Changeset 18158


Ignore:
Timestamp:
03/05/15 08:41:58 (10 years ago)
Author:
ftemme
Message:
Added an additional photon acceptance to the simulation.
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Mars/mjobs/MJSimulation.cc

    r18038 r18158  
    244244    write.AddContainer("MRawEvtHeader",       "Events");
    245245    write.AddContainer("MMcEvt",              "Events", kFALSE);
     246    write.AddContainer("MMcEvtBasic",         "Events", kFALSE);
    246247    write.AddContainer("IncidentAngle",       "Events", kFALSE);
    247248    write.AddContainer("MPointingPos",        "Events", kFALSE);
     
    439440    MParSpline splinecones("ConesAngularAcceptance");
    440441    MParSpline splinecones2("ConesTransmission");
     442    MParSpline splineAdditionalPhotonAcceptance("AdditionalPhotonAcceptance");
    441443    plist.AddToList(&splinepde);
    442444    plist.AddToList(&splinemirror);
    443445    plist.AddToList(&splinecones);
    444446    plist.AddToList(&splinecones2);
     447    plist.AddToList(&splineAdditionalPhotonAcceptance);
    445448
    446449    const TString sin2 = "sin(MCorsikaEvtHeader.fZd)*sin(MCorsikaRunHeader.fZdMin*TMath::DegToRad())";
     
    458461    MSimAbsorption cones("SimConesAngularAcceptance");
    459462    MSimAbsorption cones2("SimConesTransmission");
     463    MSimAbsorption additionalPhotonAcceptance("SimAdditionalPhotonAcceptance");
    460464    absapd.SetParName("PhotonDetectionEfficiency");
    461465    absmir.SetParName("MirrorReflectivity");
     
    463467    cones.SetUseTheta();
    464468    cones2.SetParName("ConesTransmission");
     469    additionalPhotonAcceptance.SetParName("AdditionalPhotonAcceptance");
    465470 
    466471    MSimPointingPos pointing;
     
    796801            tasks.AddToList(&absapd);
    797802            tasks.AddToList(&absmir);
     803            tasks.AddToList(&additionalPhotonAcceptance);
    798804            if (0)
    799805                tasks.AddToList(&simatm); // FASTER?
Note: See TracChangeset for help on using the changeset viewer.