| 1 | #ifndef MARS_MMuonSearchPar
|
|---|
| 2 | #define MARS_MMuonSearchPar
|
|---|
| 3 |
|
|---|
| 4 | #ifndef MARS_MParContainer
|
|---|
| 5 | #include "MParContainer.h"
|
|---|
| 6 | #endif
|
|---|
| 7 |
|
|---|
| 8 | #ifndef ROOT_TH1
|
|---|
| 9 | #include <TH1.h>
|
|---|
| 10 | #endif
|
|---|
| 11 |
|
|---|
| 12 | #ifndef ROOT_TArrayI
|
|---|
| 13 | #include <TArrayI.h>
|
|---|
| 14 | #endif
|
|---|
| 15 |
|
|---|
| 16 | class MHillas;
|
|---|
| 17 | class MGeomCam;
|
|---|
| 18 | class MCerPhotEvt;
|
|---|
| 19 |
|
|---|
| 20 | class MMuonSearchPar : public MParContainer
|
|---|
| 21 | {
|
|---|
| 22 | private:
|
|---|
| 23 | Float_t fRad; // An estimated radius of the muon ring [mm]
|
|---|
| 24 | Float_t fDev; // The standard deviation from the estimated ring [mm]
|
|---|
| 25 | Float_t fCenX; // An estimated center position in X of the muon ring [mm]
|
|---|
| 26 | Float_t fCenY; // An estimated center position in Y of the muon ring [mm]
|
|---|
| 27 | Bool_t fNoMuon; // if the radius is estimated above 600 mm (2 deg.), assumed it's not muon.
|
|---|
| 28 |
|
|---|
| 29 | public:
|
|---|
| 30 | MMuonSearchPar(const char *name=NULL, const char *title=NULL);
|
|---|
| 31 |
|
|---|
| 32 | void Reset();
|
|---|
| 33 |
|
|---|
| 34 | Float_t GetRad() const { return fRad; }
|
|---|
| 35 | Float_t GetDev() const { return fDev; }
|
|---|
| 36 | Float_t GetCenX() const { return fCenX; }
|
|---|
| 37 | Float_t GetCenY() const { return fCenY; }
|
|---|
| 38 | Bool_t IsNoMuon() const { return fNoMuon; }
|
|---|
| 39 |
|
|---|
| 40 | void Print(Option_t *opt=NULL) const;
|
|---|
| 41 | void Print(const MGeomCam &geom, Option_t *opt=NULL) const;
|
|---|
| 42 | void CalcTempCenter(const MHillas &hillas, Float_t *xtmp1,
|
|---|
| 43 | Float_t *ytmp1, Float_t *xtmp2, Float_t *ytmp2);
|
|---|
| 44 | Bool_t CalcRadius(const MGeomCam &geom, const MCerPhotEvt &evt, Float_t x,
|
|---|
| 45 | Float_t y, Float_t *r, Float_t *sigma);
|
|---|
| 46 | void CalcMinimumDev(const MGeomCam &geom, const MCerPhotEvt &evt,
|
|---|
| 47 | Float_t x, Float_t y, Float_t xcog,
|
|---|
| 48 | Float_t ycog, Float_t sigma, Float_t *opt_rad,
|
|---|
| 49 | Float_t *new_sigma, Float_t *newx, Float_t *newy);
|
|---|
| 50 | void Calc(const MGeomCam &geom, const MCerPhotEvt &evt,
|
|---|
| 51 | const MHillas &hillas);
|
|---|
| 52 |
|
|---|
| 53 | ClassDef(MMuonSearchPar, 1) // Container to hold muon search parameters
|
|---|
| 54 | };
|
|---|
| 55 |
|
|---|
| 56 | #endif
|
|---|
