The description of parameters rarely used:  &HParam data
----------------------------------------------------------

<<AnihiE>>   real*8 
 if E(positron) < Eanihi, annihilation is considered. 

<<BackAngLimit>>   real*8 
 If the cosine of the angle between a particle and the primary becomes smaller than 
 this value, the particle is discarded. See also BorderHeighH. If you give a value 
 less than -1.0, such rejection will never happen. Default is 0. 

<<BaseErg>>   real*8 
 See BaseTime. The default is 1000 (GeV). 

<<BasePower>>   real*8 
 See BaseTime. Default is 1.0 

<<BorderHeightH>>   real*8 
 If a particle goes higher than this, discard it. This should be larger than HeightOfInj or 0. 
 If 0, it is adjusted to be the same as HeightOfInj. NOTE: For upgoin primary cases, you have 
 to set this one explicitly. 

<<BorderHeightL>>   real*8 
 If a particle reaches this hight, call observation routine. No further tracking is done. 
 This is for neutrino observation. See ObsPlane. 

<<Cekaon>>   real*8 
 Obsolute 

<<Ceneuc>>   real*8 
 p -> n ; n-> p; p~->n~; n~->p~ prob. 

<<Cepic0>>   real*8 
 Obsolute 

<<Charge2heavyG>>   integer 
 charge of heavy --> heavy group index conversion array. 

<<Code2heavyG>>   integer 
 particle code --> heavy group index conversion array. 

<<Code2massN>>   integer 
 particle code --> mass number conversion array. 

<<Deltkp>>   real*8 
 k_p pi_p xsection increases as E^Deltpik(E>100GeV) 

<<Deltpip>>   real*8 
 pi_p xsection increases as E^Deltpip(E>100GeV) 

<<Deltpp>>   real*8 
 p_p xsection increases as E^Deltpp(E>100GeV) 

<<Ecrit>>   real*8 
 critical energy in GeV 

<<Efermi>>   real*8 
 If Kinetic E < Efermi, Fermi Momentum is considered for Nucleus target. 

<<Elund2>>   real*8 
 Elund < E < Elund2 ==> ad hoc model is used; Elund2 must be > ~100, if Elund2 < Elund3. 
 Defualt is Elund 

<<Elund3>>   real*8 
 Elund2< E < Elund3 ==> New Lund Fritiof is used. Elund3 can be < Elund2 
 If E > Elund3, ad hoc model is used. Default is Elund. 

<<Elund>>   real*8 
 E < Elund==> Lund M.C is used. in h-n. (IntModel='int1'). If IntModel='int2' 
 ==>Gheisha is used. Elund can be as small as 4.99. (Default is 500. 
 For NEXT 4.99). 

<<EndLevel2>>   integer 
 Don't worry. This is system use. 

<<EndLevel>>   integer 
 Used for skeleton/flesh-out job. In a normal job, system default value 0 is reset by 
 the system to be the max number of observation levels. (=NoOfSites). Its real use is in such a 
 skeleton/flesh-out job that one first follows the particles up to some high depth and later chooses 
 events and flesh them out to deeper depths. In such a skeleton-making job, the user must give the 
 depth list which is used flesh-out job, too. In the skeleton job, particle tracking is terminated 
 at the level specified by EndLevel. In such a flesh-out job, the user must give a larger value 
 or 0 to EndLevel 

<<ErrorOut>>   integer 
 Error output logical dev number. 

<<Es>>   real*8 
 Modified scattering constant. 20.e-3 GeV 

<<Eta2Pi0>>   real*8 
 eta/pi0 ratio. this is used to see the effect due to non-decay of pi0 
 at very high energies. Only soruce of h.e gamma can be eta and LPM may work 
 for them. default is 0.2 

<<EthinRatio>>   real*8 
 if ThinsSamplig 
= 0, thin sampling is performed if the energy of a particle is 
 < EthinRatio * PrimaryEnergy(/nucleon) (=Ethin) 

<<EventNo>>   integer 
 cumulative event number counter. 

<<EventsInTheRun>>   integer 
 Counter for event number in the run. Internal use. 

<<ExactThick>>   logical 
 If T, a given length is converted into thickness with best accuracy even for very 
 inclined trajectory by using numerical integration. 

<<FragmentTbl>>   real*8 
 tbl(i,j)=<Number> of frag. j when a heavy of heavy group index i breaks up at air. 

<<Generate2>>   character*16 
 don't touch this. for skeleton/flesh use. 

<<HeavyG2charge>>   integer 
 heavy group index --> charge of heavy conversion array. 

<<HeavyG2code>>   integer 
 heavy group index --> particle code conversion array. 

<<HeavyG2massN>>   integer 
 heavy group index --> mass number conversion array. 

<<HeavyG2symbol>>   character*4 
 heavy group index --> 'Fe' etc conversion array. 

<<HowGeomag>>   integer 
 if 1, no magnetic field until first coll. 
 2, mag.f always exists. If Reverse not=0, use this. 
 11, same as 1 but mag.f is const. 
 12, same as 2 but mag.f is const. 
 21, same as 1 but mag.f is const. 
 22, same as 2 but mag.f is const. 
 const value is the one at deepest observation plane. for 11,12 or should be given by 
 MagN, MagE, MagD for 21, 22. For normal applications, 11 is good. 

<<HowIntNuc>>   integer 
 If 0, 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 1, the number is determined as the total number of collisions 
 including successive interactions. Default is 1. (There is uncertaninity in 
 interpretation of the formula; value 1 gives larger number of interacting nucleons.) 

<<IncreaseXsec>>   real*8 
 how the xsection increases. 1.0--> power of E. 

<<KEminObs2>>   real*8 
 Don't touch this. skeleton/flesh use. 

<<Knockon>>   logical 
 if T, Knockon process is considered 

<<Kpicns>>   real*8 
 See Kpilog. 

<<Kpilog>>   real*8 
 k-ch/pi-ch=(Kpilog*log(ss+.069)+Kpicns)*exp(-8/s') where ss(GeV**2)= effective s. 
 s'(GeV**2)=s - 4.63. See also Kpicns. 

<<LpmBremEmin>>   real*8 
 The LPM effect is taken into account for bremsstrahlung when LpmEffect is .true. 
 and the electron energy is higher than this. 

<<LpmPairEmin>>   real*8 
 The LPM effect is taken into account for pair creation when LpmEffect is .true. 
 and the gamma energy is higher than this. 

<<LundPara>>   integer 
 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 =>Pythia message will appear. LundPara(4) >0 Fritiof message; both on ErrorOut. 
 LundPara(5) =0 ==>All kaons collisions are treated as pi- in Fritiof, else they 
 are treated by adhoc model as they are. 

<<MagBrem>>   integer 
 If 0, no magnetic bremsstrahlung is considered. 
 if 1 and Ee > MagBremEmin, energy loss due to magnetic brems is considered 
 if 2 and Ee > MagBremEmin, real sampling of gamma is performed. 
 (note, actually upsilon is referred further). 
 if generate='as' with really high energy primaries, WaitRatio 
 must be made small so that WaitRatio*E0 ~MagBremEmin 

<<MagBremEmin>>   real*8 
 E > this, magnetic bremsstrahlung by electrons may be considered. However, if 
 MagBrem = 0, not considered at all 
 MagBrem = 1, total energy loss due to brems is considered. 
 MagBrem = 2, gamma energy is sampled actually. 
 If upsilon (Ee/m * B/Bcr) is small, the effective treatment will be 
 the same as MagBrem = 0 case. 

<<MagChgDist>>   real*8 
 Distance where mag. can be seen as const.(m) at sea level 

<<MagD>>   real*8 
 See HowGeomag (in Tesla) 

<<MagE>>   real*8 
 See HowGeomag (in Tesla) 

<<MagN>>   real*8 
 See HowGeomag (in Tesla) 

<<MagPair>>   integer 
 If 0, no magnetic pair creation is considered. 
 if 1 and Eg > MagPairEmin, real sampling is tried. 
 (note, actually upsilon is referred further). To see these magnetic effects, 
 HowGeoMag=2 and HightOfInj ~ 5000 km are desirable. 

<<MagPairEmin>>   real*8 
 E > this, magnetic pair creation by gamma may be considered. However, if 
 MagPair = 0, not considered at all. 
 MagPair = 1, pair creation is sampled. 
 However, again, actual occurrence will be dependent on the angle between 
 B and photon direction. 

<<MaxComptonE>>   real*8 
 Above this energy, Compton is neglected. 

<<MaxPhotoE>>   real*8 
 Above this energy, photoelectric effect is neglected. 

<<Moliere>>   integer 
 0--> use all Gaussian approx (with air density change and 
 energy loss effect) 
 1--> use Moliere scattering for non-electrons (default) 
 2--> use Moliere scattering for all charged particles. 
 If negative, angle correlated displacement is made to be 0 
 since Moliere theory cannot give it. (if >0, we use Gaussian 
 approximation for correlation). 

<<Mudirp>>   real*8 
 DD~ # enhancement factor. D is only for prompt muon. 

<<MulLow>>   integer 
 if 1, QCD predicted multiplicity low is used in the adhoc model else UA5 
 parametalization is used. Default is 0. 
 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) 

<<NoOfSites2>>   integer 
 No of Sites for particle observation; not to be touched; for skeleton/flesh use. 

<<OffsetHeight>>   real*8 
 The vertical offset height from the deepest detector. 
 The primary is directed to this height above the detector. 

<<PathLimit>>   real*8 
 If the sum of (path/beta) of a particle exceeds this, it is judged as dead. 
 (to avoid infinite cyclotron loop). However, for normal applications, 
 this will not be effective because of BackAnglLimit. See Reverse. 
 TimeStructure should be T if Reverse 
= 0 and PathLimit is to be effective. 

<<PrevEventNo>>   integer 
 The event number already finished. System use for Cont job. 

<<PtAvFrag>>   real*8 
 <Pt> of heavy fragments. 

<<PtAvNonInteNuc>>   real*8 
 <Pt> of non interacting nucleons. 

<<RatioToE0>>   real*8 
 In the A.S generation, hadronic interactions are followed down to at 
 least RatioToE0 * E0/nucleon energy. 

<<RecoilKineMinE>>   real*8 
 Recoil Kinetic Min Energy above which the recoil (=knock-on process) 
 is treated. Below this energy, the effect is included as continuous 
 energy loss. 

<<Reverse>>   integer 
 0=> Normal tracking. 
 1=> incident is tracked to a direction opposite to the given one. 
 the incident is charge-conjugated. 
 All interactions are ignored. (Use when to make cut-off table or to see 
 a given particle (say, observed anti proton) can go out of Earth. 
 2=> same as 1 but energy gain (not loss) is taken into account 
 TimeStructure should be T if Reverse 
= 0. See BackAnglLimit. 

<<SeedFileDev>>   integer 
 logical device number of SeedFile. 

<<SucInt>>   integer 
 The number of successive interactions inside A is affected by this parameter. 
 If 0--> old one (before uv3.500) is used, which give rather smaller number 
 (<Nsuc> in Air = 1.7 for 30 mb pp), if 1--> new one <Nsuc>=2.2 for 30 mb pp). 
 Default is 1. 

<<SucPw>>   real*8 
 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. 

<<TempDev>>   integer 
 Logical Dev. number for temporary disk use. 

<<TraceDev>>   integer 
 Logical dev # for TraceDir/trace1,2,.... 

<<Truncc>>   real*8 
 coeff. for truncating path. 

<<Truncn>>   real*8 
 coeff. for truncating path. 

<<Truncx>>   real*8 
 coeff. for truncating path. 

<<UpsilonMin>>   real*8 
 Magnetic bremsstralhung is considered only if upsilon > UpsilonMin. 

<<X0>>   real*8 
 radiation length in kg/m2 for air 

<<XaxisFromSouth>>   real*8 
 Angle between the horizontal detector X-axis and the south(deg). + is counter 
 clockwise. If |XaxisFromSouth| > 360, it is computed so that the direction is 
 to the magnetic east at the deepest observation point. Default is 361. 

