source: releases/Mars.2014.05.26/mhflux/MHEffectiveOnTime.h@ 18066

Last change on this file since 18066 was 9346, checked in by tbretz, 16 years ago
*** empty log message ***
File size: 3.0 KB
Line 
1#ifndef MARS_MHEffectiveOnTime
2#define MARS_MHEffectiveOnTime
3
4#ifndef MARS_MH
5#include "MH.h"
6#endif
7#ifndef MARS_MTime
8#include "MTime.h"
9#endif
10#ifndef ROOT_TH1
11#include <TH1.h>
12#endif
13#ifndef ROOT_TH2
14#include <TH2.h>
15#endif
16#ifndef ROOT_TF1
17#include <TF1.h>
18#endif
19
20class MTime;
21class MPointingPos;
22class MParameterDerr;
23
24class MParList;
25
26class MHEffectiveOnTime : public MH
27{
28private:
29 MPointingPos *fPointPos; //! Container to get the zenith angle from
30 MTime fLastTime; //! Time-Stamp of last event
31
32 MTime *fTime; //! Time-Stamp of "effective on-time" event
33 MParameterDerr *fParam; //! Output container for effective on-time
34
35 TH2D fH2DeltaT; // Counts vs Delta T and Theta
36 TH1D fH1DeltaT; //! Counts vs Delta T (for a time interval)
37
38 TH1D fHThetaEffOn; // Effective On time versus Theta
39 TH1D fHThetaProb; // Chisq prob fit of Effective On time versus Theta
40 TH1D fHThetaNDF; // NDF vs Theta
41 TH1D fHThetaLambda; // Slope (rate) vs Theta
42
43 TH1D fHTimeEffOn; // Effective On time versus Time
44 TH1D fHTimeProb; // Chisq prob fit of Effective On time versus Time
45 TH1D fHTimeLambda; // Slope (rate) vs Time
46
47 Bool_t fIsFinalized; // Flag telling you whether fHThetaEffOn is the final result
48
49 UInt_t fNumEvents; // Number of events to be used for a bin in time
50 UInt_t fFirstBin; // For bin for fit
51
52 Double_t fTotalTime; //[s] Total time accumulated from run-header
53 Double_t fEffOnTime; //[s] On time as fitted from the DeltaT distribution
54 Double_t fEffOnTimeErr; //[s] On time error as fitted from the DeltaT distribution
55
56 // MHEffectiveOnTime
57 Bool_t FitH(TH1D *h, Double_t *res, Bool_t paint=kFALSE) const;
58 void FitThetaBins();
59 void FitTimeBin();
60 void PaintProb(TH1 &h) const;
61 void PaintText(Double_t val, Double_t error, Double_t range=-1) const;
62 void PaintText(Double_t *res) const;
63 void DrawRightAxis(const char *title);
64 void UpdateRightAxis(TH1 &h);
65 void PrintStatistics();
66
67 // MH
68 Bool_t SetupFill(const MParList *pList);
69 Bool_t ReInit(MParList *pList);
70 Int_t Fill(const MParContainer *par, const Stat_t w=1);
71 Bool_t Finalize();
72
73 // MParContainer
74 Int_t ReadEnv(const TEnv &env, TString prefix, Bool_t print);
75
76public:
77 MHEffectiveOnTime(const char *name=NULL, const char *title=NULL);
78
79 // Setter
80 void SetNumEvents(Int_t i) { fNumEvents=i; }
81 void SetFirstBin(Int_t i) { fFirstBin=i; }
82
83 // Getter
84 Double_t GetTotalTime() const { return fTotalTime; }
85 Double_t GetEffOnTime() const { return fEffOnTime; }
86
87 const TH1D &GetHEffOnTheta() const { return fHThetaEffOn; }
88 const TH1D &GetHEffOnTime() const { return fHTimeEffOn; }
89
90 const TH2D &GetH2DeltaT() const { return fH2DeltaT; }
91
92 // TObject
93 void Draw(Option_t *option="");
94 void Paint(Option_t *opt="");
95
96 ClassDef(MHEffectiveOnTime, 4) // Histogram to determin effective On-Time vs Time and Zenith Angle
97};
98
99#endif
100
Note: See TracBrowser for help on using the repository browser.