Changeset 2931 for trunk/MagicSoft/Mars/mcalib/MHCalibrationPINDiode.cc

Timestamp:

01/27/04 20:32:42 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MHCalibrationPINDiode.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mcalib/MHCalibrationPINDiode.cc

@@ -36 +36 @@
 #include <TH1.h>
 #include <TF1.h>
+#include <TPaveText.h>
 
 ClassImp(MHCalibrationPINDiode);
@@ -46 +47 @@
 //
 MHCalibrationPINDiode::MHCalibrationPINDiode(const char *name, const char *title)
-    : fVarGausFit(NULL)
+    : fChargeNbins(1000),
+      fTimeNbins(64),
+      fChargevsNbins(10000),
+      fTimeFirst(-0.25),
+      fTimeLast(31.75)
 {
 
@@ -53 +58 @@
 
     // Create a large number of bins, later we will rebin
-    fChargeFirstHiGain = -1000.;
-    fChargeLastHiGain  = gkStartPINDiodeBinNr;
-    fChargeNbinsHiGain = gkStartPINDiodeBinNr;
+    fChargeFirstHiGain = -100.5;
+    fChargeLastHiGain  = 1999.5;
 
-    fHPCharge = new TH1I("HPCharge","Distribution of Summed FADC Slices",
-                         fChargeNbinsHiGain,fChargeFirstHiGain,fChargeLastHiGain);
-    fHPCharge->SetXTitle("Sum FADC Slices");
-    fHPCharge->SetYTitle("Nr. of events");
-    fHPCharge->Sumw2();
-    fHPCharge->SetDirectory(0);
+    fHPINDiodeCharge   = new TH1F("HPINDiodeCharge","Distribution of Summed FADC Slices PINDiode",
+                                  fChargeNbins,fChargeFirstHiGain,fChargeLastHiGain);
+    fHPINDiodeCharge->SetXTitle("Sum FADC Slices");
+    fHPINDiodeCharge->SetYTitle("Nr. of events");
+    fHPINDiodeCharge->Sumw2();
+    fHPINDiodeCharge->SetDirectory(NULL);
 
-    fErrChargeFirst = 0.;
-    fErrChargeLast  = gkStartPINDiodeBinNr;
-    fErrChargeNbins = gkStartPINDiodeBinNr;
+    fHPINDiodeTime = new TH1F("HPINDiodeTime","Distribution of Mean Arrival Times PINDiode",
+                              fTimeNbins,fTimeFirst,fTimeLast);
+    fHPINDiodeTime->SetXTitle("Mean Arrival Times [FADC slice nr]");
+    fHPINDiodeTime->SetYTitle("Nr. of events");
+    fHPINDiodeTime->Sumw2();
+    fHPINDiodeTime->SetDirectory(NULL);
 
-    fHErrCharge = new TH1F("HErrCharge","Distribution of Variances of Summed FADC Slices",fErrChargeNbins,fErrChargeFirst,fErrChargeLast);
-    fHErrCharge->SetXTitle("Variance Summed FADC Slices");
-    fHErrCharge->SetYTitle("Nr. of events");
-    fHErrCharge->Sumw2();
-    fHErrCharge->SetDirectory(0);
+    fHPINDiodeChargevsN = new TH1I("HPINDiodeChargevsN","Sum of Hi Gain Charges vs. Event Number Pixel ",
+                                 fChargevsNbins,-0.5,(Axis_t)fChargevsNbins - 0.5);
+    fHPINDiodeChargevsN->SetXTitle("Event Nr.");
+    fHPINDiodeChargevsN->SetYTitle("Sum of Hi Gain FADC slices");
+    fHPINDiodeChargevsN->SetDirectory(NULL);
 
-    Int_t tfirst = 0;
-    Int_t tlast  = 31;
-    Int_t nbins   = 32;
-
-    fHPTime = new TH1F("HPTime","Distribution of Mean Arrival Times",nbins,tfirst,tlast);
-    fHPTime->SetXTitle("Mean Arrival Times [FADC slice nr]");
-    fHPTime->SetYTitle("Nr. of events");
-    fHPTime->Sumw2();
-    fHPTime->SetDirectory(0);
-
+    Clear();
 }
 
@@ -89 +87 @@
 {
 
-  delete fHPCharge;
-  delete fHErrCharge;
+  delete fHPINDiodeCharge;
+  delete fHPINDiodeTime;
+  delete fHPINDiodeChargevsN;
+
+  if (fChargeGausFit)
+    delete fChargeGausFit;
+  if (fTimeGausFit)
+    delete fTimeGausFit;
+  if (fFitLegend)
+    delete fFitLegend;
+}
+
+void MHCalibrationPINDiode::Clear(Option_t *o)
+{
   
-  if (fVarGausFit)
-    delete fVarGausFit;
-  
-  delete fHPTime;
+  fTotalEntries            = 0;
+
+  fChargeFirstHiGain       = -100.5;
+  fChargeLastHiGain        = 1999.5;
+
+  fChargeChisquare         = -1.;
+  fChargeProb              = -1.;
+  fChargeNdf               = -1;
+  fTimeChisquare           = -1.;
+  fTimeProb                = -1.;
+  fTimeNdf                 = -1;
+  fTimeMean                = -1.;
+  fTimeSigma               = -1.;
+
+  fTimeLowerFitRangeHiGain = 0;
+  fTimeUpperFitRangeHiGain = 0;
+  fTimeLowerFitRangeLoGain = 0;
+  fTimeUpperFitRangeLoGain = 0;
+
+  if (fChargeGausFit)
+    delete fChargeGausFit;
+  if (fTimeGausFit)
+    delete fTimeGausFit;
+  if (fFitLegend)
+    delete fFitLegend;
+
+  return;
 
 }
 
-const Double_t MHCalibrationPINDiode::GetTime() const
+void MHCalibrationPINDiode::Reset()
 {
-    return fVarGausFit->GetParameter(2);
+
+  Clear();
+
+  fHPINDiodeCharge->Reset();
+  fHPINDiodeTime->Reset();
+  fHPINDiodeChargevsN->Reset();
 }
 
-const Double_t MHCalibrationPINDiode::GetErrTime() const
+
+Bool_t MHCalibrationPINDiode::FillCharge(Float_t q)
 {
-    return fVarGausFit->GetParameter(3);
+    return (fHPINDiodeCharge->Fill(q) > -1);
 }
+
+Bool_t MHCalibrationPINDiode::FillTime(Float_t t)
+{
+    return (fHPINDiodeTime->Fill(t)   > -1);
+}
+
+Bool_t MHCalibrationPINDiode::FillChargevsN(Float_t q, Int_t n)
+{
+    return (fHPINDiodeChargevsN->Fill(n,q) > -1);
+}
+
+void MHCalibrationPINDiode::CutAllEdges()
+{
+
+  Int_t nbins = 30;
+
+  CutEdges(fHPINDiodeCharge,nbins);
+
+  fChargeFirstHiGain = fHPINDiodeCharge->GetBinLowEdge(fHPINDiodeCharge->GetXaxis()->GetFirst());
+  fChargeLastHiGain  = fHPINDiodeCharge->GetBinLowEdge(fHPINDiodeCharge->GetXaxis()->GetLast())
+                      +fHPINDiodeCharge->GetBinWidth(0);
+  CutEdges(fHPINDiodeChargevsN,0);
+
+}
+