| 1 | #ifndef MARS_MHAlpha
|
|---|
| 2 | #define MARS_MHAlpha
|
|---|
| 3 |
|
|---|
| 4 | #ifndef MARS_MH
|
|---|
| 5 | #include "MH.h"
|
|---|
| 6 | #endif
|
|---|
| 7 |
|
|---|
| 8 | #ifndef MARS_MTime
|
|---|
| 9 | #include "MTime.h"
|
|---|
| 10 | #endif
|
|---|
| 11 |
|
|---|
| 12 | #ifndef MARS_MAlphaFitter
|
|---|
| 13 | #include "MAlphaFitter.h"
|
|---|
| 14 | #endif
|
|---|
| 15 |
|
|---|
| 16 | #ifndef ROOT_TH3
|
|---|
| 17 | #include <TH3.h>
|
|---|
| 18 | #endif
|
|---|
| 19 |
|
|---|
| 20 | class MParList;
|
|---|
| 21 | class MParameterD;
|
|---|
| 22 | class MEnergyEst;
|
|---|
| 23 | class MHMatrix;
|
|---|
| 24 | class MPointingPos;
|
|---|
| 25 |
|
|---|
| 26 | /*
|
|---|
| 27 | class MAlphaFitter : public MParContainer
|
|---|
| 28 | {
|
|---|
| 29 | private:
|
|---|
| 30 | Float_t fSigInt;
|
|---|
| 31 | Float_t fSigMax;
|
|---|
| 32 | Float_t fBgMin;
|
|---|
| 33 | Float_t fBgMax;
|
|---|
| 34 | Int_t fPolynom;
|
|---|
| 35 |
|
|---|
| 36 | Double_t fSignificance;
|
|---|
| 37 | Double_t fExcessEvents;
|
|---|
| 38 | Double_t fChiSqSignal;
|
|---|
| 39 | Double_t fChiSqBg;
|
|---|
| 40 | Double_t fSigmaGaus;
|
|---|
| 41 | Double_t fIntegralMax;
|
|---|
| 42 |
|
|---|
| 43 | public:
|
|---|
| 44 | MAlphaFitter() : fSigInt(10), fSigMax(75), fBgMin(45), fBgMax(85), fPolynom(2)
|
|---|
| 45 | {
|
|---|
| 46 | }
|
|---|
| 47 |
|
|---|
| 48 | MAlphaFitter(const MAlphaFitter &f)
|
|---|
| 49 | {
|
|---|
| 50 | f.Copy(*this);
|
|---|
| 51 | }
|
|---|
| 52 |
|
|---|
| 53 | void Copy(TObject &o) const
|
|---|
| 54 | {
|
|---|
| 55 | MAlphaFitter &f = static_cast<MAlphaFitter&>(o);
|
|---|
| 56 |
|
|---|
| 57 | f.fSigInt = fSigInt;
|
|---|
| 58 | f.fSigMax = fSigMax;
|
|---|
| 59 | f.fBgMin = fBgMin;
|
|---|
| 60 | f.fBgMax = fBgMax;
|
|---|
| 61 | f.fPolynom = fPolynom;
|
|---|
| 62 | }
|
|---|
| 63 |
|
|---|
| 64 | void SetSignalIntegralMax(Float_t s) { fSigInt = s; }
|
|---|
| 65 | void SetSignalFitMax(Float_t s) { fSigMax = s; }
|
|---|
| 66 | void SetBackgroundFitMin(Float_t s) { fBgMin = s; }
|
|---|
| 67 | void SetBackgroundFitMax(Float_t s) { fBgMax = s; }
|
|---|
| 68 | void SetNumPolynom(Int_t s) { fPolynom = s; }
|
|---|
| 69 |
|
|---|
| 70 | Double_t GetExcessEvents() const { return fExcessEvents; }
|
|---|
| 71 | Double_t GetSignificance() const { return fSignificance; }
|
|---|
| 72 | Double_t GetChiSqSignal() const { return fChiSqSignal; }
|
|---|
| 73 | Double_t GetChiSqBg() const { return fChiSqBg; }
|
|---|
| 74 | Double_t GetSigmaGaus() const { return fSigmaGaus; }
|
|---|
| 75 |
|
|---|
| 76 | void PaintResult(Float_t x=0.04, Float_t y=0.94, Float_t size=0.035) const;
|
|---|
| 77 |
|
|---|
| 78 | Bool_t Fit(TH1D &h, Bool_t paint=kFALSE);
|
|---|
| 79 | ClassDef(MAlphaFitter, 1)
|
|---|
| 80 | };
|
|---|
| 81 | */
|
|---|
| 82 |
|
|---|
| 83 | class MHAlpha : public MH
|
|---|
| 84 | {
|
|---|
| 85 | private:
|
|---|
| 86 | MAlphaFitter fFit;
|
|---|
| 87 |
|
|---|
| 88 | TH3D fHAlpha; // Alpha vs. theta and energy
|
|---|
| 89 | TH1D fHEnergy; // excess events vs energy
|
|---|
| 90 | TH1D fHTheta; // excess events vs theta
|
|---|
| 91 | TH1D fHTime; // excess events vs time;
|
|---|
| 92 | TH1D fHAlphaTime; //! temporary histogram to get alpha vs. time
|
|---|
| 93 |
|
|---|
| 94 | MParameterD *fResult; //!
|
|---|
| 95 | MEnergyEst *fEnergy; //!
|
|---|
| 96 | MPointingPos *fPointPos; //!
|
|---|
| 97 |
|
|---|
| 98 | MTime *fTimeEffOn; //! Time to steer filling of fHTime
|
|---|
| 99 | MTime *fTime; //! Event-Time to finalize fHTime
|
|---|
| 100 | MTime fLastTime; //! Last fTimeEffOn
|
|---|
| 101 |
|
|---|
| 102 | const TString fNameProjAlpha; //! This should make sure, that gROOT doen't confuse the projection with something else
|
|---|
| 103 |
|
|---|
| 104 | MHMatrix *fMatrix; //!
|
|---|
| 105 | Int_t fMap[2]; //!
|
|---|
| 106 |
|
|---|
| 107 | void FitEnergySpec(Bool_t paint=kFALSE);
|
|---|
| 108 | void FitEnergyBins(Bool_t paint=kFALSE);
|
|---|
| 109 | void FitThetaBins(Bool_t paint=kFALSE);
|
|---|
| 110 |
|
|---|
| 111 | void UpdateAlphaTime(Bool_t final=kFALSE);
|
|---|
| 112 | void InitAlphaTime(const MTime &t);
|
|---|
| 113 | void FinalAlphaTime(MBinning &bins);
|
|---|
| 114 |
|
|---|
| 115 | void PaintText(Double_t val, Double_t error) const;
|
|---|
| 116 |
|
|---|
| 117 | public:
|
|---|
| 118 | MHAlpha(const char *name=NULL, const char *title=NULL);
|
|---|
| 119 |
|
|---|
| 120 | Bool_t SetupFill(const MParList *pl);
|
|---|
| 121 | Bool_t Fill(const MParContainer *par, const Stat_t w=1);
|
|---|
| 122 | Bool_t Finalize();
|
|---|
| 123 |
|
|---|
| 124 | const TH3D &GetHist() const { return fHAlpha; }
|
|---|
| 125 | const MAlphaFitter &GetAlphaFitter() const { return fFit; }
|
|---|
| 126 |
|
|---|
| 127 | void Paint(Option_t *opt="");
|
|---|
| 128 | void Draw(Option_t *option="");
|
|---|
| 129 |
|
|---|
| 130 | void InitMapping(MHMatrix *mat);
|
|---|
| 131 | void StopMapping();
|
|---|
| 132 |
|
|---|
| 133 | ClassDef(MHAlpha, 1) // Alpha-Plot which is fitted online
|
|---|
| 134 | };
|
|---|
| 135 |
|
|---|
| 136 | #endif
|
|---|
