Changeset 3982 for trunk/MagicSoft/Mars/mjobs

Timestamp:

05/05/04 17:26:14 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mjobs

Files:

• MJExtractCalibTest.cc (modified) (6 diffs)
• MJExtractCalibTest.h (modified) (5 diffs)
• MJPedestal.cc (modified) (1 diff)
• Makefile (modified) (1 diff)

trunk/MagicSoft/Mars/mjobs/MJExtractCalibTest.cc

@@ -48 +48 @@
 #include "MBadPixelsCam.h"
 #include "MCerPhotEvt.h"
+#include "MArrivalTime.h"
 #include "MCalibrationChargeCam.h"
+#include "MCalibrationRelTimeCam.h"
 #include "MCalibrationQECam.h"
 #include "MHCamEvent.h"
@@ -56 +58 @@
 #include "MExtractSlidingWindow.h"
 #include "MExtractor.h"
+#include "MExtractTime.h"
+#include "MExtractTimeFastSpline.h"
 #include "MFillH.h"
 #include "MCalibrate.h"
+#include "MCalibrateRelTimes.h"
 #include "MPedPhotCalc.h"
 #include "MWriteRootFile.h"
@@ -70 +75 @@
 // Default constructor. 
 //
-// Sets fRuns to 0, fExtractor to NULL
+// Sets fRuns to 0, fExtractor to NULL, fTimeExtractor to NULL
 //
 MJExtractCalibTest::MJExtractCalibTest(const char *name, const char *title) 
-    : fRuns(NULL), fExtractor(NULL)
+    : fRuns(NULL), fExtractor(NULL), fTimeExtractor(NULL)
 {
     fName  = name  ? name  : "MJExtractCalibTest";
@@ -151 +156 @@
 }
 
+TString MJExtractCalibTest::GetOutputFileT() const
+{
+  if (!fRuns)
+    return "";
+  
+  return Form("%s/%s-F4.root", (const char*)fOutputPath, (const char*)fRuns->GetRunsAsFileName());
+}
+
 TString MJExtractCalibTest::GetOutputFileP() const
 {
@@ -240 +253 @@
     }
 
-  if (TestBit(kEnableGraphicalOutput))
+  //  if (TestBit(kEnableGraphicalOutput))
     tlist.AddToList(&fill);
   tlist.AddToList(&photcalc);
   tlist.AddToList(&fillcam);
-  //  tlist.AddToList(&write);
   
   // Create and setup the eventloop
@@ -267 +279 @@
   return kTRUE;
 }
+
+Bool_t MJExtractCalibTest::ProcessT(MPedestalCam &pedcam, MCalibrationRelTimeCam &relcam)
+{
+  const TString fname = GetOutputFileT();
+  
+  if (gSystem->AccessPathName(fname, kFileExists))
+    return ProcessFileT(pedcam,relcam);
+  
+  return kTRUE;
+}
+
+Bool_t MJExtractCalibTest::ProcessFileT(MPedestalCam &pedcam, MCalibrationRelTimeCam &relcam)
+{
+
+  if (!fRuns)
+    {
+      *fLog << err << "No Runs choosen... abort." << endl;
+      return kFALSE;
+    }
+  if (fRuns->GetNumRuns() != fRuns->GetNumEntries())
+    {
+      *fLog << err << "Number of files found doesn't match number of runs... abort." << endl;
+      return kFALSE;
+    }
+  
+  *fLog << inf;
+  fLog->Separator(GetDescriptor());
+  *fLog << "Calculate MExtractedSignalCam from Runs " << fRuns->GetRunsAsString() << endl;
+  *fLog << endl;
+  
+  MArrivalTime  arrtime;
+
+  // Setup Lists
+  MParList plist;
+  plist.AddToList(&pedcam);
+  plist.AddToList(&relcam);
+  plist.AddToList(&arrtime);
+  plist.AddToList(&fTestCam);
+  plist.AddToList(&fBadPixels);
+  
+  MTaskList tlist;
+  plist.AddToList(&tlist);
+
+  // Setup Task-lists
+  MReadMarsFile read("Events");
+  read.DisableAutoScheme();
+  static_cast<MRead&>(read).AddFiles(*fRuns);
+  
+  MGeomApply             apply; // Only necessary to craete geometry
+  MExtractTimeFastSpline extract;
+  MCalibrateRelTimes     timecalc;
+
+  MHCamEvent evt("ExtTimes");
+  evt.SetType(0);
+  MFillH fill(&evt, "MArrivalTimeCam");
+  
+  MFillH fillcam("MHCalibrationTestTimeCam", "MArrivalTime");
+  fillcam.SetNameTab("TestTime");
+
+  /*
+  MWriteRootFile write(GetOutputFileD(), "RECREATE", fRuns->GetRunsAsString(), 2);
+  write.AddContainer("MExtractedSignalCam", "Events");
+  write.AddContainer("MTime",               "Events");
+  write.AddContainer("MRawEvtHeader",       "Events");
+  write.AddContainer("MPedestalCam",        "RunHeaders");
+  write.AddContainer("MRawRunHeader",       "RunHeaders");
+  write.AddContainer("MBadPixelsCam",       "RunHeaders");
+  */
+
+  tlist.AddToList(&read);
+  tlist.AddToList(&apply);
+
+  if (fTimeExtractor)
+    tlist.AddToList(fTimeExtractor);
+    else
+    {
+      *fLog << warn << GetDescriptor() 
+            << ": No extractor has been chosen, take default MExtractSlidingWindow " << endl;
+      tlist.AddToList(&extract);
+    }
+
+  //  if (TestBit(kEnableGraphicalOutput))
+  tlist.AddToList(&fill);
+  tlist.AddToList(&timecalc);
+  tlist.AddToList(&fillcam);
+  
+  // Create and setup the eventloop
+  MEvtLoop evtloop(fName);
+  evtloop.SetParList(&plist);
+  evtloop.SetDisplay(fDisplay);
+  evtloop.SetLogStream(fLog);
+  
+  // Execute first analysis
+  if (!evtloop.Eventloop())
+    {
+      *fLog << err << GetDescriptor() << ": Failed." << endl;
+      return kFALSE;
+    }
+  
+  tlist.PrintStatistics();
+  
+  DisplayResult(plist);
+
+  *fLog << inf << GetDescriptor() << ": Done." << endl;
+  
+  return kTRUE;
+}
+
 
 Bool_t MJExtractCalibTest::ReadPedPhotCam()

trunk/MagicSoft/Mars/mjobs/MJExtractCalibTest.h

@@ -18 +18 @@
 class MParList;
 class MPedestalCam;
+class MCalibrationRelTimeCam;
 class MCalibrationChargeCam;
 class MCalibrationQECam;
 class MExtractor;
+class MExtractTime;
 class MJExtractCalibTest : public MParContainer, public MGCamDisplays
 {
@@ -27 +29 @@
   TString fOutputPath;
 
-  MRunIter   *fRuns;                                       // Data files
-  MExtractor *fExtractor;                                  // Signal extractor
+  MRunIter     *fRuns;                       // Data files
+  MExtractor   *fExtractor;                  // Signal extractor
+  MExtractTime *fTimeExtractor;              // Arrival time extractor  
   
   MBadPixelsCam        fBadPixels;
@@ -40 +43 @@
   Bool_t ProcessFileD(MPedestalCam &pedcam, MCalibrationChargeCam &calcam, MCalibrationQECam &qecam);
   Bool_t ProcessFileP(MPedestalCam &pedcam, MCalibrationChargeCam &calcam, MCalibrationQECam &qecam);
+  Bool_t ProcessFileT(MPedestalCam &pedcam, MCalibrationRelTimeCam &relcam);  
 
 public:
@@ -49 +53 @@
   TString GetOutputFileP() const;
   TString GetOutputFileD() const;
+  TString GetOutputFileT() const;
   
   MHCalibrationTestCam &GetTestCam()    { return fTestCam; }
@@ -55 +60 @@
   
   void SetExtractor(MExtractor* ext)              { fExtractor = ext; }
+  void SetTimeExtractor(MExtractTime* ext)         { fTimeExtractor = ext; }
   void SetBadPixels(const MBadPixelsCam &bad) { bad.Copy(fBadPixels); }
   
   Bool_t ProcessD(MPedestalCam &pedcam, MCalibrationChargeCam &calcam, MCalibrationQECam &qecam);
   Bool_t ProcessP(MPedestalCam &pedcam, MCalibrationChargeCam &calcam, MCalibrationQECam &qecam);
+  Bool_t ProcessT(MPedestalCam &pedcam, MCalibrationRelTimeCam &relcam);  
   
   ClassDef(MJExtractCalibTest, 0) // Tool to extract, calibrate and test the signal 

trunk/MagicSoft/Mars/mjobs/MJPedestal.cc

@@ -258 +258 @@
     MGeomApply      geomapl;
     // MBadPixelsMerge merge(&fBadPixels);
-    //MExtractSignal sigcalc;
     MPedCalcPedRun  pedcalc;
 

trunk/MagicSoft/Mars/mjobs/Makefile

@@ -21 +21 @@
 INCLUDES = -I. -I../mbase -I../mgui -I../mgeom -I../mdata -I../mhbase \
            -I../mfileio -I../mfilter -I../manalysis -I../mhist -I../mcalib \
-           -I../mbadpixels -I../msignal -I../mraw -I../mpedestal
+           -I../mbadpixels -I../msignal -I../mraw -I../mpedestal -I../mtools
 
 #manalysis: MChisqEval (MParameters)