Changeset 3889 for trunk/MagicSoft/Mars/mjobs/MJPedestal.cc

Timestamp:

04/29/04 15:48:38 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mjobs/MJPedestal.cc (modified) (6 diffs)

trunk/MagicSoft/Mars/mjobs/MJPedestal.cc

@@ -45 +45 @@
 #include "MTaskList.h"
 #include "MEvtLoop.h"
+#include "MExtractor.h"
 
 #include "MStatusDisplay.h"
@@ -62 +63 @@
 
 using namespace std;
-
+// --------------------------------------------------------------------------
+//
+// Default constructor. 
+//
+// Sets fRuns to 0, fExtractor to NULL, fDataCheck to kFALSE
+//
 MJPedestal::MJPedestal(const char *name, const char *title) 
-    : fRuns(0), fDataCheck(kFALSE)
-
+    : fRuns(0), fExtractor(NULL), fDataCheck(kFALSE)
 {
     fName  = name  ? name  : "MJPedestal";
@@ -115 +120 @@
 }
 
-void MJPedestal::DrawProjection(MHCamera *obj1, Int_t fit) const
-{
-    TH1D *obj2 = (TH1D*)obj1->Projection(obj1->GetName());
-    obj2->SetDirectory(0);
-    obj2->Draw();
-    obj2->SetBit(kCanDelete);
-    gPad->SetLogy();
-
-    if (obj1->GetGeomCam().InheritsFrom("MGeomCamMagic"))
-    {
-        TArrayI s0(3);
-        s0[0] = 6;
-        s0[1] = 1;
-        s0[2] = 2;
-
-        TArrayI s1(3);
-        s1[0] = 3;
-        s1[1] = 4;
-        s1[2] = 5;
-
-        TArrayI inner(1);
-        inner[0] = 0;
-
-        TArrayI outer(1);
-        outer[0] = 1;
-
-        // Just to get the right (maximum) binning
-        TH1D *half[4];
-        half[0] = obj1->ProjectionS(s0, inner, "Sector 6-1-2 Inner");
-        half[1] = obj1->ProjectionS(s1, inner, "Sector 3-4-5 Inner");
-        half[2] = obj1->ProjectionS(s0, outer, "Sector 6-1-2 Outer");
-        half[3] = obj1->ProjectionS(s1, outer, "Sector 3-4-5 Outer");
-
-        for (int i=0; i<4; i++)
-        {
-            half[i]->SetLineColor(kRed+i);
-            half[i]->SetDirectory(0);
-            half[i]->SetBit(kCanDelete);
-            half[i]->Draw("same");
-        }
-    }
-
-    const Double_t min   = obj2->GetBinCenter(obj2->GetXaxis()->GetFirst());
-    const Double_t max   = obj2->GetBinCenter(obj2->GetXaxis()->GetLast());
-    const Double_t integ = obj2->Integral("width")/2.5066283;
-    const Double_t mean  = obj2->GetMean();
-    const Double_t rms   = obj2->GetRMS();
-    const Double_t width = max-min;
-
-    if (rms==0 || width==0)
-        return;
-
-    const TString dgausformula("([0]-[3])/[2]*exp(-0.5*(x-[1])*(x-[1])/[2]/[2])"
-                               "+[3]/[5]*exp(-0.5*(x-[4])*(x-[4])/[5]/[5])");
-
-    TF1 *f=0;
-    switch (fit)
-    {
-        // FIXME: MAYBE add function to TH1?
-    case 0:
-        f = new TF1("sgaus", "gaus(0)", min, max);
-        f->SetLineColor(kYellow);
-        f->SetBit(kCanDelete);
-        f->SetParNames("Area", "#mu", "#sigma");
-        f->SetParameters(integ/rms, mean, rms);
-        f->SetParLimits(0, 0,   integ);
-        f->SetParLimits(1, min, max);
-        f->SetParLimits(2, 0,   width/1.5);
-        obj2->Fit(f, "QLRM");
-        break;
-
-    case 1:
-        f = new TF1("dgaus", dgausformula, min, max);
-        f->SetLineColor(kYellow);
-        f->SetBit(kCanDelete);
-        f->SetParNames("A_{tot}", "#mu_{1}", "#sigma_{1}", "A_{2}", "#mu_{2}", "#sigma_{2}");
-        f->SetParameters(integ,         (min+mean)/2, width/4,
-                         integ/width/2, (max+mean)/2, width/4);
-        // The left-sided Gauss
-        f->SetParLimits(0, integ-1.5,      integ+1.5);
-        f->SetParLimits(1, min+(width/10), mean);
-        f->SetParLimits(2, 0,              width/2);
-        // The right-sided Gauss
-        f->SetParLimits(3, 0,    integ);
-        f->SetParLimits(4, mean, max-(width/10));
-        f->SetParLimits(5, 0,    width/2);
-        obj2->Fit(f, "QLRM");
-        break;
-
-    default:
-        obj2->Fit("gaus", "Q");
-        obj2->GetFunction("gaus")->SetLineColor(kYellow);
-        break;
-    }
-}
-
-void MJPedestal::CamDraw(TCanvas &c, Int_t x, Int_t y, MHCamera &cam1, Int_t fit)
-{
-    c.cd(x);
-    gPad->SetBorderMode(0);
-    MHCamera *obj1=(MHCamera*)cam1.DrawCopy("hist");
-
-    obj1->AddNotify(&fPedestalCam);
-
-    c.cd(x+y);
-    gPad->SetBorderMode(0);
-    obj1->Draw();
-
-    c.cd(x+2*y);
-    gPad->SetBorderMode(0);
-    DrawProjection(obj1, fit);
-}
 
 void MJPedestal::DisplayResult(MParList &plist)
@@ -251 +144 @@
     // Create container to display
     //
-    MHCamera disp0(geomcam, "MPedestalCam;ped", "Pedestals");
-    MHCamera disp1(geomcam, "MPedestalCam;var", "Sigma Pedestal");
+    MHCamera disp0(geomcam, "MPedestalCam;ped", "Pedestal");
+    MHCamera disp1(geomcam, "MPedestalCam;RMS", "Ped. RMS");
 
     disp0.SetCamContent(fPedestalCam, 0);
@@ -269 +162 @@
     c3.Divide(2,3);
 
-    CamDraw(c3, 1, 2, disp0, 0);
-    CamDraw(c3, 2, 2, disp1, 1);
+    CamDraw(c3, 1, 2, disp0, 1);
+    CamDraw(c3, 2, 2, disp1, 6);
 }
 
@@ -368 +261 @@
     MPedCalcPedRun  pedcalc;
 
+    if (fExtractor)
+      {
+        pedcalc.SetWindowSize((Int_t)fExtractor->GetNumHiGainSamples());
+        pedcalc.SetRange(fExtractor->GetHiGainFirst(), fExtractor->GetHiGainLast());
+      }
+    else
+      {
+        *fLog << warn << GetDescriptor() 
+            << ": No extractor has been chosen, take default number of FADC samples " << endl;
+      }
+
     tlist.AddToList(&geomapl);
     //    tlist.AddToList(&merge);