medium-vac.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file medium-vac.f

C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C++ Copyright (C) 2017 Korinna C. Zapp [Korinna.Zapp@cern.ch]       ++
C++                                                                 ++
C++ This file is part of JEWEL 2.2.0                                ++
C++                                                                 ++
C++ The JEWEL homepage is jewel.hepforge.org                        ++
C++                                                                 ++
C++ The medium model was partly implemented by Jochen Klein.        ++
C++ Raghav Kunnawalkam Elayavalli helped with the implementation    ++
C++ of the V+jet processes.                                         ++
C++                                                                 ++
C++ Please follow the MCnet GUIDELINES and cite Eur.Phys.J. C74     ++
C++ (2014) no.2, 2762 [arXiv:1311.0048] for the code and            ++
C++ JHEP 1303 (2013) 080 [arXiv:1212.1599] and                      ++
C++ optionally EPJC 60 (2009) 617 [arXiv:0804.3568] for the         ++
C++ physics. The reference for V+jet processes is EPJC 76 (2016)    ++
C++ no.12 695 [arXiv:1608.03099] and for recoil effects it is       ++
C++ arXiv:1707.01539.
C++                                                                 ++
C++ JEWEL relies heavily on PYTHIA 6 for the event generation. The  ++
C++ modified version of PYTHIA 6.4.25 that is distributed with      ++
C++ JEWEL is, however, not an official PYTHIA release and must not  ++
C++ be used for anything else. Please refer to results as           ++
C++ "JEWEL+PYTHIA".                                                 ++
C++                                                                 ++
C++ JEWEL also uses code provided by S. Zhang and J. M. Jing        ++
C++ (Computation of Special Functions, John Wiley & Sons, New York, ++
C++ 1996 and http://jin.ece.illinois.edu) for computing the         ++
C++ exponential integral Ei(x).                                     ++
C++                                                                 ++
C++                                                                 ++
C++ JEWEL  is free software; you can redistribute it and/or         ++
C++ modify it under the terms of the GNU General Public License     ++
C++ as published by the Free Software Foundation; either version 2  ++
C++ of the License, or (at your option) any later version.          ++
C++                                                                 ++
C++ JEWEL is distributed in the hope that it will be useful,        ++
C++ but WITHOUT ANY WARRANTY; without even the implied warranty of  ++
C++ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the    ++
C++ GNU General Public License for more details.                    ++
C++                                                                 ++
C++ You should have received a copy of the GNU General Public       ++  
C++ License along with this program; if not, write to the Free      ++
C++ Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, ++
C++ MA 02110-1301 USA                                               ++
C++                                                                 ++
C++ Linking JEWEL statically or dynamically with other modules is   ++
C++ making a combined work based on JEWEL. Thus, the terms and      ++
C++ conditions of the GNU General Public License cover the whole    ++
C++ combination.                                                    ++
C++                                                                 ++
C++ In addition, as a special exception, I give you permission to   ++
C++ combine JEWEL with the code for the computation of special      ++
C++ functions provided by S. Zhang and J. M. Jing. You may copy and ++
C++ distribute such a system following the terms of the GNU GPL for ++
C++ JEWEL and the licenses of the other code concerned, provided    ++
C++ that you include the source code of that other code when and as ++
C++ the GNU GPL requires distribution of source code.               ++
C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

      SUBROUTINE medinit
      IMPLICIT NONE
      common/medparam/bmin,bmax,centrmin,centrmax,breal,rau
      DOUBLE PRECISION bmin,bmax,centrmin,centrmax,breal,rau
      common/medium/medium
      LOGICAL medium
      DATA medium/.false./
      rau=5.d0
      END



      SUBROUTINE mednextevt
      IMPLICIT NONE
      END



      double precision function getcentrality()
      implicit none
      getcentrality=-1.d0
      end



      SUBROUTINE pickvtx(X,Y)
      IMPLICIT NONE
      DOUBLE PRECISION x,y
      x=0.d0
      y=0.d0
      END



        SUBROUTINE setb(BVAL)
        IMPLICIT NONE
        DOUBLE PRECISION bval
        END



      SUBROUTINE getscatterer(X,Y,Z,T,TYPE,PX,PY,PZ,E,MS,MD,TEMP)
      IMPLICIT NONE
      DOUBLE PRECISION x,y,z,t,ms,px,py,pz,e,md,temp
      INTEGER type
      WRITE(*,*)'GETSCATTERER called although in vacuum'
      END



      SUBROUTINE avscatcen(X,Y,Z,T,PX,PY,PZ,E,MS)
      IMPLICIT NONE
      DOUBLE PRECISION x,y,z,t,ms,px,py,pz,e
      WRITE(*,*)'AVSCATCEN called although in vacuum'
      END



      SUBROUTINE maxscatcen(PX,PY,PZ,E,MS)
      IMPLICIT NONE
      DOUBLE PRECISION ms,px,py,pz,e
      ms=0.5d0
        e=0.5d0
        px=0.d0
        py=0.d0
        pz=0.d0
      END



      DOUBLE PRECISION FUNCTION getneff(X,Y,Z,T)
      IMPLICIT NONE
      DOUBLE PRECISION x,y,z,t
      getneff=0.d0
      END



      DOUBLE PRECISION FUNCTION gettemp(X,Y,Z,T)
      IMPLICIT NONE
      DOUBLE PRECISION x,y,z,t
      gettemp=0.d0
      END


      
      DOUBLE PRECISION FUNCTION gettempmax()
      IMPLICIT NONE
      gettempmax=0.d0
      END



      DOUBLE PRECISION FUNCTION getmdmax()
      IMPLICIT NONE
      getmdmax=0.d0
      END



      DOUBLE PRECISION FUNCTION getmdmin()
      IMPLICIT NONE
      getmdmin=0.d0
      END



      DOUBLE PRECISION FUNCTION getmsmax()
      IMPLICIT NONE
      getmsmax=0.d0
      END



      DOUBLE PRECISION FUNCTION getneffmax()
      IMPLICIT NONE
      getneffmax=0.d0
      END



        DOUBLE PRECISION FUNCTION getnatmdmin()
        IMPLICIT NONE
      getnatmdmin=0.d0
        END



      DOUBLE PRECISION FUNCTION getltimemax()
      IMPLICIT NONE
        getltimemax=0.
        END



      DOUBLE PRECISION FUNCTION getmd(X,Y,Z,T)
      IMPLICIT NONE
      DOUBLE PRECISION x,y,z,t,gettemp
      getmd=0.d0
      END



      DOUBLE PRECISION FUNCTION getms(X,Y,Z,T)
      IMPLICIT NONE
      DOUBLE PRECISION x,y,z,t,getmd
      getms=0.d0
      END



      DOUBLE PRECISION FUNCTION medderiv(XVAL,W)
      medderiv=0.d0
      END