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)