Changeset 9425 for trunk/MagicSoft/Mars/msimcamera

Timestamp:

04/16/09 12:04:29 (15 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msimcamera

Files:

• MSimCalibrationSignal.cc (modified) (5 diffs)
• MSimCalibrationSignal.h (modified) (2 diffs)
• MSimCamera.cc (modified) (2 diffs)
• MSimCamera.h (modified) (1 diff)
• MSimGeomCam.cc (modified) (3 diffs)
• MSimGeomCam.h (modified) (2 diffs)
• MSimRandomPhotons.cc (modified) (7 diffs)
• MSimRandomPhotons.h (modified) (3 diffs)

trunk/MagicSoft/Mars/msimcamera/MSimCalibrationSignal.cc

@@ -61 +61 @@
 
 #include "MGeomCam.h"
-#include "MPulseShape.h"
+#include "MParSpline.h"
 #include "MTriggerPattern.h"
 
@@ -136 +136 @@
     }
 
-    fPulse = (MPulseShape*)pList->FindObject("MPulseShape");
+    fPulse = (MParSpline*)pList->FindObject("PulseShape", "MParSpline");
     if (!fPulse)
     {
-        *fLog << err << "MPulsShape not found... aborting." << endl;
+        *fLog << err << "PulsShape [MParSpline] not found... aborting." << endl;
         return kFALSE;
     }
@@ -173 +173 @@
     }
 
+    // FIXME: Is there a way to write them as LAST entry in the file?
     fRunHeader->SetReadyToSave();
 
@@ -195 +196 @@
         for (UInt_t idx=0; idx<fGeom->GetNumPixels(); idx++)
         {
+            // FIXME: Scale number of photons with the pixel size!
             const Int_t num = TMath::Nint(gRandom->Gaus(fNumPhotons, fNumPhotons/10));
 
@@ -229 +231 @@
     // Length (ns), Pulse position (Units ns)
     const Float_t pp = fPulsePos->GetVal();
-    const Float_t pw = fPulse->GetPulseWidth();
+    const Float_t pw = fPulse->GetWidth();
 
     const Float_t first = cnt>0 ? fEvt->GetFirst()->GetTime() : 0;

trunk/MagicSoft/Mars/msimcamera/MSimCalibrationSignal.h

@@ -8 +8 @@
 class MGeomCam;
 class MParList;
-class MPulseShape;
+class MParSpline;
 class MPhotonEvent;
 class MPhotonStatistics;
@@ -20 +20 @@
     MParList          *fParList;    //! Store pointer to MParList for initializing ReInit
     MGeomCam          *fGeom;       //! Camera geometry to know the number of expected pixels
-    MPulseShape       *fPulse;      //! Pulse Shape to get pulse width from
+    MParSpline        *fPulse;      //! Pulse Shape to get pulse width from
     MParameterD       *fPulsePos;   //! Expected position at which the pulse should be
     MParameterD       *fTrigger;    //! Position in analog channels at which the triggersignal  is raised

trunk/MagicSoft/Mars/msimcamera/MSimCamera.cc

@@ -48 +48 @@
 
 #include "MSpline3.h"
-#include "MPulseShape.h"
+#include "MParSpline.h"
 
 #include "MParList.h"
-//#include "MParameters.h"
 
 #include "MPhotonEvent.h"
@@ -134 +133 @@
         return kFALSE;
 
-    MPulseShape *pulse = (MPulseShape*)pList->FindObject("MPulseShape");
+    MParSpline *pulse = (MParSpline*)pList->FindObject("PulseShape", "MParSpline");
     if (!pulse)
     {
-        *fLog << err << "MPulseShape not found... aborting." << endl;
+        *fLog << err << "PulseShape [MParSpline] not found... aborting." << endl;
         return kFALSE;
     }

trunk/MagicSoft/Mars/msimcamera/MSimCamera.h

@@ -28 +28 @@
     MMcEvt            *fMcEvt;           //! For information stored in MMcEvt
 
-    MSpline3 *fSpline;     // Pusle Shape
+    const MSpline3    *fSpline;          // Pusle Shape
 
     Bool_t fBaselineGain;  // Should the gain be applied to baseline and electronic noise?

trunk/MagicSoft/Mars/msimcamera/MSimGeomCam.cc

@@ -60 +60 @@
 
 #include "MParameters.h"
+#include "MParSpline.h"
 #include "MRawRunHeader.h"
-
-#include "MPulseShape.h"
 
 ClassImp(MSimGeomCam);
@@ -106 +105 @@
     }
 
-    fPulse = (MPulseShape*)pList->FindObject("MPulseShape");
+    fPulse = (MParSpline*)pList->FindObject("PulseShape", "MParSpline");
 /*
     if (!fPulse)
@@ -191 +190 @@
 
     // Length (ns), Pulse position (Units ns)
-    const Float_t pp   = fPulsePos ? fPulsePos->GetVal()     : 0;
-    const Float_t pw   = fPulse    ? fPulse->GetPulseWidth() : 0;
+    const Float_t pp   = fPulsePos ? fPulsePos->GetVal() : 0;
+    const Float_t pw   = fPulse    ? fPulse->GetWidth()  : 0;
 
     fStat->SetTimeMedDev(fEvt->GetTimeMedianDev());

trunk/MagicSoft/Mars/msimcamera/MSimGeomCam.h

@@ -13 +13 @@
 class MParameterD;
 class MRawRunHeader;
-class MPulseShape;
+class MParSpline;
 
 class MSimGeomCam : public MTask
@@ -23 +23 @@
     MParameterD       *fPulsePos; //! Intended pulse position in digitization window [ns]
     MRawRunHeader     *fHeader;   //! Length of digitization window
-    MPulseShape       *fPulse;    //! 
+    MParSpline        *fPulse;    //!
 
     TString fNameGeomCam;

trunk/MagicSoft/Mars/msimcamera/MSimRandomPhotons.cc

@@ -43 +43 @@
 //  Output Containers:
 //   MPhotonEvent
+//   AccidentalPhotonRate [MPedestalCam]
 //
 //////////////////////////////////////////////////////////////////////////////
@@ -62 +63 @@
 #include "MPhotonData.h"
 
+#include "MPedestalCam.h"
+#include "MPedestalPix.h"
+
 #include "MCorsikaRunHeader.h"
+
+#include "MSpline3.h"
+#include "MParSpline.h"
+#include "MReflector.h"
 
 ClassImp(MSimRandomPhotons);
@@ -74 +82 @@
 MSimRandomPhotons::MSimRandomPhotons(const char* name, const char *title)
     : fGeom(0), fEvt(0), fStat(0), /*fEvtHeader(0),*/ fRunHeader(0),
-    fSimulateWavelength(kFALSE),  fNameGeomCam("MGeomCam")
+    fRates(0), fSimulateWavelength(kFALSE),  fNameGeomCam("MGeomCam")
 {
     fName  = name  ? name  : "MSimRandomPhotons";
@@ -113 +121 @@
     }
 
+    fRates = (MPedestalCam*)pList->FindCreateObj("MPedestalCam", "AccidentalPhotonRates");
+    if (!fRates)
+        return kFALSE;
+
     /*
     fEvtHeader = (MCorsikaEvtHeader*)pList->FindObject("MCorsikaEvtHeader");
@@ -132 +144 @@
     }
 
+    MReflector *r = (MReflector*)pList->FindObject("Reflector", "MReflector");
+    if (!r)
+    {
+        *fLog << err << "Reflector [MReflector] not found... aborting." << endl;
+        return kFALSE;
+    }
+
+    const MParSpline *s1 = (MParSpline*)pList->FindObject("PhotonDetectionEfficiency", "MParSpline");
+    const MParSpline *s2 = (MParSpline*)pList->FindObject("ConesAngularAcceptance",    "MParSpline");
+    const MParSpline *s3 = (MParSpline*)pList->FindObject("MirrorReflectivity",        "MParSpline");
+
+    const Double_t d2  = fGeom->GetCameraDist()*fGeom->GetCameraDist();
+    const Double_t pde = s1 && s1->GetSpline() ? s1->GetSpline()->Integral() : 1;
+    const Double_t sr  = s2 && s2->GetSpline() ? s2->GetSpline()->IntegralSolidAngle() : 1;
+    const Double_t mir = s3 && s3->GetSpline() ? s3->GetSpline()->Integral() : 1;
+    const Double_t Ar  = r->GetA()/1e4;
+
+    // Conversion factor to convert pixel area to steradians (because it
+    // is a rather small area we can assume it is flat)
+    const Double_t conv = fGeom->GetConvMm2Deg()*TMath::DegToRad();
+
+    // Multiply all relevant efficiencies
+    MParSpline *s4 = (MParSpline*)s1->Clone();
+    s4->Multiply(*s3->GetSpline());
+
+    const Double_t nm = s4 && s4->GetSpline() ? s4->GetSpline()->Integral() : 1;
+
+    delete s4;
+
+    // /100 to convert the pixel area from mm^2 to cm^2
+    fScale =  nm * TMath::Min(Ar, sr*d2) * conv*conv;
+
+    *fLog << inf;
+    *fLog << "Effective cone acceptance:      " << Form("%.2f", sr*d2) << "m^2" << endl;
+    *fLog << "Reflector area:                 " << Form("%.2f", Ar) << "m^2" << endl;
+    *fLog << "Resulting eff. collection area: " << Form("%.1f", TMath::Min(Ar, sr*d2)) << "m^2" << endl;
+    *fLog << "Eff. wavelength band (PDE):     " << Form("%.1f", pde) << "nm" << endl;
+    *fLog << "Eff. wavelength band (Mirror):  " << Form("%.1f", mir) << "nm" << endl;
+    *fLog << "Eff. wavelength band (PDE+MIR): " << Form("%.1f", nm) << "nm" << endl;
+    *fLog << "Pixel area of " << fNameGeomCam << "[0]: " << Form("%.2e", (*fGeom)[0].GetA()*conv*conv) << "sr" << endl;
+    //*fLog << "Effective angular acceptance:     " << sr << "sr" << endl;
+    //*fLog << "Resulting NSB frequency:           " << fFreqNSB*nm*Ar*1000 << "MHz/sr" << endl;
+    *fLog << "Resulting Freq. in " << fNameGeomCam << "[0]: " << Form("%.2f", fFreqNSB*(*fGeom)[0].GetA()*fScale*1000) << "MHz" << endl;
+
+    // const MMcRunHeader *mcrunheader = (MMcRunHeader*)pList->FindObject("MMcRunHeader");
+    // Set NumPheFromDNSB
+
+    // # Number of photons from the diffuse NSB (nphe / ns 0.1*0.1 deg^2 239 m^2) and
+    // nsb_mean 0.20
+    // Magic pixel: 0.00885361 deg
+    // dnsbpix = 0.2*50/15
+    // ampl = MMcFadcHeader->GetAmplitud()
+    // sqrt(pedrms*pedrms + dnsbpix*ampl*ampl/ratio)
+
+    return kTRUE;
+}
+
+Bool_t MSimRandomPhotons::ReInit(MParList *pList)
+{
+    // Overwrite the default set by MGeomApply
+    fRates->Init(*fGeom);
     return kTRUE;
 }
@@ -160 +233 @@
     for (UInt_t idx=0; idx<npix; idx++)
     {
-        // Scale the rate with the pixel size. The rate must
-        // always be given for the pixel with index 0.
-        const Double_t avglen = 1./(fFreqFixed+fFreqNSB*(*fGeom)[idx].GetA()/100.);
+        // Scale the rate with the pixel size.
+        const Double_t rate = fFreqFixed+fFreqNSB*(*fGeom)[idx].GetA()*fScale;
+
+        (*fRates)[idx].SetPedestal(rate);
+
+        // Calculate the average distance between two consequtive photons
+        const Double_t avglen = 1./rate;
 
         // Start producing photons at time "start"
@@ -215 +292 @@
 //    FrequencyNSB:   0.040
 //
-// The frequency is given in units fitting the units of the time.
-// Usually the time is given in nanoseconds thus, e.g., 40 means 40MHz
+// The fixed frequency is given in units fitting the units of the time.
+// Usually the time is given in nanoseconds thus, e.g., 0.040 means 40MHz.
+//
+// The FrequencyNSB is scaled by the area of the pixel in cm^2. Therefore
+// 0.040 would mean 40MHz/cm^2
 //
 Int_t MSimRandomPhotons::ReadEnv(const TEnv &env, TString prefix, Bool_t print)

trunk/MagicSoft/Mars/msimcamera/MSimRandomPhotons.h

@@ -12 +12 @@
 //class MCorsikaEvtHeader;
 class MCorsikaRunHeader;
+class MPedestalCam;
 
 class MSimRandomPhotons : public MTask
@@ -21 +22 @@
 //    MCorsikaEvtHeader *fEvtHeader;  //! Header storing event information
     MCorsikaRunHeader *fRunHeader;  //! Header storing run information
-
+    MPedestalCam      *fRates;   // Random count rate per pixel
 
     // FIXME: Make this a single number per Pixel/APD
     Double_t fFreqFixed; // [1/ns]      A fixed frequency per pixel
     Double_t fFreqNSB;   // [1/ns/cm^2] A frequency depending on area
+
+    Double_t fScale;
 
     Bool_t fSimulateWavelength;
@@ -32 +35 @@
 
     // MTask
-    Int_t PreProcess(MParList *pList);
-    Int_t Process();
+    Int_t  PreProcess(MParList *pList);
+    Bool_t ReInit(MParList *pList);
+    Int_t  Process();
 
     // MParContainer