# Output file
process.TFileService = cms.Service("TFileService",
        fileName = cms.string("test_combined_mc.root")
        )

# Import rho value for isolation correction
from ZFinder.Event.kt6_pfjets_cfi import kt6PFJetsForIsolation
process.kt6PFJetsForIsolation = kt6PFJetsForIsolation.clone()

# Particle flow isolation
from CommonTools.ParticleFlow.Tools.pfIsolation import setupPFElectronIso
process.eleIsoSequence = setupPFElectronIso(process, 'gsfElectrons')
process.pfiso = cms.Sequence(process.pfParticleSelectionSequence + process.eleIsoSequence)

# ZFinder
from ZFinder.Event.zdefinitions_cfi import zdefs_combined_mc
from ZFinder.Event.zfinder_cfi import ZFinder
process.ZFinder = ZFinder.clone(
        ZDefinitions = zdefs_combined_mc,
        is_mc = cms.bool(True),
        use_muon_acceptance = cms.bool(True),
        require_gen_z = cms.bool(False),
        )

# RUN
process.p = cms.Path(
        process.kt6PFJetsForIsolation
        * process.pfiso
        * process.ZFinder
        )
Ejemplo n.º 2
0
    "FSRWeightProducer",
    GenTag=cms.untracked.InputTag("genParticles"),
)

# ZFinder
from ZFinder.Event.zdefinitions_cfi import zdefs_combined_mc
from ZFinder.Event.zfinder_cfi import ZFinder

process.ZFinder = ZFinder.clone(
    # Use the calibrated electrons we make with process.CalibratedElectrons
    ecalElectronsInputTag=cms.InputTag("CalibratedElectrons",
                                       "calibratedGsfElectrons"),
    ZDefinitions=zdefs_combined_mc,
    is_mc=cms.bool(True),
    use_muon_acceptance=cms.bool(True),
    require_gen_z=cms.bool(False),
    run_pdf_weights=cms.bool(True),
    run_fsr_weight=cms.bool(True),
    gen_electrons=cms.string("Born"),
    # Pick the GEN level electrons to use. Select one of:
    #    - "Dressed"
    #    - "Born"
    #    - "Bare" or "Naked"
)

# RUN
process.p = cms.Path(
    process.pdfWeights * process.fsrWeight * process.kt6PFJetsForIsolation *
    process.eleRegressionEnergy * process.CalibratedElectrons * process.pfiso *
    process.ZFinder)
# Compute FSR weights for uncertainty calculation
process.fsrWeight = cms.EDProducer("FSRWeightProducer",
        GenTag = cms.untracked.InputTag("genParticles"),
        )

# ZFinder
from ZFinder.Event.zdefinitions_cfi import zdefs_combined_mc
from ZFinder.Event.zfinder_cfi import ZFinder
process.ZFinder = ZFinder.clone(
        # Use the calibrated electrons we make with process.CalibratedElectrons
        ecalElectronsInputTag = cms.InputTag("CalibratedElectrons", "calibratedGsfElectrons"),
        ZDefinitions = zdefs_combined_mc,
        is_mc = cms.bool(True),
        use_muon_acceptance = cms.bool(True),
        require_gen_z = cms.bool(False),
        run_pdf_weights = cms.bool(True),
        run_fsr_weight = cms.bool(True),
        gen_electrons = cms.string("Naked"),
	# Pick the GEN level electrons to use. Select one of:
        #    - "Dressed"
        #    - "Born"
        #    - "Bare" or "Naked"
        )

# RUN
process.p = cms.Path(
        process.pdfWeights
        *process.fsrWeight
        *process.kt6PFJetsForIsolation
        * process.eleRegressionEnergy
        * process.CalibratedElectrons
# Import rho value for isolation correction
from ZFinder.Event.kt6_pfjets_cfi import kt6PFJetsForIsolation
process.kt6PFJetsForIsolation = kt6PFJetsForIsolation.clone()

# Particle flow isolation
from CommonTools.ParticleFlow.Tools.pfIsolation import setupPFElectronIso
process.eleIsoSequence = setupPFElectronIso(process, 'CalibratedElectrons:calibratedGsfElectrons')
process.pfiso = cms.Sequence(process.pfParticleSelectionSequence + process.eleIsoSequence)

# ZFinder
from ZFinder.Event.zdefinitions_cfi import zdefs_combined_data
from ZFinder.Event.zfinder_cfi import ZFinder

process.ZFinder = ZFinder.clone(
        # Use the calibrated electrons we make with process.CalibratedElectrons
        ecalElectronsInputTag = cms.InputTag("CalibratedElectrons", "calibratedGsfElectrons"),
        ZDefinitions = zdefs_combined_data,
        use_muon_acceptance = cms.bool(True),

        )

# RUN
process.p = cms.Path(
        process.kt6PFJetsForIsolation
        * process.eleRegressionEnergy
        * process.CalibratedElectrons
        * process.pfiso
        * process.ZFinder
        )
from ZFinder.Event.electron_regression_cfi import CalibratedElectrons, RandomNumberGeneratorService, ElectronEnergyRegressions
process.RandomNumberGeneratorService = RandomNumberGeneratorService
process.CalibratedElectrons = CalibratedElectrons
process.eleRegressionEnergy = ElectronEnergyRegressions

# Import rho value for isolation correction
from ZFinder.Event.kt6_pfjets_cfi import kt6PFJetsForIsolation
process.kt6PFJetsForIsolation = kt6PFJetsForIsolation.clone()

# Particle flow isolation
from CommonTools.ParticleFlow.Tools.pfIsolation import setupPFElectronIso
process.eleIsoSequence = setupPFElectronIso(process, 'CalibratedElectrons:calibratedGsfElectrons')
process.pfiso = cms.Sequence(process.pfParticleSelectionSequence + process.eleIsoSequence)

# ZFinder
from ZFinder.Event.zfinder_cfi import ZFinder
process.ZFinder = ZFinder.clone(
        # Use the calibrated electrons we make with
        # process.CalibratedElectrons
        ecalElectronsInputTag = cms.InputTag("CalibratedElectrons", "calibratedGsfElectrons"),
        )

# RUN
process.p = cms.Path(
        process.kt6PFJetsForIsolation
        * process.eleRegressionEnergy
        * process.CalibratedElectrons
        * process.pfiso
        * process.ZFinder
        )
Ejemplo n.º 6
0
process.CalibratedElectrons = CalibratedElectrons
process.eleRegressionEnergy = ElectronEnergyRegressions

# Import rho value for isolation correction
from ZFinder.Event.kt6_pfjets_cfi import kt6PFJetsForIsolation
process.kt6PFJetsForIsolation = kt6PFJetsForIsolation.clone()

# Particle flow isolation
from CommonTools.ParticleFlow.Tools.pfIsolation import setupPFElectronIso
process.eleIsoSequence = setupPFElectronIso(
    process, 'CalibratedElectrons:calibratedGsfElectrons')
process.pfiso = cms.Sequence(process.pfParticleSelectionSequence +
                             process.eleIsoSequence)

# ZFinder
from ZFinder.Event.zdefinitions_cfi import zdefs_combined_data
from ZFinder.Event.zfinder_cfi import ZFinder

process.ZFinder = ZFinder.clone(
    # Use the calibrated electrons we make with process.CalibratedElectrons
    ecalElectronsInputTag=cms.InputTag("CalibratedElectrons",
                                       "calibratedGsfElectrons"),
    ZDefinitions=zdefs_combined_data,
    use_muon_acceptance=cms.bool(True),
)

# RUN
process.p = cms.Path(
    process.kt6PFJetsForIsolation * process.eleRegressionEnergy *
    process.CalibratedElectrons * process.pfiso * process.ZFinder)
Ejemplo n.º 7
0
        )

# Output file
process.TFileService = cms.Service("TFileService",
        fileName = cms.string("test_extended_mc.root")
        )

# Import rho value for isolation correction
from ZFinder.Event.kt6_pfjets_cfi import kt6PFJetsForIsolation
process.kt6PFJetsForIsolation = kt6PFJetsForIsolation.clone()

# Particle flow isolation
from CommonTools.ParticleFlow.Tools.pfIsolation import setupPFElectronIso
process.eleIsoSequence = setupPFElectronIso(process, 'gsfElectrons')
process.pfiso = cms.Sequence(process.pfParticleSelectionSequence + process.eleIsoSequence)

# ZFinder
from ZFinder.Event.zdefinitions_cfi import zdefs_extended_mc
from ZFinder.Event.zfinder_cfi import ZFinder
process.ZFinder = ZFinder.clone(
        ZDefinitions = zdefs_extended_mc,
        is_mc = cms.bool(True),
        )

# RUN
process.p = cms.Path(
        process.kt6PFJetsForIsolation
        * process.pfiso
        * process.ZFinder
        )
        )

# Run only on lumis specified in the lumi file
# Recipe from:
# https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuidePythonTips#Use_a_JSON_file_of_good_lumi_sec
from FWCore.ParameterSet.Types import untracked, VLuminosityBlockRange
from FWCore.PythonUtilities.LumiList import LumiList
json_file = "../../Metadata/lumi_json/Run2012ABCD.json"  # File location
run_2012abcd_lumis = LumiList(filename = json_file).getCMSSWString().split(',')
process.source.lumisToProcess = untracked(VLuminosityBlockRange(run_2012abcd_lumis))

# Import rho value for isolation correction
from ZFinder.Event.kt6_pfjets_cfi import kt6PFJetsForIsolation
process.kt6PFJetsForIsolation = kt6PFJetsForIsolation.clone()

# Particle flow isolation
from CommonTools.ParticleFlow.Tools.pfIsolation import setupPFElectronIso
process.eleIsoSequence = setupPFElectronIso(process, 'gsfElectrons')
process.pfiso = cms.Sequence(process.pfParticleSelectionSequence + process.eleIsoSequence)

# ZFinder
from ZFinder.Event.zfinder_cfi import ZFinder
process.ZFinder = ZFinder.clone()

# RUN
process.p = cms.Path(
        process.kt6PFJetsForIsolation
        * process.pfiso
        * process.ZFinder
        )
# https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuidePythonTips#Use_a_JSON_file_of_good_lumi_sec
from FWCore.ParameterSet.Types import untracked, VLuminosityBlockRange
from FWCore.PythonUtilities.LumiList import LumiList
json_file = "../../Metadata/lumi_json/Run2012ABCD.json"  # File location
run_2012abcd_lumis = LumiList(filename = json_file).getCMSSWString().split(',')
process.source.lumisToProcess = untracked(VLuminosityBlockRange(run_2012abcd_lumis))

# Import rho value for isolation correction
from ZFinder.Event.kt6_pfjets_cfi import kt6PFJetsForIsolation
process.kt6PFJetsForIsolation = kt6PFJetsForIsolation.clone()

# Particle flow isolation
from CommonTools.ParticleFlow.Tools.pfIsolation import setupPFElectronIso
process.eleIsoSequence = setupPFElectronIso(process, 'gsfElectrons')
process.pfiso = cms.Sequence(process.pfParticleSelectionSequence + process.eleIsoSequence)

# ZFinder
from ZFinder.Event.zdefinitions_cfi import zdefs_combined_data
from ZFinder.Event.zfinder_cfi import ZFinder
process.ZFinder = ZFinder.clone(
        ZDefinitions = zdefs_combined_data,
        use_muon_acceptance = cms.bool(True),
        )

# RUN
process.p = cms.Path(
        process.kt6PFJetsForIsolation
        * process.pfiso
        * process.ZFinder
        )