jet charge is processed correctly

unelg · Aug 9, 2022 · f5c246b · f5c246b
1 parent 5c4adfe
commit f5c246b
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 1 deletion.
diff --git 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: "<<jet->PT<<", eta: "<<jet->Eta<<", phi: "<<jet->Phi <<" T:"<<jet->T<<endl;
+//      cout<<"This Jet pt: "<<jet->PT<<", eta: "<<jet->Eta<<", phi: "<<jet->Phi <<" Q:"<<jet->Charge<<endl;
       adbxj= new dbxJet(alv);
       adbxj->setCharge(jet->Charge);
       adbxj->setParticleIndx(i);

diff --git a/analysis_core/FuncNode.cpp b/analysis_core/FuncNode.cpp
@@ -154,6 +154,7 @@ void FuncNode::partConstruct(AnalysisObjects *ao, std::vector<myParticle*> *inpu
                                             case jet_t: DEBUG("jet:"<<ai<<" \t");
                                                         inputPart->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() );