Merge pull request #19702 from intrepid42/ParticleLevelProducer_80X

Bugfix and c tagging for ParticleLevelProducer (80X)
cms-sw · Jul 19, 2017 · c735270 · c735270
2 parents 92b4c15 + 2114e74
commit c735270
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 1 deletion.
diff --git a/GeneratorInterface/RivetInterface/interface/RivetAnalysis.h b/GeneratorInterface/RivetInterface/interface/RivetAnalysis.h
@@ -156,7 +156,7 @@ namespace Rivet {
         _fatjets   = applyProjection<FastJets>(event, "FatJets").jetsByPt(fatjet_cut);
         _photons   = applyProjection<FinalState>(event, "Photons").particlesByPt();
         _neutrinos = applyProjection<FinalState>(event, "Neutrinos").particlesByPt();
-        _met       = -applyProjection<MissingMomentum>(event, "MET").vectorEt();
+        _met       = applyProjection<MissingMomentum>(event, "MET").missingMomentum().p3();
       };
 
       // Do nothing here

diff --git a/GeneratorInterface/RivetInterface/plugins/ParticleLevelProducer.cc b/GeneratorInterface/RivetInterface/plugins/ParticleLevelProducer.cc
@@ -46,6 +46,7 @@ void ParticleLevelProducer::addGenJet(Rivet::Jet jet, std::unique_ptr<reco::GenJ
   genJet.setP4(p4(jet));
   genJet.setVertex(genVertex_);
   if ( jet.bTagged() ) genJet.setPdgId(5);
+  else if ( jet.cTagged() ) genJet.setPdgId(4);
   genJet.setJetArea(pjet.has_area() ? pjet.area() : 0);
 
   for ( auto const & p : jet.particles()) {