- hDx ->Write() ;
- hDz ->Write() ;
- hDr ->Write() ;
- hNrp ->Write() ;
- hNrpX->Write() ;
- hNrpZ->Write() ;
-
- // Plot histograms
-
- TCanvas *cpvCanvas = new TCanvas("CPV","CPV analysis",20,20,800,400);
- gStyle->SetOptStat(111111);
- gStyle->SetOptFit(1);
- gStyle->SetOptDate(1);
- cpvCanvas->Divide(3,2);
-
- cpvCanvas->cd(1);
- gPad->SetFillColor(10);
- hNrp->SetFillColor(16);
- hNrp->Draw();
-
- cpvCanvas->cd(2);
- gPad->SetFillColor(10);
- hNrpX->SetFillColor(16);
- hNrpX->Draw();
-
- cpvCanvas->cd(3);
- gPad->SetFillColor(10);
- hNrpZ->SetFillColor(16);
- hNrpZ->Draw();
-
- cpvCanvas->cd(4);
- gPad->SetFillColor(10);
- hDx->SetFillColor(16);
- hDx->Fit("gaus");
- hDx->Draw();
-
- cpvCanvas->cd(5);
- gPad->SetFillColor(10);
- hDz->SetFillColor(16);
- hDz->Fit("gaus");
- hDz->Draw();
-
- cpvCanvas->cd(6);
- gPad->SetFillColor(10);
- hDr->SetFillColor(16);
- hDr->Draw();
-
- cpvCanvas->Print("CPV.ps");
+// for (Int_t itrack=0; itrack<ntracks; itrack++) {
+// // Get the Hits Tree for the Primary track itrack
+// gAlice->ResetHits();
+// gAlice->TreeH()->GetEvent(itrack);
+// for (Int_t iModule=0; iModule < nOfModules; iModule++) {
+// cpvModule = fPHOS->GetCPVModule(iModule);
+// cpvHits = cpvModule.Hits();
+// nCPVhits = cpvHits->GetEntriesFast();
+// for (Int_t ihit=0; ihit<nCPVhits; ihit++) {
+// cpvHit = (AliPHOSCPVHit*)cpvHits->UncheckedAt(ihit);
+// p = cpvHit->GetMomentum();
+// xzgen[0] = cpvHit->X();
+// xzgen[1] = cpvHit->Y();
+// ipart = cpvHit->GetIpart();
+// TClonesArray &lhits = *(TClonesArray *)hitsPerModule[iModule];
+// new(lhits[hitsPerModule[iModule]->GetEntriesFast()]) AliPHOSCPVHit(*cpvHit);
+// }
+// cpvModule.Clear();
+// }
+// }
+// for (iModule=0; iModule < nOfModules; iModule++) {
+// Int_t nsum = hitsPerModule[iModule]->GetEntriesFast();
+// printf("Module %d has %d hits\n",iModule,nsum);
+// }
+
+// // Then go through reconstructed points and for each find
+// // the closeset hit
+// // The distance from the rec.point to the closest hit
+// // gives the coordinate resolution of the CPV
+
+// // Get the Reconstruction Tree
+// gAlice->TreeR()->GetEvent(0) ;
+// TIter nextRP(*fPHOS->PpsdRecPoints() ) ;
+// AliPHOSCpvRecPoint *cpvRecPoint ;
+// Float_t xgen, zgen;
+// while( ( cpvRecPoint = (AliPHOSCpvRecPoint *)nextRP() ) ) {
+// TVector3 locpos;
+// cpvRecPoint->GetLocalPosition(locpos);
+// Int_t phosModule = cpvRecPoint->GetPHOSMod();
+// Int_t rpMult = cpvRecPoint->GetDigitsMultiplicity();
+// Int_t rpMultX, rpMultZ;
+// cpvRecPoint->GetClusterLengths(rpMultX,rpMultZ);
+// Float_t xrec = locpos.X();
+// Float_t zrec = locpos.Z();
+// Float_t dxmin = 1.e+10;
+// Float_t dzmin = 1.e+10;
+// Float_t r2min = 1.e+10;
+// Float_t r2;
+
+// cpvHits = hitsPerModule[phosModule-1];
+// Int_t nCPVhits = cpvHits->GetEntriesFast();
+// for (Int_t ihit=0; ihit<nCPVhits; ihit++) {
+// cpvHit = (AliPHOSCPVHit*)cpvHits->UncheckedAt(ihit);
+// xgen = cpvHit->X();
+// zgen = cpvHit->Y();
+// r2 = TMath::Power((xgen-xrec),2) + TMath::Power((zgen-zrec),2);
+// if ( r2 < r2min ) {
+// r2min = r2;
+// dxmin = xgen - xrec;
+// dzmin = zgen - zrec;
+// }
+// }
+// hDx ->Fill(dxmin);
+// hDz ->Fill(dzmin);
+// hDr ->Fill(TMath::Sqrt(r2min));
+// hNrp ->Fill(rpMult);
+// hNrpX->Fill(rpMultX);
+// hNrpZ->Fill(rpMultZ);
+// }
+// delete [] hitsPerModule;
+// }
+// // Save histograms
+
+// Text_t outputname[80] ;
+// sprintf(outputname,"%s.analyzed",fRootFile->GetName());
+// TFile output(outputname,"RECREATE");
+// output.cd();
+
+// hDx ->Write() ;
+// hDz ->Write() ;
+// hDr ->Write() ;
+// hNrp ->Write() ;
+// hNrpX->Write() ;
+// hNrpZ->Write() ;
+
+// // Plot histograms
+
+// TCanvas *cpvCanvas = new TCanvas("CPV","CPV analysis",20,20,800,400);
+// gStyle->SetOptStat(111111);
+// gStyle->SetOptFit(1);
+// gStyle->SetOptDate(1);
+// cpvCanvas->Divide(3,2);
+
+// cpvCanvas->cd(1);
+// gPad->SetFillColor(10);
+// hNrp->SetFillColor(16);
+// hNrp->Draw();
+
+// cpvCanvas->cd(2);
+// gPad->SetFillColor(10);
+// hNrpX->SetFillColor(16);
+// hNrpX->Draw();
+
+// cpvCanvas->cd(3);
+// gPad->SetFillColor(10);
+// hNrpZ->SetFillColor(16);
+// hNrpZ->Draw();
+
+// cpvCanvas->cd(4);
+// gPad->SetFillColor(10);
+// hDx->SetFillColor(16);
+// hDx->Fit("gaus");
+// hDx->Draw();
+
+// cpvCanvas->cd(5);
+// gPad->SetFillColor(10);
+// hDz->SetFillColor(16);
+// hDz->Fit("gaus");
+// hDz->Draw();
+
+// cpvCanvas->cd(6);
+// gPad->SetFillColor(10);
+// hDr->SetFillColor(16);
+// hDr->Draw();
+
+// cpvCanvas->Print("CPV.ps");