- if (fDoEoverP) {
- fHistEoverP = new TH2F("fHistEoverP","E/P vs. E", fNbins, fMinPt, fMaxPt, fNbins, 0, 10);
- fHistEoverP->GetXaxis()->SetTitle("E [GeV]");
- fHistEoverP->GetYaxis()->SetTitle("E/P [c]");
- fOutput->Add(fHistEoverP);
+ if (fAnaType == kEMCAL) {
+ if (fDoEoverP) {
+ fHistEoverP = new TH2F("fHistEoverP","E/P vs. E", fNbins, fMinBinPt, fMaxBinPt, fNbins, 0, 10);
+ fHistEoverP->GetXaxis()->SetTitle("E [GeV]");
+ fHistEoverP->GetYaxis()->SetTitle("E/P [c]");
+ fOutput->Add(fHistEoverP);
+ }
+
+ fHistCellsEnergy = new TH1F("fHistCellsEnergy","Energy spectrum of cells", fNbins, fMinBinPt, fMaxBinPt);
+ fHistCellsEnergy->GetXaxis()->SetTitle("E [GeV]");
+ fHistCellsEnergy->GetYaxis()->SetTitle("counts");
+ fOutput->Add(fHistCellsEnergy);
+
+ fHistChVSneCells = new TH2F("fHistChVSneCells","Charged energy vs. neutral (cells) energy", fNbins, fMinBinPt * 4, fMaxBinPt * 4, fNbins, fMinBinPt * 4, fMaxBinPt * 4);
+ fHistChVSneCells->GetXaxis()->SetTitle("E [GeV]");
+ fHistChVSneCells->GetYaxis()->SetTitle("P [GeV/c]");
+ fOutput->Add(fHistChVSneCells);
+
+ fHistChVSneClus = new TH2F("fHistChVSneClus","Charged energy vs. neutral (clusters) energy", fNbins, fMinBinPt * 4, fMaxBinPt * 4, fNbins, fMinBinPt * 4, fMaxBinPt * 4);
+ fHistChVSneClus->GetXaxis()->SetTitle("E [GeV]");
+ fHistChVSneClus->GetYaxis()->SetTitle("P [GeV/c]");
+ fOutput->Add(fHistChVSneClus);
+
+ fHistChVSneCorrCells = new TH2F("fHistChVSneCorrCells","Charged energy vs. neutral (corrected cells) energy", fNbins, fMinBinPt * 4, fMaxBinPt * 4, fNbins, fMinBinPt * 4, fMaxBinPt * 4);
+ fHistChVSneCorrCells->GetXaxis()->SetTitle("E [GeV]");
+ fHistChVSneCorrCells->GetYaxis()->SetTitle("P [GeV/c]");
+ fOutput->Add(fHistChVSneCorrCells);