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
  14 /////////////////////////  Ntuple analysis ///////////////////////////////
  15
  16 // ntuple is created inside the hit loop for the hits in the cluster region;
  17
  18 // ntuple1 is created after a finish of the hit loop if one or more hits
  19 // are in the cluster region;
  20
  21 // ntuple2 is created befor the hit loop for all clusters;
  22
  23 // -----------------------------------------------------------------------
  24 // lay       - number of ITS layer;
  25 // x        - coordinates in r*phi(x) direction (mm);
  26 // z        - coordinates in z direction (mm);
  27 // nx        - cluster size in the r*phi(x) direction;
  28 // nz        - cluster size in the z direction;
  29 // hitprim   - primary particle(hit) flag ( = 1 for primery particle);
  30 // dx        - difference of hit(mediate) and reconstructed (from cluster)
  31 //             coordinates in r*phi(x) direction;
  32 // dz        - difference of hit(mediate) and reconstructed (from cluster)
  33 //             coordinates in z direction;
  34 // noverlaps - number of particles overlapping in one cluster found in this
  35 //             macros in the cluster region;
  36 // ntrover   - number of particles overlapping in one cluster found in the
  37 //             AliITSClusterFinderSPD class;
  38 // noverprim - number of primary particles overlapping in one cluster;
  39 // qcl       - cluster charge in electron number
  40 // -------------------------------------------------------------------------
  41
  42 /*
  43 c2->cd(1);
  44 gPad->SetFillColor(33);
  45       ntuple->SetFillColor(42);
  46       ntuple->Draw("dx","lay == 1");
  47 c2->cd(2);
  48 gPad->SetFillColor(33);
  49       ntuple->SetFillColor(46);
  50       ntuple->Draw("dz","lay == 1");
  51 c2->cd(3);
  52 gPad->SetFillColor(33);
  53       ntuple->SetFillColor(42);
  54       ntuple->Draw("dx","lay == 2");
  55 c2->cd(4);
  56 gPad->SetFillColor(33);
  57       ntuple->SetFillColor(46);
  58       ntuple->Draw("dz","lay == 2");
  59 */
  60
  61 /*
  62 c2->cd(1);
  63 gPad->SetFillColor(33);
  64       ntuple->SetFillColor(42);
  65       ntuple->Draw("dx","lay == 1 && hitprim == 1");
  66 c2->cd(2);
  67 gPad->SetFillColor(33);
  68       ntuple->SetFillColor(46);
  69       ntuple->Draw("dz","lay == 1 && hitprim == 1");
  70 c2->cd(3);
  71 gPad->SetFillColor(33);
  72       ntuple->SetFillColor(42);
  73       ntuple->Draw("dx","lay == 2 && hitprim == 1");
  74 c2->cd(4);
  75 gPad->SetFillColor(33);
  76       ntuple->SetFillColor(46);
  77       ntuple->Draw("dz","lay == 2 && hitprim == 1");
  78 */
  79
  80 /*
  81 c2->cd(1);
  82 gPad->SetFillColor(33);
  83       ntuple->SetFillColor(42);
  84       ntuple->Draw("dx","lay == 1 && hitprim == 1 && ntrover>=1&&nx==4");
  85 c2->cd(2);
  86 gPad->SetFillColor(33);
  87       ntuple->SetFillColor(46);
  88       ntuple->Draw("dz","lay == 1 && hitprim == 1 && ntrover>=1&&nz==4");
  89 c2->cd(3);
  90 gPad->SetFillColor(33);
  91       ntuple->SetFillColor(42);
  92       ntuple->Draw("dx","lay == 2 && hitprim == 1 && ntrover>=1&&nx==4");
  93 c2->cd(4);
  94 gPad->SetFillColor(33);
  95       ntuple->SetFillColor(46);
  96       ntuple->Draw("dz","lay == 2 && hitprim == 1 && ntrover>=1&&nz==4");
  97 */
  98
  99
 100 c2->cd(1);
 101 gPad->SetFillColor(33);
 102       ntuple->SetFillColor(42);
 103       ntuple->Draw("dx","lay==1&&hitprim==1&&dx>-150&&dx<150");
 104 c2->cd(2);
 105 gPad->SetFillColor(33);
 106       ntuple->SetFillColor(46);
 107       ntuple->Draw("dz","lay==1&&hitprim==1&&dz>-500&&dz<500");
 108 c2->cd(3);
 109 gPad->SetFillColor(33);
 110       ntuple->SetFillColor(42);
 111       ntuple->Draw("dx","lay==2&&hitprim==1&&dx>-150&&dx<150");
 112 c2->cd(4);
 113 gPad->SetFillColor(33);
 114       ntuple->SetFillColor(46);
 115       ntuple->Draw("dz","lay==2&&hitprim==1&&dz>-500&&dz<500");
 116
 117
 118 /*
 119 c2->cd(1);
 120 gPad->SetFillColor(33);
 121       ntuple->SetFillColor(42);
 122       ntuple->Draw("dx","lay==1&&hitprim==1&&dx>-150&&dx<150&&ntrover==1");
 123 c2->cd(2);
 124 gPad->SetFillColor(33);
 125       ntuple->SetFillColor(46);
 126       ntuple->Draw("dz","lay==1&&hitprim==1&&dz>-500&&dz<500&&ntrover==1");
 127 c2->cd(3);
 128 gPad->SetFillColor(33);
 129       ntuple->SetFillColor(42);
 130       ntuple->Draw("dx","lay==2&&hitprim==1&&dx>-150&&dx<150&&ntrover==1");
 131 c2->cd(4);
 132 gPad->SetFillColor(33);
 133       ntuple->SetFillColor(46);
 134       ntuple->Draw("dz","lay==2&&hitprim==1&&dz>-500&&dz<500&&ntrover==1");
 135 */
 136
 137       /*
 138 c2->cd(1);
 139 gPad->SetFillColor(33);
 140       ntuple->SetFillColor(42);
 141       ntuple->Draw("dx","lay==1&&hitprim==1&&nx>4");
 142 c2->cd(2);
 143 gPad->SetFillColor(33);
 144       ntuple->SetFillColor(46);
 145       ntuple->Draw("dz","lay==1&&hitprim==1&&nz>4");
 146 c2->cd(3);
 147 gPad->SetFillColor(33);
 148       ntuple->SetFillColor(42);
 149       ntuple->Draw("dx","lay==2&&hitprim==1&&nx>4");
 150 c2->cd(4);
 151 gPad->SetFillColor(33);
 152       ntuple->SetFillColor(46);
 153       ntuple->Draw("dz","lay==2&&hitprim==1&&nz>4");
 154       */
 155
 156 /*
 157 c2->cd(1);
 158 gPad->SetFillColor(33);
 159       ntuple1->SetFillColor(42);
 160       ntuple1->Draw("noverlaps","lay == 1&&noverprim>0&&noverlaps>0");
 161 c2->cd(2);
 162 gPad->SetFillColor(33);
 163       ntuple1->SetFillColor(46);
 164       ntuple1->Draw("noverlaps","lay == 2&&noverprim>0&&noverlaps>0");
 165 c2->cd(3);
 166 gPad->SetFillColor(33);
 167       ntuple1->SetFillColor(42);
 168       ntuple1->Draw("ntrover","lay == 1&&noverprim>0&&ntrover>0");
 169 c2->cd(4);
 170 gPad->SetFillColor(33);
 171       ntuple1->SetFillColor(46);
 172       ntuple1->Draw("ntrover","lay == 2&&noverprim>0&&ntrover>0");
 173 */
 174
 175       /*
 176 c2->cd(1);
 177 gPad->SetFillColor(33);
 178       ntuple1->SetFillColor(42);
 179       ntuple1->Draw("noverprim","lay == 1");
 180 c2->cd(2);
 181 gPad->SetFillColor(33);
 182       ntuple1->SetFillColor(46);
 183       ntuple1->Draw("noverprim","lay == 2");
 184 c2->cd(3);
 185 gPad->SetFillColor(33);
 186       ntuple1->SetFillColor(42);
 187       ntuple1->Draw("noverprim","lay == 1 && noverprim < 1");
 188 c2->cd(4);
 189 gPad->SetFillColor(33);
 190       ntuple1->SetFillColor(46);
 191       ntuple1->Draw("noverprim","lay == 2 && noverprim < 1");
 192       */
 193
 194 /*
 195 c2->cd(1);
 196 gPad->SetFillColor(33);
 197       ntuple1->SetFillColor(42);
 198       ntuple1->Draw("nx:nz","lay == 1");
 199 c2->cd(2);
 200 gPad->SetFillColor(33);
 201       ntuple1->SetFillColor(46);
 202       ntuple1->Draw("nx:nz","lay == 2");
 203 */
 204
 205       /*
 206 c2->cd(1);
 207 gPad->SetFillColor(33);
 208       ntuple1->SetFillColor(42);
 209       ntuple1->Draw("qcl","lay == 1");
 210 c2->cd(2);
 211 gPad->SetFillColor(33);
 212       ntuple1->SetFillColor(46);
 213       ntuple1->Draw("qcl","lay == 2");
 214       */
 215
 216 /*
 217 c2->cd(1);
 218 gPad->SetFillColor(33);
 219       ntuple1->SetFillColor(42);
 220       ntuple1->Draw("nx","lay == 1");
 221 c2->cd(2);
 222 gPad->SetFillColor(33);
 223       ntuple1->SetFillColor(46);
 224       ntuple1->Draw("nx","lay == 2");
 225 c2->cd(3);
 226 gPad->SetFillColor(33);
 227       ntuple1->SetFillColor(42);
 228       ntuple1->Draw("nz","lay == 1");
 229 c2->cd(4);
 230 gPad->SetFillColor(33);
 231       ntuple1->SetFillColor(46);
 232       ntuple1->Draw("nz","lay == 2");
 233 */
 234
 235       /*
 236 c2->cd(1);
 237 gPad->SetFillColor(33);
 238       ntuple1->SetFillColor(42);
 239       ntuple1->Draw("nx","lay == 1&&noverprim>0");
 240 c2->cd(2);
 241 gPad->SetFillColor(33);
 242       ntuple1->SetFillColor(46);
 243       ntuple1->Draw("nx","lay == 2&&noverprim>0");
 244 c2->cd(3);
 245 gPad->SetFillColor(33);
 246       ntuple1->SetFillColor(42);
 247       ntuple1->Draw("nz","lay == 1 && noverprim>0");
 248 c2->cd(4);
 249 gPad->SetFillColor(33);
 250       ntuple1->SetFillColor(46);
 251       ntuple1->Draw("nz","lay == 2 && noverprim>0");
 252       */
 253
 254 /*
 255 c2->cd(1);
 256 gPad->SetFillColor(33);
 257       ntuple1->SetFillColor(42);
 258       ntuple1->Draw("x","lay == 1&&noverprim>=0");
 259 c2->cd(2);
 260 gPad->SetFillColor(33);
 261       ntuple1->SetFillColor(46);
 262       ntuple1->Draw("x","lay == 2&&noverprim>=0");
 263 c2->cd(3);
 264 gPad->SetFillColor(33);
 265       ntuple1->SetFillColor(42);
 266       ntuple1->Draw("z","lay == 1 && noverprim>=0");
 267 c2->cd(4);
 268 gPad->SetFillColor(33);
 269       ntuple1->SetFillColor(46);
 270       ntuple1->Draw("z","lay == 2 && noverprim>=0");
 271 */
 272
 273       /*
 274 c2->cd(1);
 275 gPad->SetFillColor(33);
 276       ntuple1->SetFillColor(42);
 277       ntuple1->Draw("nx","lay == 1&&noverprim>0 && ntrover == 1");
 278 c2->cd(2);
 279 gPad->SetFillColor(33);
 280       ntuple1->SetFillColor(46);
 281       ntuple1->Draw("nx","lay == 2&&noverprim>0 && ntrover == 1");
 282 c2->cd(3);
 283 gPad->SetFillColor(33);
 284       ntuple1->SetFillColor(42);
 285       ntuple1->Draw("nz","lay == 1 && noverprim>0 && ntrover == 1");
 286 c2->cd(4);
 287 gPad->SetFillColor(33);
 288       ntuple1->SetFillColor(46);
 289       ntuple1->Draw("nz","lay == 2 && noverprim>0 && ntrover ==1");
 290       */
 291
 292 /*
 293 c2->cd(1);
 294 gPad->SetFillColor(33);
 295       ntuple1->SetFillColor(42);
 296       ntuple1->Draw("dx","lay == 1&&noverprim>0");
 297 c2->cd(2);
 298 gPad->SetFillColor(33);
 299       ntuple1->SetFillColor(46);
 300       ntuple1->Draw("dx","lay == 2&&noverprim>0");
 301 c2->cd(3);
 302 gPad->SetFillColor(33);
 303       ntuple1->SetFillColor(42);
 304       ntuple1->Draw("dz","lay == 1 && noverprim>0");
 305 c2->cd(4);
 306 gPad->SetFillColor(33);
 307       ntuple1->SetFillColor(46);
 308       ntuple1->Draw("dz","lay == 2 && noverprim>0");
 309 */
 310
 311       /*
 312 c2->cd(1);
 313 gPad->SetFillColor(33);
 314       ntuple1->SetFillColor(42);
 315       ntuple1->Draw("dx","lay == 1");
 316 c2->cd(2);
 317 gPad->SetFillColor(33);
 318       ntuple1->SetFillColor(46);
 319       ntuple1->Draw("dx","lay == 2");
 320 c2->cd(3);
 321 gPad->SetFillColor(33);
 322       ntuple1->SetFillColor(42);
 323       ntuple1->Draw("dz","lay == 1");
 324 c2->cd(4);
 325 gPad->SetFillColor(33);
 326       ntuple1->SetFillColor(46);
 327       ntuple1->Draw("dz","lay == 2");
 328       */
 329
 330       /*
 331 c2->cd(1);
 332 gPad->SetFillColor(33);
 333       ntuple2->SetFillColor(42);
 334       ntuple2->Draw("x","lay==1");
 335 c2->cd(2);
 336 gPad->SetFillColor(33);
 337       ntuple2->SetFillColor(46);
 338       ntuple2->Draw("z","lay==1");
 339 c2->cd(3);
 340 gPad->SetFillColor(33);
 341       ntuple2->SetFillColor(42);
 342       ntuple2->Draw("x","lay==2");
 343 c2->cd(4);
 344 gPad->SetFillColor(33);
 345       ntuple2->SetFillColor(46);
 346       ntuple2->Draw("z","lay==2");
 347       */
 348
 349       /*
 350 c2->cd(1);
 351 gPad->SetFillColor(33);
 352       ntuple2->SetFillColor(42);
 353       ntuple2->Draw("nx","lay==1");
 354 c2->cd(2);
 355 gPad->SetFillColor(33);
 356       ntuple2->SetFillColor(46);
 357       ntuple2->Draw("nz","lay==1");
 358 c2->cd(3);
 359 gPad->SetFillColor(33);
 360       ntuple2->SetFillColor(42);
 361       ntuple2->Draw("nx","lay==2");
 362 c2->cd(4);
 363 gPad->SetFillColor(33);
 364       ntuple2->SetFillColor(46);
 365       ntuple2->Draw("nz","lay==2");
 366       */
 367
 368 /*
 369 c2->cd(1);
 370 gPad->SetFillColor(33);
 371       ntuple2->SetFillColor(42);
 372       ntuple2->Draw("qcl","lay==1");
 373 c2->cd(2);
 374 gPad->SetFillColor(33);
 375       ntuple2->SetFillColor(46);
 376       ntuple2->Draw("nx","lay==1");
 377 c2->cd(3);
 378 gPad->SetFillColor(33);
 379       ntuple2->SetFillColor(42);
 380       ntuple2->Draw("qcl","lay==2");
 381 c2->cd(4);
 382 gPad->SetFillColor(33);
 383       ntuple2->SetFillColor(46);
 384       ntuple2->Draw("nz","lay==2");
 385 */
 386
 387 /////////////////////   Histogramm analysis  ////////////////////////
 388
 389 /*
 390 c2->cd(1);
 391 gPad->SetFillColor(33);
 392       Xres1->SetFillColor(42);
 393       Xres1->Draw();
 394 c2->cd(2);
 395 gPad->SetFillColor(33);
 396       Zres1->SetFillColor(42);
 397       Zres1->Draw();
 398 c2->cd(3);
 399 gPad->SetFillColor(33);
 400       Xres2->SetFillColor(46);
 401       Xres2->Draw();
 402 c2->cd(4);
 403 gPad->SetFillColor(33);
 404       Zres2->SetFillColor(46);
 405       Zres2->Draw();
 406 */
 407
 408       /*
 409 c2->cd(1);
 410 gPad->SetFillColor(33);
 411       Nzpix1->SetFillColor(42);
 412       Nzpix1->Draw();
 413 c2->cd(2);
 414 gPad->SetFillColor(33);
 415       Nxpix1->SetFillColor(42);
 416       Nxpix1->Draw();
 417 c2->cd(3);
 418 gPad->SetFillColor(33);
 419       Nzpix2->SetFillColor(46);
 420       Nzpix2->Draw();
 421 c2->cd(4);
 422 gPad->SetFillColor(33);
 423       Nxpix2->SetFillColor(46);
 424       Nxpix2->Draw();
 425       */
 426
 427       /*
 428 c2->cd(1);
 429 gPad->SetFillColor(33);
 430       Zpix1->SetFillColor(42);
 431       Zpix1->Draw();
 432 c2->cd(2);
 433 gPad->SetFillColor(33);
 434       Xpix1->SetFillColor(46);
 435       Xpix1->Draw();
 436 c2->cd(3);
 437 gPad->SetFillColor(33);
 438       Zpix2->SetFillColor(42);
 439       Zpix2->Draw();
 440 c2->cd(4);
 441 gPad->SetFillColor(33);
 442       Xpix2->SetFillColor(46);
 443       Xpix2->Draw();
 444       */
 445
 446
 447 /*
 448 c2->cd(1);
 449 gPad->SetFillColor(33);
 450       Theta1->SetFillColor(42);
 451       Theta1->Draw();
 452 c2->cd(2);
 453 gPad->SetFillColor(33);
 454       Phi1->SetFillColor(46);
 455       Phi1->Draw();
 456 c2->cd(3);
 457 gPad->SetFillColor(33);
 458       Theta2->SetFillColor(42);
 459       Theta2->Draw();
 460 c2->cd(4);
 461 gPad->SetFillColor(33);
 462       Phi2->SetFillColor(46);
 463       Phi2->Draw();
 464 */
 465
 466 /*
 467 c2->cd(1);
 468 gPad->SetFillColor(33);
 469       Ptot1->SetFillColor(42);
 470       Ptot1->Draw();
 471 c2->cd(2);
 472 gPad->SetFillColor(33);
 473       Pz1->SetFillColor(46);
 474       Pz1->Draw();
 475 c2->cd(3);
 476 gPad->SetFillColor(33);
 477       Ptot2->SetFillColor(42);
 478       Ptot2->Draw();
 479 c2->cd(4);
 480 gPad->SetFillColor(33);
 481       Pz2->SetFillColor(46);
 482       Pz2->Draw();
 483 */
 484
 485 /*
 486 c2->cd(1);
 487 gPad->SetFillColor(33);
 488       Eta1->SetFillColor(42);
 489       Eta1->Draw();
 490 c2->cd(2);
 491 gPad->SetFillColor(33);
 492       Y1->SetFillColor(46);
 493       Y1->Draw();
 494 c2->cd(3);
 495 gPad->SetFillColor(33);
 496       Eta2->SetFillColor(42);
 497       Eta2->Draw();
 498 c2->cd(4);
 499 gPad->SetFillColor(33);
 500       Y2->SetFillColor(46);
 501       Y2->Draw();
 502 */
 503
 504       /*
 505 c2->cd(1);
 506 gPad->SetFillColor(33);
 507       Eta1Den->SetFillColor(42);
 508       Eta1Den->Draw();
 509 c2->cd(2);
 510 gPad->SetFillColor(33);
 511       Y1Den->SetFillColor(46);
 512       Y1Den->Draw();
 513 c2->cd(3);
 514 gPad->SetFillColor(33);
 515       Eta2Den->SetFillColor(42);
 516       Eta2Den->Draw();
 517 c2->cd(4);
 518 gPad->SetFillColor(33);
 519       Y2Den->SetFillColor(46);
 520       Y2Den->Draw();
 521       */
 522
 523       /*
 524 c2->cd(1);
 525 gPad->SetFillColor(33);
 526       Eta1DenA->SetFillColor(42);
 527       Eta1DenA->Draw();
 528 c2->cd(2);
 529 gPad->SetFillColor(33);
 530       Y1DenA->SetFillColor(46);
 531       Y1DenA->Draw();
 532 c2->cd(3);
 533 gPad->SetFillColor(33);
 534       Eta2DenA->SetFillColor(42);
 535       Eta2DenA->Draw();
 536 c2->cd(4);
 537 gPad->SetFillColor(33);
 538       Y2DenA->SetFillColor(46);
 539       Y2DenA->Draw();
 540       */
 541
 542       /*
 543 c2->Draw();
 544 c2->Print("spd_res.ps");
 545       */
 546       /*
 547 c2->Draw();
 548 c2->Print("spd_clsize.ps");
 549       */
 550
 551 }
 552
 553
 554
 555
 556
 557
 558
 559