Skip to content

Commit

Permalink
Fixed behavior in QjetsAdder when reclustering jet does not return or…
Browse files Browse the repository at this point in the history
…iginal jet
  • Loading branch information
sidnarayanan committed Aug 5, 2015
1 parent 6e97fa6 commit 1260a0a
Showing 1 changed file with 1 addition and 2 deletions.
3 changes: 1 addition & 2 deletions RecoJets/JetProducers/plugins/QjetsAdder.cc
Expand Up @@ -40,11 +40,10 @@ void QjetsAdder::produce(edm::Event & iEvent, const edm::EventSetup & iSetup) {
fastjet::ClusterSequence thisClustering_basic(allconstits, jetDef);
std::vector<fastjet::PseudoJet> out_jets_basic = sorted_by_pt(thisClustering_basic.inclusive_jets(cutoff_));
//std::cout << newCand.pt() << " " << out_jets_basic.size() <<std::endl;
if(out_jets_basic.size()==0){ // jet reclustering failed, most likely due to the higher cutoff. Use a recognizeable default value for this jet
if(out_jets_basic.size()!=1){ // jet reclustering did not return exactly 1 jet, most likely due to the higher cutoff or large cone size. Use a recognizeable default value for this jet
QjetsVolatility.push_back(-1);
continue;
}
assert( out_jets_basic.size()==1 ); // jet reclustering of one jet should yield exactly one jet at this stage

// setup objects for qjets computation
fastjet::JetDefinition qjet_def(&qjetsAlgo_);
Expand Down

0 comments on commit 1260a0a

Please sign in to comment.