void readVertex(Int_t evNumber=1) { // Dynamically link some shared libs if (gClassTable->GetID("AliRun") < 0) { gROOT->LoadMacro("loadlibs.C"); loadlibs(); } // Connect the Root Galice file containing Geometry, Kine and Hits TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root"); if (!file) file = new TFile("~/w0/START/galice.root","UPDATE"); // Get AliRun object from file or create it if not on file if (!gAlice) { gAlice = (AliRun*)file->Get("gAlice"); if (gAlice) printf("AliRun object found on file\n"); if (!gAlice) gAlice = new AliRun("gAlice","Alice test program"); } char nameTD[8],nameTR[8]; TH1F *hVertex = new TH1F("hVertex","Z position of vertex",100,-100,100); TH1F *hRealVertex = new TH1F("hRealVertex","Z position of vertex", 100,-100,100); digits = new AliSTARTdigit(); vertex = new AliSTARTvertex(); // Event ------------------------- LOOP for (Int_t j=0; jGet(nameTR); cout<<" tr "<Read(nameTR); hVertex->Fill(vertex->GetVertex()); printf(" Z position %f\n",vertex->GetVertex()); gAlice->GetEvent(j); AliHeader *header = gAlice->GetHeader(); AliGenEventHeader* genHeader = header->GenEventHeader(); TArrayF *o = new TArrayF(3); genHeader->PrimaryVertex(*o); cout<<" o "<At(2); cout<<" zRealVertex "<Fill(zRealVertex); } Hfile = new TFile("figs.root","RECREATE","Histograms for START Vertex"); printf("Writting histograms to root file \n"); Hfile->cd(); //Create a canvas, set the view range, show histograms gStyle->SetOptStat(111111); hVertex->SetXTitle("vertex position, mm"); hVertex->SetYTitle("number of events"); hVertex->Write(); hRealVertex->SetXTitle("vertex position, mm"); hRealVertex->SetYTitle("number of events"); hRealVertex->Write(); } // end of macro