source/LQGENEP/LQguser.f

0001 *-----------------
0002 * File: LQguser.f
0003 *-----------------
0004 *
0005       Program LQGUSER
0006 *******************************************
0007 * Generates 1000 events of the subprocess *
0008 *         e+ qi --> S_1/2L --> mu+ qj     *
0009 *           at the HERA energies          *
0010 *                                         *
0011 * LQ mass = 200 GeV                       *
0012 * LQ couplings = 0.3                      *
0013 * first generation quark involved         *
0014 *******************************************
0015 !      Implicit None
0016 *
0017
0018       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
0019       IMPLICIT INTEGER(I-N)
0020       INTEGER PYK,PYCHGE,PYCOMP
0021
0022
0023       double precision mass,cross_max,
0024      >  ptnt,phnt,ptt,pht,pth,phh,
0025      >  ppnt,ppt,pph,
0026      >  ptx,pty,ptz,phx,phy,phz,
0027      >  ptid,phid,ppid,pxid,pyid,pzid,lam,
0028      >  q2min
0029
0030       integer lqtype,Nevt,qi,qj,iproc,id,
0031      > Nevtdown,Nevtup
0032       logical first
0033       data first /.true./
0034       save first
0035 *
0036
0037 C...Pythia Ranom Seed Test
0038       COMMON/PYDATR/MRPY(6),RRPY(100)
0039       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
0040
0041 C...LQGENEP run setup parameters
0042       double precision BEAMPAR,LQGPAR3
0043       integer LQGPAR1,LQGPAR2
0044       CHARACTER*256 out_file
0045       character*4  j_string
0046       COMMON/LQGDATA/BEAMPAR(3),LQGPAR1(10),LQGPAR2(10),LQGPAR3(20),
0047      > out_file
0048
0049 C...LQGENEP event informations
0050       double precision LQGKPAR,LQGDDX
0051       integer LQGPID
0052       COMMON/LQGEVT/LQGKPAR(3),LQGDDX(3),LQGPID(3)
0053 *
0054       integer NeV,i,j,N_tot,N_div
0055 *
0056
0057
0058       OPEN(UNIT=22,FILE='inputfile',STATUS='OLD')
0059
0060          Read(22,*)Nevt,Mass,lqtype,beampar(2),beampar(3),
0061      >    q2min,xmin,xmax,ymin,ymax,qi,qj
0062          Read(22,*)out_file
0063          close(22)
0064
0065 c     set output file name
0066 c         out_file = 'LQGENEP_output.txt'
0067
0068 c     set random seed
0069          MRPY(1) = 1*10000
0070 c      Mass=1936.D0
0071 c      qi=1
0072 c      qj=1
0073       iproc=0
0074 c      cross_max=2.d-5
0075       cross_max=7.d-15
0076 *
0077 * Charge of lepton beam (+1/-1 positron/electron)
0078 c      beampar(1)=+1.  ! positron
0079       beampar(1)=-1.  ! electron
0080 * Energy of lepton beam
0081 c      beampar(2)=10.     ! EIC1
0082 c      beampar(2)=20.    ! EIC2
0083 c       beampar(2)=27.5  ! HERA
0084 * Energy of proton beam
0085 c      beampar(3)=250.   ! EIC1
0086 c      beampar(3)=325.   ! EIC2
0087 c      beampar(3)=820.   ! HERA
0088
0089 * Number of events to generate
0090       lqgpar1(1)=Nevt
0091 * Number of events which will be fully listed in the output
0092       lqgpar1(2)=0
0093 * Histogram flag (if 1 some histograms will be produced)
0094       lqgpar1(3)=1
0095 * Current event number (first generated event will be lqgpar1(4)+1)
0096       lqgpar1(4)=0
0097 * Pythia process number (should be > 400, if = 0 it will be set
0098 * to 401, first value available for external processes)
0099       lqgpar1(5)=iproc
0100 * LQ type
0101       lqgpar2(1)=lqtype
0102 * generation of the input quark in the s-channel process
0103       lqgpar2(2)=qi
0104 * generation of the output quark in the s-channel process
0105       lqgpar2(3)=qj
0106 * generation of the output lepton
0107       lqgpar2(4)=3
0108 * LQ mass in GeV
0109       lqgpar3(1)=Mass
0110 * Initial state s-channel coupling
0111       lqgpar3(2)=0.3
0112 * Final state s-channel coupling (in case of process eq -> LQ -> eq
0113 *                               the two couplings should be the same)
0114       lqgpar3(3)=0.3
0115 * x range low limit
0116       lqgpar3(4)=xmin
0117 * x range high limit
0118       lqgpar3(5)=xmax
0119 * y range low limit
0120       lqgpar3(6)=ymin
0121 * y range high limit
0122       lqgpar3(7)=ymax
0123 * minimum allowed Q^2 in Gev^2
0124       lqgpar3(8)=q2min
0125 * Parton distribution type according to PDFLIB
0126       lqgpar3(9)=1.
0127 * Parton distribution group according to PDFLIB
0128       lqgpar3(10)=4.
0129 * Parton distribution set according to PDFLIB
0130       lqgpar3(11)=32.
0131 * Max cross section
0132       lqgpar3(12)=cross_max
0133 * Eventually switch off initial state QCD and QED radiation
0134 * setting to 0 the following Pythia parameters
0135 c     MSTP(61)=0
0136 c     MSTP(71)=0
0137 * Eventually switch off multiple interaction
0138 c     MSTP(81)=0
0139 * Eventually switch off fragmentation and decay
0140 c     MSTP(111)=0
0141
0142 * LQGENEP Initialization
0143       call LQGENEP(Nevt,0)
0144
0145       Nev=lqgpar1(1)
0146
0147 * LQGENEP generation loop
0148
0149       do i=1,Nevt
0150 cc       if (i.NE.Nevt) then
0151 c        if ((i.GT.(Nevt-1000)).and.
0152 c     >    (i.LE.Nevt)) then
0153 cc       CALL PYEVNT
0154 cc      else
0155        call LQGENEP(Nevt,1)
0156 call LQGENEP(Nevt,0)
0157 cc      endif
0158       enddo
0159
0160
0161 * LQGENEP termination
0162        call LQGENEP(Nevt,2)
0163
0164        stop
0165       end
0166
0167
0168