COSMOS v7.655  COSMOSv7655
(AirShowerMC)
Zobsp.h
Go to the documentation of this file.
1 ! Zobsvp.h---parameters to be given by input.
2 ! This must be preceded by Zobs.h
3 
4 ! (-> ---------------------------------------------------
5 
6  real*8 HeightList !1 Height of observation levels in m. This is made from DepthList internally.
7  ! I.e., this one is usually not an input. However, if the DepthList values are
8  ! negative, this is used as input and corresponding DepthList is computed internally.
9  real*8 DepthList !1 Depth List of Observation level in kg/m$^2$. If $< 0$, HeightList has priority.
10  ! (See HeightList)
11  real*8 ASHeightList !1 This is HeightList for Air Shower observ. Used only if Generate contains
12  ! "as". See HeightList.
13  real*8 ASDepthList !1 This is DepthList for AS observation. Used only if Generate contains
14  ! "as". See DepthList.
15  real*8 LatitOfSite !1 Latitude of the deepest observation level in degree. East is positive.
16  real*8 LongitOfSite !1 Longitude of the deepest observation level in degree. North is positive.
17 
18  real*8 DtGMT !1 Difference of the local time of the observation place from GMT (hour).
19  real*8 YearOfGeomag !1 Like 1999.5. Year when Geomagnetic field is to be calculated.
20  integer ObsPlane !1 How to observe particles. \newline
21  ! 0$ \Rightarrow $ no detector plane is used for observation. BorderHeightL
22  ! and BorderHeightH are used to detect particles. This is for, say, neutrino
23  ! observation. See BorderHeight{L,H}. However, the primary is directed to
24  ! the deepest depth. \newline
25  ! 1,-1$ \Rightarrow $ detector at the observation place is horizontal. Note
26  ! that the horizontal means not tangential plane, but rather a spherical surface \newline
27  ! 2,-2$ \Rightarrow $ detector is perpendicular to the primary. \newline
28  ! 3$ \Rightarrow $ spherical observation. See text. \newline
29  ! For ObsPlane={1,2}, the user observation routine will receive coordinate values in
30  ! the corresponding detector system. However, if it is 0, 3 or negative, Exyz values
31  ! are obtained.
32  integer NoOfSites2 !2 No of Sites for particle observation; not to be touched; for skeleton/flesh use.
33  real*8 XaxisFromSouth !2 Angle between the horizontal detector X-axis and the south(deg). + is counter
34  ! clockwise. If $|$XaxisFromSouth$| > 360$, it is computed so that the direction is
35  ! to the magnetic east at the deepest observation point. Default is 361.
36 ! <-) --------------------------------------------
37 
38 
39 
40  common /Zobsc/
48 
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos depth
Definition: ZavoidUnionMap.h:1
real(8), parameter, public m
Definition: cpdgXs.f:13
! Zobsvp h parameters to be given by input ! This must be preceded by Zobs * ASHeightList(maxNoOfASSites)
float real
Definition: Zdef.h:2
int integer
Definition: Zdef.h:3
********************block data cblkHeavy ********************integer j data *HeavyG2symbol H
Definition: cblkHeavy.h:7
! Zobs h header file for observation sites definition ! integer maxNoOfSites
Definition: Zobs.h:4
! Zobsvp h parameters to be given by input ! This must be preceded by Zobs * DepthList(maxNoOfSites)
! hadronic collision parameters(-> --------------------------------------------- character *64 IntModel !1 Interaction model description. Usage was changed from v6.0. ! One may list code name and upper energy bound for the code.\newline ! E.g. IntModel='"dpmjet3"' ;to specify the dpmjet3 in the entire energy region ! IntModel='"dpmjet3" 100 "qgsjet2" to specify dpmjet3 at $<$ 100 GeV and qgsjetII ! at E$> $100 GeV \newline ! IntModel='"nucrin" 5 "fritiof1.6" 500 "adhoc" to specify Nucrin, ! fritiof1.6, and ad-hoc model in the respective energy region. This ! corresponds to the old IntModel='int1'. \newline ! IntModel='"nucrin" 5 "fritiof1.6" 10 "fritiof7.02" and \newline ! IntModel='"dpmjet3"' \newline ! are most promissing models that fit the observed ! data(muons and gamma rays) for which the primary is well known by ! BESS and AMS observations($<$ 100 GeV). character *64 XsecModel !1 Xsection model description. Noarmally should not be given. ! Defaul is to use the hadronic xsection given by each active model ! fixed by IntModel. However, for some experimental purposes, ! one may employ x-section by other interaction model. ! e.g. ! IntModel='"phits" 2 "dpmjet" ' ! XsecModel='"Phits" 2 "dpmjet" 80 "qgsjet2"' ! is one example. Default is blank and is replaced by IntModel character *100 InclusiveFile !2 The path to the inclusive data file, xdist.d. Default is ! "../Contrib/Inclusive/xdist.d" real *8 SucPw !2 In the 2nd, 3rd,.. collision of a nucleon inside a nucleus, the collision is ! made to be more elastic than normal one. The leading particle spectrum is ! sampled from x **SucPw dx. SucPw should be in 1 to 2. real *8 Eta2Pi0 !2 eta/pi0 ratio. this is used to see the effect due to non-decay of pi0 ! at very high energies. Only source of h.e gamma can be eta and LPM may work ! for them. default is 0.2 integer MulLow !2 if 1, QCD predicted multiplicity law is used in the adhoc model else UA5 ! parametalization is used. Default is 1.(from v5), ! 0.6135exp(23/18sqrt(2log(roots/0.3))) is QCD jet prediction. ! 7.2roots **0.254 -7 is UA5 data. The branch point is set at roots=900 GeV. !(I have adjusted 0.6135 so that 900 GeV is the b.p) integer LundPara !2 To control Lund program. LundPara(1) is set to kfr(7);kfr(7)=1 is for Frititof ! hard scattering. 2 is for Pythia H.S. 2 gives higher multiplicity but shape is ! strange. Default is 1. LundPara(2) is set to kfr(12):1 by for OPAL hard scattering ! parameterization. 2 by DELPHI. Default is 2.(2 gives bit higher PT). LundPara(3) ! $>$ 0 $\Rightarrow$ Pythia message will appear. LundPara(4) $>$ 0 Fritiof ! message;both on ErrorOut. LundPara(5)=0 $\Rightarrow$ All kaons collisions ! are treated as pi- in Fritiof, else they are treated by adhoc model as they are. ! below:Obsolete integer TotXSopt !2 option for total x-section. 1. PDG 3. TOTEM(COMPETE fitting) ! 2. between 1 and 3.(1 is lowest 3 is highest X-sec.) ! Diff. becomes gradually seen at roots > few x 10^2 GeV ! This for p-p case. Default is 2. ! For other collision types, 1 is used. ! If the current interaction model supplies the inelastic ! cross-section, it is used without referring to this. ! However, see XsecModel. integer SxAbySxpOpt !2 For nucleus target with mass # A, cross section is converted ! from the one for proton target. This option fixes which ! SxA/sxp table is used. This is used when the current interaction ! model dose not supply the cross-section.(See XsecModel too). ! 1. use table derived from QGSJET-II-04(default) ! 2. use table derived from dpjmet3 ! 3. use table dervied from EPOS(but can be used for A<=64). !(as of 2013/Jun). ! However, for small cross-sections(such as gamma-A, nu-A, old ! cxp2xA routine is used). real *8 Cepic0 !2 Obsolete real *8 Cekaon !2 Obsolete integer SucInt !2 The number of successive interactions inside A is affected by this parameter.\newline ! If $0\rightarrow$ old formula(before uv3.500) is used, which give rather ! smaller number($< Nsuc >$ in Air=1.7 for 30 mb pp), \newline ! if $1\rightarrow$ new one $< Nsuc >=2.2$ for 30 mb pp). \newline ! Default is 0(from V5.00 again). real *8 Ceneuc !2 \verb|p -> n
subroutine time(xxx)
Definition: chook.f:5
! This namelist data is frequently used ! Some of them should be given mandatory namelist Param HeightList
Definition: Znamelist1.h:5
! Parameters used for hadronic cascade shower is generated newline ! For you may give as as or em as(order/case/separator insensitive) is to generate EM-cascade and AS. \newline ! Generate
! Zobsvp h parameters to be given by input ! This must be preceded by Zobs LongitOfSite
Definition: Zobsp.h:45
integer Charge2heavyG charge of heavy rightarrow heavy group index conversion array integer HeavyG2massN heavy group index rightarrow mass number conversion array integer HeavyG2charge heavy group index rightarrow charge of heavy conversion array integer HeavyG2code heavy group index rightarrow particle code conversion array integer Code2massN particle code rightarrow mass number conversion array integer Code2heavyG particle code rightarrow heavy group index conversion array real *FragmentTbl the number of interacting nucleons among a projectile heavy nucleus is ! determined as the number of first collision of each interacting nucleon inside ! the nucleus If
Definition: Zheavyp.h:17
! Zobs h header file for observation sites definition ! integer horizontal
Definition: Zobs.h:4
logical, save user
Definition: cEfield0.f:32
#define BorderHeightL
Definition: Ztrackp.h:244
integer maxbin nregion c minsave drx2 ! drawing region real maxdep integer maxpos integer to
Definition: Zfit.h:15
max ptcl codes in the system(except krare) integer kphoton
! Zobsvp h parameters to be given by input ! This must be preceded by Zobs * XaxisFromSouth
Definition: Zobsp.h:45
! Parameters used for hadronic cascade shower is generated newline ! For you may give as as or em quick generation of AS for heavy primaries is tried See chookASbyH f character *Generate2 don t touch this for skeleton flesh use integer MagBrem no magnetic bremsstrahlung is considered newline ! if and Ee energy loss due to magnetic brems is considered newline ! if and Ee real sampling of gamma is performed WaitRatio ! must be made small so that WaitRatio *E0 sim MagBremEmin integer MagPair no magnetic pair creation is considered newline ! if and Eg real sampling is the LPM effect is considered when Ee LpmBremEmin for electrons and ! Eg LpmPairEmin for gamma rays real *MagBremEmin E magnetic bremsstrahlung by electrons may be considered However
Definition: Ztrackp.h:169
#define BorderHeightH
Definition: Ztrackp.h:239
! Zobs h header file for observation sites definition ! integer maxNoOfASSites
Definition: Zobs.h:4
real(4), save a
Definition: cNRLAtmos.f:20
! Zobsvp h parameters to be given by input ! This must be preceded by Zobs DtGMT
Definition: Zobsp.h:45
! Zobsvp h parameters to be given by input ! This must be preceded by Zobs * ObsPlane
Definition: Zobsp.h:45
! Parameters used for hadronic cascade shower is generated newline ! For you may give as as or em quick generation of AS for heavy primaries is tried See chookASbyH f character *Generate2 don t touch this for skeleton flesh use integer MagBrem no magnetic bremsstrahlung is considered newline ! if and Ee energy loss due to magnetic brems is considered newline ! if and Ee real sampling of gamma is performed newline(note, actually upsilon is referred further). ! if generate
! Zobs h header file for observation sites definition ! integer * perpendicular
Definition: Zobs.h:4
integer maxbin nregion c minsave drx2 ! drawing region real maxdep integer maxpos integer from
Definition: Zfit.h:15
********************block data cblkHeavy ********************integer j data *HeavyG2symbol L
Definition: cblkHeavy.h:7
! structure defining a particle at production ! Basic idea of what is to be contained in ! the particle structue is that dynamical ones should be included those derivable from the particle code ! is not included ******************************************************type fmom momentum sequence union map real e endmap map real t endmap map ! pt before pz is set real rap endmap tm endmap endunion end type fmom ! Important say
Definition: Zptcl.h:21
! Zobsvp h parameters to be given by input ! This must be preceded by Zobs * ASDepthList(maxNoOfASSites)
! Zobsvp h parameters to be given by input ! This must be preceded by Zobs YearOfGeomag
Definition: Zobsp.h:45
block data cblkIncident data *Za1ry is
Definition: cblkIncident.h:5
real zen logical degree
Definition: Zflux.h:1
! Zobsvp h parameters to be given by input ! This must be preceded by Zobs * LatitOfSite
Definition: Zobsp.h:45
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const Note
Definition: cblkElemag.h:7
! This namelist data is frequently used ! Some of them should be given mandatory namelist Param Generate
Definition: Znamelist1.h:5
! This namelist data may be given only some ! special purpose such as to see the scattering effect ! by putting NoOfSites2
Definition: Znamelist2.h:3
const int spherical
Definition: Zobs.h:19