3 C*********************************************************************
5 SUBROUTINE LUTEST_HIJING(MTEST)
7 C...Purpose: to provide a simple program (disguised as subroutine) to
8 C...run at installation as a check that the program works as intended.
9 #include "lujets_hijing.inc"
10 #include "ludat1_hijing.inc"
11 DIMENSION PSUM(5),PINI(6),PFIN(6)
13 C...Loop over events to be generated.
14 IF(MTEST.GE.1) CALL LUTABU_HIJING(20)
18 C...Reset parameter values. Switch on some nonstandard features.
33 IF(IEV.EQ.301.OR.IEV.EQ.351.OR.IEV.EQ.401) MSTJ(116)=3
35 C...Ten events each for some single jets configurations.
39 IF(ITY.EQ.3.OR.ITY.EQ.4) MSTJ(11)=2
40 IF(ITY.EQ.1) CALL LU1ENT_HIJING(1,1,15.,0.,0.)
41 IF(ITY.EQ.2) CALL LU1ENT_HIJING(1,3101,15.,0.,0.)
42 IF(ITY.EQ.3) CALL LU1ENT_HIJING(1,-2203,15.,0.,0.)
43 IF(ITY.EQ.4) CALL LU1ENT_HIJING(1,-4,30.,0.,0.)
44 IF(ITY.EQ.5) CALL LU1ENT_HIJING(1,21,15.,0.,0.)
46 C...Ten events each for some simple jet systems; string fragmentation.
47 ELSEIF(IEV.LE.130) THEN
49 IF(ITY.EQ.1) CALL LU2ENT_HIJING(1,1,-1,40.)
50 IF(ITY.EQ.2) CALL LU2ENT_HIJING(1,4,-4,30.)
51 IF(ITY.EQ.3) CALL LU2ENT_HIJING(1,2,2103,100.)
52 IF(ITY.EQ.4) CALL LU2ENT_HIJING(1,21,21,40.)
53 IF(ITY.EQ.5) CALL LU3ENT_HIJING(1,2101,21,-3203,30.,0.6,0.8)
54 IF(ITY.EQ.6) CALL LU3ENT_HIJING(1,5,21,-5,40.,0.9,0.8)
55 IF(ITY.EQ.7) CALL LU3ENT_HIJING(1,21,21,21,60.,0.7,0.5)
56 IF(ITY.EQ.8) CALL LU4ENT_HIJING(1,2,21,21,-2,40.,0.4,0.64,0.6,0
59 C...Seventy events with independent fragmentation and momentum cons.
60 ELSEIF(IEV.LE.200) THEN
62 MSTJ(2)=1+MOD(IEV-131,4)
63 MSTJ(3)=1+MOD((IEV-131)/4,4)
64 IF(ITY.EQ.1) CALL LU2ENT_HIJING(1,4,-5,40.)
65 IF(ITY.EQ.2) CALL LU3ENT_HIJING(1,3,21,-3,40.,0.9,0.4)
66 IF(ITY.EQ.3) CALL LU4ENT_HIJING(1,2,21,21,-2,40.,0.4,0.64,0.6,0
68 IF(ITY.GE.4) CALL LU4ENT_HIJING(1,2,-3,3,-2,40.,0.4,0.64,0.6,0
71 C...A hundred events with random jets (check invariant mass).
72 ELSEIF(IEV.LE.300) THEN
75 NJET=2.+6.*RLU_HIJING(0)
78 IF(I.EQ.1) KFL=INT(1.+4.*RLU_HIJING(0))
79 IF(I.EQ.NJET) KFL=-INT(1.+4.*RLU_HIJING(0))
80 EJET=5.+20.*RLU_HIJING(0)
81 THETA=ACOS(2.*RLU_HIJING(0)-1.)
82 PHI=6.2832*RLU_HIJING(0)
83 IF(I.LT.NJET) CALL LU1ENT_HIJING(-I,KFL,EJET,THETA,PHI)
84 IF(I.EQ.NJET) CALL LU1ENT_HIJING(I,KFL,EJET,THETA,PHI)
85 IF(I.EQ.1.OR.I.EQ.NJET) PSUM(5)=PSUM(5)+ULMASS_HIJING(KFL)
87 120 PSUM(J)=PSUM(J)+P(I,J)
88 IF(PSUM(4)**2-PSUM(1)**2-PSUM(2)**2-PSUM(3)**2.LT.
89 & (PSUM(5)+PARJ(32))**2) GOTO 100
91 C...Fifty e+e- continuum events with matrix elements.
92 ELSEIF(IEV.LE.350) THEN
94 CALL LUEEVT_HIJING(0,40.)
96 C...Fifty e+e- continuum event with varying shower options.
97 ELSEIF(IEV.LE.400) THEN
99 MSTJ(43)=1+MOD(IEV/2,4)
100 MSTJ(44)=MOD(IEV/8,3)
101 CALL LUEEVT_HIJING(0,90.)
103 C...Fifty e+e- continuum events with coherent shower, including top.
104 ELSEIF(IEV.LE.450) THEN
106 CALL LUEEVT_HIJING(0,500.)
108 C...Fifty Upsilon decays to ggg or gammagg with coherent shower.
109 ELSEIF(IEV.LE.500) THEN
110 CALL LUONIA_HIJING(5,9.46)
112 C...One decay each for some heavy mesons.
113 ELSEIF(IEV.LE.560) THEN
118 CALL LU1ENT_HIJING(1,100*KFLB+10*KFLC+KFLS,0.,0.,0.)
120 C...One decay each for some heavy baryons.
121 ELSEIF(IEV.LE.600) THEN
127 CALL LU1ENT_HIJING(1,1000*KFLA+100*KFLB+10*KFLC+KFLS,0.,0.,0.)
130 C...Generate event. Find total momentum, energy and charge.
132 130 PINI(J)=PLU_HIJING(0,J)
133 PINI(6)=PLU_HIJING(0,6)
136 140 PFIN(J)=PLU_HIJING(0,J)
137 PFIN(6)=PLU_HIJING(0,6)
139 C...Check conservation of energy, momentum and charge;
140 C...usually exact, but only approximate for single jets.
143 IF((PFIN(1)-PINI(1))**2+(PFIN(2)-PINI(2))**2.GE.4.) MERR=MERR+1
144 EPZREM=PINI(4)+PINI(3)-PFIN(4)-PFIN(3)
145 IF(EPZREM.LT.0..OR.EPZREM.GT.2.*PARJ(31)) MERR=MERR+1
146 IF(ABS(PFIN(6)-PINI(6)).GT.2.1) MERR=MERR+1
149 150 IF(ABS(PFIN(J)-PINI(J)).GT.0001*PINI(4)) MERR=MERR+1
150 IF(ABS(PFIN(6)-PINI(6)).GT.0.1) MERR=MERR+1
152 IF(MERR.NE.0) WRITE(MSTU(11),1000) (PINI(J),J=1,4),PINI(6),
153 &(PFIN(J),J=1,4),PFIN(6)
155 C...Check that all KF codes are known ones, and that partons/particles
156 C...satisfy energy-momentum-mass relation. Store particle statistics.
158 IF(K(I,1).GT.20) GOTO 160
159 IF(LUCOMP_HIJING(K(I,2)).EQ.0) THEN
160 WRITE(MSTU(11),1100) I
163 PD=P(I,4)**2-P(I,1)**2-P(I,2)**2-P(I,3)**2-P(I,5)**2
164 IF(ABS(PD).GT.MAX(0.1,0.001*P(I,4)**2).OR.P(I,4).LT.0.) THEN
165 WRITE(MSTU(11),1200) I
169 IF(MTEST.GE.1) CALL LUTABU_HIJING(21)
171 C...List all erroneous events and some normal ones.
172 IF(MERR.NE.0.OR.MSTU(24).NE.0.OR.MSTU(28).NE.0) THEN
173 CALL LULIST_HIJING(2)
174 ELSEIF(MTEST.GE.1.AND.MOD(IEV-5,100).EQ.0) THEN
175 CALL LULIST_HIJING(1)
178 C...Stop execution if too many errors. Endresult of run.
179 IF(MERR.NE.0) NERR=NERR+1
181 WRITE(MSTU(11),1300) IEV
185 IF(MTEST.GE.1) CALL LUTABU_HIJING(22)
186 WRITE(MSTU(11),1400) NERR
188 C...Reset commonblock variables changed during run.
197 C...Format statements for output.
198 1000 FORMAT(/' Momentum, energy and/or charge were not conserved ',
199 &'in following event'/' sum of',9X,'px',11X,'py',11X,'pz',11X,
200 &'E',8X,'charge'/' before',2X,4(1X,F12.5),1X,F8.2/' after',3X,
201 &4(1X,F12.5),1X,F8.2)
202 1100 FORMAT(/5X,'Entry no.',I4,' in following event not known code')
203 1200 FORMAT(/5X,'Entry no.',I4,' in following event has faulty ',
205 1300 FORMAT(/5X,'Ten errors experienced by event ',I3/
206 &5X,'Something is seriously wrong! Execution stopped now!')
207 1400 FORMAT(/5X,'Number of erroneous or suspect events in run:',I3/
208 &5X,'(0 fine, 1 acceptable if a single jet, ',
209 &'>=2 something is wrong)')