+ fHistJets2ZvsPt = new TH2F("fHistJets2ZvsPt", "fHistJets2ZvsPt", 120, 0, 1.2, fNbins, fMinBinPt, fMaxBinPt);
+ fHistJets2ZvsPt->GetXaxis()->SetTitle("Z");
+ fHistJets2ZvsPt->GetYaxis()->SetTitle("p_{T,2} (GeV/c)");
+ fHistJets2ZvsPt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistJets2ZvsPt);
+
+ fHistJets2NEFvsPt = new TH2F("fHistJets2NEFvsPt", "fHistJets2NEFvsPt", 120, 0, 1.2, fNbins, fMinBinPt, fMaxBinPt);
+ fHistJets2NEFvsPt->GetXaxis()->SetTitle("NEF");
+ fHistJets2NEFvsPt->GetYaxis()->SetTitle("p_{T,2} (GeV/c)");
+ fHistJets2NEFvsPt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistJets2NEFvsPt);
+
+ fHistJets2CEFvsCEFPt = new TH2F("fHistJets2CEFvsCEFPt", "fHistJets2CEFvsCEFPt", 120, 0, 1.2, fNbins, fMinBinPt, fMaxBinPt);
+ fHistJets2CEFvsCEFPt->GetXaxis()->SetTitle("1-NEF");
+ fHistJets2CEFvsCEFPt->GetYaxis()->SetTitle("(1-NEF)*p_{T,2} (GeV/c)");
+ fHistJets2CEFvsCEFPt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistJets2CEFvsCEFPt);
+
+ // Matching histograms
+
+ fHistCommonEnergy1vsJet1Pt = new TH2F("fHistCommonEnergy1vsJet1Pt", "fHistCommonEnergy1vsJet1Pt", fNbins/2, 0, 1.2, fNbins/2, fMinBinPt, fMaxBinPt);
+ fHistCommonEnergy1vsJet1Pt->GetXaxis()->SetTitle("Common energy 1 (%)");
+ fHistCommonEnergy1vsJet1Pt->GetYaxis()->SetTitle("p_{T,1}");
+ fHistCommonEnergy1vsJet1Pt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistCommonEnergy1vsJet1Pt);
+
+ fHistCommonEnergy2vsJet2Pt = new TH2F("fHistCommonEnergy2vsJet2Pt", "fHistCommonEnergy2vsJet2Pt", fNbins/2, 0, 1.2, fNbins/2, fMinBinPt, fMaxBinPt);
+ fHistCommonEnergy2vsJet2Pt->GetXaxis()->SetTitle("Common energy 2 (%)");
+ fHistCommonEnergy2vsJet2Pt->GetYaxis()->SetTitle("p_{T,2}");
+ fHistCommonEnergy2vsJet2Pt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistCommonEnergy2vsJet2Pt);
+
+ fHistDistancevsJet1Pt = new TH2F("fHistDistancevsJet1Pt", "fHistDistancevsJet1Pt", fNbins/2, 0, 1.2, fNbins/2, fMinBinPt, fMaxBinPt);
+ fHistDistancevsJet1Pt->GetXaxis()->SetTitle("Distance");
+ fHistDistancevsJet1Pt->GetYaxis()->SetTitle("p_{T,1}");
+ fHistDistancevsJet1Pt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDistancevsJet1Pt);
+
+ fHistDistancevsJet2Pt = new TH2F("fHistDistancevsJet2Pt", "fHistDistancevsJet2Pt", fNbins/2, 0, 1.2, fNbins/2, fMinBinPt, fMaxBinPt);
+ fHistDistancevsJet2Pt->GetXaxis()->SetTitle("Distance");
+ fHistDistancevsJet2Pt->GetYaxis()->SetTitle("p_{T,2}");
+ fHistDistancevsJet2Pt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDistancevsJet2Pt);
+
+ fHistDistancevsCommonEnergy1 = new TH2F("fHistDistancevsCommonEnergy1", "fHistDistancevsCommonEnergy1", fNbins/2, 0, 1.2, fNbins/2, 0, 1.2);
+ fHistDistancevsCommonEnergy1->GetXaxis()->SetTitle("Distance");
+ fHistDistancevsCommonEnergy1->GetYaxis()->SetTitle("Common energy 1 (%)");
+ fHistDistancevsCommonEnergy1->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDistancevsCommonEnergy1);
+
+ fHistDistancevsCommonEnergy2 = new TH2F("fHistDistancevsCommonEnergy2", "fHistDistancevsCommonEnergy2", fNbins/2, 0, 1.2, fNbins/2, 0, 1.2);
+ fHistDistancevsCommonEnergy2->GetXaxis()->SetTitle("Distance");
+ fHistDistancevsCommonEnergy2->GetYaxis()->SetTitle("Common energy 2 (%)");
+ fHistDistancevsCommonEnergy2->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDistancevsCommonEnergy2);
+
+ fHistCommonEnergy1vsCommonEnergy2 = new TH2F("fHistCommonEnergy1vsCommonEnergy2", "fHistCommonEnergy1vsCommonEnergy2", fNbins/2, 0, 1.2, fNbins/2, 0, 1.2);
+ fHistCommonEnergy1vsCommonEnergy2->GetXaxis()->SetTitle("Common energy 1 (%)");
+ fHistCommonEnergy1vsCommonEnergy2->GetYaxis()->SetTitle("Common energy 2 (%)");
+ fHistCommonEnergy1vsCommonEnergy2->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistCommonEnergy1vsCommonEnergy2);
+
+ fHistDeltaEtaDeltaPhi = new TH2F("fHistDeltaEtaDeltaPhi", "fHistDeltaEtaDeltaPhi", fNbins/4, -1, 1, fNbins/4, -TMath::Pi()/2, TMath::Pi()*3/2);
+ fHistDeltaEtaDeltaPhi->GetXaxis()->SetTitle("Common energy 1 (%)");
+ fHistDeltaEtaDeltaPhi->GetYaxis()->SetTitle("Common energy 2 (%)");
+ fHistDeltaEtaDeltaPhi->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaEtaDeltaPhi);
+
+ fHistJet2PtOverJet1PtvsJet2Pt = new TH2F("fHistJet2PtOverJet1PtvsJet2Pt", "fHistJet2PtOverJet1PtvsJet2Pt", fNbins, fMinBinPt, fMaxBinPt, 300, 0, 1.5);
+ fHistJet2PtOverJet1PtvsJet2Pt->GetXaxis()->SetTitle("p_{T,2}");
+ fHistJet2PtOverJet1PtvsJet2Pt->GetYaxis()->SetTitle("p_{T,2} / p_{T,1}");
+ fHistJet2PtOverJet1PtvsJet2Pt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistJet2PtOverJet1PtvsJet2Pt);
+
+ fHistJet1PtOverJet2PtvsJet1Pt = new TH2F("fHistJet1PtOverJet2PtvsJet1Pt", "fHistJet1PtOverJet2PtvsJet1Pt", fNbins, fMinBinPt, fMaxBinPt, 300, 0, 1.5);
+ fHistJet1PtOverJet2PtvsJet1Pt->GetXaxis()->SetTitle("p_{T,1}");
+ fHistJet1PtOverJet2PtvsJet1Pt->GetYaxis()->SetTitle("p_{T,1} / p_{T,2}");
+ fHistJet1PtOverJet2PtvsJet1Pt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistJet1PtOverJet2PtvsJet1Pt);
+
+ fHistDeltaPtvsJet1Pt = new TH2F("fHistDeltaPtvsJet1Pt", "fHistDeltaPtvsJet1Pt",
+ fNbins, fMinBinPt, fMaxBinPt, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaPtvsJet1Pt->GetXaxis()->SetTitle("p_{T,1}");
+ fHistDeltaPtvsJet1Pt->GetYaxis()->SetTitle("#deltap_{T} (GeV/c)");
+ fHistDeltaPtvsJet1Pt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaPtvsJet1Pt);
+
+ fHistDeltaPtvsJet2Pt = new TH2F("fHistDeltaPtvsJet2Pt", "fHistDeltaPtvsJet2Pt",
+ fNbins, fMinBinPt, fMaxBinPt, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaPtvsJet2Pt->GetXaxis()->SetTitle("p_{T,2}");
+ fHistDeltaPtvsJet2Pt->GetYaxis()->SetTitle("#deltap_{T} (GeV/c)");
+ fHistDeltaPtvsJet2Pt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaPtvsJet2Pt);
+
+ fHistDeltaPtOverJet1PtvsJet1Pt = new TH2F("fHistDeltaPtOverJet1PtvsJet1Pt", "fHistDeltaPtOverJet1PtvsJet1Pt",
+ fNbins, fMinBinPt, fMaxBinPt, fNbins, -5, 5);
+ fHistDeltaPtOverJet1PtvsJet1Pt->GetXaxis()->SetTitle("p_{T,1}");
+ fHistDeltaPtOverJet1PtvsJet1Pt->GetYaxis()->SetTitle("#deltap_{T} / p_{T,1}");
+ fHistDeltaPtOverJet1PtvsJet1Pt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaPtOverJet1PtvsJet1Pt);
+
+ fHistDeltaPtOverJet2PtvsJet2Pt = new TH2F("fHistDeltaPtOverJet2PtvsJet2Pt", "fHistDeltaPtOverJet2PtvsJet2Pt",
+ fNbins, fMinBinPt, fMaxBinPt, fNbins, -5, 5);
+ fHistDeltaPtOverJet2PtvsJet2Pt->GetXaxis()->SetTitle("p_{T,2}");
+ fHistDeltaPtOverJet2PtvsJet2Pt->GetYaxis()->SetTitle("#deltap_{T} / p_{T,2}");
+ fHistDeltaPtOverJet2PtvsJet2Pt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaPtOverJet2PtvsJet2Pt);
+
+ fHistDeltaPtvsDistance = new TH2F("fHistDeltaPtvsDistance", "fHistDeltaPtvsDistance",
+ fNbins/2, 0, 1.2, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaPtvsDistance->GetXaxis()->SetTitle("Distance");
+ fHistDeltaPtvsDistance->GetYaxis()->SetTitle("#deltap_{T} (GeV/c)");
+ fHistDeltaPtvsDistance->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaPtvsDistance);
+
+ fHistDeltaPtvsCommonEnergy1 = new TH2F("fHistDeltaPtvsCommonEnergy1", "fHistDeltaPtvsCommonEnergy1",
+ fNbins/2, 0, 1.2, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaPtvsCommonEnergy1->GetXaxis()->SetTitle("Common energy 1 (%)");
+ fHistDeltaPtvsCommonEnergy1->GetYaxis()->SetTitle("#deltap_{T} (GeV/c)");
+ fHistDeltaPtvsCommonEnergy1->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaPtvsCommonEnergy1);
+
+ fHistDeltaPtvsCommonEnergy2 = new TH2F("fHistDeltaPtvsCommonEnergy2", "fHistDeltaPtvsCommonEnergy2",
+ fNbins/2, 0, 1.2, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaPtvsCommonEnergy2->GetXaxis()->SetTitle("Common energy 2 (%)");
+ fHistDeltaPtvsCommonEnergy2->GetYaxis()->SetTitle("#deltap_{T} (GeV/c)");
+ fHistDeltaPtvsCommonEnergy2->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaPtvsCommonEnergy2);
+
+ fHistDeltaPtvsArea1 = new TH2F("fHistDeltaPtvsArea1", "fHistDeltaPtvsArea1",
+ fNbins/2, 0, 1.5, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaPtvsArea1->GetXaxis()->SetTitle("A_{jet,1}");
+ fHistDeltaPtvsArea1->GetYaxis()->SetTitle("#deltap_{T} (GeV/c)");
+ fHistDeltaPtvsArea1->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaPtvsArea1);
+
+ fHistDeltaPtvsArea2 = new TH2F("fHistDeltaPtvsArea2", "fHistDeltaPtvsArea2",
+ fNbins/2, 0, 1.5, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaPtvsArea2->GetXaxis()->SetTitle("A_{jet,2}");
+ fHistDeltaPtvsArea2->GetYaxis()->SetTitle("#deltap_{T} (GeV/c)");
+ fHistDeltaPtvsArea2->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaPtvsArea2);
+
+ fHistDeltaPtvsDeltaArea = new TH2F("fHistDeltaPtvsDeltaArea", "fHistDeltaPtvsDeltaArea",
+ fNbins, -1.+1./fNbins, 1.+1./fNbins, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaPtvsDeltaArea->GetXaxis()->SetTitle("#deltaA_{jet}");
+ fHistDeltaPtvsDeltaArea->GetYaxis()->SetTitle("#deltap_{T} (GeV/c)");
+ fHistDeltaPtvsDeltaArea->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaPtvsDeltaArea);
+
+ fHistJet1PtvsJet2Pt = new TH2F("fHistJet1PtvsJet2Pt", "fHistJet1PtvsJet2Pt", fNbins, fMinBinPt, fMaxBinPt, fNbins, fMinBinPt, fMaxBinPt);
+ fHistJet1PtvsJet2Pt->GetXaxis()->SetTitle("p_{T,1}");
+ fHistJet1PtvsJet2Pt->GetYaxis()->SetTitle("p_{T,2}");
+ fHistJet1PtvsJet2Pt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistJet1PtvsJet2Pt);
+
+ if (fIsJet1Rho || fIsJet2Rho) {
+ if (!fIsJet1Rho)
+ 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);
+
+ if (!fIsJet2Rho)
+ 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");
+ fOutput->Add(fHistDeltaCorrPtvsJet2CorrPt);
+
+ fHistDeltaCorrPtOverJet1CorrPtvsJet1CorrPt = new TH2F("fHistDeltaCorrPtOverJet1CorrPtvsJet1CorrPt", "fHistDeltaCorrPtOverJet1CorrPtvsJet1CorrPt",
+ 2*fNbins, -fMaxBinPt, fMaxBinPt, fNbins, -5, 5);
+ fHistDeltaCorrPtOverJet1CorrPtvsJet1CorrPt->GetXaxis()->SetTitle("p_{T,1}^{corr}");
+ fHistDeltaCorrPtOverJet1CorrPtvsJet1CorrPt->GetYaxis()->SetTitle("#deltap_{T}^{corr} / p_{T,1}^{corr}");
+ fHistDeltaCorrPtOverJet1CorrPtvsJet1CorrPt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaCorrPtOverJet1CorrPtvsJet1CorrPt);
+
+ fHistDeltaCorrPtOverJet2CorrPtvsJet2CorrPt = new TH2F("fHistDeltaCorrPtOverJet2CorrPtvsJet2CorrPt", "fHistDeltaCorrPtOverJet2CorrPtvsJet2CorrPt",
+ 2*fNbins, -fMaxBinPt, fMaxBinPt, fNbins, -5, 5);
+ fHistDeltaCorrPtOverJet2CorrPtvsJet2CorrPt->GetXaxis()->SetTitle("p_{T,2}^{corr}");
+ fHistDeltaCorrPtOverJet2CorrPtvsJet2CorrPt->GetYaxis()->SetTitle("#deltap_{T}^{corr} / p_{T,2}^{corr}");
+ fHistDeltaCorrPtOverJet2CorrPtvsJet2CorrPt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaCorrPtOverJet2CorrPtvsJet2CorrPt);
+
+ fHistDeltaCorrPtvsDistance = new TH2F("fHistDeltaCorrPtvsDistance", "fHistDeltaCorrPtvsDistance",
+ fNbins/2, 0, 1.2, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaCorrPtvsDistance->GetXaxis()->SetTitle("Distance");
+ fHistDeltaCorrPtvsDistance->GetYaxis()->SetTitle("#deltap_{T}^{corr} (GeV/c)");
+ fHistDeltaCorrPtvsDistance->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaCorrPtvsDistance);
+
+ fHistDeltaCorrPtvsCommonEnergy1 = new TH2F("fHistDeltaCorrPtvsCommonEnergy1", "fHistDeltaCorrPtvsCommonEnergy1",
+ fNbins/2, 0, 1.2, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaCorrPtvsCommonEnergy1->GetXaxis()->SetTitle("Common energy 1 (%)");
+ fHistDeltaCorrPtvsCommonEnergy1->GetYaxis()->SetTitle("#deltap_{T}^{corr} (GeV/c)");
+ fHistDeltaCorrPtvsCommonEnergy1->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaCorrPtvsCommonEnergy1);
+
+ fHistDeltaCorrPtvsCommonEnergy2 = new TH2F("fHistDeltaCorrPtvsCommonEnergy2", "fHistDeltaCorrPtvsCommonEnergy2",
+ fNbins/2, 0, 1.2, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaCorrPtvsCommonEnergy2->GetXaxis()->SetTitle("Common energy 2 (%)");
+ fHistDeltaCorrPtvsCommonEnergy2->GetYaxis()->SetTitle("#deltap_{T}^{corr} (GeV/c)");
+ fHistDeltaCorrPtvsCommonEnergy2->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaCorrPtvsCommonEnergy2);
+
+ fHistDeltaCorrPtvsArea1 = new TH2F("fHistDeltaCorrPtvsArea1", "fHistDeltaCorrPtvsArea1",
+ fNbins/2, 0, 1.5, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaCorrPtvsArea1->GetXaxis()->SetTitle("A_{jet,1}");
+ fHistDeltaCorrPtvsArea1->GetYaxis()->SetTitle("#deltap_{T}^{corr} (GeV/c)");
+ fHistDeltaCorrPtvsArea1->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaCorrPtvsArea1);
+
+ fHistDeltaCorrPtvsArea2 = new TH2F("fHistDeltaCorrPtvsArea2", "fHistDeltaCorrPtvsArea2",
+ fNbins/2, 0, 1.5, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaCorrPtvsArea2->GetXaxis()->SetTitle("A_{jet,2}");
+ fHistDeltaCorrPtvsArea2->GetYaxis()->SetTitle("#deltap_{T}^{corr} (GeV/c)");
+ fHistDeltaCorrPtvsArea2->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaCorrPtvsArea2);
+
+ fHistDeltaCorrPtvsDeltaArea = new TH2F("fHistDeltaCorrPtvsDeltaArea", "fHistDeltaCorrPtvsDeltaArea",
+ fNbins, -1.+1./fNbins, 1.+1./fNbins, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistDeltaCorrPtvsDeltaArea->GetXaxis()->SetTitle("#deltaA_{jet}");
+ fHistDeltaCorrPtvsDeltaArea->GetYaxis()->SetTitle("#deltap_{T}^{corr} (GeV/c)");
+ fHistDeltaCorrPtvsDeltaArea->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistDeltaCorrPtvsDeltaArea);
+
+ if (!fIsJet1Rho)
+ fHistJet1CorrPtvsJet2CorrPt = new TH2F("fHistJet1CorrPtvsJet2CorrPt", "fHistJet1CorrPtvsJet2CorrPt",
+ fNbins, fMinBinPt, fMaxBinPt, 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ else if (!fIsJet2Rho)
+ fHistJet1CorrPtvsJet2CorrPt = new TH2F("fHistJet1CorrPtvsJet2CorrPt", "fHistJet1CorrPtvsJet2CorrPt",
+ 2*fNbins, -fMaxBinPt, fMaxBinPt, fNbins, fMinBinPt, fMaxBinPt);
+ else
+ fHistJet1CorrPtvsJet2CorrPt = new TH2F("fHistJet1CorrPtvsJet2CorrPt", "fHistJet1CorrPtvsJet2CorrPt",
+ 2*fNbins, -fMaxBinPt, fMaxBinPt,
+ 2*fNbins, -fMaxBinPt, fMaxBinPt);
+ fHistJet1CorrPtvsJet2CorrPt->GetXaxis()->SetTitle("p_{T,1}^{corr}");
+ fHistJet1CorrPtvsJet2CorrPt->GetYaxis()->SetTitle("p_{T,2}^{corr}");
+ fHistJet1CorrPtvsJet2CorrPt->GetZaxis()->SetTitle("counts");
+ fOutput->Add(fHistJet1CorrPtvsJet2CorrPt);
+ }