AliRoot Core  v5-06-30 (35d6c57)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AnalyzeLaser.C
Go to the documentation of this file.
1 
3 
4 gSystem->Load("libSTAT");
5 TStatToolkit stat;
6 Int_t npoints;
7 Double_t chi2;
8 TVectorD vec;
9 TMatrixD mat;
10 
11 
14 TTree * tree = viewer->GetViewer()->GetTree();
15 TFile fp("/data/calib/CalibTreePulser_run33834_Cside.root");
16 tree->AddFriend(treePulser,"P.");
17 
18 
19 tree->SetAlias("dt","(sector%36>30)*2+(sector<36)*0.3");
20 tree->SetAlias("T","T0_100_220.fElements-P..StandardTime0.fElements");
21 tree->SetAlias("Tm","T0_100_220_Median.fElements");
22 tree->SetAlias("Q","Q_100_220.fElements");
23 tree->SetAlias("Qm","Q_100_220_Median.fElements");
24 
25 tree->SetAlias("Qcut","abs(Q/Qm-2)<1.4&&Q>6&&Q<200");
26 tree->SetAlias("Tcut","abs(T-Tm)<2");
27 
28 
29 
30 TString strSector="";
31 {
32  for (Int_t isec=54;isec<71;isec+=1){
33  if (isec!=64) {
34  strSector+="(sector==";
35  strSector+=isec;
36  strSector+=")++";
37  strSector+="(lx.fElements-195.)*(sector==";
38  strSector+=isec;
39  strSector+=")++";
40  strSector+="((lx.fElements-195)^2)*(sector==";
41  strSector+=isec;
42  strSector+=")++";
43  }
44  }
45 }
46 
47 
48 TCut cutA="Tcut&&Qcut&&sector%36>17";
49 
50 TString *strFitG =stat.FitPlane(tree,"T+dt","gx.fElements++gy.fElements",cutA,chi2,npoints,vec,mat);
51 
52 TString *strFitGL =stat.FitPlane(tree,"T+dt","lx.fElements++ly.fElements++gx.fElements++gy.fElements",cutA,chi2,npoints,vec,mat);
53 
54 TString *strFitGL2 =stat.FitPlane(tree,"T+dt","lx.fElements++ly.fElements++gx.fElements++gy.fElements++lx.fElements^2++ly.fElements^2",cutA,chi2,npoints,vec,mat);
55 
56 TString *strFitGLA =stat.FitPlane(tree,"T+dt",strSector+"lx.fElements++ly.fElements++gx.fElements++gy.fElements++lx.fElements^2++ly.fElements^2",cutA,chi2,npoints,vec,mat);
57 
58 
59 tree->SetAlias("tfitG",strFitG->Data())
60 tree->SetAlias("tfitGL",strFitGL->Data())
61 tree->SetAlias("tfitGL2",strFitGL2->Data())
62 tree->SetAlias("tfitGLA",strFitGLA->Data())
63 
64 
TString * strFitGLA
Definition: AnalyzeLaser.C:56
#define TObjArray
TVectorD vec
Definition: AnalyzeLaser.C:8
GUI for the AliTPCCalibViewer used for the calibration monitor All functionalities of the AliTPCCalib...
TMatrixD mat
Definition: AnalyzeLaser.C:9
TFile fp("/data/calib/CalibTreePulser_run33834_Cside.root")
TCut cutA
Definition: AnalyzeLaser.C:48
AliTPCCalibViewer * GetViewer()
AliTPCCalibViewerGUI * viewer
Definition: AnalyzeLaser.C:13
TStatToolkit stat
Definition: AnalyzeLaser.C:5
TObjArray * array
Definition: AnalyzeLaser.C:12
Int_t npoints
Definition: AnalyzeLaser.C:6
TString strSector
Definition: AnalyzeLaser.C:30
Double_t chi2
Definition: AnalyzeLaser.C:7
static TObjArray * ShowGUI(const char *fileName=0)
TTree * tree
Definition: AnalyzeLaser.C:14
TString * strFitG
Definition: AnalyzeLaser.C:50
TString * strFitGL
Definition: AnalyzeLaser.C:52
TTree * GetTree() const
TString * strFitGL2
Definition: AnalyzeLaser.C:54