Skip to content

Commit

Permalink
Bugfix to get all genjets
Browse files Browse the repository at this point in the history
  • Loading branch information
nickmccoll committed Mar 27, 2018
1 parent f205677 commit b2633fa
Show file tree
Hide file tree
Showing 2 changed files with 6 additions and 6 deletions.
6 changes: 3 additions & 3 deletions TreeFillers/src/FatJetFiller.cc
Original file line number Diff line number Diff line change
Expand Up @@ -82,16 +82,16 @@ void FatJetFiller::load(const edm::Event& iEvent, const edm::EventSetup& iSetup)
std::vector<size8> FatJetFiller::processGenJets(){
std::vector<size8> fillGJ(han_genJets->size(),0);
for(const auto& jet : (*han_jets)){
// if(jet.pt() < minJetPT) continue;
if(jet.pt() < minJetPT) continue;
auto genRef = jet.genJetFwdRef().backRef();
if(!genRef.isNull()) fillGJ[genRef.key()] = true;
// if(!genRef.isNull()) fillGJ[genRef.key()] = true;
}
std::vector<size8> genInTreeIDX(han_genJets->size(),255);
size8 curIDX = 0;
const size8 maxIDX = 255;
for(unsigned int iG = 0; iG < han_genJets->size(); ++iG){
const auto& jet = han_genJets->at(iG);
if(jet.pt() < minJetPT && !fillGJ[iG]) continue;
// if(jet.pt() < minJetPT && !fillGJ[iG]) continue;
genInTreeIDX[iG] = std::min(curIDX,maxIDX);
curIDX++;
data.fillMulti(i_gen_pt ,float(jet.pt() ));
Expand Down
6 changes: 3 additions & 3 deletions TreeFillers/src/JetFiller.cc
Original file line number Diff line number Diff line change
Expand Up @@ -109,16 +109,16 @@ bool JetFiller::passTightID(const pat::Jet& jet) {
std::vector<size8> JetFiller::processGenJets(){
std::vector<size8> fillGJ(han_genJets->size(),0);
for(const auto& jet : (*han_jets)){
// if(jet.pt() < minJetPT) continue;
if(jet.pt() < minJetPT) continue;
auto genRef = jet.genJetFwdRef().backRef();
if(!genRef.isNull()) fillGJ[genRef.key()] = true;
// if(!genRef.isNull()) fillGJ[genRef.key()] = true;
}
std::vector<size8> genInTreeIDX(han_genJets->size(),255);
size8 curIDX = 0;
const size8 maxIDX = 255;
for(unsigned int iG = 0; iG < han_genJets->size(); ++iG){
const auto& jet = han_genJets->at(iG);
if(jet.pt() < minJetPT && !fillGJ[iG]) continue;
// if(jet.pt() < minJetPT && !fillGJ[iG]) continue;
genInTreeIDX[iG] = std::min(curIDX,maxIDX);
curIDX++;
data.fillMulti(i_gen_pt ,float(jet.pt() ));
Expand Down

0 comments on commit b2633fa

Please sign in to comment.