Test macro to analysis the SPD simulation made with the Dubna model

diff --git a/ITS/SPD_ntuple_Dubna.C b/ITS/SPD_ntuple_Dubna.C
new file mode 100644 (file)
index 0000000..4294ab7
--- /dev/null
+++ b/ITS/SPD_ntuple_Dubna.C
@@ -0,0 +1,560 @@
+void SPD_ntuple_Dubna()
+{
+
+TFile *f = new TFile("SPD_his_dubna.root");
+//TH1F *Pix = (TH1F*)f->Get("Pix");
+   
+//gStyle->SetOptStat(1111111);
+//gStyle->SetOptLogy();
+//TCanvas *c1 = new TCanvas("c1","SPD clusters",400,10,600,700);
+TCanvas *c2 = new TCanvas("c2","SPD clusters",400,10,600,700);
+//c1->Divide(2,2);
+c2->Divide(2,2);
+
+/////////////////////////  Ntuple analysis ///////////////////////////////
+
+// ntuple is created inside the hit loop for the hits in the cluster region;
+
+// ntuple1 is created after a finish of the hit loop if one or more hits
+// are in the cluster region;
+
+// ntuple2 is created befor the hit loop for all clusters;
+
+// -----------------------------------------------------------------------
+// lay       - number of ITS layer;
+// x        - coordinates in r*phi(x) direction (mm);
+// z        - coordinates in z direction (mm);
+// nx        - cluster size in the r*phi(x) direction;
+// nz        - cluster size in the z direction;
+// hitprim   - primary particle(hit) flag ( = 1 for primery particle);     
+// dx        - difference of hit(mediate) and reconstructed (from cluster)
+//             coordinates in r*phi(x) direction;
+// dz        - difference of hit(mediate) and reconstructed (from cluster)
+//             coordinates in z direction;
+// noverlaps - number of particles overlapping in one cluster found in this
+//             macros in the cluster region;
+// ntrover   - number of particles overlapping in one cluster found in the
+//             AliITSClusterFinderSPD class; 
+// noverprim - number of primary particles overlapping in one cluster;
+// qcl       - cluster charge in electron number  
+// -------------------------------------------------------------------------
+
+/*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(42);
+      ntuple->Draw("dx","lay == 1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(46);
+      ntuple->Draw("dz","lay == 1");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(42);
+      ntuple->Draw("dx","lay == 2");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(46);
+      ntuple->Draw("dz","lay == 2");
+*/
+
+/*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(42);
+      ntuple->Draw("dx","lay == 1 && hitprim == 1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(46);
+      ntuple->Draw("dz","lay == 1 && hitprim == 1");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(42);
+      ntuple->Draw("dx","lay == 2 && hitprim == 1");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(46);
+      ntuple->Draw("dz","lay == 2 && hitprim == 1");
+*/
+
+/*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(42);
+      ntuple->Draw("dx","lay == 1 && hitprim == 1 && ntrover>=1&&nx==1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(46);
+      ntuple->Draw("dz","lay == 1 && hitprim == 1 && ntrover>=1&&nz==1");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(42);
+      ntuple->Draw("dx","lay == 2 && hitprim == 1 && ntrover>=1&&nx==1");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(46);
+      ntuple->Draw("dz","lay == 2 && hitprim == 1 && ntrover>=1&&nz==1");
+*/
+
+
+/*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(42);
+      ntuple->Draw("dx","lay==1&&hitprim==1&&abs(dx)<150&&nx==1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(46);
+      ntuple->Draw("dz","lay==1&&hitprim==1&&abs(dz)<500&&nz==1");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(42);
+      ntuple->Draw("dx","lay==2&&hitprim==1&&abs(dx)<150&&nx==1");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(46);
+      ntuple->Draw("dz","lay==2&&hitprim==1&&abs(dz)<500&&nz==1");
+*/
+
+
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(42);
+      ntuple->Draw("dx","lay==1&&hitprim==1&&abs(dx)<150&&ntrover==1&&nx>0");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(46);
+      ntuple->Draw("dz","lay==1&&hitprim==1&&abs(dz)<500&&ntrover==1&&nz>0");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(42);
+      ntuple->Draw("dx","lay==2&&hitprim==1&&abs(dx)<150&&ntrover==1&&nx>0");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(46);
+      ntuple->Draw("dz","lay==2&&hitprim==1&&abs(dz)<500&&ntrover==1&&nz>0");
+
+
+      /*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(42);
+      ntuple->Draw("dx","lay==1&&hitprim==1&&nx==1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(46);
+      ntuple->Draw("dz","lay==1&&hitprim==1&&nz==1");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(42);
+      ntuple->Draw("dx","lay==2&&hitprim==1&&nx==1");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple->SetFillColor(46);
+      ntuple->Draw("dz","lay==2&&hitprim==1&&nz==1");
+      */
+
+/*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("noverlaps","lay == 1&&noverprim>0&&noverlaps>1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("noverlaps","lay == 2&&noverprim>0&&noverlaps>1");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("ntrover","lay == 1&&noverprim>0&&ntrover>1");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("ntrover","lay == 2&&noverprim>0&&ntrover>1");
+*/
+
+      /*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("noverprim","lay == 1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("noverprim","lay == 2");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("noverprim","lay == 1 && noverprim < 1");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("noverprim","lay == 2 && noverprim < 1");
+      */
+
+/*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("nx:nz","lay == 1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("nx:nz","lay == 2");
+*/
+
+      /*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("qcl","lay == 1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("qcl","lay == 2");
+      */
+
+      /*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("nx","lay == 1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("nx","lay == 2");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("nz","lay == 1");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("nz","lay == 2");
+      */
+
+      /*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("nx","lay == 1&&noverprim>0");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("nx","lay == 2&&noverprim>0");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("nz","lay == 1 && noverprim>0");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("nz","lay == 2 && noverprim>0");
+      */
+
+/*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("x","lay == 1&&noverprim>=0");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("x","lay == 2&&noverprim>=0");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("z","lay == 1 && noverprim>=0");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("z","lay == 2 && noverprim>=0");
+*/
+
+      /*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("nx","lay == 1&&noverprim>0 && ntrover == 1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("nx","lay == 2&&noverprim>0 && ntrover == 1");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("nz","lay == 1 && noverprim>0 && ntrover == 1");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("nz","lay == 2 && noverprim>0 && ntrover ==1");
+      */
+            
+/*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("dx","lay == 1&&noverprim>0");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("dx","lay == 2&&noverprim>0");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("dz","lay == 1 && noverprim>0");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("dz","lay == 2 && noverprim>0");
+*/
+
+      /*            
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("dx","lay == 1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("dx","lay == 2");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(42);
+      ntuple1->Draw("dz","lay == 1");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple1->SetFillColor(46);
+      ntuple1->Draw("dz","lay == 2");
+      */
+
+      /*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple2->SetFillColor(42);
+      ntuple2->Draw("x","lay==1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple2->SetFillColor(46);
+      ntuple2->Draw("z","lay==1");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple2->SetFillColor(42);
+      ntuple2->Draw("x","lay==2");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple2->SetFillColor(46);
+      ntuple2->Draw("z","lay==2");
+      */
+
+      /*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple2->SetFillColor(42);
+      ntuple2->Draw("nx","lay==1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple2->SetFillColor(46);
+      ntuple2->Draw("nz","lay==1");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple2->SetFillColor(42);
+      ntuple2->Draw("nx","lay==2");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple2->SetFillColor(46);
+      ntuple2->Draw("nz","lay==2");
+      */
+
+/*
+c2->cd(1);
+gPad->SetFillColor(33);
+      ntuple2->SetFillColor(42);
+      ntuple2->Draw("qcl","lay==1");
+c2->cd(2);
+gPad->SetFillColor(33);
+      ntuple2->SetFillColor(46);
+      ntuple2->Draw("nx","lay==1");
+c2->cd(3);
+gPad->SetFillColor(33);
+      ntuple2->SetFillColor(42);
+      ntuple2->Draw("qcl","lay==2");
+c2->cd(4);
+gPad->SetFillColor(33);
+      ntuple2->SetFillColor(46);
+      ntuple2->Draw("nz","lay==2");
+*/
+      
+/////////////////////   Histogramm analysis  ////////////////////////
+
+/*
+c2->cd(1);
+gPad->SetFillColor(33);
+      Xres1->SetFillColor(42);
+      Xres1->Draw();
+c2->cd(2);
+gPad->SetFillColor(33);
+      Zres1->SetFillColor(42);
+      Zres1->Draw();
+c2->cd(3);
+gPad->SetFillColor(33);
+      Xres2->SetFillColor(46);
+      Xres2->Draw();
+c2->cd(4);
+gPad->SetFillColor(33);
+      Zres2->SetFillColor(46);
+      Zres2->Draw();
+*/
+
+      /*
+c2->cd(1);
+gPad->SetFillColor(33);
+      Nzpix1->SetFillColor(42);
+      Nzpix1->Draw();
+c2->cd(2);
+gPad->SetFillColor(33);
+      Nxpix1->SetFillColor(42);
+      Nxpix1->Draw();
+c2->cd(3);
+gPad->SetFillColor(33);
+      Nzpix2->SetFillColor(46);
+      Nzpix2->Draw();
+c2->cd(4);
+gPad->SetFillColor(33);
+      Nxpix2->SetFillColor(46);
+      Nxpix2->Draw();
+      */                         
+      
+      /*                          
+c2->cd(1);
+gPad->SetFillColor(33);
+      Zpix1->SetFillColor(42);
+      Zpix1->Draw();
+c2->cd(2);
+gPad->SetFillColor(33);
+      Xpix1->SetFillColor(46);
+      Xpix1->Draw();
+c2->cd(3);
+gPad->SetFillColor(33);
+      Zpix2->SetFillColor(42);
+      Zpix2->Draw();
+c2->cd(4);
+gPad->SetFillColor(33);
+      Xpix2->SetFillColor(46);
+      Xpix2->Draw();
+      */
+
+
+/*
+c2->cd(1);
+gPad->SetFillColor(33);
+      Theta1->SetFillColor(42);
+      Theta1->Draw();
+c2->cd(2);
+gPad->SetFillColor(33);
+      Phi1->SetFillColor(46);
+      Phi1->Draw();
+c2->cd(3);
+gPad->SetFillColor(33);
+      Theta2->SetFillColor(42);
+      Theta2->Draw();
+c2->cd(4);
+gPad->SetFillColor(33);
+      Phi2->SetFillColor(46);
+      Phi2->Draw();
+*/
+
+/*
+c2->cd(1);
+gPad->SetFillColor(33);
+      Ptot1->SetFillColor(42);
+      Ptot1->Draw();
+c2->cd(2);
+gPad->SetFillColor(33);
+      Pz1->SetFillColor(46);
+      Pz1->Draw();
+c2->cd(3);
+gPad->SetFillColor(33);
+      Ptot2->SetFillColor(42);
+      Ptot2->Draw();
+c2->cd(4);
+gPad->SetFillColor(33);
+      Pz2->SetFillColor(46);
+      Pz2->Draw();
+*/
+
+/*
+c2->cd(1);
+gPad->SetFillColor(33);
+      Eta1->SetFillColor(42);
+      Eta1->Draw();
+c2->cd(2);
+gPad->SetFillColor(33);
+      Y1->SetFillColor(46);
+      Y1->Draw();
+c2->cd(3);
+gPad->SetFillColor(33);
+      Eta2->SetFillColor(42);
+      Eta2->Draw();
+c2->cd(4);
+gPad->SetFillColor(33);
+      Y2->SetFillColor(46);
+      Y2->Draw();
+*/
+
+      /*
+c2->cd(1);
+gPad->SetFillColor(33);
+      Eta1Den->SetFillColor(42);
+      Eta1Den->Draw();
+c2->cd(2);
+gPad->SetFillColor(33);
+      Y1Den->SetFillColor(46);
+      Y1Den->Draw();
+c2->cd(3);
+gPad->SetFillColor(33);
+      Eta2Den->SetFillColor(42);
+      Eta2Den->Draw();
+c2->cd(4);
+gPad->SetFillColor(33);
+      Y2Den->SetFillColor(46);
+      Y2Den->Draw();
+      */
+
+      /*
+c2->cd(1);
+gPad->SetFillColor(33);
+      Eta1DenA->SetFillColor(42);
+      Eta1DenA->Draw();
+c2->cd(2);
+gPad->SetFillColor(33);
+      Y1DenA->SetFillColor(46);
+      Y1DenA->Draw();
+c2->cd(3);
+gPad->SetFillColor(33);
+      Eta2DenA->SetFillColor(42);
+      Eta2DenA->Draw();
+c2->cd(4);
+gPad->SetFillColor(33);
+      Y2DenA->SetFillColor(46);
+      Y2DenA->Draw();
+      */
+     
+      /*                       
+c2->Draw();   
+c2->Print("spd_res.ps");
+      */
+      /*                        
+c2->Draw();   
+c2->Print("spd_clsize.ps");
+      */    
+
+}
+
+
+
+
+
+
+
+