ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
medium-vac.f
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file medium-vac.f
1 C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2 C++ Copyright (C) 2017 Korinna C. Zapp [Korinna.Zapp@cern.ch] ++
3 C++ ++
4 C++ This file is part of JEWEL 2.2.0 ++
5 C++ ++
6 C++ The JEWEL homepage is jewel.hepforge.org ++
7 C++ ++
8 C++ The medium model was partly implemented by Jochen Klein. ++
9 C++ Raghav Kunnawalkam Elayavalli helped with the implementation ++
10 C++ of the V+jet processes. ++
11 C++ ++
12 C++ Please follow the MCnet GUIDELINES and cite Eur.Phys.J. C74 ++
13 C++ (2014) no.2, 2762 [arXiv:1311.0048] for the code and ++
14 C++ JHEP 1303 (2013) 080 [arXiv:1212.1599] and ++
15 C++ optionally EPJC 60 (2009) 617 [arXiv:0804.3568] for the ++
16 C++ physics. The reference for V+jet processes is EPJC 76 (2016) ++
17 C++ no.12 695 [arXiv:1608.03099] and for recoil effects it is ++
18 C++ arXiv:1707.01539.
19 C++ ++
20 C++ JEWEL relies heavily on PYTHIA 6 for the event generation. The ++
21 C++ modified version of PYTHIA 6.4.25 that is distributed with ++
22 C++ JEWEL is, however, not an official PYTHIA release and must not ++
23 C++ be used for anything else. Please refer to results as ++
24 C++ "JEWEL+PYTHIA". ++
25 C++ ++
26 C++ JEWEL also uses code provided by S. Zhang and J. M. Jing ++
27 C++ (Computation of Special Functions, John Wiley & Sons, New York, ++
28 C++ 1996 and http://jin.ece.illinois.edu) for computing the ++
29 C++ exponential integral Ei(x). ++
30 C++ ++
31 C++ ++
32 C++ JEWEL is free software; you can redistribute it and/or ++
33 C++ modify it under the terms of the GNU General Public License ++
34 C++ as published by the Free Software Foundation; either version 2 ++
35 C++ of the License, or (at your option) any later version. ++
36 C++ ++
37 C++ JEWEL is distributed in the hope that it will be useful, ++
38 C++ but WITHOUT ANY WARRANTY; without even the implied warranty of ++
39 C++ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the ++
40 C++ GNU General Public License for more details. ++
41 C++ ++
42 C++ You should have received a copy of the GNU General Public ++
43 C++ License along with this program; if not, write to the Free ++
44 C++ Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, ++
45 C++ MA 02110-1301 USA ++
46 C++ ++
47 C++ Linking JEWEL statically or dynamically with other modules is ++
48 C++ making a combined work based on JEWEL. Thus, the terms and ++
49 C++ conditions of the GNU General Public License cover the whole ++
50 C++ combination. ++
51 C++ ++
52 C++ In addition, as a special exception, I give you permission to ++
53 C++ combine JEWEL with the code for the computation of special ++
54 C++ functions provided by S. Zhang and J. M. Jing. You may copy and ++
55 C++ distribute such a system following the terms of the GNU GPL for ++
56 C++ JEWEL and the licenses of the other code concerned, provided ++
57 C++ that you include the source code of that other code when and as ++
58 C++ the GNU GPL requires distribution of source code. ++
59 C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
60 
61  SUBROUTINE medinit
62  IMPLICIT NONE
63  common/medparam/bmin,bmax,centrmin,centrmax,breal,rau
64  DOUBLE PRECISION bmin,bmax,centrmin,centrmax,breal,rau
66  LOGICAL medium
67  DATA medium/.false./
68  rau=5.d0
69  END
70 
71 
72 
73  SUBROUTINE mednextevt
74  IMPLICIT NONE
75  END
76 
77 
78 
79  double precision function getcentrality()
80  implicit none
81  getcentrality=-1.d0
82  end
83 
84 
85 
86  SUBROUTINE pickvtx(X,Y)
87  IMPLICIT NONE
88  DOUBLE PRECISION x,y
89  x=0.d0
90  y=0.d0
91  END
92 
93 
94 
95  SUBROUTINE setb(BVAL)
96  IMPLICIT NONE
97  DOUBLE PRECISION bval
98  END
99 
100 
101 
102  SUBROUTINE getscatterer(X,Y,Z,T,TYPE,PX,PY,PZ,E,MS,MD,TEMP)
103  IMPLICIT NONE
104  DOUBLE PRECISION x,y,z,t,ms,px,py,pz,e,md,temp
105  INTEGER type
106  WRITE(*,*)'GETSCATTERER called although in vacuum'
107  END
108 
109 
110 
111  SUBROUTINE avscatcen(X,Y,Z,T,PX,PY,PZ,E,MS)
112  IMPLICIT NONE
113  DOUBLE PRECISION x,y,z,t,ms,px,py,pz,e
114  WRITE(*,*)'AVSCATCEN called although in vacuum'
115  END
116 
117 
118 
119  SUBROUTINE maxscatcen(PX,PY,PZ,E,MS)
120  IMPLICIT NONE
121  DOUBLE PRECISION ms,px,py,pz,e
122  ms=0.5d0
123  e=0.5d0
124  px=0.d0
125  py=0.d0
126  pz=0.d0
127  END
128 
129 
130 
131  DOUBLE PRECISION FUNCTION getneff(X,Y,Z,T)
132  IMPLICIT NONE
133  DOUBLE PRECISION x,y,z,t
134  getneff=0.d0
135  END
136 
137 
138 
139  DOUBLE PRECISION FUNCTION gettemp(X,Y,Z,T)
140  IMPLICIT NONE
141  DOUBLE PRECISION x,y,z,t
142  gettemp=0.d0
143  END
144 
145 
146 
147  DOUBLE PRECISION FUNCTION gettempmax()
148  IMPLICIT NONE
149  gettempmax=0.d0
150  END
151 
152 
153 
154  DOUBLE PRECISION FUNCTION getmdmax()
155  IMPLICIT NONE
156  getmdmax=0.d0
157  END
158 
159 
160 
161  DOUBLE PRECISION FUNCTION getmdmin()
162  IMPLICIT NONE
163  getmdmin=0.d0
164  END
165 
166 
167 
168  DOUBLE PRECISION FUNCTION getmsmax()
169  IMPLICIT NONE
170  getmsmax=0.d0
171  END
172 
173 
174 
175  DOUBLE PRECISION FUNCTION getneffmax()
176  IMPLICIT NONE
177  getneffmax=0.d0
178  END
179 
180 
181 
182  DOUBLE PRECISION FUNCTION getnatmdmin()
183  IMPLICIT NONE
184  getnatmdmin=0.d0
185  END
186 
187 
188 
189  DOUBLE PRECISION FUNCTION getltimemax()
190  IMPLICIT NONE
191  getltimemax=0.
192  END
193 
194 
195 
196  DOUBLE PRECISION FUNCTION getmd(X,Y,Z,T)
197  IMPLICIT NONE
198  DOUBLE PRECISION x,y,z,t,gettemp
199  getmd=0.d0
200  END
201 
202 
203 
204  DOUBLE PRECISION FUNCTION getms(X,Y,Z,T)
205  IMPLICIT NONE
206  DOUBLE PRECISION x,y,z,t,getmd
207  getms=0.d0
208  END
209 
210 
211 
212  DOUBLE PRECISION FUNCTION medderiv(XVAL,W)
213  medderiv=0.d0
214  END