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
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
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
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)
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)
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)
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. #####
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
'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
# '/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.
#!/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