def parse_args():
options = VarParsing("analysis")
options.register(
'vtxsmear',
'NoSmear',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Vertex smearing options ('Flat' or 'NoSmear')",
)
options.register(
'genseed',
456789,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Generator initial seed",
)
options.register(
'vtxseed',
98765432,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Vertex smear initial seed",
)
options.maxEvents = 10000
options.outputFile = "single_neutrino.root"
options.parseArguments()
return options
def parse_options():
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing("analysis")
# Shortcuts for options
register_float = lambda name, default, help="": options.register(
name, default, VarParsing.multiplicity.singleton, VarParsing.varType.float, help
)
register_bool = lambda name, default, help="": options.register(
name, default, VarParsing.multiplicity.singleton, VarParsing.varType.bool, help
)
register_int = lambda name, default, help="": options.register(
name, default, VarParsing.multiplicity.singleton, VarParsing.varType.int, help
)
# Actual options
register_float('pt', 100., 'Pt in GeV (only for single particle guns)')
register_float('minpt', -1., 'Pt in GeV (only for single particle guns)')
register_float('maxpt', -1., 'Pt in GeV (only for single particle guns)')
register_int('pdgid', -13, 'PdgID of the single particle gun (default -13 muon)')
register_int('seed', 1001, 'RNG seed')
register_bool('debug', False, 'Enable debug logging')
register_bool('dofinecalo', True, 'Enables fine calo for HGCAL')
register_bool('ntuple', False, 'Runs ntupler instead of saving SIM level')
register_bool('profiling', False, 'Enables memory and timing profiling')
options.maxEvents = 1
options.outputFile = 'default.root'
options.parseArguments()
return options
def options() :
options = VP('standard')
options.output = "topTuple.root"
options.maxEvents = 100
options.register('isData', default = True, mytype = VP.varType.bool)
options.register('isFastSim', default = False, mytype = VP.varType.bool)
options.register('quiet', default = True, mytype = VP.varType.bool)
options.register('requireLepton', default = True, mytype = VP.varType.bool)
options.register('globalTag', mytype = VP.varType.string )
options.register('postfix','TR', mytype = VP.varType.string )
options.register('btags', mytype = VP.varType.string, mult = VP.multiplicity.list )
options.register('doElectronEA', default = True, mytype = VP.varType.bool)
options.parseArguments()
options._tagOrder =[]
defaultGT = ('FT53_V21A_AN6' if options.isData else 'START53_V15')
sync53 = '/store/mc/Summer12_DR53X/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola/AODSIM/PU_S10_START53_V7A-v1/0000/FE4C2F81-D0E1-E111-9080-0030487E0A2D.root'
SingleElRun2012A = "root://cmsxrootd.fnal.gov///store/data/Run2012A/SingleElectron/AOD/13Jul2012-v1/0000/001A2EB8-47D4-E111-B527-003048679070.root"
defaultFiles = ( [sync53] if not options.isData else [SingleElRun2012A] )
# pickRelValInputFiles( cmsswVersion = 'CMSSW_5_3_6',
# dataTier = 'RECO',
# relVal = 'SingleMu',
# globalTag = 'GR_R_53_V15_RelVal_mu2012A',
# numberOfFiles = 10
# )
options.files = options.files if options.files else defaultFiles
if not options.globalTag : options.globalTag = defaultGT
if not options.quiet : print options
options.btags = ['combinedSecondaryVertex','jetProbability']
return options
def options() :
options = VP('standard')
options.output = "topTuple.root"
options.maxEvents = -1
options.register('procName', default = "PAT2", mytype = VP.varType.string)
options.register('isData', default = False, mytype = VP.varType.bool)
options.register('skim', default = True, mytype = VP.varType.bool)
options.register('btag', default = False, mytype = VP.varType.bool)
# options.register('outputModule', default = True, mytype = VP.varType.bool)
options.register('quiet', default = False, mytype = VP.varType.bool)
options.register('requireLepton', default = True, mytype = VP.varType.bool)
options.register('globalTag', mytype = VP.varType.string )
options.register('postfix','TR', mytype = VP.varType.string )
options.register('btags', mytype = VP.varType.string, mult = VP.multiplicity.list )
options.register('doElectronEA', default = True, mytype = VP.varType.bool)
options.register('noJetSmearing', default = False, mytype = VP.varType.bool)
try:
for key in options._register.keys():
setattr(options, key, cms_var.get(key, getattr(options, key)))
except NameError:
print "cms_var is not in __builtin__"
print sys.argv
if not "-create" in sys.argv:
options.parseArguments()
else:
print "Found '-create' in sys.argv! Omitting options.parseArguments()"
options._tagOrder =[]
defaultGT = ('GR_R_53_V18' if options.isData else 'START53_V15')
# FT_53_V21_AN4
# START53_V23
sync53 = 'file:/user/kuessel/files/Synchfiles/syncExercise53.root'
diffStep1="file:/user/kuessel/CMSSW/Synch/CMSSW_5_3_8_patch3/src/pickevents.root"
subset="file:/user/kuessel/CMSSW/Synch/CMSSW_5_3_8_patch3/src/testSkim.root"
skim="/store/user/kuessel/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola/YK_MC_MARCH13_TTbar/accebdcbeca2810af478d8af2493d41f/SynchSelMuJets_103_1_lgS.root"
defaultFiles = [skim]
# defaultFiles = [subset]
# defaultFiles = [diffStep1]
options.files = defaultFiles
if not options.globalTag :
options.globalTag = defaultGT
options.btags = ['combinedSecondaryVertex','jetProbability']
if not options.quiet : print options
return options
def options() :
from FWCore.ParameterSet.VarParsing import VarParsing as VP
options = VP('standard')
options.output = "DJ_Tree.root"
options.maxEvents = 100
options.register('GlobalTag', mytype = VP.varType.string)
options.register('quiet', default = True )
options.register('isData', default = True )
options.register('signal', default = False )
options.register('jetCorrections', default = ['L1FastJet','L2Relative','L3Absolute','L2L3Residual'], #L2L3Residual removed from options for simulation (below)
info = "jet correction levels to apply", mult = VP.multiplicity.list, mytype = VP.varType.string)
__hack_ListVarparsingBug__( options, 'jetCorrections')
options.parseArguments()
options._tagOrder =[] # weird, but something to do with options.output
defaultGT,defaultFile = {
"44X" : [('START44_V9B::All','/store/mc/Fall11/QCD_Pt-120to170_TuneZ2_7TeV_pythia6/AODSIM/PU_S6_START42_V14B-v1/0000/0A8E5F68-73F4-E011-AB1C-00E08178C0D5.root'),
('START44_V9B::All','/store/mc/Fall11/HTo2LongLivedTo4F_MH-1000_MFF-350_CTau-350_7TeV-pythia6/GEN-SIM-RECODEBUG/DEBUG-PU_S6_START44_V9B-v4/0000/E83B1DB8-158B-E111-9DB4-0017A477003C.root'),
('FT_R_44_V11::All','/store/data/Run2011B/HT/RECO/19Nov2011-v1/10001/B2E314BC-4746-E111-9BC1-00261894390A.root')],
"53X" : [
#('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-1000_MFF-350_CTau35To3500_8TeV-pythia6/AODSIM/DEBUG_PU_S10_START53_V7A-v1/0000/9C61CA50-0BF3-E111-91E1-003048D2BB22.root'),
#('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-400_MFF-50_CTau8To800_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/00000/80A0705C-7518-E211-A17D-0030487E4D11.root'),
#('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-200_MFF-50_CTau20To2000_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/00000/7A2CB2BB-6C18-E211-AB1B-003048C69416.root'),
('START53_V7F::All','/store/mc/Summer12_DR53X/LongLivedChi0ToMuQQ_MSquark-1500_MChi-494_TuneZ2Star_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/00000/B69AB21F-D50F-E211-B2A7-00215E21DCA2.root'),
('START53_V7F::All','/store/mc/Summer12_DR53X/QCD_Pt-170to300_TuneZ2star_8TeV_pythia6/AODSIM/PU_S10_START53_V7A-v2/00000/0CECB336-B70E-E211-B22D-0018F3D096CE.root'),
('FT_53_V6_AN2::All','/store/data/Run2012B/JetHT/AOD/13Jul2012-v1/00000/04989369-BFD2-E111-B190-002618943933.root'),
('FT_53_V10_AN2::All','/store/data/Run2012C/JetHT/AOD/24Aug2012-v2/00000/22AFA398-D600-E211-AD2E-00261894385A.root'),
('GR_P_V41_AN2::All','/store/data/Run2012C/JetHT/AOD/PromptReco-v2/000/198/955/144C3A0C-D8CF-E111-8A6C-001D09F27067.root'),
('GR_P_V42_AN2::All','/store/data/Run2012D/JetHT/AOD/PromptReco-v1/000/204/506/F8F3A290-8C10-E211-B2D3-001D09F2525D.root'),],
"53XCtrl" : [('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-1000_MFF-350_CTau35To3500_8TeV-pythia6/AODSIM/DEBUG_PU_S10_START53_V7A-v1/0000/9C61CA50-0BF3-E111-91E1-003048D2BB22.root'),
('START53_V7F::All','/store/mc/Summer12_DR53X/QCD_Pt-170to300_TuneZ2star_8TeV_pythia6/AODSIM/PU_S10_START53_V7A-v2/00000/0CECB336-B70E-E211-B22D-0018F3D096CE.root'),
('FT_53_V6_AN2::All','/store/data/Run2012B/SingleMu/AOD/13Jul2012-v1/0000/24DF4BD2-0DD4-E111-BB00-E0CB4E4408DF.root'),
('FT_53_V10_AN2::All','/store/data/Run2012C/SingleMu/AOD/24Aug2012-v1/00000/10FD4F26-FAEE-E111-9AB5-00259073E4FC.root'),
('GR_P_V41_AN2::All','/store/data/Run2012C/SingleMu/AOD/PromptReco-v2/000/198/955/247727CE-98CF-E111-B2EF-5404A63886A0.root'),
('GR_P_V42_AN2::All','/store/data/Run2012D/SingleMu/AOD/PromptReco-v1/000/204/506/8078B6DD-9210-E211-80BE-BCAEC5329700.root'),]
}["53X"][int(not options.signal)+ options.isData]
options.files = options.files if options.files else defaultFile
options.GlobalTag = options.GlobalTag if options.GlobalTag else defaultGT
if not options.isData : #remove L2L3Residual correction from simulation options
jecs = [jc for jc in options.jetCorrections if jc!='L2L3Residual']
options.clearList('jetCorrections')
options.jetCorrections = jecs
return options
from PhysicsTools.PatAlgos.tools.helpers import cloneProcessingSnippet
from PhysicsTools.PatAlgos.tools.helpers import massSearchReplaceAnyInputTag
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
inputFiles = 'root://cmseos.fnal.gov//eos/uscms/store/user/zhangj/events/ALP/RunIISummer17DR94Premix/ALP_m10_w1_htjmin400_RunIISummer17DR94Premix_AODSIM_1.root'
options.maxEvents = 200
options.register('skipEvents', 0, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Events to skip")
options.register('reportEvery', 1, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Report every")
options.register('isMC', 1, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Sample is MC")
options.register('doSlimming', 1, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Output content is reduced")
options.register('isStandard', 0, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Output content is reduced")
options.register('isReMINIAOD', 1, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Output content is reduced")
options.register('numThreads', 8, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Set number of threads")
if options.isStandard:
options.doSlimming = 0
outputFile = 'root://cmseos.fnal.gov//store/user/mwulansa/DIS/TCPAnalysis/RunIISummer17DR94Premix/ALP_m10_w1_htjmin400_RunIISummer17DR94Premix_MINIAODSIM_Standard_1.root'
elif options.doSlimming:
# outputFile = 'root://cmseos.fnal.gov//store/user/mwulansa/DIS/TCPAnalysis/RunIISummer17DR94Premix/ALP_m10_w1_htjmin400_RunIISummer17DR94Premix_MINIAODSIM_Slimmed_1.root'
outputFile = 'test.root'
else:
VarParsing.multiplicity.singleton, VarParsing.varType.string,
"'POWHEG','aMC', 'MadGraph' or 'pythia8'")
options.register("analysisType", "SL", VarParsing.multiplicity.singleton,
VarParsing.varType.string, "'SL' or 'DL'")
options.register("globalTag", "74X_mcRun2_asymptotic_v2",
VarParsing.multiplicity.singleton, VarParsing.varType.string,
"global tag")
options.register(
"useJson", False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"apply the json filter (on the grid there are better ways to do this)")
options.parseArguments()
# re-set some defaults
if options.maxEvents is -1: # maxEvents is set in VarParsing class by default to -1
options.maxEvents = 100 # reset to 100 for testing
if not options.inputFiles:
options.inputFiles = [
'file:/pnfs/desy.de/cms/tier2/store/user/hmildner/ttHTobb_M125_13TeV_powheg_pythia8/Boostedv2MiniAOD/151017_154254/0000/BoostedTTH_MiniAOD_1.root'
]
# checks for correct values and consistency
if options.analysisType not in ["SL", "DL"]:
print "\n\nConfig ERROR: unknown analysisType '" + options.analysisType + "'"
print "Options are 'SL' or 'DL'\n\n"
sys.exit()
if "data" in options.globalTag.lower() and not options.isData:
print "\n\nConfig ERROR: GT contains seems to be for data but isData==False\n\n"
sys.exit()
if "mc" in options.globalTag.lower() and options.isData:
print "\n\nConfig ERROR: GT contains seems to be for MC but isData==True\n\n"
def options() :
from FWCore.ParameterSet.VarParsing import VarParsing as VP
options = VP('standard')
options.output = "SusyCAF_Tree.root"
options.secondaryOutput = "" #switch PAT-tuple output off by default
options.maxEvents = 100
options.register('GlobalTag', mytype = VP.varType.string)
options.register('quiet', default = True )
options.register('isData', default = True )
options.register('dropMore', default = False , info = "drop more for a smaller event size")
options.register('leptonPtMin', default = 0 )
options.register('patify', default = True, info = "run SUSYPAT on the fly")
options.register('AllTracks', default = False, info = "include all tracks")
options.register('SourceName', info = "'S:stream' or 'DS:dataset' to store those HLT paths", mytype = VP.varType.string,)
options.register('jetCollections', default = ['ak5calo','ak5pf'], # + ['ak7calo','ak7pf','ak5pf2pat','ak7pf2pat']
info = "jet types to store", mult = VP.multiplicity.list, mytype = VP.varType.string)
options.register('jetCorrections', default = ['L1FastJet','L2Relative','L3Absolute','L2L3Residual'], #L2L3Residual removed from options for simulation (below)
info = "jet correction levels to apply", mult = VP.multiplicity.list, mytype = VP.varType.string)
options.register('scan', default = "", info = "code for CMSSM or SMS scan", mytype = VP.varType.string)
options.register('triggers', default = -1, info = "store trigger information in the ntuple") # defaults to isData (below)
options.register('taus', default = True, info = "store tau information in the ntuple")
options.register('dqm', default = False, info = "store dqm information in the ntuple")
options.register('doTauReco', default = False, info = "redo tau recostruction (needed in order to read 51X files with 52X recipe)")
options.register('doTypeIMetReco', default = False, info = "do type I MET recostruction")
options.register('doTypeIMetPat', default = True, info = "do type I MET Pat (RA4 Recipe + also caloMET)")
options.register('doPfMetPhiCorrections', default = False, info = "do phi corrections in pfMET")
options.register('beamHaloVars', default = True, info = "store beam halo information in the ntuple")
options.register('doPhotonPFIso', default = True, info = "include simple cut based Photon IDs 2012 in the ntuple")
options.register('doPoolOutput', default = False, info = "Output edm file with full event content")
options.register('pdfSets', default = "", info = "pdf Sets to use for pdf Uncertainty weights", mytype = VP.varType.string)
options.register('defPdf', default = "cteq6ll.LHpdf", info = "pdf set used when producing sample", mytype = VP.varType.string)
__hack_ListVarparsingBug__( options, 'jetCollections')
__hack_ListVarparsingBug__( options, 'jetCorrections')
options.parseArguments()
options._tagOrder =[] # weird, but something to do with options.output
defaultGT,defaultFile = {
"53X" : [('START53_V10::All','/store/mc/Summer12_DR53X/ZZ_TuneZ2star_8TeV_pythia6_tauola/AODSIM/PU_S10_START53_V7A-v1/0000/FEEB0275-24FA-E111-B9A6-00266CF9B970.root'),
('GR_P_V41_AN1::All' ,'/store/data/Run2012C/SingleMu/AOD/PromptReco-v2/000/201/624/F40EFD20-9CF0-E111-8F16-5404A63886C4.root')],
"52X" : [('START52_V4::All','/store/relval/CMSSW_5_2_5_cand1/RelValTTbar/GEN-SIM-RECO/START52_V9-v1/0261/249F81B9-8C91-E111-AE72-003048679236.root'),
('GR_R_52_V4::All' ,'/store/relval/CMSSW_5_2_5_cand1/SingleMu/RECO/GR_R_52_V7_RelVal_mu2011B-v1/0262/7425B24B-9E91-E111-BA46-002618943800.root')],
"51X" : [('START50_V15::All','/store/relval/CMSSW_5_1_2/RelValTTbar/GEN-SIM-RECO/PU_START50_V15A-v1/0003/0AB02E79-7C61-E111-8594-002481E0D73C.root'),
('GR_R_50_V12::All' ,'/store/relval/CMSSW_5_1_2/Jet/RECO/GR_R_50_V12_RelVal_jet2010B-v1/0237/326E6DA4-B460-E111-AAFA-0026189438D6.root')],
"44X" : [('START44_V12::All','/../user/b/bbetchar//CMSSW_4_4_2_patch10/RelValProdTTbar/AODSIM/START44_V7_special_120119-v1/0088/0ECC7FA5-7A43-E111-9B36-002618943861.root'),
('GR_R_44_V13::All','/../user/b/bbetchar/data/Run2011A/SingleMu/AOD/08Nov2011-v1/0001/80D61727-200D-E111-B7F4-1CC1DE0571C8.root')],
"42X" : [('START42_V17::All','/../user/b/bbetchar/CMSSW_4_2_8/RelValProdTTbar/GEN-SIM-RECO/MC_42_V12-v1/0026/9C0E8835-9ABB-E011-95B0-0026189438BA.root'),
('GR_R_42_V24::All','/store/data/Run2011B/MultiJet/AOD/PromptReco-v1/000/175/832/485ABBC8-10DC-E011-980B-BCAEC518FF8E.root')]
}["53X"][options.isData]
options.files = options.files if options.files else defaultFile
options.GlobalTag = options.GlobalTag if options.GlobalTag else defaultGT
if options.triggers < 0 : options.triggers = int(options.isData)
if not options.isData : #remove L2L3Residual correction from simulation options
jecs = [jc for jc in options.jetCorrections if jc!='L2L3Residual']
options.clearList('jetCorrections')
options.jetCorrections = jecs
return options
def options():
from FWCore.ParameterSet.VarParsing import VarParsing as VP
options = VP('standard')
options.output = "DJ_Tree.root"
options.maxEvents = 100
options.register('GlobalTag', mytype=VP.varType.string)
options.register('quiet', default=True)
options.register('isData', default=True)
options.register('signal', default=False)
options.register(
'jetCorrections',
default=['L1FastJet', 'L2Relative', 'L3Absolute', 'L2L3Residual'
], #L2L3Residual removed from options for simulation (below)
info="jet correction levels to apply",
mult=VP.multiplicity.list,
mytype=VP.varType.string)
__hack_ListVarparsingBug__(options, 'jetCorrections')
options.parseArguments()
options._tagOrder = [] # weird, but something to do with options.output
defaultGT, defaultFile = {
"44X":
[('START44_V9B::All',
'/store/mc/Fall11/QCD_Pt-120to170_TuneZ2_7TeV_pythia6/AODSIM/PU_S6_START42_V14B-v1/0000/0A8E5F68-73F4-E011-AB1C-00E08178C0D5.root'
),
('START44_V9B::All',
'/store/mc/Fall11/HTo2LongLivedTo4F_MH-1000_MFF-350_CTau-350_7TeV-pythia6/GEN-SIM-RECODEBUG/DEBUG-PU_S6_START44_V9B-v4/0000/E83B1DB8-158B-E111-9DB4-0017A477003C.root'
),
('FT_R_44_V11::All',
'/store/data/Run2011B/HT/RECO/19Nov2011-v1/10001/B2E314BC-4746-E111-9BC1-00261894390A.root'
)],
"53X": [
#('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-1000_MFF-350_CTau35To3500_8TeV-pythia6/AODSIM/DEBUG_PU_S10_START53_V7A-v1/0000/9C61CA50-0BF3-E111-91E1-003048D2BB22.root'),
#('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-400_MFF-50_CTau8To800_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/00000/80A0705C-7518-E211-A17D-0030487E4D11.root'),
#('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-200_MFF-50_CTau20To2000_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/00000/7A2CB2BB-6C18-E211-AB1B-003048C69416.root'),
('START53_V7F::All',
'/store/mc/Summer12_DR53X/LongLivedChi0ToMuQQ_MSquark-1500_MChi-494_TuneZ2Star_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/00000/B69AB21F-D50F-E211-B2A7-00215E21DCA2.root'
),
('START53_V7F::All',
'/store/mc/Summer12_DR53X/QCD_Pt-170to300_TuneZ2star_8TeV_pythia6/AODSIM/PU_S10_START53_V7A-v2/00000/0CECB336-B70E-E211-B22D-0018F3D096CE.root'
),
('FT_53_V6_AN2::All',
'/store/data/Run2012B/JetHT/AOD/13Jul2012-v1/00000/04989369-BFD2-E111-B190-002618943933.root'
),
('FT_53_V10_AN2::All',
'/store/data/Run2012C/JetHT/AOD/24Aug2012-v2/00000/22AFA398-D600-E211-AD2E-00261894385A.root'
),
('GR_P_V41_AN2::All',
'/store/data/Run2012C/JetHT/AOD/PromptReco-v2/000/198/955/144C3A0C-D8CF-E111-8A6C-001D09F27067.root'
),
('GR_P_V42_AN2::All',
'/store/data/Run2012D/JetHT/AOD/PromptReco-v1/000/204/506/F8F3A290-8C10-E211-B2D3-001D09F2525D.root'
),
],
"53XCtrl": [
('START53_V7F::All',
'/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-1000_MFF-350_CTau35To3500_8TeV-pythia6/AODSIM/DEBUG_PU_S10_START53_V7A-v1/0000/9C61CA50-0BF3-E111-91E1-003048D2BB22.root'
),
('START53_V7F::All',
'/store/mc/Summer12_DR53X/QCD_Pt-170to300_TuneZ2star_8TeV_pythia6/AODSIM/PU_S10_START53_V7A-v2/00000/0CECB336-B70E-E211-B22D-0018F3D096CE.root'
),
('FT_53_V6_AN2::All',
'/store/data/Run2012B/SingleMu/AOD/13Jul2012-v1/0000/24DF4BD2-0DD4-E111-BB00-E0CB4E4408DF.root'
),
('FT_53_V10_AN2::All',
'/store/data/Run2012C/SingleMu/AOD/24Aug2012-v1/00000/10FD4F26-FAEE-E111-9AB5-00259073E4FC.root'
),
('GR_P_V41_AN2::All',
'/store/data/Run2012C/SingleMu/AOD/PromptReco-v2/000/198/955/247727CE-98CF-E111-B2EF-5404A63886A0.root'
),
('GR_P_V42_AN2::All',
'/store/data/Run2012D/SingleMu/AOD/PromptReco-v1/000/204/506/8078B6DD-9210-E211-80BE-BCAEC5329700.root'
),
]
}["53X"][int(not options.signal) + options.isData]
options.files = options.files if options.files else defaultFile
options.GlobalTag = options.GlobalTag if options.GlobalTag else defaultGT
if not options.isData: #remove L2L3Residual correction from simulation options
jecs = [jc for jc in options.jetCorrections if jc != 'L2L3Residual']
options.clearList('jetCorrections')
options.jetCorrections = jecs
return options
process.load('Configuration.StandardSequences.Services_cff')
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.StandardSequences.Generator_cff')
process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic8TeVCollision_cfi')
process.load("PhysicsTools.HepMCCandAlgos.genParticles_cfi")
process.load("RecoJets.Configuration.GenJetParticles_cff")
process.load("RecoJets.Configuration.RecoGenJets_cff")
# ================= Input Variable Parsing ======================
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('standard')
options.output = 'QCDAna.root'
options.maxEvents = 1000
options.register('tune',
"4C",
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Pythia tune")
options.register('processType',
"HardQCDSelection",
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Pythia process type with pT_hat range")
options.register('sqrtS',
'''
Job option for the B-initiated HNL generation - specialised for the Bc case
'''
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
# define the defaults here, changed from command line
options.maxEvents = -1 # -1 means all events, maxEvents considers the total over files considered
options.outputFile = 'BPH-test.root'
# add costum parameters
options.register(
'severityLevel',
'ERROR', # default value
VarParsing.multiplicity.singleton, # singleton or list
VarParsing.varType.string, # string, int, or float
'severity level for log messages, DEBUG, INFO, WARNING, ERROR')
options.register('nThr', 1, VarParsing.multiplicity.singleton,
VarParsing.varType.int, 'Number of threads')
options.register('seedOffset', 1, VarParsing.multiplicity.singleton,
VarParsing.varType.int, 'Seed offset')
options.register('mass', 1, VarParsing.multiplicity.singleton,
VarParsing.varType.float, 'mass of the HNL')
options.register('ctau', 100, VarParsing.multiplicity.singleton,
VarParsing.varType.float, 'ctau of the HNL [mm]')
options.register('doSkipMuonFilter', False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, 'Skip the muon filter')
options.register('doDisplFilter', False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
'In muon filter, add a cut on the HNL displacement')
options.register('doMajorana', False, VarParsing.multiplicity.singleton,
process.load('Configuration.StandardSequences.Services_cff')
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.EventContent.EventContent_cff')
process.load('SimGeneral.MixingModule.mixNoPU_cfi')
process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
process.load('Configuration.StandardSequences.MagneticField_38T_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff')
# using VarParsing module to setup analysis option
from FWCore.ParameterSet.VarParsing import VarParsing
myopts = VarParsing ('analysis')
##### set default value #####
myopts.maxEvents = 10
myopts.outputFile = 'miniaod.root'
myopts.register('basketSize',
16384,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
'Default ROOT basket size in output file.')
myopts.register('compressionLevel',
4,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
'ROOT compression level of output file.')
myopts.register('compressionAlgorithm',
'LZMA',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
## Configure message logger
process.load("FWCore.MessageLogger.MessageLogger_cfi")
process.MessageLogger.cerr.threshold = 'INFO'
process.MessageLogger.cerr.FwkReport.reportEvery = 100
## Define input
if not options.inputFiles:
if options.runOnGenOrAODsim:
options.inputFiles=['/store/mc/RunIISummer15GS/TTToSemiLeptonic_TuneCUETP8M1_alphaS01273_13TeV-powheg-scaledown-pythia8/GEN-SIM/MCRUN2_71_V1-v2/40000/DE7952A2-6E2F-E611-A803-001E673D1B21.root']
else:
options.inputFiles=['/store/mc/RunIIFall17MiniAOD/TTToSemiLeptonic_TuneCP5_PSweights_13TeV-powheg-pythia8/MINIAODSIM/94X_mc2017_realistic_v10-v1/50000/DC5D3109-F2E1-E711-A26E-A0369FC5FC9C.root']
## Define maximum number of events to loop over
if options.maxEvents is -1: # maxEvents is set in VarParsing class by default to -1
options.maxEvents = 1000 # reset for testing
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(int(options.maxEvents)))
process.source = cms.Source( "PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles),
skipEvents=cms.untracked.uint32(int(options.skipEvents)),
)
## Set input particle collections to be used by the tools
genParticleCollection = ''
genJetInputParticleCollection = ''
genJetCollection = 'ak4GenJetsCustom'
if options.runOnGenOrAODsim:
genParticleCollection = 'genParticles'
genJetInputParticleCollection = genParticleCollection
# import ROOT in batch mode
import sys
oldargv = sys.argv[:]
sys.argv = [ '-b-' ]
import ROOT
ROOT.gROOT.SetBatch(True)
sys.argv = oldargv
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('python')
#default options
options.inputFiles="/eos/cms/store/relval/CMSSW_9_4_0_pre3/RelValTTbar_13/MINIAODSIM/PU25ns_94X_mc2017_realistic_PixFailScenario_Run305081_FIXED_HS_AVE50-v1/10000/02B605A1-86C2-E711-A445-4C79BA28012B.root"
options.outputFile="jetmetNtuples.root"
options.maxEvents=-1
#overwrite if given any command line arguments
options.parseArguments()
# define deltaR
from math import hypot, pi, sqrt, fabs
import numpy as n
from jetmet_tree import *
from functions import *
# create an oput tree.
#fout = ROOT.TFile ("jetmet.root","recreate")
fout = ROOT.TFile (options.outputFile,"recreate")
t = ROOT.TTree ("events","events")
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.EventContent.EventContent_cff')
process.load('SimGeneral.MixingModule.mixNoPU_cfi')
process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
process.load('Configuration.StandardSequences.MagneticField_38T_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load(
'Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff'
)
# using VarParsing module to setup analysis option
from FWCore.ParameterSet.VarParsing import VarParsing
myopts = VarParsing('analysis')
##### set default value #####
myopts.maxEvents = 10
myopts.outputFile = 'miniaod.root'
myopts.register('basketSize', 16384, VarParsing.multiplicity.singleton,
VarParsing.varType.int,
'Default ROOT basket size in output file.')
myopts.register('compressionLevel', 4, VarParsing.multiplicity.singleton,
VarParsing.varType.int,
'ROOT compression level of output file.')
myopts.register(
'compressionAlgorithm', 'LZMA', VarParsing.multiplicity.singleton,
VarParsing.varType.string,
'Algorithm used to compress data in the ROOT output file, allowed values are ZLIB and LZMA.'
)
myopts.register(
'fastCloning', False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
## Configure message logger
process.load("FWCore.MessageLogger.MessageLogger_cfi")
process.MessageLogger.cerr.threshold = 'INFO'
process.MessageLogger.cerr.FwkReport.reportEvery = 100
## Define input
if not options.inputFiles:
if options.runOnGenOrAODsim:
options.inputFiles=['/store/mc/RunIISummer15GS/TTToSemiLeptonic_TuneCUETP8M1_alphaS01273_13TeV-powheg-scaledown-pythia8/GEN-SIM/MCRUN2_71_V1-v2/40000/DE7952A2-6E2F-E611-A803-001E673D1B21.root']
else:
options.inputFiles=['/store/mc/RunIISpring15MiniAODv2/ttbb_4FS_ckm_amcatnlo_madspin_pythia8/MINIAODSIM/74X_mcRun2_asymptotic_v2-v1/40000/06D46D97-C66D-E511-9ABF-00266CFAE20C.root']
## Define maximum number of events to loop over
if options.maxEvents is -1: # maxEvents is set in VarParsing class by default to -1
options.maxEvents = 1000 # reset for testing
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(int(options.maxEvents)))
process.source = cms.Source( "PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles),
skipEvents=cms.untracked.uint32(int(options.skipEvents)),
)
## Set input particle collections to be used by the tools
genParticleCollection = ''
genJetInputParticleCollection = ''
genJetCollection = 'ak4GenJetsCustom'
if options.runOnGenOrAODsim:
genParticleCollection = 'genParticles'
genJetInputParticleCollection = genParticleCollection
process.load('Configuration.StandardSequences.MagneticField_38T_PostLS1_cff')
process.load('Configuration.StandardSequences.Reconstruction_cff')
#process.load('Configuration.StandardSequences.MagneticField_38T_PostLS1_cff')
#process.load('RecoTracker.CkfPattern.CkfTrackCandidates_cff')
###
# Defaults
###
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
#options.inputFiles = "files.list"
#options.inputFiles = "file:GEM2019Upg14.root"
options.outputFile = "file:rereco_fail_me21_me31_me41_0p0000.root"
options.maxEvents = "-1"
#options.histFile = "validation_hists_0p0.root"
#options.register ('failureRate', 0.15, VarParsing.multiplicity.singleton, VarParsing.varType.float, "choose failure probability")
#options.register ('txt_input', "files.list", VarParsing.multiplicity.singleton, VarParsing.varType.string, "List of input files")
#options.register ('histFile', "validation_hists_me21_me32_me41_0p0.root", VarParsing.multiplicity.singleton, VarParsing.varType.string, "Output validation histograms file")
options.register ('failRateME21', 0.0, VarParsing.multiplicity.singleton, VarParsing.varType.float, "ME2/1 failure probability")
options.register ('failRateME31', 0.0, VarParsing.multiplicity.singleton, VarParsing.varType.float, "ME3/1 failure probability")
options.register ('failRateME41', 0.0, VarParsing.multiplicity.singleton, VarParsing.varType.float, "ME4/1 failure probability")
options.register ('doUniformFailure', True, VarParsing.multiplicity.singleton, VarParsing.varType.bool, "run failure in all chambers evenly")
options.register ('doCFEBFailure', True, VarParsing.multiplicity.singleton, VarParsing.varType.bool, "explicitly set CFEB failure")
#options.register ('maxEvents', -1, VarParsing.multiplicity.singleton, VarParsing.varType.int, "max number of events to process") -- this is already registered by default
options.parseArguments()
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(options.maxEvents)
)
def main():
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.inputFiles = [
"/afs/cern.ch/work/l/lata/HH_bbyy/CMSSW_9_4_2/src/flashgg/root_files/myMicroAODOutput_6.root"
]
options.maxEvents = -1
options.outputFile = "fout.root"
options.parseArguments()
###### save hisots in root file ######
fout = ROOT.TFile("output_3.root", 'RECREATE')
fout.cd()
##### define Histograms ######
h_H_pt = ROOT.TH1F("gen_H_pT", ";Higgs p_{T} [GeV];Events;;", 50, 0., 500.)
h_H_eta = ROOT.TH1F("gen_H_eta", ";Higgs #eta;Events;;", 50, -5., 5.)
h_b_pt = ROOT.TH1F("gen_b_pT", ";p_{T}of b [GeV];Events;;", 50, 0., 500.)
h_b_eta = ROOT.TH1F("gen_b_eta", ";#eta of b;Events;;", 50, -5., 5.)
h_y_pt = ROOT.TH1F("gen_y_pT", ";p_{T}of #gamma[GeV];Events;;", 50, 0.,
500)
h_y_eta = ROOT.TH1F("gen_y_eta", ";#eta of #gamma;Events;;", 50, -5., 5.)
h_jet_pt = ROOT.TH1F("jet_pt", ";jet p_{T} [GeV];Events;;", 50, 0., 500.)
h_jet_eta = ROOT.TH1F("jet_eta", ";jet #eta;Events;;", 50, -5., 5.)
h_jet_deepcsv = ROOT.TH1F("jet_deepcsv", ";jet DeepCSV [GeV];Events;;", 50,
0., 1.)
h_btaggedjet0_pt = ROOT.TH1F("btagged_jet0_pt",
";b-tagged leading jet p_{T} [GeV];Events;;",
50, 0., 500.)
h_btaggedjet0_eta = ROOT.TH1F("btagged_jet0_eta",
";b-tagged leading jet #eta;Events;;", 50,
-5., 5.)
h_btaggedjet1_pt = ROOT.TH1F(
"btagged_jet1_pt", ";b-tagged subleading jet p_{T} [GeV];Events;;", 50,
0, 500)
h_btaggedjet1_eta = ROOT.TH1F("btagged_jet1_eta",
";b-tagged subleading jet #eta;Events;;", 50,
-5, 5)
h_photon_pt = ROOT.TH1F("photon_pt", ";#gamma p_{T} [GeV];Events;;", 50,
0., 500.)
h_photon_eta = ROOT.TH1F("photon_eta", ";#gamma #eta [GeV];Events;;", 50,
-5, 5)
h_photon0_pt = ROOT.TH1F("photon0_pt",
";leading #gamma p_{T} [GeV];Events;;", 50, 0.,
500.)
h_photon0_eta = ROOT.TH1F("photon0_eta",
";leading #gamma #eta [GeV];Events;;", 50, -5, 5)
h_photon1_pt = ROOT.TH1F("photon1_pt",
";subleading #gamma p_{T} [GeV];Events;;", 50, 0.,
500)
h_photon1_eta = ROOT.TH1F("photon1_eta",
";subleading #gamma #eta [GeV];Events;;", 50, -5,
5)
h_njet = ROOT.TH1D("njet", ";Njet;Events;;", 15, 0, 15)
h_nphoton = ROOT.TH1D("nphoton", ";Nphoton;Events;;", 5, 0, 5)
h_dipho_mass = ROOT.TH1F("dipho_mass ", ";diphoton mass [GeV];Events;;",
50, 0., 200.)
h_reconsM_photon = ROOT.TH1F("photon_recomass",
";M_{#gamma#gamma} [GeV];Events;;", 50, 0.,
200.)
h_reconsM_bb = ROOT.TH1F("bb_recomass", ";M_{bb} [GeV];Events;;", 50, 0.,
200)
events = Events(options.inputFiles)
nevents = 0
for event in events:
if options.maxEvents > 0 and nevents > options.maxEvents: break
### Find gen particles
h_prunedGenpar = Handle("std::vector<reco::GenParticle>")
event.getByLabel("flashggPrunedGenParticles", h_prunedGenpar)
particles = h_prunedGenpar.product()
for p in particles:
print('prunedgenparticles P pid = {} status = {}'.format(
p.pdgId(), p.status()))
if p.pdgId() == 25:
h_H_pt.Fill(p.pt())
h_H_eta.Fill(p.eta())
if abs(p.pdgId()) == 5:
h_b_pt.Fill(p.pt())
h_b_eta.Fill(p.eta())
if p.pdgId() == 22:
h_y_pt.Fill(p.pt())
h_y_eta.Fill(p.eta())
h_photons = Handle("std::vector<flashgg::Photon>")
event.getByLabel("flashggRandomizedPhotons", h_photons)
print("N(photons) = {}".format(len(h_photons.product())))
photons = h_photons.product()
h_nphoton.Fill(len(photons))
if (len(photons) >= 2):
p4_g0 = ROOT.TLorentzVector()
p4_g1 = ROOT.TLorentzVector()
p4_g0.SetPtEtaPhiM(photons[0].pt(), photons[0].eta(),
photons[0].phi(), photons[0].mass())
p4_g1.SetPtEtaPhiM(photons[1].pt(), photons[1].eta(),
photons[1].phi(), photons[1].mass())
h_photon0_pt.Fill(photons[0].pt())
h_photon1_pt.Fill(photons[1].pt())
h_photon0_eta.Fill(photons[0].eta())
h_photon1_eta.Fill(photons[1].eta())
if photons[0].pt() > 20 and photons[1].pt() > 20:
h_reconsM_photon.Fill((p4_g0 + p4_g1).M())
for gam in photons:
print('photon pt = {}'.format(gam.pt()))
h_photon_eta.Fill(gam.eta())
h_photon_pt.Fill(gam.pt())
h_packedGenpar = Handle("std::vector<pat::PackedGenParticle>")
event.getByLabel("flashggGenPhotons", h_packedGenpar)
particles = h_packedGenpar.product()
for p in particles:
print('packedgenparticles P pid = {} status = {}'.format(
p.pdgId(), p.status()))
h_jets = Handle("std::vector<std::vector<flashgg::Jet> >")
event.getByLabel("flashggFinalJets", h_jets)
print("N(jets) = %i" % len(h_jets.product()))
if len(h_jets.product()) <= 0: continue
jets = h_jets.product()[0]
h_njet.Fill(len(jets))
if (len(jets) >= 2):
disc1 = jets[0].bDiscriminator("pfDeepCSVJetTags:probb") + jets[
0].bDiscriminator("pfDeepCSVJetTags:probbb")
disc2 = jets[1].bDiscriminator("pfDeepCSVJetTags:probb") + jets[
1].bDiscriminator("pfDeepCSVJetTags:probbb")
pt0 = jets[0].pt()
pt1 = jets[1].pt()
eta0 = jets[0].eta()
eta1 = jets[1].eta()
p4_b0 = ROOT.TLorentzVector()
p4_b1 = ROOT.TLorentzVector()
p4_b0.SetPtEtaPhiM(jets[0].pt(), jets[0].eta(), jets[0].phi(),
jets[0].mass())
p4_b1.SetPtEtaPhiM(jets[1].pt(), jets[1].eta(), jets[1].phi(),
jets[1].mass())
if disc1 > 0.4941 and disc2 > 0.4941:
h_btaggedjet0_pt.Fill(pt0)
h_btaggedjet1_pt.Fill(pt1)
h_btaggedjet0_eta.Fill(eta0)
h_btaggedjet1_eta.Fill(eta1)
if pt0 > 25 and pt1 > 25 and abs(eta0) <= 2.5 and abs(
eta1) <= 2.5:
h_reconsM_bb.Fill((p4_b0 + p4_b1).M())
for jet in jets:
print("jet pt = %f DeepCSVBDisc = %f" %
(jet.pt(), jet.bDiscriminator("pfDeepCSVJetTags:probb") +
jet.bDiscriminator("pfDeepCSVJetTags:probbb")))
deepCSV = jet.bDiscriminator(
"pfDeepCSVJetTags:probb") + jet.bDiscriminator(
"pfDeepCSVJetTags:probbb")
h_jet_pt.Fill(jet.pt())
h_jet_eta.Fill(jet.eta())
h_jet_deepcsv.Fill(deepCSV)
h_dipho = Handle("std::vector<flashgg::DiPhotonCandidate>")
event.getByLabel("flashggDiPhotons", h_dipho)
diphos = h_dipho.product()
print("N(diphotons) = %i" % len(diphos))
for dipho in diphos:
print("Diphoton pt = %f" % dipho.pt())
h_dipho_mass.Fill(dipho.mass())
nevents += 1
fout.Write()
fout.Close()
Authors: L. Dodd, N. Woods, I. Ojalvo, S. Dasu, M. Cepeda, E. Friis (UW Madison)
'''
import FWCore.ParameterSet.Config as cms
import os
# Get command line options
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
# Set useful defaults
options.inputFiles = ('/store/data/Run2012C/ZeroBias1/RAW/v1/000/198/609/04F0C0E3-72CA-E111-A802-003048F117EC.root',
'/store/data/Run2012C/ZeroBias1/RAW/v1/000/198/609/08361F06-71CA-E111-B50E-003048F1C424.root')
options.outputFile = "test_pu_tree.root"
options.maxEvents = 555
options.register(
"puMax",
7,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Threshold for considering a region to be pileup only")
options.parseArguments()
process = cms.Process("puMethodComparison")
process.load('Configuration/StandardSequences/FrontierConditions_GlobalTag_cff')
process.GlobalTag.globaltag = 'GR_H_V28::All'
inputFiles=myFile.root outputFile=outputFile.root maxEvents=-1
Authors: L. Dodd, N. Woods, I. Ojalvo, S. Dasu, M. Cepeda, E. Friis (UW Madison)
'''
import FWCore.ParameterSet.Config as cms
import os
# Get command line options
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
# Set useful defaults
options.inputFiles = $inputFileNames
options.outputFile = "$outputFileName"
options.maxEvents = -1
options.register(
'eicIsolationThreshold',
3,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"EIC Isolation threshold")
options.register(
'hActivityCut',
0.5,
VarParsing.multiplicity.singleton,
VarParsing.varType.float,
"HCAL activity threshold")
options.register(
"stage1B",
0,
options = VarParsing ('analysis')
options.register ('globalTag',
'MCRUN2_74_V9A',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
'GlobalTag')
options.register ('TriggerSet',
'SingleMu', # SingleMu or DoubleMu
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
'TirggerSet')
#-------------------------------------------------------------------------------
# defaults
options.maxEvents = -1 #all events ---> already the default
#options.maxEvents = 500 #all events ---> already the default
options.outputFile = 'tnpZ_MC.root'
print " options.maxEvents = ",options.maxEvents
print " options.outputFile = ",options.outputFile
#-------------------------------------------------------------------------------
options.parseArguments()
# load configurations
#label = options.label
#globalTag = options.globalTag + "::All"
globalTag = options.globalTag
TriggerSet = options.TriggerSet
print " globalTag = ",globalTag
print " TriggerSet = ", TriggerSet
# maxEvents: singleton, int; default = -1
# inputFiles: (comma separated, no spaces!) list, string: default empty
options.register( "outName", "testrun", VarParsing.multiplicity.singleton, VarParsing.varType.string, "name and path of the output files (without extension)" )
options.register( "weight", 0.01, VarParsing.multiplicity.singleton, VarParsing.varType.float, "xs*lumi/(nPosEvents-nNegEvents)" )
options.register( "isData", False, VarParsing.multiplicity.singleton, VarParsing.varType.bool, "is it data or MC?" )
options.register( "isBoostedMiniAOD", True, VarParsing.multiplicity.singleton, VarParsing.varType.bool, "has the file been prepared with the BoostedProducer ('custom' MiniAOD)?" )
options.register( "makeSystematicsTrees", True, VarParsing.multiplicity.singleton, VarParsing.varType.bool, "do you need all systematics (e.g. to calculate limits)?" )
options.register( "generatorName", "notSpecified", VarParsing.multiplicity.singleton, VarParsing.varType.string, "'POWHEG','aMC', 'MadGraph' or 'pythia8'" )
options.register( "analysisType", "SL", VarParsing.multiplicity.singleton, VarParsing.varType.string, "'SL' or 'DL'" )
options.register( "globalTag", "74X_mcRun2_asymptotic_v2", VarParsing.multiplicity.singleton, VarParsing.varType.string, "global tag" )
options.register( "useJson",False, VarParsing.multiplicity.singleton, VarParsing.varType.bool, "apply the json filter (on the grid there are better ways to do this)" )
options.parseArguments()
# re-set some defaults
if options.maxEvents is -1: # maxEvents is set in VarParsing class by default to -1
options.maxEvents = 100 # reset to 100 for testing
if not options.inputFiles:
options.inputFiles=['file:/pnfs/desy.de/cms/tier2/store/user/hmildner/ttHTobb_M125_13TeV_powheg_pythia8/Boostedv2MiniAOD/151017_154254/0000/BoostedTTH_MiniAOD_1.root']
# checks for correct values and consistency
if options.analysisType not in ["SL","DL"]:
print "\n\nConfig ERROR: unknown analysisType '"+options.analysisType+"'"
print "Options are 'SL' or 'DL'\n\n"
sys.exit()
if "data" in options.globalTag.lower() and not options.isData:
print "\n\nConfig ERROR: GT contains seems to be for data but isData==False\n\n"
sys.exit()
if "mc" in options.globalTag.lower() and options.isData:
print "\n\nConfig ERROR: GT contains seems to be for MC but isData==True\n\n"
sys.exit()
if not options.inputFiles:
process.load('Configuration.StandardSequences.Services_cff')
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.StandardSequences.Generator_cff')
process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic8TeVCollision_cfi')
process.load("PhysicsTools.HepMCCandAlgos.genParticles_cfi")
process.load("RecoJets.Configuration.GenJetParticles_cff")
process.load("RecoJets.Configuration.RecoGenJets_cff")
# ================= Input Variable Parsing ======================
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('standard')
options.output = 'QCDAna.root'
options.maxEvents = 200000
options.register('tune',
"Z2",
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Pythia tune")
options.register('processType',
"NSD",
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Pythia process type with pT_hat range")
options.register('sqrtS',
2760.0,
import FWCore.ParameterSet.Config as cms
from PhysicsTools.PatAlgos.tools.helpers import getPatAlgosToolsTask
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
#options.inputFiles = '/store/mc/RunIIFall17MiniAODv2/TTToHadronic_TuneCP5_13TeV-powheg-pythia8/MINIAODSIM/PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/90000/DCFE3F5F-AE42-E811-B6DB-008CFAF72A64.root'
options.inputFiles = 'file:/storage/local/data1/home/jduarte1/forPatrick/FFA0194D-1BBC-EF4F-9B8F-8FBED2C62FC8.root'
#options.inputFiles = '/store/mc/RunIISummer19UL17MiniAOD/TTToHadronic_TuneCP5_13TeV-powheg-pythia8/MINIAODSIM/106X_mc2017_realistic_v6-v4/30000/FFA0194D-1BBC-EF4F-9B8F-8FBED2C62FC8.root'
#options.inputFiles = 'file:FFA0194D-1BBC-EF4F-9B8F-8FBED2C62FC8.root'
options.maxEvents = 100
options.parseArguments()
from Configuration.ProcessModifiers.enableSonicTriton_cff import enableSonicTriton
process = cms.Process('PATtest', enableSonicTriton)
## MessageLogger
process.load("FWCore.MessageLogger.MessageLogger_cfi")
process.MessageLogger.cerr.FwkReport.reportEvery = 1
## Options and Output Report
process.options = cms.untracked.PSet(wantSummary=cms.untracked.bool(True))
## Source
process.source = cms.Source("PoolSource",
fileNames=cms.untracked.vstring(
options.inputFiles))
## Maximal Number of Events
process.maxEvents = cms.untracked.PSet(
input=cms.untracked.int32(options.maxEvents))
process.load("HeterogeneousCore.SonicTriton.TritonService_cff")
## 9
# '/store/user/icali/HIMinBiasUPC/HIMinBiasUPC_Skim_HLT_HIMinBiasHfOrBSC_Centrality0-10//64ca16868e481177958780733023cfa2/SD_MB_Cen0_10_123_1_soE.root',
# '/store/user/icali/HIMinBiasUPC/HIMinBiasUPC_Skim_HLT_HIMinBiasHfOrBSC_Centrality0-10//64ca16868e481177958780733023cfa2/SD_MB_Cen0_10_124_1_p5Y.root',
## 10
#'/store/user/icali/HIMinBiasUPC/HIMinBiasUPC_Skim_HLT_HIMinBiasHfOrBSC_Centrality0-10//64ca16868e481177958780733023cfa2/SD_MB_Cen0_10_126_1_pCl.root',
#'/store/user/icali/HIMinBiasUPC/HIMinBiasUPC_Skim_HLT_HIMinBiasHfOrBSC_Centrality0-10//64ca16868e481177958780733023cfa2/SD_MB_Cen0_10_128_1_FuJ.root',
#"file:/afs/cern.ch/work/r/richard/public/HI_L1_FirmwareTesting/Mismatched_HI_Jets_RAW.root"
#"file:/afs/cern.ch/work/r/richard/public/HI_L1_FirmwareTesting/RUN_258694/l1tCalo_EDM_258694_0.root"
#"root://xrootd.cmsaf.mit.edu//store/user/mnguyen/Hydjet_Quenched_MinBias_5020GeV/HydjetMB_740pre8_MCHI2_74_V3_53XBS_DIGI-RAW/6da45e4e90741bc03dbd9aec5f36c050/step2_DIGI_L1_DIGI2RAW_HLT_RAW2DIGI_L1Reco_100_1_nRy.root"
#"file:/afs/cern.ch/work/r/richard/public/HI_L1_FirmwareTesting/RUN_259637/RUN_259637_centralityMismatches.root"
"file:/afs/cern.ch/work/r/richard/public/HI_L1_FirmwareTesting/RUN_259637/RUN_259637_egammaMismatches.root"
)
options.maxEvents = 119
options.parseArguments()
process = cms.Process("ORSCPATTERNS")
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff')
from Configuration.AlCa.GlobalTag_condDBv2 import GlobalTag
process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_data', '')
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(options.maxEvents)
)
process.source = cms.Source(
"PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles)
