payloads = [
    'Jec11_V3_L1FastJet_AK5PFchs.txt', 'Jec11_V3_L2Relative_AK5PFchs.txt',
    'Jec11_V3_L3Absolute_AK5PFchs.txt', 'Jec11_V3_Uncertainty_AK5PFchs.txt'
]

## Maximal Number of Events
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(-1))

from Analysis.SHyFT.shyftPFSelection_cfi import shyftPFSelection as shyftPFSelectionInput

process.load('Analysis.PdfWeights.pdfWeightProducer_cfi')

process.pfShyftProducer = cms.EDFilter(
    'EDSHyFTPFSelector',
    shyftPFSelection=shyftPFSelectionInput.clone(
        jetSrc=cms.InputTag('goodPatJetsPFlow'),
        jecPayloads=cms.vstring(payloads)))

## std sequence to produce the kinematic fit for semi-leptonic events
process.load('TopQuarkAnalysis.TopKinFitter.TtSemiLepKinFitProducer_Muons_cfi')

process.kinFitTtSemiLepEvent.jets = cms.InputTag("pfShyftProducer", "jets")
process.kinFitTtSemiLepEvent.mets = cms.InputTag("pfShyftProducer", "MET")
process.kinFitTtSemiLepEvent.leps = cms.InputTag("pfShyftProducer", "muons")

process.load("TopQuarkAnalysis.TopKinFitter.TtSemiLepKinFitProducer_Muons_cfi")

process.pfShyftTupleJets = cms.EDProducer(
    "CandViewNtpProducer",
    src=cms.InputTag("pfShyftProducer", "jets"),
    lazyParser=cms.untracked.bool(True),
Beispiel #2
0

## Maximum number of events ("-1" means all events)
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )

from Analysis.SHyFT.shyftPFSelection_cfi import shyftPFSelection as shyftPFSelectionInput

process.pileup = cms.EDFilter('PileUpProducer',
                              pvSrc = cms.InputTag('goodOfflinePrimaryVertices')
                              )

## Define various producers
process.pfShyftProducerAK5 = cms.EDFilter('EDSHyFTPFSelector',
                                          shyftPFSelection = shyftPFSelectionInput.clone(
                                          jetSrc = cms.InputTag('goodPatJetsPFlow'),
                                          rhoSrc = cms.InputTag('kt6PFJets', 'rho'),
                                          doElectrons = cms.bool(inputMuOrEle), # 0 = mu, 1 = ele
                                          jecPayloads = cms.vstring( payloads )
                                          )
                                          )
#!!!Change below doElectrons from False to True
process.pfShyftProducerAK5Loose = cms.EDFilter('EDSHyFTPFSelector',
                                               shyftPFSelection = shyftPFSelectionInput.clone(
                                               jetSrc = cms.InputTag('goodPatJetsPFlow'),
                                               muonSrc = cms.InputTag('selectedPatMuonsPFlowLoose'),
                                               electronSrc = cms.InputTag('selectedPatElectronsPFlowLoose'),
                                               doElectrons = cms.bool(inputMuOrEle), # 0 = mu, 1 = ele
                                               rhoSrc = cms.InputTag('kt6PFJets', 'rho'),
                                               jecPayloads = cms.vstring( payloads ),
                                               removeLooseLep = cms.bool(True)
                                               )
                                               )
        process.hltSelection = cms.Sequence( ~process.hltSelectionMu * process.hltSelectionEle )
    else :
        process.hltSelection = cms.Sequence( process.hltSelectionMu )


## Maximal Number of Events
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )

from Analysis.SHyFT.shyftPFSelection_cfi import shyftPFSelection as shyftPFSelectionInput



process.pfShyftProducer = cms.EDFilter('EDSHyFTPFSelector',
                                    shyftPFSelection = shyftPFSelectionInput.clone(
                                           jetSrc = cms.InputTag('goodPatJetsPFlow'),
                                           muonSrc = cms.InputTag('selectedPatMuonsLoosePFlow'),
                                           removeLooseLep = cms.bool(True),
                                           jecPayloads = cms.vstring( payloads )
                                        )
                                    )

process.pfShyftProducerPrunedCA8 = cms.EDFilter('EDSHyFTPFSelector',
                                            shyftPFSelection = shyftPFSelectionInput.clone(
                                                jetSrc = cms.InputTag('goodPatJetsCA8PrunedPF'),
                                                jecPayloads = cms.vstring( payloads )
                                                )
                                            )

process.pfShyftProducer.shyftPFSelection.cutsToIgnore.append('== 1 Tight Lepton')
process.pfShyftProducer.shyftPFSelection.cutsToIgnore.append('0 other lepton')
process.pfShyftProducer.shyftPFSelection.cutsToIgnore.append('>=1 Jets')

## Maximal Number of Events
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )

from Analysis.SHyFT.shyftPFSelection_cfi import shyftPFSelection as shyftPFSelectionInput



process.pfShyftProducer = cms.EDFilter('EDSHyFTPFSelector',
                                    shyftPFSelection = shyftPFSelectionInput.clone(
                                           muonSrc = cms.InputTag('selectedPatMuonsLoosePFlow'),
                                           electronSrc = cms.InputTag('selectedPatElectronsLoosePFlow'),
                                           muonIdPFTight = shyftPFSelectionInput.muonIdPFTight.clone(
                                               cutsToIgnore = cms.vstring('PFIso')
                                               ),
                                           electronIdPFTight = shyftPFSelectionInput.electronIdPFTight.clone(
                                               cutsToIgnore = cms.vstring('PFIso')
                                               ),
                                           removeLooseLep = cms.bool(True)
                                             # replace things here if you want
                                        )
                                    )


## std sequence to produce the kinematic fit for semi-leptonic events
process.load('TopQuarkAnalysis.TopKinFitter.TtSemiLepKinFitProducer_Muons_cfi')

process.kinFitTtSemiLepEvent.jets = cms.InputTag("pfShyftProducer", "jets")
process.kinFitTtSemiLepEvent.mets = cms.InputTag("pfShyftProducer", "MET")
process.kinFitTtSemiLepEvent.leps = cms.InputTag("pfShyftProducer", "muons")
    'Jec11_V3_Uncertainty_AK5PFchs.txt'
]


## Maximal Number of Events
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )

from Analysis.SHyFT.shyftPFSelection_cfi import shyftPFSelection as shyftPFSelectionInput


process.load('Analysis.PdfWeights.pdfWeightProducer_cfi')


process.pfShyftProducer = cms.EDFilter('EDSHyFTPFSelector',
                                    shyftPFSelection = shyftPFSelectionInput.clone(
                                           jetSrc = cms.InputTag('goodPatJetsPFlow'),
                                           jecPayloads = cms.vstring( payloads )
                                        )
                                    )


## std sequence to produce the kinematic fit for semi-leptonic events
process.load('TopQuarkAnalysis.TopKinFitter.TtSemiLepKinFitProducer_Muons_cfi')

process.kinFitTtSemiLepEvent.jets = cms.InputTag("pfShyftProducer", "jets")
process.kinFitTtSemiLepEvent.mets = cms.InputTag("pfShyftProducer", "MET")
process.kinFitTtSemiLepEvent.leps = cms.InputTag("pfShyftProducer", "muons")

process.load("TopQuarkAnalysis.TopKinFitter.TtSemiLepKinFitProducer_Muons_cfi")



## Maximal Number of Events
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )

from Analysis.SHyFT.shyftPFSelection_cfi import shyftPFSelection as shyftPFSelectionInput

process.pileup = cms.EDFilter('PileUpProducer',
                                pvSrc = cms.InputTag('goodOfflinePrimaryVertices')

)
#!!!Change below doElectrons .. False to True
process.pfShyftProducerAK5 = cms.EDFilter('EDSHyFTPFSelector',
                                    shyftPFSelection = shyftPFSelectionInput.clone(
                                           jetSrc = cms.InputTag('goodPatJetsPFlow'),
    					   rhoSrc = cms.InputTag('kt6PFJets', 'rho'),
					   doElectrons = cms.bool(False),
                                           jecPayloads = cms.vstring( payloads )
                                        )
                                    )
#!!!Change below doElectrons from False to True
process.pfShyftProducerAK5Loose = cms.EDFilter('EDSHyFTPFSelector',
                                            shyftPFSelection = shyftPFSelectionInput.clone(
                                                jetSrc = cms.InputTag('goodPatJetsPFlow'),
                                                muonSrc = cms.InputTag('selectedPatMuonsPFlowLoose'),
                                                electronSrc = cms.InputTag('selectedPatElectronsPFlowLoose'),
						doElectrons = cms.bool(False),
    					   	rhoSrc = cms.InputTag('kt6PFJets', 'rho'),
                                                jecPayloads = cms.vstring( payloads ),
                                                removeLooseLep = cms.bool(True)
                                                )
                                            )
        process.hltSelection = cms.Sequence(process.hltSelectionMu)


## Maximum number of events ("-1" means all events)
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(-1))

from Analysis.SHyFT.shyftPFSelection_cfi import shyftPFSelection as shyftPFSelectionInput

process.pileup = cms.EDFilter("PileUpProducer", pvSrc=cms.InputTag("goodOfflinePrimaryVertices"))

## Define various producers
process.pfShyftProducerAK5 = cms.EDFilter(
    "EDSHyFTPFSelector",
    shyftPFSelection=shyftPFSelectionInput.clone(
        jetSrc=cms.InputTag("goodPatJetsPFlow"),
        rhoSrc=cms.InputTag("kt6PFJets", "rho"),
        doElectrons=cms.bool(inputMuOrEle),  # 0 = mu, 1 = ele
        jecPayloads=cms.vstring(payloads),
    ),
)

process.pfShyftProducerAK5Loose = cms.EDFilter(
    "EDSHyFTPFSelector",
    shyftPFSelection=shyftPFSelectionInput.clone(
        jetSrc=cms.InputTag("goodPatJetsPFlow"),
        muonSrc=cms.InputTag("selectedPatMuonsPFlowLoose"),
        electronSrc=cms.InputTag("selectedPatElectronsPFlowLoose"),
        doElectrons=cms.bool(inputMuOrEle),  # 0 = mu, 1 = ele
        rhoSrc=cms.InputTag("kt6PFJets", "rho"),
        jecPayloads=cms.vstring(payloads),
        removeLooseLep=cms.bool(True),
    ),
Beispiel #8
0
    if options.muOrEle:
        process.hltSelection = cms.Sequence(~process.hltSelectionMu *
                                            process.hltSelectionEle)
    else:
        process.hltSelection = cms.Sequence(process.hltSelectionMu)

## Maximal Number of Events
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(-1))

from Analysis.SHyFT.shyftPFSelection_cfi import shyftPFSelection as shyftPFSelectionInput

process.pfShyftProducer = cms.EDFilter(
    'EDSHyFTPFSelector',
    shyftPFSelection=shyftPFSelectionInput.clone(
        jetSrc=cms.InputTag('goodPatJetsPFlow'),
        muonSrc=cms.InputTag('selectedPatMuonsLoosePFlow'),
        removeLooseLep=cms.bool(True),
        jecPayloads=cms.vstring(payloads)))

process.pfShyftProducerPrunedCA8 = cms.EDFilter(
    'EDSHyFTPFSelector',
    shyftPFSelection=shyftPFSelectionInput.clone(
        jetSrc=cms.InputTag('goodPatJetsCA8PrunedPF'),
        jecPayloads=cms.vstring(payloads)))

process.pfShyftProducer.shyftPFSelection.cutsToIgnore.append(
    '== 1 Tight Lepton')
process.pfShyftProducer.shyftPFSelection.cutsToIgnore.append('0 other lepton')
process.pfShyftProducer.shyftPFSelection.cutsToIgnore.append('>=1 Jets')

process.pfShyftProducerPrunedCA8.shyftPFSelection.cutsToIgnore.append(
Beispiel #9
0
        'dcap:///pnfs/cms/WAX/11/store/user/lpctlbsm/srappocc/SingleMu/ttbsm_v8_Run2011-May10ReReco/0d3d9a54f3a29af186ad87df2a0c3ce1/ttbsm_42x_data_120_1_dbd.root',
        'dcap:///pnfs/cms/WAX/11/store/user/lpctlbsm/srappocc/SingleMu/ttbsm_v8_Run2011-May10ReReco/0d3d9a54f3a29af186ad87df2a0c3ce1/ttbsm_42x_data_121_1_kcS.root',
    ))

## Maximal Number of Events
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(-1))

from Analysis.SHyFT.shyftPFSelection_cfi import shyftPFSelection as shyftPFSelectionInput

process.pfShyftProducer = cms.EDFilter(
    'EDSHyFTPFSelector',
    shyftPFSelection=shyftPFSelectionInput.clone(
        muonSrc=cms.InputTag('selectedPatMuonsLoosePFlow'),
        electronSrc=cms.InputTag('selectedPatElectronsLoosePFlow'),
        muonIdPFTight=shyftPFSelectionInput.muonIdPFTight.clone(
            cutsToIgnore=cms.vstring('PFIso')),
        electronIdPFTight=shyftPFSelectionInput.electronIdPFTight.clone(
            cutsToIgnore=cms.vstring('PFIso')),
        removeLooseLep=cms.bool(True)
        # replace things here if you want
    ))

## std sequence to produce the kinematic fit for semi-leptonic events
process.load('TopQuarkAnalysis.TopKinFitter.TtSemiLepKinFitProducer_Muons_cfi')

process.kinFitTtSemiLepEvent.jets = cms.InputTag("pfShyftProducer", "jets")
process.kinFitTtSemiLepEvent.mets = cms.InputTag("pfShyftProducer", "MET")
process.kinFitTtSemiLepEvent.leps = cms.InputTag("pfShyftProducer", "muons")

process.load("TopQuarkAnalysis.TopKinFitter.TtSemiLepKinFitProducer_Muons_cfi")

process.pfShyftTupleJets = cms.EDProducer(

## Maximal Number of Events
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )

from Analysis.SHyFT.shyftPFSelection_cfi import shyftPFSelection as shyftPFSelectionInput

process.pileup = cms.EDFilter('PileUpProducer',
                                pvSrc = cms.InputTag('goodOfflinePrimaryVertices')

)

process.pfShyftProducer = cms.EDFilter('EDSHyFTPFSelector',
                                    shyftPFSelection = shyftPFSelectionInput.clone(
                                           jetSrc = cms.InputTag('goodPatJetsCA8PrunedPFPacked'),
    					   rhoSrc = cms.InputTag('kt6PFJets', 'rho'),
					   doElectrons = cms.bool(False),
                                           jecPayloads = cms.vstring( payloads )
                                        )
                                    )

process.pfShyftProducerLoose = cms.EDFilter('EDSHyFTPFSelector',
                                            shyftPFSelection = shyftPFSelectionInput.clone(
                                                jetSrc = cms.InputTag('goodPatJetsCA8PrunedPFPacked'),
                                                muonSrc = cms.InputTag('selectedPatMuonsPFlowLoose'),
                                                electronSrc = cms.InputTag('selectedPatElectronsPFlowLoose'),
						doElectrons = cms.bool(False),
    					   	rhoSrc = cms.InputTag('kt6PFJets', 'rho'),
                                                jecPayloads = cms.vstring( payloads ),
                                                removeLooseLep = cms.bool(True)
                                                )
                                            )
import FWCore.ParameterSet.Config as cms

from Analysis.BoostedTopAnalysis.BoostedTopWTagParams_cfi import boostedTopWTagParams
from Analysis.SHyFT.shyftPFSelection_cfi import shyftPFSelection


boostedTopMassAnalysis = cms.PSet(
    hadronicWParams = boostedTopWTagParams.clone(),
    shyftSelection = shyftPFSelection.clone()
)