--- /dev/null
+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");
+ */
+
+}
+
+
+
+
+
+
+
+