don't sort clusters after local reco, do this in AliITSUTrackerGlo
[u/mrichter/AliRoot.git] / ITS / UPGRADE / testDetector.C
1
2 // You have to load the class before ... ;-)
3 // .L DetectorK.cxx++
4
5 void standardPlots() {
6
7   DetectorK its("ALICE","ITS");
8
9   its.MakeAliceCurrent(0,0); 
10
11   // its.SetRadius("bpipe",2.1);
12   // its.AddLayer("spd0",2.2,0.001,0.0012,0.0012);
13
14   its.PrintLayout();
15   its.SolveViaBilloir();
16  
17   its.MakeStandardPlots(0,2);
18
19   its.AddTPC(0.1,0.1);
20   its.SolveViaBilloir();
21  
22   its.MakeStandardPlots(1,1);
23
24 }
25
26
27 void particleDependendResolution() { 
28 // particle dependency on resolution
29
30   // .L Detector.cxx++
31
32   Detector its("ALICE","ITS");
33
34   its.MakeAliceCurrent(); 
35   its.PrintLayout();
36   its.PlotLayout();
37
38   its.SolveViaBilloir(0);
39  
40   its.SetRadius("bpipe",2.1);
41   its.AddLayer("spd0",2.2,0.001,0.0012,0.0012);
42
43   TCanvas *c1 = new TCanvas("c1","c1");
44  
45   c1->Divide(2,1);
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();
50
51   
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");
55
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");
60  
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");
65
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");
70
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");
75
76
77
78 }