1 void eval(char *inFile)
5 // Connect the Root Galice file containing Geometry, Kine and Hits
7 file = new TFile(inFile);
9 gAlice = (AliRun*)file->Get("gAlice");
10 if (gAlice) printf("AliRun object found on file\n");
11 if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
14 AliTPC *TPC = (AliTPC*)gAlice->GetDetector("TPC");
16 AliTPCParam *param = (AliTPCParam*)file->Get("75x40_100x60");
18 TFile *nfile = new TFile("/prog/alice/data/exact_clusters.root","RECREATE");
21 //setup AliTPCClustersArray
22 AliTPCClustersArray * arr=new AliTPCClustersArray;
23 arr->SetClusterType("AliComplexCluster");
28 TPC->SetClustersArray(arr);
29 TPC->Hits2ExactClustersSector(1);
30 TPC->Hits2ExactClustersSector(37);
35 sprintf(treeName,"TreeCExact_%s",param->GetTitle());
36 TPC->GetClustersArray()->GetTree()->Write(treeName);
37 param->Write(param->GetTitle());
44 void calc(char *rootfile,char *digitsfile,char *cfile)
48 a = new AliL3Evaluate(rootfile,digitsfile,slice);
49 a->SetupSlow("tracks.raw","./");
50 ntuppel = (TNtuple*)a->EvaluatePoints(cfile);
52 f = new TFile("results.root","RECREATE");
61 gStyle->SetOptFit(0110);
64 TCanvas *c1 = new TCanvas("c1","",2);
65 TH1F *hist = new TH1F("hist","",100,-2,2);
67 ntuppel->Draw("resy>>hist","pt>1");
69 float rms = hist->GetRMS();
70 hist->SetXTitle("#delta_{T} [cm]");
72 TF1 *f1 = new TF1("f1","gaus",-rms,rms);