2 // You have to load the class before ... ;-)
7 DetectorK its("ALICE","ITS");
9 its.MakeAliceCurrent(0,0);
11 // its.SetRadius("bpipe",2.1);
12 // its.AddLayer("spd0",2.2,0.001,0.0012,0.0012);
15 its.SolveViaBilloir();
17 its.MakeStandardPlots(0,2);
20 its.SolveViaBilloir();
22 its.MakeStandardPlots(1,1);
27 void particleDependendResolution() {
28 // particle dependency on resolution
32 Detector its("ALICE","ITS");
34 its.MakeAliceCurrent();
38 its.SolveViaBilloir(0);
40 its.SetRadius("bpipe",2.1);
41 its.AddLayer("spd0",2.2,0.001,0.0012,0.0012);
43 TCanvas *c1 = new TCanvas("c1","c1");
46 c1->cd(1); gPad->SetGridx(); gPad->SetGridy();
47 gPad->SetLogx(); //gPad->SetLogy();
48 c1->cd(2); gPad->SetGridx(); gPad->SetGridy();
49 gPad->SetLogx(); //gPad->SetLogy();
52 // compare to telescope equation ?
53 // c1->cd(1); its.GetGraphPointingResolutionTeleEqu(0,1)->Draw("AC");
54 // c1->cd(2); its.GetGraphPointingResolutionTeleEqu(1,1)->Draw("AC");
56 its.SetParticleMass(0.140); // pion
57 its.SolveViaBilloir(0,0);
58 c1->cd(1); its.GetGraphPointingResolution(0,1)->Draw("AC");
59 c1->cd(2); its.GetGraphPointingResolution(1,1)->Draw("AC");
61 its.SetParticleMass(0.498); // kaon
62 its.SolveViaBilloir(0,0);
63 c1->cd(1); its.GetGraphPointingResolution(0,2)->Draw("C");
64 c1->cd(2); its.GetGraphPointingResolution(1,2)->Draw("C");
66 its.SetParticleMass(0.00051); // electron
67 its.SolveViaBilloir(0,0);
68 c1->cd(1); its.GetGraphPointingResolution(0,3)->Draw("C");
69 c1->cd(2); its.GetGraphPointingResolution(1,3)->Draw("C");
71 its.SetParticleMass(0.938); // proton
72 its.SolveViaBilloir(0,0);
73 c1->cd(1); its.GetGraphPointingResolution(0,4)->Draw("C");
74 c1->cd(2); its.GetGraphPointingResolution(1,4)->Draw("C");