JetAnalysis/python/PatAna_cff.py

import FWCore.ParameterSet.Config as cms

# Pat Jets
from PhysicsTools.PatAlgos.patHeavyIonSequences_cff import *
patJets.jetSource  = cms.InputTag("iterativeConePu5CaloJets")
patJets.addBTagInfo         = False
patJets.addTagInfos         = False
patJets.addDiscriminators   = False
patJets.addAssociatedTracks = False
patJets.addJetCharge        = False
patJets.addJetID            = False
patJets.getJetMCFlavour     = False
patJets.addGenPartonMatch   = True
patJets.addGenJetMatch      = True
patJets.embedGenJetMatch    = True
patJets.embedGenPartonMatch = True
patJets.embedCaloTowers     = False


patJetCorrFactors.useNPV = False
# full reco
icPu5corr = patJetCorrFactors.clone(
  src = cms.InputTag("iterativeConePu5CaloJets"),
  levels = cms.vstring('L2Relative','L3Absolute'),
  payload = cms.string('IC5Calo_2760GeV')
  )

#icPu5patJets = patJets.clone(
#  jetSource = cms.InputTag("iterativeConePu5CaloJets"),
#  genJetMatch = cms.InputTag("icPu5match"),
#  genPartonMatch = cms.InputTag("icPu5parton"),
#  jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr"))
#  )

akPu5PFcorr = icPu5corr.clone(
  src = cms.InputTag("akPu5PFJets"),
  payload = cms.string('AK5PFTowers_hiGoodTightTracks')
  )
akPu5PFpatJets = patJets.clone(
  jetSource = cms.InputTag("akPu5PFJets"),
  genJetMatch = cms.InputTag("akPu5PFmatch"),
  genPartonMatch = cms.InputTag("akPu5PFparton"),
  jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5PFcorr"))
  )

akPu3PFcorr = icPu5corr.clone(
  src = cms.InputTag("akPu3PFJets"),
  payload = cms.string('AK3PFTowers_hiGoodTightTracks')
  )

#akPu3PFpatJets = patJets.clone(
#  jetSource = cms.InputTag("akPu3PFJets"),
#  genJetMatch = cms.InputTag("akPu3PFmatch"),
#  genPartonMatch = cms.InputTag("akPu3PFparton"),
#  jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3PFcorr"))
#  )

akPu5corr = icPu5corr.clone(
    src = cms.InputTag("akPu5CaloJets"),
    payload = cms.string('AK5Calo')
    )

akPu5patJets = patJets.clone(
    jetSource = cms.InputTag("akPu5CaloJets"),
    genJetMatch = cms.InputTag("akPu5match"),
    genPartonMatch = cms.InputTag("akPu5parton"),
    jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5corr"))
    )

akPu3corr = icPu5corr.clone(
    src = cms.InputTag("akPu3CaloJets"),
    payload = cms.string('AK5Calo')
    )

akPu3patJets = patJets.clone(
    jetSource = cms.InputTag("akPu3CaloJets"),
    genJetMatch = cms.InputTag("akPu3match"),
    genPartonMatch = cms.InputTag("akPu3parton"),
    jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3corr"))
    )


# mc matching
patJetPartonMatch.matched = cms.InputTag("hiPartons")

icPu5clean = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone5HiGenJets') )
icPu5match = patJetGenJetMatch.clone(
  src = cms.InputTag("iterativeConePu5CaloJets"),
  matched = cms.InputTag("icPu5clean")
  )
icPu5parton = patJetPartonMatch.clone(
  src = cms.InputTag("iterativeConePu5CaloJets")
        )


akPu5PFclean = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak5HiGenJets') ) 
akPu5PFmatch = patJetGenJetMatch.clone(
  src = cms.InputTag("akPu5PFJets"),
  matched = cms.InputTag("akPu5PFclean")
  )
akPu5PFparton = patJetPartonMatch.clone(
  src = cms.InputTag("akPu5PFJets")
        )

akPu3PFclean = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak3HiGenJets') ) 
akPu3PFmatch = patJetGenJetMatch.clone(
  src = cms.InputTag("akPu3PFJets"),
  matched = cms.InputTag("akPu3PFclean")
  )
akPu3PFparton = patJetPartonMatch.clone(
  src = cms.InputTag("akPu3PFJets")
        )

akPu5match = patJetGenJetMatch.clone(
    src = cms.InputTag("akPu5CaloJets"),
    matched = cms.InputTag("akPu5PFclean")
    )
akPu5parton = patJetPartonMatch.clone(
    src = cms.InputTag("akPu5CaloJets")
    )

akPu3match = patJetGenJetMatch.clone(
    src = cms.InputTag("akPu3CaloJets"),
    matched = cms.InputTag("akPu3PFclean")
    )
akPu3parton = patJetPartonMatch.clone(
    src = cms.InputTag("akPu3CaloJets")
    )

from RecoJets.JetAssociationProducers.ak5JTA_cff import *
from RecoBTag.Configuration.RecoBTag_cff import *
# @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ calo jet
#### b-tagging:
# b-tagging general configuration
# @@@@@@@@@ calo jet
#### b-tagging:
# b-tagging general configuration
icPu5JetTracksAssociatorAtVertex        = ak5JetTracksAssociatorAtVertex.clone()
icPu5JetTracksAssociatorAtVertex.jets   = cms.InputTag("iterativeConePu5CaloJets")
icPu5JetTracksAssociatorAtVertex.tracks = cms.InputTag("hiSelectedTracks")

# impact parameter b-tag
icPu5ImpactParameterTagInfos                = impactParameterTagInfos.clone()
icPu5ImpactParameterTagInfos.jetTracks      = cms.InputTag("icPu5JetTracksAssociatorAtVertex")
icPu5ImpactParameterTagInfos.primaryVertex  = cms.InputTag("hiSelectedVertex")

icPu5TrackCountingHighEffBJetTags          = trackCountingHighEffBJetTags.clone()
icPu5TrackCountingHighEffBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"))
icPu5TrackCountingHighPurBJetTags          = trackCountingHighPurBJetTags.clone()
icPu5TrackCountingHighPurBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"))
icPu5JetProbabilityBJetTags                = jetProbabilityBJetTags.clone()
icPu5JetProbabilityBJetTags.tagInfos       = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"))
icPu5JetBProbabilityBJetTags               = jetBProbabilityBJetTags.clone()
icPu5JetBProbabilityBJetTags.tagInfos      = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"))

# secondary vertex b-tag
icPu5SecondaryVertexTagInfos                     = secondaryVertexTagInfos.clone()
icPu5SecondaryVertexTagInfos.trackIPTagInfos     = cms.InputTag("icPu5ImpactParameterTagInfos")
icPu5SimpleSecondaryVertexBJetTags               = simpleSecondaryVertexBJetTags.clone()
icPu5SimpleSecondaryVertexBJetTags.tagInfos      = cms.VInputTag(cms.InputTag("icPu5SecondaryVertexTagInfos"))
icPu5CombinedSecondaryVertexBJetTags             = combinedSecondaryVertexBJetTags.clone()
icPu5CombinedSecondaryVertexBJetTags.tagInfos    = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"),
                                                                         cms.InputTag("icPu5SecondaryVertexTagInfos"))
icPu5CombinedSecondaryVertexMVABJetTags          = combinedSecondaryVertexMVABJetTags.clone()
icPu5CombinedSecondaryVertexMVABJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"),
                                                                         cms.InputTag("icPu5SecondaryVertexTagInfos"))

# soft muon b-tag
icPu5SoftMuonTagInfos                = softMuonTagInfos.clone()
icPu5SoftMuonTagInfos.jets           = cms.InputTag("iterativeConePu5CaloJets")
icPu5SoftMuonTagInfos.primaryVertex  = cms.InputTag("hiSelectedVertex")
icPu5SoftMuonBJetTags                = softMuonBJetTags.clone()
icPu5SoftMuonBJetTags.tagInfos       = cms.VInputTag(cms.InputTag("icPu5SoftMuonTagInfos"))
icPu5SoftMuonByIP3dBJetTags          = softMuonByIP3dBJetTags.clone()
icPu5SoftMuonByIP3dBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5SoftMuonTagInfos"))
icPu5SoftMuonByPtBJetTags            = softMuonByPtBJetTags.clone()
icPu5SoftMuonByPtBJetTags.tagInfos   = cms.VInputTag(cms.InputTag("icPu5SoftMuonTagInfos"))

# ghost tracks
icPu5GhostTrackVertexTagInfos                 = ghostTrackVertexTagInfos.clone()
icPu5GhostTrackVertexTagInfos.trackIPTagInfos = cms.InputTag("icPu5ImpactParameterTagInfos")
icPu5GhostTrackBJetTags                       = ghostTrackBJetTags.clone()
icPu5GhostTrackBJetTags.tagInfos              = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"),
                                                                      cms.InputTag("icPu5GhostTrackVertexTagInfos"))
# prepare a path running the new modules
icPu5JetTracksAssociator = cms.Sequence(icPu5JetTracksAssociatorAtVertex)
icPu5JetBtaggingIP       = cms.Sequence(icPu5ImpactParameterTagInfos * (icPu5TrackCountingHighEffBJetTags +
                                                                                        icPu5TrackCountingHighPurBJetTags +
                                                                                        icPu5JetProbabilityBJetTags +
                                                                                        icPu5JetBProbabilityBJetTags
                                                                                        )
                                                )

icPu5JetBtaggingSV = cms.Sequence(icPu5ImpactParameterTagInfos *
                                          icPu5SecondaryVertexTagInfos * (icPu5SimpleSecondaryVertexBJetTags +
                                                                                  icPu5CombinedSecondaryVertexBJetTags +
                                                                                  icPu5CombinedSecondaryVertexMVABJetTags
                                                                                  )
                                          +icPu5GhostTrackVertexTagInfos
                                          *icPu5GhostTrackBJetTags
                                          )


icPu5JetBtaggingMu = cms.Sequence(icPu5SoftMuonTagInfos * (icPu5SoftMuonBJetTags +
                                                                           icPu5SoftMuonByIP3dBJetTags +
                                                                           icPu5SoftMuonByPtBJetTags
                                                                           )
                                          )

icPu5JetBtagging = cms.Sequence(icPu5JetBtaggingIP 
                                        *icPu5JetBtaggingSV 
                                        *icPu5JetBtaggingMu
                                        )
#----------------------


#----------------------
icPu5corr    = patJetCorrFactors.clone(src      = cms.InputTag("iterativeConePu5CaloJets"),
                                                       levels   = cms.vstring('L2Relative','L3Absolute'),
                                                       payload  = cms.string('IC5Calo_2760GeV'))
icPu5clean   = heavyIonCleanedGenJets.clone(src = cms.InputTag('iterativeCone5HiGenJets')) # cleans the jets, but NOT the partons
icPu5match   = patJetGenJetMatch.clone(src      = cms.InputTag("iterativeConePu5CaloJets"),
                                                       matched  = cms.InputTag("icPu5clean"))

icPu5parton  = patJetPartonMatch.clone(src      = cms.InputTag("iterativeConePu5CaloJets"))

# ----- flavour bit
icPu5PatJetPartonAssociation       = patJetPartonAssociation.clone(jets    = cms.InputTag("iterativeConePu5CaloJets"),
                                                                                   partons = cms.InputTag("genPartons"),
                                                                                   coneSizeToAssociate = cms.double(0.4))
icPu5PatJetFlavourAssociation      = patJetFlavourAssociation.clone(srcByReference = cms.InputTag("icPu5PatJetPartonAssociation"))

icPu5PatJetFlavourId               = cms.Sequence(icPu5PatJetPartonAssociation*icPu5PatJetFlavourAssociation)

#-------

icPu5patJets = patJets.clone(jetSource            = cms.InputTag("iterativeConePu5CaloJets"),
                                             genJetMatch          = cms.InputTag("icPu5match"),
                                             genPartonMatch       = cms.InputTag("icPu5parton"),
                                             jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr")),
                                             JetPartonMapSource   = cms.InputTag("icPu5PatJetFlavourAssociation"),
                                             trackAssociationSource = cms.InputTag("icPu5JetTracksAssociatorAtVertex"),
                                             discriminatorSources = cms.VInputTag(cms.InputTag("icPu5CombinedSecondaryVertexBJetTags"),
                                                                                  cms.InputTag("icPu5CombinedSecondaryVertexMVABJetTags"),
                                                                                  cms.InputTag("icPu5JetBProbabilityBJetTags"),
                                                                                  cms.InputTag("icPu5JetProbabilityBJetTags"),
                                                                                  cms.InputTag("icPu5SoftMuonByPtBJetTags"),                
                                                                                  cms.InputTag("icPu5SoftMuonByIP3dBJetTags"),
                                                                                  cms.InputTag("icPu5TrackCountingHighEffBJetTags"),
                                                                                  cms.InputTag("icPu5TrackCountingHighPurBJetTags"),
                                                                                  ),
                                             )


#### B-tagging for this bit:
# b-tagging general configuration
akPu3PFJetTracksAssociatorAtVertex        = ak5JetTracksAssociatorAtVertex.clone()
akPu3PFJetTracksAssociatorAtVertex.jets   = cms.InputTag("akPu3PFJets")
akPu3PFJetTracksAssociatorAtVertex.tracks = cms.InputTag("hiSelectedTracks")
# impact parameter b-tag
akPu3PFImpactParameterTagInfos               = impactParameterTagInfos.clone()
akPu3PFImpactParameterTagInfos.jetTracks     = cms.InputTag("akPu3PFJetTracksAssociatorAtVertex")
akPu3PFImpactParameterTagInfos.primaryVertex = cms.InputTag("hiSelectedVertex")
akPu3PFTrackCountingHighEffBJetTags          = trackCountingHighEffBJetTags.clone()
akPu3PFTrackCountingHighEffBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"))
akPu3PFTrackCountingHighPurBJetTags          = trackCountingHighPurBJetTags.clone()
akPu3PFTrackCountingHighPurBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"))
akPu3PFJetProbabilityBJetTags                = jetProbabilityBJetTags.clone()
akPu3PFJetProbabilityBJetTags.tagInfos       = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"))
akPu3PFJetBProbabilityBJetTags               = jetBProbabilityBJetTags.clone()
akPu3PFJetBProbabilityBJetTags.tagInfos      = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"))

# secondary vertex b-tag
akPu3PFSecondaryVertexTagInfos                     = secondaryVertexTagInfos.clone()
akPu3PFSecondaryVertexTagInfos.trackIPTagInfos     = cms.InputTag("akPu3PFImpactParameterTagInfos")
akPu3PFSimpleSecondaryVertexBJetTags               = simpleSecondaryVertexBJetTags.clone()
akPu3PFSimpleSecondaryVertexBJetTags.tagInfos      = cms.VInputTag(cms.InputTag("akPu3PFSecondaryVertexTagInfos"))
akPu3PFCombinedSecondaryVertexBJetTags             = combinedSecondaryVertexBJetTags.clone()
akPu3PFCombinedSecondaryVertexBJetTags.tagInfos    = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"), 
                                                                           cms.InputTag("akPu3PFSecondaryVertexTagInfos"))
akPu3PFCombinedSecondaryVertexMVABJetTags          = combinedSecondaryVertexMVABJetTags.clone()
akPu3PFCombinedSecondaryVertexMVABJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"), 
                                                                           cms.InputTag("akPu3PFSecondaryVertexTagInfos"))

# soft muon b-tag
akPu3PFSoftMuonTagInfos                = softMuonTagInfos.clone()
akPu3PFSoftMuonTagInfos.jets           = cms.InputTag("akPu3PFJets")
akPu3PFSoftMuonTagInfos.primaryVertex  = cms.InputTag("hiSelectedVertex")
akPu3PFSoftMuonBJetTags                = softMuonBJetTags.clone()
akPu3PFSoftMuonBJetTags.tagInfos       = cms.VInputTag(cms.InputTag("akPu3PFSoftMuonTagInfos"))
akPu3PFSoftMuonByIP3dBJetTags          = softMuonByIP3dBJetTags.clone()
akPu3PFSoftMuonByIP3dBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFSoftMuonTagInfos"))
akPu3PFSoftMuonByPtBJetTags            = softMuonByPtBJetTags.clone()
akPu3PFSoftMuonByPtBJetTags.tagInfos   = cms.VInputTag(cms.InputTag("akPu3PFSoftMuonTagInfos"))

# ghost tracks
akPu3PFGhostTrackVertexTagInfos                 = ghostTrackVertexTagInfos.clone()
akPu3PFGhostTrackVertexTagInfos.trackIPTagInfos = cms.InputTag("akPu3PFImpactParameterTagInfos")
akPu3PFGhostTrackBJetTags                       = ghostTrackBJetTags.clone()
akPu3PFGhostTrackBJetTags.tagInfos              = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"),
                                                                        cms.InputTag("akPu3PFGhostTrackVertexTagInfos"))
# prepare a path running the new modules
akPu3PFJetTracksAssociator = cms.Sequence(akPu3PFJetTracksAssociatorAtVertex)
akPu3PFJetBtaggingIP       = cms.Sequence(akPu3PFImpactParameterTagInfos * (akPu3PFTrackCountingHighEffBJetTags +
                                                                                        akPu3PFTrackCountingHighPurBJetTags +
                                                                                        akPu3PFJetProbabilityBJetTags +
                                                                                        akPu3PFJetBProbabilityBJetTags
                                                                                        )
                                                )

akPu3PFJetBtaggingSV = cms.Sequence(akPu3PFImpactParameterTagInfos *
                                          akPu3PFSecondaryVertexTagInfos * (akPu3PFSimpleSecondaryVertexBJetTags +
                                                                                  akPu3PFCombinedSecondaryVertexBJetTags +
                                                                                  akPu3PFCombinedSecondaryVertexMVABJetTags
                                                                                  )
                                          +akPu3PFGhostTrackVertexTagInfos
                                          *akPu3PFGhostTrackBJetTags
                                          )


akPu3PFJetBtaggingMu = cms.Sequence(akPu3PFSoftMuonTagInfos * (akPu3PFSoftMuonBJetTags +
                                                                           akPu3PFSoftMuonByIP3dBJetTags +
                                                                           akPu3PFSoftMuonByPtBJetTags
                                                                           )
                                          )

akPu3PFJetBtagging = cms.Sequence(akPu3PFJetBtaggingIP 
                                        *akPu3PFJetBtaggingSV 
                                        *akPu3PFJetBtaggingMu
                                        )

#__________________________________________________________
# ----- flavour bit
akPu3PFPatJetPartonAssociation       = patJetPartonAssociation.clone(jets    = cms.InputTag("akPu3PFJets"),
                                                                                     partons = cms.InputTag("genPartons"),
                                                                                     coneSizeToAssociate = cms.double(0.4))
akPu3PFPatJetFlavourAssociation      = patJetFlavourAssociation.clone(srcByReference = cms.InputTag("akPu3PFPatJetPartonAssociation"))

akPu3PFPatJetFlavourId               = cms.Sequence(akPu3PFPatJetPartonAssociation*akPu3PFPatJetFlavourAssociation)
#
#-------
akPu3PFcorr = patJetCorrFactors.clone(src     = cms.InputTag("akPu3PFJets"),
                                      levels  = cms.vstring('L2Relative','L3Absolute'),
                                      payload = cms.string('AK3PFTowers_hiGoodTightTracks')
                                                      )
akPu3PFclean   = heavyIonCleanedGenJets.clone(src = cms.InputTag('ak3HiGenJets'))
akPu3PFmatch   = patJetGenJetMatch.clone(src      = cms.InputTag("akPu3PFJets"),
                                                         matched  = cms.InputTag("akPu3PFclean"))
akPu3PFparton  = patJetPartonMatch.clone(src      = cms.InputTag("akPu3PFJets"))
akPu3PFpatJets = patJets.clone(jetSource            = cms.InputTag("akPu3PFJets"),
                                               genJetMatch          = cms.InputTag("akPu3PFmatch"),
                                               genPartonMatch       = cms.InputTag("akPu3PFparton"),
                                               jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3PFcorr")),
                                               JetPartonMapSource   = cms.InputTag("akPu3PFPatJetFlavourAssociation"),
                                               trackAssociationSource = cms.InputTag("akPu3PFJetTracksAssociatorAtVertex"),
                                               discriminatorSources = cms.VInputTag(cms.InputTag("akPu3PFCombinedSecondaryVertexBJetTags"),
                                                                                    cms.InputTag("akPu3PFCombinedSecondaryVertexMVABJetTags"),
                                                                                    cms.InputTag("akPu3PFJetBProbabilityBJetTags"),
                                                                                    cms.InputTag("akPu3PFJetProbabilityBJetTags"),
                                                                                    cms.InputTag("akPu3PFSoftMuonByPtBJetTags"),                
                                                                                    cms.InputTag("akPu3PFSoftMuonByIP3dBJetTags"),
                                                                                    cms.InputTag("akPu3PFTrackCountingHighEffBJetTags"),
                                                                                    cms.InputTag("akPu3PFTrackCountingHighPurBJetTags"),
                                                                                    ),
                                               )

# end batagging


# === mc sequences ===
icPu5patSequence = cms.Sequence(icPu5corr * icPu5clean * icPu5match * icPu5parton  *  icPu5patJets)
akPu5PFpatSequence = cms.Sequence(akPu5PFcorr * akPu5PFclean * akPu5PFmatch * akPu5PFparton * akPu5PFpatJets)
akPu3PFpatSequence = cms.Sequence(akPu3PFcorr * akPu3PFclean * akPu3PFmatch * akPu3PFparton * akPu3PFpatJets)
akPu5patSequence = cms.Sequence(akPu5corr * akPu5PFclean * akPu5match * akPu5parton * akPu5patJets)
akPu3patSequence = cms.Sequence(akPu3corr * akPu3PFclean * akPu3match * akPu3parton * akPu3patJets)

icPu5patSequence_withBtagging = cms.Sequence(icPu5corr * icPu5clean * icPu5match * icPu5parton  * icPu5PatJetFlavourId * icPu5JetTracksAssociator * icPu5JetBtagging * icPu5patJets)
akPu3PFpatSequence_withBtagging = cms.Sequence(akPu3PFcorr * akPu3PFclean * akPu3PFmatch * akPu3PFparton * akPu3PFPatJetFlavourId * akPu3PFJetTracksAssociator *akPu3PFJetBtagging * akPu3PFpatJets)


# === data sequences ===
# Note still need to use enableData function in cfg to remove mc dep of patjet
icPu5patSequence_data = cms.Sequence( icPu5corr * icPu5patJets )
akPu5PFpatSequence_data = cms.Sequence(akPu5PFcorr * akPu5PFpatJets )
akPu3PFpatSequence_data = cms.Sequence( akPu3PFcorr * akPu3PFpatJets )
akPu5patSequence_data = cms.Sequence( akPu5corr * akPu5patJets )
akPu3patSequence_data = cms.Sequence( akPu3corr * akPu3patJets )

icPu5patSequence_withBtagging_data = cms.Sequence( icPu5corr * icPu5JetTracksAssociator * icPu5JetBtagging * icPu5JetTracksAssociator * icPu5JetBtagging * icPu5patJets )
akPu3PFpatSequence_withBtagging_data = cms.Sequence( akPu3PFcorr * akPu3PFJetTracksAssociator *akPu3PFJetBtagging * akPu3PFJetTracksAssociator *akPu3PFJetBtagging * akPu3PFpatJets )
Revision:	1.8
Committed:	Mon Apr 16 22:08:38 2012 UTC (13 years ago) by mnguyen
Content type:	text/x-python
Branch:	MAIN
CVS Tags:	HiForest_V02_00
Changes since 1.7:	+269 -20 lines
Log Message:	Added b-tagging (turned off by default)
#	User	Rev	Content
1	frankma	1.1	import FWCore.ParameterSet.Config as cms
2
3			# Pat Jets
4			from PhysicsTools.PatAlgos.patHeavyIonSequences_cff import *
5			patJets.jetSource = cms.InputTag("iterativeConePu5CaloJets")
6			patJets.addBTagInfo = False
7			patJets.addTagInfos = False
8			patJets.addDiscriminators = False
9			patJets.addAssociatedTracks = False
10			patJets.addJetCharge = False
11			patJets.addJetID = False
12			patJets.getJetMCFlavour = False
13			patJets.addGenPartonMatch = True
14			patJets.addGenJetMatch = True
15			patJets.embedGenJetMatch = True
16			patJets.embedGenPartonMatch = True
17			patJets.embedCaloTowers = False
18
19	mnguyen	1.8
20	yilmaz	1.5	patJetCorrFactors.useNPV = False
21	frankma	1.1	# full reco
22			icPu5corr = patJetCorrFactors.clone(
23			src = cms.InputTag("iterativeConePu5CaloJets"),
24			levels = cms.vstring('L2Relative','L3Absolute'),
25	yilmaz	1.7	payload = cms.string('IC5Calo_2760GeV')
26	frankma	1.1	)
27	yilmaz	1.6
28	mnguyen	1.8	#icPu5patJets = patJets.clone(
29			# jetSource = cms.InputTag("iterativeConePu5CaloJets"),
30			# genJetMatch = cms.InputTag("icPu5match"),
31			# genPartonMatch = cms.InputTag("icPu5parton"),
32			# jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr"))
33			# )
34	frankma	1.1
35			akPu5PFcorr = icPu5corr.clone(
36			src = cms.InputTag("akPu5PFJets"),
37	yilmaz	1.7	payload = cms.string('AK5PFTowers_hiGoodTightTracks')
38	frankma	1.1	)
39			akPu5PFpatJets = patJets.clone(
40			jetSource = cms.InputTag("akPu5PFJets"),
41			genJetMatch = cms.InputTag("akPu5PFmatch"),
42			genPartonMatch = cms.InputTag("akPu5PFparton"),
43			jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5PFcorr"))
44			)
45
46			akPu3PFcorr = icPu5corr.clone(
47			src = cms.InputTag("akPu3PFJets"),
48	yilmaz	1.7	payload = cms.string('AK3PFTowers_hiGoodTightTracks')
49	frankma	1.1	)
50	yilmaz	1.6
51	mnguyen	1.8	#akPu3PFpatJets = patJets.clone(
52			# jetSource = cms.InputTag("akPu3PFJets"),
53			# genJetMatch = cms.InputTag("akPu3PFmatch"),
54			# genPartonMatch = cms.InputTag("akPu3PFparton"),
55			# jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3PFcorr"))
56			# )
57	frankma	1.1
58	yilmaz	1.6	akPu5corr = icPu5corr.clone(
59			src = cms.InputTag("akPu5CaloJets"),
60			payload = cms.string('AK5Calo')
61			)
62
63			akPu5patJets = patJets.clone(
64			jetSource = cms.InputTag("akPu5CaloJets"),
65			genJetMatch = cms.InputTag("akPu5match"),
66			genPartonMatch = cms.InputTag("akPu5parton"),
67			jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu5corr"))
68			)
69
70			akPu3corr = icPu5corr.clone(
71			src = cms.InputTag("akPu3CaloJets"),
72	yilmaz	1.7	payload = cms.string('AK5Calo')
73	yilmaz	1.6	)
74
75			akPu3patJets = patJets.clone(
76			jetSource = cms.InputTag("akPu3CaloJets"),
77			genJetMatch = cms.InputTag("akPu3match"),
78			genPartonMatch = cms.InputTag("akPu3parton"),
79			jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3corr"))
80			)
81
82
83	frankma	1.1
84			# mc matching
85			patJetPartonMatch.matched = cms.InputTag("hiPartons")
86
87			icPu5clean = heavyIonCleanedGenJets.clone( src = cms.InputTag('iterativeCone5HiGenJets') )
88			icPu5match = patJetGenJetMatch.clone(
89			src = cms.InputTag("iterativeConePu5CaloJets"),
90			matched = cms.InputTag("icPu5clean")
91			)
92			icPu5parton = patJetPartonMatch.clone(
93			src = cms.InputTag("iterativeConePu5CaloJets")
94			)
95
96	yilmaz	1.3
97	frankma	1.1	akPu5PFclean = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak5HiGenJets') )
98			akPu5PFmatch = patJetGenJetMatch.clone(
99			src = cms.InputTag("akPu5PFJets"),
100			matched = cms.InputTag("akPu5PFclean")
101			)
102			akPu5PFparton = patJetPartonMatch.clone(
103			src = cms.InputTag("akPu5PFJets")
104			)
105
106			akPu3PFclean = heavyIonCleanedGenJets.clone( src = cms.InputTag('ak3HiGenJets') )
107			akPu3PFmatch = patJetGenJetMatch.clone(
108			src = cms.InputTag("akPu3PFJets"),
109			matched = cms.InputTag("akPu3PFclean")
110			)
111			akPu3PFparton = patJetPartonMatch.clone(
112			src = cms.InputTag("akPu3PFJets")
113			)
114
115	yilmaz	1.6	akPu5match = patJetGenJetMatch.clone(
116			src = cms.InputTag("akPu5CaloJets"),
117			matched = cms.InputTag("akPu5PFclean")
118			)
119			akPu5parton = patJetPartonMatch.clone(
120			src = cms.InputTag("akPu5CaloJets")
121			)
122
123			akPu3match = patJetGenJetMatch.clone(
124			src = cms.InputTag("akPu3CaloJets"),
125			matched = cms.InputTag("akPu3PFclean")
126			)
127			akPu3parton = patJetPartonMatch.clone(
128			src = cms.InputTag("akPu3CaloJets")
129			)
130
131	mnguyen	1.8	from RecoJets.JetAssociationProducers.ak5JTA_cff import *
132			from RecoBTag.Configuration.RecoBTag_cff import *
133			# @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ calo jet
134			#### b-tagging:
135			# b-tagging general configuration
136			# @@@@@@@@@ calo jet
137			#### b-tagging:
138			# b-tagging general configuration
139			icPu5JetTracksAssociatorAtVertex = ak5JetTracksAssociatorAtVertex.clone()
140			icPu5JetTracksAssociatorAtVertex.jets = cms.InputTag("iterativeConePu5CaloJets")
141			icPu5JetTracksAssociatorAtVertex.tracks = cms.InputTag("hiSelectedTracks")
142
143			# impact parameter b-tag
144			icPu5ImpactParameterTagInfos = impactParameterTagInfos.clone()
145			icPu5ImpactParameterTagInfos.jetTracks = cms.InputTag("icPu5JetTracksAssociatorAtVertex")
146			icPu5ImpactParameterTagInfos.primaryVertex = cms.InputTag("hiSelectedVertex")
147
148			icPu5TrackCountingHighEffBJetTags = trackCountingHighEffBJetTags.clone()
149			icPu5TrackCountingHighEffBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"))
150			icPu5TrackCountingHighPurBJetTags = trackCountingHighPurBJetTags.clone()
151			icPu5TrackCountingHighPurBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"))
152			icPu5JetProbabilityBJetTags = jetProbabilityBJetTags.clone()
153			icPu5JetProbabilityBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"))
154			icPu5JetBProbabilityBJetTags = jetBProbabilityBJetTags.clone()
155			icPu5JetBProbabilityBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"))
156
157			# secondary vertex b-tag
158			icPu5SecondaryVertexTagInfos = secondaryVertexTagInfos.clone()
159			icPu5SecondaryVertexTagInfos.trackIPTagInfos = cms.InputTag("icPu5ImpactParameterTagInfos")
160			icPu5SimpleSecondaryVertexBJetTags = simpleSecondaryVertexBJetTags.clone()
161			icPu5SimpleSecondaryVertexBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5SecondaryVertexTagInfos"))
162			icPu5CombinedSecondaryVertexBJetTags = combinedSecondaryVertexBJetTags.clone()
163			icPu5CombinedSecondaryVertexBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"),
164			cms.InputTag("icPu5SecondaryVertexTagInfos"))
165			icPu5CombinedSecondaryVertexMVABJetTags = combinedSecondaryVertexMVABJetTags.clone()
166			icPu5CombinedSecondaryVertexMVABJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"),
167			cms.InputTag("icPu5SecondaryVertexTagInfos"))
168
169			# soft muon b-tag
170			icPu5SoftMuonTagInfos = softMuonTagInfos.clone()
171			icPu5SoftMuonTagInfos.jets = cms.InputTag("iterativeConePu5CaloJets")
172			icPu5SoftMuonTagInfos.primaryVertex = cms.InputTag("hiSelectedVertex")
173			icPu5SoftMuonBJetTags = softMuonBJetTags.clone()
174			icPu5SoftMuonBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5SoftMuonTagInfos"))
175			icPu5SoftMuonByIP3dBJetTags = softMuonByIP3dBJetTags.clone()
176			icPu5SoftMuonByIP3dBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5SoftMuonTagInfos"))
177			icPu5SoftMuonByPtBJetTags = softMuonByPtBJetTags.clone()
178			icPu5SoftMuonByPtBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5SoftMuonTagInfos"))
179
180			# ghost tracks
181			icPu5GhostTrackVertexTagInfos = ghostTrackVertexTagInfos.clone()
182			icPu5GhostTrackVertexTagInfos.trackIPTagInfos = cms.InputTag("icPu5ImpactParameterTagInfos")
183			icPu5GhostTrackBJetTags = ghostTrackBJetTags.clone()
184			icPu5GhostTrackBJetTags.tagInfos = cms.VInputTag(cms.InputTag("icPu5ImpactParameterTagInfos"),
185			cms.InputTag("icPu5GhostTrackVertexTagInfos"))
186			# prepare a path running the new modules
187			icPu5JetTracksAssociator = cms.Sequence(icPu5JetTracksAssociatorAtVertex)
188			icPu5JetBtaggingIP = cms.Sequence(icPu5ImpactParameterTagInfos * (icPu5TrackCountingHighEffBJetTags +
189			icPu5TrackCountingHighPurBJetTags +
190			icPu5JetProbabilityBJetTags +
191			icPu5JetBProbabilityBJetTags
192			)
193			)
194
195			icPu5JetBtaggingSV = cms.Sequence(icPu5ImpactParameterTagInfos *
196			icPu5SecondaryVertexTagInfos * (icPu5SimpleSecondaryVertexBJetTags +
197			icPu5CombinedSecondaryVertexBJetTags +
198			icPu5CombinedSecondaryVertexMVABJetTags
199			)
200			+icPu5GhostTrackVertexTagInfos
201			*icPu5GhostTrackBJetTags
202			)
203
204
205			icPu5JetBtaggingMu = cms.Sequence(icPu5SoftMuonTagInfos * (icPu5SoftMuonBJetTags +
206			icPu5SoftMuonByIP3dBJetTags +
207			icPu5SoftMuonByPtBJetTags
208			)
209			)
210
211			icPu5JetBtagging = cms.Sequence(icPu5JetBtaggingIP
212			*icPu5JetBtaggingSV
213			*icPu5JetBtaggingMu
214			)
215			#----------------------
216
217
218			#----------------------
219			icPu5corr = patJetCorrFactors.clone(src = cms.InputTag("iterativeConePu5CaloJets"),
220			levels = cms.vstring('L2Relative','L3Absolute'),
221			payload = cms.string('IC5Calo_2760GeV'))
222			icPu5clean = heavyIonCleanedGenJets.clone(src = cms.InputTag('iterativeCone5HiGenJets')) # cleans the jets, but NOT the partons
223			icPu5match = patJetGenJetMatch.clone(src = cms.InputTag("iterativeConePu5CaloJets"),
224			matched = cms.InputTag("icPu5clean"))
225
226			icPu5parton = patJetPartonMatch.clone(src = cms.InputTag("iterativeConePu5CaloJets"))
227
228			# ----- flavour bit
229			icPu5PatJetPartonAssociation = patJetPartonAssociation.clone(jets = cms.InputTag("iterativeConePu5CaloJets"),
230			partons = cms.InputTag("genPartons"),
231			coneSizeToAssociate = cms.double(0.4))
232			icPu5PatJetFlavourAssociation = patJetFlavourAssociation.clone(srcByReference = cms.InputTag("icPu5PatJetPartonAssociation"))
233
234			icPu5PatJetFlavourId = cms.Sequence(icPu5PatJetPartonAssociation*icPu5PatJetFlavourAssociation)
235
236			#-------
237
238			icPu5patJets = patJets.clone(jetSource = cms.InputTag("iterativeConePu5CaloJets"),
239			genJetMatch = cms.InputTag("icPu5match"),
240			genPartonMatch = cms.InputTag("icPu5parton"),
241			jetCorrFactorsSource = cms.VInputTag(cms.InputTag("icPu5corr")),
242			JetPartonMapSource = cms.InputTag("icPu5PatJetFlavourAssociation"),
243			trackAssociationSource = cms.InputTag("icPu5JetTracksAssociatorAtVertex"),
244			discriminatorSources = cms.VInputTag(cms.InputTag("icPu5CombinedSecondaryVertexBJetTags"),
245			cms.InputTag("icPu5CombinedSecondaryVertexMVABJetTags"),
246			cms.InputTag("icPu5JetBProbabilityBJetTags"),
247			cms.InputTag("icPu5JetProbabilityBJetTags"),
248			cms.InputTag("icPu5SoftMuonByPtBJetTags"),
249			cms.InputTag("icPu5SoftMuonByIP3dBJetTags"),
250			cms.InputTag("icPu5TrackCountingHighEffBJetTags"),
251			cms.InputTag("icPu5TrackCountingHighPurBJetTags"),
252			),
253			)
254
255
256			#### B-tagging for this bit:
257			# b-tagging general configuration
258			akPu3PFJetTracksAssociatorAtVertex = ak5JetTracksAssociatorAtVertex.clone()
259			akPu3PFJetTracksAssociatorAtVertex.jets = cms.InputTag("akPu3PFJets")
260			akPu3PFJetTracksAssociatorAtVertex.tracks = cms.InputTag("hiSelectedTracks")
261			# impact parameter b-tag
262			akPu3PFImpactParameterTagInfos = impactParameterTagInfos.clone()
263			akPu3PFImpactParameterTagInfos.jetTracks = cms.InputTag("akPu3PFJetTracksAssociatorAtVertex")
264			akPu3PFImpactParameterTagInfos.primaryVertex = cms.InputTag("hiSelectedVertex")
265			akPu3PFTrackCountingHighEffBJetTags = trackCountingHighEffBJetTags.clone()
266			akPu3PFTrackCountingHighEffBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"))
267			akPu3PFTrackCountingHighPurBJetTags = trackCountingHighPurBJetTags.clone()
268			akPu3PFTrackCountingHighPurBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"))
269			akPu3PFJetProbabilityBJetTags = jetProbabilityBJetTags.clone()
270			akPu3PFJetProbabilityBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"))
271			akPu3PFJetBProbabilityBJetTags = jetBProbabilityBJetTags.clone()
272			akPu3PFJetBProbabilityBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"))
273
274			# secondary vertex b-tag
275			akPu3PFSecondaryVertexTagInfos = secondaryVertexTagInfos.clone()
276			akPu3PFSecondaryVertexTagInfos.trackIPTagInfos = cms.InputTag("akPu3PFImpactParameterTagInfos")
277			akPu3PFSimpleSecondaryVertexBJetTags = simpleSecondaryVertexBJetTags.clone()
278			akPu3PFSimpleSecondaryVertexBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFSecondaryVertexTagInfos"))
279			akPu3PFCombinedSecondaryVertexBJetTags = combinedSecondaryVertexBJetTags.clone()
280			akPu3PFCombinedSecondaryVertexBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"),
281			cms.InputTag("akPu3PFSecondaryVertexTagInfos"))
282			akPu3PFCombinedSecondaryVertexMVABJetTags = combinedSecondaryVertexMVABJetTags.clone()
283			akPu3PFCombinedSecondaryVertexMVABJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"),
284			cms.InputTag("akPu3PFSecondaryVertexTagInfos"))
285
286			# soft muon b-tag
287			akPu3PFSoftMuonTagInfos = softMuonTagInfos.clone()
288			akPu3PFSoftMuonTagInfos.jets = cms.InputTag("akPu3PFJets")
289			akPu3PFSoftMuonTagInfos.primaryVertex = cms.InputTag("hiSelectedVertex")
290			akPu3PFSoftMuonBJetTags = softMuonBJetTags.clone()
291			akPu3PFSoftMuonBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFSoftMuonTagInfos"))
292			akPu3PFSoftMuonByIP3dBJetTags = softMuonByIP3dBJetTags.clone()
293			akPu3PFSoftMuonByIP3dBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFSoftMuonTagInfos"))
294			akPu3PFSoftMuonByPtBJetTags = softMuonByPtBJetTags.clone()
295			akPu3PFSoftMuonByPtBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFSoftMuonTagInfos"))
296
297			# ghost tracks
298			akPu3PFGhostTrackVertexTagInfos = ghostTrackVertexTagInfos.clone()
299			akPu3PFGhostTrackVertexTagInfos.trackIPTagInfos = cms.InputTag("akPu3PFImpactParameterTagInfos")
300			akPu3PFGhostTrackBJetTags = ghostTrackBJetTags.clone()
301			akPu3PFGhostTrackBJetTags.tagInfos = cms.VInputTag(cms.InputTag("akPu3PFImpactParameterTagInfos"),
302			cms.InputTag("akPu3PFGhostTrackVertexTagInfos"))
303			# prepare a path running the new modules
304			akPu3PFJetTracksAssociator = cms.Sequence(akPu3PFJetTracksAssociatorAtVertex)
305			akPu3PFJetBtaggingIP = cms.Sequence(akPu3PFImpactParameterTagInfos * (akPu3PFTrackCountingHighEffBJetTags +
306			akPu3PFTrackCountingHighPurBJetTags +
307			akPu3PFJetProbabilityBJetTags +
308			akPu3PFJetBProbabilityBJetTags
309			)
310			)
311
312			akPu3PFJetBtaggingSV = cms.Sequence(akPu3PFImpactParameterTagInfos *
313			akPu3PFSecondaryVertexTagInfos * (akPu3PFSimpleSecondaryVertexBJetTags +
314			akPu3PFCombinedSecondaryVertexBJetTags +
315			akPu3PFCombinedSecondaryVertexMVABJetTags
316			)
317			+akPu3PFGhostTrackVertexTagInfos
318			*akPu3PFGhostTrackBJetTags
319			)
320
321
322			akPu3PFJetBtaggingMu = cms.Sequence(akPu3PFSoftMuonTagInfos * (akPu3PFSoftMuonBJetTags +
323			akPu3PFSoftMuonByIP3dBJetTags +
324			akPu3PFSoftMuonByPtBJetTags
325			)
326			)
327
328			akPu3PFJetBtagging = cms.Sequence(akPu3PFJetBtaggingIP
329			*akPu3PFJetBtaggingSV
330			*akPu3PFJetBtaggingMu
331			)
332
333			#__________________________________________________________
334			# ----- flavour bit
335			akPu3PFPatJetPartonAssociation = patJetPartonAssociation.clone(jets = cms.InputTag("akPu3PFJets"),
336			partons = cms.InputTag("genPartons"),
337			coneSizeToAssociate = cms.double(0.4))
338			akPu3PFPatJetFlavourAssociation = patJetFlavourAssociation.clone(srcByReference = cms.InputTag("akPu3PFPatJetPartonAssociation"))
339
340			akPu3PFPatJetFlavourId = cms.Sequence(akPu3PFPatJetPartonAssociation*akPu3PFPatJetFlavourAssociation)
341			#
342			#-------
343			akPu3PFcorr = patJetCorrFactors.clone(src = cms.InputTag("akPu3PFJets"),
344			levels = cms.vstring('L2Relative','L3Absolute'),
345			payload = cms.string('AK3PFTowers_hiGoodTightTracks')
346			)
347			akPu3PFclean = heavyIonCleanedGenJets.clone(src = cms.InputTag('ak3HiGenJets'))
348			akPu3PFmatch = patJetGenJetMatch.clone(src = cms.InputTag("akPu3PFJets"),
349			matched = cms.InputTag("akPu3PFclean"))
350			akPu3PFparton = patJetPartonMatch.clone(src = cms.InputTag("akPu3PFJets"))
351			akPu3PFpatJets = patJets.clone(jetSource = cms.InputTag("akPu3PFJets"),
352			genJetMatch = cms.InputTag("akPu3PFmatch"),
353			genPartonMatch = cms.InputTag("akPu3PFparton"),
354			jetCorrFactorsSource = cms.VInputTag(cms.InputTag("akPu3PFcorr")),
355			JetPartonMapSource = cms.InputTag("akPu3PFPatJetFlavourAssociation"),
356			trackAssociationSource = cms.InputTag("akPu3PFJetTracksAssociatorAtVertex"),
357			discriminatorSources = cms.VInputTag(cms.InputTag("akPu3PFCombinedSecondaryVertexBJetTags"),
358			cms.InputTag("akPu3PFCombinedSecondaryVertexMVABJetTags"),
359			cms.InputTag("akPu3PFJetBProbabilityBJetTags"),
360			cms.InputTag("akPu3PFJetProbabilityBJetTags"),
361			cms.InputTag("akPu3PFSoftMuonByPtBJetTags"),
362			cms.InputTag("akPu3PFSoftMuonByIP3dBJetTags"),
363			cms.InputTag("akPu3PFTrackCountingHighEffBJetTags"),
364			cms.InputTag("akPu3PFTrackCountingHighPurBJetTags"),
365			),
366			)
367
368			# end batagging
369
370
371	frankma	1.1	# === mc sequences ===
372	mnguyen	1.8	icPu5patSequence = cms.Sequence(icPu5corr * icPu5clean * icPu5match * icPu5parton * icPu5patJets)
373	frankma	1.1	akPu5PFpatSequence = cms.Sequence(akPu5PFcorr * akPu5PFclean * akPu5PFmatch * akPu5PFparton * akPu5PFpatJets)
374			akPu3PFpatSequence = cms.Sequence(akPu3PFcorr * akPu3PFclean * akPu3PFmatch * akPu3PFparton * akPu3PFpatJets)
375	yilmaz	1.6	akPu5patSequence = cms.Sequence(akPu5corr * akPu5PFclean * akPu5match * akPu5parton * akPu5patJets)
376			akPu3patSequence = cms.Sequence(akPu3corr * akPu3PFclean * akPu3match * akPu3parton * akPu3patJets)
377	mnguyen	1.8
378			icPu5patSequence_withBtagging = cms.Sequence(icPu5corr * icPu5clean * icPu5match * icPu5parton * icPu5PatJetFlavourId * icPu5JetTracksAssociator * icPu5JetBtagging * icPu5patJets)
379			akPu3PFpatSequence_withBtagging = cms.Sequence(akPu3PFcorr * akPu3PFclean * akPu3PFmatch * akPu3PFparton * akPu3PFPatJetFlavourId * akPu3PFJetTracksAssociator akPu3PFJetBtagging akPu3PFpatJets)
380
381
382			# === data sequences ===
383			# Note still need to use enableData function in cfg to remove mc dep of patjet
384			icPu5patSequence_data = cms.Sequence( icPu5corr * icPu5patJets )
385			akPu5PFpatSequence_data = cms.Sequence(akPu5PFcorr * akPu5PFpatJets )
386			akPu3PFpatSequence_data = cms.Sequence( akPu3PFcorr * akPu3PFpatJets )
387			akPu5patSequence_data = cms.Sequence( akPu5corr * akPu5patJets )
388			akPu3patSequence_data = cms.Sequence( akPu3corr * akPu3patJets )
389
390			icPu5patSequence_withBtagging_data = cms.Sequence( icPu5corr * icPu5JetTracksAssociator * icPu5JetBtagging * icPu5JetTracksAssociator * icPu5JetBtagging * icPu5patJets )
391			akPu3PFpatSequence_withBtagging_data = cms.Sequence( akPu3PFcorr * akPu3PFJetTracksAssociator akPu3PFJetBtagging akPu3PFJetTracksAssociator akPu3PFJetBtagging akPu3PFpatJets )