]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TPC/AliTPCGoofieValues.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
bugfix
[u/mrichter/AliRoot.git] / TPC / AliTPCGoofieValues.cxx
CommitLineData
d4450cf4 1
2/**************************************************************************
3 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * *
5 * Author: The ALICE Off-line Project. *
6 * Contributors are mentioned in the code where appropriate. *
7 * *
8 * Permission to use, copy, modify and distribute this software and its *
9 * documentation strictly for non-commercial purposes is hereby granted *
10 * without fee, provided that the above copyright notice appears in all *
11 * copies and that both the copyright notice and this permission notice *
12 * appear in the supporting documentation. The authors make no claims *
13 * about the suitability of this software for any purpose. It is *
14 * provided "as is" without express or implied warranty. *
15 **************************************************************************/
16
17/****************************************************************************
18 * TPC Calibration Class for GOOFIE values. Drift velocity, gas composition *
19 * and the gain. *
20 ****************************************************************************/
21
22#include "AliTPCGoofieValues.h"
23#include <iostream>
24
25/*****************************************************************************
26* The class AliTPCGoofieValues allows the access to GoofieValues. The only *
27* construtor is load a data from ASCI file. The methods make Tgraphs and *
28* TSplines of the time dependace of th values. One method allows save the *
29* the graphs and spline togather with tree of allvalues into file. *
30
31Current example usage:
32
33AliTPCGoofieValues *goofieVal = new AliTPCGoofieValues("Goofie_data_january_run_01_08.txt");
34TFile f("goofieValues.root","recreate");
35goofieVal->Write("goofie");
36TBrowser b;
37And now you can browse
38
39
40*****************************************************************************/
41
42ClassImp(AliTPCGoofieValues)
43
44//____________________________________________________________________________
45AliTPCGoofieValues::AliTPCGoofieValues():
46 TNamed(),
47 fLinesInFile(0),///< lines in ASCII file
48 fStartTime(0),///< StartTime
49 fEndTime(0),///< EndTime
50 fTimeOfRun(0),///<TimeOfRun
51 fTempGrad(0),///<TempGrad
52 fAverageTempGraph(0), ///<graph of average temperatures
53 fTempGradGraph(0),///<graph of temperature gradients
54 fPressGraph(0),///<graph of pressures
55 fVdriftGraph(0),///<graph of drift velocities
56 fVdriftcorGraph(0),///<graph of drift velocities corrected
57 fGainFGraph(0),///<graph of near gain
58 fGainNGraph(0),///<graph of far gain
59 fCO2Graph(0),///<graph of CO2 content
60 fN2Graph(0),///<graph of N2 content
61
62 fAverageTempSpline(0),///< spline of average temperatures
63 fTempGradSpline(0),///<spline of temperature gradients
64 fPressSpline(0),///<spline of pressures
65 fVdriftSpline(0),///<spline of drift velocities
66 fVdriftcorSpline(0),///<spline of drift velocities corrected
67 fGainFSpline(0),///<spline of near gain
68 fGainNSpline(0),///<spline of far gain
69 fCO2Spline(0),///<spline of CO2 content
70 fN2Spline(0),
71 fGoofieValues(0)///<spline of N2 content
72{
73 //
74 // Default constructor
75 //
76}
77
78 //____________________________________________________________________________
79AliTPCGoofieValues::AliTPCGoofieValues(const char *fname):
80 TNamed(),
81 fLinesInFile(0),///< lines in ASCII file
82 fStartTime(0),///< StartTime
83 fEndTime(0),///< EndTime
84 fTimeOfRun(0),///<TimeOfRun
85 fTempGrad(0),///<TempGrad
86 fAverageTempGraph(0), ///<graph of average temperatures
87 fTempGradGraph(0),///<graph of temperature gradients
88 fPressGraph(0),///<graph of pressures
89 fVdriftGraph(0),///<graph of drift velocities
90 fVdriftcorGraph(0),///<graph of drift velocities corrected
91 fGainFGraph(0),///<graph of near gain
92 fGainNGraph(0),///<graph of far gain
93 fCO2Graph(0),///<graph of CO2 content
94 fN2Graph(0),///<graph of N2 content
95
96 fAverageTempSpline(0),///< spline of average temperatures
97 fTempGradSpline(0),///<spline of temperature gradients
98 fPressSpline(0),///<spline of pressures
99 fVdriftSpline(0),///<spline of drift velocities
100 fVdriftcorSpline(0),///<spline of drift velocities corrected
101 fGainFSpline(0),///<spline of near gain
102 fGainNSpline(0),///<spline of far gain
103 fCO2Spline(0),///<spline of CO2 content
104 fN2Spline(0),
105 fGoofieValues(0)///<spline of N2 content
106{
107 /**
108 Constructor take a values from ASCI file with raw values. <br>
109 example: <b>fname = AliTPCGoofie_run_001.txt </b><br>
110 Read the next values, in this order (each line of the file is a data point) <br>
111 fGoofieTime fGoofieTempF fGoofieTempN fGoofiePress fGoofieVdrift fGoofieVdriftcor fGoofieAreaF fGoofieAreaN fGoofieCO2 fGoofieN2
112
113 */
114
115 fLinesInFile= 0 ;
116 fStartTime= 0 ;
117 fEndTime= 0 ;
118 fTimeOfRun= 0 ;
119 fTempGrad= 0 ;
120
121 fGoofieValues = new TTree("tree","goofie values");
122 fLinesInFile = fGoofieValues->ReadFile(fname,"fGoofieTime/D:fGoofieTempF/D:fGoofieTempN/D:fGoofiePress/D:fGoofieVdrift/D:fGoofieVdriftcor/D:fGoofieAreaF/D:fGoofieAreaN/D:fGoofieCO2/D:fGoofieN2/D");
123
124 fAverageTempGraph = new TGraph();
125 fTempGradGraph = new TGraph();
126 fPressGraph = new TGraph();
127 fVdriftGraph = new TGraph();
128 fVdriftcorGraph = new TGraph();
129 fGainFGraph = new TGraph();
130 fGainNGraph = new TGraph();
131 fCO2Graph = new TGraph();
132 fN2Graph = new TGraph();
133
134 FillAllGraphs();
135 //Splines are allocated via FillAllSplines
136 FillAllSplines();
137 }
138 //________________________________________________________________________________
139 AliTPCGoofieValues::~AliTPCGoofieValues()
140 {
141 /**
142 AliTPCGoofieValues destructor
143 */
144
145 }
146 //_________________________________________________________________________________
147 Long64_t AliTPCGoofieValues::GetLinesInFile(){
148 return fLinesInFile;
149 }
150 //_________________________________________________________________________________
151 Double_t AliTPCGoofieValues::GetStartTime()
152 {
153 /**
154 take a time in beging of run from Goofie in seconds
155 */
156
157 return fStartTime;
158 }
159 //__________________________________________________________________________________
160 Double_t AliTPCGoofieValues::GetEndTime()
161 {
162 /**
163 take the time in the end of run from Goofie in seconds
164 */
165
166 return fEndTime;
167 }
168
169 //__________________________________________________________________________________
170 Double_t AliTPCGoofieValues::GetTimeOfRun()
171 {
172 /**
173 return time of run in seconds
174 */
175
176 Double_t time = GetEndTime(/*const char *fname*/) - GetStartTime(/*const char *fname*/);
177 return time;
178 }
179 //__________________________________________________________________________________
180 Double_t AliTPCGoofieValues::GetTempGrad(Double_t timeSec)
181 {
182
183 /**
184 gradient of temperature in the chosen time in run
185 */
186
187 Double_t tempGrad;
188 FillTempGradGraph();
189 tempGrad = EvalTempGrad(timeSec);
190 return tempGrad;
191 }
192 //__________________________________________________________________________________
193 void AliTPCGoofieValues::FillAllGraphs(){
194 /**
195 Fill ALL the Graphs.<br>
196 There are individual methods, to do the same.<br>
197 This function is anyway called in the ctor.
198 Selection has to be implemented:
199 //
200 Example selection :
201
202 AliTPCGoofieValues *a = new AliTPCGoofieValues("Goofie_data_january_run_01_08.txt");
203 TEventList list("listGood","listGood");
204 a->GetTree()->Draw(">>listGood","abs(fGoofieAreaN-2600)<500&&abs(fGoofieAreaF-2900)<300&&abs(fGoofieAreaF/fGoofieAreaN-1.1)<0.2","");
205 a->GetTree()->Draw("fGoofieVdrift:fGoofieTime");
206 a->GetTree()->SetEventList(&list);
207 a->GetTree()->SetMarkerColor(2);
208 //
209 a->GetTree()->Draw("fGoofieVdrift:fGoofieTime","","same*");
210 TGraph gr(200,a.GetTree()->GetV2(),a.GetTree()->GetV1());
211 //example residuals
212 AliSplineFit fit;
213 fit.SetGraph(&gr)
214 fit->SetMinPoints(200);
215 fit->InitKnots(&gr,15,0,0.2)
216 fit.SplineFit(1)
217 fit.MakeDiffHisto(&gr)->Draw();
218
219
220 */
221 Long64_t nevent = fGoofieValues->GetEntries();
222
223 //temporal for reading the branches
224 Double_t fGoofieTime = 0;
225 Double_t fGoofieTempN = 0;
226 Double_t fGoofieTempF = 0;
227 Double_t fGoofiePress = 0;
228 Double_t fGoofieAreaF = 0;
229 Double_t fGoofieAreaN = 0;
230 Double_t fGoofieVdriftcor = 0;
231 Double_t fGoofieVdrift = 0;
232 Double_t fGoofieCO2 = 0;
233 Double_t fGoofieN2 = 0;
234 fGoofieValues->SetBranchAddress("fGoofieTime",&fGoofieTime);
235 fGoofieValues->SetBranchAddress("fGoofieTempN",&fGoofieTempN);
236 fGoofieValues->SetBranchAddress("fGoofieTempF",&fGoofieTempF);
237 fGoofieValues->SetBranchAddress("fGoofiePress",&fGoofiePress);
238 fGoofieValues->SetBranchAddress("fGoofieAreaF",&fGoofieAreaF);
239 fGoofieValues->SetBranchAddress("fGoofieAreaN",&fGoofieAreaN);
240 fGoofieValues->SetBranchAddress("fGoofieVdriftcor",&fGoofieVdriftcor);
241 fGoofieValues->SetBranchAddress("fGoofieVdrift",&fGoofieVdrift);
242 fGoofieValues->SetBranchAddress("fGoofieCO2",&fGoofieCO2);
243 fGoofieValues->SetBranchAddress("fGoofieN2",&fGoofieN2);
244 fStartTime = 0; //fGoofieValues->GetBranch("fGoofieTime")->GetEntry(0);
245
246 for (int i = 0; i < nevent; ++i){
247 if (fGoofieValues->GetEvent(i)<0){
248 // cout<< "you're done, man !" << endl;
249 continue;
250 }
251 if (i==0) fStartTime = fGoofieTime;
252 else if (i == (nevent-1)) fEndTime = fGoofieTime;
253 //// cout<< " reading " << i << endl;
254 fAverageTempGraph->SetPoint(i,(fGoofieTime - fStartTime),(fGoofieTempF+fGoofieTempN)/2);
255 fTempGradGraph->SetPoint(i,(fGoofieTime - fStartTime),(fGoofieTempF-fGoofieTempN)/25);
256 fPressGraph->SetPoint(i,(fGoofieTime - fStartTime),fGoofiePress);
257 fGainFGraph->SetPoint(i,(fGoofieTime - fStartTime),fGoofieAreaF);
258 fGainNGraph->SetPoint(i,(fGoofieTime - fStartTime),fGoofieAreaN);
259 fVdriftGraph->SetPoint(i,(fGoofieTime - fStartTime),fGoofieVdrift);
260 fVdriftcorGraph->SetPoint(i,(fGoofieTime - fStartTime),fGoofieVdriftcor);
261 fCO2Graph->SetPoint(i,(fGoofieTime - fStartTime),fGoofieCO2);
262 fN2Graph->SetPoint(i,(fGoofieTime - fStartTime),fGoofieN2);
263 }
264
265 }
266
267 void AliTPCGoofieValues::FillAllSplines(){
268 /**
269 Fill ALL the Splines.<br>
270 There are individual methods, to do the same.<br>
271 This function is anyway called in the ctor.
272 */
273
274 fAverageTempSpline = new TSpline3("temperature",fAverageTempGraph);
275 fTempGradSpline = new TSpline3("temperature gradient",fTempGradGraph);
276 fPressSpline = new TSpline3("pressure",fPressGraph);
277 fVdriftcorSpline = new TSpline3("vdriftcor",fVdriftcorGraph);
278 fVdriftSpline = new TSpline3("vdrift",fVdriftGraph);
279 fGainFSpline = new TSpline3("gainF",fGainFGraph);
280 fGainNSpline = new TSpline3("gainN",fGainNGraph);
281 fCO2Spline = new TSpline3("co2",fCO2Graph);
282 fN2Spline = new TSpline3("n2",fN2Graph);
283
284 }
285
286 //___________________________________________________________________________________
287 void AliTPCGoofieValues::FillAverageTempGraph()
288 {
289 /**
290 graph of temperature vs time. <br>
291 If the graph is already filled, it's only putting a title.
292 */
293
294 fAverageTempGraph->SetTitle("time dependence of temperature");
295
296 if (!fAverageTempGraph->GetN()==0){
297 // cout<< " graph already filled !!! " << endl;
298 }
299 else {
300 Double_t fGoofieTime = 0;
301 Double_t fGoofieTempN = 0;
302 Double_t fGoofieTempF = 0;
303 fGoofieValues->SetBranchAddress("fGoofieTime",&fGoofieTime);
304 fGoofieValues->SetBranchAddress("fGoofieTempN",&fGoofieTempN);
305 fGoofieValues->SetBranchAddress("fGoofieTempF",&fGoofieTempF);
306
307 Long64_t nevent = fGoofieValues->GetEntries();
308
309 for (int i = 0; i < nevent; ++i){
310 if (fGoofieValues->GetEvent(i)<0){
311 // cout<< "you're done, man !" << endl;
312 continue;
313 }
314 if (i==0) fStartTime = fGoofieTime;
315 else if (i == (nevent-1)) fEndTime = fGoofieTime;
316 fAverageTempGraph->SetPoint(i,(fGoofieTime - fStartTime),(fGoofieTempF+fGoofieTempN)/2);
317 }
318 }
319 }
320 //__________________________________________________________________________________
321 void AliTPCGoofieValues::FillTempGradGraph()
322 {
323 /**
324 graph of temperature gradient [K/cm]<br>
325 If the graph is already filled, it's only putting a title.
326 */
327
328 fTempGradGraph->SetTitle("time dpendance of Temperature`s gradient");
329
330 if(! fTempGradGraph->GetN()==0){
331 // cout<< " graph already filled !!! " << endl;
332 }
333 else{
334 Double_t fGoofieTime = 0;
335 Double_t fGoofieTempN = 0;
336 Double_t fGoofieTempF = 0;
337 fGoofieValues->SetBranchAddress("fGoofieTime",&fGoofieTime);
338 fGoofieValues->SetBranchAddress("fGoofieTempN",&fGoofieTempN);
339 fGoofieValues->SetBranchAddress("fGoofieTempF",&fGoofieTempF);
340
341 Long64_t nevent = fGoofieValues->GetEntries();
342 for (int i = 0; i < nevent; ++i){
343 if (fGoofieValues->GetEvent(i)<0){
344 // cout<< "you're done, man !" << endl;
345 continue;
346 }
347 if (i==0) fStartTime = fGoofieTime;
348 else if (i == (nevent-1)) fEndTime = fGoofieTime;
349 fTempGradGraph->SetPoint(i,(fGoofieTime - fStartTime),(fGoofieTempF-fGoofieTempN)/25);
350 }
351 }
352 }
353
354 //____________________________________________________________________________________
355 void AliTPCGoofieValues::FillPressGraph()
356 {
357 /**
358 Graph of pressure<br>
359 If the graph is already filled, it's only putting a title.
360 */
361
362 fPressGraph->SetTitle("time dpendance of Pressure");
363 if(!fPressGraph->GetN()==0){
364 // cout<< " graph already filled !!! " << endl;
365 }
366 else{
367 Double_t fGoofieTime = 0;
368 Double_t fGoofiePress = 0;
369 fGoofieValues->SetBranchAddress("fGoofieTime",&fGoofieTime);
370 fGoofieValues->SetBranchAddress("fGoofiePress",&fGoofiePress);
371
372 Long64_t nevent = fGoofieValues->GetEntries();
373 for (int i = 0; i < nevent; ++i){
374 if (fGoofieValues->GetEvent(i)<0){
375 // cout<< "you're done, man !" << endl;
376 continue;
377 }
378 if (i==0) fStartTime = fGoofieTime;
379 else if (i == (nevent-1)) fEndTime = fGoofieTime;
380 fPressGraph->SetPoint(i,(fGoofieTime - fStartTime),fGoofiePress);
381 }
382 }
383 }
384
385
386//________________________________________________________________________________________________
387 void AliTPCGoofieValues::FillGainFGraph()
388 {
389
390 /**
391 return graph of Gain in the far point of Goofie`s detector<br>
392 If the graph is already filled, it's only putting a title.
393 */
394
395 fGainFGraph->SetTitle("time dpendance of Gain in the far point in natural unit");
396
397 if(!fGainFGraph->GetN()==0){
398 // cout<< " graph already filled !!! " << endl;
399 }
400 else{
401 Double_t fGoofieTime = 0;
402 Double_t fGoofieAreaF = 0;
403 fGoofieValues->SetBranchAddress("fGoofieTime",&fGoofieTime);
404 fGoofieValues->SetBranchAddress("fGoofieAreaF",&fGoofieAreaF);
405 Long64_t nevent = fGoofieValues->GetEntries();
406 for (int i = 0; i < nevent; ++i){
407 if (fGoofieValues->GetEvent(i)<0){
408 // cout<< "you're done, man !" << endl;
409 continue;
410 }
411 if (i==0) fStartTime = fGoofieTime;
412 else if (i == (nevent-1)) fEndTime = fGoofieTime;
413 fGainFGraph->SetPoint(i,(fGoofieTime - fStartTime),fGoofieAreaF);
414 }
415 }
416 }
417
418//_____________________________________________________________________________________
419 void AliTPCGoofieValues::FillGainNGraph()
420 {
421 /**
422 Graph of Gain in the near point of Goofie`s detector<br>
423 If the graph is already filled, it's only putting a title.
424 */
425
426 fGainNGraph->SetTitle("time dependance of Gain in the far poin in natural unit");
427 if(!fGainNGraph->GetN()==0){
428 // cout<< " graph already filled !!! " << endl;
429 }
430 else{
431
432 Double_t fGoofieTime = 0;
433 Double_t fGoofieAreaN = 0;
434 fGoofieValues->SetBranchAddress("fGoofieTime",&fGoofieTime);
435 fGoofieValues->SetBranchAddress("fGoofieAreaN",&fGoofieAreaN);
436 Long64_t nevent = fGoofieValues->GetEntries();
437 for (int i = 0; i < nevent; ++i){
438 if (fGoofieValues->GetEvent(i)<0){
439 // cout<< "you're done, man !" << endl;
440 continue;
441 }
442 if (i==0) fStartTime = fGoofieTime;
443 else if (i == (nevent-1)) fEndTime = fGoofieTime;
444 fGainNGraph->SetPoint(i,(fGoofieTime - fStartTime),fGoofieAreaN);
445 }
446 }
447
448 }
449//_______________________________________________________________________________________
450void AliTPCGoofieValues::FillVdriftGraph()
451 {
452 /**
453 Graph of raw Vdrift<br>
454 If the graph is already filled, it's only putting a title.
455 */
456
457 fVdriftGraph->SetTitle("time dpendance of raw Vdrift");
458 if(!fVdriftGraph->GetN()==0){
459 // cout<< " graph already filled !!! " << endl;
460 }
461 else{
462 Double_t fGoofieTime = 0;
463 Double_t fGoofieVdrift = 0;
464 fGoofieValues->SetBranchAddress("fGoofieTime",&fGoofieTime);
465 fGoofieValues->SetBranchAddress("fGoofieVdrift",&fGoofieVdrift);
466 Long64_t nevent = fGoofieValues->GetEntries();
467 for (int i = 0; i < nevent; ++i){
468 if (fGoofieValues->GetEvent(i)<0){
469 // cout<< "you're done, man !" << endl;
470 continue;
471 }
472 if (i==0) fStartTime = fGoofieTime;
473 else if (i == (nevent-1)) fEndTime = fGoofieTime;
474 fVdriftGraph->SetPoint(i,(fGoofieTime - fStartTime),fGoofieVdrift);
475 }
476 }
477 }
478//_____________________________________________________________________________________________
479void AliTPCGoofieValues::FillVdriftcorGraph()
480 {
481 /**
482 Graph of Vdrift corrected by temperature and pressure<br>
483 If the graph is already filled, it's only putting a title.
484 */
485
486 fVdriftcorGraph->SetTitle("time dpendance of raw Vdriftcor");
487 if(!fVdriftcorGraph->GetN()==0){
488 // cout<< " graph already filled !!! " << endl;
489 }
490 else{
491 Double_t fGoofieTime = 0;
492 Double_t fGoofieVdriftcor = 0;
493 fGoofieValues->SetBranchAddress("fGoofieTime",&fGoofieTime);
494 fGoofieValues->SetBranchAddress("fGoofieVdriftcor",&fGoofieVdriftcor);
495 Long64_t nevent = fGoofieValues->GetEntries();
496 for (int i = 0; i < nevent; ++i){
497 if (fGoofieValues->GetEvent(i)<0){
498 // cout<< "you're done, man !" << endl;
499 continue;
500 }
501 if (i==0) fStartTime = fGoofieTime;
502 else if (i == (nevent-1)) fEndTime = fGoofieTime;
503 fVdriftcorGraph->SetPoint(i,(fGoofieTime - fStartTime),fGoofieVdriftcor);
504 }
505 }
506 }
507//____________________________________________________________________________________________
508void AliTPCGoofieValues::FillCO2Graph()
509 {
510 /**
511 Graph of CO2<br>
512 If the graph is already filled, it's only putting a title.
513 */
514
515 fCO2Graph->SetTitle("time dependance of CO2");
516 if(!fCO2Graph->GetN()==0){
517 // cout<< " graph already filled !!! " << endl;
518 }
519 else{
520 Double_t fGoofieTime = 0;
521 fGoofieValues->SetBranchAddress("fGoofieTime",&fGoofieTime);
522 Double_t fGoofieCO2 = 0;
523 fGoofieValues->SetBranchAddress("fGoofieCO2",&fGoofieCO2);
524 Long64_t nevent = fGoofieValues->GetEntries();
525 for (int i = 0; i < nevent; ++i){
526 if (fGoofieValues->GetEvent(i)<0){
527 // cout<< "you're done, man !" << endl;
528 continue;
529 }
530 if (i==0) fStartTime = fGoofieTime;
531 else if (i == (nevent-1)) fEndTime = fGoofieTime;
532 fCO2Graph->SetPoint(i,(fGoofieTime - fStartTime),fGoofieCO2);
533 }
534 }
535 }
536//_____________________________________________________________________________________________
537void AliTPCGoofieValues::FillN2Graph()
538 {
539 /**
540 Graph of N2<br>
541 If the graph is already filled, it's only putting a title.
542 */
543
544
545 fN2Graph->SetTitle("time dependance of N2");
546 if(!fN2Graph->GetN()==0){
547 // cout<< " graph already filled !!! " << endl;
548 }
549 else{
550 Double_t fGoofieTime = 0;
551 fGoofieValues->SetBranchAddress("fGoofieTime",&fGoofieTime);
552 Double_t fGoofieN2 = 0;
553 fGoofieValues->SetBranchAddress("fGoofieN2",&fGoofieN2);
554 Long64_t nevent = fGoofieValues->GetEntries();
555 for (int i = 0; i < nevent; ++i){
556 if (fGoofieValues->GetEvent(i)<0){
557 // cout<< "you're done, man !" << endl;
558 continue;
559 }
560 if (i==0) fStartTime = fGoofieTime;
561 else if (i == (nevent-1)) fEndTime = fGoofieTime;
562 fN2Graph->SetPoint(i,(fGoofieTime - fStartTime),fGoofieN2);
563 }
564 }
565 }
566
567//_________________________________________________________________________________________
568void AliTPCGoofieValues::FillAverageTempSpline()
569{
570 /**
571 TSpline object of average temperature<br>
572 If the spline is already filled, it's printing a message.
573 */
574 if(!fAverageTempSpline)
575 fAverageTempSpline = new TSpline3("temperature",fAverageTempGraph);
576 //else
577 // cout<< " spline already filled !!! " << endl;
578
579}
580
581//_________________________________________________________________________________________
582void AliTPCGoofieValues::FillTempGradSpline()
583{
584
585 /**
586 TSpline object of temperature gradient<br>
587 If the spline is already filled, it's printing a message.
588 */
589 if(!fTempGradSpline)
590 fTempGradSpline = new TSpline3("temperature gradient",fTempGradGraph);
591 //else
592 // cout<< " spline already filled !!! " << endl;
593}
594
595//_________________________________________________________________________________________
596void AliTPCGoofieValues::FillPressSpline()
597{
598
599 /**
600 TSpline object of pressure<br>
601 If the spline is already filled, it's printing a message.
602 */
603 if(!fPressSpline)
604 fPressSpline = new TSpline3("pressure",fPressGraph);
605 //else
606 // cout<< " spline already filled !!! " << endl;
607
608}
609
610//_________________________________________________________________________________________
611void AliTPCGoofieValues::FillVdriftSpline()
612{
613
614 /**
615 TSpline object of drift velocity<br>
616 If the spline is already filled, it's printing a message.
617 */
618 if(!fVdriftSpline)
619 fVdriftSpline = new TSpline3("vdrift",fVdriftGraph);
620 //else{
621 // cout<< " spline already filled !!! " << endl;
622 //}
623}
624
625//_________________________________________________________________________________________
626void AliTPCGoofieValues::FillVdriftcorSpline()
627{
628
629
630 /**
631 TSpline object of drift velocit corrected<br>
632 If the spline is already filled, it's printing a message.
633 */
634 if(!fVdriftcorSpline)
635 fVdriftcorSpline = new TSpline3("vdriftcor",fVdriftcorGraph);
636 //else
637 // cout<< " spline already filled !!! " << endl;
638
639}
640
641//_________________________________________________________________________________________
642void AliTPCGoofieValues::FillGainFSpline()
643{
644
645 /**
646 return TSpline object of gain in far point<br>
647 If the spline is already filled, it's printing a message.
648 */
649 if(!fGainFSpline)
650 fGainFSpline = new TSpline3("gainF",fGainFGraph);
651 //else
652 // cout<< " spline already filled !!! " << endl;
653
654}
655
656//_________________________________________________________________________________________
657void AliTPCGoofieValues::FillGainNSpline()
658{
659
660 /**
661 TSpline object of gain in near point<br>
662 If the spline is already filled, it's printing a message.
663 */
664 if(!fGainNSpline)
665 fGainNSpline = new TSpline3("gainN",fGainNGraph);
666 //else
667 // cout<< " spline already filled !!! " << endl;
668
669}
670
671//_________________________________________________________________________________________
672void AliTPCGoofieValues::FillCO2Spline()
673{
674
675 /**
676 TSpline object of average temperature<br>
677 If the spline is already filled, it's printing a message.
678 */
679 if(!fCO2Spline)
680 fCO2Spline = new TSpline3("co2",fCO2Graph);
681 //else
682 // cout<< " spline already filled !!! " << endl;
683}
684
685//_________________________________________________________________________________________
686void AliTPCGoofieValues::FillN2Spline()
687{
688
689 /**
690 TSpline object of average temperature<br>
691 If the spline is already filled, it's printing a message.
692 */
693 if(!fN2Spline)
694 fN2Spline = new TSpline3("n2",fN2Graph);
695 //else
696 // cout<< " spline already filled !!! " << endl;
697}
698
699//__________________________________________________________________________________________
700Double_t AliTPCGoofieValues::EvalTempGrad(Double_t timeSec)
701{
702 if(!fTempGradSpline)
703 fTempGradSpline = new TSpline3("temperature gradient",fTempGradGraph);
704 Double_t a = fTempGradSpline->Eval(timeSec);
705 return a;
706}
707
708//___________________________________________________________________________________________
709Double_t AliTPCGoofieValues::EvalAverageTemp(Double_t timeSec)
710{
711 if(!fAverageTempSpline)
712 fAverageTempSpline = new TSpline3("temperature",fAverageTempGraph);
713 Double_t a = fAverageTempSpline->Eval(timeSec);
714 return a;
715}
716
717//___________________________________________________________________________________________
718Double_t AliTPCGoofieValues::EvalPress(Double_t timeSec)
719{
720 if(!fPressSpline)
721 fPressSpline = new TSpline3("pressure",fPressGraph);
722 Double_t a = fPressSpline->Eval(timeSec);
723 return a;
724}
725
726//___________________________________________________________________________________________
727Double_t AliTPCGoofieValues::EvalVdrift(Double_t timeSec)
728{
729 if(!fVdriftSpline)
730 fVdriftSpline = new TSpline3("vdrift",fVdriftGraph);
731 Double_t a = fVdriftSpline->Eval(timeSec);
732 return a;
733}
734
735//___________________________________________________________________________________________
736Double_t AliTPCGoofieValues::EvalVdriftcor(Double_t timeSec)
737{
738 if(!fVdriftcorSpline)
739 fVdriftcorSpline = new TSpline3("vdriftcor",fVdriftcorGraph);
740 Double_t a = fVdriftcorSpline->Eval(timeSec);
741 return a;
742}
743
744//___________________________________________________________________________________________
745Double_t AliTPCGoofieValues::EvalGainF(Double_t timeSec)
746{
747 if(!fGainFSpline)
748 fGainFSpline = new TSpline3("gainF",fGainFGraph);
749 Double_t a = fGainFSpline->Eval(timeSec);
750 return a;
751}
752
753//___________________________________________________________________________________________
754Double_t AliTPCGoofieValues::EvalGainN(Double_t timeSec)
755{
756 if(!fGainNSpline)
757 fGainNSpline = new TSpline3("gainN",fGainNGraph);
758 Double_t a = fGainNSpline->Eval(timeSec);
759 return a;
760}
761
762//___________________________________________________________________________________________
763Double_t AliTPCGoofieValues::EvalCO2(Double_t timeSec)
764{
765 if(!fCO2Spline)
766 fCO2Spline = new TSpline3("co2",fCO2Graph);
767 Double_t a = fCO2Spline->Eval(timeSec);
768 return a;
769}
770
771//___________________________________________________________________________________________
772Double_t AliTPCGoofieValues::EvalN2(Double_t timeSec)
773{
774 if(!fN2Spline)
775 fN2Spline = new TSpline3("n2",fN2Graph);
776 Double_t a = fN2Spline->Eval(timeSec);
777 return a;
778}
779
780void AliTPCGoofieValues::PrintTree(){
781 /**
782 Testing function: it prints all the information stored in the tree
783 */
784
785 Long64_t nevent = fGoofieValues->GetEntries();
786 // cout<< " number of entries: " << nevent << endl;
787 //temporal for reading the branches
788 Double_t fGoofieTime = 0;
789 Double_t fGoofieTempN = 0;
790 Double_t fGoofieTempF = 0;
791 Double_t fGoofiePress = 0;
792 Double_t fGoofieAreaF = 0;
793 Double_t fGoofieAreaN = 0;
794
795 Double_t fGoofieVdriftcor = 0;
796 Double_t fGoofieVdrift = 0;
797 Double_t fGoofieCO2 = 0;
798 Double_t fGoofieN2 = 0;
799
800 fGoofieValues->SetBranchAddress("fGoofieTime",&fGoofieTime);
801 fGoofieValues->SetBranchAddress("fGoofieTempN",&fGoofieTempN);
802 fGoofieValues->SetBranchAddress("fGoofieTempF",&fGoofieTempF);
803 fGoofieValues->SetBranchAddress("fGoofiePress",&fGoofiePress);
804
805 fGoofieValues->SetBranchAddress("fGoofieAreaF",&fGoofieAreaF);
806 fGoofieValues->SetBranchAddress("fGoofieAreaN",&fGoofieAreaN);
807
808 fGoofieValues->SetBranchAddress("fGoofieVdriftcor",&fGoofieVdriftcor);
809 fGoofieValues->SetBranchAddress("fGoofieVdrift",&fGoofieVdrift);
810
811 fGoofieValues->SetBranchAddress("fGoofieCO2",&fGoofieCO2);
812 fGoofieValues->SetBranchAddress("fGoofieN2",&fGoofieN2);
813
814 for (int j = 0; j < nevent; ++j){
815 if (fGoofieValues->GetEvent(j)<0){
816 // cout<< "you're done, man !" << endl;
817 continue;
818 }
819// // cout<< " reading event: " << j << " " << fGoofieTime << " "
820// << fGoofieTempN << " " << fGoofieTempF << " "
821// << fGoofiePress << " " << fGoofiePress << " "
822// << fGoofieAreaF << " " << fGoofieAreaN << " "
823// << fGoofieVdriftcor << " " << fGoofieVdrift << " "
824// << fGoofieCO2 << " " << fGoofieN2 << endl;
825 }
826}
827
828
829
830
831
832
833
834