KaliVeda
Toolkit for HIC analysis
Loading...
Searching...
No Matches
KVINDRADB_e475s.cpp
1/*
2$Id: KVINDRADB_e475s.cpp,v 1.6 2009/04/09 09:23:20 ebonnet Exp $
3$Revision: 1.6 $
4$Date: 2009/04/09 09:23:20 $
5*/
6
7//Created by KVClassFactory on Tue Sep 18 11:17:53 2007
8//Author: Eric Bonnet
9
10#include "KVINDRADB_e475s.h"
11#include "Riostream.h"
12#include "TF1.h"
13#include "TString.h"
14#include "TObjString.h"
15#include "TList.h"
16
17#include "KVDBParameterSet.h"
18#include "KVINDRA.h"
19#include "KVMultiDetArray.h"
20#include "KVDataSet.h"
21#include "KVINDRADB.h"
22#include "KVNumberList.h"
23
24using namespace std;
25
26ClassImp(KVINDRADB_e475s)
27
28
29
30//------------------------------
31void KVINDRADB_e475s::init()
32//------------------------------
33
36
37{
38 //Define duplicate run, create the tables of the database.
39 cout << "init de KVINDRADB_e475s" << endl;
40
41 fCalibrations = AddTable("Calibrations", "List of available calibrations");
42
43}
44
45
46
47//------------------------------
48KVINDRADB_e475s::KVINDRADB_e475s()
49//------------------------------
50
51
53
54{
55 //Default constructor
56 init();
57}
58
59
60//------------------------------
61KVINDRADB_e475s::~KVINDRADB_e475s()
62//------------------------------
63
64
66
67{
68 //Destructor
69}
70
71
72//------------------------------
73KVINDRADB_e475s::KVINDRADB_e475s(const Char_t* name): KVINDRADB(name)
74//------------------------------
75
78
79{
80 //ctor.
81 init();
82}
83
84
85//------------------------------
86
89
90void KVINDRADB_e475s::Build()
91{
92//------------------------------
93 KVINDRADB::Build();
94
95 ReadCalibrations();
96
97}
98
99
100//------------------------------
101void KVINDRADB_e475s::ReadPedestalList()
102//------------------------------
103
104
105
106{
107
108 ifstream finput;
109 if (!OpenCalibFile("Pedestals", finput)) {
110 Error("ReadPedestalList()", "Could not open file %s",
111 GetCalibFileName("Pedestals"));
112 return;
113 }
114
115 Info("ReadPedestalList()",
116 "Reading calibration parameters...");
117
118 TString sline;
119 TString stit;
120 TObjArray* toks = NULL;
121 while (finput.good()) {
122 TString calib_det = "", calib_acq = "", calib_file = "";
123 sline.ReadLine(finput);
124 if (!sline.IsNull()) {
125 cout << sline << endl;
126 toks = sline.Tokenize(" ");
127 if (toks->GetEntries() >= 1) {
128 calib_file = ((TObjString*)(*toks)[0])->GetString();
129 if (toks->GetEntries() >= 2) {
130 calib_det = ((TObjString*)(*toks)[1])->GetString();
131 if (toks->GetEntries() >= 3) {
132 calib_acq = ((TObjString*)(*toks)[2])->GetString();
133 }
134 }
135 }
136
137 stit.Form("%s/%s", gDataSet->GetDataSetDir(), calib_file.Data());
138 ifstream fin(stit.Data());
139 ReadPedestalFile(fin, calib_det, calib_acq);
140 fin.close();
141 }
142 }
143 finput.close();
144
145
146}
147
148
149//------------------------------
150void KVINDRADB_e475s::ReadCalibrations()
151//------------------------------
152
153
154
155{
156
157 ifstream finput;
158 if (!OpenCalibFile("CalibFile", finput)) {
159 Error("ReadCalibrations()", "Could not open file %s",
160 GetCalibFileName("CalibFile"));
161 return;
162 }
163
164 Info("ReadCalibrations()",
165 "Reading calibration parameters...");
166
167 TString sline;
168 TString stit;
169 TString calib_det, calib_gain, calib_file;
170 TObjArray* toks = NULL;
171 while (finput.good()) {
172 sline.ReadLine(finput);
173 if (!sline.IsNull()) {
174 cout << sline << endl;
175 toks = sline.Tokenize(" ");
176 calib_det = ((TObjString*)(*toks)[0])->GetString();
177 calib_gain = ((TObjString*)(*toks)[1])->GetString();
178 calib_file = ((TObjString*)(*toks)[2])->GetString();
179 stit.Form("%s/%s", gDataSet->GetDataSetDir(), calib_file.Data());
180 ifstream fin(stit.Data());
181 ReadCalibFile(fin, calib_det, calib_gain);
182 fin.close();
183 }
184 }
185 finput.close();
186
187}
188
189
190//------------------------------
191
194
195void KVINDRADB_e475s::ReadCalibFile(ifstream& fin, TString dettype, TString detgain)
196{
197//------------------------------
198 Info("ReadCalibFile()",
199 "Reading calibration parameters of file %s%s", dettype.Data(), detgain.Data());
200
201 KVDBParameterSet* parset = nullptr;
202 TList list_record;
203
204 TString detname;
205 TString sline;
206 TString calib_type, calib_signal, calib_runs, calib_ring, calib_mods;
207 unique_ptr<TObjArray> toks;
208 unique_ptr<TObjArray> tokspara;
209 while (fin.good()) { //reading the file
210 sline.ReadLine(fin);
211 while (sline.BeginsWith("+KVCalibrator")) {
212 sline.ReadLine(fin);
213 while (sline.BeginsWith("<LIST> Type=")) {
214 sline.ReplaceAll("<LIST> Type=", "");
215 if (sline.Contains("(")) {
216 toks.reset(sline.Tokenize("("));
217 sline = ((TObjString*)toks->At(0))->GetString();
218 }
219 calib_type = sline;
220 sline.ReadLine(fin);
221 while (sline.BeginsWith("<LIST> Signal=")) {
222 sline.ReplaceAll("<LIST> Signal=", "");
223 calib_signal = sline;
224 cout << calib_signal << endl;
225 if (calib_signal != dettype + detgain) Error("ReadCalibFile()", "Erreur de synchro ...");
226 sline.ReadLine(fin);
227 while (sline.BeginsWith("++<LIST> Runs=")) {
228 sline.ReplaceAll("++<LIST> Runs=", "");
229 calib_runs = sline;
230 cout << calib_runs << endl;
231 sline.ReadLine(fin);
232 list_record.RemoveAll();
233 while (sline.BeginsWith("<LIST> Rings=")) {
234 sline.ReplaceAll("<LIST> Rings=", "");
235 calib_ring = sline;
236 cout << calib_ring << endl;
237 sline.ReadLine(fin);
238 while (sline.BeginsWith("<LIST> Modules=") || sline.BeginsWith("<PARAMETER> Function=")) {
239 if (sline.BeginsWith("<LIST> Modules=")) {
240 sline.ReplaceAll("<LIST> Modules=", "");
241 calib_mods = sline;
242 cout << calib_mods << endl;
243 //sprintf(detname,"%s_%02d%02d_%s",dettype.Data(),calib_ring.Atoi(),calib_mods.Atoi(),detgain.Data());
244 detname.Form("%s_%02d%02d", dettype.Data(), calib_ring.Atoi(), calib_mods.Atoi());
245 }
246 else {
247 sline.ReplaceAll("<PARAMETER> Function=", "");
248 cout << sline << endl;
249 toks.reset(sline.Tokenize(":"));
250 TF1 ff(Form("f%s", detname.Data()), ((TObjString*)toks->At(0))->GetString().Data());
251
252 TString sparaline = ((TObjString*)toks->At(1))->GetString();
253
254 tokspara.reset(sparaline.Tokenize(" "));
255 if (tokspara->GetEntries() != ff.GetNpar()) cout << tokspara->GetEntries() << " " << ff.GetNpar() << endl;
256 for (Int_t pp = 0; pp < ff.GetNpar(); pp += 1) {
257 ff.SetParameter(pp, ((TObjString*)tokspara->At(pp))->GetString().Atof());
258 printf("%d %g\n", pp, ff.GetParameter(pp));
259 }
260
261 parset = new KVDBParameterSet(detname, Form("%s(%s)", calib_type.Data(), detgain.Data()), ff.GetNpar() + 3);
262 parset->SetParamName(0, ff.GetExpFormula().Data());
263 for (Int_t pp = 0; pp < ff.GetNpar(); pp += 1) {
264 parset->SetParameter(pp + 1, ff.GetParameter(pp));
265 parset->SetParamName(pp + 1, ff.GetParName(pp));
266 }
267 parset->SetParameter(ff.GetNpar() + 1, 0.);
268 parset->SetParamName(ff.GetNpar() + 1, "xmin");
269 parset->SetParameter(ff.GetNpar() + 2, 4096.);
270 parset->SetParamName(ff.GetNpar() + 2, "xmax");
271
272 fCalibrations->AddRecord(parset);
273 list_record.Add(parset);
274
275 }
276 sline.ReadLine(fin);
277 }
278 }
279 KVNumberList nl(calib_runs);
280 nl.Begin();
281 while (!nl.End()) {
282 Int_t run = nl.Next();
283 for (Int_t ll = 0; ll < list_record.GetEntries(); ll += 1)
284 if (!(((KVDBRecord*)list_record.At(ll))->AddLink("Runs", GetRun(run))))
285 Error("ReadCalibFile()", "Pb de chargements ...");
286 }
287 }
288 }
289 }
290 }
291 }
292
293 fin.close();
294}
295
296
297//------------------------------
298
301
302void KVINDRADB_e475s::ReadPedestalFile(ifstream& fin, TString dettype, TString detgain)
303{
304//------------------------------
305 Info("ReadPedestalFile()",
306 "Reading pedestal value of file %s - %s", dettype.Data(), detgain.Data());
307
308 KVDBParameterSet* parset = NULL;
309 TList* list_record = new TList();
310
311 Char_t detname[50];
312 TString sline;
313 TString calib_type, calib_signal, calib_det, calib_acq;
314 TString calib_runs, calib_ring, calib_mods;
315 TObjArray* toks = NULL;
316 TObjArray* tokspara = NULL;
317 while (fin.good()) { //reading the file
318 sline.ReadLine(fin);
319 while (sline.BeginsWith("+KVCalibrator")) {
320 sline.ReadLine(fin);
321 while (sline.BeginsWith("<LIST> Type=")) {
322 sline.ReplaceAll("<LIST> Type=", "");
323 calib_type = sline;
324 cout << calib_type << endl;
325 sline.ReadLine(fin);
326 while (sline.BeginsWith("<LIST> KVDetector=")) {
327 sline.ReplaceAll("<LIST> KVDetector=", "");
328 calib_det = sline;
329 cout << calib_det << endl;
330 sline.ReadLine(fin);
331 while (sline.BeginsWith("<LIST> KVACQParam=")) {
332 sline.ReplaceAll("<LIST> KVACQParam=", "");
333 calib_acq = sline;
334 cout << calib_acq << endl;
335 sline.ReadLine(fin);
336 while (sline.BeginsWith("++<LIST> Runs=")) {
337 sline.ReplaceAll("++<LIST> Runs=", "");
338 calib_runs = sline;
339 cout << calib_runs << endl;
340 sline.ReadLine(fin);
341
342 list_record->RemoveAll();
343 while (sline.BeginsWith("<LIST> Rings=")) {
344 sline.ReplaceAll("<LIST> Rings=", "");
345 calib_ring = sline;
346 cout << calib_ring << endl;
347 sline.ReadLine(fin);
348 while (sline.BeginsWith("<LIST> Modules=") || sline.BeginsWith("<PARAMETER> Function=")) {
349 if (sline.BeginsWith("<LIST> Modules=")) {
350 sline.ReplaceAll("<LIST> Modules=", "");
351 calib_mods = sline;
352 cout << calib_mods << endl;
353 //sprintf(detname,"%s_%02d%02d_%s",dettype.Data(),calib_ring.Atoi(),calib_mods.Atoi(),detgain.Data());
354 sprintf(detname, "%s_%02d%02d", calib_det.Data(), calib_ring.Atoi(), calib_mods.Atoi());
355 }
356 else {
357 sline.ReplaceAll("<PARAMETER> Function=", "");
358 cout << sline << endl;
359 toks = sline.Tokenize(":");
360 TF1 ff(Form("f%s", detname), ((TObjString*)(*toks)[0])->GetString().Data());
361
362 TString sparaline = ((TObjString*)(*toks)[1])->GetString();
363 delete toks;
364
365 tokspara = sparaline.Tokenize(" ");
366 if (tokspara->GetEntries() != ff.GetNpar()) cout << tokspara->GetEntries() << " " << ff.GetNpar() << endl;
367 for (Int_t pp = 0; pp < ff.GetNpar(); pp += 1) {
368 ff.SetParameter(pp, ((TObjString*)(*tokspara)[pp])->GetString().Atof());
369 printf("%d %lf\n", pp, ff.GetParameter(pp));
370 }
371 delete tokspara;
372
373 parset = new KVDBParameterSet(detname, Form("%s", calib_acq.Data()), ff.GetNpar() + 1);
374 parset->SetParamName(0, ff.GetExpFormula().Data());
375 for (Int_t pp = 0; pp < ff.GetNpar(); pp += 1) {
376 parset->SetParameter(pp + 1, ff.GetParameter(pp));
377 parset->SetParamName(pp + 1, ff.GetParName(pp));
378 }
379
380 fPedestals->AddRecord(parset);
381 list_record->Add(parset);
382
383 }
384 sline.ReadLine(fin);
385 }
386 }
387 KVNumberList nl(calib_runs);
388 nl.Begin();
389 while (!nl.End()) {
390 Int_t run = nl.Next();
391 for (Int_t ll = 0; ll < list_record->GetEntries(); ll += 1)
392 if (!(((KVDBRecord*)list_record->At(ll))->AddLink("Runs", GetRun(run))))
393 Error("ReadCalibFile()", "Pb de chargements ...");
394 }
395 }
396 }
397 }
398 }
399 }
400 }
401 fin.close();
402
403 list_record->Clear();
404 delete list_record;
405
406}
407
408
Int_t
int Int_t
Char_t
char Char_t
name
char name[80]
Form
char * Form(const char *fmt,...)
KVDBParameterSet
To store calibration parameters in a database ,.
Definition KVDBParameterSet.h:34
KVDBParameterSet::SetParamName
void SetParamName(UShort_t i, const Char_t *name)
Definition KVDBParameterSet.cpp:116
KVDBParameterSet::SetParameter
void SetParameter(UShort_t i, Double_t val)
Definition KVDBParameterSet.h:89
KVDBRecord
Record folder for the database.
Definition KVDBRecord.h:43
KVDBTable::AddRecord
virtual Bool_t AddRecord(KVDBRecord *add)
Definition KVDBTable.cpp:74
KVDataSet::GetDataSetDir
const Char_t * GetDataSetDir() const
Definition KVDataSet.cpp:729
KVExpDB::GetCalibFileName
TString GetCalibFileName(const Char_t *type) const
Definition KVExpDB.h:109
KVExpDB::OpenCalibFile
Bool_t OpenCalibFile(const Char_t *type, std::ifstream &fs) const
Definition KVExpDB.cpp:433
KVINDRADB_e475s
Database for E475s experiment (2006)
Definition KVINDRADB_e475s.h:26
KVINDRADB_e475s::init
void init()
Define duplicate run, create the tables of the database.
Definition KVINDRADB_e475s.cpp:31
KVINDRADB_e475s::ReadPedestalFile
void ReadPedestalFile(std::ifstream &fin, TString, TString)
Definition KVINDRADB_e475s.cpp:302
KVINDRADB_e475s::~KVINDRADB_e475s
virtual ~KVINDRADB_e475s()
Destructor.
Definition KVINDRADB_e475s.cpp:61
KVINDRADB_e475s::fCalibrations
KVDBTable * fCalibrations
Definition KVINDRADB_e475s.h:31
KVINDRADB_e475s::ReadPedestalList
void ReadPedestalList()
Definition KVINDRADB_e475s.cpp:101
KVINDRADB_e475s::Build
virtual void Build()
Definition KVINDRADB_e475s.cpp:90
KVINDRADB_e475s::ReadCalibrations
void ReadCalibrations()
Definition KVINDRADB_e475s.cpp:150
KVINDRADB_e475s::KVINDRADB_e475s
KVINDRADB_e475s()
Default constructor.
Definition KVINDRADB_e475s.cpp:48
KVINDRADB_e475s::ReadCalibFile
void ReadCalibFile(std::ifstream &fin, TString, TString)
Definition KVINDRADB_e475s.cpp:195
KVINDRADB
DataBase of parameters for an INDRA campaign.
Definition KVINDRADB.h:59
KVINDRADB::GetRun
KVINDRADBRun * GetRun(Int_t run) const
Definition KVINDRADB.h:133
KVINDRADB::fPedestals
KVDBTable * fPedestals
table of pedestal files
Definition KVINDRADB.h:73
KVINDRADB::Build
virtual void Build()
Definition KVINDRADB.cpp:893
KVNumberList
Strings used to represent a set of ranges of values.
Definition KVNumberList.h:85
KVNumberList::End
Bool_t End(void) const
Definition KVNumberList.h:199
KVNumberList::Begin
void Begin(void) const
Definition KVNumberList.cpp:873
KVNumberList::Next
Int_t Next(void) const
Definition KVNumberList.cpp:841
TCollection::RemoveAll
void RemoveAll()
TCollection::GetEntries
virtual Int_t GetEntries() const
TF1::GetParameter
virtual Double_t GetParameter(const TString &name) const
TF1::GetNpar
virtual Int_t GetNpar() const
TF1::GetExpFormula
virtual TString GetExpFormula(Option_t *option="") const
TF1::GetParName
virtual const char * GetParName(Int_t ipar) const
TF1::SetParameter
virtual void SetParameter(const TString &name, Double_t value)
TList::Clear
void Clear(Option_t *option="") override
TList::Add
void Add(TObject *obj) override
TList::At
TObject * At(Int_t idx) const override
TObjArray::GetEntries
Int_t GetEntries() const override
TObject::Error
virtual void Error(const char *method, const char *msgfmt,...) const
TObject::Info
virtual void Info(const char *method, const char *msgfmt,...) const
TString::Atoi
Int_t Atoi() const
TString::Data
const char * Data() const
TString::Tokenize
TObjArray * Tokenize(const TString &delim) const
TString::BeginsWith
Bool_t BeginsWith(const char *s, ECaseCompare cmp=kExact) const
TString::IsNull
Bool_t IsNull() const
TString::Form
void Form(const char *fmt,...)
TString::Contains
Bool_t Contains(const char *pat, ECaseCompare cmp=kExact) const
TString::ReplaceAll
TString & ReplaceAll(const char *s1, const char *s2)
TString::ReadLine
std::istream & ReadLine(std::istream &str, Bool_t skipWhite=kTRUE)
TString
ClassImp
ClassImp(TPyArg)