| 1 | #ifndef MARS_MCT1SupercutsCalc
|
|---|
| 2 | #define MARS_MCT1SupercutsCalc
|
|---|
| 3 |
|
|---|
| 4 | #ifndef MARS_MFilter
|
|---|
| 5 | #include "MFilter.h"
|
|---|
| 6 | #endif
|
|---|
| 7 |
|
|---|
| 8 | #ifndef ROOT_TArrayD
|
|---|
| 9 | #include <TArrayD.h>
|
|---|
| 10 | #endif
|
|---|
| 11 |
|
|---|
| 12 | class MParList;
|
|---|
| 13 | class MHillas;
|
|---|
| 14 | class MHillasSrc;
|
|---|
| 15 | class MMcEvt;
|
|---|
| 16 | class MCerPhotEvt;
|
|---|
| 17 | class MGeomCam;
|
|---|
| 18 | class MHadronness;
|
|---|
| 19 | class MHMatrix;
|
|---|
| 20 |
|
|---|
| 21 |
|
|---|
| 22 | class MCT1SupercutsCalc : public MTask
|
|---|
| 23 | {
|
|---|
| 24 | private:
|
|---|
| 25 | MHillas *fHil;
|
|---|
| 26 | MHillasSrc *fHilSrc;
|
|---|
| 27 | MMcEvt *fMcEvt;
|
|---|
| 28 | MHadronness *fHadronness; //! output container for hadronness
|
|---|
| 29 |
|
|---|
| 30 | TString fHadronnessName; // name of container to store hadronness
|
|---|
| 31 | TString fHilName;
|
|---|
| 32 | TString fHilSrcName;
|
|---|
| 33 |
|
|---|
| 34 | Double_t fMm2Deg;
|
|---|
| 35 |
|
|---|
| 36 | Int_t fMap[8];
|
|---|
| 37 | MHMatrix *fMatrix;
|
|---|
| 38 |
|
|---|
| 39 | //---------------------------------
|
|---|
| 40 | // cut parameters
|
|---|
| 41 | TArrayD fLengthUp;
|
|---|
| 42 | TArrayD fLengthLo;
|
|---|
| 43 | TArrayD fWidthUp;
|
|---|
| 44 | TArrayD fWidthLo;
|
|---|
| 45 | TArrayD fDistUp;
|
|---|
| 46 | TArrayD fDistLo;
|
|---|
| 47 | TArrayD fAsymUp;
|
|---|
| 48 | TArrayD fAsymLo;
|
|---|
| 49 | TArrayD fAlphaUp;
|
|---|
| 50 | //---------------------------------
|
|---|
| 51 |
|
|---|
| 52 | void InitParams();
|
|---|
| 53 |
|
|---|
| 54 | Int_t PreProcess(MParList *pList);
|
|---|
| 55 | Int_t Process();
|
|---|
| 56 |
|
|---|
| 57 | void Set(TArrayD &a, const TArrayD &b) { if (a.GetSize()==b.GetSize()) a=b; }
|
|---|
| 58 | Double_t GetVal(Int_t i) const;
|
|---|
| 59 |
|
|---|
| 60 | Double_t CtsMCut(
|
|---|
| 61 | #if ROOT_VERSION_CODE > ROOT_VERSION(3,05,00)
|
|---|
| 62 | const
|
|---|
| 63 | #endif
|
|---|
| 64 | TArrayD &a, Double_t ls, Double_t ct,
|
|---|
| 65 | Double_t ls2, Double_t dd2);
|
|---|
| 66 |
|
|---|
| 67 | public:
|
|---|
| 68 | MCT1SupercutsCalc(const char *hilname="MHillas",
|
|---|
| 69 | const char *hilsrcname="MHillasSrc",
|
|---|
| 70 | const char *name=NULL, const char *title=NULL);
|
|---|
| 71 |
|
|---|
| 72 | void SetParams(Double_t *par);
|
|---|
| 73 | void GetParams(Double_t *par);
|
|---|
| 74 |
|
|---|
| 75 | void SetHadronnessName(const TString name) { fHadronnessName = name; }
|
|---|
| 76 | TString GetHadronnessName() const { return fHadronnessName; }
|
|---|
| 77 |
|
|---|
| 78 | void InitMapping(MHMatrix *mat);
|
|---|
| 79 | void StopMapping() { InitMapping(NULL); }
|
|---|
| 80 |
|
|---|
| 81 | void SetLengthUp(const TArrayD &d) { Set(fLengthUp, d); }
|
|---|
| 82 | void SetLengthLo(const TArrayD &d) { Set(fLengthLo, d); }
|
|---|
| 83 | void SetWidthUp (const TArrayD &d) { Set(fWidthUp, d); }
|
|---|
| 84 | void SetWidthLo (const TArrayD &d) { Set(fWidthLo, d); }
|
|---|
| 85 | void SetDistUp (const TArrayD &d) { Set(fDistUp, d); }
|
|---|
| 86 | void SetDistLo (const TArrayD &d) { Set(fDistLo, d); }
|
|---|
| 87 | void SetAsymUp (const TArrayD &d) { Set(fAsymUp, d); }
|
|---|
| 88 | void SetAsymLo (const TArrayD &d) { Set(fAsymLo, d); }
|
|---|
| 89 | void SetAlphaUp (const TArrayD &d) { Set(fAlphaUp, d); }
|
|---|
| 90 |
|
|---|
| 91 | const TArrayD &GetLengthUp() const { return fLengthUp; }
|
|---|
| 92 | const TArrayD &GetLengthLo() const { return fLengthLo; }
|
|---|
| 93 | const TArrayD &GetWidthUp() const { return fWidthUp; }
|
|---|
| 94 | const TArrayD &GetWithLo() const { return fWidthLo; }
|
|---|
| 95 | const TArrayD &GetDistUp() const { return fDistUp; }
|
|---|
| 96 | const TArrayD &GetDistLo() const { return fDistLo; }
|
|---|
| 97 | const TArrayD &GetAsymUp() const { return fAsymUp; }
|
|---|
| 98 | const TArrayD &GetAsymLo() const { return fAsymLo; }
|
|---|
| 99 | const TArrayD &GetAlphaUp() const { return fAlphaUp; }
|
|---|
| 100 |
|
|---|
| 101 | ClassDef(MCT1SupercutsCalc, 0) // A class to evaluate the Supercuts
|
|---|
| 102 | };
|
|---|
| 103 |
|
|---|
| 104 | #endif
|
|---|
| 105 |
|
|---|
| 106 |
|
|---|
| 107 |
|
|---|
| 108 |
|
|---|
| 109 |
|
|---|
| 110 |
|
|---|
| 111 |
|
|---|
| 112 |
|
|---|
