[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

gmu.in



Dear Peter,

this is the gmu.in which I belive having used latetly,
it corresponds to the one in the CVS directory on kaon.

	VLG
           Bernhard


-----------------------------------------------------------------------------


LIST
TRIG 9999999999
DEBU  1   2   10000
CUTS      0.000100 0.000100 0.0001 0.0001 0.0001 0.000100 0.000100
          0.000100 0.000100 0.01 10000. 0.
HADR      1
LOSS      1
DRAY      1
MUNU      1
COMP      1
PAIR      1
BREM      1
ANNI      1
AUTO      0
PHOT      1  
MULS      1
C-->  run 1=formatted iswit(3)=(b1=1:nowall,b2=1:noframe,b3=1:nomag)
SWIT  108 1  0  0  0  0  0   0  0  0
C-->	  ipart=ipart+100 p,x,y random  
C  KINE:    ipart    p       dp/p    dY     z	sigmax	sigmay divX	divY
KINE        106    0.03263   0.016  -0.7   -45.  0.40	0.90	 11.	 38.
C  MATE  21
BINS  20 IN Z  20 IN R
TIME  10. 10. 1.
MXEV  100000
CYCL  21
END
STOP

--------------------------------------------------------------------------------

The beam divergence in gukine.f is created as:



C add beam divergence -> create x and y momentum

C          CALL GCHEAP(RNDM(1))
C     x-Divergence PKINE(7)
C     y-Divergence PKINE(8)

          PLAB(1) = tan(PKINE(7)/1000*(RNDM(4)*2-1))*PLAB(3)
          PLAB(3) = sqrt(PLAB(3)**2 - PLAB(1)**2)

          PLAB(2) = tan(PKINE(8)/1000*(RNDM(5)*2-1))*PLAB(3)
          PLAB(3) = sqrt(PLAB(3)**2 - PLAB(2)**2)












     _/   _/         _/                 Dr. Bernhard Lauss
    _/   _/         _/                  Institut Laue-Langevin
   _/   _/         _/                   B.P.156 - 6 rue Jules Horowitz
  _/   _/         _/                    38042 Grenoble - Cedex 9
 _/   _/_/_/_/   _/_/_/_/               France

					Tel:+33-(0)4-7620 7134   
					Fax:+33-(0)4-7620 7777

					Email: lauss@ill.fr