* 
* runs Pythia+GSTAR for minb events
* 
 
 MACRO aaaa  neventk=10   file=minb200 seed=13543676

* echo seed=[seed] neventk=[nevents]
 rndm  9543713 [seed]
 detp  geom year2003 gcalor_on debug_on split_off help 

* sigma of the beam in X=Y & Z
 vsig 0.05 7.5  
* offset added to the vertex (X,Y,Z)
 gvert 0.58 0.34 0.0

 make  geometry 
 gclose all

 gfile    o [file].fzd
 
 gexec $STAF/lib/apythia.sl 

* ....................  NICE DRAWING .............. START 
 ne; 
 dcut cave x   1 10  5 .01 .01; 
 daxis 0.0 0.0 0.0 90.0

 dcut cave z   1 3 14 .01 .01; * BEMC
 daxis 0.0 0.0 0.0 90.0

 dcut cave z 270 17 14 .01 .01; *EEMC
 daxis 0.0 0.0 0.0 90.0

 dcut cave z 376 10 15 .04 .04; *BBC in  x4
 dcut cave z 373 10 15 .04 .04;   *BBC out x4
 daxis 0.0 0.0 0.0 30.0
* ....................  NICE DRAWING .............. END 


ener 200.
PMA1 6=175.
PMA1 23=91.2   24=80.5


MSEL 1      ! 1:nondifractive only, 2:all

MSTP 51=7   ! choose pdf. 7:CTEQ5L LO in Pythia 6.1
CKIN 3=-1.0 ! pt min     <<<<<<<<<< Pt MIN cut here
CKIN 8=10.  ! y max

mdcy (111,1)=0   !=0: make pi0 stable, =1 let it decay
mdcy (221,1)=0   !make eta stable
MDCY (310,1)=0   !MAKE K_SHORT STABLE
*MDCY (3112,1)=0  !MAKE SIGMA- STABLE
*MDCY (3122,1)=0  !MAKE LAMBDA0 STABLE
*MDCY (3222,1)=0  !MAKE SIGMA+ STABLE
*MDCY (3312,1)=0  !MAKE CASCADE- STABLE
*MDCY (3322,1)=0  !MAKE CASCADE0 STABLE
*MDCY (3334,1)=0  !MAKE OMEGA- STABLE


*MSTP 2=2 33=3 81=1 82=4
*MSTJ  22=2


*---------------- S I M U L A T I O N  -----------------
    
 do i = 1, [neventk]
    trig 1
*   gprint hits ecal
 enddo
 
 mess  STAF .......................................................... end
  
*  dhits;
return