Changeset 2405 for trunk/MagicSoft/Mars/macros

Timestamp:

10/20/03 15:09:44 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/status.C (modified) (10 diffs)

trunk/MagicSoft/Mars/macros/status.C

@@ -31 +31 @@
 /////////////////////////////////////////////////////////////////////////////
 
-void status(const char *fname="*1947*cosmic*.root")
+void status()
 {
     //
@@ -57 +57 @@
     plist.AddToList(&tlist);
 
-    /*
-     MSrcPosCam src;
-     src.SetXY(1./geomcam.GetConvMm2Deg(), 0);
-     plist.AddToList(&src);
-     */
+    // ------------- user change -----------------
+    MDirIter files("~MAGIC/online_data/rootdata", "*Mkn501-I*.root", -1);
+    //files.Print("all");
 
     //
@@ -69 +67 @@
     MReadMarsFile read("Events");
     read.DisableAutoScheme();
-
-    // ------------- user change -----------------
-    read.AddFile(fname);
+    read.AddFiles(files);
 
     MGeomApply        geomapl;
@@ -79 +75 @@
     MCerPhotCalc      ncalc;
 
-    TArrayS blinds(2);
-    blinds[0] = 271;
-    blinds[1] = 291;
+    TArrayS blinds(0);
 
     MBlindPixelCalc blind;
-    //blind.SetPixelIndices(blinds);
+    blind.SetPixelIndices(blinds);
     blind.SetUseInterpolation();
 
@@ -91 +85 @@
     MHillasCalc       hcalc;
     MHillasSrcCalc    scalc; // !!Preliminary!! Will be removed later!
-    MCT1SupercutsCalc calc1;
 
     MHCamEvent hist("PedestalRms");
@@ -98 +91 @@
 
     // -------------------------------------------
+    MHTriggerLvl0 trighi(254, "SaturationHi", "Saturation Rate of Hi Gains");
+    trighi.SetType(1);
+
+    MFillH fillhi(&trighi, "MRawEvtData");
     MFillH hfill0("Uncleaned [MHCamEvent]", "MCerPhotEvt");
     MFillH hfill1("Pedestals [MHCamEvent]", "MPedestalCam");
@@ -108 +105 @@
     MFillH hfill8("MHStarMap", "MHillas");
     MFillH hfill9("Cleaned [MHCamEvent]", "MCerPhotEvt");
-    MFillH hfill10("MHHadronness", "MHadronness");
-    MFillH hfill11("MHSigmaTheta");
+
+    MContinue cont1("MNewImagePar.fNumCorePixels<0");
 
     tlist.AddToList(&read);
     tlist.AddToList(&geomapl);
+    tlist.AddToList(&fillhi);
     tlist.AddToList(&pcopy);
     tlist.AddToList(&pnsb);
@@ -121 +119 @@
     tlist.AddToList(&clean);
     tlist.AddToList(&hcalc);
+    tlist.AddToList(&cont1);
     tlist.AddToList(&scalc);
-    //tlist.AddToList(&calc1);
     tlist.AddToList(&hfill1);
     tlist.AddToList(&hfill2);
@@ -132 +130 @@
     tlist.AddToList(&hfill8);
     tlist.AddToList(&hfill9);
-    //tlist.AddToList(&hfill10);
 
     MEvtLoop evtloop;
@@ -146 +143 @@
     tlist.PrintStatistics();
 
-    MHCamera &uncl = *((MHCamEvent*)plist.FindObject("Uncleaned"))->GetHistByName();
-    MHCamera &hped = *((MHCamEvent*)plist.FindObject("Pedestals"))->GetHistByName();
-    MHCamera &hrms = *((MHCamEvent*)plist.FindObject("PedestalRms"))->GetHistByName();
-
-    for (int i=0; i<577; i++)
-    {
-        if (uncl.IsUsed(i))
-        {
-            if (uncl[i+1]>uncl.GetMean()+3*uncl.GetRMS())
-                cout << "Mean Charge of Pixel-Index #" << i << ": " << uncl[i+1] << " > 3*rms" << endl;
-            if (uncl[i+1]==0)
-                cout << "Mean Charge of Pixel-Index #" << i << ": " << uncl[i+1] << " <= 0" << endl;
-            if (uncl[i+1]<uncl.GetMean()-3*uncl.GetRMS())
-                cout << "Mean Charge of Pixel-Index #" << i << ": " << uncl[i+1] << " < 3*rms" << endl;
-        }
-        if (hped.IsUsed(i))
-        {
-            if (hped[i+1]>hped.GetMean()+1.5*hped.GetRMS())
-                cout << "Mean Pedestal of Pixel-Index #" << i << ": " << hped[i+1] << " > 1.5*rms" << endl;
-            if (hped[i+1]==0)
-                cout << "Mean Pedestal of Pixel-Index #" << i << ": " << hped[i+1] << " <= 0" << endl;
-            if (hped[i+1]<hped.GetMean()-1.5*hped.GetRMS())
-                cout << "Mean Pedestal of Pixel-Index #" << i << ": " << hped[i+1] << " < 1.5*rms" << endl;
-        }
-        if (hrms.IsUsed(i))
-        {
-            if (hrms[i+1]>hrms.GetMean()+4*hrms.GetRMS())
-                cout << "Mean PedestalRMS of Pixel-Index #" << i << ": " << hrms[i+1] << " > 4*rms" << endl;
-            if (hrms[i+1]==0)
-                cout << "Mean PedestalRMS of Pixel-Index #" << i << ": " << hrms[i+1] << " <= 0" << endl;
-            if (hrms[i+1]<hrms.GetMean()-4*hrms.GetRMS())
-                cout << "Mean PedestalRMS of Pixel-Index #" << i << ": " << hrms[i+1] << " < 4*rms" << endl;
-        }
-    }
+    MHCamEvent *uncl = (MHCamEvent*)plist.FindObject("Uncleaned");
+    MHCamEvent *hped = (MHCamEvent*)plist.FindObject("Pedestals");
+    MHCamEvent *hrms = (MHCamEvent*)plist.FindObject("PedestalRms");
+    uncl->PrintOutliers(3);
+    hped->PrintOutliers(1.5);
+    hrms->PrintOutliers(4);
+    trighi.PrintOutliers(2.5);
 
     //