#include "AliRun.h"
#include "AliRunLoader.h"
#include "AliPDG.h"
-#include "STEER/AliESDEvent.h"
-#include "STEER/AliESDZDC.h"
+#include "STEER/ESD/AliESDEvent.h"
+#include "STEER/ESD/AliESDZDC.h"
#endif
-void CheckAlienZDCESD(Int_t year=2010, const char* period="10a",
+void CheckAlienZDCESD(Int_t year=2010, const char* period="10f",
Int_t nRun=0, Int_t recoPass=1, Int_t nMaxFiles=1,
- Bool_t plot=kTRUE, Bool_t esdWordCut=kFALSE)
+ Bool_t plot=kFALSE, Bool_t esdWordCut=kFALSE)
{
if(nRun==0){
// Histogram definition
// ----------------------------------------------------------------
-/* TH2F *centroidZNsideC = new TH2F("centroidZNsideC","Impact point over ZNC",100,-5.,5.,100,-5.,5.);
+ TH2F *centroidZNsideC = new TH2F("centroidZNsideC","Impact point over ZNC",100,-3.5,3.5,100,-3.5,3.5);
centroidZNsideC->SetXTitle("X_{ZNC} (cm)");
centroidZNsideC->SetYTitle("Y_{ZNC} (cm)");
- TH2F * centroidZNsideA = new TH2F("centroidZNsideA","Impact point over ZNA",100,-5.,5.,100,-5.,5.);
+ TH2F * centroidZNsideA = new TH2F("centroidZNsideA","Impact point over ZNA",100,-3.5,3.5,100,-3.5,3.5);
centroidZNsideA->SetXTitle("X_{ZNA} (cm)");
centroidZNsideA->SetYTitle("Y_{ZNA} (cm)");
-*/
- TH1F * enZNC = new TH1F("enZNC", "Energy in ZNC",100,0.,5000.);
+ TH1F * enZNC = new TH1F("enZNC", "ZNC signal",100,0.,2000.);
enZNC->SetXTitle("E (GeV)");
- TH1F * enZPC = new TH1F("enZPC", "Energy in ZPC",100,0.,5000.);
+ TH1F * enZPC = new TH1F("enZPC", "ZPC signal",100,0.,2000.);
enZPC->SetXTitle("E (GeV)");
- TH1F * enZNA = new TH1F("enZNA", "Energy in ZNA",100,0.,5000.);
+ TH1F * enZNA = new TH1F("enZNA", "ZNA signal",100,0.,2000.);
enZNA->SetXTitle("E (GeV)");
- TH1F * enZPA = new TH1F("enZPA", "Energy in ZPA",100,0.,5000.);
+ TH1F * enZPA = new TH1F("enZPA", "ZPA signal",100,0.,2000.);
enZPA->SetXTitle("E (GeV)");
- TH1D * enZEM1 = new TH1D("enZEM1", "Energy in ZEM1",100,0.,100.);
+ TH1D * enZEM1 = new TH1D("enZEM1", "Energy in ZEM1",100,0.,2000.);
enZEM1->SetXTitle("E (GeV)");
- TH1D * enZEM2 = new TH1D("enZEM2", "Energy in ZEM2",100,0.,100.);
+ TH1D * enZEM2 = new TH1D("enZEM2", "Energy in ZEM2",100,0.,2000.);
enZEM2->SetXTitle("E (GeV)");
// ----------------------------------------------------------------
TH1D * hZNCTow[5]; TH1D * hZPCTow[5];
char nomehistznc[30], nomehistzpc[30], nomehistzna[30], nomehistzpa[30];
for(Int_t i=0; i<5; i++){
sprintf(nomehistznc,"ZNC-pm%d",i);
- hZNCTow[i] = new TH1D(nomehistznc, nomehistznc, 100, 0.,100.);
+ hZNCTow[i] = new TH1D(nomehistznc, nomehistznc, 100, 0.,1200.);
sprintf(nomehistzpc,"ZPC-pm%d",i);
- hZPCTow[i] = new TH1D(nomehistzpc, nomehistzpc, 100, 0.,100.);
+ hZPCTow[i] = new TH1D(nomehistzpc, nomehistzpc, 100, 0.,1200.);
sprintf(nomehistzna,"ZNA-pm%d",i);
- hZNATow[i] = new TH1D(nomehistzna, nomehistzna, 100, 0.,100.);
+ hZNATow[i] = new TH1D(nomehistzna, nomehistzna, 100, 0.,1200.);
sprintf(nomehistzpa,"ZPA-pm%d",i);
- hZPATow[i] = new TH1D(nomehistzpa, nomehistzpa, 100, 0.,100.);
+ hZPATow[i] = new TH1D(nomehistzpa, nomehistzpa, 100, 0.,1200.);
}
-
- TH1D *hSumQZNC = new TH1D("hSumQZNC", "hSumQZNC", 100, 0., 1000.);
- TH1D *hSumQZPC = new TH1D("hSumQZPC", "hSumQZPC", 100, 0., 1000.);
- TH1D *hSumQZNA = new TH1D("hSumQZNA", "hSumQZNA", 100, 0., 1000.);
- TH1D *hSumQZPA = new TH1D("hSumQZPA", "hSumQZPA", 100, 0., 1000.);
//
TH1F *hESDword = new TH1F("hESDword","hESDword",6,0.5,6.5);
hESDword->SetXTitle("ZDC trigger pattern");
+
+ // ----------------------------------------------------------------
+ // %%%%% TDCs
+ // ----------------------------------------------------------------
+ TH1F * hTDC[11];
+ for(Int_t itdc=0; itdc<11; itdc++){
+
+ if(itdc==0) hTDC[itdc] = new TH1F("TDCZNC", "TDC ZNC", 200, -150., 50.);
+ else if(itdc==1) hTDC[itdc] = new TH1F("TDCZNA", "TDC ZNA", 200, -150., 50.);
+ else if(itdc==2) hTDC[itdc] = new TH1F("TDCL0", "TDC L0", 100, 250., 350.);
+ else if(itdc==3) hTDC[itdc] = new TH1F("TDCADCgate", "TDC ADCgate", 100, 500., 600.);
+ else if(itdc==4) hTDC[itdc] = new TH1F("TDCZNsum", "TDC ZNsum", 60, -100., -40.);
+ else if(itdc==5) hTDC[itdc] = new TH1F("TDCZNdiff", "TDC ZNdiff", 60, -30., 30.);
+ else if(itdc==6) hTDC[itdc] = new TH1F("TDCZPC", "TDC ZPC", 200, -150., 50.);
+ else if(itdc==7) hTDC[itdc] = new TH1F("TDCZPA", "TDC ZPA", 200, -150., 50.);
+ else if(itdc==8) hTDC[itdc] = new TH1F("TDCZEM1", "TDC ZEM1", 200, -150., 50.);
+ else if(itdc==9) hTDC[itdc] = new TH1F("TDCZEM2", "TDC ZEM2", 200, -150., 50.);
+ else if(itdc==10) hTDC[itdc] = new TH1F("TDCorZN", "TDC OR ZN", 400, -150., 50.);
+
+ hTDC[itdc]->SetXTitle("TDC (ns)");
+ }
+
+ TH2F *hDebunch = new TH2F("hDebunch","hDebunch",240,-30.,30.,240,-30.,30.);
+ hDebunch->SetYTitle("t_{ZNC}+t_{ZNA} (ns)");
+ hDebunch->SetXTitle("t_{ZNC}-t_{ZNA} (ns)");
+
+ TH1F * hTDCcorr[4];
+ for(Int_t itdc=0; itdc<4; itdc++){
+ if(itdc==0) hTDCcorr[itdc] = new TH1F("TDCZNCcorr", "TDC ZNC corr", 200, -100., 100.);
+ else if(itdc==1) hTDCcorr[itdc] = new TH1F("TDCZNAcorr", "TDC ZNA corr", 200, -100., 100.);
+ else if(itdc==2) hTDCcorr[itdc] = new TH1F("TDCDiffcorr", "TDC Diff corr", 100, -50., 50.);
+ else if(itdc==3) hTDCcorr[itdc] = new TH1F("TDCSumcorr", "TDC Sum corr", 100, -50., 50.);
+
+ hTDCcorr[itdc]->SetXTitle("TDC (ns)");
+ }
+
TGrid::Connect("alien:",0,0,"t");
gSystem->Exec(Form("gbbox find \"/alice/data/%d/LHC%s/000%d/ESDs/pass%d\" \"AliESDs.root\" > ESDFiles.txt",
year, period, nRun, recoPass));
+ //gSystem->Exec(Form("gbbox find \"/alice/data/%d/LHC%s/000%d/ESDs/hlt_clustering\" \"AliESDs.root\" > ESDFiles.txt",
+ // year, period, nRun));
FILE* listruns=fopen("ESDFiles.txt","r");
char esdFileName[200], filnamalien[200];
sprintf(directory,"/alice/data/%d",year);
if(!strstr(esdFileName,directory)) continue;
sscanf(esdFileName,"/alice/data/%d/LHC%s/000%d/ESDs/pass%d/%s/AliESDs.root",&year,&yperiod,&nRun,&recoPass,&dirESD);
+ //sscanf(esdFileName,"/alice/data/%d/LHC%s/000%d/ESDs/hlt_clustering/%s/AliESDs.root",&year,&yperiod,&nRun,&dirESD);
sprintf(filnamalien,"alien://%s",esdFileName);
printf("\n Opening file: %s\n",filnamalien);
// get the event summary data
tree->GetEvent(iEvent);
- //printf(" ev. type %d\n",esd->GetEventType());
//
if(esd->GetEventType() == 7){
- nevPhys++;
if(!esd) {
Error("CheckESD", "no ESD object found for event %d", iEvent);
return;
Double_t sumQznc=0., sumQzpc=0., sumQzna=0., sumQzpa=0.;
AliESDZDC *esdZDC = esd->GetESDZDC();
- //Double_t centrZNC={-999.,-999.}, centrZNA={-999.,-999.};
- //esdZDC->GetZNCentroidInpp(centrZNC, centrZNA);
- //Short_t npart = esdZDC->GetZDCParticipants();
+ //printf(" ev.%d \n",nevPhys);
+ //esdZDC->Print("");
+ nevPhys++;
+
+
+ Double_t centrZNC[2]={-999.,-999.}, centrZNA[2]={-999.,-999.};
+ esdZDC->GetZNCentroidInPbPb(2760.,centrZNC, centrZNA);
+ //
Double_t energyZNC = esdZDC->GetZDCN1Energy();
Double_t energyZPC = esdZDC->GetZDCP1Energy();
Double_t energyZNA = esdZDC->GetZDCN2Energy();
nevZPC++;
hESDword->Fill(6.);
}
- //if(centrZNC[0]!=-999. && centrZNC[1]!=-999) centroidZNsideC->Fill(centrZNC[0], centrZNC[1]);
- //if(centrZNA[0]!=-999. && centrZNA[1]!=-999) centroidZNsideA->Fill(centrZNA[0], centrZNA[1]);
enZNC->Fill(energyZNC);
enZPC->Fill(energyZPC);
enZNA->Fill(energyZNA);
if((iWord & 0x00000020) == 0x00000020) sumQzpc += towZPC[jj];
if((iWord & 0x00000001) == 0x00000001) sumQzna += towZNA[jj];
if((iWord & 0x00000002) == 0x00000002) sumQzpa += towZPA[jj];
- //
- /*w = TMath::Power(towZNC[jj], alpha);
- numX += x[jj]*w;
- numY += y[jj]*w;
- den += w;*/
}
else{
sumQznc += towZNC[jj];
}
}
//
- if(esdWordCut){
+
+ if(centrZNA[0]!=-999. && centrZNA[1]!=-999) centroidZNsideA->Fill(centrZNA[0], centrZNA[1], towZNA[0]);
+ if(centrZNC[0]!=-999. && centrZNC[1]!=-999) centroidZNsideC->Fill(centrZNC[0], centrZNC[1], towZNC[0]);
+/* if(esdWordCut){
if((iWord & 0x00000010) == 0x00000010) {
hSumQZNC->Fill(sumQznc);
- //
- //if(centrZNA[0]!=-999. && centrZNA[1]!=-999) centroidZNsideC->Fill(centrZNA[0], centrZNA[1]);
}
//
if((iWord & 0x00000020) == 0x00000020) {
hSumQZPC->Fill(sumQzpc);
}
+
}
else{
hSumQZNC->Fill(sumQznc);
//
hSumQZNA->Fill(sumQzna);
hSumQZPA->Fill(sumQzpa);
- }
+*/
+
+ Int_t tdcData[32][4];
+ Float_t tdcCorr[32][4];
+ for(int ij=0; ij<32; ij++){
+ for(int ik=0; ik<4;ik++){
+ tdcData[ij][ik] = 9999.;
+ tdcCorr[ij][ik] = -9999./0.025;
+ }
+ }
+ Float_t tdcC=-9999./0.025, tdcA=-9999./0.025, tdcGate=-9999./0.025, tdcL0=-9999./0.025;
+ Float_t tdcDiff=-9999./0.025, tdcSum=-9999./0.025;
+ Float_t tdcZPC=-9999./0.025, tdcZPA=-9999./0.025;
+ Float_t tdcZEM1=-9999./0.025, tdcZEM2=-9999./0.025;
+ Float_t tdcZNCcorr=-9999./0.025, tdcZNAcorr=-9999./0.025;
+ //
+ for(Int_t itdc=0; itdc<21; itdc++){
+ for(Int_t j=0; j<4; j++){
+
+ tdcData[itdc][j] = esdZDC->GetZDCTDCData(itdc, j);
+ if(tdcData[itdc][j] == 0.) continue; //change in 9999
+
+ //printf(" TDC ch.%d hit %d: %f ns\n",itdc,j,0.025*esdZDC->GetZDCTDCData(itdc, j));
+ tdcCorr[itdc][j] = esdZDC->GetZDCTDCCorrected(itdc, j);
+ //printf(" %f ns\n",esdZDC->GetZDCTDCCorrected(itdc, j));
+ if(itdc==12 && j==0 && (tdcData[10][0]!=0. && tdcData[12][0]!=0.)){ //change in 9999!!!
+ tdcDiff = esdZDC->GetZDCTimeDiff();
+ tdcSum = esdZDC->GetZDCTimeSum();
+ }
+
+ if(itdc==8 && tdcData[itdc][j]!=0.){
+ tdcZEM1 = 0.025*(tdcData[itdc][j]);
+ }
+ else if(itdc==9 && tdcData[itdc][j]!=0.){
+ tdcZEM2 = 0.025*(tdcData[itdc][j]);
+ }
+ else if(itdc==10 && tdcData[itdc][j]!=0.){
+ tdcC = 0.025*(tdcData[itdc][j]);
+ tdcZNCcorr = tdcCorr[itdc][j];
+ hTDCcorr[0]->Fill(tdcZNCcorr);
+ }
+ else if(itdc==11 && tdcData[itdc][j]!=0.){
+ tdcZPC = 0.025*(tdcData[itdc][j]);
+ }
+ else if(itdc==12 && tdcData[itdc][j]!=0.){
+ tdcA = 0.025*(tdcData[itdc][j]);
+ tdcZNAcorr = tdcCorr[itdc][j];
+ hTDCcorr[1]->Fill(tdcZNAcorr);
+ }
+ else if(itdc==13 && tdcData[itdc][j]!=0.){
+ tdcZPA = 0.025*(tdcData[itdc][j]);
+ }
+ else if(itdc==14 && tdcData[itdc][j]!=0.){
+ tdcGate = 0.025*tdcData[itdc][j];
+ hTDC[3]->Fill(tdcGate);
+ }
+ else if(itdc==15 && tdcData[itdc][j]!=0.){
+ tdcL0 = 0.025*tdcData[itdc][j];
+ hTDC[2]->Fill(tdcL0);
+ //
+ if(tdcC!=-9999./0.025) hTDC[0]->Fill(tdcC-tdcL0);
+ if(tdcA!=-9999./0.025) hTDC[1]->Fill(tdcA-tdcL0);
+ if(tdcZPC!=-9999./0.025) hTDC[6]->Fill(tdcZPC-tdcL0);
+ if(tdcZPA!=-9999./0.025) hTDC[7]->Fill(tdcZPA-tdcL0);
+ if(tdcZEM1!=-9999./0.025) hTDC[8]->Fill(tdcZEM1-tdcL0);
+ if(tdcZEM2!=-9999./0.025) hTDC[9]->Fill(tdcZEM2-tdcL0);
+ if(tdcC!=-9999./0.025 && tdcA!=-9999./0.025 && tdcL0!=-9999./0.025){
+ tdcC = tdcC-tdcL0;
+ tdcA = tdcA-tdcL0;
+ hTDC[4]->Fill(tdcC+tdcA);
+ hTDC[5]->Fill(tdcC-tdcA);
+ //printf(" TDCdata: sum %f diff %f ", tdcC+tdcA, tdcC-tdcA);
+ //hDebunch->Fill(tdcC-tdcA+2.1, tdcC+tdcA+65.5);
+ }
+ if(tdcZNAcorr!=-9999./0.025 && tdcZNCcorr!=-9999./0.025 && tdcL0!=-9999./0.025){
+ hTDCcorr[2]->Fill(tdcZNCcorr-tdcZNAcorr);
+ hTDCcorr[3]->Fill(tdcZNCcorr+tdcZNAcorr);
+ hDebunch->Fill(tdcDiff, tdcSum);
+ //printf(" TDCcorr: sum %f diff %f (ns) \n", tdcSum, tdcDiff);
+ }
+ }
+ else if(itdc==18 && tdcData[itdc][j]!=0.){
+ hTDC[10]->Fill(0.025*tdcData[itdc][j]-tdcL0);
+ }
+ }
+ }//tdc scan
+
+
+ }//physics event
+
- }
+ }//event loop
nAnalyzedFiles++;
esdFile->Close();
}//if(nAnalyzedFiles<=nMaxFiles)
else{
- printf("\t %d files analyzed\n\n",nMaxFiles);
- break;
+ printf("\t %d files analyzed\n\n",nMaxFiles);
+ break;
}
} // while closing
- printf(" No. of events over threshold: ZNA %d ZPA %d ZEM1 %d ZEM2 %d ZNC %d ZPC %d\n\n",
- nevZNA, nevZPA, nevZEM1, nevZEM2, nevZNC, nevZPC);
+/* printf(" No. of events over threshold: ZNA: %d ZPA: %d ZEM1: %d "
+ " ZEM2: %d ZNC: %d ZPC: %d\n\n",
+ nevZNA, nevZPA, nevZEM1, nevZEM2, nevZNC, nevZPC);*/
+
+ TFile * fileout = new TFile("ESDhistos.root","recreate");
+ fileout->cd();
+ centroidZNsideC->Write();
+ centroidZNsideA->Write();
+ enZNC->Write();
+ enZNA->Write();
+ enZPC->Write();
+ enZPA->Write();
+ enZEM1->Write();
+ enZEM2->Write();
+ for(Int_t jj=0; jj<5; jj++){
+ hZNCTow[jj]->Write();
+ hZPCTow[jj]->Write();
+ hZNATow[jj]->Write();
+ hZPATow[jj]->Write();
+ }
+ /*hSumQZNC->Write();
+ hSumQZPC->Write();
+ hSumQZNA->Write();
+ hSumQZPA->Write();*/
+ //
+ hESDword->Write();
+ for(Int_t jj=0; jj<11; jj++) hTDC[jj]->Write();
+ for(Int_t jj=0; jj<4; jj++) hTDCcorr[jj]->Write();
+ //
+ fileout->Close();
+
if(plot){
//***********************************************************
// #### ROOT initialization
gStyle->SetFrameFillColor(10);
gStyle->SetOptTitle(1);
if(esdWordCut) gStyle->SetOptStat(1111111);
- else gStyle->SetOptStat(11111);
+ else gStyle->SetOptStat(1111);
gStyle->SetOptFit(0);
gStyle->SetTitleTextColor(4);
gStyle->SetStatTextColor(4);
//-------------------------------------------------
- TCanvas *c1 = new TCanvas("c1","Hadronic PMCs + ZEMs",400,0,500,800);
+ TCanvas *c1 = new TCanvas("c1","ZDCs + ZEMs signals",400,0,500,800);
c1->Divide(2,3);
c1->cd(1);
- hZNATow[0]->SetLineColor(kBlue);
- hZNATow[0]->SetFillColor(kBlue);
- hZNATow[0]->Draw("");
+ gPad->SetLogy(1);
+ enZNC->Draw("");
+ enZNC->SetLineColor(kBlue);
+ enZNC->SetFillColor(kBlue);
c1->cd(2);
- hZPATow[0]->SetLineColor(kBlue);
- hZPATow[0]->SetFillColor(kBlue);
- hZPATow[0]->Draw("");
+ gPad->SetLogy(1);
+ enZPC->Draw("");
+ enZPC->SetLineColor(kBlue+3);
+ enZPC->SetFillColor(kBlue+3);
c1->cd(3);
+ gPad->SetLogy(1);
enZEM1->SetLineColor(kRed);
enZEM1->SetFillColor(kRed);
enZEM1->Draw("");
c1->cd(4);
+ gPad->SetLogy(1);
enZEM2->SetLineColor(kRed);
enZEM2->SetFillColor(kRed);
enZEM2->Draw("");
c1->cd(5);
- hZNCTow[0]->Draw("");
- c1->cd(6);
- hZPCTow[0]->Draw("");
-
- //-------------------------------------------------
- TCanvas *c2 = new TCanvas("c2","Total signal in hadronic ZDCs",0,0,500,500);
- c2->Divide(2,2);
- c2->cd(1);
- enZNC->Draw("");
- enZNC->SetLineColor(kBlue);
- enZNC->SetFillColor(kBlue);
- c2->cd(2);
- enZPC->Draw("");
- enZPC->SetLineColor(kBlue+3);
- enZPC->SetFillColor(kBlue+3);
- c2->cd(3);
+ gPad->SetLogy(1);
enZNA->Draw("");
enZNA->SetLineColor(kRed);
enZNA->SetFillColor(kRed);
- c2->cd(4);
+ c1->cd(6);
+ gPad->SetLogy(1);
enZPA->Draw("");
enZPA->SetLineColor(kRed+1);
- enZPA->SetFillColor(kRed+1);
-
+ enZPA->SetFillColor(kRed+1);
//-------------------------------------------------
TCanvas *c3 = new TCanvas("c3","Side C ZDCs",0,0,800,400);
c3->Divide(5,2);
c3->cd(1);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZNCTow[0]->SetLineColor(kBlue);
hZNCTow[0]->SetFillColor(kBlue);
hZNCTow[0]->Draw("");
hZNCTow[1]->SetFillColor(kBlue);
hZNCTow[1]->Draw("");
c3->cd(3);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZNCTow[2]->SetLineColor(kBlue);
hZNCTow[2]->SetFillColor(kBlue);
hZNCTow[2]->Draw("");
c3->cd(4);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZNCTow[3]->SetLineColor(kBlue);
hZNCTow[3]->SetFillColor(kBlue);
hZNCTow[3]->Draw("");
c3->cd(5);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZNCTow[4]->SetLineColor(kBlue);
hZNCTow[4]->SetFillColor(kBlue);
hZNCTow[4]->Draw("");
//
c3->cd(6);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZPCTow[0]->SetLineColor(kBlue+3);
hZPCTow[0]->SetFillColor(kBlue+3);
hZPCTow[0]->Draw("");
c3->cd(7);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZPCTow[1]->SetLineColor(kBlue+3);
hZPCTow[1]->SetFillColor(kBlue+3);
hZPCTow[1]->Draw("");
c3->cd(8);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZPCTow[2]->SetLineColor(kBlue+3);
hZPCTow[2]->SetFillColor(kBlue+3);
hZPCTow[2]->Draw("");
c3->cd(9);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZPCTow[3]->SetLineColor(kBlue+3);
hZPCTow[3]->SetFillColor(kBlue+3);
hZPCTow[3]->Draw("");
c3->cd(10);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZPCTow[4]->SetLineColor(kBlue+3);
hZPCTow[4]->SetFillColor(kBlue+3);
hZPCTow[4]->Draw("");
TCanvas *c32 = new TCanvas("c32","side A ZDCs",700,0,800,400);
c32->Divide(5,2);
c32->cd(1);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZNATow[0]->SetLineColor(kRed);
hZNATow[0]->SetFillColor(kRed);
hZNATow[0]->Draw("");
hZNATow[1]->SetFillColor(kRed);
hZNATow[1]->Draw("");
c32->cd(3);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZNATow[2]->SetLineColor(kRed);
hZNATow[2]->SetFillColor(kRed);
hZNATow[2]->Draw("");
c32->cd(4);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZNATow[3]->SetLineColor(kRed);
hZNATow[3]->SetFillColor(kRed);
hZNATow[3]->Draw("");
c32->cd(5);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZNATow[4]->SetLineColor(kRed);
hZNATow[4]->SetFillColor(kRed);
hZNATow[4]->Draw("");
//
c32->cd(6);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZPATow[0]->SetLineColor(kRed+1);
hZPATow[0]->SetFillColor(kRed+1);
hZPATow[0]->Draw("");
c32->cd(7);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZPATow[1]->SetLineColor(kRed+1);
hZPATow[1]->SetFillColor(kRed+1);
hZPATow[1]->Draw("");
c32->cd(8);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZPATow[2]->SetLineColor(kRed+1);
hZPATow[2]->SetFillColor(kRed+1);
hZPATow[2]->Draw("");
c32->cd(9);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZPATow[3]->SetLineColor(kRed+1);
hZPATow[3]->SetFillColor(kRed+1);
hZPATow[3]->Draw("");
c32->cd(10);
- //gPad->SetLogy(1);
+ gPad->SetLogy(1);
hZPATow[4]->SetLineColor(kRed+1);
hZPATow[4]->SetFillColor(kRed+1);
hZPATow[4]->Draw("");
+
+ //-------------------------------------------------
+ TCanvas *c4 = new TCanvas("c4","TDCs",400,0,700,500);
+ c4->Divide(3,2);
+ c4->cd(1);
+ hTDC[1]->SetFillColor(kPink); hTDC[1]->SetLineColor(kPink);
+ hTDC[1]->Draw("");
+ c4->cd(2);
+ hTDC[3]->SetFillColor(kPink+4); hTDC[3]->SetLineColor(kPink+4);
+ hTDC[3]->Draw("");
+ c4->cd(3);
+ hTDC[5]->SetFillColor(kGreen); hTDC[5]->SetLineColor(kGreen);
+ hTDC[5]->Draw("");
+ c4->cd(4);
+ hTDC[7]->SetFillColor(kGreen+2); hTDC[7]->SetLineColor(kGreen+2);
+ hTDC[7]->Draw("");
+ c4->cd(5);
+ hTDC[8]->SetFillColor(kAzure); hTDC[8]->SetLineColor(kAzure);
+ hTDC[8]->Draw("");
+ c4->cd(6);
+ hTDC[9]->SetFillColor(kAzure+5); hTDC[9]->SetLineColor(kAzure+5);
+ hTDC[9]->Draw("");
}
-
- TFile * fileout = new TFile("ESDhistos.root","recreate");
- fileout->cd();
- //centroidZNsideC->Write();
- //centroidZNsideA->Write();
- enZNC->Write();
- enZNA->Write();
- enZPC->Write();
- enZPA->Write();
- enZEM1->Write();
- enZEM2->Write();
- for(Int_t jj=0; jj<5; jj++){
- hZNCTow[jj]->Write();
- hZPCTow[jj]->Write();
- hZNATow[jj]->Write();
- hZPATow[jj]->Write();
- }
- hSumQZNC->Write();
- hSumQZPC->Write();
- hSumQZNA->Write();
- hSumQZPA->Write();
- //
- hESDword->Write();
- //
- fileout->Close();
}