HepMCCandAlgos/python/flavorHistoryFilter_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 #
0004 # Create sequence to separate HF composition samples.
0005 #
0006 # These are exclusive priorities, so sample "i" will not overlap with "i+1".
0007 # Note that the "dr" values below correspond to the dr between the
0008 # matched genjet, and the sister genjet.
0009 #
0010 # This filter will write out a single unsigned integer:
0011 #
0012 # 1) W+bb with >= 2 jets from the ME (dr > 0.5)
0013 # 2) W+b or W+bb with 1 jet from the ME
0014 # 3) W+cc from the ME (dr > 0.5)
0015 # 4) W+c or W+cc with 1 jet from the ME
0016 # 5) W+bb with 1 jet from the parton shower (dr == 0.0)
0017 # 6) W+cc with 1 jet from the parton shower (dr == 0.0)
0018 #
0019 # These are the "trash bin" samples that we're throwing away:
0020 #
0021 # 7) W+bb with >= 2 partons but 1 jet from the ME (dr == 0.0)
0022 # 8) W+cc with >= 2 partons but 1 jet from the ME (dr == 0.0)
0023 # 9) W+bb with >= 2 partons but 2 jets from the PS (dr > 0.5)
0024 # 10)W+cc with >= 2 partons but 2 jets from the PS (dr > 0.5)
0025 #
0026 # And here is the true "light flavor" sample:
0027 #
0028 # 11) Veto of all the previous (W+ light jets)
0029 #
0030 # The user can select one of these or another to filter on
0031 # via the "pathToSelect" member, which can be 1-11.
0032
0033 flavorHistoryFilter = cms.EDFilter("FlavorHistoryFilter",
0034                                    bsrc = cms.InputTag("bFlavorHistoryProducer", "bPartonFlavorHistory"),
0035                                    csrc = cms.InputTag("cFlavorHistoryProducer", "cPartonFlavorHistory"),
0036                                    pathToSelect = cms.int32(-1),   # no path selected by default
0037                                    dr = cms.double(0.5),           # cutoff delta r
0038                                    verbose = cms.bool(False)
0039                                    )