4 TFile *f = new TFile("SSD_his.root");
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);
13 ///////////////////////// Ntuple analysis ///////////////////////////////
15 // ntuple is created inside the hit loop for the hits corresponding to the
18 // ntuple1 is created after a finish of the hit loop if one or more hits
19 // correspond to the recpoint;
21 // ntuple2 is created befor the hit loop for all recpoints;
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 // -------------------------------------------------------------------------
46 gPad->SetFillColor(33);
47 ntuple->SetFillColor(42);
48 ntuple->Draw("dx","lay == 5 && hitprim == 1&&abs(dx)<200");
50 gPad->SetFillColor(33);
51 ntuple->SetFillColor(46);
52 ntuple->Draw("dz","lay == 5 && hitprim == 1&&abs(dz)<5000");
54 gPad->SetFillColor(33);
55 ntuple->SetFillColor(42);
56 ntuple->Draw("dx","lay == 6 && hitprim == 1&&abs(dx)<200");
58 gPad->SetFillColor(33);
59 ntuple->SetFillColor(46);
60 ntuple->Draw("dz","lay == 6 && hitprim == 1&&abs(dz)<5000");
65 gPad->SetFillColor(33);
66 ntuple1->SetFillColor(42);
67 ntuple1->Draw("nxP","lay == 5&&noverprim>=0");
69 gPad->SetFillColor(33);
70 ntuple1->SetFillColor(46);
71 ntuple1->Draw("nxN","lay == 5&&noverprim>=0");
73 gPad->SetFillColor(33);
74 ntuple1->SetFillColor(42);
75 ntuple1->Draw("nxP","lay == 6 && noverprim>=0");
77 gPad->SetFillColor(33);
78 ntuple1->SetFillColor(46);
79 ntuple1->Draw("nxN","lay == 6 && noverprim>=0");
84 gPad->SetFillColor(33);
85 ntuple1->SetFillColor(42);
86 ntuple1->Draw("qclP","lay == 5&&noverprim>=0");
88 gPad->SetFillColor(33);
89 ntuple1->SetFillColor(42);
90 ntuple1->Draw("qclN","lay == 5&&noverprim>=0");
92 gPad->SetFillColor(33);
93 ntuple1->SetFillColor(46);
94 ntuple1->Draw("qclP","lay == 6&&noverprim>=0");
96 gPad->SetFillColor(33);
97 ntuple1->SetFillColor(46);
98 ntuple1->Draw("qclN","lay == 6&&noverprim>=0");
102 ///////////////////// Histogramm/ntuple analysis ////////////////////////
106 gPad->SetFillColor(33);
107 ntuple1->SetFillColor(42);
108 ntuple1->Draw("qclP","noverprim>=0");
110 gPad->SetFillColor(33);
111 ntuple1->SetFillColor(42);
112 ntuple1->Draw("qclN","noverprim>=0");
114 gPad->SetFillColor(33);
115 adcPadcN5cut->SetFillColor(42);
116 adcPadcN5cut->Draw();
118 gPad->SetFillColor(33);
119 adcPadcN6cut->SetFillColor(46);
120 adcPadcN6cut->Draw();
126 c2->Print("ssd_res.ps");
131 c2->Print("spd_clsize.ps");