From f5c246b3699412d81b0db2f134b6267f305b8d28 Mon Sep 17 00:00:00 2001 From: unelg Date: Tue, 9 Aug 2022 10:11:40 +0200 Subject: [PATCH] jet charge is processed correctly --- CLA/delphes.C | 2 +- analysis_core/FuncNode.cpp | 1 + 2 files changed, 2 insertions(+), 1 deletion(-) diff --git a/CLA/delphes.C b/CLA/delphes.C index 90800842..24148ad5 100644 --- a/CLA/delphes.C +++ b/CLA/delphes.C @@ -234,7 +234,7 @@ void delphes::Loop(analy_struct aselect, char *extname) for(i = 0; i < branchJet->GetEntriesFast(); ++i) { jet = (Jet*) branchJet->At(i); alv.SetPtEtaPhiM( jet->PT, jet->Eta, jet->Phi, jet->Mass ); // all in GeV -// cout<<"This Jet pt: "<PT<<", eta: "<Eta<<", phi: "<Phi <<" T:"<T<PT<<", eta: "<Eta<<", phi: "<Phi <<" Q:"<Charge<setCharge(jet->Charge); adbxj->setParticleIndx(i); diff --git a/analysis_core/FuncNode.cpp b/analysis_core/FuncNode.cpp index 8f8c1a3f..db7e3d9b 100644 --- a/analysis_core/FuncNode.cpp +++ b/analysis_core/FuncNode.cpp @@ -154,6 +154,7 @@ void FuncNode::partConstruct(AnalysisObjects *ao, std::vector *inpu case jet_t: DEBUG("jet:"<setTlv(inputPart->lv()+sgn*ao->jets[ac].at(ai).lv() ); // any jet inputPart->setFlavor(inputPart->Flavor() +ao->jets[ac].at(ai).Flavor() ); + inputPart->setCharge(inputPart->q()+ao->jets[ac].at(ai).q() ); // inputPart->setIsTight( inputPart->isTight() // add to the existing one // + ao->jets[ac].at(ai).isbtagged_77() +100* ao->jets[ac].at(ai).isTautagged() ); inputPart->setIsTight (ao->jets[ac].at(ai).isTight() );