1 void Eval(char *rootfile="")
5 // l.UnSet(AliL3Logger::kDebug);
6 // l.UnSet(AliL3Logger::kAll);
7 l.Set(AliL3Logger::kError);
12 e = new AliL3Evaluate(rootfile,slice);
13 e->SetupSlow("tracks.raw",".");
14 //e->SetupFast("tracks.raw","/nfs/david/subatom/alice/data/V3.04/fast/clusters/hg_8k_v0_s1-3_e0_cl.root",".");
16 ntuppel = (TNtuple*)e->EvaluatePoints();
19 file = new TFile("CFeval_nodeconv.root","RECREATE");
26 void plotPt(char *rootfile)
28 gStyle->SetStatColor(10);
29 gStyle->SetOptStat(0);
30 gStyle->SetOptFit(1100);
33 c = new TCanvas("c","",2);
37 f1 = new TFile(rootfile);
39 hist = new TH1F("hist","",50,-10,10);
41 fNtuppel->Draw("(pt_found-pt_gen)/pt_gen*100>>hist","nHits>30");
42 hist->GetXaxis()->SetTitle("%");
43 hist->GetYaxis()->SetTitle("Counts");
45 float rms = hist->GetRMS();
46 printf("Rms value : %f\n",rms);
48 TF1 *f = new TF1("f","gaus",-rms,rms);
54 void plot(char *rootfile)
57 gStyle->SetStatColor(10);
58 gStyle->SetOptStat(1);
61 file = new TFile(rootfile);
62 hist = new TH1F("hist","",100,-0.6,0.6);
66 can = new TCanvas("can","Residuals",900,600);
68 SetCanvasOptions(can);
70 //ntuppel->Draw("residual_trans>>hist","zHit < 50 && padrow > 55");//beta < 10*4.1515/180");
71 ntuppel->Draw("resy>>hist","ptgen > 1.0");
73 float rms = hist->GetRMS();
74 printf("Rms value : %f\n",rms);
76 TF1 *f = new TF1("f","gaus",-rms,rms);
79 hist->GetXaxis()->SetTitle("#delta_{T} [cm]");
80 hist->GetYaxis()->SetTitle("Counts");
82 f2 = new TFile("results_fast_oldparams.root");
83 hist2 = new TH1F("hist2","",100,-0.6,0.6);
86 ntuppel_fast->Draw("residual_trans>>hist2","nHits>100 && pt>1.0 && padrow > 0 && zHit < 50");//dipangle < 20*3.1415/180");