EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
demo.F
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file demo.F
1 C***********************************************************************
2 C...
3  PROGRAM demo
4  IMPLICIT NONE
5  INTEGER nset
6  parameter(nset=20)
7 C...Demonstration job for LEPTO
8  COMMON /leptou/ cut(14),lst(40),parl(30),x,y,w2,q2,u
9  REAL cut,parl,x,y,w2,q2,u
10  INTEGER lst
11  common/lujets/n,k(4000,5),p(4000,5),v(4000,5)
12  REAL p,v
13  INTEGER n,k
14  INTEGER nevent
15  DATA nevent/100/
16 
17  INTEGER ne,set
18  REAL t1,t2,t3
19 
20  INTEGER lst7(nset),lst8(nset),lst11(nset),lst15(nset),lst19(nset),
21  & lst20(nset),
22  & lst23(nset),lst34(nset),lst35(nset)
23  INTEGER ltype(nset)
24 
25  DATA lst7/1,0,18*1/
26  DATA lst8/2*12,1,17*12/
27  DATA lst11/3*0,111,16*0/
28  DATA lst15/17*9,-4,-4,9/
29  DATA lst19/4*-1,10,-10,14*-1/
30  DATA lst20/6*5,1,6,12*5/
31  DATA lst23/8*4,1,2,3,9*4/
32  DATA lst34/11*1,0,0,7*1/
33  DATA lst35/12*1,0,0,3*1,0,0,1/
34  DATA ltype/9*-11,13,10*-11/
35 
36  REAL parl7(nset),parl8(nset),parl9(nset)
37  REAL pe(nset),le(nset)
38 
39  DATA parl7/20*0.5/
40  DATA parl8/6*0.04,2*0.005,6*0.04,0.001,5*0.04/
41  DATA parl9/15*4.0,10.0,100.0,3*4.0/
42  DATA pe/20*820.0/
43  DATA le/20*-27.5/
44 
45 C...Set cuts
46  lst(3)=0
47  cut(7)=100.
48 
49  DO 600 set=1,nset
50  lst(7)=lst7(set)
51  lst(8)=lst8(set)
52  lst(11)=lst11(set)
53  lst(15)=lst15(set)
54  lst(19)=lst19(set)
55  lst(20)=lst20(set)
56  lst(23)=lst23(set)
57  lst(34)=lst34(set)
58  lst(35)=lst35(set)
59  parl(7)=parl7(set)
60  parl(8)=parl8(set)
61  parl(9)=parl9(set)
62 
63  WRITE(6,*) 'LST(7)=',lst(7),'hadronisation'
64  WRITE(6,*) 'LST(8)=',lst(8),'QCD'
65  WRITE(6,*) 'LST(11)=',lst(11),'LQCD+10*LTM+100*LHT'
66  WRITE(6,*) 'LST(19)=',lst(19),'Grid'
67  WRITE(6,*) 'LST(20)=',lst(20),'Cut-off scheme'
68  WRITE(6,*) 'LST(34)=',lst(34),'SCI'
69  WRITE(6,*) 'LST(35)=',lst(35),'SQT'
70  WRITE(6,*) 'PARL(7)=',parl(7),'SCI prob'
71  WRITE(6,*) 'PARL(8)=',parl(8),'ycut or zmin'
72  WRITE(6,*) 'PARL(9)=',parl(9),'smin'
73  WRITE(6,*) 'LEPTON TYPE=',ltype(set)
74  WRITE(6,*) 'LEPTON ENERGY=',le(set)
75  WRITE(6,*) 'PROTON ENERGY=',pe(set)
76 
77  CALL timex(t1)
78  CALL linit(0,ltype(set),le(set),pe(set),lst23(set))
79  CALL timex(t2)
80 
81  DO 500 ne=1,nevent
82  100 CALL lepto
83  IF(lst(21).NE.0) THEN
84  WRITE(6,*) ' Warning: LST(21)=',lst(21),' event skipped'
85  goto 100
86  ENDIF
87  IF(mod(ne,1000).EQ.0) THEN
88 C CALL LULIST(1)
89 C WRITE(6,*) 'N,x,y,Q2,W= ',NE,X,Y,Q2,SQRT(W2)
90  ENDIF
91  500 CONTINUE
92  CALL timex(t3)
93  WRITE(6,1000) parl(23),parl(24),nevent,
94  & t2-t1,t3-t2,(t3-t2)/max(1,nevent)
95  600 CONTINUE
96  1000 FORMAT(/,' Cross section from numerical integration:',g12.3,' pb',
97  & /,' Cross section from MonteCarlo simulation:',g12.3,' pb',
98  & /,' Number of events generated: ',i12,
99  & /,' Time for initialization (LINIT): ',f12.3,' s',
100  & /,' Time for event generation, total: ',f12.3,' s',
101  & /,' per event: ',f12.3,' s')
102  END