Changes between Version 3 and Version 4 of InstallingCorsika


Ignore:
Timestamp:
Oct 20, 2018, 12:02:19 PM (5 months ago)
Author:
tbretz
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • InstallingCorsika

    v3 v4  
    5757
    5858CEFFIC/CERWLEN: The 8-th number (production height) of each photon is replaced by the wavelength
     59
     60{{{#!fortran
     61C  LIMITING FACTOR FOR STEP SIZE OF ELECTRON IN MAGNETIC FIELD
     62#if __CERENKOV__ && __IACT__
     63C  LIMIT IN DEFLECTION ANGLE IS 2.5 MILLIRADIAN = 0.143 DEG
     64C  WE USE A LIMIT OF ABOUT 0.05 DEG (APPROX. 1 MILLIRAD)
     65      BLIMIT   = 0.001D0 / BNORM
     66#else
     67C  WE USE A LIMIT OF ABOUT 11.4 DEG (0.2 RAD)
     68      BLIMIT   = 0.2D0 / BNORM
     69#endif
     70}}}
     71
     72{{{#!fortran
     73#if __CERENKOV__ && __IACT__
     74C  SCATTERING ANGLES OF MUONS SHOULD BE SMALLER THAN THE PIXEL SIZE.
     75*     AUX = MIN( 1.D0, 0.015D0*GAMMA )
     76C  THE SAME SHOULD HOLD FOR DEFLECTION IN THE GEOMAGNETIC FIELD.
     77C  HERE USING A MAXIMUM RMS SCATTERING / DEFLECTION ANGLE OF 0.05 DEG
     78C  AND APPROXIMATE ALL BETA*GAMMA TERMS BY GAMMA.
     79Cxx      Write(*,*) 'mu step old-style step=',MIN( 1.D0,0.015D0*GAMMA )
     80C  FOR A MEAN SCATTERING ANGLE THETA WE HAVE A STEP LENGTH OF ABOUT
     81C  (THETA / (13.6 MEV/(BETA*C*P))**2 RADIATION LENGTHS (PDG),
     82C  NOT TAKING INTO ACCOUNT THE NON-GAUSSIAN PART OF THE DISTIBUTION.
     83C  FOR THE MOMENT DON''T CARE ABOUT THE DIFFERENCE BETWEEN THE
     84C  'COULOMB SCATTERING LENGTH' 37.7 G/CM**2 (=C(21)) AND THE
     85C  RADIATION LENGTH OF 36.66 OR 36.62 G/CM**2 IN AIR.
     86C  NOTE: PI/180/(13.6 MEV/(BETA*C*P)) APPROX 0.136*GAMMA FOR MUONS.
     87      AUX = MIN( 1.D0, ((0.05*0.136)*GAMMA)**2 )
     88      IF ( BNORMC .GT. 0.D0 ) THEN
     89C  NOTE: PI/180*PAMA(5)*BETA*GAMMA APPROX 0.00185*GAMMA
     90         AUX = MIN( AUX, (0.05*0.00185)*GAMMA*RHOF(H)/(BNORMC*C(21)) )
     91      ENDIF
     92Cxx      Write(*,*) 'mu step new-style step=',AUX
     93#else
     94      AUX = MIN( 10.D0, 0.015D0*GAMMA )
     95#endif
     96}}}