Changeset 17641 for trunk/Mars/fact

Timestamp:

04/01/14 14:38:58 (11 years ago)

Author:

ftemme

Message:

Added the calculation of Disp and ThetaSq in the MC-Star macro and commented the lightpulser calibration in the MC-Callisto macro out

Location:

trunk/Mars/fact/analysis/mc

Files:

• callisto.C (modified) (2 diffs)
• star.C (modified) (4 diffs)

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

@@ -92 +92 @@
     Bool_t maximum = kTRUE;
 
-    const char *lp_template    = maximum ?
-        "/cm/shared/apps/fact/Mars_svn_LP/template-lp-extractor-maximum.root" :
-        "/cm/shared/apps/fact/Mars_svn_LP/template-lp-extractor-leading-edge.root";
-
-    const char *pulse_template = "/cm/shared/apps/fact/Mars_svn_LP/template-pulse.root";
+    //const char *lp_template    = maximum ?
+    //    "/cm/shared/apps/fact/Mars_svn_LP/template-lp-extractor-maximum.root" :
+    //    "/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";
 
     // ------------------------------------------------------
@@ -182 +186 @@
     lpref->SetDirectory(0);
 
-    MHCamera *gain = (MHCamera*)arrt.FindObjectInCanvas("gain", "MHCamera", "Gain");
-    if (!gain)
-    {
-        gLog << err << "ERROR - Gain not found in " << lp_template << endl;
-        return 101;
-    }
-    gain->SetDirectory(0);
+    // MHCamera *gain = (MHCamera*)arrt.FindObjectInCanvas("gain", "MHCamera", "Gain");
+    // if (!gain)
+    // {
+    //     gLog << err << "ERROR - Gain not found in " << lp_template << endl;
+    //     return 101;
+    // }
+    // gain->SetDirectory(0);
 
     TFile fp(pulse_template);

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

@@ -42 +42 @@
 #include "MParameters.h"
 #include "MWriteAsciiFile.h"
+#include "MFMagicCuts.h"
 
 #include "MMuonSetup.h"
@@ -288 +289 @@
     write.AddContainer("MNewImagePar",    "Events");
     write.AddContainer("MRawEvtHeader",   "Events");
+    write.AddContainer("ThetaSquared",   "Events");
+    write.AddContainer("Disp",           "Events");
     write.AddContainer("MRawRunHeader",   "RunHeaders");
     write.AddContainer("MGeomCam",        "RunHeaders");
@@ -294 +297 @@
     fillmhn2.SetFilter(&fmuonhn);
 
+    TArrayD param(12);
+    // Parametrization of Disp
+    param[0]  =  1.15136;     // constant
+    param[8]  =  0.0681437;   // slope
+    param[9]  =  2.62932;     // leak
+    param[10] =  1.51279;     // size-offset
+    param[11] =  0.0507821;   // size-slope
+    // Parametrization for sign of disp (m3long, slope)
+    param[5]  = -0.07;
+    param[6]  =  7.2;
+    param[7]  =  0.5;
+    // ThetaSq-Cut
+    param[1]  =  0.11;  // 0.215
+    // Area-Cut
+    param[2]  =  0.215468;
+    param[3]  =  5.63973;
+    param[4]  =  0.0836169;
+    MFMagicCuts cuts;
+    cuts.SetHadronnessCut(MFMagicCuts::kNoCut);
+    cuts.SetVariables(param);
+
+
     tlist.AddToList(&read);
     tlist.AddToList(&cont);
@@ -299 +324 @@
     tlist.AddToList(&clean);
     tlist.AddToList(&hcalc);
+    tlist.AddToList(&cuts);
     tlist.AddToList(&fillC1);
     tlist.AddToList(&fillC2);