-
- deltaZ = 4.91617e-01+2.67567e-02*vzero[2]+1.49626e-02*TMath::Power(vzero[2],2);
- Float_t multFactorZ = 28000./(Float_t)nopoints1;
- Int_t nbin=(Int_t)((deltaZ/0.005)/multFactorZ);
- Int_t nbinxy=250;
- Int_t *vectorBinZ,*vectorBinXY;
- vectorBinZ=new Int_t[nbin];
- vectorBinXY=new Int_t[nbinxy];
- Float_t f1= 0;
- Float_t f2= 0;
- Double_t sigma,averagebg;
-
- TH1D *hITSZv = new TH1D("hITSZv","",nbin,vzero[2]-deltaZ,vzero[2]+deltaZ);
- TH1D *hITSXv = new TH1D("hITSXv","",nbinxy,-3,3);
- TH1D *hITSYv = new TH1D("hITSYv","",nbinxy,-3,3);
-
- // cout << "deltaZeta: " << deltaZ << " cm" << "\n";
-
-
- start:
-
- hITSZv->Add(hITSZv,-1);
- hITSXv->Add(hITSXv,-1);
- hITSYv->Add(hITSYv,-1);
-
- np1=np2=0;
-
-
- for(i=g2->GetStartSPD();i<=g2->GetLastSPD();i++)
- {
- fITS->ResetRecPoints();
- itsloader->TreeR()->GetEvent(i);
- npoints = recpoints->GetEntries();
- for (ipoint=0;ipoint<npoints;ipoint++) {
-
- pnt = (AliITSRecPoint*)recpoints->UncheckedAt(ipoint);
- l[0]=pnt->GetX();
- l[1]=0;
- l[2]=pnt->GetZ();
- g2->LtoG(i, l, p);
-
- if(i<80 && TMath::Abs(p[2])<14.35) {
- p[0]=p[0]-vzero[0];
- p[1]=p[1]-vzero[1];
- r=TMath::Sqrt(TMath::Power(p[0],2)+TMath::Power(p[1],2));
- y1[np1]=p[1];
- x1[np1]=p[0];
- z1[np1]=p[2];
- r1[np1]=r;
- phi1[np1]=PhiFunc(p);
- np1++;
- }