Changeset 2661 for trunk/MagicSoft/Mars/mhist

Timestamp:

12/12/03 09:04:01 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhist

Files:

• MHCalibrationPixel.cc (modified) (5 diffs)
• MHCalibrationPixel.h (modified) (2 diffs)

trunk/MagicSoft/Mars/mhist/MHCalibrationPixel.cc

@@ -45 +45 @@
 #include <TF1.h>
 #include <TH2.h>
+#include <TProfile.h>
 #include <TCanvas.h>
 #include <TPad.h>
@@ -85 +86 @@
         fTimeProb(-1.),
         fTimeNdf(-1),
-        fUseLoGain(kFALSE)
+        fUseLoGain(kFALSE),
+        fOffset(0.),
+        fSlope(0.)
 { 
 
@@ -391 +394 @@
     }
 
-    gStyle->SetOptFit(0);
-    gStyle->SetOptStat(1111111);
+  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);
+
+  gStyle->SetOptFit(0);
+  gStyle->SetOptStat(1111111);
     
     TCanvas *c = MakeDefCanvas(this,600,900); 
@@ -432 +440 @@
         gROOT->SetSelectedPad(NULL);
         gStyle->SetOptFit();
-        fHiGainvsLoGain->DrawClone("Apq")->SetBit(kCanDelete);
-        fHiGainvsLoGain->Fit("p1","q");
+        //      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;
+
+
         gPad->Modified();
         gPad->Update();
@@ -467 +484 @@
         gROOT->SetSelectedPad(NULL);
         gStyle->SetOptFit();
-        fHiGainvsLoGain->DrawClone("Apq")->SetBit(kCanDelete);
-        fHiGainvsLoGain->Fit("p1","q");
+        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");
         gPad->Modified();
         gPad->Update();

trunk/MagicSoft/Mars/mhist/MHCalibrationPixel.h

@@ -85 +85 @@
 
   Bool_t fUseLoGain;
+
+  Double_t fOffset;
+  Double_t fSlope;
   
   virtual void DrawLegend();
@@ -137 +140 @@
   const TH1I *GetHChargevsN()               { return fHChargevsNHiGain; }
   const TH1I *GetHChargevsN()         const { return fHChargevsNHiGain; }
+
+  Double_t GetOffset()  { return fOffset; }
+  Double_t GetSlope()   { return fSlope;  }
   
   Bool_t FitChargeHiGain(Option_t *option="RQ0");