GenParticles2HepMC for MiniAOD #18115
Conversation
|
A new Pull Request was created by @jhgoh (Junghwan John Goh) for master. It involves the following packages: GeneratorInterface/RivetInterface @cmsbuild, @govoni, @perrozzi, @thuer, @davidlange6 can you please review it and eventually sign? Thanks. cms-bot commands are listed here #13028 |
|
please test |
|
The tests are being triggered in jenkins. |
|
(to be fair, most of the credit goes to David Sperka @dsperka and Andrew Gilbert) |
|
Comparison job queued. |
|
Main body is modified to include main contributors. |
| @@ -3,4 +3,6 @@ | |||
| genParticles2HepMC = cms.EDProducer("GenParticles2HepMCConverter", | |||
| genParticles = cms.InputTag("genParticles"), | |||
There was a problem hiding this comment.
would it make sense to add another commented line like
genParticles = cms.InputTag("mergedGenParticles"), # in case mergedGenParticles are created
or add an example in the test directory?
There was a problem hiding this comment.
Done, added commented line just below of line 4, also added a cfg file under the test directory.
cmsbuild
added
comparison-pending
and removed
comparison-available
tests-approved
labels
cmsbuild
added
comparison-pending
generators-pending
pending-signatures
tests-pending
and removed
generators-approved
tests-approved
labels
|
@davidlange6 simplified code by using simple vector and basic std::find. |
|
please test |
|
The tests are being triggered in jenkins. |
|
Comparison job queued. |
|
+1 |
cmsbuild
added
fully-signed
generators-approved
and removed
generators-pending
pending-signatures
labels
|
This pull request is fully signed and it will be integrated in one of the next master IBs (tests are also fine). This pull request requires discussion in the ORP meeting before it's merged. @Muzaffar, @davidlange6, @smuzaffar |
|
Comparison is ready Comparison Summary:
|
|
+1 |
This PR enables to reproduce HepMC from the prunedGenParticles+packedGenParticles in the MiniAOD, to be used in the Rivet analyses.
For the RECO/AOD's, the workflow remains same: genParticles -> genParticle2HepMC ->rivetAnalyzer
and for the MiniAOD users: prunedGenParticles+packedGenParticles -> mergedGenParticles -> genParticles2HepMC -> RivetAnalyzer
The MergedGenParticleProducer is developed by @dsperka, @ajgilbert, @perrozzi for the Higgs group, with modification to the GenParticles2HepMCConverter to assign the "signal vertex".
(https://indico.cern.ch/event/596884/contributions/2436141/attachments/1396171/2128423/HTXS_GENmeeting_16Jan.pdf#page=10 https://twiki.cern.ch/twiki/bin/view/CMS/HiggsTemplateCrossSection)
The signal vertex assignment can be done by adding PDG id of the interested particle in the configuration, e.g.,
process.genParticles2HepMC.signalParticlePdgIds = cms.vint32(25). Default value assigns the first one to be the signal vtx.The suggested codes are modified by @intrepid42 - inserts dummy pions to make Rivet to distinguish particles from hadron decay which are dropped during the pruning. This trick is necessary to construct pseudo-top. (https://indico.cern.ch/event/607918/contributions/2450369/attachments/1405142/2146446/Elvire_170201.pdf#page=5 https://indico.cern.ch/event/610889/contributions/2463094/attachments/1412258/2160811/PseudoTopRivet.pdf#page=10)
This is supposed to run at the user analysis level, therefore it does not change the central workflows.
(contents are modified to leave original contributers of the codes)