Changeset 9991 for trunk/Mars

Timestamp:

09/30/10 15:34:56 (14 years ago)

Author:

tbretz

Message:

Changed the checks in MSimRandomPhotons which didn't allow pedestal and calibration runs to be produced

Location:

trunk/Mars

Files:

• Changelog (modified) (1 diff)
• msimcamera/MSimRandomPhotons.cc (modified) (7 diffs)
• msimcamera/MSimRandomPhotons.h (modified) (1 diff)

trunk/Mars/Changelog

@@ -35 +35 @@
    * msim/MPhotonEvent.cc:
      - don't loop over all (mostly empty) slots in Print()
+
+   * msimcamera/MSimRandomPhotons.[h,cc]:
+     - changed again the checks to also allow pedestal and calibration
+       runs to be produced without having a MCorsikaRunHeader available
 
 

trunk/Mars/msimcamera/MSimRandomPhotons.cc

@@ -150 +150 @@
 }
 
-Bool_t MSimRandomPhotons::CheckWavelengthRange(const MParSpline &sp, const char *txt) const
-{
-    const Float_t min = sp.GetXmin();
-    const Float_t max = sp.GetXmax();
-
-    if (min>fRunHeader->GetWavelengthMin())
-    {
-        *fLog << err << "ERROR - Minimum wavelength (" << min << "nm)";
-        *fLog << " defined for " << txt;
-        *fLog << " exceeds minimum wavelength simulated (";
-        *fLog << fRunHeader->GetWavelengthMin() << "nm)." << endl;
-        return kFALSE;
-    }
-    if (max<fRunHeader->GetWavelengthMax())
-    {
-        *fLog << err << "ERROR - Maximum wavelength (" << max << "nm)";
-        *fLog << " defined for " << txt;
-        *fLog << " undershoots maximum wavelength simulated (";
-        *fLog << fRunHeader->GetWavelengthMax() << "nm)." << endl;
-        return kFALSE;
-    }
-
-    return kTRUE;
-}
-
 // --------------------------------------------------------------------------
 //
@@ -221 +196 @@
 
     fRunHeader = (MCorsikaRunHeader*)pList->FindObject("MCorsikaRunHeader");
-    if (!fRunHeader)
+    if (fSimulateWavelength && !fRunHeader)
     {
         *fLog << err << "MCorsikaRunHeader not found... aborting." << endl;
@@ -239 +214 @@
     const MParSpline *s4 = (MParSpline*)pList->FindObject("ConesAngularAcceptance",    "MParSpline");
 
-    // Check if all splines are defined in the relevant range
-    if (!CheckWavelengthRange(*s1, "PhotonDetectionEfficiency [MParSpline]"))
-        return kFALSE;
-    if (!CheckWavelengthRange(*s2, "ConesTransmission [MParSpline]"))
-        return kFALSE;
-    if (!CheckWavelengthRange(*s3, "MirrorReflectivity [MParSpline]"))
-        return kFALSE;
-
-    const Float_t wmin = fRunHeader->GetWavelengthMin();
-    const Float_t wmax = fRunHeader->GetWavelengthMax();
+    // Ensure that all efficiencies are at least defined in the raneg of the
+    // photon detection efficiency. We assume that this is the limiting factor
+    // and has to be zero at both ends.
+    if (s2->GetXmin()>s1->GetXmin())
+    {
+        *fLog << err << "ERROR - ConeTransmission range must be defined down to " << s1->GetXmin() << "nm (PhotonDetectionEffciency)." << endl;
+        return kFALSE;
+    }
+    if (s2->GetXmax()<s1->GetXmax())
+    {
+        *fLog << err << "ERROR - ConeTransmission range must be defined up to " << s1->GetXmax() << "nm (PhotonDetectionEffciency)." << endl;
+        return kFALSE;
+    }
+    if (s3->GetXmin()>s1->GetXmin())
+    {
+        *fLog << err << "ERROR - MirrorReflectivity range must be defined down to " << s1->GetXmin() << "nm (PhotonDetectionEffciency)." << endl;
+        return kFALSE;
+    }
+    if (s3->GetXmax()<s1->GetXmax())
+    {
+        *fLog << err << "ERROR - MirrorReflectivity range must be defined up to " << s1->GetXmax() << "nm (PhotonDetectionEffciency)." << endl;
+        return kFALSE;
+    }
+
+    // If the simulated wavelenth range exists and is smaller reduce the
+    // range to it. Later it is checked that at both edges the transmission
+    // is 0. This must be true in both cases: The simulated wavelength range
+    // exceed the PDE or the PDE range exceeds the simulated waveband.
+    const Float_t wmin = fRunHeader && fRunHeader->GetWavelengthMin()>s1->GetXmin() ? fRunHeader->GetWavelengthMin() : s1->GetXmin();
+    const Float_t wmax = fRunHeader && fRunHeader->GetWavelengthMax()<s1->GetXmax() ? fRunHeader->GetWavelengthMax() : s1->GetXmax();
 
     const Int_t min = TMath::FloorNint(wmin);
@@ -316 +311 @@
         return kFALSE;
 
-    if (!CheckWavelengthRange(flux, TString("night sky background flux from ")+fFileNameNSB))
-        return kFALSE;
+    if (flux.GetXmin()>wmin)
+    {
+        *fLog << err << "ERROR - NSB flux from " << fFileNameNSB << " must be defined down to " << wmin << "nm." << endl;
+        return kFALSE;
+    }
+    if (flux.GetXmax()<wmax)
+    {
+        *fLog << err << "ERROR - NSB flux from " << fFileNameNSB << " must be defined up to " << wmax << "nm." << endl;
+        return kFALSE;
+    }
 
     MParSpline nsb;
@@ -334 +337 @@
     if (eff.GetSpline()->Eval(min)>1e-5)
     {
-        *fLog << warn << "WARNING - Total transmission efficiency at detector at ";
+        *fLog << warn << "WARNING - Total transmission efficiency at ";
         *fLog << min << "nm is not zero, but " << nsb.GetSpline()->Eval(min) << "... abort." << endl;
     }
     if (eff.GetSpline()->Eval(max)>1e-5)
     {
-        *fLog << warn << "WARNING - Total transmission efficiency at detector at ";
+        *fLog << warn << "WARNING - Total transmission efficiency at ";
         *fLog << max << "nm is not zero, but " << nsb.GetSpline()->Eval(max) << "... abort." << endl;
     }
@@ -346 +349 @@
     if (eff.GetSpline()->Eval(wmin)>1e-5)
     {
-        *fLog << err << "ERROR - Total transmission efficiency at detector at minimum simulated wavelength ";
+        *fLog << err << "ERROR - Total transmission efficiency at ";
         *fLog << wmin << "nm is not zero... abort." << endl;
         *fLog << "        PhotonDetectionEfficency: " << s1->GetSpline()->Eval(wmin) << endl;
@@ -356 +359 @@
     if (eff.GetSpline()->Eval(wmax)>1e-5)
     {
-        *fLog << err << "ERROR - Total transmission efficiency at detector at maximum simulated wavelength ";
+        *fLog << err << "ERROR - Total transmission efficiency at ";
         *fLog << wmax << "nm is not zero... abort." << endl;
         *fLog << "        PhotonDetectionEfficency: " << s1->GetSpline()->Eval(wmax) << endl;

trunk/Mars/msimcamera/MSimRandomPhotons.h

@@ -43 +43 @@
     Int_t ReadEnv(const TEnv &env, TString prefix, Bool_t print=kFALSE);
 
-    // MSimRandomPhotons
-    Bool_t CheckWavelengthRange(const MParSpline &sp, const char *txt) const;
-
 public:
     MSimRandomPhotons(const char *name=NULL, const char *title=NULL);