Changeset 2662 for trunk/MagicSoft/Mars/mhist

Timestamp:

12/12/03 09:41:11 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhist

Files:

• MHCalibrationPixel.cc (modified) (8 diffs)
• MHCalibrationPixel.h (modified) (4 diffs)

trunk/MagicSoft/Mars/mhist/MHCalibrationPixel.cc

@@ -65 +65 @@
 MHCalibrationPixel::MHCalibrationPixel(const char *name, const char *title)
       : fPixId(-1),
-        fHiGainvsLoGain(NULL),
         fTotalEntries(0),
         fHiGains(NULL),
@@ -71 +70 @@
         fChargeGausFit(NULL), 
         fTimeGausFit(NULL), 
+        fHivsLoGain(NULL),
         fFitLegend(NULL),
         fLowerFitRange(-2000.),
@@ -166 +166 @@
   delete fHChargevsNLoGain;
 
-  if (fHiGainvsLoGain)
-    delete fHiGainvsLoGain;
-
   delete fHiGains;
   delete fLoGains;
@@ -178 +175 @@
   if (fFitLegend)
     delete fFitLegend;
+  if (fHivsLoGain)
+    delete fHivsLoGain;
+
 
 }
@@ -379 +379 @@
 {
 
-  if (!fHiGainvsLoGain)
-    {
-
-    // Create TGraph, we set the points later, get hold of the number of points with fTotalEntries
-    fHiGainvsLoGain = new TGraph(fTotalEntries,fHiGains->GetArray(),fLoGains->GetArray());
-    fHiGainvsLoGain->SetName("HiGainvsLoGain");
-    fHiGainvsLoGain->SetTitle("Plot the High Gain vs. Low Gain");
-    fHiGainvsLoGain->GetXaxis()->Set(300,0.,1500.);
-    fHiGainvsLoGain->GetYaxis()->Set(400,0.,2000.);
-    fHiGainvsLoGain->GetXaxis()->SetTitle("Sum of Charges High Gain");
-    fHiGainvsLoGain->GetYaxis()->SetTitle("Sum of Charges Low Gain");
-    fHiGainvsLoGain->SetMarkerStyle(7);
-
-    }
-
-  TProfile *hist = new TProfile("PePe","Shit is still moving",100,0.,1000.,0.,1000.);
-
-  for (Int_t i=0;i<fTotalEntries;i++)
-      hist->Fill(fHiGains->At(i),fLoGains->At(i),1);
-
+  if (!fHivsLoGain)
+    FitHiGainvsLoGain();
+      
+  
   gStyle->SetOptFit(0);
   gStyle->SetOptStat(1111111);
-    
-    TCanvas *c = MakeDefCanvas(this,600,900); 
-
-    gROOT->SetSelectedPad(NULL);
-
+  
+  TCanvas *c = MakeDefCanvas(this,600,900); 
+  
+  gROOT->SetSelectedPad(NULL);
+  
     c->Divide(2,4);
 
@@ -440 +424 @@
         gROOT->SetSelectedPad(NULL);
         gStyle->SetOptFit();
-        //      fHiGainvsLoGain->DrawClone("Apq")->SetBit(kCanDelete);
-        //      fHiGainvsLoGain->Fit("p1","q");
-        hist->Draw("prof");
-        hist->Fit("pol1","pepe","markus",hist->GetMean()-3.*hist->GetRMS(),hist->GetMean()+2.*hist->GetRMS());
-
-        fOffset = hist->GetFunction("pepe")->GetParameter(0);
-        fSlope  = hist->GetFunction("pepe")->GetParameter(1);
-
-        cout << "Pepe: " << fOffset << " Markus: " << fSlope << endl;
+        fHivsLoGain->Draw("prof");
 
 
@@ -484 +460 @@
         gROOT->SetSelectedPad(NULL);
         gStyle->SetOptFit();
-        hist->Draw("prof");
-        hist->Fit("p1","pepe","markus",hist->GetMean()-3.*hist->GetRMS(),hist->GetMean()+3.*hist->GetRMS());
-        //      fHiGainvsLoGain->DrawClone("Apq")->SetBit(kCanDelete);
-        //      fHiGainvsLoGain->Fit("p1","q");
+        fHivsLoGain->Draw("prof");
         gPad->Modified();
         gPad->Update();
@@ -563 +536 @@
 }
 
+void MHCalibrationPixel::FitHiGainvsLoGain()
+{
+
+  if (fHivsLoGain)
+    return;
+
+  gStyle->SetOptFit();
+
+  fHivsLoGain = new TProfile("HiGainvsLoGain","Plot the High Gain vs. Low Gain",100,0.,1000.,0.,1000.);
+  fHivsLoGain->GetXaxis()->SetTitle("Sum of Charges High Gain");
+  fHivsLoGain->GetYaxis()->SetTitle("Sum of Charges Low Gain");
+  
+  for (Int_t i=0;i<fTotalEntries;i++)
+    fHivsLoGain->Fill(fHiGains->At(i),fLoGains->At(i),1);
+  
+  fHivsLoGain->Fit("pol1","rq","",fHivsLoGain->GetMean()-2.5*fHivsLoGain->GetRMS(),fHivsLoGain->GetMean()+2.*fHivsLoGain->GetRMS());
+  
+  fOffset = fHivsLoGain->GetFunction("pol1")->GetParameter(0);
+  fSlope  = fHivsLoGain->GetFunction("pol1")->GetParameter(1);
+  
+}
 
 

trunk/MagicSoft/Mars/mhist/MHCalibrationPixel.h

@@ -18 +18 @@
 #endif
 
-#ifndef ROOT_TGraph
-#include "TGraph.h"
+#ifndef ROOT_TProfile
+#include "TProfile.h"
 #endif
 
@@ -36 +36 @@
 
   Int_t fPixId;                  // Pixel Nr
-  TGraph *fHiGainvsLoGain;       //! Plot the HiGain vs. the LoGain
   Int_t fTotalEntries;           // Number of entries
 
@@ -55 +54 @@
   TF1* fTimeGausFit;
   
+  TProfile* fHivsLoGain;
+
   TPaveText *fFitLegend;  
   
@@ -102 +103 @@
 
   void   SetPointInGraph(Float_t qhi, Float_t qlo);
+  void   FitHiGainvsLoGain();
 
   Bool_t FillChargeLoGain(Float_t q)             { return (fHChargeLoGain->Fill(q) > -1); }