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","SPD 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&&abs(dx)<200");
  49 c2->cd(2);
  50 gPad->SetFillColor(33);
  51       ntuple->SetFillColor(46);
  52       ntuple->Draw("dz","lay == 5 && hitprim == 1&&abs(dz)<5000");
  53 c2->cd(3);
  54 gPad->SetFillColor(33);
  55       ntuple->SetFillColor(42);
  56       ntuple->Draw("dx","lay == 6 && hitprim == 1&&abs(dx)<200");
  57 c2->cd(4);
  58 gPad->SetFillColor(33);
  59       ntuple->SetFillColor(46);
  60       ntuple->Draw("dz","lay == 6 && hitprim == 1&&abs(dz)<5000");
  61
  62
  63 /*
  64 c2->cd(1);
  65 gPad->SetFillColor(33);
  66       ntuple1->SetFillColor(42);
  67       ntuple1->Draw("nxP","lay == 5&&noverprim>=0");
  68 c2->cd(2);
  69 gPad->SetFillColor(33);
  70       ntuple1->SetFillColor(46);
  71       ntuple1->Draw("nxN","lay == 5&&noverprim>=0");
  72 c2->cd(3);
  73 gPad->SetFillColor(33);
  74       ntuple1->SetFillColor(42);
  75       ntuple1->Draw("nxP","lay == 6 && noverprim>=0");
  76 c2->cd(4);
  77 gPad->SetFillColor(33);
  78       ntuple1->SetFillColor(46);
  79       ntuple1->Draw("nxN","lay == 6 && noverprim>=0");
  80 */
  81
  82 /*
  83 c2->cd(1);
  84 gPad->SetFillColor(33);
  85       ntuple1->SetFillColor(42);
  86       ntuple1->Draw("qclP","lay == 5&&noverprim>=0");
  87 c2->cd(2);
  88 gPad->SetFillColor(33);
  89       ntuple1->SetFillColor(42);
  90       ntuple1->Draw("qclN","lay == 5&&noverprim>=0");
  91 c2->cd(3);
  92 gPad->SetFillColor(33);
  93       ntuple1->SetFillColor(46);
  94       ntuple1->Draw("qclP","lay == 6&&noverprim>=0");
  95 c2->cd(4);
  96 gPad->SetFillColor(33);
  97       ntuple1->SetFillColor(46);
  98       ntuple1->Draw("qclN","lay == 6&&noverprim>=0");
  99 */
 100
 101
 102 /////////////////////   Histogramm/ntuple  analysis  ////////////////////////
 103
 104       /*
 105 c2->cd(1);
 106 gPad->SetFillColor(33);
 107       ntuple1->SetFillColor(42);
 108       ntuple1->Draw("qclP","noverprim>=0");
 109 c2->cd(2);
 110 gPad->SetFillColor(33);
 111       ntuple1->SetFillColor(42);
 112       ntuple1->Draw("qclN","noverprim>=0");
 113 c2->cd(3);
 114 gPad->SetFillColor(33);
 115       adcPadcN5cut->SetFillColor(42);
 116       adcPadcN5cut->Draw();
 117 c2->cd(4);
 118 gPad->SetFillColor(33);
 119       adcPadcN6cut->SetFillColor(46);
 120       adcPadcN6cut->Draw();
 121       */
 122
 123
 124       /*
 125 c2->Draw();
 126 c2->Print("ssd_res.ps");
 127       */
 128
 129       /*
 130 c2->Draw();
 131 c2->Print("spd_clsize.ps");
 132       */
 133
 134 }
 135
 136
 137
 138
 139
 140
 141
 142
 143