Changeset 17641 for trunk/Mars/fact/analysis/mc/star.C

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

File:

• trunk/Mars/fact/analysis/mc/star.C (modified) (4 diffs)

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);