Changeset 18548 for trunk/Mars/msim/MSimAtmosphere.cc

Timestamp:

08/26/16 13:01:43 (8 years ago)

Author:

tbretz

Message:

If the CEFFIC option is enables and the user has not requested execution explicitly, it is skipped.

File:

• trunk/Mars/msim/MSimAtmosphere.cc (modified) (6 diffs)

trunk/Mars/msim/MSimAtmosphere.cc

@@ -698 +698 @@
 //
 MSimAtmosphere::MSimAtmosphere(const char* name, const char *title)
-    : fEvt(0), fAtmosphere(0),
+    : fRunHeader(0), fEvt(0), fAtmosphere(0),
     fFileAerosols("resmc/atmosphere-aerosols.txt"),
-    fFileOzone("resmc/atmosphere-ozone.txt")
+    fFileOzone("resmc/atmosphere-ozone.txt"),
+    fForce(kFALSE)
 {
     fName  = name  ? name  : "MSimAtmosphere";
@@ -739 +740 @@
 Bool_t MSimAtmosphere::ReInit(MParList *pList)
 {
-    MCorsikaRunHeader *h = (MCorsikaRunHeader*)pList->FindObject("MCorsikaRunHeader");
-    if (!h)
+    fRunHeader = (MCorsikaRunHeader*)pList->FindObject("MCorsikaRunHeader");
+    if (!fRunHeader)
     {
         *fLog << err << "MCorsikaRunHeader not found... aborting." << endl;
         return kFALSE;
+    }
+
+    if (fRunHeader->Has(MCorsikaRunHeader::kCeffic) && !fForce)
+    {
+        *fLog << inf << "CEFFIC enabled... task will be skipped." << endl;
+        return kTRUE;
     }
 
@@ -759 +766 @@
     */
 
-    fAtmosphere->Init(*h, fFileOzone, fFileAerosols);
+    fAtmosphere->Init(*fRunHeader, fFileOzone, fFileAerosols);
 
     if (!fAtmosphere->IsAllValid())
@@ -767 +774 @@
     }
 
-    if (h->GetWavelengthMin()<fAtmosphere->GetWavelengthMin())
+    if (fRunHeader->GetWavelengthMin()<fAtmosphere->GetWavelengthMin())
         *fLog << warn << "WARNING - Lower bound of wavelength bandwidth exceeds valid range of atmosphere." << endl;
 
-    if (h->GetWavelengthMax()>fAtmosphere->GetWavelengthMax())
+    if (fRunHeader->GetWavelengthMax()>fAtmosphere->GetWavelengthMax())
         *fLog << warn << "WARNING - Lower bound of wavelength bandwidth exceeds  valid range of atmosphere." << endl;
 
-    if (!h->Has(MCorsikaRunHeader::kAtmext))
+    if (!fRunHeader->Has(MCorsikaRunHeader::kAtmext))
         *fLog << warn << "WARNING - ATMEXT option not used for Corsika data." << endl;
 
-    if (!h->Has(MCorsikaRunHeader::kRefraction))
+    if (!fRunHeader->Has(MCorsikaRunHeader::kRefraction))
         *fLog << warn << "WARNING - Refraction calculation disabled for Corsika data." << endl;
 
@@ -786 +793 @@
 Int_t MSimAtmosphere::Process()
 {
+    // Skip the task if the CEFFIC option has been enabled and
+    // its excution has not been forced by the user
+    if (!fForce && fRunHeader->Has(MCorsikaRunHeader::kCeffic))
+        return kTRUE;
+
     // Get the number of photons in the list
     const Int_t num = fEvt->GetNumPhotons();
@@ -872 +884 @@
     }
 
+    if (IsEnvDefined(env, prefix, "Force", print))
+    {
+        rc = kTRUE;
+        fForce = GetEnvValue(env, prefix, "Force", fForce);
+    }
+
     return rc;
 }