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