From: loizides Date: Mon, 29 Apr 2013 14:04:42 +0000 (+0000) Subject: remove unnecessary bins from histograms depending on given parameters (salvatore) X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=5c045cdec490d9cef9d9bb947e044c20bba2c33a;p=u%2Fmrichter%2FAliRoot.git remove unnecessary bins from histograms depending on given parameters (salvatore) --- diff --git a/PWGJE/EMCALJetTasks/AliJetResponseMaker.cxx b/PWGJE/EMCALJetTasks/AliJetResponseMaker.cxx index f6dbe01b5c5..2ce1f420123 100644 --- a/PWGJE/EMCALJetTasks/AliJetResponseMaker.cxx +++ b/PWGJE/EMCALJetTasks/AliJetResponseMaker.cxx @@ -770,16 +770,26 @@ void AliJetResponseMaker::UserCreateOutputObjects() fHistJet1PtvsJet2Pt->GetZaxis()->SetTitle("counts"); fOutput->Add(fHistJet1PtvsJet2Pt); - if (!fRhoName.IsNull() || !fRho2Name.IsNull()) { - fHistDeltaCorrPtvsJet1CorrPt = new TH2F("fHistDeltaCorrPtvsJet1CorrPt", "fHistDeltaCorrPtvsJet1CorrPt", - 2*fNbins, -fMaxBinPt, fMaxBinPt, 2*fNbins, -fMaxBinPt, fMaxBinPt); + if (!fRhoName.IsNull() || !fRho2Name.IsNull()) { + if (fRhoName.IsNull()) + fHistDeltaCorrPtvsJet1CorrPt = new TH2F("fHistDeltaCorrPtvsJet1CorrPt", "fHistDeltaCorrPtvsJet1CorrPt", + fNbins, fMinBinPt, fMaxBinPt, 2*fNbins, -fMaxBinPt, fMaxBinPt); + else + fHistDeltaCorrPtvsJet1CorrPt = new TH2F("fHistDeltaCorrPtvsJet1CorrPt", "fHistDeltaCorrPtvsJet1CorrPt", + 2*fNbins, -fMaxBinPt, fMaxBinPt, 2*fNbins, -fMaxBinPt, fMaxBinPt); + fHistDeltaCorrPtvsJet1CorrPt->GetXaxis()->SetTitle("p_{T,1}^{corr}"); fHistDeltaCorrPtvsJet1CorrPt->GetYaxis()->SetTitle("#deltap_{T}^{corr} (GeV/c)"); fHistDeltaCorrPtvsJet1CorrPt->GetZaxis()->SetTitle("counts"); fOutput->Add(fHistDeltaCorrPtvsJet1CorrPt); - fHistDeltaCorrPtvsJet2CorrPt = new TH2F("fHistDeltaCorrPtvsJet2CorrPt", "fHistDeltaCorrPtvsJet2CorrPt", - 2*fNbins, -fMaxBinPt, fMaxBinPt, 2*fNbins, -fMaxBinPt, fMaxBinPt); + if (fRho2Name.IsNull()) + fHistDeltaCorrPtvsJet2CorrPt = new TH2F("fHistDeltaCorrPtvsJet2CorrPt", "fHistDeltaCorrPtvsJet2CorrPt", + fNbins, fMinBinPt, fMaxBinPt, 2*fNbins, -fMaxBinPt, fMaxBinPt); + else + fHistDeltaCorrPtvsJet2CorrPt = new TH2F("fHistDeltaCorrPtvsJet2CorrPt", "fHistDeltaCorrPtvsJet2CorrPt", + 2*fNbins, -fMaxBinPt, fMaxBinPt, 2*fNbins, -fMaxBinPt, fMaxBinPt); + fHistDeltaCorrPtvsJet2CorrPt->GetXaxis()->SetTitle("p_{T,2}^{corr}"); fHistDeltaCorrPtvsJet2CorrPt->GetYaxis()->SetTitle("#deltap_{T}^{corr} (GeV/c)"); fHistDeltaCorrPtvsJet2CorrPt->GetZaxis()->SetTitle("counts");