ITS/SSD_ntuple.C

   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","SSD 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");
  49 c2->cd(2);
  50 gPad->SetFillColor(33);
  51       ntuple->SetFillColor(46);
  52       ntuple->Draw("dz","lay == 5 && hitprim == 1");
  53 c2->cd(3);
  54 gPad->SetFillColor(33);
  55       ntuple->SetFillColor(42);
  56       ntuple->Draw("dx","lay == 6 && hitprim == 1");
  57 c2->cd(4);
  58 gPad->SetFillColor(33);
  59       ntuple->SetFillColor(46);
  60       ntuple->Draw("dz","lay == 6 && hitprim == 1");
  61 */
  62
  63 /*
  64 c2->cd(1);
  65 gPad->SetFillColor(33);
  66       ntuple->SetFillColor(42);
  67       ntuple->Draw("dx","lay == 5 && hitprim == 1&&ntrover==1");
  68 c2->cd(2);
  69 gPad->SetFillColor(33);
  70       ntuple->SetFillColor(46);
  71       ntuple->Draw("dz","lay == 5 && hitprim == 1&&ntrover==1");
  72 c2->cd(3);
  73 gPad->SetFillColor(33);
  74       ntuple->SetFillColor(42);
  75       ntuple->Draw("dx","lay == 6 && hitprim == 1&&ntrover==1");
  76 c2->cd(4);
  77 gPad->SetFillColor(33);
  78       ntuple->SetFillColor(46);
  79       ntuple->Draw("dz","lay == 6 && hitprim == 1&&ntrover==1");
  80 */
  81
  82
  83 /*
  84 c2->cd(1);
  85 gPad->SetFillColor(33);
  86       ntuple->SetFillColor(42);
  87       ntuple->Draw("dx","lay == 5 && hitprim == 1&&abs(dx)<200");
  88 c2->cd(2);
  89 gPad->SetFillColor(33);
  90       ntuple->SetFillColor(46);
  91       ntuple->Draw("dz","lay == 5 && hitprim == 1&&abs(dz)<5000");
  92 c2->cd(3);
  93 gPad->SetFillColor(33);
  94       ntuple->SetFillColor(42);
  95       ntuple->Draw("dx","lay == 6 && hitprim == 1&&abs(dx)<200");
  96 c2->cd(4);
  97 gPad->SetFillColor(33);
  98       ntuple->SetFillColor(46);
  99       ntuple->Draw("dz","lay == 6 && hitprim == 1&&abs(dz)<5000");
 100 */
 101
 102 /*
 103 c2->cd(1);
 104 gPad->SetFillColor(33);
 105       ntuple->SetFillColor(42);
 106       ntuple->Draw("dx","lay == 5 && hitprim == 1&&abs(dx)<200&&ntrover==1");
 107 c2->cd(2);
 108 gPad->SetFillColor(33);
 109       ntuple->SetFillColor(46);
 110       ntuple->Draw("dz","lay == 5 && hitprim == 1&&abs(dz)<5000&&ntrover==1");
 111 c2->cd(3);
 112 gPad->SetFillColor(33);
 113       ntuple->SetFillColor(42);
 114       ntuple->Draw("dx","lay == 6 && hitprim == 1&&abs(dx)<200&&ntrover==1");
 115 c2->cd(4);
 116 gPad->SetFillColor(33);
 117       ntuple->SetFillColor(46);
 118       ntuple->Draw("dz","lay == 6 && hitprim == 1&&abs(dz)<5000&&ntrover==1");
 119 */
 120
 121 /*
 122 c2->cd(1);
 123 gPad->SetFillColor(33);
 124       ntuple1->SetFillColor(42);
 125       ntuple1->Draw("nxP","lay == 5&&noverprim>=0");
 126 c2->cd(2);
 127 gPad->SetFillColor(33);
 128       ntuple1->SetFillColor(46);
 129       ntuple1->Draw("nxN","lay == 5&&noverprim>=0");
 130 c2->cd(3);
 131 gPad->SetFillColor(33);
 132       ntuple1->SetFillColor(42);
 133       ntuple1->Draw("nxP","lay == 6 && noverprim>=0");
 134 c2->cd(4);
 135 gPad->SetFillColor(33);
 136       ntuple1->SetFillColor(46);
 137       ntuple1->Draw("nxN","lay == 6 && noverprim>=0");
 138 */
 139
 140 /*
 141 c2->cd(1);
 142 gPad->SetFillColor(33);
 143       ntuple1->SetFillColor(42);
 144       ntuple1->Draw("dx","lay == 5&&noverprim>=0&&abs(dx)<50&&nxP==2&&nxN==2");
 145 c2->cd(2);
 146 gPad->SetFillColor(33);
 147       ntuple1->SetFillColor(46);
 148       ntuple1->Draw("dz","lay == 5&&noverprim>=0&&abs(dz)<2000&&nxP==2&&nxN==2");
 149 c2->cd(3);
 150 gPad->SetFillColor(33);
 151       ntuple1->SetFillColor(42);
 152       ntuple1->Draw("dx","lay == 6 && noverprim>=0&&abs(dx)<200");
 153 c2->cd(4);
 154 gPad->SetFillColor(33);
 155       ntuple1->SetFillColor(46);
 156       ntuple1->Draw("dz","lay == 6 && noverprim>=0&&abs(dz)<5000");
 157 */
 158
 159 /*
 160 c2->cd(1);
 161 gPad->SetFillColor(33);
 162       ntuple1->SetFillColor(42);
 163       ntuple1->Draw("dx","lay == 5&&noverprim>=0&&noverlaps==0&&qcut<0.18");
 164 c2->cd(2);
 165 gPad->SetFillColor(33);
 166       ntuple1->SetFillColor(46);
 167       ntuple1->Draw("dz","lay == 5&&noverprim>=0&&noverlaps==0&&qcut<0.18");
 168 c2->cd(3);
 169 gPad->SetFillColor(33);
 170       ntuple1->SetFillColor(42);
 171       ntuple1->Draw("dx","lay == 6 && noverprim>=0&&noverlaps==0&&qcut<0.18");
 172 c2->cd(4);
 173 gPad->SetFillColor(33);
 174       ntuple1->SetFillColor(46);
 175       ntuple1->Draw("dz","lay == 6 && noverprim>=0&&noverlaps==0&&qcut<0.18");
 176 */
 177
 178 /*
 179 c2->cd(1);
 180 gPad->SetFillColor(33);
 181       ntuple1->SetFillColor(42);
 182       ntuple1->Draw("dx","lay == 5&&noverprim>=0&&noverlaps==0&&abs(dx)<100&&qcut<0.18");
 183 c2->cd(2);
 184 gPad->SetFillColor(33);
 185       ntuple1->SetFillColor(46);
 186       ntuple1->Draw("dz","lay == 5&&noverprim>=0&&noverlaps==0&&abs(dz)<5000&&qcut<0.18");
 187 c2->cd(3);
 188 gPad->SetFillColor(33);
 189       ntuple1->SetFillColor(42);
 190       ntuple1->Draw("dx","lay == 6&&noverprim>=0&&noverlaps==0&&abs(dx)<100&&qcut<0.18");
 191 c2->cd(4);
 192 gPad->SetFillColor(33);
 193       ntuple1->SetFillColor(46);
 194       ntuple1->Draw("dz","lay == 6&&noverprim>=0&&noverlaps==0&&abs(dz)<5000&&qcut<0.18");
 195 */
 196
 197 /*
 198 c2->cd(1);
 199 gPad->SetFillColor(33);
 200       ntuple1->SetFillColor(42);
 201       ntuple1->Draw("noverprim","lay == 5&&noverprim>=-1");
 202 c2->cd(2);
 203 gPad->SetFillColor(33);
 204       ntuple1->SetFillColor(46);
 205       ntuple1->Draw("noverlaps","lay == 5&&noverprim>=-1");
 206 c2->cd(3);
 207 gPad->SetFillColor(33);
 208       ntuple1->SetFillColor(42);
 209       ntuple1->Draw("noverprim","lay == 6 && noverprim>=-1");
 210 c2->cd(4);
 211 gPad->SetFillColor(33);
 212       ntuple1->SetFillColor(46);
 213       ntuple1->Draw("noverlaps","lay == 6 && noverprim>=-1");
 214 */
 215
 216       /*
 217 c2->cd(1);
 218 gPad->SetFillColor(33);
 219       ntuple1->SetFillColor(42);
 220       ntuple1->Draw("ntrover","lay == 5&&noverprim>=0");
 221 c2->cd(2);
 222 gPad->SetFillColor(33);
 223       ntuple1->SetFillColor(46);
 224       ntuple1->Draw("noverlaps","lay == 5&&noverprim>=0");
 225 c2->cd(3);
 226 gPad->SetFillColor(33);
 227       ntuple1->SetFillColor(42);
 228       ntuple1->Draw("ntrover","lay == 6 && noverprim>=0");
 229 c2->cd(4);
 230 gPad->SetFillColor(33);
 231       ntuple1->SetFillColor(46);
 232       ntuple1->Draw("noverlaps","lay == 6 && noverprim>=0");
 233       */
 234
 235 /*
 236 c2->cd(1);
 237 gPad->SetFillColor(33);
 238       ntuple1->SetFillColor(42);
 239       ntuple1->Draw("qclP","lay == 5&&noverprim>=0");
 240 c2->cd(2);
 241 gPad->SetFillColor(33);
 242       ntuple1->SetFillColor(42);
 243       ntuple1->Draw("qclN","lay == 5&&noverprim>=0");
 244 c2->cd(3);
 245 gPad->SetFillColor(33);
 246       ntuple1->SetFillColor(46);
 247       ntuple1->Draw("qclP","lay == 6&&noverprim>=0");
 248 c2->cd(4);
 249 gPad->SetFillColor(33);
 250       ntuple1->SetFillColor(46);
 251       ntuple1->Draw("qclN","lay == 6&&noverprim>=0");
 252 */
 253
 254
 255 /////////////////////   Histogramm/ntuple  analysis  ////////////////////////
 256
 257
 258 c2->cd(1);
 259 gPad->SetFillColor(33);
 260       adcPadcN5all->SetFillColor(42);
 261       adcPadcN5all->Draw();
 262 c2->cd(2);
 263 gPad->SetFillColor(33);
 264       adcPadcN6all->SetFillColor(42);
 265       adcPadcN6all->Draw();
 266 c2->cd(3);
 267 gPad->SetFillColor(33);
 268       adcPadcN5cut->SetFillColor(42);
 269       adcPadcN5cut->Draw();
 270 c2->cd(4);
 271 gPad->SetFillColor(33);
 272       adcPadcN6cut->SetFillColor(46);
 273       adcPadcN6cut->Draw();
 274
 275
 276       /*
 277 c2->cd(1);
 278 gPad->SetFillColor(33);
 279       ntuple1->SetFillColor(42);
 280       ntuple1->Draw("qclP","noverprim>=0&&qclP<500");
 281 c2->cd(2);
 282 gPad->SetFillColor(33);
 283       ntuple1->SetFillColor(42);
 284       ntuple1->Draw("qclN","noverprim>=0&&qclN<500");
 285 c2->cd(3);
 286 gPad->SetFillColor(33);
 287       adcPadcN5cut->SetFillColor(42);
 288       adcPadcN5cut->Draw();
 289 c2->cd(4);
 290 gPad->SetFillColor(33);
 291       adcPadcN6cut->SetFillColor(46);
 292       adcPadcN6cut->Draw();
 293       */
 294
 295
 296       /*
 297 c2->Draw();
 298 c2->Print("ssd_res.ps");
 299       */
 300
 301       /*
 302 c2->Draw();
 303 c2->Print("spd_clsize.ps");
 304       */
 305
 306 }
 307
 308
 309
 310
 311
 312
 313
 314
 315