Changeset 19551


Ignore:
Timestamp:
Jun 5, 2019, 1:23:38 PM (5 months ago)
Author:
tbretz
Message:
When the bunch size is 1 then the photons still carry their original probability (that is my interpretation) which is stored as 'bunch size'. Therefore, all bunch sizes <= 1 are ok.
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Mars/msim/MPhotonData.cc

    r19533 r19551  
    260260        return kCONTINUE;
    261261
    262     if (f[0]!=1)
    263         cout << "MPhotonData::FillCorsika: WARNING - Bunch size !=1 (" << f[0] << ")" <<endl;
     262    // My understanding is that each photon internally gets a weight (e.g. lambda^-2)
     263    // and according to this weight a dice is thrown. The weights are still written
     264    // to the output file but only for the surviving photons
     265    if (f[0]>1)
     266        cout << "MPhotonData::FillCorsika: WARNING - Bunch size > 1 (" << f[0] << ")" <<endl;
    264267
    265268#ifdef __MMCS__
     
    308311Int_t MPhotonData::FillCorsikaThin(const Float_t f[8], Int_t i)
    309312{
     313    // DATAB2(LHCER(IBUF)+1,IBUF) = PHOTCM*WTCER/MAX(1.D-10,PROBTH)
     314
     315    // For the THIN option the photon bunch size is multiplied with the
     316    // thinning weight. See Sect. 4.89 page 99.
     317
     318    /*
     319     #if __THIN__
     320     C      (CONTAINING UP TO 39 BUNCHES, 8 WORDS EACH)
     321     C
     322     C    8*(N-1)+1  NUMBER OF PHOTONS IN BUNCH
     323     C                  (FOR STANDARD PARTICLE OUTPUT FILE:
     324     C                   99.E5 + NINT(NUMBER OF PHOTONS IN BUNCH)*10 + 1
     325     C    8*(N-1)+2  X- COORDINATE IN CM
     326     C    8*(N-1)+3  Y- COORDINATE IN CM
     327     C    8*(N-1)+4  DIRECTION COSINUS TO X AXIS
     328     C    8*(N-1)+5  DIRECTION COSINUS TO Y AXIS
     329     C    8*(N-1)+6  T TIME SINCE FIRST INTERACTION (OR ENTRANCE INTO
     330     C               ATMOSPHERE) IN NSEC
     331     C    8*(N-1)+7  PRODUCTION HEIGHT OF BUNCH IN CM
     332     C    8*(N-1)+8  WEIGHT OF BUNCH
     333     #else
     334     */
     335
    310336    // From the Corsika manual:
    311337    //
     
    346372        return kCONTINUE;
    347373
    348     if (f[0]!=1)
    349         cout << "MPhotonData::FillCorsikaThin: WARNING - Bunch size !=1 (" << f[0] << ")" <<endl;
    350 
     374    // My understanding is that each photon internally gets a weight (e.g. lambda^-2)
     375    // and according to this weight a dice is thrown. The weights are still written
     376    // to the output file but only for the surviving photons
     377    if (f[0]>1)
     378        cout << "MPhotonData::FillCorsikaThin: WARNING - Bunch size > 1 (" << f[0] << ")" <<endl;
    351379
    352380    // x=north, y=west
     
    365393
    366394    fProductionHeight = f[6]; // [cm]
     395    // f[7]<0: Photoelectron bunches of specific wavelength (if __CEFFIC__)
     396    // f[7]>0: __CERWLEN__ (but __CEFFIC takes priority)
    367397
    368398    fWavelength = TMath::Abs(f[7]);
Note: See TracChangeset for help on using the changeset viewer.