def printify(process):
    process.MessageLogger.categories.append('PrintEvent')

    process.load('HLTrigger.HLTcore.triggerSummaryAnalyzerAOD_cfi')
    process.triggerSummaryAnalyzerAOD.inputTag = cms.InputTag(
        'hltTriggerSummaryAOD', '', 'HLT')
    if hasattr(process, 'pathSimple'):
        process.pathSimple *= process.triggerSummaryAnalyzerAOD

    process.PrintOriginalMuons = cms.EDAnalyzer(
        'PrintEvent',
        muon_src=cms.InputTag('cleanPatMuonsTriggerMatch'),
        trigger_results_src=cms.InputTag('TriggerResults', '', 'HLT'))
    process.pathSimple *= process.PrintOriginalMuons

    pe = process.PrintEventSimple = cms.EDAnalyzer(
        'PrintEvent', dilepton_src=cms.InputTag('SimpleMuonsPlusMuonsMinus'))
    if hasattr(process, 'pathSimple'):
        process.pathSimple *= process.PrintEventSimple

    #- 2011-2012 selection (Nlayers > 8)
    #process.PrintEventOurNew = pe.clone(dilepton_src = cms.InputTag('OurNewMuonsPlusMuonsMinus'))
    #process.PrintEventOurNewSS = pe.clone(dilepton_src = cms.InputTag('OurNewMuonsSameSign'))
    #process.PrintEventOurNewEmu = pe.clone(dilepton_src = cms.InputTag('OurNewMuonsElectronsOppSign'))
    #process.pathOurNew *= process.PrintEventOurNew * process.PrintEventOurNewSS * process.PrintEventOurNewEmu

    #- December 2012 selection (Nlayers > 5, re-tuned TuneP, dpT/pT < 0.3)
    if hasattr(process, 'pathOur2012'):
        process.PrintEventOur2012 = pe.clone(
            dilepton_src=cms.InputTag('Our2012MuonsPlusMuonsMinus'))
        process.PrintEventOur2012SS = pe.clone(
            dilepton_src=cms.InputTag('Our2012MuonsSameSign'))
        process.PrintEventOur2012Emu = pe.clone(
            dilepton_src=cms.InputTag('Our2012MuonsElectronsOppSign'))
        process.pathOur2012 *= process.PrintEventOur2012 * process.PrintEventOur2012SS * process.PrintEventOur2012Emu
예제 #2
0
def ntuplify(process, fill_gen_info=True):

    if miniAOD:
        process.load('SUSYBSMAnalysis.Zprime2muAnalysis.PrunedMCLeptons_cfi')
        obj = process.prunedMCLeptons
        obj.src = cms.InputTag('prunedGenParticles')

        process.SimpleNtupler = cms.EDAnalyzer(
            'SimpleNtupler_miniAOD',
            dimu_src=cms.InputTag('SimpleMuonsAllSigns'),
            met_src=cms.InputTag("slimmedMETs"),
            jet_src=cms.InputTag("slimmedJets"),
            beamspot_src=cms.InputTag('offlineBeamSpot'),
            vertices_src=cms.InputTag('offlineSlimmedPrimaryVertices'),
            TriggerResults_src=cms.InputTag('TriggerResults', '', 'PAT'),  #mc
            # 								TriggerResults_src = cms.InputTag('TriggerResults', '', 'RECO'),	#data
            genEventInfo=cms.untracked.InputTag('generator'),
            metFilter=cms.VInputTag(
                cms.InputTag("Flag_HBHENoiseFilter"),
                cms.InputTag("Flag_HBHENoiseIsoFilter"),
                cms.InputTag("Flag_EcalDeadCellTriggerPrimitiveFilter"),
                cms.InputTag("Flag_eeBadScFilter"),
                cms.InputTag("Flag_globalTightHalo2016Filter")))

    else:
        process.SimpleNtupler = cms.EDAnalyzer(
            'SimpleNtupler',
            dimu_src=cms.InputTag('SimpleMuonsAllSigns'),
            met_src=cms.InputTag("patMETsPF"),
            jet_src=cms.InputTag("cleanPatJets"),
            beamspot_src=cms.InputTag('offlineBeamSpot'),
            vertices_src=cms.InputTag('offlinePrimaryVertices'),
            TriggerResults_src=cms.InputTag('TriggerResults', '', 'PAT'),
            genEventInfo=cms.untracked.InputTag('generator'))
    if miniAOD and Electrons:
        process.SimpleNtuplerEmu = process.SimpleNtupler.clone(
            dimu_src=cms.InputTag('SimpleMuonsElectronsAllSigns'))

    if fill_gen_info:
        from SUSYBSMAnalysis.Zprime2muAnalysis.HardInteraction_cff import hardInteraction
        process.SimpleNtupler.hardInteraction = hardInteraction

    if hasattr(process, 'pathSimple'):
        if miniAOD and fill_gen_info:
            process.pathSimple *= obj * process.SimpleNtupler
            if Electrons:
                process.pathSimple *= obj * process.SimpleNtupler * process.SimpleNtuplerEmu
        else:
            process.pathSimple *= process.SimpleNtupler
            if Electrons:
                process.pathSimple *= process.SimpleNtupler * process.SimpleNtuplerEmu
예제 #3
0
def printify(process):
    process.MessageLogger.categories.append('PrintEvent')

    process.load('HLTrigger.HLTcore.triggerSummaryAnalyzerAOD_cfi')
    process.triggerSummaryAnalyzerAOD.inputTag = cms.InputTag(
        'hltTriggerSummaryAOD', '', 'HLT')
    if hasattr(process, 'pathSimple'):
        process.pathSimple *= process.triggerSummaryAnalyzerAOD

    process.PrintOriginalMuons = cms.EDAnalyzer(
        'PrintEvent',
        muon_src=cms.InputTag('cleanPatMuonsTriggerMatch'),
        trigger_results_src=cms.InputTag('TriggerResults', '', 'HLT'))
    process.pathSimple *= process.PrintOriginalMuons

    pe = process.PrintEventSimple = cms.EDAnalyzer(
        'PrintEvent', dilepton_src=cms.InputTag('SimpleMuonsPlusMuonsMinus'))
    if hasattr(process, 'pathSimple'):
        process.pathSimple *= process.PrintEventSimple
예제 #4
0
def printify(process):
    process.MessageLogger.categories.append('PrintEvent')

    process.load('HLTrigger.HLTcore.triggerSummaryAnalyzerAOD_cfi')
    process.triggerSummaryAnalyzerAOD.inputTag = cms.InputTag('hltTriggerSummaryAOD','','HLT')
    process.pathSimple *= process.triggerSummaryAnalyzerAOD

    process.PrintOriginalMuons = cms.EDAnalyzer('PrintEvent', muon_src = cms.InputTag('cleanPatMuonsTriggerMatch'), trigger_results_src = cms.InputTag('TriggerResults','','HLT'))
    process.pathSimple *= process.PrintOriginalMuons

    pe = process.PrintEventSimple = cms.EDAnalyzer('PrintEvent', dilepton_src = cms.InputTag('SimpleMuonsPlusMuonsMinus'))
    process.pathSimple *= process.PrintEventSimple

    process.PrintEventOurNew = pe.clone(dilepton_src = cms.InputTag('OurNewMuonsPlusMuonsMinus'))
    process.PrintEventOurNewSS = pe.clone(dilepton_src = cms.InputTag('OurNewMuonsSameSign'))
    process.PrintEventOurNewEmu = pe.clone(dilepton_src = cms.InputTag('OurNewMuonsElectronsOppSign'))
    process.pathOurNew *= process.PrintEventOurNew * process.PrintEventOurNewSS * process.PrintEventOurNewEmu

    process.PrintEventVBTF = pe.clone(dilepton_src = cms.InputTag('VBTFMuonsPlusMuonsMinus'))
    process.pathVBTF *= process.PrintEventVBTF
def printify(process):
    process.MessageLogger.categories.append('PrintEvent')

    process.load('HLTrigger.HLTcore.triggerSummaryAnalyzerAOD_cfi')
    process.triggerSummaryAnalyzerAOD.inputTag = cms.InputTag('hltTriggerSummaryAOD','','HLT')
    process.pathSimple *= process.triggerSummaryAnalyzerAOD

    process.PrintOriginalMuons = cms.EDAnalyzer('PrintEvent', muon_src = cms.InputTag('cleanPatMuonsTriggerMatch'), trigger_results_src = cms.InputTag('TriggerResults','','HLT'))
    process.pathSimple *= process.PrintOriginalMuons

    pe = process.PrintEventSimple = cms.EDAnalyzer('PrintEvent', dilepton_src = cms.InputTag('SimpleMuonsPlusMuonsMinus'))
    process.pathSimple *= process.PrintEventSimple

    #- 2011-2012 selection (Nlayers > 8)
    #process.PrintEventOurNew = pe.clone(dilepton_src = cms.InputTag('OurNewMuonsPlusMuonsMinus'))
    #process.PrintEventOurNewSS = pe.clone(dilepton_src = cms.InputTag('OurNewMuonsSameSign'))
    #process.PrintEventOurNewEmu = pe.clone(dilepton_src = cms.InputTag('OurNewMuonsElectronsOppSign'))
    #process.pathOurNew *= process.PrintEventOurNew * process.PrintEventOurNewSS * process.PrintEventOurNewEmu

    #- December 2012 selection (Nlayers > 5, re-tuned TuneP, dpT/pT < 0.3)
    process.PrintEventOur2012    = pe.clone(dilepton_src = cms.InputTag('Our2012MuonsPlusMuonsMinus'))
    process.PrintEventOur2012SS  = pe.clone(dilepton_src = cms.InputTag('Our2012MuonsSameSign'))
    process.PrintEventOur2012Emu = pe.clone(dilepton_src = cms.InputTag('Our2012MuonsElectronsOppSign'))
    process.pathOur2012 *= process.PrintEventOur2012 * process.PrintEventOur2012SS * process.PrintEventOur2012Emu
import sys, os
from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cfg import cms, process

process.source.fileNames = ['/store/data/Run2012C/SingleMu/AOD/PromptReco-v2/000/202/272/94E42528-89F9-E111-BE44-BCAEC53296F8.root']
process.GlobalTag.globaltag = 'GR_P_V42_AN2::All'
process.maxEvents.input = 5000
process.options.wantSummary = True
process.MessageLogger.cerr.FwkReport.reportEvery = 100000

process.load('SUSYBSMAnalysis.Zprime2muAnalysis.CheckPrescale_cfi')
process.CheckPrescale.dump_prescales = True

process.Mu17       = process.CheckPrescale.clone(trigger_paths=cms.vstring('HLT_Mu17_v3'))
process.Mu15eta2p1 = process.CheckPrescale.clone(trigger_paths=cms.vstring('HLT_Mu15_eta2p1_v3', 'HLT_Mu15_eta2p1_v4'))
process.Mu24eta2p1 = process.CheckPrescale.clone(trigger_paths=cms.vstring('HLT_Mu24_eta2p1_v3', 'HLT_Mu24_eta2p1_v4', 'HLT_Mu24_eta2p1_v5'))

process.MessageLogger.suppressWarning = cms.untracked.vstring('Mu17', 'Mu15eta2p1', 'Mu24eta2p1')

#process.p = cms.Path(process.Mu17 * process.Mu15eta2p1 * process.Mu24eta2p1)
process.p = cms.Path(process.Mu24eta2p1)

if __name__ == '__main__' and 'submit' in sys.argv:
    crab_cfg = '''
[CRAB]
jobtype = cmssw
scheduler = condor

[CMSSW]
datasetpath = %(dataset)s
pset = GetPrescales.py
total_number_of_lumis = -1
import FWCore.ParameterSet.Config as cms

from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cfg import process

process.load('SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAsymmetry_cfi')

process.genDimuons = cms.EDProducer('GenDilCandViewShallowCloneCombiner',
                                    decay = cms.string('genParticles@+ genParticles@-'),
                                    cut = cms.string('')
                                    )

process.source.fileNames = ['file:/uscms_data/d1/tucker/pattuplesofeventsabove200gev/2010a.root', 'file:/uscms_data/d1/tucker/pattuplesofeventsabove200gev/2010b.root']
process.p = cms.Path(process.Zprime2muAnalysisSequence * process.Zprime2muAsymmetry)

process.Zprime2muAsymmetry.numFits = 4

process.Zprime2muAsymmetry.calculateMistag = False
process.Zprime2muAsymmetry.useGen = False
#process.Zprime2muAsymmetry.noFit = True

#process.Zprime2muAsymmetry.peakMass = cms.double(91.2)
#process.Zprime2muAsymmetry.massDistType = cms.int32(3),
#process.Zprime2muAsymmetry.maxPt = cms.double(200.0),
#process.Zprime2muAsymmetry.maxRapidity = cms.double(3.5),
#process.Zprime2muAsymmetry.fitWindow = cms.vdouble(80.0, 100.0),
##process.Zprime2muAsymmetry.fitWindow = cms.vdouble(200.0, 500.0),
##   order: cos_cs, rapidity, pT, phi, mass, phi_cs
#process.Zprime2muAsymmetry.recSigma = cms.vdouble(5.8e-4, 8.6e-3, 3.9, 8.6e-2, 6.14, 9.5e-2)


예제 #8
0
from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cfg import cms, process

process.maxEvents.input = -1
process.source.fileNames = [
    '/store/mc/RunIISummer17MiniAOD/ZToMuMu_NNPDF30_13TeV-powheg_M_120_200/MINIAODSIM/NZSFlatPU28to62_92X_upgrade2017_realistic_v10-v1/150000/C08A23CE-A5AC-E711-9142-002590D9D98E.root'
    # '/store/mc/PhaseIFall16MiniAOD/ZToMuMu_NNPDF30_13TeV-powheg_M_200_400/MINIAODSIM/FlatPU28to62HcalNZSRAW_PhaseIFall16_exo52_90X_upgrade2017_realistic_v6_C1-v1/120000/304E419F-CC13-E711-93E9-FA163E0231A1.root'
]
process.options.wantSummary = True

ex = ''

from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cff import leptons, leptonsMini, muonPhotonMatchMiniAOD, muonPhotonMatch, allDimuons, dimuons, rec_level_module
#tracks = ['global', 'inner', 'tpfms', 'picky', 'tunep', 'tmr', 'tunepnew']
tracks = ['tunepnew']

process.load('SUSYBSMAnalysis.Zprime2muAnalysis.HardInteractionFilter_cfi')
process.HardInteractionFilterRes = process.HardInteractionFilter.clone(
    use_resonance_mass=True)

process.load('SUSYBSMAnalysis.Zprime2muAnalysis.EfficiencyFromMC_cfi')

import SUSYBSMAnalysis.Zprime2muAnalysis.VBTFSelection_cff as VBTFSelection
process.allDimuonsVBTF = VBTFSelection.allDimuons.clone()
process.dimuonsVBTF = VBTFSelection.dimuons.clone(src='allDimuonsVBTF')
process.VBTFEfficiencyFromMCMini = process.EfficiencyFromMCMini.clone(
    dimuon_src='dimuonsVBTF', acceptance_max_eta_2=2.4)
process.VBTFEfficiencyFromMC = process.EfficiencyFromMC.clone(
    dimuon_src='dimuonsVBTF', acceptance_max_eta_2=2.4)
process.VBTFEfficiencyFromMCnoTrigger = process.EfficiencyFromMCnoTrigger.clone(
    dimuon_src='dimuonsVBTF', acceptance_max_eta_2=2.4)  ### NO TRIGGER PROCESS
def check_prescale(process, trigger_paths, hlt_process_name='HLT'):
    process.load('SUSYBSMAnalysis.Zprime2muAnalysis.CheckPrescale_cfi')
    process.CheckPrescale.trigger_paths = cms.vstring(*trigger_paths)
    process.pCheckPrescale = cms.Path(process.CheckPrescale)
예제 #10
0
from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cfg import cms, process

process.source.fileNames = ['/store/data/Run2011A/SingleMu/AOD/PromptReco-v4/000/167/913/F2B6FB32-7AA3-E011-BCAF-BCAEC5329710.root']
process.GlobalTag.globaltag = 'GR_R_42_V13::All'
process.maxEvents.input = -1
process.options.wantSummary = True

process.load('SUSYBSMAnalysis.Zprime2muAnalysis.CheckPrescale_cfi')
process.CheckPrescale.trigger_paths = cms.vstring('HLT_Mu30_v1', 'HLT_Mu30_v2', 'HLT_Mu30_v3', 'HLT_Mu30_v4', 'HLT_Mu30_v5')
process.p = cms.Path(process.CheckPrescale)
def check_prescale(process, trigger_paths, hlt_process_name='HLT'):
    process.load('SUSYBSMAnalysis.Zprime2muAnalysis.CheckPrescale_cfi')
    process.CheckPrescale.trigger_paths = cms.vstring(*trigger_paths)
    process.pCheckPrescale = cms.Path(process.CheckPrescale)
#!/usr/bin/env python

import sys, os, FWCore.ParameterSet.Config as cms
from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cff import switch_hlt_process_name
from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cfg import process
process.getPrescales = cms.EDProducer('GetPrescaleL1HLT',
    hlt_src = cms.InputTag('TriggerResults','','HLT'),
)       
process.insertPrescales= cms.Path(process.getPrescales)
from SUSYBSMAnalysis.Zprime2muAnalysis.hltTriggerMatch_cfi import trigger_match, prescaled_trigger_match, trigger_paths, prescaled_trigger_paths, overall_prescale, offline_pt_threshold, prescaled_offline_pt_threshold

# Since the prescaled trigger comes with different prescales in
# different runs/lumis, this filter prescales it to a common factor to
# make things simpler.
process.load('SUSYBSMAnalysis.Zprime2muAnalysis.PrescaleToCommon_cff')
process.PrescaleToCommon.trigger_paths = prescaled_trigger_paths
process.PrescaleToCommon.overall_prescale = overall_prescale
process.PrescaleToCommon.disable = cms.untracked.bool(True)

# The histogramming module that will be cloned multiple times below
# for making histograms with different cut/dilepton combinations.
from SUSYBSMAnalysis.Zprime2muAnalysis.HistosFromPAT_cfi import HistosFromPAT
HistosFromPAT.leptonsFromDileptons = True
HistosFromPAT.usePrescaleWeight= cms.untracked.bool(True)

# These modules define the basic selection cuts. For the monitoring
# sets below, we don't need to define a whole new module, since they
# just change one or two cuts -- see below.
import SUSYBSMAnalysis.Zprime2muAnalysis.VBTFSelection_cff as VBTFSelection
import SUSYBSMAnalysis.Zprime2muAnalysis.OurSelectionOld_cff as OurSelectionOld
import SUSYBSMAnalysis.Zprime2muAnalysis.OurSelection2011EPS_cff as OurSelection2011EPS
#process.source.fileNames = ['file:crab1/crab_effres_zp1000/res/pat_2_1_375.root']
process.source.fileNames = ['/store/user/slava/DYToMuMu_M_20_TuneZ2star_8TeV_pythia6/effres_dy20/20941d9c676d6826327c8223aa3d20e0/pat_9_1_5zY.root']
process.options.wantSummary = True

ex = ''

if use_old_selection:
    ex += 'oldsel'
    from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cff import switch_to_old_selection
    switch_to_old_selection(process)

from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cff import rec_levels, rec_level_module
tracks = ['global', 'inner', 'tpfms', 'picky', 'tunep', 'tmr', 'tunepnew']
rec_levels(process, tracks)

process.load('SUSYBSMAnalysis.Zprime2muAnalysis.HardInteractionFilter_cfi')
process.HardInteractionFilterRes = process.HardInteractionFilter.clone(use_resonance_mass=True)

process.load('SUSYBSMAnalysis.Zprime2muAnalysis.EfficiencyFromMC_cfi')

#process.HardInteractionFilter.use_resonance_mass = True
#process.EfficiencyFromMC.use_resonance_mass_denom = True

# Since LooseTightPairSelector ignores the cutFor that
# Zprime2muLeptonProducer sets, don't need to redo leptons for the
# VBTF path.
import SUSYBSMAnalysis.Zprime2muAnalysis.VBTFSelection_cff as VBTFSelection
process.allDimuonsVBTF = VBTFSelection.allDimuons.clone()
process.dimuonsVBTF = VBTFSelection.dimuons.clone(src = 'allDimuonsVBTF')
process.VBTFEfficiencyFromMC = process.EfficiencyFromMC.clone(dimuon_src = 'dimuonsVBTF', acceptance_max_eta_2 = 2.1)
예제 #14
0
    if "Run2016H" in fileName:
        process.GlobalTag.globaltag = '80X_dataRun2_Prompt_v14'
    elif "Run2016" in fileName:
        process.GlobalTag.globaltag = '80X_dataRun2_2016SeptRepro_v6'
    else:
        # https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideFrontierConditions#Global_Tags_for_PdmVMCcampai_AN3
        process.GlobalTag.globaltag = '80X_mcRun2_asymptotic_2016_TrancheIV_v8'
        isMC = True
process.MessageLogger.cerr.FwkReport.reportEvery = 1000  # default 1000

from SUSYBSMAnalysis.Zprime2muAnalysis.hltTriggerMatch_cfi import trigger_match, prescaled_trigger_match, trigger_paths, prescaled_trigger_paths, overall_prescale, offline_pt_threshold, prescaled_offline_pt_threshold

# Since the prescaled trigger comes with different prescales in
# different runs/lumis, this filter prescales it to a common factor to
# make things simpler.
process.load('SUSYBSMAnalysis.Zprime2muAnalysis.PrescaleToCommon_cff')
process.PrescaleToCommon.trigger_paths = prescaled_trigger_paths
process.PrescaleToCommon.overall_prescale = overall_prescale

process.PrescaleToCommonMiniAOD.trigger_paths = prescaled_trigger_paths
process.PrescaleToCommonMiniAOD.overall_prescale = overall_prescale

# The histogramming module that will be cloned multiple times below
# for making histograms with different cut/dilepton combinations.

from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cff import electrons_miniAOD
electrons_miniAOD(process)

from SUSYBSMAnalysis.Zprime2muAnalysis.HistosFromPAT_cfi import HistosFromPAT_MiniAOD as HistosFromPAT
HistosFromPAT.leptonsFromDileptons = True
HistosFromPAT.usekFactor = False  #### Set TRUE to use K Factor on DY. If used, the k factor will be applied to ALL samples submitted. #####
예제 #15
0
import FWCore.ParameterSet.Config as cms

from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cfg import process

process.load('SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAsymmetry_cfi')

process.genDimuons = cms.EDProducer(
    'GenDilCandViewShallowCloneCombiner',
    decay=cms.string('genParticles@+ genParticles@-'),
    cut=cms.string(''))

process.source.fileNames = [
    'file:/uscms_data/d1/tucker/pattuplesofeventsabove200gev/2010a.root',
    'file:/uscms_data/d1/tucker/pattuplesofeventsabove200gev/2010b.root'
]
process.p = cms.Path(process.Zprime2muAnalysisSequence *
                     process.Zprime2muAsymmetry)

process.Zprime2muAsymmetry.numFits = 4

process.Zprime2muAsymmetry.calculateMistag = False
process.Zprime2muAsymmetry.useGen = False
#process.Zprime2muAsymmetry.noFit = True

#process.Zprime2muAsymmetry.peakMass = cms.double(91.2)
#process.Zprime2muAsymmetry.massDistType = cms.int32(3),
#process.Zprime2muAsymmetry.maxPt = cms.double(200.0),
#process.Zprime2muAsymmetry.maxRapidity = cms.double(3.5),
#process.Zprime2muAsymmetry.fitWindow = cms.vdouble(80.0, 100.0),
##process.Zprime2muAsymmetry.fitWindow = cms.vdouble(200.0, 500.0),
##   order: cos_cs, rapidity, pT, phi, mass, phi_cs
예제 #16
0
    'file:/u/user/msoh/Samples/MC/DY/M6000/B86337F0-16C7-E611-BA6A-24BE05C44B91.root',
    'file:/u/user/msoh/Samples/MC/DY/M6000/EA325C83-17C7-E611-AFA4-001E67E5E8B6.root',
    'file:/u/user/msoh/Samples/MC/DY/M6000/F214122A-17C7-E611-9134-001E674FC800.root',

    #'file:/u/user/msoh/Data/MC/E412A19B-DCCF-E611-8050-0CC47A546E5E.root',  ##Moriond17 M50-120
    #'file:/u/user/msoh/Data/MC/824C363B-0AC8-E611-B4A5-20CF3027A580.root', ##Moriond17
    #'file:/u/user/msoh/Data/MC/52FEE967-17C7-E611-9A69-FA163E2EED86.root', ##Moriond17 M4500-6000

    #'file:/u/user/msoh/data/mc/m50_120/1A6B76DF-153B-E611-BEC5-0CC47A4DEDF8.root',
]
process.options.wantSummary = True
process.MessageLogger.cerr.FwkReport.reportEvery = 1000  # default 1000

from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cff import leptons, leptonsMini, muonPhotonMatchMiniAOD, muonPhotonMatch, allDimuons, dimuons, rec_level_module

process.load('SUSYBSMAnalysis.Zprime2muAnalysis.HardInteractionFilter_cfi')
process.HardInteractionFilterRes = process.HardInteractionFilter.clone(
    use_resonance_mass=True)

process.load('SUSYBSMAnalysis.Zprime2muAnalysis.EfficiencyFromMC_cfi')

import SUSYBSMAnalysis.Zprime2muAnalysis.OurSelection2016_cff as OurSelection
process.allDimuonsOur = OurSelection.allDimuons.clone()
process.dimuonsOur = OurSelection.dimuons.clone(src='allDimuonsOur')
process.OurEfficiencyFromMCMini = process.EfficiencyFromMCMini.clone(
    dimuon_src='dimuonsOur', acceptance_max_eta_2=2.4)
process.OurEfficiencyFromMC = process.EfficiencyFromMC.clone(
    dimuon_src='dimuonsOur', acceptance_max_eta_2=2.4)
process.OurEfficiencyFromMCnoTrigger = process.EfficiencyFromMCnoTrigger.clone(
    dimuon_src='dimuonsOur', acceptance_max_eta_2=2.4)  ### NO TRIGGER PROCESS
예제 #17
0
        # '/store/mc/RunIISummer16MiniAODv2/WWTo2L2Nu_13TeV-powheg/MINIAODSIM/PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/80000/08E155A9-FAB6-E611-92BF-00259073E45E.root',
        # '/store/mc/RunIISummer16MiniAODv2/WWTo2L2Nu_13TeV-powheg/MINIAODSIM/PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/80000/0CB16E4E-F0B6-E611-9D13-0090FAA58294.root',
        # '/store/mc/RunIISummer16MiniAODv2/WWTo2L2Nu_13TeV-powheg/MINIAODSIM/PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/80000/1AA5E3E6-F4B6-E611-8A54-0090FAA575E0.root',
        #                                                                 '/store/data/Run2016B/SingleMuon/MINIAOD/23Sep2016-v3/00000/162AD1DB-1E98-E611-9893-008CFA56D58C.root',
        #'/store/data/Run2016B/SingleMuon/MINIAOD/23Sep2016-v3/00000/1A1F07FF-2698-E611-915C-0242AC130004.root'
        #                                                                '/store/mc/RunIISpring16MiniAODv2/ZToMuMu_NNPDF30_13TeV-powheg_M_120_200/MINIAODSIM/PUSpring16RAWAODSIM_reHLT_80X_mcRun2_asymptotic_v14-v1/90000/18C80393-613A-E611-86DF-0090FAA573E0.root'
    ),
    secondaryFileNames=secFiles)

secFiles.extend([])

process.maxEvents.input = -1
process.GlobalTag.globaltag = '80X_dataRun2_2016SeptRepro_v6'
#process.MessageLogger.cerr.FwkReport.reportEvery = 1 # default 1000

process.load('SUSYBSMAnalysis.Zprime2muAnalysis.PrunedMCLeptons_cfi')
process.DYGenMassFilter = cms.EDFilter(
    'DibosonGenMass',
    src=cms.InputTag('prunedGenParticles'),
    min_mass=cms.double(50),
    max_mass=cms.double(200),
)

# Define the numerators and denominators, removing cuts from the
# allDimuons maker. "NoX" means remove cut X entirely (i.e. the
# loose_cut denominators), "TiX" means move cut X from the loose_cut
# to the tight_cut (meaning only one muon instead of two has to pass
# the cut).  "NoNo" means remove nothing (i.e. the numerator). This
# will break if loose_, tight_cut strings are changed upstream, so we
# try to check those with a simple string test below.
예제 #18
0
#!/usr/bin/env python

import sys, os, FWCore.ParameterSet.Config as cms
from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cff import switch_hlt_process_name
from SUSYBSMAnalysis.Zprime2muAnalysis.Zprime2muAnalysis_cfg import process
#process.source.fileNames = ['/store/user/slava/DYToMuMu_M-2000_CT10_TuneZ2star_8TeV-powheg-pythia6/datamc_dy2000/ecac376f8fa7ccc229aaa06d757d785a/pat_1_1_G72.root']
#process.maxEvents.input = 100
from SUSYBSMAnalysis.Zprime2muAnalysis.hltTriggerMatch_cfi import trigger_match, prescaled_trigger_match, trigger_paths, prescaled_trigger_paths, overall_prescale, offline_pt_threshold, prescaled_offline_pt_threshold

# Since the prescaled trigger comes with different prescales in
# different runs/lumis, this filter prescales it to a common factor to
# make things simpler.
process.load('SUSYBSMAnalysis.Zprime2muAnalysis.PrescaleToCommon_cff')
process.PrescaleToCommon.trigger_paths = prescaled_trigger_paths
process.PrescaleToCommon.overall_prescale = overall_prescale

# The histogramming module that will be cloned multiple times below
# for making histograms with different cut/dilepton combinations.
from SUSYBSMAnalysis.Zprime2muAnalysis.HistosFromPAT_cfi import HistosFromPAT
HistosFromPAT.leptonsFromDileptons = True

# These modules define the basic selection cuts. For the monitoring
# sets below, we don't need to define a whole new module, since they
# just change one or two cuts -- see below.
#import SUSYBSMAnalysis.Zprime2muAnalysis.VBTFSelection_cff as VBTFSelection
#import SUSYBSMAnalysis.Zprime2muAnalysis.OurSelectionOld_cff as OurSelectionOld
#import SUSYBSMAnalysis.Zprime2muAnalysis.OurSelection2011EPS_cff as OurSelection2011EPS
import SUSYBSMAnalysis.Zprime2muAnalysis.OurSelectionNew_cff as OurSelectionNew
import SUSYBSMAnalysis.Zprime2muAnalysis.OurSelectionDec2012_cff as OurSelectionDec2012

# CandCombiner includes charge-conjugate decays with no way to turn it