Changeset 8929 for trunk/MagicSoft/Mars/mjobs/MJStar.cc

Timestamp:

06/08/08 13:55:36 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mjobs/MJStar.cc (modified) (10 diffs)

trunk/MagicSoft/Mars/mjobs/MJStar.cc

@@ -18 +18 @@
 !   Author(s): Thomas Bretz, 1/2004 <mailto:tbretz@astro.uni-wuerzburg.de>
 !
-!   Copyright: MAGIC Software Development, 2000-2007
+!   Copyright: MAGIC Software Development, 2000-2008
 !
 !
@@ -35 +35 @@
 #include <TFile.h>
 
+// Core
 #include "MLog.h"
 #include "MLogManip.h"
@@ -45 +46 @@
 #include "MStatusDisplay.h"
 
+// Histograms
 #include "MHSectorVsTime.h"
 #include "MHEffectiveOnTime.h"
@@ -50 +52 @@
 #include "MBinning.h"
 
+///NEEWWWW
+#include "MH3.h"
+
+// Tasks
 #include "MReadReports.h"
 #include "MReadMarsFile.h"
-#include "MFDataPhrase.h"
-#include "MFilterList.h"
-#include "MFDataMember.h"
-#include "MFDeltaT.h"
-#include "MFSoftwareTrigger.h"
 #include "MContinue.h"
 #include "MGeomApply.h"
@@ -69 +70 @@
 #include "MWriteRootFile.h"
 
+// Filter
+#include "MFDataPhrase.h"
+#include "MFTriggerPattern.h"
+#include "MFilterList.h"
+#include "MFDataMember.h"
+#include "MFDeltaT.h"
+#include "MFSoftwareTrigger.h"
+
+// Parameter container
 #include "MMuonSetup.h"
 #include "MObservatory.h"
+//#include "MTriggerPattern.h"
 #include "MPointingPosCalc.h"
 
@@ -242 +253 @@
     fill0d.SetNameTab("Sparks");
 
+    // FIXME: This should be replaced by a cut working for all triggers!
+    MTriggerPattern par;   // An empty trigger pattern for older files
+    plist.AddToList(&par); // without MTriggerPattern stored in the file
+
+    // For the effective on-time calculation we don't want SUM-only events
+    MFTriggerPattern fsum;
+    fsum.SetDefault(kTRUE);
+    fsum.DenySumTrigger();
+    fsum.AllowTriggerLvl1();
+    fill9.SetFilter(&fsum);
+
+    // Plot the trigger pattern rates vs. run-number
+    MH3 hrate("MRawRunHeader.fRunNumber", "MTriggerPattern.GetUnprescaled");
+    hrate.SetWeight("1./MRawRunHeader.GetRunLength");
+    hrate.SetName("Rate");
+    hrate.SetTitle("Event rate after cleaning [Hz];Run Number;Trigger Type;");
+    hrate.InitLabels(MH3::kLabelsX);
+    hrate.DefaultLabelY("ERROR");
+    hrate.DefineLabelY( 0, "0");        //  0: No pattern
+    hrate.DefineLabelY( 1, "Trig");     //  1: Lvl1
+    hrate.DefineLabelY( 4, "Trig");     //  4: Lvl2
+    hrate.DefineLabelY( 5, "Trig");     //     Lvl2+Lvl1
+    hrate.DefineLabelY(32, "Sum");      // 32: Sum
+    hrate.DefineLabelY(33, "Trig");     //     Sum+Lvl1
+    hrate.DefineLabelY(36, "Trig");     //     Sum+Lvl2
+    hrate.DefineLabelY(37, "Trig");     //     Sum+Lvl1+Lvl2
+
+    MFillH frate(&hrate, "", "FillRate");
+    frate.SetDrawOption("box");
+
     // ------------------ Setup write task ----------------
 
@@ -260 +301 @@
     write.AddContainer("MRawEvtHeader",             "Events");
     write.AddContainer("MPointingPos",              "Events");
+    write.AddContainer("MTriggerPattern",           "Events");
 
     // Run Header
@@ -333 +375 @@
         tlist2.AddToList(&rate);
         tlist2.AddToList(&fillvs);
+        tlist2.AddToList(&fsum);
         tlist2.AddToList(&fill9);
         tlist2.AddToList(&writet);
@@ -358 +401 @@
     tlist2.AddToList(&fill0d);
     tlist2.AddToList(&fill0b);
+    tlist2.AddToList(&frate);
     tlist2.AddToList(&fill1);
     tlist2.AddToList(&fill2);
@@ -409 +453 @@
     MHSectorVsTime histdc, histrms;
     histdc.SetNameTime("MTimeCurrents");
-    histdc.SetTitle("Mean of all DC Currents;;<I> [nA]");
+    histdc.SetTitle("Average DC currents of all pixels vs time;;<I> [nA]");
     histdc.SetMinimum(0);
     histdc.SetMaximum(10);
     histrms.SetNameTime("MTimeCurrents");
-    histrms.SetTitle("Mean pedestal rms of all pixels;;<\\sigma_{p}> [phe]");
+    histrms.SetTitle("Average pedestal rms of all pixels vs time;;<\\sigma_{p}> [phe]");
     histrms.SetType(5);
     histrms.SetMinimum(0);