1 | /* ======================================================================== *\
|
---|
2 | !
|
---|
3 | ! *
|
---|
4 | ! * This file is part of MARS, the MAGIC Analysis and Reconstruction
|
---|
5 | ! * Software. It is distributed to you in the hope that it can be a useful
|
---|
6 | ! * and timesaving tool in analysing Data of imaging Cerenkov telescopes.
|
---|
7 | ! * It is distributed WITHOUT ANY WARRANTY.
|
---|
8 | ! *
|
---|
9 | ! * Permission to use, copy, modify and distribute this software and its
|
---|
10 | ! * documentation for any purpose is hereby granted without fee,
|
---|
11 | ! * provided that the above copyright notice appear in all copies and
|
---|
12 | ! * that both that copyright notice and this permission notice appear
|
---|
13 | ! * in supporting documentation. It is provided "as is" without express
|
---|
14 | ! * or implied warranty.
|
---|
15 | ! *
|
---|
16 | !
|
---|
17 | !
|
---|
18 | ! Author(s): Thomas Bretz, 06/2001 <mailto:tbretz@astro.uni-wuerzburg.de>
|
---|
19 | !
|
---|
20 | ! Copyright: MAGIC Software Development, 2000-2003
|
---|
21 | !
|
---|
22 | !
|
---|
23 | \* ======================================================================== */
|
---|
24 |
|
---|
25 | /////////////////////////////////////////////////////////////////////////////
|
---|
26 | //
|
---|
27 | // MWriteAsciiFile
|
---|
28 | //
|
---|
29 | // If you want to store a single container into an Ascii file you have
|
---|
30 | // to use this class. You must know the name of the file you wanne write
|
---|
31 | // (you should know it) and the name of the container you want to write.
|
---|
32 | // This can be the name of the class or a given name, which identifies
|
---|
33 | // the container in a parameter container list (MParList).
|
---|
34 | // The container is written to the ascii file if its ReadyToSave flag is
|
---|
35 | // set (MParContainer)
|
---|
36 | //
|
---|
37 | // You can write more than one container in one line of the file, see
|
---|
38 | // AddContainer.
|
---|
39 | //
|
---|
40 | // You can also write single data members of a container (like fWidth
|
---|
41 | // of MHillas). For more details see AddContainer. Make sure, that a
|
---|
42 | // getter method for the data member exist. The name of the method
|
---|
43 | // must be the same than the data member itself, but the f must be
|
---|
44 | // replaced by a Get.
|
---|
45 | //
|
---|
46 | /////////////////////////////////////////////////////////////////////////////
|
---|
47 | #include "MWriteAsciiFile.h"
|
---|
48 |
|
---|
49 | #include <fstream>
|
---|
50 |
|
---|
51 | #include <TMethodCall.h> // TMethodCall, AsciiWrite
|
---|
52 |
|
---|
53 | #include "MIter.h"
|
---|
54 |
|
---|
55 | #include "MDataList.h" // MDataList
|
---|
56 | #include "MDataChain.h" // MDataChain
|
---|
57 | #include "MDataValue.h" // MDataValue
|
---|
58 | #include "MDataMember.h" // MDataMember
|
---|
59 |
|
---|
60 | #include "MLog.h"
|
---|
61 | #include "MLogManip.h"
|
---|
62 |
|
---|
63 | #include "MParList.h"
|
---|
64 |
|
---|
65 | ClassImp(MWriteAsciiFile);
|
---|
66 |
|
---|
67 | using namespace std;
|
---|
68 |
|
---|
69 | // --------------------------------------------------------------------------
|
---|
70 | //
|
---|
71 | // Init. Replaces the same code used in all constructors.
|
---|
72 | //
|
---|
73 | void MWriteAsciiFile::Init(const char *filename, const char *name, const char *title)
|
---|
74 | {
|
---|
75 | fName = name ? name : "MWriteAsciiFile";
|
---|
76 | fTitle = title ? title : "Task to write one container to an ascii file";
|
---|
77 |
|
---|
78 | fNameFile = filename;
|
---|
79 |
|
---|
80 | fOut = new ofstream(fNameFile);
|
---|
81 | }
|
---|
82 |
|
---|
83 | // --------------------------------------------------------------------------
|
---|
84 | //
|
---|
85 | // Specify the name of the ascii output file 'filename' and the name
|
---|
86 | // of the container you want to write. (Normally this is the name
|
---|
87 | // of the class (eg. MHillas) but it can also be a different name which
|
---|
88 | // identifies the container in the parameter list.
|
---|
89 | // Because you cannot write more than one container there is no Add-function
|
---|
90 | // like in MWriteRootFile.
|
---|
91 | //
|
---|
92 | // For Example: MWriteAsciiFile("file.txt", "MHillas");
|
---|
93 | //
|
---|
94 | MWriteAsciiFile::MWriteAsciiFile(const char *filename, const char *contname,
|
---|
95 | const char *name, const char *title)
|
---|
96 | {
|
---|
97 | Init(filename, name, title);
|
---|
98 |
|
---|
99 | if (contname)
|
---|
100 | AddColumns(contname);
|
---|
101 | }
|
---|
102 |
|
---|
103 | // --------------------------------------------------------------------------
|
---|
104 | //
|
---|
105 | // Specify a the name of the ascii output file 'filename' and a pointer to
|
---|
106 | // the container you want to write.
|
---|
107 | // Because you cannot write more than one container there is no Add-function
|
---|
108 | // like in MWriteRootFile.
|
---|
109 | //
|
---|
110 | // For Example: MHillas hillas;
|
---|
111 | // MWriteAsciiFile("file.txt", &hillas);
|
---|
112 | //
|
---|
113 | //
|
---|
114 | MWriteAsciiFile::MWriteAsciiFile(const char *filename, MParContainer *cont,
|
---|
115 | const char *name, const char *title)
|
---|
116 | {
|
---|
117 | Init(filename, name, title);
|
---|
118 |
|
---|
119 | if (cont)
|
---|
120 | AddColumns(cont);
|
---|
121 | }
|
---|
122 |
|
---|
123 | // --------------------------------------------------------------------------
|
---|
124 | //
|
---|
125 | // Destructor. Delete the output file if necessary (it is closed
|
---|
126 | // automatically by its destructor.
|
---|
127 | //
|
---|
128 | MWriteAsciiFile::~MWriteAsciiFile()
|
---|
129 | {
|
---|
130 | fAutoDel.SetOwner();
|
---|
131 | delete fOut;
|
---|
132 | }
|
---|
133 |
|
---|
134 | // --------------------------------------------------------------------------
|
---|
135 | //
|
---|
136 | // Return open state of the file
|
---|
137 | //
|
---|
138 | Bool_t MWriteAsciiFile::IsFileOpen() const
|
---|
139 | {
|
---|
140 | return (bool)(*fOut);
|
---|
141 | }
|
---|
142 |
|
---|
143 | // --------------------------------------------------------------------------
|
---|
144 | //
|
---|
145 | // Tries to get all containers from the ParList which were given by name
|
---|
146 | // adds them to the list of pointers to the container which should be
|
---|
147 | // written to the ascii file.
|
---|
148 | //
|
---|
149 | Bool_t MWriteAsciiFile::GetContainer(MParList *plist)
|
---|
150 | {
|
---|
151 | TObject *obj=NULL;
|
---|
152 |
|
---|
153 | TIter Next(&fList);
|
---|
154 | while ((obj=Next()))
|
---|
155 | {
|
---|
156 | //
|
---|
157 | // MData is the highest class in the inheritance tree
|
---|
158 | //
|
---|
159 | if (obj->InheritsFrom(MData::Class()))
|
---|
160 | {
|
---|
161 | if (!((MData*)obj)->PreProcess(plist))
|
---|
162 | return kFALSE;
|
---|
163 | continue;
|
---|
164 | }
|
---|
165 |
|
---|
166 | //
|
---|
167 | // MParContainer is the next class in the inheritance tree
|
---|
168 | //
|
---|
169 | if (obj->InheritsFrom(MParContainer::Class()))
|
---|
170 | continue;
|
---|
171 |
|
---|
172 | //
|
---|
173 | // It is neither a MData nor a MParContainer, it must be a TNamed
|
---|
174 | //
|
---|
175 | TObject *o = plist->FindObject(obj->GetName());
|
---|
176 | if (!o)
|
---|
177 | return kFALSE;
|
---|
178 |
|
---|
179 | fList[fList.IndexOf(obj)] = o;
|
---|
180 | }
|
---|
181 | return kTRUE;
|
---|
182 | }
|
---|
183 |
|
---|
184 | // --------------------------------------------------------------------------
|
---|
185 | //
|
---|
186 | // Check if the containers are ready for writing. If so write them.
|
---|
187 | // The containers are written in one line, one after each other.
|
---|
188 | // If not all containers are written (because of the IsReadyToSave-flag)
|
---|
189 | // a warning message is print.
|
---|
190 | //
|
---|
191 | Bool_t MWriteAsciiFile::CheckAndWrite()
|
---|
192 | {
|
---|
193 | MParContainer *obj = NULL;
|
---|
194 |
|
---|
195 | //
|
---|
196 | // Check for the Write flag
|
---|
197 | //
|
---|
198 | Bool_t write = kFALSE;
|
---|
199 | MIter Next(&fList);
|
---|
200 | while ((obj=Next()))
|
---|
201 | if (obj->IsReadyToSave())
|
---|
202 | {
|
---|
203 | write = kTRUE;
|
---|
204 | break;
|
---|
205 | }
|
---|
206 |
|
---|
207 | //
|
---|
208 | // Do not write if not at least one ReadyToSave-flag set
|
---|
209 | //
|
---|
210 | if (!write)
|
---|
211 | return kTRUE;
|
---|
212 |
|
---|
213 | Bool_t written = kFALSE;
|
---|
214 | Int_t num = fList.GetEntries();
|
---|
215 |
|
---|
216 | Next.Reset();
|
---|
217 | while ((obj=Next()))
|
---|
218 | {
|
---|
219 | if (written)
|
---|
220 | *fOut << " ";
|
---|
221 |
|
---|
222 | //
|
---|
223 | // Write container to file
|
---|
224 | //
|
---|
225 | if (!obj->AsciiWrite(*fOut))
|
---|
226 | continue;
|
---|
227 |
|
---|
228 | //
|
---|
229 | // Store that at least one container was written
|
---|
230 | //
|
---|
231 | written = kTRUE;
|
---|
232 |
|
---|
233 | num--;
|
---|
234 | }
|
---|
235 |
|
---|
236 | if (written)
|
---|
237 | {
|
---|
238 | *fOut << endl;
|
---|
239 |
|
---|
240 | if (num!=0)
|
---|
241 | *fLog << warn << "WARNING - Number of objects written mismatch!" << endl;
|
---|
242 | }
|
---|
243 | return kTRUE;
|
---|
244 | }
|
---|
245 |
|
---|
246 | // --------------------------------------------------------------------------
|
---|
247 | //
|
---|
248 | // Add a rule to be written as a column to the ascii file.
|
---|
249 | // For more information about rules see MDataChain.
|
---|
250 | //
|
---|
251 | // eg: MWriteAsciiFile::AddColumn("log10(MHillas.fEnergy)/2")
|
---|
252 | //
|
---|
253 | void MWriteAsciiFile::AddColumn(const TString rule)
|
---|
254 | {
|
---|
255 | MDataChain *chain = new MDataChain(rule);
|
---|
256 | fList.Add(chain);
|
---|
257 | }
|
---|
258 |
|
---|
259 | // --------------------------------------------------------------------------
|
---|
260 | //
|
---|
261 | // Add another column to be written to the ascii file. The columns will be
|
---|
262 | // output one after each other in one line.
|
---|
263 | // Specify the name of the data member to be written (eg fWidth) and
|
---|
264 | // a possible scale factor (eg. to transform millimeters to degrees)
|
---|
265 | //
|
---|
266 | // eg:
|
---|
267 | // MMcEvt evt;
|
---|
268 | // MWriteAsciiFile::AddColumn(&evt, "fImpact", 0.01);
|
---|
269 | //
|
---|
270 | void MWriteAsciiFile::AddColumn(MParContainer *cont, const TString member, Double_t scale)
|
---|
271 | {
|
---|
272 | MData *data = new MDataMember(cont, member);
|
---|
273 |
|
---|
274 | if (scale!=1)
|
---|
275 | {
|
---|
276 | MDataList *list = new MDataList('*');
|
---|
277 | MDataValue *val = new MDataValue(scale);
|
---|
278 |
|
---|
279 | list->SetOwner();
|
---|
280 | list->AddToList(data);
|
---|
281 | list->AddToList(val);
|
---|
282 |
|
---|
283 | data = list;
|
---|
284 | }
|
---|
285 | fList.Add(data);
|
---|
286 | fAutoDel.Add(data);
|
---|
287 | }
|
---|
288 |
|
---|
289 | // --------------------------------------------------------------------------
|
---|
290 | //
|
---|
291 | // Add another container (by name) to be written to the ascii file.
|
---|
292 | // The container will be output one after each other in one line.
|
---|
293 | // The output will be done either by MParContainer::AsciiWrite or
|
---|
294 | // by the corresponding overloaded function.
|
---|
295 | //
|
---|
296 | // eg: MWriteAsciiFile::AddColumns("MMcEvt")
|
---|
297 | //
|
---|
298 | void MWriteAsciiFile::AddColumns(const TString cont)
|
---|
299 | {
|
---|
300 | TNamed *name = new TNamed(cont, cont);
|
---|
301 | fList.Add(name);
|
---|
302 | fAutoDel.Add(name);
|
---|
303 | }
|
---|
304 |
|
---|
305 | // --------------------------------------------------------------------------
|
---|
306 | //
|
---|
307 | // Add another container (by pointer) to be written to the ascii file.
|
---|
308 | // The container will be output one after each other in one line.
|
---|
309 | // The output will be done either by MParContainer::AsciiWrite or
|
---|
310 | // by the corresponding overloaded function.
|
---|
311 | //
|
---|
312 | // eg:
|
---|
313 | // MMcEvt evt;
|
---|
314 | // MWriteAsciiFile::AddColumns(&evt);
|
---|
315 | //
|
---|
316 | void MWriteAsciiFile::AddColumns(MParContainer *cont)
|
---|
317 | {
|
---|
318 | fList.Add(cont);
|
---|
319 | }
|
---|