]>
Commit | Line | Data |
---|---|---|
6b8f55ce | 1 | void SSD_ntuple() |
2 | { | |
3 | ||
4 | TFile *f = new TFile("SSD_his.root"); | |
5 | ||
6 | //gStyle->SetOptStat(1111111); | |
7 | //gStyle->SetOptLogy(); | |
8 | //TCanvas *c1 = new TCanvas("c1","SPD clusters",400,10,600,700); | |
9 | TCanvas *c2 = new TCanvas("c2","SPD clusters",400,10,600,700); | |
10 | //c1->Divide(2,2); | |
11 | c2->Divide(2,2); | |
12 | ||
13 | ///////////////////////// Ntuple analysis /////////////////////////////// | |
14 | ||
15 | // ntuple is created inside the hit loop for the hits corresponding to the | |
16 | // recpoint; | |
17 | ||
18 | // ntuple1 is created after a finish of the hit loop if one or more hits | |
19 | // correspond to the recpoint; | |
20 | ||
21 | // ntuple2 is created befor the hit loop for all recpoints; | |
22 | ||
23 | // ----------------------------------------------------------------------- | |
24 | // lay - number of ITS layer; | |
25 | // lad - number of ITS ladder; | |
26 | // det - number of ITS detector; | |
27 | // nxP/N - cluster size in the r*phi(x) direction for P/N sides; | |
28 | // hitprim - primary particle(hit) flag ( = 1 for primery particle); | |
29 | // x - x local coordinate in mm; | |
30 | // z - z local coordinate in mm; | |
31 | // dx - difference of hit(mediate) and reconstructed (from cluster) | |
32 | // coordinates in r*phi(x) direction in microns; | |
33 | // dz - difference of hit(mediate) and reconstructed (from cluster) | |
34 | // coordinates in z direction in microns; | |
35 | // noverlaps - number of particles overlapping in one cluster; | |
36 | // noverprim - number of primary particles overlapping in one cluster; | |
37 | // qclP/N - cluster signals in ADC normalized to the path length in Si | |
38 | // qrec - recpoint signal (maximum from qclP and qclN) | |
39 | // pmod - particle momentum at the vertex in MeV/c | |
40 | // partcode - particle code | |
41 | // ------------------------------------------------------------------------- | |
42 | ||
43 | ||
44 | ||
45 | c2->cd(1); | |
46 | gPad->SetFillColor(33); | |
47 | ntuple->SetFillColor(42); | |
48 | ntuple->Draw("dx","lay == 5 && hitprim == 1&&abs(dx)<200"); | |
49 | c2->cd(2); | |
50 | gPad->SetFillColor(33); | |
51 | ntuple->SetFillColor(46); | |
52 | ntuple->Draw("dz","lay == 5 && hitprim == 1&&abs(dz)<5000"); | |
53 | c2->cd(3); | |
54 | gPad->SetFillColor(33); | |
55 | ntuple->SetFillColor(42); | |
56 | ntuple->Draw("dx","lay == 6 && hitprim == 1&&abs(dx)<200"); | |
57 | c2->cd(4); | |
58 | gPad->SetFillColor(33); | |
59 | ntuple->SetFillColor(46); | |
60 | ntuple->Draw("dz","lay == 6 && hitprim == 1&&abs(dz)<5000"); | |
61 | ||
62 | ||
0315d466 | 63 | /* |
6b8f55ce | 64 | c2->cd(1); |
65 | gPad->SetFillColor(33); | |
66 | ntuple1->SetFillColor(42); | |
67 | ntuple1->Draw("nxP","lay == 5&&noverprim>=0"); | |
68 | c2->cd(2); | |
69 | gPad->SetFillColor(33); | |
70 | ntuple1->SetFillColor(46); | |
71 | ntuple1->Draw("nxN","lay == 5&&noverprim>=0"); | |
72 | c2->cd(3); | |
73 | gPad->SetFillColor(33); | |
74 | ntuple1->SetFillColor(42); | |
75 | ntuple1->Draw("nxP","lay == 6 && noverprim>=0"); | |
76 | c2->cd(4); | |
77 | gPad->SetFillColor(33); | |
78 | ntuple1->SetFillColor(46); | |
79 | ntuple1->Draw("nxN","lay == 6 && noverprim>=0"); | |
0315d466 | 80 | */ |
6b8f55ce | 81 | |
82 | /* | |
83 | c2->cd(1); | |
84 | gPad->SetFillColor(33); | |
85 | ntuple1->SetFillColor(42); | |
86 | ntuple1->Draw("qclP","lay == 5&&noverprim>=0"); | |
87 | c2->cd(2); | |
88 | gPad->SetFillColor(33); | |
89 | ntuple1->SetFillColor(42); | |
90 | ntuple1->Draw("qclN","lay == 5&&noverprim>=0"); | |
91 | c2->cd(3); | |
92 | gPad->SetFillColor(33); | |
93 | ntuple1->SetFillColor(46); | |
94 | ntuple1->Draw("qclP","lay == 6&&noverprim>=0"); | |
95 | c2->cd(4); | |
96 | gPad->SetFillColor(33); | |
97 | ntuple1->SetFillColor(46); | |
98 | ntuple1->Draw("qclN","lay == 6&&noverprim>=0"); | |
99 | */ | |
100 | ||
101 | ||
102 | ///////////////////// Histogramm/ntuple analysis //////////////////////// | |
103 | ||
104 | /* | |
105 | c2->cd(1); | |
106 | gPad->SetFillColor(33); | |
0315d466 | 107 | ntuple1->SetFillColor(42); |
108 | ntuple1->Draw("qclP","noverprim>=0"); | |
6b8f55ce | 109 | c2->cd(2); |
110 | gPad->SetFillColor(33); | |
0315d466 | 111 | ntuple1->SetFillColor(42); |
112 | ntuple1->Draw("qclN","noverprim>=0"); | |
6b8f55ce | 113 | c2->cd(3); |
114 | gPad->SetFillColor(33); | |
115 | adcPadcN5cut->SetFillColor(42); | |
116 | adcPadcN5cut->Draw(); | |
117 | c2->cd(4); | |
118 | gPad->SetFillColor(33); | |
119 | adcPadcN6cut->SetFillColor(46); | |
120 | adcPadcN6cut->Draw(); | |
121 | */ | |
122 | ||
123 | ||
124 | /* | |
125 | c2->Draw(); | |
126 | c2->Print("ssd_res.ps"); | |
127 | */ | |
128 | ||
129 | /* | |
130 | c2->Draw(); | |
131 | c2->Print("spd_clsize.ps"); | |
132 | */ | |
133 | ||
134 | } | |
135 | ||
136 | ||
137 | ||
138 | ||
139 | ||
140 | ||
141 | ||
142 | ||
143 |