| Line | |
|---|
| 1 | #ifndef MARS_MCT1SupercutsCalc
|
|---|
| 2 | #define MARS_MCT1SupercutsCalc
|
|---|
| 3 |
|
|---|
| 4 | /////////////////////////////////////////////////////////////////////////////
|
|---|
| 5 | // //
|
|---|
| 6 | // MCT1SupercutsCalc //
|
|---|
| 7 | // //
|
|---|
| 8 | /////////////////////////////////////////////////////////////////////////////
|
|---|
| 9 |
|
|---|
| 10 | #ifndef MARS_MFilter
|
|---|
| 11 | #include "MFilter.h"
|
|---|
| 12 | #endif
|
|---|
| 13 |
|
|---|
| 14 | class MParList;
|
|---|
| 15 | class MHillas;
|
|---|
| 16 | class MHillasSrc;
|
|---|
| 17 | class MMcEvt;
|
|---|
| 18 | class MCerPhotEvt;
|
|---|
| 19 | class MGeomCam;
|
|---|
| 20 | class MHadronness;
|
|---|
| 21 |
|
|---|
| 22 |
|
|---|
| 23 | class MCT1SupercutsCalc : public MTask
|
|---|
| 24 | {
|
|---|
| 25 | private:
|
|---|
| 26 | MHillas *fHil;
|
|---|
| 27 | MHillasSrc *fHilSrc;
|
|---|
| 28 | MMcEvt *fMcEvt;
|
|---|
| 29 | MGeomCam *fCam;
|
|---|
| 30 | MHadronness *fHadronness; //! output container for hadronness
|
|---|
| 31 |
|
|---|
| 32 | TString fHadronnessName; // name of container to store hadronness
|
|---|
| 33 | TString fHilName;
|
|---|
| 34 | TString fHilSrcName;
|
|---|
| 35 |
|
|---|
| 36 | Double_t fMm2Deg;
|
|---|
| 37 |
|
|---|
| 38 | //---------------------------------
|
|---|
| 39 | // cut parameters
|
|---|
| 40 | Double_t fLengthUp[8];
|
|---|
| 41 | Double_t fWidthUp[8];
|
|---|
| 42 | Double_t fDistUp[8];
|
|---|
| 43 | Double_t fLengthLo[8];
|
|---|
| 44 | Double_t fWidthLo[8];
|
|---|
| 45 | Double_t fDistLo[8];
|
|---|
| 46 | Double_t fAsymUp[8];
|
|---|
| 47 | Double_t fAsymLo[8];
|
|---|
| 48 | Double_t fAlphaUp[8];
|
|---|
| 49 | //---------------------------------
|
|---|
| 50 |
|
|---|
| 51 | void InitParams();
|
|---|
| 52 |
|
|---|
| 53 | public:
|
|---|
| 54 | MCT1SupercutsCalc(const char *hilname="MHillas",
|
|---|
| 55 | const char *hilsrcname="MHillasSrc",
|
|---|
| 56 | const char *name=NULL, const char *title=NULL);
|
|---|
| 57 | ~MCT1SupercutsCalc();
|
|---|
| 58 |
|
|---|
| 59 | Double_t CtsMCut(Double_t *a, Double_t ls, Double_t ct,
|
|---|
| 60 | Double_t ls2, Double_t dd2);
|
|---|
| 61 |
|
|---|
| 62 | void SetHadronnessName(const TString name) { fHadronnessName = name; }
|
|---|
| 63 | TString GetHadronnessName() const { return fHadronnessName; }
|
|---|
| 64 |
|
|---|
| 65 | Bool_t PreProcess(MParList *pList);
|
|---|
| 66 | Bool_t Process();
|
|---|
| 67 |
|
|---|
| 68 | ClassDef(MCT1SupercutsCalc, 0) // A class to evaluate the Supercuts
|
|---|
| 69 | };
|
|---|
| 70 |
|
|---|
| 71 | #endif
|
|---|
| 72 |
|
|---|
| 73 |
|
|---|
| 74 |
|
|---|
| 75 |
|
|---|
| 76 |
|
|---|
Note:
See
TracBrowser
for help on using the repository browser.
