Changes between Version 25 and Version 26 of InstallingCorsika


Ignore:
Timestamp:
Oct 20, 2018, 2:33:11 PM (8 months ago)
Author:
tbretz
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • InstallingCorsika

    v25 v26  
    123123> **IACTEXT** The interface to theTELOUTfunction is extended by parameters describing theemitting particle.   This extended information is stored as an additional photon bunch(after the normal one) with mass, charge, energy, and emission time replacing thecx,cy,photons, andzemfields, respectively, and are identified by a wavelength of 9999.The compact output format is disabled for making that possible. In addition, all particlesarriving at the CORSIKA observation level are included in the eventio format output file,in a photon-bunch like block identified by array and detector numbers 999.
    124124
     125== The 700nm problem ==
     126
     127This is were the data is read from the files
     128{{{#!
     129C=======================================================================
     130
     131      SUBROUTINE TPDINI( OBLECE )
     132
     133C-----------------------------------------------------------------------
     134[...]
     135C  READ THE TABLE OF ATMOSPHERIC EXTINCTION AND STORE IT IN THE
     136C  ATMABS MATRIX.
     137[...]
     138        DO  I = 1, 105
     139          READ(MCERABS,21) WLT
     140 21       FORMAT(I4)
     141          READ(MCERABS,*,ERR=995,END=995) (ATMABS(I,J), J=0,50)
     142        ENDDO
     143[...]
     144C  CALCULATE THE ATM. EXTINCTION FOR OBS. LEVEL INTERPOLATING THE TABLE
     145        DO  IWL = 1, 105
     146          FX0 = ATMABS(IWL,X0)
     147          FX1 = ATMABS(IWL,X1)
     148          ATEOBS(IWL) = LINEAR(X,X0,X1,FX0,FX1)
     149        ENDDO
     150      ENDIF
     151[...]
     152C  READ THE QUANTUM EFFECIENCY OF THE PMT
     153        OPEN(UNIT=MCERQEF,FILE=DATDIR(1:INDEX(DATDIR,' ')-1)//
     154     *          'quanteff.dat',STATUS='OLD',FORM='FORMATTED',ERR=996)
     155        READ(MCERQEF,20) TEXTQEF
     156        READ(MCERQEF,23,ERR=997,END=997) (QUAEFF(I),I = 1,105)
     157 23     FORMAT(8F6.3)
     158        CLOSE( MCERQEF )
     159        WRITE(MONIOU,22) TEXTQEF
     160      ENDIF
     161[...]
     162C  READ THE MIRROR REFLECTIVITY
     163        OPEN(UNIT=MCERMIR,FILE=DATDIR(1:INDEX(DATDIR,' ')-1)//
     164     *          'mirreff.dat',STATUS='OLD',FORM='FORMATTED',ERR=998)
     165        READ(MCERMIR,20) TEXTREF
     166        READ(MCERMIR,23,END=999,ERR=999) (MIRREF(I), I=1,105)
     167        CLOSE( MCERMIR )
     168        WRITE(MONIOU,22) TEXTREF
     169      ENDIF
     170}}}
     171
     172
     173This is where the index or the y-value and the corresponding x-values are calculated for mirror and pde efficiency.
     174{{{#!fortran
     175C=====================================================================
     176
     177      SUBROUTINE TELEFF( ABSORB )
     178
     179C-----------------------------------------------------------------------
     180[...]
     181C  CALCULATE THE REFERENCE WL AND INDEX OF WL FOR THE INTERPOLATIONS
     182      RIWL = 1 + INT( (WL-180.D0)/5.D0 )
     183      WLI0 = RIWL*5 + 175
     184      WLI1 = RIWL*5 + 175 + 5
     185}}}
     186
     187This is the atmospheric absorption
     188{{{#!fortran
     189*-- Author : V. de Souza Filho, Uni. Campinas   22/06/1999
     190C=======================================================================
     191
     192      SUBROUTINE ATABSO( ABSORB )
     193
     194C-----------------------------------------------------------------------
     195[...]
     196C  CALCULATE THE REFERENCE WL AND INDEX OF WL FOR THE INTERPOLATIONS
     197      RIWL = 1 + INT( (WL-180.D0)/5.D0 )
     198      WLI0 = RIWL*5 + 175
     199      WLI1 = RIWL*5 + 175 + 5
     200}}}
     201
     202
    125203== Conclusion ==
    126204