3 Macro to perform fits of the Laser Central electrode data
4 Several fit methods implemented
6 0. RebuildCE("ce.root","pul.root"); - rebuild data from the scratch
7 - the data will be storered in file inname
9 1. RebuildData() - transform arbitrary layout of the Input data to the internal format
10 StoreData(); - The data tree expected in file inname (see variable bellow)
11 StoreTree(); - Modify inname and xxside and tcor in order to transform data
14 2. MakeFit(); - Make a fit of the data - already in internal format
18 3. MakeRes(); - Make the final calibration + combination of different components
20 4. LoadViewer(); - Browse the fit parameters
24 gSystem->AddIncludePath("-I$ALICE_ROOT/TPC -I$ALICE_ROOT/STAT");
25 gSystem->Load("libSTAT.so");
26 .L $ALICE_ROOT/TPC/CalibMacros/AnalyzeLaserCE.C+
39 Calibration viewer variables:
41 Result - resulting correction
42 out - outlyers not used for fit
43 tcor - offset specified by user before fitting
44 timeF1 - sector local fit - plane
45 timeF2 - sector local fit - parabola
48 out - outlyers not used for fit
49 tcor - offset specified by user before fitting
50 timeF1 - sector time local fit - plane
51 timeF2 - sector time local fit - parabola
52 qF1 - sector q local fit - plane
53 qF2 - sector q local fit - parabola
57 ffit1 - adding common shifts - alpha dependendent
58 ffit2 - adding opposite shifts - alpha dependent
60 fGXY - global fit parameter - XY
61 fInOut - global fit parameter - inner-outer sector matching
62 fLX - global LX dependence
64 Gloabl fit o consist of
65 -fGXY~-fLX~-fTL~-fOff~:ffit0~
68 // Control variable - check results
71 ffit2~-(timeIn~):lx~ - fit value minus input time
74 (timeF2~-ffit2~+fTL~+fInOut~):Result~
76 timeF2~-Result~:ffit2~-fTL~-fInOut~
84 #include "TEntryList.h"
86 #include "TStatToolkit.h"
87 #include "AliTPCCalibViewer.h"
88 #include "AliTPCCalibViewerGUI.h"
89 #include "AliTPCPreprocessorOnline.h"
90 #include "AliTPCCalibCE.h"
91 #include "AliTPCCalibPulser.h"
92 #include "TStopwatch.h"
94 //Define interesting variables - file names
96 char * inname = "treeCE.root"; // input file with tree
98 // variable name definition in input tree - change it according the input
100 TString qaside("CE_Q");
101 TString taside("CE_T");
102 TString raside("CE_RMS");
103 TString qcside("CE_Q");
104 TString tcside("CE_T");
105 TString rcside("CE_RMS");
108 // correction variable - usually Pulser time
110 TString tcor("-pulCorr");
113 char * fname = "treefitCE.root"; // output file with tree
114 char * oname = "fitCE.root"; // output file with CalPads fit
120 AliTPCCalPad *calPadIn = 0; // original time pad
121 AliTPCCalPad *calPadF1 = 0; // original time pad - fit plane
122 AliTPCCalPad *calPadF2 = 0; // original time pad - fit parabola
123 AliTPCCalPad *calPadQIn = 0; // original Q pad
124 AliTPCCalPad *calPadQF1 = 0; // original Q pad
125 AliTPCCalPad *calPadQF2 = 0; // original Q pad
127 AliTPCCalPad *calPadCor = 0; // base correction CalPad
128 AliTPCCalPad *calPadOut = 0; // outlyer CalPad
132 const Float_t tThr=0.3; // max diff in the sector
133 const Float_t qThr0=0.5; // max Q diff in the sector
134 const Float_t qThr1=2; // max Q diff in the sector
139 AliTPCCalPad *calPad0 = 0; // global fit 0 - base
140 AliTPCCalPad *calPad1 = 0; // global fit 1 - common behavior rotation -A-C
141 AliTPCCalPad *calPad2 = 0; // gloabl fit 2 - CE missalign rotation A-C
143 AliTPCCalPad *calPadInOut = 0; // misaalign in-out
144 AliTPCCalPad *calPadLX = 0; // local x missalign
145 AliTPCCalPad *calPadTL = 0; // tan
146 AliTPCCalPad *calPadQ = 0; // time (q) correction
147 AliTPCCalPad *calPadGXY = 0; // global XY missalign (drift velocity grad)
148 AliTPCCalPad *calPadOff = 0; // normalization offset fit
149 AliTPCCalPad *calPadRes = 0; // final calibration
151 TObjString strFit0=""; // string fit 0
152 TObjString strFit1=""; // string fit 1
153 TObjString strFit2=""; // string fit 2
154 TVectorD vecFit0; // parameters fit 0
155 TVectorD vecFit1; // parameters fit 1
156 TVectorD vecFit2; // parameters fit 2
157 TEntryList *elist=0; // event list -slection criterias
161 AliTPCCalibViewerGUI * viewer=0; //viewerGUI
162 AliTPCCalibViewer *makePad=0; //viewer
163 TTree * tree=0; // working tree
166 void RebuildData(); // transform the input data to the fit format
167 void MakeFit(); // make fits
168 void MakeFitPulser(); // make fit for pulser correction data
170 // internal functions
172 void MakeAliases0(); // Make Aliases 0 - for user tree
173 void MakeAliases1(); // Make Aliases 1 - for default tree
174 void LoadData(); // Load data
175 void StoreData(); // store current data
176 void StoreTree(); // store fit data in the output tree
188 void MakeFitPulser(){
192 TVectorD vec0,vec1,vec2;
194 TString fitvar="P_T.fElements";
195 TCut out("abs(P_T.fElements/P_T_Mean.fElements-1.001)<.002");
196 TCut fitadd("P_T.fElements>446");
198 fstring+="(sector>=0&§or<9)++";
199 fstring+="(sector>=9&§or<18)++";
200 fstring+="(sector>=18&§or<27)++";
201 fstring+="(sector>=27&§or<36)++";
202 fstring+="(sector>=36&§or<45)++";
203 fstring+="(sector>=45&§or<54)++";
204 fstring+="(sector>=54&§or<63)++";
205 fstring+="(sector>=63&§or<72)";
207 TString *pulOff =stat.FitPlane(tree,fitvar.Data(),fstring.Data(),(out+fitadd).GetTitle(),chi2,npoints,vec0,mat);
209 tree->SetAlias("pul0",pulOff->Data());
210 tree->SetAlias("pulCorr","P_T.fElements-pul0");
211 tree->SetAlias("pulOut",out.GetTitle());
220 makePad = new AliTPCCalibViewer(fname);
221 tree = makePad->GetTree();
233 fstring+="side++"; // offset on 2 different sides //1
234 //fstring+="(1/qp)++"; // Q -threshold effect correction //2
235 //fstring+="(qp)++"; // Q -threshold effect correction //3
236 fstring+="(inn)++"; // inner outer misalign - common //4
237 fstring+="(side*inn)++"; // - opposite //5
239 fstring+="(gyr)++"; // drift velocity gradient - common //6
240 fstring+="(side*gyr)++"; // - opposite //7
241 fstring+="(gxr)++"; // global x tilt - common //8
242 fstring+="(side*gxr)++"; // - opposite //9
244 fstring+="tl^2++"; // local phi correction //10
246 fstring+="(lxr)++"; // zr angle - common //11
247 fstring+="(side*lxr)++"; // - opposite //12
248 fstring+="(inn*lxr)++"; // inner outer angle - common //13
249 fstring+="(side*inn*lxr)++";// - opposite //14
250 fstring+="(lxr^2)++"; // zr second - common //15
251 fstring+="(side*lxr^2)++"; // - opposite //16
252 fstring+="(inn*lxr^2)++"; // zr second - common //15
253 fstring+="(inn*side*lxr^2)++"; // - opposite //16
255 printf("Fit0\t start\n");
256 TString *fit0 =stat.FitPlane(tree,"dt",fstring.Data(),"cutF&&cutCE",chi2,npoints,vecFit0,mat,0.9);
257 tree->SetAlias("f0",fit0->Data());
258 strFit0 = fit0->Data();
259 printf("Fit0\t end\n");
260 printf("Chi2/npoints\t=%f\n",TMath::Sqrt(chi2/npoints));
261 fit0->Tokenize("++")->Print();
263 printf("Global tendencies extraction\n");
267 TString tmpstr = fstring;
268 TObjArray *arr = tmpstr.Tokenize("++");
269 TString fitQ("0"); // q correction
270 TString fitLX("0"); // lx correction
271 TString fitInOut("0"); // inner-outer - match
272 TString fitGXY("0"); // global xy fit
273 TString fitOff("0"); // side offsets
274 TString fitTL("0"); // side offsets
281 for(Int_t i=0;i<arr->GetEntriesFast();i++){
282 if (!arr->At(i)) continue;
283 TString *fitstr = new TString(arr->At(i)->GetName());
285 //Bool_t isQ = fitstr->Contains("qp)");
286 Bool_t isRot = fitstr->Contains("sin(")+fitstr->Contains("cos(");
287 Bool_t isLX = fitstr->Contains("lxr");
288 Bool_t isIn = fitstr->Contains("inn");
289 Bool_t isGXY = fitstr->Contains("gxr")+fitstr->Contains("gyr");
290 if (fitstr->Contains("tl^2")){
292 fitTL+=(*fitstr)+"*";
297 fitGXY+=(*fitstr)+"*";
298 fitGXY+=vecFit0[i+1];
301 if (isLX&&!isRot&&!isIn){
303 fitLX+=(*fitstr)+"*";
309 fitInOut+=(*fitstr)+"*";
310 fitInOut+=vecFit0[i+1];
315 tree->SetAlias("fInOut",fitInOut.Data());
316 tree->SetAlias("fLX",fitLX.Data());
317 tree->SetAlias("fGXY",fitGXY.Data());
318 tree->SetAlias("fOff",fitOff.Data());
319 //tree->SetAlias("fQ",fitQ.Data());
320 tree->SetAlias("fTL",fitTL.Data());
324 // Common "deformation" tendencies
326 fstring+="(sin(atan2(gy.fElements,gx.fElements)))++";
327 fstring+="(cos(atan2(gy.fElements,gx.fElements)))++";
329 fstring+="(sin(atan2(gy.fElements,gx.fElements)*2))++";
330 fstring+="(cos(atan2(gy.fElements,gx.fElements)*2))++";
331 fstring+="(sin(atan2(gy.fElements,gx.fElements)*3))++";
332 fstring+="(cos(atan2(gy.fElements,gx.fElements)*3))++";
334 fstring+="(sin(atan2(gy.fElements,gx.fElements)*2))*lxr++";
335 fstring+="(cos(atan2(gy.fElements,gx.fElements)*2))*lxr++";
336 fstring+="(sin(atan2(gy.fElements,gx.fElements)*3))*lxr++";
337 fstring+="(cos(atan2(gy.fElements,gx.fElements)*3))*lxr++";
340 TString *fit1 =stat.FitPlane(tree,"dt",fstring.Data(),"cutF&&cutCE",chi2,npoints,vecFit1,mat,0.9);
341 tree->SetAlias("f1",fit1->Data());
342 strFit1 = fit1->Data();
343 printf("Fit1\t end\n");
344 printf("Chi2/npoints\t=%f\n",TMath::Sqrt(chi2/npoints));
345 fit1->Tokenize("++")->Print();
348 // Central electrode "deformation"
350 fstring+="(side*sin(atan2(gy.fElements,gx.fElements)))++";
351 fstring+="(side*cos(atan2(gy.fElements,gx.fElements)))++";
353 fstring+="(side*sin(atan2(gy.fElements,gx.fElements)*2))++";
354 fstring+="(side*cos(atan2(gy.fElements,gx.fElements)*2))++";
355 fstring+="(side*sin(atan2(gy.fElements,gx.fElements)*3))++";
356 fstring+="(side*cos(atan2(gy.fElements,gx.fElements)*3))++";
358 fstring+="(side*sin(atan2(gy.fElements,gx.fElements)*2))*lxr++";
359 fstring+="(side*cos(atan2(gy.fElements,gx.fElements)*2))*lxr++";
360 fstring+="(side*sin(atan2(gy.fElements,gx.fElements)*3))*lxr++";
361 fstring+="(side*cos(atan2(gy.fElements,gx.fElements)*3))*lxr++";
363 TString *fit2 =stat.FitPlane(tree,"dt",fstring.Data(),"cutF&&abs(dt-f0)<0.7&&cutCE",chi2,npoints,vecFit2,mat,0.9);
364 strFit2 = fit2->Data();
365 printf("Fit2\t end\n");
366 printf("Chi2/npoints\t=%f\n",TMath::Sqrt(chi2/npoints));
367 fit2->Tokenize("++")->Print();
368 tree->SetAlias("f2",fit2->Data());
372 calPad0 = makePad->GetCalPad("f0","1", "ffit0");
373 calPad1 = makePad->GetCalPad("f1","1", "ffit1");
374 calPad2 = makePad->GetCalPad("f2","1", "ffit2");
375 calPadInOut = makePad->GetCalPad("fInOut","1", "fInOut");
376 calPadLX = makePad->GetCalPad("fLX","1", "fLX");
377 calPadTL = makePad->GetCalPad("fTL","1", "fTL");
378 //calPadQ = makePad->GetCalPad("fQ","1", "fQ");
379 calPadGXY = makePad->GetCalPad("fGXY","1", "fGXY");
380 calPadOff = makePad->GetCalPad("fOff","1", "fOff");
389 TObjArray * array = AliTPCCalibViewerGUI::ShowGUI(fname);
390 viewer = (AliTPCCalibViewerGUI*)array->At(0);
391 makePad = viewer->GetViewer();
392 tree = viewer->GetViewer()->GetTree();
403 // transform the input data to the fit format
406 makePad = new AliTPCCalibViewer(inname);
407 tree = makePad->GetTree();
412 printf("-1. MaQkeFitPulser\n");
418 printf("0. GetCalPads\n");
420 calPadCor = makePad->GetCalPad("tcor","1", "tcor");
421 calPadOut = makePad->GetCalPad("1","!((cutA||cutC)&&abs(ly.fElements/lx.fElements)<0.16)", "out");
422 calPadIn = makePad->GetCalPad("dt","1","timeIn");
424 calPadF1 = calPadIn->GlobalFit("timeF1", calPadOut,kFALSE,0,0,0.8); // use robust fit
425 calPadQIn = makePad->GetCalPad("qa*(sector%36<18)+qc*(sector%36>17)","1","qIn");
428 printf("1. Create outlyer map\n");
432 // Update outlyer maps
434 for (Int_t isector=0;isector<72; isector++){
435 for (UInt_t ich=0;ich<calPadIn->GetCalROC(isector)->GetNchannels();ich++){
436 Float_t val0= calPadIn->GetCalROC(isector)->GetValue(ich);
437 Float_t val1= calPadF1->GetCalROC(isector)->GetValue(ich);
438 if (TMath::Abs(val0-val1)>tThr) calPadOut->GetCalROC(isector)->SetValue(ich,1);
441 printf("2. Fit sector parabolas\n");
444 calPadF1 = calPadIn->GlobalFit("timeF1", calPadOut,kFALSE,0,0,0.9);
445 calPadF2 = calPadIn->GlobalFit("timeF2", calPadOut,kFALSE,1,0,0.9);
446 calPadQF1 = calPadQIn->GlobalFit("qF1", calPadOut,kFALSE,1);
447 calPadQF2 = calPadQIn->GlobalFit("qF2", calPadOut,kFALSE,1);
450 // Update outlyer maps
452 printf("3. Update outlyer map\n");
453 for (Int_t isector=0;isector<72; isector++){
454 for (UInt_t ich=0;ich<calPadIn->GetCalROC(isector)->GetNchannels();ich++){
455 Float_t val0= calPadQIn->GetCalROC(isector)->GetValue(ich);
456 Float_t val1= calPadQF2->GetCalROC(isector)->GetValue(ich);
460 if (TMath::Abs(val0/val1)<qThr0) calPadOut->GetCalROC(isector)->SetValue(ich,1);
461 if (TMath::Abs(val0/val1)>qThr1) calPadOut->GetCalROC(isector)->SetValue(ich,1);
464 printf("4. Redo fit of the of parabola \n");
469 AliTPCCalPad *calPadInCor = makePad->GetCalPad("dt","1","timeIn");
470 calPadF1 = calPadInCor->GlobalFit("timeF1", calPadOut,kFALSE,0,0.9);
471 calPadF2 = calPadInCor->GlobalFit("timeF2", calPadOut,kFALSE,1,0.9);
472 calPadQF1 = calPadQIn->GlobalFit("qF1", calPadOut,kFALSE,0,0,0.9);
473 calPadQF2 = calPadQIn->GlobalFit("qF2", calPadOut,kFALSE,1,0,0.9);
484 calPadIn = (AliTPCCalPad*)f.Get("timeIn"); // original time pad
485 calPadF1 = (AliTPCCalPad*)f.Get("timeF1"); // original time pad - fit plane
486 calPadF2 = (AliTPCCalPad*)f.Get("timeF2"); // original time pad - fit parabola
488 calPadQIn = (AliTPCCalPad*)f.Get("qIn"); // original time pad
489 calPadQF1 = (AliTPCCalPad*)f.Get("qF1"); // original time pad - fit plane
490 calPadQF2 = (AliTPCCalPad*)f.Get("qF2"); // original time pad - fit parabola
492 calPadCor = (AliTPCCalPad*)f.Get("tcor"); // base correction CalPad
493 calPadOut = (AliTPCCalPad*)f.Get("out"); // outlyer CalPad
495 calPad0 = (AliTPCCalPad*)f.Get("ffit0"); // global fit 0 - base
496 calPad1 = (AliTPCCalPad*)f.Get("ffit1"); // global fit 1 - common behavior rotation -A-C
497 calPad2 = (AliTPCCalPad*)f.Get("ffit2"); // gloabl fit 2 - CE missalign rotation A-C
498 calPadInOut = (AliTPCCalPad*)f.Get("fInOut");// misaalign in-out
499 calPadLX = (AliTPCCalPad*)f.Get("fLX"); // local x missalign
500 calPadTL = (AliTPCCalPad*)f.Get("fTL"); // local y/x missalign
501 calPadQ = (AliTPCCalPad*)f.Get("fQ"); // time (q) correction
502 calPadGXY = (AliTPCCalPad*)f.Get("fGXY"); // global XY missalign (drift velocity grad)
503 calPadOff = (AliTPCCalPad*)f.Get("fOff"); // normalization offset fit
504 calPadRes = (AliTPCCalPad*)f.Get("Result"); //resulting calibration
511 TFile * fstore = new TFile(oname,"recreate");
512 if (calPadIn) calPadIn->Write("timeIn"); // original time pad
513 if (calPadF1) calPadF1->Write("timeF1"); // original time pad - fit plane
514 if (calPadF2) calPadF2->Write("timeF2"); // original time pad - fit parabola
516 if (calPadQIn) calPadQIn->Write("qIn"); // original time pad
517 if (calPadQF1) calPadQF1->Write("qF1"); // original time pad - fit plane
518 if (calPadQF2) calPadQF2->Write("qF2"); // original time pad - fit parabola
520 if (calPadCor) calPadCor->Write("tcor"); // base correction CalPad
521 if (calPadOut) calPadOut->Write("out"); // outlyer CalPad
523 if (calPad0) calPad0->Write("ffit0"); // global fit 0 - base
524 if (calPad1) calPad1->Write("ffit1"); // global fit 1 - common behavior rotation -A-C
525 if (calPad2) calPad2->Write("ffit2"); // gloabl fit 2 - CE missalign rotation A-C
526 if (calPadInOut)calPadInOut->Write("fInOut"); // misaalign in-out
527 if (calPadLX) calPadLX->Write("fLX"); // local x missalign
528 if (calPadTL) calPadTL->Write("fTL"); // local y/x missalign
529 if (calPadQ) calPadQ->Write("fQ"); // time (q) correction
530 if (calPadGXY) calPadGXY->Write("fGXY"); // global XY missalign (drift velocity grad)
531 if (calPadOff) calPadOff->Write("fOff"); // normalization offset fit
532 if (calPadRes) calPadRes->Write("Result"); //resulting calibration
541 AliTPCPreprocessorOnline * preprocesor = new AliTPCPreprocessorOnline;
543 if (calPadIn) preprocesor->AddComponent(calPadIn->Clone());
544 if (calPadF1) preprocesor->AddComponent(calPadF1->Clone());
545 if (calPadF2) preprocesor->AddComponent(calPadF2->Clone());
547 if (calPadQIn) preprocesor->AddComponent(calPadQIn->Clone());
548 if (calPadQF1) preprocesor->AddComponent(calPadQF1->Clone());
549 if (calPadQF2) preprocesor->AddComponent(calPadQF2->Clone());
551 if (calPadCor) preprocesor->AddComponent(calPadCor->Clone());
552 if (calPadOut) preprocesor->AddComponent(calPadOut->Clone());
554 if (calPad0) preprocesor->AddComponent(calPad0->Clone());
555 if (calPad1) preprocesor->AddComponent(calPad1->Clone());
556 if (calPad2) preprocesor->AddComponent(calPad2->Clone());
557 if (calPadInOut)preprocesor->AddComponent(calPadInOut->Clone());
558 if (calPadLX) preprocesor->AddComponent(calPadLX->Clone());
559 if (calPadTL) preprocesor->AddComponent(calPadTL->Clone());
560 if (calPadQ) preprocesor->AddComponent(calPadQ->Clone());
561 if (calPadGXY) preprocesor->AddComponent(calPadGXY->Clone());
562 if (calPadOff) preprocesor->AddComponent(calPadOff->Clone());
563 if (calPadRes) preprocesor->AddComponent(calPadRes->Clone());
564 preprocesor->DumpToFile(fname);
571 // Define variables and selection of outliers - for user defined tree
573 tree->SetAlias("tcor",tcor.Data()); // correction variable
574 tree->SetAlias("ta",taside+".fElements");
575 tree->SetAlias("tc",tcside+".fElements");
576 tree->SetAlias("qa",qaside+".fElements");
577 tree->SetAlias("qc",qcside+".fElements");
578 tree->SetAlias("ra",raside+".fElements");
579 tree->SetAlias("rc",rcside+".fElements");
580 tree->SetAlias("side","1-(sector%36>17)*2");
581 tree->SetAlias("dt","(ta)*(sector%36<18)+(tc)*(sector%36>17)+tcor");
582 tree->SetAlias("cutA","qa>30&&qa<400&&abs(ta)<2&&ra>0.5&&ra<2");
583 tree->SetAlias("cutC","qc>30&&qc<400&&abs(tc)<2&&rc>0.5&&rc<2");
584 tree->SetAlias("cutF","((row.fElements+pad.fElements+sector)%19==0)"); // use just part of the statistic - 5 %
585 tree->SetAlias("cutCE","V.out.fElements");
589 tree->SetAlias("inn","sector<36");
590 tree->SetAlias("gxr","(gx.fElements/250.)"); //
591 tree->SetAlias("gyr","(gy.fElements/250.)"); //
592 tree->SetAlias("lxr","(lx.fElements-133.41)/250.");
593 tree->SetAlias("qp","((sector%36<18)*sqrt(qa)/10.+(sector%36>17)*sqrt(qc)/10.)"); //
594 tree->SetAlias("tl","(ly.fElements/lx.fElements)/0.17");
600 // Define variables and selection of outliers -for default usage
602 tree->SetAlias("tcor","tcor.fElements"); // correction variable
603 tree->SetAlias("side","1-(sector%36>17)*2");
604 tree->SetAlias("dt","timeIn.fElements");
606 tree->SetAlias("cutA","out.fElements==1");
607 tree->SetAlias("cutC","out.fElements==1");
608 tree->SetAlias("cutF","((row.fElements+pad.fElements+sector.fElements)%19==0)"); // use just part of the statistic - 5 %
610 tree->SetAlias("cutCE","out.fElements<0.5");
614 tree->SetAlias("inn","sector<36");
615 tree->SetAlias("gxr","(gx.fElements/250.)"); //
616 tree->SetAlias("gyr","(gy.fElements/250.)"); //
617 tree->SetAlias("lxr","(lx.fElements-133.41)/250.");
618 tree->SetAlias("qp","(sqrt(qIn.fElements)/10.+(out.fElements>0.5))"); //
619 tree->SetAlias("tl","(ly.fElements/lx.fElements)/0.17");
626 // make final calibration
628 AliTPCCalPad * calPadRes0 =new AliTPCCalPad(*calPadIn);
629 calPadRes0->Add(calPad2,-1); // remove global fit
630 calPadRes = calPadRes0->GlobalFit("Result", calPadOut,kTRUE,1,0.9);
634 Float_t tlmedian = calPadTL->GetMedian();
635 for (Int_t isector=0;isector<72; isector++){
636 for (UInt_t ich=0;ich<calPadIn->GetCalROC(isector)->GetNchannels();ich++){
639 Float_t val0 = calPadRes->GetCalROC(isector)->GetValue(ich);
640 if (TMath::Abs(val0)>0.5) calPadRes->GetCalROC(isector)->SetValue(ich,0);
641 Float_t tl = calPadTL->GetCalROC(isector)->GetValue(ich);
642 Float_t inOut = calPadInOut->GetCalROC(isector)->GetValue(ich);
643 calPadRes->GetCalROC(isector)->SetValue(ich,calPadRes->GetCalROC(isector)->GetValue(ich)+tl+inOut);
648 calPadRes->Add(calPadCor,-1); // remove back correction (e.g Pulser time 0)
654 void RebuildCE(char *finname, char *pulname){
656 // Transformation from the CE to the visualization-analisys output
658 // finname = CE_Vscan_Run_61684-50_170.root;
660 AliTPCCalibCE * ce = (AliTPCCalibCE*)f.Get("AliTPCCalibCE");
662 AliTPCCalPad *padtime = new AliTPCCalPad((TObjArray*)ce->GetCalPadT0());
663 AliTPCCalPad *padRMS = new AliTPCCalPad((TObjArray*)ce->GetCalPadRMS());
664 AliTPCCalPad *padq = new AliTPCCalPad((TObjArray*)ce->GetCalPadQ());
665 padtime->SetName("CE_T");
666 padRMS->SetName("CE_RMS");
667 padq->SetName("CE_Q");
670 AliTPCCalibPulser *pul = (AliTPCCalibPulser*)f2.Get("AliTPCCalibPulser");
671 AliTPCCalPad *pultime = new AliTPCCalPad((TObjArray*)pul->GetCalPadT0());
672 AliTPCCalPad *pulRMS = new AliTPCCalPad((TObjArray*)pul->GetCalPadRMS());
673 AliTPCCalPad *pulq = new AliTPCCalPad((TObjArray*)pul->GetCalPadQ());
674 pultime->SetName("P_T");
675 pulRMS->SetName("P_RMS");
676 pulq->SetName("P_Q");
679 AliTPCPreprocessorOnline * preprocesor = new AliTPCPreprocessorOnline;
680 preprocesor->AddComponent(padtime);
681 preprocesor->AddComponent(padq);
682 preprocesor->AddComponent(padRMS);
683 preprocesor->AddComponent(pultime);
684 preprocesor->AddComponent(pulq);
685 preprocesor->AddComponent(pulRMS);
686 preprocesor->DumpToFile(inname);
690 void AnalyzeLaserCE(){
691 RebuildCE("ce.root","pul.root");
705 void AddFiles(char *list){
707 // prepare viewer for data sets
710 finput.open(list, ios_base::in);
715 elist = new TEntryList("elist","elist");
716 tree->Draw(Form(">>elist%d",counter),"1","entrylist");
718 while(finput.good()) {
719 finput >> currentFile;
720 printf("Getting file%s\n",currentFile.Data());
721 TFile * fC = new TFile(currentFile.Data());
722 TTree * treeC = (TTree*) fC->Get("calPads");
723 makePad->AddFriend(treeC,Form("T%d",counter));
725 TString cutStrF1=Form("abs(T%d.timeIn.fElements-T%d.timeF1.fElements)<1.5",counter,counter);
726 TEntryList *celist = new TEntryList(Form("listF1%d",counter),Form("listF1%d",counter));
727 tree->Draw(Form(">>listF1%d",counter),cutStrF1,"entrylist");
728 if (celist->GetN()>0) elist->Add(celist);
734 //void AnalyzeSectors(){
736 // // Analyze sector by sector
738 // TTreeSRedirector cstream("fits.root");