[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
gmu.in
- To: Peter Kammel <kammel@npl.uiuc.edu>
- Subject: gmu.in
- From: Bernhard LAUSS <lauss@ill.fr>
- Date: Fri, 20 Oct 2006 12:34:09 +0200
- In-reply-to: <Pine.LNX.4.21.0609281409580.27419-100000@three.npl.uiuc.edu>
- References: <Pine.LNX.4.21.0609281409580.27419-100000@three.npl.uiuc.edu>
- User-agent: Thunderbird 1.5.0.4 (Windows/20060516)
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