void readRec() { //read T0 RecPoints and plots histos // Dynamically link some shared libs /* if (gClassTable->GetID("AliRun") < 0) { gROOT->LoadMacro("loadlibs.C"); loadlibs(); } */ char filename[100]; sprintf(filename,"galice.root"); AliRunLoader* rl = AliRunLoader::Open("galice.root",AliConfig::GetDefaultEventFolderName(),"read"); if (rl == 0x0) { cerr<<"Can not open session for file galice.root\n"; return; } rl->LoadgAlice(); gAlice = rl->GetAliRun(); AliT0* T0 = (AliT0 *)gAlice->GetDetector("T0"); rl->LoadHeader(); Int_t retval; AliLoader* lstart = rl->GetLoader("T0Loader"); Int_t iNevents=rl->GetNumberOfEvents(); cout<<" nevents "<GetEvent(ievent); AliHeader *header = gAlice->GetHeader(); AliGenEventHeader* genHeader = header->GenEventHeader(); TArrayF *o = new TArrayF(3); genHeader->PrimaryVertex(*o); Float_t zRealVertex=o->At(2); hRealVertex->Fill( zRealVertex); lstart->LoadRecPoints("READ"); TTree *recTree = lstart->TreeR(); TBranch *brRec=recTree->GetBranch("T0"); AliT0RecPoint *fRec = new AliT0RecPoint(); if (brRec) { brRec->SetAddress(&fRec); }else{ cerr<<"EXEC Branch T0 Recpoints not found"<GetEntry(0); Int_t mean = fRec->GetMeanTime(); hMean->Fill(mean/1000.); Int_t besttimeleft = fRec->GetBestTimeLeft(); Int_t besttimeright = fRec->GetBestTimeRight(); hBestTimeC->Fill(0.001 * besttimeright); hBestTimeA->Fill(0.001 * besttimeleft ); Float_t vertex= fRec->GetVertex(); if(vertex<99){ cout<Fill(zRealVertex-vertex); hVertex->Fill(vertex); for (Int_t i=0; i<24; i++){ hAmp->Fill(fRec->GetAmp(i)); hTime->Fill(fRec->GetTime(i)); // cout<<"time "<GetTime(i)<<" amp "<GetAmp(i)<cd(); //Create a canvas, set the view range, show histograms gStyle->SetOptStat(111111); // TCanvas *c1 = new TCanvas("c1","Alice T0 Time ",400,10,600,600); // hTimediff->Write(); hBestTimeC->Write(); hBestTimeA ->Write(); hVertex->Write(); hRealVertex->Write(); hAcc->Write(); hMean->Write(); hAmp->Write(); hTime->Write(); } // end of macro