Changeset 4925

Timestamp:

09/10/04 15:31:34 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhist

Files:

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

trunk/MagicSoft/Mars/mhist/MHEffectiveOnTime.cc

@@ -28 +28 @@
 //  MHEffectiveOnTime
 //
-//  Fills a 2-D histogram Time-difference vs. Theta
+//  Filling this you will get the effective on-time versus theta and
+//  observation time.
 //
 //  From this histogram the effective on-time is determined by a fit.
@@ -34 +35 @@
 //  are stored in corresponding histograms
 //
+//  To determin the efective on time a poisson fit is done. For more details
+//  please have a look into the source code of FitH() it should be simple
+//  to understand. In this function a Delta-T distribution is fitted, while
+//  Delta-T is the time between two consecutive events.
+//
+//  The fit is done for projections of a 2D histogram in Theta and Delta T.
+//  So you get the effective on time versus theta.
+//
+//  To get the effective on-time versus time a histogram is filled with
+//  the Delta-T distribution of a number of events set by SetNumEvents().
+//  The default is 12000 (roughly 1min at 200Hz)
+//
+//  For each "time-bin" the histogram is fitted and the resulting effective
+//  on-time is stored in the fHEffOnTime histogram. Each entry in this
+//  histogram is the effective observation time between the upper and
+//  lower edges of the bins.
+//
+//  In addition the calculated effective on time is stored in a
+//  "MEffectiveOnTime [MParameterDerr]" and the corresponding time-stamp
+//  (the upper edge of the bin) "MTimeEffectiveOnTime [MTime]"
+//
+//  The class takes two binnings from the Parameter list; if these binnings
+//  are not available the defaultbinning is used:
+//      MBinning("BinningDeltaT"); // Units of seconds
+//      MBinning("BinningTheta");  // Units of degrees
+//
+//
+//  Usage:
+//  ------
+//    MFillH fill("MHEffectiveOnTime", "MTime");
+//    tlist.AddToList(&fill);
+//
+//
+//  Input Container:
+//    MPointingPos
+//
+//  Output Container:
+//    MEffectiveOnTime [MParameterDerr]
+//    MTimeEffectiveOnTime [MTime]
+//
 //////////////////////////////////////////////////////////////////////////////
-
 #include "MHEffectiveOnTime.h"
 
@@ -73 +113 @@
     //
     fName  = name  ? name  : "MHEffectiveOnTime";
-    fTitle = title ? title : "2-D histogram in Theta and time difference";
+    fTitle = title ? title : "Histogram to determin effective On-Time vs Time and Zenith Angle";
 
     // Main histogram
@@ -264 +304 @@
     h->Fit(&func, "0IRQ");
 
-    const Double_t chi2   = func.GetChisquare();
-    const Int_t    NDF    = func.GetNDF();
+    const Double_t chi2 = func.GetChisquare();
+    const Int_t    NDF  = func.GetNDF();
 
     // was fit successful ?
@@ -293 +333 @@
     gMinuit->mnemat((Double_t*)emat, 2);
 
-    const Double_t dldl   = emat[0][0];
+    const Double_t dldl  = emat[0][0];
     //const Double_t dN0dN0 = emat[1][1];
 
-    const Double_t teff   = Nm/lambda;
-    const Double_t dteff  = teff * TMath::Sqrt(dldl/(lambda*lambda) + 1.0/Nm);
-
-    const Double_t dl     = TMath::Sqrt(dldl);
+    const Double_t teff  = Nm/lambda;
+    const Double_t dteff = teff * TMath::Sqrt(dldl/(lambda*lambda) + 1.0/Nm);
+    const Double_t dl    = TMath::Sqrt(dldl);
 
     //const Double_t kappa  = Nm/N0;

trunk/MagicSoft/Mars/mhist/MHEffectiveOnTime.h

@@ -24 +24 @@
 {
 private:
-    MPointingPos   *fPointPos;   //!
-    MTime           fLastTime;   //!
+    MPointingPos   *fPointPos; //! Container to get the zenith angle from
+    MTime           fLastTime; //! Time-Stamp of last event
 
-    MTime          *fTime;       //!
-    MParameterDerr *fParam;      //!
+    MTime          *fTime;     //! Time-Stamp of "effective on-time" event
+    MParameterDerr *fParam;    //! Output container for effective on-time
 
-    TH2D fH2DeltaT;              // Counts vs Delta T and Theta
-    TH1D fH1DeltaT;              //! Counts vs Delta T (for a time interval)
+    TH2D fH2DeltaT;      // Counts vs Delta T and Theta
+    TH1D fH1DeltaT;      //! Counts vs Delta T (for a time interval)
 
-    TH1D fHEffOnTheta;           // Effective On time versus Theta
-    TH1D fHEffOnTime;            // Effective On time versus Time
+    TH1D fHEffOnTheta;   // Effective On time versus Theta
+    TH1D fHEffOnTime;    // Effective On time versus Time
 
-    TH1D fHProbTheta;            // Chisq prob fit of Effective On time versus Theta
-    TH1D fHProbTime;             // Chisq prob fit of Effective On time versus Time
+    TH1D fHProbTheta;    // Chisq prob fit of Effective On time versus Theta
+    TH1D fHProbTime;     // Chisq prob fit of Effective On time versus Time
 
     TH1D fHN0;
     TH1D fHLambda;
 
-    Bool_t fIsFinalized;
+    Bool_t fIsFinalized; // Flag telling you whether fHEffOnTheta is the final result
 
-    Int_t fNumEvents;
+    Int_t fNumEvents;    // Number of events to be used for a bin in time
 
-    const TString fNameProjDeltaT;  //! This should make sure, that gROOT doen't confuse the projection with something else
-    const TString fNameProjTheta;   //! This should make sure, that gROOT doen't confuse the projection with something else
+    const TString fNameProjDeltaT; //! This should make sure, that gROOT doen't confuse the projection with something else
+    const TString fNameProjTheta;  //! This should make sure, that gROOT doen't confuse the projection with something else
 
     Bool_t FitH(TH1D *h, Double_t *res, Bool_t paint=kFALSE) const;
@@ -63 +63 @@
     Bool_t Finalize();
 
+    const TH1D &GetHEffOnTheta() const { return fHEffOnTheta; }
+    const TH1D &GetHEffOnTime() const { return fHEffOnTime; }
+
     void Draw(Option_t *option="");
     void Paint(Option_t *opt="");
 
-    ClassDef(MHEffectiveOnTime, 1) // 2D-histogram to determin effective on-time vs. theta
+    ClassDef(MHEffectiveOnTime, 1) // Histogram to determin effective On-Time vs Time and Zenith Angle
 };