Changeset 2440 for trunk/MagicSoft/Mars/manalysis

Timestamp:

10/29/03 13:35:42 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/manalysis

Files:

• AnalysisLinkDef.h (modified) (2 diffs)
• MMcPedestalCopy.cc (modified) (7 diffs)
• MMcPedestalCopy.h (modified) (2 diffs)
• MMcPedestalNSBAdd.cc (modified) (5 diffs)
• Makefile (modified) (2 diffs)

trunk/MagicSoft/Mars/manalysis/AnalysisLinkDef.h

@@ -8 +8 @@
 #pragma link C++ class MCerPhotEvt+;
 #pragma link C++ class MCerPhotAnal+;
+#pragma link C++ class MCerPhotAnal2+;
 #pragma link C++ class MCerPhotCalc+;
+
+//#pragma link C++ class MRelSignal+;
 
 #pragma link C++ class MBlindPixels+;
@@ -57 +60 @@
 #pragma link C++ class MMinuitInterface+;
 #pragma link C++ class MFiltercutsCalc+;
-#pragma link C++ class MMcPedestalRead+;
-
 
 #endif
-
-
-
-
-
-
-
-
-
-
-
-

trunk/MagicSoft/Mars/manalysis/MMcPedestalCopy.cc

@@ -16 +16 @@
 !
 !
-!   Author(s): Thomas Bretz  12/2000 <mailto:tbretz@uni-sw.gwdg.de>
+!   Author(s): Thomas Bretz, 12/2000 <mailto:tbretz@astro.uni-wuerzburg.de>
 !
-!   Copyright: MAGIC Software Development, 2000-2001
+!   Copyright: MAGIC Software Development, 2000-2003
 !
 !
@@ -24 +24 @@
 
 /////////////////////////////////////////////////////////////////////////////
-//                                                                         //
-//  MMcPedestalCopy                                                        //
-//                                                                         //
-//  This task looks for the ìnformation about FADC pedestals in            //
-//  MMcFadcHeader and translates it to the pedestal values in              //
-//  MPedestalCam                                                           //
-//                                                                         //
-//  Input Containers:                                                      //
-//   MMcFadcHeader                                                         //
-//                                                                         //
-//  Output Containers:                                                     //
-//   MPedestalCam                                                          //
-//                                                                         //
+//
+//  MMcPedestalCopy
+//
+//  This task looks for the ìnformation about FADC pedestals in
+//  MMcFadcHeader and translates it to the pedestal values in
+//  MPedestalCam
+//
+//  Input Containers:
+//   MMcFadcHeader
+//
+//  Output Containers:
+//   MPedestalCam
+//
 /////////////////////////////////////////////////////////////////////////////
-
 #include "MMcPedestalCopy.h"
 
@@ -50 +49 @@
 
 #include "MRawRunHeader.h"
+#include "MMcRunHeader.hxx"
 #include "MMcFadcHeader.hxx"
 
@@ -89 +89 @@
 // --------------------------------------------------------------------------
 //
-// Check runtype and search for MPedestalCam and MMcFadcHeader.
-// If the runtype check fails the task is removed from the task list.
+// Search for MPedestalCam and MMcFadcHeader.
 //
 Int_t MMcPedestalCopy::PreProcess(MParList *pList)
 {
-    if (!CheckRunType(pList))
-    {
-        *fLog << warn << dbginf << " MMcPedestalCopy is for Monte Carlo files only... ";
-        *fLog << "removing task from list." << endl;
-        return kSKIP;
-    }
+    fMcRun = (MMcRunHeader*)pList->FindObject("MMcRunHeader");
+    if (!fMcRun)
+        *fLog << warn << dbginf << "MMcRunHeader not found... assuming camera<0.7" << endl;
 
     fPedCam = (MPedestalCam*)pList->FindCreateObj("MPedestalCam");
@@ -108 +104 @@
     if (!fMcPed)
     {
-        *fLog << warn << dbginf << "MMcFadcHeader not found... aborting." << endl;
+        *fLog << err << dbginf << "MMcFadcHeader not found... aborting." << endl;
         return kFALSE;
     }
@@ -128 +124 @@
     const int num = fPedCam->GetSize();
 
+    const Bool_t camver70 = fMcRun && fMcRun->GetCamVersion()>=70;
+
     for (int i=0; i<num; i++)
     {
@@ -136 +134 @@
 
         const Float_t pedest = fMcPed->GetPedestal(i);
-        //const Float_t pedrms = pedest/sqrt((float)num);
+        const Float_t sigma  = camver70 ? fMcPed->GetPedestalRmsHigh(i) : fMcPed->GetElecNoise(i);
 
-        const Float_t sigma  = fMcPed->GetElecNoise(i);
-        //const Float_t sigrms = sigma/sqrt(num*2.);
-
-        //pix.SetPedestalRms(pedrms, sigrms);
         pix.Set(pedest, sigma);
     }
 
+    if (camver70)
+        fPedCam->SetReadyToSave();
+
     return kTRUE;
 }
-

trunk/MagicSoft/Mars/manalysis/MMcPedestalCopy.h

@@ -7 +7 @@
 
 class MMcFadcHeader;
+class MMcRunHeader;
 class MPedestalCam;
 
@@ -13 +14 @@
 private:
     const MMcFadcHeader *fMcPed;
+    const MMcRunHeader  *fMcRun;
+
     MPedestalCam *fPedCam;
 
     Bool_t CheckRunType(MParList *pList) const;
 
-    Int_t PreProcess(MParList *pList);
+    Int_t  PreProcess(MParList *pList);
     Bool_t ReInit(MParList *pList);
 

trunk/MagicSoft/Mars/manalysis/MMcPedestalNSBAdd.cc

@@ -111 +111 @@
 // --------------------------------------------------------------------------
 //
+// Check for the camera version. This class should not be used with
+// camera files >= 0.7
+//
+Bool_t MMcPedestalNSBAdd::CheckCamVersion(MParList *pList) const
+{
+    const MMcRunHeader *run = (MMcRunHeader*)pList->FindObject("MMcRunHeader");
+    if (!run)
+        *fLog << warn << dbginf << "MMcRunHeader not found... assuming camera<0.7" << endl;
+
+    return !run || run->GetCamVersion()<70;
+}
+
+// --------------------------------------------------------------------------
+//
 // - check whether we have a monte carlo file. if not skip this task
 // - try to get MMcFadcHeader, MGeomCam and MPedestalCam from the parameter
@@ -118 +132 @@
 Int_t MMcPedestalNSBAdd::PreProcess(MParList *pList)
 {
-    if (!CheckRunType(pList))
-    {
-        *fLog << warn << dbginf << "Warning - MMcPedestalNSB is for Monte Carlo files only... removing this task from list." << endl;
-        return kSKIP;
-    }
-
     fFadc = (MMcFadcHeader*)pList->FindObject("MMcFadcHeader");
     if (!fFadc)
@@ -141 +149 @@
     if (!fPedCam)
         return kFALSE;
-
-    const MMcRunHeader *mcrunheader = (MMcRunHeader*)pList->FindObject("MMcRunHeader");
-    if (!mcrunheader && fDnsbPixel<0)
-    {
-        *fLog << err << dbginf << "Using the default argument only ";
-        *fLog << "allowed if MMcRunHeader is available... aborting." << endl;
-        return kFALSE;
-    }
 
     return kTRUE;
@@ -202 +202 @@
         return kFALSE;
 
+    if (!CheckCamVersion(pList))
+        return kTRUE;
+
     const Float_t dnsbpix = GetDnsb(pList) * 50.0/15.0;
 
@@ -220 +223 @@
     }
 
+    fPedCam->SetReadyToSave();
+
     return kTRUE;
 }

trunk/MagicSoft/Mars/manalysis/Makefile

@@ -49 +49 @@
            MCerPhotEvt.cc \
            MCerPhotAnal.cc \
+           MCerPhotAnal2.cc \
            MCerPhotCalc.cc \
            MPedCalcPedRun.cc \
@@ -66 +67 @@
            MMinuitInterface.cc \
            MFiltercutsCalc.cc \
-           MCT1PadONOFF.cc \
-           MMcPedestalRead.cc
+           MCT1PadONOFF.cc
 
+#           MRelSignal.cc \
 
 SRCS    = $(SRCFILES)