source: branches/Corsika7500Compatibility/msim/MPhotonEvent.h@ 19851

Last change on this file since 19851 was 18526, checked in by dbaack, 8 years ago
Bugfixes and gcc-5 compability changes
File size: 3.1 KB
Line 
1#ifndef MARS_MPhotonEvent
2#define MARS_MPhotonEvent
3
4#ifndef MARS_MParContainer
5#include "MParContainer.h"
6#endif
7
8#ifndef ROOT_TClonesArray
9#include <TClonesArray.h>
10#endif
11
12#include <iosfwd>
13
14using namespace std;
15
16class MPhotonData;
17class MCorsikaRunHeader;
18class MCorsikaFormat;
19
20class MPhotonEvent : public MParContainer
21{
22private:
23 TClonesArray fData;
24
25public:
26 MPhotonEvent(const char *name=NULL, const char *title=NULL);
27
28 void Sort(Bool_t force=kFALSE);
29 Bool_t IsSorted() const { return fData.IsSorted(); }
30
31 // Getter/Setter
32 Int_t GetNumPhotons() const { return fData.GetEntriesFast(); }
33 Int_t GetNumExternal() const;
34
35 Float_t GetTimeFirst() const;
36 Float_t GetTimeLast() const;
37 Double_t GetTimeMedianDev() const;
38
39 Double_t GetMeanX() const;
40 Double_t GetMeanY() const;
41 Double_t GetMeanT() const;
42
43 TClonesArray &GetArray() { return fData; }
44 const TClonesArray &GetArray() const { return fData; }
45
46 Bool_t AsciiWrite(ostream &out) const;
47
48 MPhotonData &Add(Int_t n);
49 MPhotonData &Add();
50
51 MPhotonData *GetFirst() const;
52 MPhotonData *GetLast() const;
53
54 MPhotonData &operator[](UInt_t idx);
55 const MPhotonData &operator[](UInt_t idx) const;
56
57 Int_t Shrink(Int_t n);
58 void Resize(Int_t n);
59
60 void AddXY(Double_t x, Double_t y);
61 void SimWavelength(Float_t wmin, Float_t wmax);
62
63 // I/O
64 Int_t ReadEventIoEvt(MCorsikaFormat *fInFormat);
65 Int_t ReadEventIoEvtCompact(MCorsikaFormat *fInFormat);
66 Int_t ReadCorsikaEvt(Float_t * data, Int_t numEvents, Int_t arrayIdx);
67
68 // TObject
69 void Paint(Option_t *o="");
70 void Print(Option_t * = NULL) const;
71 //void Clear(Option_t * = NULL);
72
73 void Reset();
74
75 ClassDef(MPhotonEvent, 1) //Container to store the raw Event Data
76};
77
78// FIXME: Should we merge this into MPhotonEvent?
79class MPhotonStatistics : public MParContainer
80{
81private:
82 Float_t fTimeFirst; //! Start of (simulated) sampling window
83 Float_t fTimeLast; //! Start of (simulated) sampling window
84
85 Float_t fLength; // Time between first and last photon
86 Float_t fTimeMedDev; // Median deviation
87
88// Float_t fOffset;
89// Float_t fWindow;
90
91 Int_t fMaxIndex; //!
92
93public:
94 MPhotonStatistics(const char *name=NULL, const char *title=NULL) : fMaxIndex(-1)
95 {
96 fName = name ? name : "MPhotonStatistics";
97 fTitle = title ? title : "Corsika Event Data Information";
98 }
99
100 void SetTime(Float_t first, Float_t last) { fTimeFirst=first; fTimeLast=last; }
101 void SetLength(Float_t len) { fLength=len; }
102 void SetMaxIndex(UInt_t idx) { fMaxIndex=idx; }
103 void SetTimeMedDev(Float_t dev) { fTimeMedDev=dev; }
104
105// Float_t GetRawTimeFirst() const { return fTimeFirst; }
106// Float_t GetRawTimeLast() const { return fTimeLast; }
107
108 Float_t GetTimeFirst() const { return fTimeFirst; }
109 Float_t GetTimeLast() const { return fTimeLast; }
110
111 Float_t GetLength() const { return fLength; }
112 Float_t GetTimeMedDev() const { return fTimeMedDev; }
113
114 Int_t GetMaxIndex() const { return fMaxIndex; }
115
116// Bool_t IsValid() const { return fTimeLast>=fTimeFirst; }
117
118 ClassDef(MPhotonStatistics, 1)
119};
120#endif
Note: See TracBrowser for help on using the repository browser.