ITS/SPD_ntuple_Dubna.C

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