Changeset 4906 for trunk/MagicSoft/Mars/mhist

Timestamp:

09/09/04 17:51:57 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHEffectiveOnTime.cc (modified) (19 diffs)

trunk/MagicSoft/Mars/mhist/MHEffectiveOnTime.cc

@@ -63 +63 @@
 //
 MHEffectiveOnTime::MHEffectiveOnTime(const char *name, const char *title)
-    : fPointPos(0), fTime(0), fParam(0), fIsFinalized(kFALSE), fInterval(60)
+    : fPointPos(0), fTime(0), fParam(0), fIsFinalized(kFALSE), fInterval(60),
+    fNameProjDeltaT(Form("DeltaT_%p", this)), fNameProjTheta(Form("Theta_%p", this))
 {
     //
@@ -72 +73 @@
 
     // Main histogram
-    fHTimeDiff.SetXTitle("\\Delta t [s]");
-    fHTimeDiff.SetYTitle("\\Theta [\\circ]");
-    fHTimeDiff.UseCurrentStyle();
-    fHTimeDiff.SetDirectory(NULL);
+    fHDeltaT.SetName("DeltaT");
+    fHDeltaT.SetXTitle("\\Delta t [s]");
+    fHDeltaT.SetYTitle("\\Theta [\\circ]");
+    fHDeltaT.UseCurrentStyle();
+    fHDeltaT.SetDirectory(NULL);
 
     // effective on time versus theta
-    fHEffOn.SetName("EffOnTime");
-    fHEffOn.SetTitle("Effective On Time T_{eff}");
-    fHEffOn.SetXTitle("\\Theta [\\circ]");
-    fHEffOn.SetYTitle("T_{eff} [s]");
-    fHEffOn.UseCurrentStyle();
-    fHEffOn.SetDirectory(NULL);
-    fHEffOn.GetYaxis()->SetTitleOffset(1.2);
+    fHEffOnTheta.SetName("EffOnTime");
+    fHEffOnTheta.SetTitle("Effective On Time T_{eff}");
+    fHEffOnTheta.SetXTitle("\\Theta [\\circ]");
+    fHEffOnTheta.SetYTitle("T_{eff} [s]");
+    fHEffOnTheta.UseCurrentStyle();
+    fHEffOnTheta.SetDirectory(NULL);
+    fHEffOnTheta.GetYaxis()->SetTitleOffset(1.2);
     //fHEffOn.Sumw2();
+
+    // effective on time versus time
+    fHEffOnTime.SetName("EffOnTime");
+    fHEffOnTime.SetTitle("Effective On Time T_{eff}");
+    fHEffOnTime.SetXTitle("Time");
+    fHEffOnTime.SetYTitle("T_{eff} [s]");
+    fHEffOnTime.UseCurrentStyle();
+    fHEffOnTime.SetDirectory(NULL);
+    fHEffOnTime.GetYaxis()->SetTitleOffset(1.2);
+    fHEffOnTime.GetXaxis()->SetLabelSize(0.033);
+    fHEffOnTime.GetXaxis()->SetTimeFormat("%H:%M:%S %F1995-01-01 00:00:00 GMT");
+    fHEffOnTime.GetXaxis()->SetTimeDisplay(1);
+    fHEffOnTime.Sumw2();
 
     // chi2/NDF versus theta
@@ -128 +143 @@
     // setup binning
     MBinning btheta("BinningTheta");
-    btheta.SetEdgesCos(51, 0, 60);
-
-    MBinning btime("BinningTimeDiff");
+    btheta.SetEdgesCos(101, 0, 60);
+
+    MBinning btime("BinningDeltaT");
     btime.SetEdges(50, 0, 0.1);
 
-    MH::SetBinning(&fHTimeDiff, &btime, &btheta);
-
-    btheta.Apply(fHEffOn);
+    MH::SetBinning(&fHDeltaT, &btime, &btheta);
+
+    btheta.Apply(fHEffOnTheta);
     btheta.Apply(fHChi2);
     btheta.Apply(fHLambda);
@@ -141 +156 @@
     btheta.Apply(fHProb);
 }
-
-// FIXME: Just for a preliminary check
-static Double_t testval = 0;
-static Double_t testerr = 0;
 
 // --------------------------------------------------------------------------
@@ -167 +178 @@
        return kFALSE;
 
-   const MBinning* binsdtime = (MBinning*)plist->FindObject("BinningTimeDiff");
+   const MBinning* binsdtime = (MBinning*)plist->FindObject("BinningDeltaT");
    const MBinning* binstheta = (MBinning*)plist->FindObject("BinningTheta");
    if (!binstheta || !binsdtime)
@@ -173 +184 @@
    else
    {
-       SetBinning(&fHTimeDiff, binsdtime, binstheta);
-
-       binstheta->Apply(fHEffOn);
+       SetBinning(&fHDeltaT, binsdtime, binstheta);
+
+       binstheta->Apply(fHEffOnTheta);
        binstheta->Apply(fHChi2);
        binstheta->Apply(fHLambda);
@@ -200 +211 @@
 
     // nbins = number of Theta bins
-    const Int_t nbins = fHTimeDiff.GetNbinsY();
+    const Int_t nbins = fHDeltaT.GetNbinsY();
 
     TH1D *h=0;
@@ -206 +217 @@
     {
         //        TH1D &h = *hist->ProjectionX("Calc-theta", i, i);
-        h = fHTimeDiff.ProjectionX(name, i, i, "E");
+        h = fHDeltaT.ProjectionX(name, i, i, "E");
 
         Double_t res[7];
@@ -213 +224 @@
 
         // the effective on time is Nm/lambda
-        fHEffOn.SetBinContent(i, res[0]);
-        fHEffOn.SetBinError  (i, res[1]);
+        fHEffOnTheta.SetBinContent(i, res[0]);
+        fHEffOnTheta.SetBinError  (i, res[1]);
 
         // plot chi2-probability of fit
@@ -346 +357 @@
         TH1D *h0=0;
 
-        padsave->cd(1);
-        if ((h0 = (TH1D*)gPad->FindObject("TimeDiff")))
+        padsave->GetPad(1)->cd(1);
+        if ((h0 = (TH1D*)gPad->FindObject(fNameProjDeltaT)))
         {
-            h0 = fHTimeDiff.ProjectionX("TimeDiff", -1, 9999, "E");
+            h0 = fHDeltaT.ProjectionX(fNameProjDeltaT, -1, 9999, "E");
             if (h0->GetEntries()>0)
                 gPad->SetLogy();
         }
 
-        padsave->cd(2);
-        if ((h0 = (TH1D*)gPad->FindObject("Theta")))
-            fHTimeDiff.ProjectionY("Theta", -1, 9999, "E");
+        padsave->GetPad(2)->cd(1);
+        if ((h0 = (TH1D*)gPad->FindObject(fNameProjTheta)))
+            fHDeltaT.ProjectionY(fNameProjTheta, -1, 9999, "E");
 
         if (!fIsFinalized)
             FitThetaBins();
     }
-    if (o==(TString)"paint")
+
+    TH1D *h=0;
+    if (o==(TString)"theta")
+        h = &fHEffOnTheta;
+    if (o==(TString)"time")
+        h = &fHEffOnTime;
+
+    if (h)
     {
         Double_t error = 0;
-        for (int i=0; i<fHEffOn.GetXaxis()->GetNbins(); i++)
-            error += fHEffOn.GetBinError(i);
-
-        TLatex text(0.45, 0.94, Form("T_{eff} = %.1fs \\pm %.1fs", fHEffOn.Integral(), error));
+        for (int i=0; i<h->GetXaxis()->GetNbins(); i++)
+            error += h->GetBinError(i);
+
+        TLatex text(0.45, 0.94, Form("T_{eff} = %.1fs \\pm %.1fs", h->Integral(), error));
         text.SetBit(TLatex::kTextNDC);
         text.SetTextSize(0.04);
@@ -386 +404 @@
     AppendPad("fit");
 
-    pad->Divide(2,2);
+    pad->Divide(2, 1, 0, 0);
 
     TH1 *h;
@@ -392 +410 @@
     pad->cd(1);
     gPad->SetBorderMode(0);
-    h = fHTimeDiff.ProjectionX("TimeDiff", -1, 9999, "E");
+    gPad->Divide(1, 2, 0, 0);
+    pad->GetPad(1)->cd(1);
+    gPad->SetBorderMode(0);
+    h = fHDeltaT.ProjectionX(fNameProjDeltaT, -1, 9999, "E");
     h->SetTitle("Distribution of \\Delta t [s]");
     h->SetXTitle("\\Delta t [s]");
@@ -401 +422 @@
     h->Draw();
 
+    pad->GetPad(1)->cd(2);
+    gPad->SetBorderMode(0);
+    fHEffOnTime.Draw();
+    AppendPad("time");
+
     pad->cd(2);
     gPad->SetBorderMode(0);
-    h = fHTimeDiff.ProjectionY("Theta", -1, 9999, "E");
+    gPad->Divide(1, 3, 0, 0);
+
+    pad->GetPad(2)->cd(1);
+    gPad->SetBorderMode(0);
+    h = fHDeltaT.ProjectionY(fNameProjTheta, -1, 9999, "E");
     h->SetTitle("Distribution of  \\Theta [\\circ]");
     h->SetXTitle("\\Theta [\\circ]");
@@ -413 +443 @@
     h->Draw();
 
-    pad->cd(3);
-    gPad->SetBorderMode(0);
-    fHEffOn.Draw();
-    AppendPad("paint");
-
-    pad->cd(4);
+    pad->GetPad(2)->cd(2);
     gPad->SetBorderMode(0);
     fHProb.Draw();
+
+    pad->GetPad(2)->cd(3);
+    gPad->SetBorderMode(0);
+    fHEffOnTheta.Draw();
+    AppendPad("theta");
+}
+
+void MHEffectiveOnTime::FillTimeBin(Double_t val, Double_t err)
+{
+    // Get x-axis
+    TAxis &x = *fHEffOnTime.GetXaxis();
+
+    // Get number of bins
+    const Int_t n = x.GetNbins();
+
+    // Fill last bin with new values
+    fHEffOnTime.SetBinContent(n, fParam->GetVal());
+    fHEffOnTime.SetBinError(n, fParam->GetErr());
+
+    // Get time range of histogram, get bin width
+    const Double_t lo = fHEffOnTime.GetXaxis()->GetXmin();
+    const Double_t up = fHEffOnTime.GetXaxis()->GetXmax();
+    const Double_t w  = fHEffOnTime.GetXaxis()->GetBinWidth(1);
+
+    // Enhance histogram by one bin
+    MBinning bins;
+    bins.SetEdges(n+2, lo, up+w);
+    bins.Apply(fHEffOnTime);
+
+    // Transform overflow bin
+    // fHEffOnTime.SetBinContent(n+2, fHEffOnTime.GetBinContent(n+1));
+    // fHEffOnTime.SetBinError(n+2, fHEffOnTime.GetBinError(n+1));
 }
 
 void MHEffectiveOnTime::Calc()
 {
-    TH1D *h = fHTimeDiff.ProjectionX("", -1, 99999, "E");
+    TH1D *h = fHDeltaT.ProjectionX("", -1, 99999, "E");
     h->SetDirectory(0);
 
@@ -438 +495 @@
     const Double_t error = res[1];
 
-    fParam->SetVal(val-fEffOnTime0, error-fEffOnErr0);
+    fParam->SetVal(val>fEffOnTime0?val-fEffOnTime0:0, error-fEffOnErr0);
     fParam->SetReadyToSave();
 
-    testval += fParam->GetVal();
-    testerr += fParam->GetErr();
+    FillTimeBin(fParam->GetVal(), fParam->GetErr());
 
     fEffOnTime0 = val;
@@ -453 +509 @@
     *fLog << Form("T_{eff} = %.1fs \\pm %.1fs",
                   fParam->GetVal(), fParam->GetErr());
-    *fLog << Form("   %.1f %.1f   %.1f %.1f",
-                  val, testval, error, testerr) << endl;
+    *fLog << Form("   %.1f   %.1f",
+                  val, error) << endl;
 
     fTime->AddMilliSeconds(fInterval*1000);
@@ -468 +524 @@
     if (!time)
     {
-        *fLog << err << "ERROR - MHEffectiveOnTimeTime::Fill without argument or container doesn't inherit from MTime... abort." << endl;
+        *fLog << err << "ERROR - MHEffectiveOnTime::Fill without argument or container doesn't inherit from MTime... abort." << endl;
         return kFALSE;
     }
@@ -483 +539 @@
         *fTime = *time;
 
+        MBinning bins;
+        bins.SetEdges(2, fTime->GetAxisTime(), fTime->GetAxisTime()+fInterval);
+        bins.Apply(fHEffOnTime);
+
         // Make this just a ns before the first event
         fTime->Minus1ns();
@@ -493 +553 @@
     }
 
-    fHTimeDiff.Fill(*time-fLastTime, fPointPos->GetZd(), w);
+    fHDeltaT.Fill(*time-fLastTime, fPointPos->GetZd(), w);
     fLastTime = *time;