customiseL1CustomReco.py

Go to the documentation of this file.

import FWCore.ParameterSet.Config as cms

## from CondCore.DBCommon.CondDBSetup_cfi import *
from CondCore.CondDB.CondDB_cfi import *
from PhysicsTools.SelectorUtils.tools.vid_id_tools import *
from PhysicsTools.SelectorUtils.centralIDRegistry import central_id_registry
from JetMETCorrections.Configuration.JetCorrectors_cff import *

def L1NtupleCustomReco(process):


####  Custom Jet reco ####

    # load JEC from SQLite file
    process.load("CondCore.DBCommon.CondDBCommon_cfi")

    ## process.jec = cms.ESSource(
    ##     "PoolDBESSource",
    ##     DBParameters = cms.PSet(
    ##         messageLevel = cms.untracked.int32(0)
    ##         ),
    ##     timetype = cms.string('runnumber'),
    ##     toGet = cms.VPSet(
    ##         cms.PSet(
    ##             record = cms.string('JetCorrectionsRecord'),
    ##             tag    = cms.string('JetCorrectorParametersCollection_Summer15_25nsV6_DATA_AK4PFchs'),
    ##             label  = cms.untracked.string('AK4PFCHS')
    ##             ),
    ##         ), 
    ##     connect = cms.string('sqlite:Summer15_25nsV6_DATA.db')
    ##     # uncomment above tag lines and this comment to use MC JEC
    ##     # connect = cms.string('sqlite:Summer12_V7_MC.db')
    ##     )
    ## 
    ## process.es_prefer_jec = cms.ESPrefer('PoolDBESSource','jec')
    
    # re-apply JEC for AK4 CHS PF jets
    process.load('JetMETCorrections.Configuration.JetCorrectors_cff')
    #process.load('JetMETCorrections.Configuration.JetCorrectionProducers_cff')
    #process.load('JetMETCorrections.Configuration.CorrectedJetProducers_cff')

    #process.ak4PFCHSJetsL1FastL2L3Residual = process.ak4PFCHSJetsL1.clone(correctors = ['ak4PFCHSL1FastL2L3ResidualCorrector'])

####  Custom Met Filter reco

    # load hbhe noise filter result producer
    process.load('CommonTools/RecoAlgos/HBHENoiseFilterResultProducer_cfi')


####  Custom E/Gamma reco ####

    # turn on VID producer, indicate data format  to be
    # DataFormat.AOD or DataFormat.MiniAOD, as appropriate 
    dataFormat = DataFormat.AOD
    switchOnVIDElectronIdProducer(process, dataFormat)
    process.load("RecoEgamma.ElectronIdentification.egmGsfElectronIDs_cfi")
    process.egmGsfElectronIDSequence = cms.Sequence(process.egmGsfElectronIDs)
    # define which IDs we want to produce
    idmod = 'RecoEgamma.ElectronIdentification.Identification.cutBasedElectronID_Spring15_25ns_V1_cff'  
    setupAllVIDIdsInModule(process,idmod,setupVIDElectronSelection)




    process.l1CustomReco = cms.Path(
        process.ak4PFCHSL1FastL2L3ResidualCorrectorChain
        +process.HBHENoiseFilterResultProducer
        +process.egmGsfElectronIDSequence
        )
    
    process.schedule.append(process.l1CustomReco)

    return process