| 1 | #ifndef MARS_MMySuperCutsCalc
|
|---|
| 2 | #define MARS_MMySuperCutsCalc
|
|---|
| 3 |
|
|---|
| 4 | #ifndef MARS_MTask
|
|---|
| 5 | #include "MTask.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 MHillasExt;
|
|---|
| 16 | class MNewImagePar;
|
|---|
| 17 | class MMcEvt;
|
|---|
| 18 | class MCerPhotEvt;
|
|---|
| 19 | class MGeomCam;
|
|---|
| 20 | class MHadronness;
|
|---|
| 21 | class MHMatrix;
|
|---|
| 22 | class MMySuperCuts;
|
|---|
| 23 | class MPointingPos;
|
|---|
| 24 | class MNewImagePar;
|
|---|
| 25 | class MHillasExt;
|
|---|
| 26 |
|
|---|
| 27 | class MMySuperCutsCalc : public MTask
|
|---|
| 28 | {
|
|---|
| 29 | private:
|
|---|
| 30 |
|
|---|
| 31 | MHillas *fHil;
|
|---|
| 32 | MHillasSrc *fHilSrc;
|
|---|
| 33 | MHillasExt *fHilExt;
|
|---|
| 34 | MNewImagePar *fNewPar;
|
|---|
| 35 | MMcEvt *fMcEvt;
|
|---|
| 36 | MHadronness *fHadronness; //! output container for hadronness
|
|---|
| 37 | MMySuperCuts *fSuper; // container for supercut parameters
|
|---|
| 38 | MPointingPos *fPointingPos;
|
|---|
| 39 | MNewImagePar *fNewImagePar;
|
|---|
| 40 | MHillasExt *fHillasExt;
|
|---|
| 41 |
|
|---|
| 42 | TString fHadronnessName; // name of container to store hadronness
|
|---|
| 43 | TString fHilName;
|
|---|
| 44 | TString fHilSrcName;
|
|---|
| 45 | TString fHilExtName;
|
|---|
| 46 | TString fNewParName;
|
|---|
| 47 | TString fSuperName; // name of container for supercut parameters
|
|---|
| 48 |
|
|---|
| 49 | Double_t fMm2Deg; //!
|
|---|
| 50 |
|
|---|
| 51 | Int_t fMap[11]; //!
|
|---|
| 52 | MHMatrix *fMatrix; //!
|
|---|
| 53 |
|
|---|
| 54 | Bool_t fNoDistCut;
|
|---|
| 55 | Float_t fSizeCutLow,fSizeCutUp ;
|
|---|
| 56 |
|
|---|
| 57 |
|
|---|
| 58 | Int_t PreProcess(MParList *pList);
|
|---|
| 59 | Int_t Process();
|
|---|
| 60 |
|
|---|
| 61 | Double_t GetVal(Int_t i) const;
|
|---|
| 62 |
|
|---|
| 63 | Double_t CtsMCut(const Double_t* a, Double_t ls, Double_t ls2, Double_t ct, Double_t dd2) const;
|
|---|
| 64 |
|
|---|
| 65 |
|
|---|
| 66 | public:
|
|---|
| 67 |
|
|---|
| 68 | MMySuperCutsCalc(const char *hilname="MHillas",
|
|---|
| 69 | const char *hilsrcname="MHillasSrc",
|
|---|
| 70 | const char *name=NULL, const char *title=NULL);
|
|---|
| 71 |
|
|---|
| 72 | void SetHadronnessName(const TString name) { fHadronnessName = name; }
|
|---|
| 73 | TString GetHadronnessName() const { return fHadronnessName; }
|
|---|
| 74 |
|
|---|
| 75 | void InitMapping(MHMatrix *mat);
|
|---|
| 76 | void StopMapping() { InitMapping(NULL); }
|
|---|
| 77 |
|
|---|
| 78 |
|
|---|
| 79 | /* Double_t Calc(MMySuperCuts* super, MHillas* hillas, MHillasSrc* hillasSrc, MHadronness* hadronness); */
|
|---|
| 80 |
|
|---|
| 81 | Bool_t CalcDistCut(MHillasSrc* hillasSrc);
|
|---|
| 82 |
|
|---|
| 83 |
|
|---|
| 84 | void SetNoDistCut(Bool_t flag) { fNoDistCut = flag; }
|
|---|
| 85 | void SetSizeCuts(Float_t sizeLow, Float_t sizeUp)
|
|---|
| 86 | {
|
|---|
| 87 | fSizeCutLow = sizeLow;
|
|---|
| 88 | fSizeCutUp = sizeUp;
|
|---|
| 89 | }
|
|---|
| 90 |
|
|---|
| 91 |
|
|---|
| 92 | ClassDef(MMySuperCutsCalc, 0) // A class to evaluate the Supercuts
|
|---|
| 93 | };
|
|---|
| 94 |
|
|---|
| 95 | #endif
|
|---|
| 96 |
|
|---|
| 97 |
|
|---|
| 98 |
|
|---|
| 99 |
|
|---|
| 100 |
|
|---|
| 101 |
|
|---|
| 102 |
|
|---|
| 103 |
|
|---|
| 104 |
|
|---|
| 105 |
|
|---|
| 106 |
|
|---|
| 107 |
|
|---|
| 108 |
|
|---|
| 109 |
|
|---|
| 110 |
|
|---|
| 111 |
|
|---|
| 112 |
|
|---|
| 113 |
|
|---|
| 114 |
|
|---|
