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 6/2005 <mailto:tbretz@astro.uni-wuerzburg.de>
|
---|
19 | !
|
---|
20 | ! Copyright: MAGIC Software Development, 2000-2005
|
---|
21 | !
|
---|
22 | !
|
---|
23 | \* ======================================================================== */
|
---|
24 |
|
---|
25 | /////////////////////////////////////////////////////////////////////////////
|
---|
26 | //
|
---|
27 | // MParameterCalc
|
---|
28 | //
|
---|
29 | // Task to calculate a parameter stored ina MParametrD by a rule, eg:
|
---|
30 | //
|
---|
31 | // MParameterCalc est;
|
---|
32 | // est.SetRule("0.5 + (1.1*MHillas.fLength) + (2.2*MHillasSrc.fDist) + (3.3*MHillas.fSize) +"
|
---|
33 | // "(4.4*MHillas.fSize*MHillas.fLength) + (5.5*MHillasSrc.fDist*MHillas.fLength)");
|
---|
34 | //
|
---|
35 | // For description of rules, see MDataPhrase.
|
---|
36 | //
|
---|
37 | // Output:
|
---|
38 | // MEnergyEst
|
---|
39 | //
|
---|
40 | /////////////////////////////////////////////////////////////////////////////
|
---|
41 | #include "MParameterCalc.h"
|
---|
42 |
|
---|
43 | #include "MParList.h"
|
---|
44 |
|
---|
45 | #include "MDataPhrase.h"
|
---|
46 | #include "MParameters.h"
|
---|
47 |
|
---|
48 | #include "MLog.h"
|
---|
49 | #include "MLogManip.h"
|
---|
50 |
|
---|
51 | ClassImp(MParameterCalc);
|
---|
52 |
|
---|
53 | using namespace std;
|
---|
54 |
|
---|
55 | // --------------------------------------------------------------------------
|
---|
56 | //
|
---|
57 | // Default constructor. Initialize fData with default rule "MMcEvt.fEnergy"
|
---|
58 | //
|
---|
59 | MParameterCalc::MParameterCalc(const char *def, const char *name, const char *title)
|
---|
60 | : fData(0), fParameter(0), fNameParameter("MParameterD")
|
---|
61 | {
|
---|
62 | fName = name ? name : "MParameterCalc";
|
---|
63 | fTitle = title ? title : "Task to calculate a MParameterD";
|
---|
64 |
|
---|
65 | fData = new MDataPhrase(def);
|
---|
66 | }
|
---|
67 |
|
---|
68 | // --------------------------------------------------------------------------
|
---|
69 | //
|
---|
70 | // delete fData
|
---|
71 | //
|
---|
72 | MParameterCalc::~MParameterCalc()
|
---|
73 | {
|
---|
74 | delete fData;
|
---|
75 | }
|
---|
76 |
|
---|
77 | // --------------------------------------------------------------------------
|
---|
78 | //
|
---|
79 | // Delete fData. Initialize a new MDataPhrase with rule.
|
---|
80 | // Returns if fData->IsValid()
|
---|
81 | //
|
---|
82 | Bool_t MParameterCalc::SetRule(const char *rule)
|
---|
83 | {
|
---|
84 | delete fData;
|
---|
85 | fData = new MDataPhrase(rule);
|
---|
86 |
|
---|
87 | return fData->IsValid();
|
---|
88 | }
|
---|
89 |
|
---|
90 | // --------------------------------------------------------------------------
|
---|
91 | //
|
---|
92 | // Forwards SetVariables to fData to allow optimizations.
|
---|
93 | //
|
---|
94 | void MParameterCalc::SetVariables(const TArrayD &arr)
|
---|
95 | {
|
---|
96 | fData->SetVariables(arr);
|
---|
97 | }
|
---|
98 |
|
---|
99 | // --------------------------------------------------------------------------
|
---|
100 | //
|
---|
101 | // FindCreate "MEnergyEst"
|
---|
102 | // PreProcess fData.
|
---|
103 | //
|
---|
104 | Int_t MParameterCalc::PreProcess(MParList *plist)
|
---|
105 | {
|
---|
106 | memset(fCounter, 0, sizeof(ULong_t)*2);
|
---|
107 |
|
---|
108 | fParameter = (MParameterD*)plist->FindCreateObj("MParameterD", fNameParameter);
|
---|
109 | if (!fParameter)
|
---|
110 | return kFALSE;
|
---|
111 |
|
---|
112 | *fLog << inf << "Rule: " << fData->GetRule() << endl;
|
---|
113 |
|
---|
114 | if (!fData->PreProcess(plist))
|
---|
115 | return kFALSE;
|
---|
116 |
|
---|
117 | return kTRUE;
|
---|
118 | }
|
---|
119 |
|
---|
120 | // --------------------------------------------------------------------------
|
---|
121 | //
|
---|
122 | // Get value from fData and set it to fEnergy. SetReadyToSave for fEnergy.
|
---|
123 | // Return kCONTINUE if value is NaN (Not a Number)
|
---|
124 | //
|
---|
125 | Int_t MParameterCalc::Process()
|
---|
126 | {
|
---|
127 | const Double_t val = fData->GetValue();
|
---|
128 | if (TMath::IsNaN(val))
|
---|
129 | {
|
---|
130 | fCounter[0]++;
|
---|
131 | return kCONTINUE;
|
---|
132 | }
|
---|
133 | if (!TMath::Finite(val))
|
---|
134 | {
|
---|
135 | fCounter[1]++;
|
---|
136 | return kCONTINUE;
|
---|
137 | }
|
---|
138 |
|
---|
139 | fParameter->SetVal(val);
|
---|
140 | fParameter->SetReadyToSave();
|
---|
141 | return kTRUE;
|
---|
142 | }
|
---|
143 |
|
---|
144 | Int_t MParameterCalc::PostProcess()
|
---|
145 | {
|
---|
146 | if (GetNumExecutions()==0)
|
---|
147 | return kTRUE;
|
---|
148 |
|
---|
149 | *fLog << inf << endl;
|
---|
150 | *fLog << GetDescriptor() << " execution statistics:" << endl;
|
---|
151 | PrintSkipped(fCounter[0], "Skipped due to NaN-result (Not a Number)");
|
---|
152 | PrintSkipped(fCounter[1], "Skipped due to inf-result (infinite)");
|
---|
153 | *fLog << endl;
|
---|
154 |
|
---|
155 | return kTRUE;
|
---|
156 | }
|
---|
157 |
|
---|
158 | // --------------------------------------------------------------------------
|
---|
159 | //
|
---|
160 | // Print the rule used for energy estimation
|
---|
161 | //
|
---|
162 | void MParameterCalc::Print(Option_t *o) const
|
---|
163 | {
|
---|
164 | *fLog << all << GetDescriptor() << ":";
|
---|
165 | if (!fData)
|
---|
166 | *fLog << " <n/a>" << endl;
|
---|
167 | else
|
---|
168 | *fLog << endl << fData->GetRule() << endl;
|
---|
169 | }
|
---|
170 |
|
---|
171 | // --------------------------------------------------------------------------
|
---|
172 | //
|
---|
173 | // Check for corresponding entries in resource file and setup filters.
|
---|
174 | // Avoid trailing 0's!
|
---|
175 | //
|
---|
176 | // Example:
|
---|
177 | // test.C:
|
---|
178 | // MParameterCalc est("MyEstimator");
|
---|
179 | //
|
---|
180 | // test.rc:
|
---|
181 | // MyEstimator.Rule: {0} + {1}
|
---|
182 | // MyEstimator.0: log10(MHillas.fSize)
|
---|
183 | // MyEstimator.1: 5.5
|
---|
184 | //
|
---|
185 | // For more details see MDataPhrase::ReadEnv
|
---|
186 | //
|
---|
187 | Int_t MParameterCalc::ReadEnv(const TEnv &env, TString prefix, Bool_t print)
|
---|
188 | {
|
---|
189 | MDataPhrase *f = new MDataPhrase;
|
---|
190 | f->SetName(fName);
|
---|
191 |
|
---|
192 | const Bool_t rc = f->ReadEnv(env, prefix, print);
|
---|
193 | if (rc!=kTRUE)
|
---|
194 | {
|
---|
195 | delete f;
|
---|
196 | return rc;
|
---|
197 | }
|
---|
198 |
|
---|
199 | delete fData;
|
---|
200 | fData = f;
|
---|
201 |
|
---|
202 | if (!fData->IsValid())
|
---|
203 | {
|
---|
204 | *fLog << err << "MParameterCalc::ReadEnv - ERROR: Inavlid rule from resource file." << endl;
|
---|
205 | return kERROR;
|
---|
206 | }
|
---|
207 |
|
---|
208 | return kTRUE;
|
---|
209 | }
|
---|