void readESD() { //read T0 ESD from reconstructed cosmic data char filename[100]; sprintf(filename,""); TFile *file = new TFile("AliESDs.root"); // TH1F *hMean= new TH1F("hMean"," T0 ", 10000,1.249e6,1.259e6); TH1F *hMean= new TH1F("hMean"," T0 ",6000,0,6000); TH1F *hVertex = new TH1F("hVertex","Z position of vertex", 100,-120,-110); // TH1F *hAmp = new TH1F("hAmp"," amplitude", 2500, 10,2510); // TH1F *hTime = new TH1F("hTime"," time", 500,1250000,1250500); TH1F *hAmp[24]; TH1F * hTime[24]; Char_t buf1[20], buf2[20]; for(Int_t ic=0; ic<24; ic++) { sprintf(buf1,"Amp%i",ic+1); hAmp[ic] = new TH1F(buf1,"Amp",50, 0,50 ); sprintf(buf2,"Time%i",ic+1); // hTime[ic] = new TH1F("hTime"," time", 10000,1.249e6,1.259e6); hTime[ic] = new TH1F("hTime"," time", 6000,0,6000); } Double_t *amp, *time; TTree *esdTree = (TTree*)file->Get("esdTree"); esdTree->ls(); TBranch *brRec=esdTree->GetBranch("AliESDTZERO."); AliESDTZERO *fesdT0 = new AliESDTZERO(); if (brRec) { brRec->SetAddress(&fesdT0); }else{ cerr<<"EXEC Branch T0 not found"<GetEntry(0); // Event ------------------------- LOOP for (Int_t ievent=0; ieventGetEntries(); ievent++){ brRec->GetEntry(ievent); Int_t mean = fesdT0->GetT0(); hMean->Fill(mean); Float_t vertex= fesdT0->GetT0zVertex(); if (ievent <100) cout<Fill(vertex); amp=fesdT0->GetT0amplitude(); time=fesdT0->GetT0time(); for (Int_t i=0; i<24; i++){ hAmp[i]->Fill(amp[i]); hTime[i]->Fill(time[i]); if (ievent <100) cout<cd(); //Create a canvas, set the view range, show histograms // g // gStyle->SetOptStat(111111); hVertex->Write(); hMean->Write(); for (Int_t i=0; i<24; i++){ hAmp[i]->Write(); hTime[i]->Write(); } TCanvas *c1 = new TCanvas("c1", "Time C side",0,48,1280,951); c1->Divide(4,3); // gStyle->SetOptFit(1111); for (Int_t i=0; i<12; i++) { c1->cd(i+1); hTime[i]->Draw(); } TCanvas *c2 = new TCanvas("c2", "Time A side",0,48,1280,951); c2->Divide(4,3); gStyle->SetOptFit(1111); for (Int_t i=12; i<24; i++) { c2->cd(i+1-12); hTime[i]->Draw(); } TCanvas *c3 = new TCanvas("c3", "Time mean & vertex",0,48,1280,951); c3->Divide(1,2); c3->cd(1); hVertex->Draw(); c3->cd(2); hMean->Draw(); } // end of macro