5 gSystem->Load("libSTAT.so");
7 .L $ALICE_ROOT/TPC/fastSimul/AliTPCclusterFast.cxx+
8 .L $ALICE_ROOT/TPC/fastSimul/AliTPCclusterAnalysis.C
19 TChain * treeCluster=0;
21 AliTPCfastTrack * track =0;
23 void LoadTree(const char* fname="cluterSimul.root"){
27 treeCluster = new TChain("simul","simul");
28 treeCluster->AddFile(fname);
31 void LoadTrack(const char* fname="trackerSimul.root"){
35 treeTrack = new TChain("simulTrack","simulTrack");
36 treeTrack->AddFile(fname);
46 void MakeQNormalization(Int_t maxPoints){
48 // Normalize Q to the diffusion and angular effect
58 fstringQ+="abs(fAngleY)++";
59 fstringQ+="abs(fAngleZ)++";
63 strQmax = toolkit.FitPlane(treeCluster,"s.GetQmax(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 strQmax->Tokenize("++")->Print();
66 treeCluster->SetAlias("qMaxCorr",strQmax->Data());
68 strQtot = toolkit.FitPlane(treeCluster,"s.GetQtot(0.33,2.5,1,1,1)/s.fNtot",fstringQ.Data(),"", chi2,npoints,fitParam,covMatrix,-1,0,50000);
69 printf("Qmax norm\t%f\n",TMath::Sqrt(chi2/npoints));
70 strQtot->Tokenize("++")->Print();
71 treeCluster->SetAlias("qTotCorr",strQtot->Data());
80 // Correction tests - normalization of response functions
84 TF2 f2("f2","AliTPCclusterFast::GaussConvolution(x,y,2,0.5,0.5,0.5)",-3,3,-3,3);
86 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);
92 Double_t testSumGaus(Float_t k0,Float_t k1, Float_t s0, Float_t s1){
93 TF2 f2("f2",Form("AliTPCclusterFast::GaussConvolution(x,y,%f,%f,%f,%f)",k0,k1,s0,s1),-2,2,-2,2);
95 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);
97 printf("%f\t%f\t%f\t%f\t%f\n",k0,k1,s0,s1, sumg);
101 Double_t testSumExp(Float_t s0, Float_t k1){
102 TF1 f1("f1",Form("AliTPCclusterFast::GaussExpConvolution(x,%f,%f)",s0,k1),-2,2);
104 for (Float_t x=-5; x<5; x+=0.2) sumg+=f1.Eval(x);
106 printf("%f\t%f\t%f\t%f\t%f\n",s0,k1, sumg);