Changes between Version 3 and Version 4 of InstallingCorsika

Timestamp:

10/20/18 12:02:19 (6 years ago)

Author:

tbretz

Comment:

--

InstallingCorsika

@@ -57 +57 @@
 
 CEFFIC/CERWLEN: The 8-th number (production height) of each photon is replaced by the wavelength
+
+{{{#!fortran
+C  LIMITING FACTOR FOR STEP SIZE OF ELECTRON IN MAGNETIC FIELD
+#if __CERENKOV__ && __IACT__
+C  LIMIT IN DEFLECTION ANGLE IS 2.5 MILLIRADIAN = 0.143 DEG
+C  WE USE A LIMIT OF ABOUT 0.05 DEG (APPROX. 1 MILLIRAD)
+      BLIMIT   = 0.001D0 / BNORM
+#else
+C  WE USE A LIMIT OF ABOUT 11.4 DEG (0.2 RAD)
+      BLIMIT   = 0.2D0 / BNORM
+#endif
+}}}
+
+{{{#!fortran
+#if __CERENKOV__ && __IACT__
+C  SCATTERING ANGLES OF MUONS SHOULD BE SMALLER THAN THE PIXEL SIZE.
+*     AUX = MIN( 1.D0, 0.015D0*GAMMA )
+C  THE SAME SHOULD HOLD FOR DEFLECTION IN THE GEOMAGNETIC FIELD.
+C  HERE USING A MAXIMUM RMS SCATTERING / DEFLECTION ANGLE OF 0.05 DEG
+C  AND APPROXIMATE ALL BETA*GAMMA TERMS BY GAMMA.
+Cxx      Write(*,*) 'mu step old-style step=',MIN( 1.D0,0.015D0*GAMMA )
+C  FOR A MEAN SCATTERING ANGLE THETA WE HAVE A STEP LENGTH OF ABOUT
+C  (THETA / (13.6 MEV/(BETA*C*P))**2 RADIATION LENGTHS (PDG),
+C  NOT TAKING INTO ACCOUNT THE NON-GAUSSIAN PART OF THE DISTIBUTION.
+C  FOR THE MOMENT DON''T CARE ABOUT THE DIFFERENCE BETWEEN THE
+C  'COULOMB SCATTERING LENGTH' 37.7 G/CM**2 (=C(21)) AND THE
+C  RADIATION LENGTH OF 36.66 OR 36.62 G/CM**2 IN AIR.
+C  NOTE: PI/180/(13.6 MEV/(BETA*C*P)) APPROX 0.136*GAMMA FOR MUONS.
+      AUX = MIN( 1.D0, ((0.05*0.136)*GAMMA)**2 )
+      IF ( BNORMC .GT. 0.D0 ) THEN
+C  NOTE: PI/180*PAMA(5)*BETA*GAMMA APPROX 0.00185*GAMMA
+         AUX = MIN( AUX, (0.05*0.00185)*GAMMA*RHOF(H)/(BNORMC*C(21)) )
+      ENDIF
+Cxx      Write(*,*) 'mu step new-style step=',AUX
+#else
+      AUX = MIN( 10.D0, 0.015D0*GAMMA )
+#endif
+}}}