Changeset 9343 for trunk/MagicSoft/Mars/mhflux

Timestamp:

02/15/09 23:12:04 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhflux

Files:

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

trunk/MagicSoft/Mars/mhflux/MHEffectiveOnTime.cc

@@ -88 +88 @@
 //   +  Double_t fTotalTime;
 //
+//  Class version 4:
+//  ----------------
+//   +  Double_t fEffOnTime;    //[s] On time as fitted from the DeltaT distribution
+//   +  Double_t fEffOnTimeErr; //[s] On time error as fitted from the DeltaT distribution
 //
 // ==========================================================================
@@ -301 +305 @@
 MHEffectiveOnTime::MHEffectiveOnTime(const char *name, const char *title)
     : fPointPos(0), fTime(0), fParam(0), fIsFinalized(kFALSE), 
-    fNumEvents(200*60), fFirstBin(3), fTotalTime(-1)
+    fNumEvents(200*60), fFirstBin(3), fTotalTime(-1),
+    fEffOnTime(-1), fEffOnTimeErr(0)
     //fNumEvents(2*60), fFirstBin(1)
 {
@@ -479 +484 @@
    fTotalTime = 0;
 
+   fEffOnTime    = -1;
+   fEffOnTimeErr =  0;
+
    return kTRUE;
 }
@@ -714 +722 @@
     fTime->Plus1ns();
 
-    *fLog << fLastTime << ":  Val=" << res[0] << "  Err=" << res[1] << endl;
+    *fLog << inf2 << fLastTime << ":  Val=" << res[0] << "  Err=" << res[1] << endl;
 }
 
@@ -780 +788 @@
     FitThetaBins();
     FitTimeBin();
+
+
+    TH1D *h = fH2DeltaT.ProjectionX("FinalizeProjDeltaT", -1, -1, "E");
+    Double_t res[7];
+    if (FitH(h, res))
+    {
+        fEffOnTime    = res[0];
+        fEffOnTimeErr = res[1];
+    }
+    delete h;
+
 
     fIsFinalized = kTRUE;

trunk/MagicSoft/Mars/mhflux/MHEffectiveOnTime.h

@@ -50 +50 @@
     UInt_t fFirstBin;    // For bin for fit
 
-    Double_t fTotalTime; //[s] Total time accumulated from run-header
+    Double_t fTotalTime;    //[s] Total time accumulated from run-header
+    Double_t fEffOnTime;    //[s] On time as fitted from the DeltaT distribution
+    Double_t fEffOnTimeErr; //[s] On time error as fitted from the DeltaT distribution
 
     // MHEffectiveOnTime
@@ -81 +83 @@
     // Getter
     Double_t GetTotalTime() const { return fTotalTime; }
+    Double_t GetEffOnTime() const { return fEffOnTime; }
 
     const TH1D &GetHEffOnTheta() const { return fHThetaEffOn; }