Changeset 17699 for trunk/Mars/fact/analysis/mc/callisto.C

Timestamp:

04/26/14 15:35:48 (11 years ago)

Author:

tbretz

Message:

Improved the calling convention. There is really no need to give so many arguments. Pipe some Monte Carlo values through. I wonder what special MC macros are for which omit all the Monte Carlo data in their output.

File:

• trunk/Mars/fact/analysis/mc/callisto.C (modified) (10 diffs)

trunk/Mars/fact/analysis/mc/callisto.C

@@ -1 @@
-using namespace std;
-
 #include <sstream>
 #include <iostream>
+
+#include "MLog.h"
+#include "MLogManip.h"
 
 #include "TH1F.h"
@@ -10 +11 @@
 #include "TLine.h"
 
+#include "../mcore/DrsCalib.h"
 #include "MDrsCalibration.h"
-//#include "MLogManip.h"
 #include "MExtralgoSpline.h"
 #include "MSequence.h"
@@ -45 +46 @@
 #include "MWriteAsciiFile.h"
 
+using namespace std;
+
 /* Maybe you wanna use this macro like this:
  * 
@@ -67 +70 @@
  */
 
-//using namespace std;
-
-int callisto(
-            const char *drs_calib_300_path="/fhgfs/groups/app/fact/mc_test/testingdrsfile/test300samples.drs.fits",
-            const char *pedestal_file_path="/fhgfs/groups/app/fact/mc_test/ceresfitstest/mcFilesForTests/mcNsbPedestal/00000001.001_P_MonteCarlo000_Events.fits",
-            const char *data_file_path="/fhgfs/groups/app/fact/mc_test/ceresfitstest/mcFilesForTests/mcProton/00000003.387_D_MonteCarlo010_Events.fits",
-            
-            const char *root_file_output_path = "/fhgfs/groups/app/callisto_star_test/callisto_for_mc_test_output/callisto.root",
-            const char *status_display_output_path = "/fhgfs/groups/app/callisto_star_test/callisto_for_mc_test_output/callisto_status_display.root", 
-            const char *status_display_title = "callisto_status_display")
+int callisto(const TString drsfile="test300samples.drs.fits",
+             const TString pedfile="00000001.001_P_MonteCarlo000_Events.fits",
+             const TString datfile="00000003.387_D_MonteCarlo010_Events.fits",
+             TString outfile = "",
+             TString displayfile = "", TString displaytitle = "")
 {
-    
+
+    // ======================================================
+
+    if (displaytitle.IsNull())
+        displaytitle = gSystem->BaseName(datfile);
+
+    FileStat_t fstat;
+    int rc = gSystem->GetPathInfo(outfile, fstat);
+    bool isdir = !rc || R_ISDIR(fstat.fMode);
+
+    const char *buf = gSystem->ConcatFileName(outfile, "callisto.root");
+    outfile = buf;
+    delete [] buf;
+
+    if (displayfile.IsNull())
+    {
+        displayfile = outfile;
+        displayfile.Insert(displayfile.Last('.'), "-display");
+    }
+    else
+    {
+        if (isdir && gSystem->DirName(displayfile)==TString("."))
+        {
+            buf = gSystem->ConcatFileName(outfile, displayfile);
+            displayfile = buf;
+            delete [] buf;
+        }
+    }
+
     // ======================================================
 
@@ -88 +114 @@
 
     // map file to use (get that from La Palma!)
-    const char *map = usemap ? "/fhgfs/groups/app/fact/resources/monte_carlo_FACTmap.txt" : NULL;
+    const char *map = usemap ? "TestForThomas/FACT/FACTmap111030.txt" : NULL;
 
     Bool_t maximum = kTRUE;
@@ -96 +122 @@
     //    "/cm/shared/apps/fact/Mars_svn_LP/template-lp-extractor-leading-edge.root";
 
-    const char *lp_template    = "/cm/shared/apps/fact/Mars_Trunk/template-lp.root";
-
-    //const char *pulse_template = "/cm/shared/apps/fact/Mars_svn_LP/template-pulse.root";
-
-    const char *pulse_template = "/cm/shared/apps/fact/Mars_Trunk/template-pulse.root";
+    const char *pulse_template = "TestForThomas/FACT/template-pulse.root";
 
     // ------------------------------------------------------
@@ -147 +169 @@
     }
 
-    TString datfile = TString(data_file_path);
-    TString drsfile = TString(drs_calib_300_path);
-    TString pedfile = TString(pedestal_file_path);
-
     gLog.Separator("Callisto");
-    gLog << all;    
-    gLog << "Data File        : " << datfile << "\n";
-    gLog << "DRS calib     300: " << drsfile << '\n';
+    gLog << all;
+    gLog << "Data File:     " << datfile << '\n';
+    gLog << "DRS calib 300: " << drsfile << endl;;
 
     MDrsCalibration drscalib300;
@@ -162 +180 @@
     }
     gLog << all;
-    gLog << "Pedestal     file: " << pedfile << '\n';
-
-    gLog << "root_file_output_path: " << root_file_output_path << endl;
-    gLog << "status_display_output_path: " << status_display_output_path << endl;
-    gLog << "status_display_title: " << status_display_title << endl;
+    gLog << "Pedestal file: " << pedfile << '\n';
+    gLog << "Output file:   " << outfile << '\n';
+    gLog << "Display file:  " << displayfile << '\n';
+    gLog << "Display title: " << displaytitle << endl;
 
     // ------------------------------------------------------
     MStatusArray arrt, arrp;
 
-    TFile ft(lp_template);
-    if (arrt.Read()<=0)
-    {
-        gLog << err << "ERROR - Reading LP template from " << lp_template << endl;
-        return 100;
-    }
-
-    MHCamera *lpref = (MHCamera*)arrt.FindObjectInCanvas("ExtCalSig;avg", "MHCamera", "Cam");
-    if (!lpref)
-    {
-        gLog << err << "ERROR - LP Template not found in " << lp_template << endl;
-        return 101;
-    }
-    lpref->SetDirectory(0);
+    // TFile ft(lp_template);
+    // if (arrt.Read()<=0)
+    // {
+    //     gLog << err << "ERROR - Reading LP template from " << lp_template << endl;
+    //     return 100;
+    // }
+
+    // MHCamera *lpref = (MHCamera*)arrt.FindObjectInCanvas("ExtCalSig;avg", "MHCamera", "Cam");
+    // if (!lpref)
+    // {
+    //     gLog << err << "ERROR - LP Template not found in " << lp_template << endl;
+    //     return 101;
+    // }
+    // lpref->SetDirectory(0);
 
     // MHCamera *gain = (MHCamera*)arrt.FindObjectInCanvas("gain", "MHCamera", "Gain");
@@ -619 +636 @@
     // The second rule is for the case reading raw-files!
     
-    MWriteRootFile write5(root_file_output_path, "RECREATE", "Calibrated Data", 2);
-    write5.AddContainer("MRawRunHeader", "RunHeaders");
-    write5.AddContainer("MGeomCam",      "RunHeaders");
-    write5.AddContainer("MSignalCam",    "Events");
-    write5.AddContainer("MTime",         "Events");
-    write5.AddContainer("MRawEvtHeader", "Events");
+    MWriteRootFile write5(outfile, "RECREATE", "Calibrated Data", 2);
+    write5.AddContainer("MRawRunHeader",       "RunHeaders");
+    write5.AddContainer("MGeomCam",            "RunHeaders");
+    write5.AddContainer("MMcCorsikaRunHeader", "RunHeaders", kFALSE);
+    write5.AddContainer("MCorsikaRunHeader",   "RunHeaders", kFALSE);
+    write5.AddContainer("MMcRunHeader",        "RunHeaders", kFALSE);
+
+    // Common events
+    write5.AddContainer("MCorsikaEvtHeader",   "Events", kFALSE);
+    write5.AddContainer("MMcEvt",              "Events", kFALSE);
+    write5.AddContainer("IncidentAngle",       "Events", kFALSE);
+    write5.AddContainer("MPointingPos",        "Events", kFALSE);
+    write5.AddContainer("MSignalCam",          "Events");
+    write5.AddContainer("MTime",               "Events");
+    write5.AddContainer("MRawEvtHeader",       "Events");
     //write.AddContainer("MTriggerPattern", "Events");
-    
     
     // ------------------ Setup histograms and fill tasks ----------------
@@ -669 +694 @@
         return 19;
 
-    d->SetTitle(status_display_title, kFALSE);
-    d->SaveAs(status_display_output_path);
+    d->SetTitle(displaytitle, kFALSE);
+    d->SaveAs(displayfile);
 
     return 0;