From c0ec1c443ed886d7a73bf52fd40ae7cfdcfc00d9 Mon Sep 17 00:00:00 2001 From: kleinb Date: Fri, 14 Jan 2011 16:51:11 +0000 Subject: [PATCH] added protection against missing 4 vector area --- .../AliAnalysisTaskJetBackgroundSubtract.cxx | 31 ++++++++++--------- 1 file changed, 16 insertions(+), 15 deletions(-) diff --git a/JETAN/AliAnalysisTaskJetBackgroundSubtract.cxx b/JETAN/AliAnalysisTaskJetBackgroundSubtract.cxx index ed0868d3960..b4b32bb7d33 100644 --- a/JETAN/AliAnalysisTaskJetBackgroundSubtract.cxx +++ b/JETAN/AliAnalysisTaskJetBackgroundSubtract.cxx @@ -453,22 +453,23 @@ void AliAnalysisTaskJetBackgroundSubtract::UserExec(Option_t */*option*/) }//kRhoRC - else if(fSubtraction==k4Area){ - - backgroundv.SetPxPyPzE(rho*(jet->VectorAreaCharged())->Px(),rho*(jet->VectorAreaCharged())->Py(),rho*(jet->VectorAreaCharged())->Pz(),rho*(jet->VectorAreaCharged())->E()); - if((backgroundv.E()>jet->E())&&(backgroundv.Pt()>jet->Pt())){ - - // optionally rescale it and keep?? - bAdd = RescaleJetMomentum(&tmpNewJet,0.1); - if(h2PtInOut)h2PtInOut->Fill(jet->Pt(),0.1); - } - else{ - bAdd = RescaleJet4vector(&tmpNewJet,backgroundv); - } - // add background estimates to the new jet object - // allows to recover old p_T and rho... - tmpNewJet.SetBgEnergy(backgroundv.P(),0); + else if(fSubtraction==k4Area&&jet->VectorAreaCharged()){ + + backgroundv.SetPxPyPzE(rho*(jet->VectorAreaCharged())->Px(),rho*(jet->VectorAreaCharged())->Py(),rho*(jet->VectorAreaCharged())->Pz(),rho*(jet->VectorAreaCharged())->E()); + if((backgroundv.E()>jet->E())&&(backgroundv.Pt()>jet->Pt())){ + + // optionally rescale it and keep?? + bAdd = RescaleJetMomentum(&tmpNewJet,0.1); + if(h2PtInOut)h2PtInOut->Fill(jet->Pt(),0.1); + } + else{ + bAdd = RescaleJet4vector(&tmpNewJet,backgroundv); + } + // add background estimates to the new jet object + // allows to recover old p_T and rho... + tmpNewJet.SetBgEnergy(backgroundv.P(),0); + }//kArea4vector -- 2.43.0