AliRoot Core  a565103 (a565103)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliTPCclusterAnalysis.C
Go to the documentation of this file.
1 
18 TChain * treeCluster=0;
19 TChain * treeTrack=0;
20 AliTPCfastTrack * track =0;
21 
22 void LoadTree(const char* fname="cluterSimul.root"){
24 
25  treeCluster = new TChain("simul","simul");
26  treeCluster->AddFile(fname);
27 }
28 
29 void LoadTrack(const char* fname="trackerSimul.root"){
31 
32  treeTrack = new TChain("simulTrack","simulTrack");
33  treeTrack->AddFile(fname);
34  TFile f(fname);
35 
36 }
37 
38 
39 
40 
41 
42 
43 void MakeQNormalization(Int_t maxPoints){
45 
46  TStatToolkit toolkit;
47  Double_t chi2;
48  TVectorD fitParam;
49  TMatrixD covMatrix;
50  Int_t npoints;
51  TString fstringQ="";
52  fstringQ+="fDiff++";
53  fstringQ+="abs(fAngleY)++";
54  fstringQ+="abs(fAngleZ)++";
55  TString *strQmax=0;
56  TString *strQtot=0;
57  //
58  strQmax = toolkit.FitPlane(treeCluster,"s.GetQmax(0.33,2.5,1,1,1)/s.fNtot",fstringQ.Data(), "", chi2,npoints,fitParam,covMatrix,-1,0,50000);
59  printf("Qmax norm\t%f\n",TMath::Sqrt(chi2/npoints));
60  strQmax->Tokenize("++")->Print();
61  treeCluster->SetAlias("qMaxCorr",strQmax->Data());
62  //
63  strQtot = toolkit.FitPlane(treeCluster,"s.GetQtot(0.33,2.5,1,1,1)/s.fNtot",fstringQ.Data(),"", chi2,npoints,fitParam,covMatrix,-1,0,50000);
64  printf("Qmax norm\t%f\n",TMath::Sqrt(chi2/npoints));
65  strQtot->Tokenize("++")->Print();
66  treeCluster->SetAlias("qTotCorr",strQtot->Data());
67  //
68  //
69 }
70 
71 
72 
73 
74 //
75 // Correction tests - normalization of response functions
76 //
77 
78 void sum(){
79  TF2 f2("f2","AliTPCclusterFast::GaussConvolution(x,y,2,0.5,0.5,0.5)",-3,3,-3,3);
80  Float_t sumg=0;
81  for (Float_t x=-5; x<5; x+=0.5) for (Float_t y=-5; y<5; y+=0.5) sumg+=f2->Eval(x,y);
82  printf("%f\n", sumg);
83 }
84 
85 
86 
87 Double_t testSumGaus(Float_t k0,Float_t k1, Float_t s0, Float_t s1){
88  TF2 f2("f2",Form("AliTPCclusterFast::GaussConvolution(x,y,%f,%f,%f,%f)",k0,k1,s0,s1),-2,2,-2,2);
89  Float_t sumg=0;
90  for (Float_t x=-5; x<5; x+=0.2) for (Float_t y=-5; y<5; y+=0.2) sumg+=f2.Eval(x,y);
91  sumg*=0.2*0.2;
92  printf("%f\t%f\t%f\t%f\t%f\n",k0,k1,s0,s1, sumg);
93  return sumg;
94 }
95 
96 Double_t testSumExp(Float_t s0, Float_t k1){
97  TF1 f1("f1",Form("AliTPCclusterFast::GaussExpConvolution(x,%f,%f)",s0,k1),-2,2);
98  Float_t sumg=0;
99  for (Float_t x=-5; x<5; x+=0.2) sumg+=f1.Eval(x);
100  sumg*=0.2;
101  printf("%f\t%f\t%f\t%f\t%f\n",s0,k1, sumg);
102  return sumg;
103 }
104 
105 
106 
107 
printf("Chi2/npoints = %f\n", TMath::Sqrt(chi2/npoints))
void LoadTrack(const char *fname="trackerSimul.root")
TFile f("CalibObjects.root")
AliTPCfastTrack * track
npoints
Definition: driftITSTPC.C:85
TChain * treeTrack
void sum()
Double_t testSumGaus(Float_t k0, Float_t k1, Float_t s0, Float_t s1)
Double_t chi2
Definition: AnalyzeLaser.C:7
Double_t testSumExp(Float_t s0, Float_t k1)
void MakeQNormalization(Int_t maxPoints)
TStatToolkit toolkit
Definition: gainCalib.C:18
TChain * treeCluster
char * fname
void LoadTree(const char *fname="cluterSimul.root")