Back to home page

sPhenix code displayed by LXR

 
 

    


File indexing completed on 2025-08-03 08:16:23

0001 C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
0002 C++ Copyright (C) 2017 Korinna C. Zapp [Korinna.Zapp@cern.ch]       ++
0003 C++                                                                 ++
0004 C++ This file is part of JEWEL 2.2.0                                ++
0005 C++                                                                 ++
0006 C++ The JEWEL homepage is jewel.hepforge.org                        ++
0007 C++                                                                 ++
0008 C++ The medium model was partly implemented by Jochen Klein.        ++
0009 C++ Raghav Kunnawalkam Elayavalli helped with the implementation    ++
0010 C++ of the V+jet processes.                                         ++
0011 C++                                                                 ++
0012 C++ Please follow the MCnet GUIDELINES and cite Eur.Phys.J. C74     ++
0013 C++ (2014) no.2, 2762 [arXiv:1311.0048] for the code and            ++
0014 C++ JHEP 1303 (2013) 080 [arXiv:1212.1599] and                      ++
0015 C++ optionally EPJC 60 (2009) 617 [arXiv:0804.3568] for the         ++
0016 C++ physics. The reference for V+jet processes is EPJC 76 (2016)    ++
0017 C++ no.12 695 [arXiv:1608.03099] and for recoil effects it is       ++
0018 C++ arXiv:1707.01539.
0019 C++                                                                 ++
0020 C++ JEWEL relies heavily on PYTHIA 6 for the event generation. The  ++
0021 C++ modified version of PYTHIA 6.4.25 that is distributed with      ++
0022 C++ JEWEL is, however, not an official PYTHIA release and must not  ++
0023 C++ be used for anything else. Please refer to results as           ++
0024 C++ "JEWEL+PYTHIA".                                                 ++
0025 C++                                                                 ++
0026 C++ JEWEL also uses code provided by S. Zhang and J. M. Jing        ++
0027 C++ (Computation of Special Functions, John Wiley & Sons, New York, ++
0028 C++ 1996 and http://jin.ece.illinois.edu) for computing the         ++
0029 C++ exponential integral Ei(x).                                     ++
0030 C++                                                                 ++
0031 C++                                                                 ++
0032 C++ JEWEL  is free software; you can redistribute it and/or         ++
0033 C++ modify it under the terms of the GNU General Public License     ++
0034 C++ as published by the Free Software Foundation; either version 2  ++
0035 C++ of the License, or (at your option) any later version.          ++
0036 C++                                                                 ++
0037 C++ JEWEL is distributed in the hope that it will be useful,        ++
0038 C++ but WITHOUT ANY WARRANTY; without even the implied warranty of  ++
0039 C++ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the    ++
0040 C++ GNU General Public License for more details.                    ++
0041 C++                                                                 ++
0042 C++ You should have received a copy of the GNU General Public       ++  
0043 C++ License along with this program; if not, write to the Free      ++
0044 C++ Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, ++
0045 C++ MA 02110-1301 USA                                               ++
0046 C++                                                                 ++
0047 C++ Linking JEWEL statically or dynamically with other modules is   ++
0048 C++ making a combined work based on JEWEL. Thus, the terms and      ++
0049 C++ conditions of the GNU General Public License cover the whole    ++
0050 C++ combination.                                                    ++
0051 C++                                                                 ++
0052 C++ In addition, as a special exception, I give you permission to   ++
0053 C++ combine JEWEL with the code for the computation of special      ++
0054 C++ functions provided by S. Zhang and J. M. Jing. You may copy and ++
0055 C++ distribute such a system following the terms of the GNU GPL for ++
0056 C++ JEWEL and the licenses of the other code concerned, provided    ++
0057 C++ that you include the source code of that other code when and as ++
0058 C++ the GNU GPL requires distribution of source code.               ++
0059 C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
0060 
0061       SUBROUTINE MEDINIT
0062       IMPLICIT NONE
0063       COMMON/MEDPARAM/BMIN,BMAX,CENTRMIN,CENTRMAX,BREAL,RAU
0064       DOUBLE PRECISION BMIN,BMAX,CENTRMIN,CENTRMAX,BREAL,RAU
0065       COMMON/MEDIUM/MEDIUM
0066       LOGICAL MEDIUM
0067       DATA MEDIUM/.FALSE./
0068       RAU=5.d0
0069       END
0070 
0071 
0072 
0073       SUBROUTINE MEDNEXTEVT
0074       IMPLICIT NONE
0075       END
0076 
0077 
0078 
0079       double precision function getcentrality()
0080       implicit none
0081       getcentrality=-1.d0
0082       end
0083 
0084 
0085 
0086       SUBROUTINE PICKVTX(X,Y)
0087       IMPLICIT NONE
0088       DOUBLE PRECISION X,Y
0089       X=0.d0
0090       Y=0.d0
0091       END
0092 
0093 
0094 
0095         SUBROUTINE SETB(BVAL)
0096         IMPLICIT NONE
0097         DOUBLE PRECISION BVAL
0098         END
0099 
0100 
0101 
0102       SUBROUTINE GETSCATTERER(X,Y,Z,T,TYPE,PX,PY,PZ,E,MS,MD,TEMP)
0103       IMPLICIT NONE
0104       DOUBLE PRECISION X,Y,Z,T,MS,PX,PY,PZ,E,MD,TEMP
0105       INTEGER TYPE
0106       WRITE(*,*)'GETSCATTERER called although in vacuum'
0107       END
0108 
0109 
0110 
0111       SUBROUTINE AVSCATCEN(X,Y,Z,T,PX,PY,PZ,E,MS)
0112       IMPLICIT NONE
0113       DOUBLE PRECISION X,Y,Z,T,MS,PX,PY,PZ,E
0114       WRITE(*,*)'AVSCATCEN called although in vacuum'
0115       END
0116 
0117 
0118 
0119       SUBROUTINE MAXSCATCEN(PX,PY,PZ,E,MS)
0120       IMPLICIT NONE
0121       DOUBLE PRECISION MS,PX,PY,PZ,E
0122       ms=0.5d0
0123         e=0.5d0
0124         px=0.d0
0125         py=0.d0
0126         pz=0.d0
0127       END
0128 
0129 
0130 
0131       DOUBLE PRECISION FUNCTION GETNEFF(X,Y,Z,T)
0132       IMPLICIT NONE
0133       DOUBLE PRECISION X,Y,Z,T
0134       GETNEFF=0.d0
0135       END
0136 
0137 
0138 
0139       DOUBLE PRECISION FUNCTION GETTEMP(X,Y,Z,T)
0140       IMPLICIT NONE
0141       DOUBLE PRECISION X,Y,Z,T
0142       GETTEMP=0.d0
0143       END
0144 
0145 
0146       
0147       DOUBLE PRECISION FUNCTION GETTEMPMAX()
0148       IMPLICIT NONE
0149       GETTEMPMAX=0.d0
0150       END
0151 
0152 
0153 
0154       DOUBLE PRECISION FUNCTION GETMDMAX()
0155       IMPLICIT NONE
0156       GETMDMAX=0.d0
0157       END
0158 
0159 
0160 
0161       DOUBLE PRECISION FUNCTION GETMDMIN()
0162       IMPLICIT NONE
0163       GETMDMIN=0.d0
0164       END
0165 
0166 
0167 
0168       DOUBLE PRECISION FUNCTION GETMSMAX()
0169       IMPLICIT NONE
0170       GETMSMAX=0.d0
0171       END
0172 
0173 
0174 
0175       DOUBLE PRECISION FUNCTION GETNEFFMAX()
0176       IMPLICIT NONE
0177       GETNEFFMAX=0.d0
0178       END
0179 
0180 
0181 
0182         DOUBLE PRECISION FUNCTION GETNATMDMIN()
0183         IMPLICIT NONE
0184       GETNATMDMIN=0.d0
0185         END
0186 
0187 
0188 
0189       DOUBLE PRECISION FUNCTION GETLTIMEMAX()
0190       IMPLICIT NONE
0191         GETLTIMEMAX=0.
0192         END
0193 
0194 
0195 
0196       DOUBLE PRECISION FUNCTION GETMD(X,Y,Z,T)
0197       IMPLICIT NONE
0198       DOUBLE PRECISION X,Y,Z,T,GETTEMP
0199       GETMD=0.D0
0200       END
0201 
0202 
0203 
0204       DOUBLE PRECISION FUNCTION GETMS(X,Y,Z,T)
0205       IMPLICIT NONE
0206       DOUBLE PRECISION X,Y,Z,T,GETMD
0207       GETMS=0.D0
0208       END
0209 
0210 
0211 
0212       DOUBLE PRECISION FUNCTION MEDDERIV(XVAL,W)
0213       MEDDERIV=0.D0
0214       END