Changeset 19343 for trunk/Mars/mcorsika


Ignore:
Timestamp:
10/30/18 22:12:37 (6 years ago)
Author:
tbretz
Message:
Reading and printing some more options from the file headers
Location:
trunk/Mars/mcorsika
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/Mars/mcorsika/MCorsikaRunHeader.cc

    r19341 r19343  
    5050// ----------------
    5151//  + UInt_t  fCerenkovFileOption
     52//  + UInt_t fHadronModelLowEnergy
     53//  + UInt_t fHadronModelHighEnergy
     54//  + Float_t fTransitionEnergy
     55//  + Bool_t fCurvedAtmosphere
    5256//  + Float_t fEnergyCutoffHadrons
    5357//  + Float_t fEnergyCutoffMuons
     
    194198    fMagneticFieldAz =  g[91];  // Azimuth angle of magnetic north expressed in telescope coordinates
    195199
     200    fHadronModelLowEnergy  = TMath::Nint(g[73]);
     201    fHadronModelHighEnergy = TMath::Nint(g[74]);
     202
    196203    // WITH rounding: unbelievable!
    197204    fCerenkovFlag = TMath::Nint(g[75]);
    198205    fCerenkovFileOption = TMath::Nint(g[90]);
     206
     207    fCurvedAtmosphere = TMath::Nint(g[77])==2;
    199208
    200209    fZdMin = g[79];                // lower edge of theta in °
     
    235244    fViewConeOuterAngle = g[152];  // outer angle of view cone (°)
    236245
     246    fTransitionEnergy = g[153];
     247
    237248    return kTRUE;
    238249}
     
    328339    else
    329340        *fLog << "<off>" << endl;
     341
     342    *fLog << "Interact.model: ";
     343    switch (fHadronModelLowEnergy)
     344    {
     345    case 1:  *fLog << "GEISHA"; break;
     346    case 2:  *fLog << "UrQMD";  break;
     347    case 3:  *fLog << "FLUKA";  break;
     348    default: *fLog << "<n/a>";  break;
     349    }
     350    *fLog << " / ";
     351    switch (fHadronModelHighEnergy)
     352    {
     353    case 0:  *fLog << "HDPM";   break;
     354    case 1:  *fLog << "VENUS";  break;
     355    case 2:  *fLog << "SIBYLL"; break;
     356    case 3:  *fLog << "QGSJET"; break;
     357    case 4:  *fLog << "DPMJET"; break;
     358    case 5:  *fLog << "neXus";  break;
     359    case 6:  *fLog << "EPOS";   break;
     360    default: *fLog << "<n/a>";  break;
     361    }
     362    *fLog << " [lo/hi], Transition at " << fTransitionEnergy << " GeV" << endl;
    330363
    331364    *fLog << "Options used:  ";
     
    343376        *fLog << " VOLUMEDET";
    344377    if (Has(kCurved))
    345         *fLog << " CURVED";
     378        *fLog << " CURVED" << (fCurvedAtmosphere?"<on>":"<off>");
    346379    if (Has(kSlant))
    347380        *fLog << " SLANT";
  • trunk/Mars/mcorsika/MCorsikaRunHeader.h

    r19332 r19343  
    7171    UInt_t fCerenkovFileOption;       // MCERFI
    7272
     73    UInt_t fHadronModelLowEnergy;
     74    UInt_t fHadronModelHighEnergy;
     75
     76    Float_t fTransitionEnergy;        // [GeV]
     77
     78    Bool_t fCurvedAtmosphere;
     79
    7380    Float_t fEnergyCutoffHadrons;     // [GeV]
    7481    Float_t fEnergyCutoffMuons;       // [GeV]
     
    8188    Float_t fThinningWeightLimitEM;    // EM weight limit: WMAX*WEITRAT
    8289    Float_t fThinningMaxRadius;        // [cm] Max radial raius for thinning
     90
    8391
    8492public:
Note: See TracChangeset for help on using the changeset viewer.