2 ! (-> ---------------------------------------------------
4 complex*16
CosZenith !1 Range of
cos(zenith angle). Say, (0.5, 1.0). Used when Za2ry
is 'is' 6 ! the zenith angle range at the incident point (not in Exyz
system).
7 ! In that
case, azimuth
is 0
to 2
pi.
8 complex*16
Azimuth !1 Range of azimuthal angle in deg. Say, (0, 45). Default
is (0,360).
9 ! Can be such
as (300., 390.). Used when Za1ry
is 'is'\newline 13 ! In
this case,
for the upper opening angle, min( Imag(
Azimuth),180.)
is used.
14 character*4 Za1ry !1 Specify the primary angle
sampling method by one of '
is', 'ps' or 'aps'.\
newline 22 !
If this is $<$ deepest obs. level and zeinth angle of primary is $< 0$,
23 ! the primary is assumed
to be upgoing even
if Reverse =0.
26 ! The primary is directed
to this height above the detector.
27 !
If ObsPlane is 3 (spherical), not used.
29 ! <-) ----------------------------------------------------
33 common /Zincidentc/ Za1ry
real(8), parameter, public m
! parameters for primary angle Ddelta
! for length to thickness conversion or v v ! integer maxnodes real Hinf ! This is used when making table for dim simulation ! The slant length for vertical height to km is ! divided by LenStep m steps ! It can cover the slant length of about km for cos
! 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
! This namelist data may be given only some ! special purpose such as to see the scattering effect ! by putting OffsetHeight
! parameters for primary angle CosZenith
! This namelist data may be given only some ! special purpose such as to see the scattering effect ! by putting Reverse
! 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
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
integer maxbin nregion c minsave drx2 ! drawing region real maxdep integer maxpos integer to
max ptcl codes in the system(except krare) integer kphoton
! parameters for primary angle sampling(-> --------------------------------------------------- complex *16 CosZenith !1 Range of cos(zenith angle). Say,(0.5, 1.0). Used when Za2ry is 'is' ! If ObsPlane=3(spherical), real(CosZenith) must be $> $0, and means ! the zenith angle range at the incident point(not in Exyz system). ! In that case, azimuth is 0 to 2pi. complex *16 Azimuth !1 Range of azimuthal angle in deg. Say,(0, 45). Default is(0, 360). ! Can be such as(300., 390.). Used when Za1ry is 'is'\newline ! If ObsPlane=3(spherical), this is used to show the half opening angle ! range where the primary injection position is uniformly distributed ! on a sphere. The center of the opening angle is(Latit, Longit, HeightOfInj). ! In this case, for the upper opening angle, min(Imag(Azimuth), 180.) is used. character *4 Za1ry !1 Specify the primary angle sampling method by one of 'is', 'ps' or 'aps'.\newline ! "is" is isotropic. The range is by CosZenith.\newline ! "ps" is for point source(See also SourceDec)\newline ! "aps" is around point source(See also SourceDec and Ddelta) \newline ! If ObsPlane=3(spherical), this must be "is". real *8 SourceDec !1 Source declination of point source.(deg) real *8 Ddelta !1 SourceDec $\pm$ Ddelta is the region for 'aps'(deg). real *8 HeightOfInj !1 The vertical height of primary injection point(m). ! If this is $<$ deepest obs. level and zeinth angle of primary is $< 0$, ! the primary is assumed to be upgoing even if Reverse=0. ! NOTE:BorderHeightH must be given explicitly in this case. real *8 OffsetHeight !2 The vertical offset height from the deepest detector. ! The primary is directed to this height above the detector. ! If ObsPlane is 3(spherical), not used. !<-) ---------------------------------------------------- common/Zincident/Azimuth
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos height
! constants thru Cosmos real * pi
! rigidity cut related commons *******************************integer maxDim2 parameter maxAzm real Latit
integer lengid integer lengdir character *dir integer kgetenv2 character *numb character *execid character *msg logical takehist save do nsites tkarspec nsites if(histdep(i) .gt. 0) then takehist
! parameters for primary angle SourceDec
! parameters for primary angle * HeightOfInj
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos *Zfirst pos *Zfirst Zfirst vec w *Zfirst vec w Zfirst vec *Zfirst Zfirst where
! 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
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below this
! This namelist data is frequently used ! Some of them should be given mandatory namelist Param Azimuth
integer maxbin nregion c minsave drx2 ! drawing region real maxdep integer maxpos integer from
block data cblkIncident data *Za1ry is
1.8.13