1 | #ifndef MARS_MSequence
|
---|
2 | #define MARS_MSequence
|
---|
3 |
|
---|
4 | #ifndef ROOT_TArrayI
|
---|
5 | #include <TArrayI.h>
|
---|
6 | #endif
|
---|
7 |
|
---|
8 | #ifndef MARS_MTime
|
---|
9 | #include "MTime.h"
|
---|
10 | #endif
|
---|
11 |
|
---|
12 | class MDirIter;
|
---|
13 |
|
---|
14 | class MSequence : public MParContainer
|
---|
15 | {
|
---|
16 | friend class MSequenceSQL;
|
---|
17 |
|
---|
18 | public:
|
---|
19 | enum LightCondition_t { kNA, kNoMoon, kTwilight, kMoon, kDay };
|
---|
20 |
|
---|
21 | enum FileExt_t {
|
---|
22 | kRaw =BIT(0),
|
---|
23 | kRoot =BIT(1),
|
---|
24 | kReport=BIT(2),
|
---|
25 | };
|
---|
26 | enum RunType_t {
|
---|
27 | kPed =BIT(3),
|
---|
28 | kCal =BIT(4),
|
---|
29 | kDat =BIT(5),
|
---|
30 | kExcl =BIT(6),
|
---|
31 | kAll =kPed|kDat|kCal,
|
---|
32 | kImage =BIT(7),
|
---|
33 | kCalib =BIT(8)
|
---|
34 | };
|
---|
35 |
|
---|
36 | enum { kRunType = kPed|kCal|kDat|kImage|kCalib };
|
---|
37 | enum { kFileType = kRaw|kRoot|kReport };
|
---|
38 | enum { kFilePaths = kFileType|kImage|kCalib };
|
---|
39 |
|
---|
40 | enum FileType_t {
|
---|
41 | kRawDat = kRaw |kDat,
|
---|
42 | kRawPed = kRaw |kPed,
|
---|
43 | kRawCal = kRaw |kCal,
|
---|
44 | kRawAll = kRaw |kAll,
|
---|
45 | kRootDat = kRoot |kDat,
|
---|
46 | kRootPed = kRoot |kPed,
|
---|
47 | kRootCal = kRoot |kCal,
|
---|
48 | kRootAll = kRoot |kAll,
|
---|
49 | kReportPed = kReport|kPed,
|
---|
50 | kReportCal = kReport|kCal,
|
---|
51 | kReportDat = kReport|kDat,
|
---|
52 | kReportAll = kReport|kAll,
|
---|
53 | kCalibrated = kRoot |kCalib,
|
---|
54 | kImages = kRoot |kImage,
|
---|
55 | };
|
---|
56 |
|
---|
57 | private:
|
---|
58 | TString fFileName; // Expanded file name
|
---|
59 | TString fDataPath; // Path to data files
|
---|
60 |
|
---|
61 | UShort_t fTelescope; // Telescope number
|
---|
62 | UInt_t fSequence; // Sequence number
|
---|
63 |
|
---|
64 | MTime fStart; // Start time of sequence
|
---|
65 |
|
---|
66 | UInt_t fLastRun; // Last run in sequence (necessary?)
|
---|
67 | UInt_t fNumEvents; // Number of events in sequence
|
---|
68 |
|
---|
69 | UInt_t fPeriod; // Observation period of the sequence
|
---|
70 | MTime fNight; // Night (day of sunrise) of the sequence
|
---|
71 |
|
---|
72 | LightCondition_t fLightCondition;
|
---|
73 |
|
---|
74 | TString fProject;
|
---|
75 | TString fSource;
|
---|
76 | TString fTriggerTable;
|
---|
77 | TString fHvSettings;
|
---|
78 | TString fComment;
|
---|
79 |
|
---|
80 | TArrayI fRuns; // Run numbers
|
---|
81 | TArrayI fRunsSub; // Sub runs (files)
|
---|
82 |
|
---|
83 | TArrayI fCalRuns; // Numbers of calibration runs/files
|
---|
84 | TArrayI fCalRunsSub; // Sub runs (files) of calibration runs
|
---|
85 |
|
---|
86 | TArrayI fPedRuns; // Numbers of pedestal runs/files
|
---|
87 | TArrayI fPedRunsSub; // Sub runs (files) of pedestal runs
|
---|
88 |
|
---|
89 | TArrayI fDatRuns; // Numbers of data runs/files
|
---|
90 | TArrayI fDatRunsSub; // Sub runs (files) of data runs
|
---|
91 |
|
---|
92 | TArrayI fExclRuns; // Numbers of excluded runs/files
|
---|
93 | TArrayI fExclRunsSub; // Sub runs (files) of excluded runs
|
---|
94 |
|
---|
95 | Bool_t fMonteCarlo; // Flag for Monte Carlo sequences
|
---|
96 |
|
---|
97 | // Helper for interpretation
|
---|
98 | void AddEntry(Int_t run, Int_t file, TArrayI &data, TArrayI &sub) const;
|
---|
99 | void EvalEntry(const TEnv *env, const TString &prefix, const TString &num, TArrayI &data, TArrayI &sub) const;
|
---|
100 |
|
---|
101 | void Split(TString &runs, TArrayI &data, TArrayI &sub, const TEnv *env=0, const TString prefix="") const;
|
---|
102 | void Split(const TEnv &env, const TString &prefix, const char *name, TArrayI &data, TArrayI &sub) const;
|
---|
103 |
|
---|
104 | LightCondition_t GetLightCondition(const TString &str) const;
|
---|
105 |
|
---|
106 | // Helper for file setup
|
---|
107 | void GetArrays(const TArrayI* &r, const TArrayI* &f, Int_t type) const;
|
---|
108 | Bool_t FindPath(TString &path, FileType_t type) const;
|
---|
109 | TString FindPath(FileType_t type, TString path="") const { FindPath(path, type); return path; }
|
---|
110 | TString FindFile(const TString &d, const TString &f) const;
|
---|
111 | TString GetPathName(TString d, FileType_t type) const;
|
---|
112 | TString GetFileRegexp(UInt_t num, const TArrayI &arr, const TArrayI &sub, FileType_t type) const;
|
---|
113 | TString GetFileName(UInt_t i, const TArrayI &arr, const TArrayI &sub, FileType_t type, const char *path) const;
|
---|
114 | UInt_t GetRuns(MDirIter &iter, const TArrayI &arr, const TArrayI &file, FileType_t type, TString path="") const;
|
---|
115 |
|
---|
116 | // Helper for Print()
|
---|
117 | TString GetNumSequence(Int_t pos, Int_t n, const TArrayI &f) const;
|
---|
118 | TString GetNumSequence(Int_t &pos, const TArrayI &n, const TArrayI &f) const;
|
---|
119 |
|
---|
120 | void PrintRunsClassic(ostream &out, const char *pre, const char *name, const TArrayI &r) const;
|
---|
121 | TString PrintRuns(ostream &out, const char *pre, const char *name, const TArrayI &r, const TArrayI &f) const;
|
---|
122 |
|
---|
123 | // General helper
|
---|
124 | Bool_t HasSubRuns() const { return fRunsSub.GetSize()!=0 || fDatRunsSub.GetSize()!=0 || fPedRunsSub.GetSize()!=0 || fCalRunsSub.GetSize()!=0 || fExclRunsSub.GetSize()!=0; }
|
---|
125 | Bool_t IsContained(const TArrayI &arr, const TArrayI &sub, UInt_t num, UInt_t file) const;
|
---|
126 |
|
---|
127 | // Some helpers to handle the arrays
|
---|
128 | static Int_t SortArraySub(Int_t p, Int_t n, Int_t *arr1, Int_t *arr2);
|
---|
129 | static void SortArrays(TArrayI &arr1, TArrayI &arr2);
|
---|
130 | static Int_t GetSubArray(Int_t p, Int_t n, Int_t *arr1);
|
---|
131 |
|
---|
132 | public:
|
---|
133 | MSequence() : fTelescope(1), fSequence((UInt_t)-1), fLastRun((UInt_t)-1),
|
---|
134 | fNumEvents((UInt_t)-1), fPeriod((UInt_t)-1), fLightCondition(kNA), fMonteCarlo(kFALSE)
|
---|
135 | {
|
---|
136 | fName = "MSequence";
|
---|
137 | fTitle = "Sequence file";
|
---|
138 | }
|
---|
139 | MSequence(const char *fname, const char *path="", UInt_t id=(UInt_t)-1);
|
---|
140 | MSequence(const char *fname, UInt_t seq, UShort_t tel=0);
|
---|
141 | MSequence(const MSequence &s) : MParContainer(s),
|
---|
142 | fFileName(s.fFileName), fDataPath(s.fDataPath),
|
---|
143 | fTelescope(s.fTelescope),
|
---|
144 | fSequence(s.fSequence), fStart(s.fStart),
|
---|
145 | fLastRun(s.fLastRun), fNumEvents(s.fNumEvents), fPeriod(s.fPeriod),
|
---|
146 | fNight(s.fNight), fLightCondition(s.fLightCondition),
|
---|
147 | fProject(s.fProject), fSource(s.fSource),
|
---|
148 | fTriggerTable(s.fTriggerTable), fHvSettings(s.fHvSettings),
|
---|
149 | fRuns(s.fRuns), fRunsSub(s.fRunsSub), fCalRuns(s.fCalRuns),
|
---|
150 | fCalRunsSub(s.fCalRunsSub), fPedRuns(s.fPedRuns),
|
---|
151 | fPedRunsSub(s.fPedRunsSub), fDatRuns(s.fDatRuns),
|
---|
152 | fDatRunsSub(s.fDatRunsSub), fExclRuns(s.fExclRuns),
|
---|
153 | fExclRunsSub(s.fExclRunsSub), fMonteCarlo(s.fMonteCarlo) { }
|
---|
154 |
|
---|
155 | Bool_t IsCompatible(const MSequence &s) const;
|
---|
156 | Bool_t IsIdentical(const MSequence &s) const;
|
---|
157 | Bool_t operator==(const MSequence &s) const;
|
---|
158 |
|
---|
159 | // I/O
|
---|
160 | Int_t Read(const char *name=0);
|
---|
161 | Bool_t ReadDatabase(TString sql, UInt_t seq=0, UShort_t tel=0);
|
---|
162 | Bool_t ReadFile(const char *filename, UInt_t id=(UInt_t)-1);
|
---|
163 | Bool_t WriteFile(const char *filename, const Option_t *o) const;
|
---|
164 | Bool_t WriteFile(const char *filename) const { return WriteFile(filename,""); } //*MENU *ARGS={filename=>fBaseName}
|
---|
165 | Bool_t SetDataPath(const char *path="");
|
---|
166 |
|
---|
167 | // TObject
|
---|
168 | void Print(ostream &out, Option_t *o) const;
|
---|
169 | void Print(Option_t *o) const;
|
---|
170 | void Print() const { Print(""); } //*MENU*
|
---|
171 |
|
---|
172 | // Genaral interface
|
---|
173 | Bool_t IsValid() const { return fTelescope>0 && fSequence!=(UInt_t)-1; }
|
---|
174 | Bool_t IsMonteCarlo() const { return fMonteCarlo; }
|
---|
175 | Bool_t IsExcluded(UInt_t run, UInt_t file) const { return IsContained(fExclRuns, fExclRunsSub, run, file); }
|
---|
176 | Bool_t IsContained(UInt_t run, UInt_t file) const { return IsContained(fCalRuns, fCalRunsSub, run, file) || IsContained(fPedRuns, fPedRunsSub, run, file) || IsContained(fDatRuns, fDatRunsSub, run, file); }
|
---|
177 |
|
---|
178 | // Setter
|
---|
179 | void SetNight(const char*night);
|
---|
180 |
|
---|
181 | void AddRun(UInt_t run, RunType_t type=kAll) { AddFile(run, 0, type); }
|
---|
182 | void AddRuns(UInt_t run1, UInt_t run2, RunType_t type=kAll) { for (UInt_t i=run1; i<=run2; i++) AddFile(i, 0, type); }
|
---|
183 | void AddFile(UInt_t run, UInt_t file, RunType_t type=kAll);
|
---|
184 | void AddFiles(UInt_t run, UInt_t f1, UInt_t f2, RunType_t type=kAll) { for (UInt_t i=f1; i<=f2; i++) AddFile(run, i, type); }
|
---|
185 |
|
---|
186 | void ExcludeFile(UInt_t run, UInt_t file=0) { /*if (!force || !IsExcluded(run, file))*/ AddFile(run, file, kExcl); }
|
---|
187 | void ExcludeRuns(TString runs);
|
---|
188 |
|
---|
189 | void SetMonteCarlo(Bool_t ismc=kTRUE) { fMonteCarlo=ismc; }
|
---|
190 |
|
---|
191 | // Getter
|
---|
192 | UInt_t GetNumExclRuns() const { return fExclRuns.GetSize(); }
|
---|
193 |
|
---|
194 | UShort_t GetTelescope() const { return fTelescope; }
|
---|
195 | UInt_t GetSequence() const { return fSequence; }
|
---|
196 | UInt_t GetLastRun() const { return fLastRun; }
|
---|
197 | UInt_t GetPeriod() const { return fPeriod; }
|
---|
198 | Bool_t HasMoon() const { return fLightCondition==kMoon; }
|
---|
199 |
|
---|
200 | Int_t GetFile(UInt_t i, RunType_t type, UInt_t &run, UInt_t &file) const;
|
---|
201 |
|
---|
202 | LightCondition_t GetLightCondition() const { return fLightCondition; }
|
---|
203 |
|
---|
204 | const MTime &GetStart() const { return fStart; }
|
---|
205 | const MTime &GetNight() const { return fNight; }
|
---|
206 | const TString &GetSource() const { return fSource; }
|
---|
207 |
|
---|
208 | const TString GetExcludedRuns() const;
|
---|
209 |
|
---|
210 | // Filesystem interface
|
---|
211 | UInt_t GetRuns(MDirIter &iter, FileType_t type, const char *path=0) const;
|
---|
212 | TString GetFileRegexp(UInt_t num, FileType_t type) const;
|
---|
213 | TString GetFileName(UInt_t i, FileType_t type, const char *path=0) const;
|
---|
214 | UInt_t GetNumEntries(RunType_t type=kAll) const;
|
---|
215 |
|
---|
216 | // Filesystem getter
|
---|
217 | const char *GetBaseName() const;
|
---|
218 | const char *GetFilePath() const;
|
---|
219 | const TString &GetFileName() const { return fFileName; }
|
---|
220 | const TString &GetDataPath() const { return fDataPath; }
|
---|
221 |
|
---|
222 | const TString GetStandardPath() const { return GetStandardPath(fMonteCarlo); }
|
---|
223 |
|
---|
224 | // Static function for external access
|
---|
225 | static TString GetStandardPath(Bool_t mc) { return mc?"/magic/montecarlo/":"/magic/data/"; }
|
---|
226 | static Bool_t InflateSeq(TString &seq, /*FIXME: TELNUM*/ Bool_t ismc=kFALSE);
|
---|
227 |
|
---|
228 | static TString InflateFileName(const MTime &night, UShort_t tel, Int_t run, Int_t file, Int_t type=kRawAll);
|
---|
229 | static TString InflateFileName(const MTime &night, Int_t run, Int_t type=kRawAll) { return InflateFileName(night, 1, run, 0, type); }
|
---|
230 |
|
---|
231 | static TString InflateFilePath(const MTime &night, UShort_t tel, Int_t run, Int_t file=0, Int_t type=kRawAll);
|
---|
232 | static TString InflateFilePath(const MTime &night, Int_t run, Int_t type=kRawAll) { return InflateFilePath(night, 1, run, 0, type); }
|
---|
233 |
|
---|
234 | static TString InflateFilePath(const MTime &night, Bool_t mc=kFALSE);
|
---|
235 |
|
---|
236 | ClassDef(MSequence, 6) // Describes a sequences, reads and writes sequence files
|
---|
237 | };
|
---|
238 |
|
---|
239 | #endif
|
---|