--- /dev/null
+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"<<endl;
+ return;
+ }
+ // brRec->GetEntry(0);
+
+ // Event ------------------------- LOOP
+ for (Int_t ievent=0; ievent<brRec->GetEntries(); ievent++){
+ brRec->GetEntry(ievent);
+ Int_t mean = fesdT0->GetT0();
+ hMean->Fill(mean);
+ Float_t vertex= fesdT0->GetT0zVertex();
+ if (ievent <100) cout<<ievent<<" "<<mean<<" vertex "<<vertex<<endl;
+ if(vertex<999){
+ hVertex->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<<ievent<<" "<<i<<" time "<<time[i]<<" amp "<<amp[i]<<endl;
+ }
+ }
+ }
+
+ Hfile = new TFile("FigESD_25081.root","RECREATE","Histograms for T0
+digits");
+ printf("Writting histograms to root file \n");
+ Hfile->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