ITS/SPD_ntuple.C

   1 void SPD_ntuple()
   2 {
   3
   4 TFile *f = new TFile("SPD_his.root");
   5 //TH1F *Pix = (TH1F*)f->Get("Pix");
   6
   7 //gStyle->SetOptStat(1111111);
   8 gStyle->SetOptLogy();
   9 //TCanvas *c1 = new TCanvas("c1","SPD clusters",400,10,600,700);
  10 TCanvas *c2 = new TCanvas("c2","SPD clusters",400,10,600,700);
  11 //c1->Divide(2,2);
  12 c2->Divide(2,2);
  13 //c2->SetLogy();
  14
  15 /////////////////////////  Ntuple analysis ///////////////////////////////
  16
  17 // ntuple is created inside the hit loop for the hits in the cluster region;
  18
  19 // ntuple1 is created after a finish of the hit loop if one or more hits
  20 // are in the cluster region;
  21
  22 // ntuple2 is created befor the hit loop for all clusters;
  23
  24 // -----------------------------------------------------------------------
  25 // lay       - number of ITS layer;
  26 // nx        - cluster size in the r*phi(x) direction;
  27 // nz        - cluster size in the z direction;
  28 // hitprim   - primary particle(hit) flag ( = 1 for primery particle);
  29 // dx        - difference of hit(mediate) and reconstructed (from cluster)
  30 //             coordinates in r*phi(x) direction;
  31 // dz        - difference of hit(mediate) and reconstructed (from cluster)
  32 //             coordinates in z direction;
  33 // noverlaps - number of particles overlapping in one cluster;
  34 // noverprim - number of primary particles overlapping in one cluster;
  35 // qcl       - cluster charge in electron number
  36 // -------------------------------------------------------------------------
  37
  38 /*
  39 c2->cd(1);
  40 gPad->SetFillColor(33);
  41       ntuple->SetFillColor(42);
  42       ntuple->Draw("dx","lay == 1 && hitprim == 1");
  43 c2->cd(2);
  44 gPad->SetFillColor(33);
  45       ntuple->SetFillColor(46);
  46       ntuple->Draw("dz","lay == 1 && hitprim == 1");
  47 c2->cd(3);
  48 gPad->SetFillColor(33);
  49       ntuple->SetFillColor(42);
  50       ntuple->Draw("dx","lay == 2 && hitprim == 1");
  51 c2->cd(4);
  52 gPad->SetFillColor(33);
  53       ntuple->SetFillColor(46);
  54       ntuple->Draw("dz","lay == 2 && hitprim == 1");
  55 */
  56
  57 /*
  58 c2->cd(1);
  59 gPad->SetFillColor(33);
  60       ntuple->SetFillColor(42);
  61       ntuple->Draw("dx","lay == 1 && hitprim == 1");
  62 c2->cd(2);
  63 gPad->SetFillColor(33);
  64       ntuple->SetFillColor(46);
  65       ntuple->Draw("dx","lay == 1 && hitprim == 1 && nx < 15");
  66 c2->cd(3);
  67 gPad->SetFillColor(33);
  68       ntuple->SetFillColor(42);
  69       ntuple->Draw("dx","lay == 1 && hitprim == 1 && nx < 5");
  70 c2->cd(4);
  71 gPad->SetFillColor(33);
  72       ntuple->SetFillColor(46);
  73       ntuple->Draw("dx","lay == 1 && hitprim == 1 && nx < 4");
  74 */
  75
  76       /*
  77 c2->cd(1);
  78 gPad->SetFillColor(33);
  79       ntuple->SetFillColor(42);
  80       ntuple->Draw("dx","lay==1&&hitprim==1&&dx>-150&&dx<150");
  81 c2->cd(2);
  82 gPad->SetFillColor(33);
  83       ntuple->SetFillColor(46);
  84       ntuple->Draw("dz","lay==1&&hitprim==1&&dz>-500&&dz<500");
  85 c2->cd(3);
  86 gPad->SetFillColor(33);
  87       ntuple->SetFillColor(42);
  88       ntuple->Draw("dx","lay==2&&hitprim==1&&dx>-150&&dx<150");
  89 c2->cd(4);
  90 gPad->SetFillColor(33);
  91       ntuple->SetFillColor(46);
  92       ntuple->Draw("dz","lay==2&&hitprim==1&&dz>-500&&dz<500");
  93       */
  94 /*
  95 c2->cd(1);
  96 gPad->SetFillColor(33);
  97       ntuple->SetFillColor(42);
  98       ntuple->Draw("dx","lay==1&&hitprim==1&&dx>-150&&dx<150&&nx<5");
  99 c2->cd(2);
 100 gPad->SetFillColor(33);
 101       ntuple->SetFillColor(46);
 102       ntuple->Draw("dz","lay==1&&hitprim==1&&dz>-500&&dz<500&&nz<5");
 103 c2->cd(3);
 104 gPad->SetFillColor(33);
 105       ntuple->SetFillColor(42);
 106       ntuple->Draw("dx","lay==2&&hitprim==1&&dx>-150&&dx<150&&nx<5");
 107 c2->cd(4);
 108 gPad->SetFillColor(33);
 109       ntuple->SetFillColor(46);
 110       ntuple->Draw("dz","lay==2&&hitprim==1&&dz>-500&&dz<500&&nz<5");
 111 */
 112 /*
 113 c2->cd(1);
 114 gPad->SetFillColor(33);
 115       ntuple->SetFillColor(42);
 116       ntuple->Draw("dx","lay==1&&hitprim==1&&nx<4");
 117 c2->cd(2);
 118 gPad->SetFillColor(33);
 119       ntuple->SetFillColor(46);
 120       ntuple->Draw("dz","lay==1&&hitprim==1&&nz<4");
 121 c2->cd(3);
 122 gPad->SetFillColor(33);
 123       ntuple->SetFillColor(42);
 124       ntuple->Draw("dx","lay==2&&hitprim==1&&nx<4");
 125 c2->cd(4);
 126 gPad->SetFillColor(33);
 127       ntuple->SetFillColor(46);
 128       ntuple->Draw("dz","lay==2&&hitprim==1&&nz<4");
 129 */
 130
 131 /*
 132 c2->cd(1);
 133 gPad->SetFillColor(33);
 134       ntuple1->SetFillColor(42);
 135       ntuple1->Draw("noverlaps","lay == 1");
 136 c2->cd(2);
 137 gPad->SetFillColor(33);
 138       ntuple1->SetFillColor(46);
 139       ntuple1->Draw("noverlaps","lay == 2");
 140
 141 c2->cd(3);
 142 gPad->SetFillColor(33);
 143       ntuple1->SetFillColor(42);
 144       ntuple1->Draw("noverlaps","lay == 1 && noverprim > 0");
 145 c2->cd(4);
 146 gPad->SetFillColor(33);
 147       ntuple1->SetFillColor(46);
 148       ntuple1->Draw("noverlaps","lay == 2 && noverprim > 0");
 149 */
 150 /*
 151 c2->cd(1);
 152 gPad->SetFillColor(33);
 153       ntuple1->SetFillColor(42);
 154       ntuple1->Draw("noverprim","lay == 1");
 155 c2->cd(2);
 156 gPad->SetFillColor(33);
 157       ntuple1->SetFillColor(46);
 158       ntuple1->Draw("noverprim","lay == 2");
 159 c2->cd(3);
 160 gPad->SetFillColor(33);
 161       ntuple1->SetFillColor(42);
 162       ntuple1->Draw("noverprim","lay == 1 && noverprim > 1");
 163 c2->cd(4);
 164 gPad->SetFillColor(33);
 165       ntuple1->SetFillColor(46);
 166       ntuple1->Draw("noverprim","lay == 2 && noverprim > 1");
 167 */
 168 /*
 169 c2->cd(1);
 170 gPad->SetFillColor(33);
 171       ntuple1->SetFillColor(42);
 172       ntuple1->Draw("qcl","lay == 1&&noverprim>0");
 173 c2->cd(2);
 174 gPad->SetFillColor(33);
 175       ntuple1->SetFillColor(46);
 176       ntuple1->Draw("qcl","lay == 2&&noverprim>0");
 177 */
 178
 179
 180 c2->cd(1);
 181 gPad->SetFillColor(33);
 182       ntuple1->SetFillColor(42);
 183       ntuple1->Draw("nx","lay == 1&&noverprim>0");
 184 c2->cd(2);
 185 gPad->SetFillColor(33);
 186       ntuple1->SetFillColor(46);
 187       ntuple1->Draw("nx","lay == 2&&noverprim>0");
 188 c2->cd(3);
 189 gPad->SetFillColor(33);
 190       ntuple1->SetFillColor(42);
 191       ntuple1->Draw("nz","lay == 1 && noverprim>0");
 192 c2->cd(4);
 193 gPad->SetFillColor(33);
 194       ntuple1->SetFillColor(46);
 195       ntuple1->Draw("nz","lay == 2 && noverprim>0");
 196
 197
 198 /*
 199 c2->cd(1);
 200 gPad->SetFillColor(33);
 201       ntuple1->SetFillColor(42);
 202       ntuple1->Draw("dx","lay == 1&&noverprim>0");
 203 c2->cd(2);
 204 gPad->SetFillColor(33);
 205       ntuple1->SetFillColor(46);
 206       ntuple1->Draw("dx","lay == 2&&noverprim>0");
 207 c2->cd(3);
 208 gPad->SetFillColor(33);
 209       ntuple1->SetFillColor(42);
 210       ntuple1->Draw("dz","lay == 1 && noverprim>0");
 211 c2->cd(4);
 212 gPad->SetFillColor(33);
 213       ntuple1->SetFillColor(46);
 214       ntuple1->Draw("dz","lay == 2 && noverprim>0");
 215 */
 216
 217       /*
 218 c2->cd(1);
 219 gPad->SetFillColor(33);
 220       ntuple1->SetFillColor(42);
 221       ntuple1->Draw("dx","lay == 1");
 222 c2->cd(2);
 223 gPad->SetFillColor(33);
 224       ntuple1->SetFillColor(46);
 225       ntuple1->Draw("dx","lay == 2");
 226 c2->cd(3);
 227 gPad->SetFillColor(33);
 228       ntuple1->SetFillColor(42);
 229       ntuple1->Draw("dz","lay == 1");
 230 c2->cd(4);
 231 gPad->SetFillColor(33);
 232       ntuple1->SetFillColor(46);
 233       ntuple1->Draw("dz","lay == 2");
 234       */
 235       /*
 236 c2->cd(1);
 237 gPad->SetFillColor(33);
 238       ntuple2->SetFillColor(42);
 239       ntuple2->Draw("nx","lay == 1");
 240 c2->cd(2);
 241 gPad->SetFillColor(33);
 242       ntuple2->SetFillColor(46);
 243       ntuple2->Draw("nz","lay == 1");
 244 c2->cd(3);
 245 gPad->SetFillColor(33);
 246       ntuple2->SetFillColor(42);
 247       ntuple2->Draw("nx","lay == 2");
 248 c2->cd(4);
 249 gPad->SetFillColor(33);
 250       ntuple2->SetFillColor(46);
 251       ntuple2->Draw("nz","lay == 2");
 252       */
 253
 254 /////////////////////   Histogramm analysis  ////////////////////////
 255
 256       /*
 257 c2->cd(1);
 258 gPad->SetFillColor(33);
 259       Xres1->SetFillColor(42);
 260       Xres1->Draw();
 261 c2->cd(2);
 262 gPad->SetFillColor(33);
 263       Zres1->SetFillColor(42);
 264       Zres1->Draw();
 265 c2->cd(3);
 266 gPad->SetFillColor(33);
 267       Xres2->SetFillColor(46);
 268       Xres2->Draw();
 269 c2->cd(4);
 270 gPad->SetFillColor(33);
 271       Zres2->SetFillColor(46);
 272       Zres2->Draw();
 273       */
 274       /*
 275 c2->cd(1);
 276 gPad->SetFillColor(33);
 277       Nzpix1->SetFillColor(42);
 278       Nzpix1->Draw();
 279 c2->cd(2);
 280 gPad->SetFillColor(33);
 281       Nxpix1->SetFillColor(42);
 282       Nxpix1->Draw();
 283 c2->cd(3);
 284 gPad->SetFillColor(33);
 285       Nzpix2->SetFillColor(46);
 286       Nzpix2->Draw();
 287 c2->cd(4);
 288 gPad->SetFillColor(33);
 289       Nxpix2->SetFillColor(46);
 290       Nxpix2->Draw();
 291       */
 292
 293       /*
 294 c2->cd(1);
 295 gPad->SetFillColor(33);
 296       Zpix1->SetFillColor(42);
 297       Zpix1->Draw();
 298 c2->cd(2);
 299 gPad->SetFillColor(33);
 300       Xpix1->SetFillColor(46);
 301       Xpix1->Draw();
 302 c2->cd(3);
 303 gPad->SetFillColor(33);
 304       Zpix2->SetFillColor(42);
 305       Zpix2->Draw();
 306 c2->cd(4);
 307 gPad->SetFillColor(33);
 308       Xpix2->SetFillColor(46);
 309       Xpix2->Draw();
 310       */
 311
 312       /*
 313 c2->Draw();
 314 c2->Print("spd_res.ps");
 315       */
 316       /*
 317 c2->Draw();
 318 c2->Print("spd_clsize.ps");
 319       */
 320
 321 }
 322
 323
 324
 325
 326
 327
 328
 329