def main():
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.inputFiles = "/afs/cern.ch/work/d/devdatta/CMSREL/Flashgg/CMSSW_9_4_2/src/flashgg/myMicroAODOutputFile.root"
options.parseArguments()
events = Events(options.inputFiles)
nevents = 0
for event in events:
h_jets = Handle("std::vector<std::vector<flashgg::Jet> >")
event.getByLabel("flashggFinalJets", h_jets)
print("N(jets) = %i" % len(h_jets.product()))
jets = h_jets.product()[0]
for jet in jets:
print("jet pt = %f DeepCSVBDisc = %f" %
(jet.pt(), jet.bDiscriminator("pfDeepCSVJetTags:probb") +
jet.bDiscriminator("pfDeepCSVJetTags:probbb")))
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())
def main():
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
options.inputFiles = "/afs/cern.ch/work/d/devdatta/CMSREL/Flashgg/CMSSW_9_4_2/src/flashgg/myMicroAODOutputFile.root"
options.parseArguments()
events = Events(options.inputFiles)
nevents = 0
for event in events:
h_jets = Handle("std::vector<std::vector<flashgg::Jet> >")
event.getByLabel("flashggFinalJets", h_jets)
print( "N(jets) = %i" % len(h_jets.product()) )
jets = h_jets.product()[0]
for jet in jets:
print( "jet pt = %f DeepCSVBDisc = %f" % (jet.pt(), jet.bDiscriminator("pfDeepCSVJetTags:probb")+jet.bDiscriminator("pfDeepCSVJetTags:probbb")) )
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() )
relBase = os.environ['CMSSW_BASE']
## PARSE ARGUMENTS
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.register('isMC', '', VarParsing.multiplicity.singleton,
VarParsing.varType.bool, 'Is MC')
options.register('isTTbar', '', VarParsing.multiplicity.singleton,
VarParsing.varType.bool, 'Is TTbar')
options.isMC = True
options.isTTbar = False
options.inputFiles = [
#matched with ~jmanagan/nobackup/LJMet94X_1lep_013019_logs/nominal/TTTT_TuneCP5_PSweights_13TeV-amcatnlo-pythia8/producer_TTTT_TuneCP5_PSweights_13TeV-amcatnlo-pythia8_1.py
'root://cmsxrootd.fnal.gov//store/mc/RunIIFall17MiniAODv2/TTTT_TuneCP5_PSweights_13TeV-amcatnlo-pythia8/MINIAODSIM/PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/70000/0ED34A55-DD52-E811-91CC-E0071B73B6B0.root',
'root://cmsxrootd.fnal.gov//store/mc/RunIIFall17MiniAODv2/TTTT_TuneCP5_PSweights_13TeV-amcatnlo-pythia8/MINIAODSIM/PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/70000/E2A4455F-FA53-E811-8017-E0071B7A8560.root',
#'root://cmsxrootd.fnal.gov//store/data/Run2017B/SingleMuon/MINIAOD/31Mar2018-v1/80000/4462B46E-653C-E811-BA97-0025905A6064.root'
#'root://cmsxrootd.fnal.gov//store/mc/RunIIFall17MiniAODv2/TprimeTprime_M-1800_TuneCP5_13TeV-madgraph-pythia8/MINIAODSIM/PU2017_12Apr2018_94X_mc2017_realistic_v14-v2/80000/12A585B9-F46B-E811-A775-FA163EFD0C51.root'
#'root://cmsxrootd.fnal.gov//store/mc/RunIIFall17MiniAODv2/TTToSemiLeptonic_TuneCP5_PSweights_13TeV-powheg-pythia8/MINIAODSIM/PU2017_12Apr2018_94X_mc2017_realistic_v14-v1/50000/5E7E4AA9-0743-E811-999A-0CC47A7C35A8.root'
]
options.maxEvents = 100
options.parseArguments()
isMC = options.isMC
isTTbar = options.isTTbar
## LJMET
process = cms.Process("LJMET")
## MessageLogger
process.load("FWCore.MessageService.MessageLogger_cfi")
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")
options.register ('user',
'',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Name the user. If not set by crab, this script will determine it.")
# defaults
#options.inputFiles = 'root:///user/jschulz/CMSSW_8_0_20/src/TreeWriter/ZNuNuGJets_MonoPhoton_PtG-130_TuneCUETP8M1_13TeV-madgraph_PUSpring16RAWAODSIM_80X_mcRun2_asymptotic_2016_miniAODv2_v0-v1_MINIAODSIM.root'
#options.inputFiles = 'root://cms-xrd-global.cern.ch//store/mc/RunIISpring16MiniAODv2/GJets_HT-600ToInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_miniAODv2_v0-v1/00000/264A540A-571A-E611-8C5E-0025904E3FCE.root'
#options.inputFiles = 'root://cms-xrd-global.cern.ch//store/mc/RunIISummer16MiniAODv2/TTJets_HT-600to800_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/MINIAODSIM/PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6_ext1-v1/60000/00E7F059-0BD5-E611-9267-001E67397CB5.root'
#options.inputFiles = 'root://cms-xrd-global.cern.ch//store/mc/RunIISummer16MiniAODv2/WJetsToLNu_HT-800To1200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/MINIAODSIM/PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6_ext1-v1/120000/02DD5E46-7ABE-E611-8F20-0025905B8582.root'
#options.inputFiles = 'root://cms-xrd-global.cern.ch//store/mc/RunIISummer16MiniAODv2/WGJets_MonoPhoton_PtG-130_TuneCUETP8M1_13TeV-madgraph/MINIAODSIM/PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/50000/BC527183-C0B7-E611-BC15-001E67348055.root'
#options.inputFiles = 'root://cms-xrd-global.cern.ch//store/mc/RunIISpring16MiniAODv2/SMS-T5Wg_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/MINIAODSIM/PUSpring16Fast_80X_mcRun2_asymptotic_2016_miniAODv2_v0-v1/80000/F227DD10-813E-E611-A722-6C3BE5B5C460.root'
#options.inputFiles = 'root://cms-xrd-global.cern.ch//store/data/Run2016B/SinglePhoton/MINIAOD/03Feb2017_ver1-v1/100000/04CFB75E-12EE-E611-918B-02163E0125C4.root'
#options.inputFiles = 'root:///user/kiesel/root-files/johannes/SinglePhoton_Run2016partBtoH-03Feb2017_MINIAOD.root'
options.inputFiles = 'root:///user/kiesel/root-files/johannes/SinglePhoton_Run2016C-03Feb2017-v1_MINIAOD.root'
options.inputFiles = 'root:///user/kiesel/root-files/johannes/SinglePhoton_Run2016D-03Feb2017-v1_MINIAOD.root'
options.inputFiles = 'root:///user/kiesel/root-files/johannes/SinglePhoton_Run2016E-03Feb2017-v1_MINIAOD.root'
options.inputFiles = 'root:///user/kiesel/root-files/johannes/SinglePhoton_Run2016F-03Feb2017-v1_MINIAOD.root'
options.inputFiles = 'root:///user/kiesel/root-files/johannes/SinglePhoton_Run2016G-03Feb2017-v1_MINIAOD.root'
#options.inputFiles = 'root:///user/kiesel/root-files/johannes/SinglePhoton_Run2016H-03Feb2017_ver2-v1_MINIAOD.root'
#options.inputFiles = 'root:///user/kiesel/root-files/johannes/SinglePhoton_Run2016H-03Feb2017_ver3-v1_MINIAOD.root'
options.outputFile = 'photonTree.root'
options.maxEvents = -1
# get and parse the command line arguments
options.parseArguments()
dataset=options.dataset or guessDatasetFromFileName(options.inputFiles[0])
print "Assumed dataset:", dataset
isRealData=not dataset.endswith("SIM")
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
# input cmsRun options
options = VarParsing("analysis")
options.inputFiles = (
"/store/data/Run2015B/SingleMuon/MINIAOD/PromptReco-v1/000/251/244/00000/68275270-7C27-E511-B1F0-02163E011A46.root"
)
#'/store/mc/RunIISpring15DR74/ZZTo4L_13TeV_powheg_pythia8/MINIAODSIM/Asympt50ns_MCRUN2_74_V9A-v1/00000/2C0279E7-9B16-E511-B5FF-90B11C2ACF20.root'
#'/store/data/Run2015B/SingleMuon/MINIAOD/PromptReco-v1/000/251/162/00000/160C08A3-4227-E511-B829-02163E01259F.root'
#'/store/data/Run2015B/SingleMuon/MINIAOD/PromptReco-v1/000/251/162/00000/160C08A3-4227-E511-B829-02163E01259F.root'
#'/store/mc/RunIISpring15DR74/ZZTo4L_13TeV_powheg_pythia8/MINIAODSIM/Asympt50ns_MCRUN2_74_V9A-v1/00000/2C0279E7-9B16-E511-B5FF-90B11C2ACF20.root'
#'/store/mc/RunIISpring15DR74/ZprimeToTauTau_M_500_TuneCUETP8M1_tauola_13TeV_pythia8/MINIAODSIM/Asympt25ns_MCRUN2_74_V9-v1/60000/CAE0DF38-7D19-E511-9DA4-0025904A8EC8.root'
# options.inputFiles = '/store/mc/RunIISpring15DR74/SUSYGluGluToHToTauTau_M-250_TuneCUETP8M1_13TeV-pythia8/MINIAODSIM/Asympt25ns_MCRUN2_74_V9-v1/80000/6EBE4F55-4E03-E511-9FE6-0CC47A13D416.root'
options.outputFile = "dataMCcomp.root"
# whether running MC or data
options.register(
"isMC",
0,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Set to 1 for simulated samples - updates GT, emulates HCAL TPGs.",
)
options.setDefault('maxEvents', options.maxE)
options.setDefault('tag', options.tg)
extension = {False: 'data', True: 'mc'}
outputFileNANO = cms.untracked.string(
'_'.join(['BParkNANO', extension[options.isMC], options.tag]) + '.root')
outputFileFEVT = cms.untracked.string(
'_'.join(['BParkFullEvt', extension[options.isMC], options.tag]) + '.root')
if not options.inputFiles:
if options.isMC:
options.inputFiles = [
'/store/mc/RunIIAutumn18MiniAOD/BuToKMuMu_probefilter_SoftQCDnonD_TuneCP5_13TeV-pythia8-evtgen/MINIAODSIM/PUPoissonAve20_BParking_102X_upgrade2018_realistic_v15-v2/40000/E5179C93-BFAA-6240-AAE7-A0A4DC410E4C.root',
'/store/mc/RunIIAutumn18MiniAOD/BuToKMuMu_probefilter_SoftQCDnonD_TuneCP5_13TeV-pythia8-evtgen/MINIAODSIM/PUPoissonAve20_BParking_102X_upgrade2018_realistic_v15-v2/40000/7F92D828-F1D0-2244-8B21-5224B26E974B.root',
'/store/mc/RunIIAutumn18MiniAOD/BuToKMuMu_probefilter_SoftQCDnonD_TuneCP5_13TeV-pythia8-evtgen/MINIAODSIM/PUPoissonAve20_BParking_102X_upgrade2018_realistic_v15-v2/40000/7B6B1A4D-B9E0-7245-988D-FB7BF863E454.root',
'/store/mc/RunIIAutumn18MiniAOD/BuToKMuMu_probefilter_SoftQCDnonD_TuneCP5_13TeV-pythia8-evtgen/MINIAODSIM/PUPoissonAve20_BParking_102X_upgrade2018_realistic_v15-v2/40000/79B0D428-E66B-1F4C-B05D-70EA12884B46.root',
'/store/mc/RunIIAutumn18MiniAOD/BuToKMuMu_probefilter_SoftQCDnonD_TuneCP5_13TeV-pythia8-evtgen/MINIAODSIM/PUPoissonAve20_BParking_102X_upgrade2018_realistic_v15-v2/40000/6EAEE0F7-3D9B-3341-AB81-54CECA85CEC0.root',
'/store/mc/RunIIAutumn18MiniAOD/BuToKMuMu_probefilter_SoftQCDnonD_TuneCP5_13TeV-pythia8-evtgen/MINIAODSIM/PUPoissonAve20_BParking_102X_upgrade2018_realistic_v15-v2/40000/46BE0C7D-4B50-6F42-8451-418494B481D5.root',
'/store/mc/RunIIAutumn18MiniAOD/BuToKMuMu_probefilter_SoftQCDnonD_TuneCP5_13TeV-pythia8-evtgen/MINIAODSIM/PUPoissonAve20_BParking_102X_upgrade2018_realistic_v15-v2/40000/43F82832-A9A2-5D45-950F-BC7D8DAC9C9B.root',
]
elif 'local' in options.tag:
options.inputFiles = [
'file:/afs/cern.ch/user/c/castilla/private/BParking/CMSSW_10_2_15/src/PhysicsTools/BParkingNano/input.root'
]
elif 'debug' in options.tag:
options.inputFiles = [
'file:/afs/cern.ch/user/c/castilla/private/BParking/CMSSW_10_2_15/src/PhysicsTools/BParkingNano/BtoKmm/pickevents.root'
]
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/RunIIFall17MiniAOD/ZprimeToWWToWlepWhad_narrow_M-3000_TuneCP5_13TeV-madgraph/MINIAODSIM/94X_mc2017_realistic_v10-v1/20000/3E25D208-8205-E811-8858-3417EBE64426.root'
options.maxEvents = -1
options.parseArguments()
process = cms.Process("PATtest")
## MessageLogger
process.load("FWCore.MessageLogger.MessageLogger_cfi")
process.MessageLogger.cerr.FwkReport.reportEvery = 100
## 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))
## Geometry and Detector Conditions (needed for a few patTuple production steps)
process.load("Configuration.Geometry.GeometryRecoDB_cff")
process.load(
"Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.outputFile = 'events.root'
options.inputFiles = '/store/data/Run2016G/DoubleMuon/MINIAOD/23Sep2016-v1/100000/0A30F7A9-ED8F-E611-91F1-008CFA1C6564.root' # ReReco
options.maxEvents = 1
options.parseArguments()
process = cms.Process("PickEvent")
process.source = cms.Source ("PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles)
)
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32 (options.maxEvents)
)
process.Out = cms.OutputModule("PoolOutputModule",
fileName = cms.untracked.string (options.outputFile)
)
process.end = cms.EndPath(process.Out)
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.outputFile = 'genTree.root'
options.inputFiles= '/store/user/dntaylor/HPlusPlusHMinusMinusHTo4L_M-500_13TeV-pythia8/RunIIFall15MiniAODv2_MINIAODSIM/160210_132739/0000/dblh_m500_13tev_miniAODv2_1.root'
options.maxEvents = -1
options.parseArguments()
process = cms.Process("GenNtuple")
process.load("FWCore.MessageService.MessageLogger_cfi")
process.MessageLogger.cerr.FwkReport.reportEvery = 100
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(options.maxEvents) )
process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles)
)
process.TFileService = cms.Service("TFileService",
fileName = cms.string(options.outputFile),
)
process.load("DevTools.GenNtuplizer.GenTree_cfi")
import ROOT
ROOT.gROOT.SetBatch(True)
sys.argv = oldargv
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('python')
# RooFit
ROOT.gSystem.Load("libRooFit.so")
ROOT.gSystem.Load("libRooFitCore.so")
ROOT.gROOT.SetStyle("Plain")
ROOT.gSystem.SetIncludePath("-I$ROOFITSYS/include/")
#default options
#options.inputFiles="../step3_inMINIAODSIM.root"
options.inputFiles = "files.txt"
options.outputFile = "result/pf_eta_in_jets.root"
options.maxEvents = -1
#overwrite if given any command line arguments
options.parseArguments()
#in case of txt input file, read the information from txt
li_f = []
for iF, F in enumerate(options.inputFiles):
print F
if F.split('.')[-1] == "txt":
options.inputFiles.pop(iF)
with open(F) as f:
li_f += f.read().splitlines()
options.inputFiles += li_f
print "analyzing files:"
if not options.inputFiles:
if options.inputDirectory is not None:
listDirectory(options.inputDirectory, options.inputFiles, options.maxFiles)
if True:
names = [f for f in options.datasetName.split('/') if f]
name = '%s-%s-%s.root' % (names[0],names[1],names[-1])
options.outputFile = os.path.join(options.outputDirectory,name)
files = getFiles(
[options.datasetName],
'wreece',
'susy_tree_CMG_[0-9]+.root'
)
if options.maxFiles > 0:
options.inputFiles = files[0:options.maxFiles]
else:
options.inputFiles = files
rt.gROOT.ProcessLine("""
struct Variables{\
Double_t mR;\
Double_t R;\
Double_t Rsq;\
Double_t maxTCHE;\
Double_t maxSSVP;\
Double_t nextTCHE;\
Double_t nextSSVP;\
Double_t met;\
Double_t jet1Pt;\
Double_t jet2Pt;\
import FWCore.ParameterSet.Config as cms
process = cms.Process("TEST")
runOnMC = True
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.parseArguments()
# Modify the defaults
if not options.inputFiles:
options.inputFiles = ['file:rawsim_numEvent100.root']
if options.outputFile == "output.root":
options.outputFile = "test_ntuple.root"
## 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 ),
#SkipEvent = cms.untracked.vstring('ProductNotFound')
)
## 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) )
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:
print "\n\nConfig ERROR: no inputFiles specified\n\n"
import FWCore.Utilities.FileUtils as FileUtils
from FWCore.ParameterSet.VarParsing import VarParsing
from flashgg.MetaData.samples_utils import SamplesManager
from flashgg.Systematics.SystematicDumperDefaultVariables import defaultVariables,defaultHistograms
## CMD LINE OPTIONS ##
options = VarParsing('analysis')
print options
processId = "tth"
targetMass = 130.
processId = "%s_%i" % (processId,int(targetMass))
# maxEvents is the max number of events processed of each file, not globally
options.maxEvents = -1
options.inputFiles = "file:myTagOutputFile_%s.root" % processId
options.outputFile = "ValidationTagsDump_%s.root" % processId
options.parseArguments()
## I/O SETUP ##
process = cms.Process("ValidationTagsDumper")
process.load("FWCore.MessageService.MessageLogger_cfi")
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(options.maxEvents))
process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32( 1 )
process.source = cms.Source ("PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles))
process.flashggSystTagMerger = cms.EDProducer("TagMerger",src=cms.VInputTag("flashggTagSorter"))
process.TFileService = cms.Service("TFileService",
fileName = cms.string(options.outputFile))
# Auto generated configuration file
# using:
# Revision: 1.20
# Source: /local/reps/CMSSW/CMSSW/Configuration/Applications/python/ConfigBuilder.py,v
# with command line options: Configuration/GenProduction/python/ThirteenTeV/Hadronizer_Tune4C_13TeV_generic_LHE_pythia8_Tauola_cff.py --filein file:pLHE_EDM_files/wpz_all_10000_000_Fall13pLHE.root --fileout file:pLHE_EDM_files/wpz_all_10000_000_Fall13pLHE_BG2-Fall13.root --mc --eventcontent RAWSIM --customise SLHCUpgradeSimulations/Configuration/postLS1Customs.customisePostLS1,Configuration/DataProcessing/Utils.addMonitoring --datatier GEN-SIM --conditions POSTLS162_V1::All --step GEN,SIM --magField 38T_PostLS1 --geometry Extended2015 --python_filename B2G-Fall13-00061_1_cfg.py --no_exec -n -1
import FWCore.ParameterSet.Config as cms
process = cms.Process('SIM')
# parse variables from cmsRun
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
options.inputFiles = "file:input.root"
options.outputFile = "step1.root"
options.parseArguments()
# import of standard configurations
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.Geometry.GeometryExtended2015Reco_cff')
process.load('Configuration.Geometry.GeometryExtended2015_cff')
process.load('Configuration.StandardSequences.MagneticField_38T_PostLS1_cff')
process.load('Configuration.StandardSequences.Generator_cff')
process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic8TeVCollision_cfi')
process.load('GeneratorInterface.Core.genFilterSummary_cff')
process.load('Configuration.StandardSequences.SimIdeal_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
Optional arguments:
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 = '/store/data/Run2012C/ZeroBias1/RAW/v1/000/198/609/04F0C0E3-72CA-E111-A802-003048F117EC.root'
#'/store/user/tapas/ETauSkim/skim_12_1_erV.root'
#options.inputFiles = '/store/user/tapas/2012-08-01-CRAB_ZEESkim/skim_10_1_wd2.root'
options.outputFile = "tAvgPU_rate_tree.root"
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")
inputFilesData = [
'/store/data/Run2015D/SingleMuon/MINIAOD/05Oct2015-v1/30000/DAF6F553-986F-E511-AB07-002618943982.root',
'/store/data/Run2015D/SingleMuon/MINIAOD/05Oct2015-v1/30000/DC242338-986F-E511-9196-00261894394B.root',
'/store/data/Run2015D/SingleMuon/MINIAOD/05Oct2015-v1/30000/DC4D3D93-906F-E511-BEAD-0026189438F9.root',
]
# file dataset=/JPsiToMuMu_Pt20to100-pythia8-gun/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v2/MINIAODSIM
inputFilesMC = [
'/store/mc/RunIISpring15MiniAODv2/JPsiToMuMu_Pt20to100-pythia8-gun/MINIAODSIM/74X_mcRun2_asymptotic_v2-v2/30000/085258CF-CA76-E511-AF11-00237DF20360.root',
'/store/mc/RunIISpring15MiniAODv2/JPsiToMuMu_Pt20to100-pythia8-gun/MINIAODSIM/74X_mcRun2_asymptotic_v2-v2/30000/2C190717-CB76-E511-911C-00237DF25430.root',
'/store/mc/RunIISpring15MiniAODv2/JPsiToMuMu_Pt20to100-pythia8-gun/MINIAODSIM/74X_mcRun2_asymptotic_v2-v2/30000/48EAA0D6-CA76-E511-82F0-00237DF25430.root',
]
if len(options.inputFiles) is 0:
if options.isMC:
options.inputFiles = inputFilesMC
else:
options.inputFiles = inputFilesData
if len(options.outputFile) is 0:
if options.isMC:
options.outputFile = 'TnPTree_mc.root'
else:
options.outputFile = 'TnPTree_data.root'
process.source = cms.Source(
"PoolSource",
fileNames=cms.untracked.vstring(options.inputFiles),
)
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.Geometry.GeometrySimDB_cff')
process.load('Configuration.StandardSequences.MagneticField_38T_PostLS1_cff')
process.load('Configuration.StandardSequences.Generator_cff')
process.load('IOMC.EventVertexGenerators.VtxSmearedNominalCollision2015_cfi')
process.load('GeneratorInterface.Core.genFilterSummary_cff')
process.load('Configuration.StandardSequences.SimIdeal_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
# parse variables from cmsRun
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
options.inputFiles = "file:test.root"
options.outputFile = "step1.root"
options.parseArguments()
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(-1)
)
# Input source
process.source = cms.Source("PoolSource",
secondaryFileNames = cms.untracked.vstring(),
fileNames = cms.untracked.vstring(options.inputFiles),
inputCommands = cms.untracked.vstring('keep *',
'drop LHEXMLStringProduct_*_*_*'),
dropDescendantsOfDroppedBranches = cms.untracked.bool(False)
)
process.options = cms.untracked.PSet(
# cmsRun Python Configuration file
# for more info on the syntax, see https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideAboutPythonConfigFile
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
# default arguments
options.outputFile = 'test0.root'
options.inputFiles = 'file:/scratch/mratti/ECALDPG_test_samples/CMSSW_10_0_0_RelValZEE_13_GEN-SIM-DIGI-RAW-RECO_PU25ns_100X_upgrade2018_realistic_v7_HS-v1__A_FILE.root'
# if not using file: it will look in the local storage element
options.maxEvents = -1 # -1 means all events, maxEvents considers the total over files considered
options.register(
'anaName',
'DoublePhoton', # default value # DoubleElectron
VarParsing.multiplicity.singleton, # singleton or list
VarParsing.varType.string, # string, int, or float
'Analysis channel') # description
options.register(
'conditions',
'102X_upgrade2018_realistic_EcalAging_mid2021_235fb_v1', #
VarParsing.multiplicity.singleton, # singleton or list
VarParsing.varType.string, # string, int, or float
'global tag of input sample') # description
options.parseArguments()
# Define the process and set some options for it
process = cms.Process('Validation')
# initialize MessageLogger and output report
## Set up command line options
options = VarParsing ('analysis')
options.register('runOnGenOrAODsim', False, VarParsing.multiplicity.singleton, VarParsing.varType.bool, "GEN SIM")
options.register( "skipEvents", 0, VarParsing.multiplicity.singleton, VarParsing.varType.int, "Number of events to skip" )
options.parseArguments()
## 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 = ''
)
#needs to be disabled for crab to work, otherwise get configuration errors
#validate using edmConfigHash
if hasattr(sys, "argv") and not "-cfg" in sys.argv:
options.parseArguments()
process = cms.Process("TTH")
process.load("Configuration.Geometry.GeometryIdeal_cff")
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
process.load("Configuration.StandardSequences.MagneticField_cff")
process.load("PhysicsTools.PatAlgos.patSequences_cff")
process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
if len(options.inputFiles) == 0:
options.inputFiles = ["/store/data/Run2012A/SingleElectron/AOD/22Jan2013-v1/20000/00738DE4-8673-E211-8741-002481E15204.root"]
process.source = cms.Source(
"PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles),
#required for taus: https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuidePFTauID#53X_recommendation_for_Run_I_ana
dropDescendantsOfDroppedBranches = cms.untracked.bool(False),
inputCommands = cms.untracked.vstring(
'keep *',
'drop recoPFTaus_*_*_*'
)
)
process.out = cms.OutputModule(
"PoolOutputModule",
#validate using edmConfigHash
if hasattr(sys, "argv") and not "-cfg" in sys.argv:
options.parseArguments()
process = cms.Process("TTH")
process.load("Configuration.Geometry.GeometryIdeal_cff")
process.load(
"Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
process.load("Configuration.StandardSequences.MagneticField_cff")
process.load("PhysicsTools.PatAlgos.patSequences_cff")
process.options = cms.untracked.PSet(wantSummary=cms.untracked.bool(True))
if len(options.inputFiles) == 0:
#options.inputFiles = ["/store/mc/Summer12_DR53X/TTH_HToBB_M-125_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/0000/00FA9388-81FC-E111-A80D-00215E2217BE.root"]
options.inputFiles = [
"/store/mc/Summer12_DR53X/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball/AODSIM/PU_S10_START53_V7A-v1/0000/00037C53-AAD1-E111-B1BE-003048D45F38.root"
]
process.source = cms.Source(
"PoolSource",
fileNames=cms.untracked.vstring(options.inputFiles),
#required for taus: https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuidePFTauID#53X_recommendation_for_Run_I_ana
dropDescendantsOfDroppedBranches=cms.untracked.bool(False),
inputCommands=cms.untracked.vstring('keep *', 'drop recoPFTaus_*_*_*'))
process.out = cms.OutputModule(
"PoolOutputModule",
fileName=cms.untracked.string(options.outputFile),
SelectEvents=cms.untracked.PSet(SelectEvents=cms.vstring()),
outputCommands=cms.untracked.vstring(["drop *"]))
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
#input cmsRun options
options = VarParsing ('analysis')
with open('files2') as f:
options.inputFiles = f.readlines()
#options.inputFiles ="/store/mc/RunIISpring16MiniAODv2/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_miniAODv2_v0_ext1-v1/50000/127657A5-4E1C-E611-A5B2-001E672486B0.root"
#options.inputFiles = "/store/mc/RunIISpring16MiniAODv1/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_v3_ext1-v1/20000/0852310B-B6FD-E511-96B1-0002C94CD0B8.root" #smaller set for testing
#options.inputFiles="/store/mc/RunIIFall15MiniAODv2/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/MINIAODSIM/PU25nsData2015v1_76X_mcRun2_asymptotic_v12-v1/70000/068CC5B6-DDB8-E511-82EB-003048FFD756.root"
options.outputFile = "MiniAOD_effi_80x_ZtoTT.root"
options.parseArguments()
#name the process
process = cms.Process("TreeProducerFromMiniAOD")
process.load('FWCore/MessageService/MessageLogger_cfi')
process.MessageLogger.cerr.FwkReport.reportEvery = 1000;
process.MessageLogger.cerr.threshold = cms.untracked.string('INFO')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
#50 ns global tag for MC replace with 'GR_P_V56' for prompt reco. https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideFrontierConditions#Prompt_reconstruction_Global_Tag
from Configuration.AlCa.GlobalTag import GlobalTag
#Make sure Global Tag mathes input file type
#process.GlobalTag = GlobalTag(process.GlobalTag, '76X_mcRun2_asymptotic_RunIIFall15DR76_v1', '')
process.GlobalTag = GlobalTag(process.GlobalTag, '80X_mcRun2_asymptotic_v6', '')
#how many events to run over
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(1000000)
)
import FWCore.ParameterSet.Config as cms
###
# Defaults
###
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
options.inputFiles = "/store/mc/Summer12/GluGluToHToZZTo4L_M-125_14TeV-powheg-pythia6/GEN-SIM-RAW-RECO/PU50bx25_POSTLS161_V12-v1/10000/0004A91B-754C-E211-9FEA-00304867BEDE.root"
options.outputFile = "POSTLS1_PU50bx25_Buffer.root"
options.parseArguments()
process = cms.Process("ReRecoBuffer")
# import of standard configurations
process.load('Configuration.StandardSequences.Services_cff')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.EventContent.EventContent_cff')
process.load('Configuration.Geometry.GeometryIdeal_cff.py')
process.load('Configuration.StandardSequences.MagneticField_38T_PostLS1_cff')
process.load('Configuration.StandardSequences.Reconstruction_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
process.source = cms.Source("PoolSource",
# replace 'myfile.root' with the source file you want to use
fileNames = cms.untracked.vstring(
'file:/afs/cern.ch/work/g/gauzinge/public/BeamHalo/',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"The FLUKA file directory")
options.parseArguments()
#specify input
inputPath = "/afs/cern.ch/work/g/gauzinge/public/BeamHalo"
#inputPath = "/afs/cern.ch/work/g/gauzinge/public/BeamGas/cms_conditioned_carbon"
#inputPath = "/afs/cern.ch/work/g/gauzinge/public/BeamGas/cms_conditioned_hydrogen"
#inputPath = "/afs/cern.ch/work/g/gauzinge/public/BeamGas/cms_conditioned_oxygen"
#inputPath = "/afs/cern.ch/work/p/pkicsiny/private/cmssw/CMSSW_11_2_0_pre6/src/GeneratorInterface/BeamHaloGenerator/input/"
options.inputFiles = [
inputPath + "/" + f for f in os.listdir(inputPath) if f[:3] == "run"
][0]
print(len(options.inputFiles))
print(options.inputFiles)
#count number of events in all input files
print("Number of events in input files: {}".format(
sum([
len(
set([
int(line.split()[0]) for line in open(f)
if line.split()[0].isdigit()
])) for f in options.inputFiles
])))
print("Number of files: {}".format(len(options.inputFiles)))
#specify output name
import FWCore.ParameterSet.Config as cms
###
# Defaults
###
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
#options.inputFiles = "file:GEM2019Upg14_CFEBBufferOverload_0p01.root"
options.inputFiles = "file:/hdfs/store/user/dntaylor/BufferOverload_HZZ_0p10_CFEB_v8-CFEBBufferOverload_cfg/CFEBBufferOverload_cfg-GEM2019Upg14_9.root"
options.outputFile = "CFEBBufferOverload_StandAloneMuon.root"
options.parseArguments()
###
# Run Skim
###
process = cms.Process("StandAloneMuonSkim")
process.load('Configuration.StandardSequences.Services_cff')
process.load('FWCore.MessageService.MessageLogger_cfi')
#process.load('Configuration.Geometry.GeometryExtended2019Reco_cff')
#process.load('Configuration.StandardSequences.MagneticField_38T_PostLS1_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
from Configuration.AlCa.GlobalTag import GlobalTag
process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:upgrade2019', '')
process.options = cms.untracked.PSet(
wantSummary = cms.untracked.bool(True),
# SkipEvent = cms.untracked.vstring('ProductNotFound')
)
import FWCore.ParameterSet.Config as cms
import os
from FWCore.ParameterSet.VarParsing import VarParsing
#input cmsRun options
options = VarParsing ('analysis')
options.inputFiles = '/store/mc/RunIISpring15DR74/WJetsToLNu_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/MINIAODSIM/Asympt50ns_MCRUN2_74_V9A-v1/80000/06293CB4-CEFD-E411-8BE4-0CC47A13D16A.root'
options.outputFile = "MiniAOD_FR_WJets.root"
options.parseArguments()
#name the process
process = cms.Process("TreeProducerFromMiniAOD")
#Make the framework shutup
process.load('FWCore/MessageService/MessageLogger_cfi')
process.MessageLogger.cerr.FwkReport.reportEvery = 100;
process.MessageLogger.cerr.threshold = cms.untracked.string('INFO')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff')
#whether running MC or data
options.register(
'isMC',
1,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
'Set to 1 for simulated samples - updates GT, emulates HCAL TPGs.'
)
process = cms.Process("Demo")
# Necessary so we can import Main_cfg.py from ordinary python
# scripts and extract the fatjets
try:
options.parseArguments()
except:
pass
if len(options.inputFiles)==0:
#options.inputFiles = cms.untracked.vstring(["/store/mc/Spring14miniaod/TTbarH_HToBB_M-125_13TeV_pythia6/MINIAODSIM/PU20bx25_POSTLS170_V5-v1/00000/0E97DD3E-2209-E411-8A04-003048945312.root"])
options.inputFiles = cms.untracked.vstring([
'/store/mc/Spring14miniaod/TTbarH_HToBB_M-125_13TeV_pythia6/MINIAODSIM/PU20bx25_POSTLS170_V5-v1/00000/0E97DD3E-2209-E411-8A04-003048945312.root',
'/store/mc/Spring14miniaod/TTbarH_HToBB_M-125_13TeV_pythia6/MINIAODSIM/PU20bx25_POSTLS170_V5-v1/00000/D0242854-2209-E411-B361-003048C559CE.root'
])
#enable debugging printout
if "TTH_DEBUG" in os.environ:
process.load("FWCore.MessageLogger.MessageLogger_cfi")
process.MessageLogger = cms.Service("MessageLogger",
destinations=cms.untracked.vstring('cout', 'debug'),
debugModules=cms.untracked.vstring('*'),
cout=cms.untracked.PSet(threshold=cms.untracked.string('INFO')),
debug=cms.untracked.PSet(threshold=cms.untracked.string('DEBUG')),
)
else:
process.load("FWCore.MessageService.MessageLogger_cfi")
import FWCore.ParameterSet.Config as cms
import os
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.inputFiles = 'root://cmsxrootd.fnal.gov//store/mc/RunIISpring15DR74/GluGluHToTauTau_M125_13TeV_powheg_pythia8/AODSIM/Asympt25ns_MCRUN2_74_V9-v1/30000/0C4BFABC-362E-E511-8591-0CC47A4DED1A.root'
options.outputFile = "testTau_fakeRate.root"
options.parseArguments()
process = cms.Process("TreeProducerFromAOD")
#from FWCore.ParameterSet.VarParsing import VarParsing
#options = VarParsing ('analysis')
#options.inputFiles ='file:/hdfs/store/mc/Phys14DR/WJetsToLNu_13TeV-madgraph-pythia8-tauola/AODSIM/PU20bx25_PHYS14_25_V1-v1/00000/28C4E0C1-7F6F-E411-AE20-0025905B85EE.root'
#options.outputFiles ='outputReRunNonStdStrip.root'
process.load('FWCore/MessageService/MessageLogger_cfi')
#process.MessageLogger.cerr.FwkReport.reportEvery = 100
#process.MessageLogger.cerr.threshold = cms.untracked.string('INFO')
process.load("Configuration.Geometry.GeometryDB_cff")
process.load(
'Configuration/StandardSequences/FrontierConditions_GlobalTag_cff')
process.load(
'Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff'
)
process.load("Configuration.StandardSequences.MagneticField_cff")
from Configuration.AlCa.GlobalTag_condDBv2 import GlobalTag as customiseGlobalTag
process.GlobalTag = customiseGlobalTag(
process.GlobalTag,
globaltag='MCRUN2_74_V9A',
Optional arguments:
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 = '/store/user/tapas/ETauSkimRun2012D/skim_101_1_MBW.root'
options.outputFile = "uct_efficiency_tree.root"
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, VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"If 1, enable production of Stage1B trees")
options.register("stage1", 1, VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"If 0, *disable* production of Stage1 trees")
options.register('ecalCalib', 'CALIB_V4', VarParsing.multiplicity.singleton,
VarParsing.varType.string,
'Can be CALIB_V1, CALIB_V3, or CALIB_V4')
options.register('eicCardHcalOnly', 0, VarParsing.multiplicity.singleton,
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
from Configuration.Eras.Era_Run3_cff import Run3
from Configuration.Eras.Era_Run2_2018_cff import Run2_2018
options = VarParsing('analysis')
options.register ("run3", True, VarParsing.multiplicity.singleton, VarParsing.varType.bool)
options.register ("mc", False, VarParsing.multiplicity.singleton, VarParsing.varType.bool)
options.register ("useB904Data", False, VarParsing.multiplicity.singleton, VarParsing.varType.bool)
options.register ("useGEMs", False, VarParsing.multiplicity.singleton, VarParsing.varType.bool)
options.parseArguments()
options.inputFiles = "file:step_DQM.root"
process_era = Run3
if not options.run3:
process_era = Run2_2018
process = cms.Process("DQMClient", process_era)
process.load("Configuration/StandardSequences/GeometryRecoDB_cff")
process.load("Configuration/StandardSequences/MagneticField_cff")
process.load("Configuration/StandardSequences/FrontierConditions_GlobalTag_cff")
process.load('Configuration.StandardSequences.DQMSaverAtRunEnd_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load('Configuration.EventContent.EventContent_cff')
process.load("DQM.L1TMonitorClient.L1TdeCSCTPGClient_cfi")
process.load("DQM.L1TMonitorClient.L1TdeGEMTPGClient_cfi")
process.load('DQMOffline.Configuration.DQMOfflineMC_cff')
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(1)
)
"The maximum number of files to read")
options.parseArguments()
if options.inputDirectory is not None:
listDirectory(options.inputDirectory, options.inputFiles,
options.maxFiles)
if True:
names = [f for f in options.datasetName.split('/') if f]
name = '%s-%s-%s.root' % (names[0], names[1], names[-1])
options.outputFile = os.path.join(options.outputDirectory, name)
files = getFiles([options.datasetName], 'wreece',
'susy_tree_CMG_[0-9]+.root')
if options.maxFiles > 0:
options.inputFiles = files[0:options.maxFiles]
else:
options.inputFiles = files
rt.gROOT.ProcessLine("""
struct Variables{\
Double_t mR;\
Double_t R;\
Double_t Rsq;\
Double_t maxTCHE;\
Double_t maxSSVP;\
Double_t nextTCHE;\
Double_t nextSSVP;\
Double_t met;\
Double_t jet1Pt;\
Double_t jet2Pt;\
import FWCore.ParameterSet.Config as cms
import os
from FWCore.ParameterSet.VarParsing import VarParsing
#input cmsRun options
options = VarParsing ('analysis')
options.inputFiles =" /store/mc/RunIISpring16MiniAODv2/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/MINIAODSIM/PUSpring16_80X_mcRun2_asymptotic_2016_miniAODv2_v0_ext1-v1/50000/127657A5-4E1C-E611-A5B2-001E672486B0.root"
#options.inputFiles="/store/mc/RunIIFall15MiniAODv2/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/MINIAODSIM/PU25nsData2015v1_76X_mcRun2_asymptotic_v12-v1/70000/068CC5B6-DDB8-E511-82EB-003048FFD756.root"
options.outputFile = "MiniAOD_FR_ZToEE.root"
options.parseArguments()
#name the process
process = cms.Process("TreeProducerFromMiniAOD")
#Make the framework shutup
process.load('FWCore/MessageService/MessageLogger_cfi')
process.MessageLogger.cerr.FwkReport.reportEvery = 100;
process.MessageLogger.cerr.threshold = cms.untracked.string('INFO')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
from Configuration.AlCa.GlobalTag import GlobalTag
#50 ns global tag for MC replace with 'GR_P_V56' for prompt reco. https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideFrontierConditions#Prompt_reconstruction_Global_Tag
#process.GlobalTag = GlobalTag(process.GlobalTag, '76X_mcRun2_asymptotic_RunIIFall15DR76_v1', '')
process.GlobalTag = GlobalTag(process.GlobalTag, '80X_mcRun2_asymptotic_2016_v3', '')
#how many events to run over
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(-1)
)
#input files
Optional arguments:
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 = '/store/data/Run2012C/ZeroBias1/RAW/v1/000/198/609/04F0C0E3-72CA-E111-A802-003048F117EC.root'
#'/store/user/tapas/ETauSkim/skim_12_1_erV.root'
#options.inputFiles = '/store/user/tapas/2012-08-01-CRAB_ZEESkim/skim_10_1_wd2.root'
options.outputFile = "tAvgPU_rate_tree.root"
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, VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"If 1, enable production of Stage1B trees")
options.register("stage1", 1, VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"If 0, *disable* production of Stage1 trees")
options.register('ecalCalib', 'CALIB_V4', VarParsing.multiplicity.singleton,
VarParsing.varType.string,
#!/usr/bin/env cmsRun
import FWCore.ParameterSet.Config as cms
import FWCore.Utilities.FileUtils as FileUtils
from FWCore.ParameterSet.VarParsing import VarParsing
## CMD LINE OPTIONS ##
options = VarParsing('analysis')
# maxEvents is the max number of events processed of each file, not globally
options.maxEvents = -1
options.inputFiles = "file:diphotonsMicroAOD.root"
options.outputFile = "quickDump.root"
options.parseArguments()
process = cms.Process("Analysis")
process.load("FWCore.MessageService.MessageLogger_cfi")
# process.source = cms.Source ("PoolSource",
# fileNames = cms.untracked.vstring(options.inputFiles))
readFiles = cms.untracked.vstring()
process.source = cms.Source ("PoolSource",fileNames = readFiles)
readFiles.extend( [
'/store/user/spigazzi/flashgg/diphotonsIsoStudies0T/1_2_0-31-g16f7065/GJet_Pt-15to6000_TuneCUETP8M1_Flat_13TeV_pythia8/diphotonsIsoStudies0T-1_2_0-31-g16f7065-v0-RunIIFall15MiniAODv2-PU25nsData2015v1_magnetOff_76X_mcRun2_asymptotic_v12-v1/160119_091710/0000/myMicroAODOutputFile_1.root',
'/store/user/spigazzi/flashgg/diphotonsIsoStudies0T/1_2_0-31-g16f7065/GJet_Pt-15to6000_TuneCUETP8M1_Flat_13TeV_pythia8/diphotonsIsoStudies0T-1_2_0-31-g16f7065-v0-RunIIFall15MiniAODv2-PU25nsData2015v1_magnetOff_76X_mcRun2_asymptotic_v12-v1/160119_091710/0000/myMicroAODOutputFile_105.root',
'/store/user/spigazzi/flashgg/diphotonsIsoStudies0T/1_2_0-31-g16f7065/GJet_Pt-15to6000_TuneCUETP8M1_Flat_13TeV_pythia8/diphotonsIsoStudies0T-1_2_0-31-g16f7065-v0-RunIIFall15MiniAODv2-PU25nsData2015v1_magnetOff_76X_mcRun2_asymptotic_v12-v1/160119_091710/0000/myMicroAODOutputFile_106.root',
'/store/user/spigazzi/flashgg/diphotonsIsoStudies0T/1_2_0-31-g16f7065/GJet_Pt-15to6000_TuneCUETP8M1_Flat_13TeV_pythia8/diphotonsIsoStudies0T-1_2_0-31-g16f7065-v0-RunIIFall15MiniAODv2-PU25nsData2015v1_magnetOff_76X_mcRun2_asymptotic_v12-v1/160119_091710/0000/myMicroAODOutputFile_107.root',
'/store/user/spigazzi/flashgg/diphotonsIsoStudies0T/1_2_0-31-g16f7065/GJet_Pt-15to6000_TuneCUETP8M1_Flat_13TeV_pythia8/diphotonsIsoStudies0T-1_2_0-31-g16f7065-v0-RunIIFall15MiniAODv2-PU25nsData2015v1_magnetOff_76X_mcRun2_asymptotic_v12-v1/160119_091710/0000/myMicroAODOutputFile_109.root',
'/store/user/spigazzi/flashgg/diphotonsIsoStudies0T/1_2_0-31-g16f7065/GJet_Pt-15to6000_TuneCUETP8M1_Flat_13TeV_pythia8/diphotonsIsoStudies0T-1_2_0-31-g16f7065-v0-RunIIFall15MiniAODv2-PU25nsData2015v1_magnetOff_76X_mcRun2_asymptotic_v12-v1/160119_091710/0000/myMicroAODOutputFile_11.root',
#!/usr/bin/env cmsRun
import FWCore.ParameterSet.Config as cms
import FWCore.Utilities.FileUtils as FileUtils
from FWCore.ParameterSet.VarParsing import VarParsing
from flashgg.MetaData.samples_utils import SamplesManager
## CMD LINE OPTIONS ##
options = VarParsing('analysis')
print options
# maxEvents is the max number of events processed of each file, not globally
options.maxEvents = -1
options.inputFiles = "file:myTagOutputFile.root"
options.outputFile = "VBFTagsDump.root"
options.parseArguments()
## I/O SETUP ##
process = cms.Process("VBFTagsDumper")
process.load("FWCore.MessageService.MessageLogger_cfi")
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(-1))
process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32( 1 )
process.source = cms.Source ("PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles))
if options.maxEvents > 0:
process.source.eventsToProcess = cms.untracked.VEventRange('1:1-1:'+str(options.maxEvents))
process.TFileService = cms.Service("TFileService",
#!/usr/bin/env cmsRun
import FWCore.ParameterSet.Config as cms
import FWCore.Utilities.FileUtils as FileUtils
from FWCore.ParameterSet.VarParsing import VarParsing
from flashgg.MetaData.samples_utils import SamplesManager
## CMD LINE OPTIONS ##
options = VarParsing('analysis')
# maxEvents is the max number of events processed of each file, not globally
options.maxEvents = -1
options.inputFiles = "../../myTagOutputFile.root"
options.outputFile = "ValidationTagsDump.root"
options.parseArguments()
## I/O SETUP ##
process = cms.Process("ValidationTagsDumper")
process.load("FWCore.MessageService.MessageLogger_cfi")
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(-1))
process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32( 1 )
process.source = cms.Source ("PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles),
eventsToProcess = cms.untracked.VEventRange('1:1-1:'+str(options.maxEvents)))
process.TFileService = cms.Service("TFileService",
fileName = cms.string(options.outputFile))
## TAGS DUMPERS ##
# Auto generated configuration file
# using:
# Revision: 1.19
# Source: /local/reps/CMSSW/CMSSW/Configuration/Applications/python/ConfigBuilder.py,v
# with command line options: step1 --filein file:paceholder_Fall13_root.root --fileout file:placeholder_Sprint14dr_1.root.root --mc --eventcontent MINIAODSIM --runUnscheduled --datatier MINIAODSIM --conditions PLS170_V7AN1::All --step PAT --python_filename EXO-Spring14miniaod_cfg.py --no_exec --customise Configuration/DataProcessing/Utils.addMonitoring -n -1
import FWCore.ParameterSet.Config as cms
process = cms.Process('PAT')
# parse variables from cmsRun
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
options.inputFiles = "file:input_placeholder.root"
options.outputFile = "file:output_placeholder.root"
options.parseArguments()
# import of standard configurations
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.PATMC_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(-1)
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing("analysis")
options.inputFiles = 'file:N4000_0000_seed1001.root', 'file:N4000_0001_seed1002.root'
options.outputFile = 'substructurepacks.root'
options.maxEvents = -1
options.register(
'coneSize',
0.8,
VarParsing.multiplicity.singleton,
VarParsing.varType.float,
"Cone size with which to cluster the genJets"
)
options.parseArguments()
from Configuration.StandardSequences.Eras import eras
process = cms.Process('SoftDropGenJets',eras.Run2_2016)
# import of standard configurations
process.load('Configuration.StandardSequences.Services_cff')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
process.load('Configuration.StandardSequences.MagneticField_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.load('Configuration.EventContent.EventContent_cff')
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(options.maxEvents)
)
)
#needs to be disabled for crab to work, otherwise get configuration errors
#validate using edmConfigHash
if hasattr(sys, "argv") and not "-cfg" in sys.argv:
options.parseArguments()
process = cms.Process("TTH")
process.load("Configuration.Geometry.GeometryIdeal_cff")
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
process.load("Configuration.StandardSequences.MagneticField_cff")
process.load("PhysicsTools.PatAlgos.patSequences_cff")
process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
if len(options.inputFiles) == 0:
options.inputFiles = ["/store/mc/Summer12_DR53X/TTH_HToBB_M-125_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/0000/00FA9388-81FC-E111-A80D-00215E2217BE.root"]
process.source = cms.Source(
"PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles),
#required for taus: https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuidePFTauID#53X_recommendation_for_Run_I_ana
dropDescendantsOfDroppedBranches = cms.untracked.bool(False),
inputCommands = cms.untracked.vstring(
'keep *',
'drop recoPFTaus_*_*_*'
)
)
process.out = cms.OutputModule(
"PoolOutputModule",
# import ROOT in batch mode
import sys
import time
oldargv = sys.argv[:]
sys.argv = ['-b-']
import ROOT
ROOT.gROOT.SetBatch(True)
sys.argv = oldargv
PLOTRAW = True #Set True if you want only raw pt from met objects
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('python')
#default options
options.inputFiles = "../step3_inMINIAODSIM.root"
options.outputFile = "jetmetNtuples.root"
options.maxEvents = -1
#overwrite if given any command line arguments
options.parseArguments()
#in case of txt input file, read the information from txt
li_f = []
for iF, F in enumerate(options.inputFiles):
print F
if F.split('.')[-1] == "txt":
options.inputFiles.pop(iF)
with open(F) as f:
li_f += f.read().splitlines()
options.inputFiles += li_f
print "analyzing files:"
for F in options.inputFiles:
#!/usr/bin/env cmsRun
import FWCore.ParameterSet.Config as cms
import FWCore.Utilities.FileUtils as FileUtils
from FWCore.ParameterSet.VarParsing import VarParsing
from flashgg.MetaData.samples_utils import SamplesManager
## CMD LINE OPTIONS ##
options = VarParsing('analysis')
print options
# maxEvents is the max number of events processed of each file, not globally
options.maxEvents = -1
options.inputFiles = "file:/tmp/nancy/vhEtTagOutputFile_VHsignal.root"
options.outputFile = "VHEtTagsDump_signal.root"
options.parseArguments()
## I/O SETUP ##
process = cms.Process("VHEtTagsDumper")
process.load("FWCore.MessageService.MessageLogger_cfi")
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(-1))
process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32( 1 )
process.source = cms.Source ("PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles))
if options.maxEvents > 0:
process.source.eventsToProcess = cms.untracked.VEventRange('1:1-1:'+str(options.maxEvents))
process.TFileService = cms.Service("TFileService",
# import ROOT in batch mode
import sys
oldargv = sys.argv[:]
sys.argv = [ '-b-' ]
import ROOT
ROOT.gROOT.SetBatch(True)
sys.argv = oldargv
from collections import defaultdict
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('python')
#default options
options.inputFiles="root://cmsxrootd-site.fnal.gov//store/mc/RunIIAutumn18MiniAOD/QCD_Pt-15to7000_TuneCP5_Flat_13TeV_pythia8/MINIAODSIM/102X_upgrade2018_realistic_v15_ext1-v1/110000/8EF57574-29F7-DA4C-A5EC-8002BF2DEB96.root"
options.outputFile="jet.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 np
from jet_tree import *
# create an oput tree.
fout = ROOT.TFile (options.outputFile,"recreate")
t = ROOT.TTree ("Events","Events")
declare_branches(t)
Optional arguments:
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(
inputFilesData = [
'/store/data/Run2015D/DoubleMuon/MINIAOD/05Oct2015-v1/40000/529E5F9F-5F6F-E511-83C6-0026189438F4.root',
'/store/data/Run2015D/DoubleMuon/MINIAOD/05Oct2015-v1/40000/52A8BE9F-5F6F-E511-86A9-0025905AA9CC.root',
'/store/data/Run2015D/DoubleMuon/MINIAOD/05Oct2015-v1/40000/56CC0262-5E6F-E511-B9BD-0025905A613C.root',
]
# file dataset=/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/RunIISpring15MiniAODv2-74X_mcRun2_asymptotic_v2-v1/MINIAODSIM
inputFilesMC = [
'/store/mc/RunIISpring15MiniAODv2/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/MINIAODSIM/74X_mcRun2_asymptotic_v2-v1/50000/00759690-D16E-E511-B29E-00261894382D.root',
'/store/mc/RunIISpring15MiniAODv2/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/MINIAODSIM/74X_mcRun2_asymptotic_v2-v1/50000/00E88378-6F6F-E511-9D54-001E6757EAA4.root',
'/store/mc/RunIISpring15MiniAODv2/DYJetsToLL_M-50_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8/MINIAODSIM/74X_mcRun2_asymptotic_v2-v1/50000/02CD8A95-736F-E511-B76E-00266CFFBF34.root',
]
if len(options.inputFiles) is 0 :
if options.isMC :
options.inputFiles = inputFilesMC
else :
options.inputFiles = inputFilesData
if len(options.outputFile) is 0 :
if options.isMC :
options.outputFile = 'TnPTreeDZ_mc.root'
else :
options.outputFile = 'TnPTreeDZ_data.root'
process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles),
)
process.load('FWCore.MessageService.MessageLogger_cfi')
process.MessageLogger.cerr.FwkReport.reportEvery = 1000
# 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")
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
#input cmsRun options
options = VarParsing ('analysis')
options.inputFiles = '/store/mc/RunIISpring15DR74/ZZTo4L_13TeV_powheg_pythia8/MINIAODSIM/Asympt50ns_MCRUN2_74_V9A-v1/00000/2C0279E7-9B16-E511-B5FF-90B11C2ACF20.root'
#'/store/mc/RunIISpring15DR74/ZZTo4L_13TeV_powheg_pythia8/MINIAODSIM/Asympt50ns_MCRUN2_74_V9A-v1/00000/2C0279E7-9B16-E511-B5FF-90B11C2ACF20.root'
#'/store/mc/RunIISpring15DR74/ZprimeToTauTau_M_500_TuneCUETP8M1_tauola_13TeV_pythia8/MINIAODSIM/Asympt25ns_MCRUN2_74_V9-v1/60000/CAE0DF38-7D19-E511-9DA4-0025904A8EC8.root'
#options.inputFiles = '/store/mc/RunIISpring15DR74/SUSYGluGluToHToTauTau_M-250_TuneCUETP8M1_13TeV-pythia8/MINIAODSIM/Asympt25ns_MCRUN2_74_V9-v1/80000/6EBE4F55-4E03-E511-9FE6-0CC47A13D416.root'
options.outputFile = "MiniAOD_effi.root"
options.parseArguments()
#name the process
process = cms.Process("TreeProducerFromMiniAOD")
process.load('FWCore/MessageService/MessageLogger_cfi')
process.MessageLogger.cerr.FwkReport.reportEvery = 100;
process.MessageLogger.cerr.threshold = cms.untracked.string('INFO')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff')
#50 ns global tag for MC replace with 'GR_P_V56' for prompt reco. https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideFrontierConditions#Prompt_reconstruction_Global_Tag
process.GlobalTag.globaltag = 'MCRUN2_74_V9A'
#how many events to run over
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(-1)
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:
print "\n\nConfig ERROR: no inputFiles specified\n\n"
sys.exit()
options.register(
'nEventsSample', 100, VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Total number of events in dataset for event weight calculation.")
## 'maxEvents' is already registered by the Framework, changing default value
options.setDefault('maxEvents', 100)
options.parseArguments()
outName = options.outputFile
print "Default output name: " + outName
if "output" in outName: # if an input file name is specified, event weights can be determined
outName = "out_" + basename(options.inputFiles[0])
print "Output root file name: " + outName
else:
options.inputFiles = "file:GGJets_M-1000To2000_Pt-50_13TeV-sherpa.root"
# outName = "ExoDiphotonAnalyzer.root"
isMC = True
# data always has "Run201" in its filename
if "Run201" in outName:
isMC = False
# to avoid processing with an incorrect global tag, don't set a valid default
globalTag = 'notset'
jetLabel = "selectedUpdatedPatJetsUpdatedJEC"
# options for data
JEC = cms.vstring(['L1FastJet', 'L2Relative', 'L3Absolute', 'L2L3Residual'])
# necessary because in re-MINIAOD the process labels are "PAT", not "RECO"
isReMINIAOD = False
options.setDefault('maxEvents', 500)
options.setDefault('tag', '10215')
options.parseArguments()
globaltag = '102X_dataRun2_v11' if not options.isMC else '102X_upgrade2018_realistic_v15'
if options._beenSet['globalTag']:
globaltag = options.globalTag
extension = {False: 'data', True: 'mc'}
outputFileNANO = cms.untracked.string(
'_'.join(['BParkNANO', extension[options.isMC], options.tag]) + '.root')
outputFileFEVT = cms.untracked.string(
'_'.join(['BParkFullEvt', extension[options.isMC], options.tag]) + '.root')
if not options.inputFiles:
options.inputFiles = [
'file:/afs/cern.ch/user/d/dryu/BFrag/CMSSW_10_2_15_skim/src/PhysicsTools/BParkingNano/test/BdToKstarJpsi_ToKPiMuMu_probefilter_SoftQCDnonD_TuneCP5_13TeV-pythia8-evtgen_BF7A19E8-ED80-B841-A6EB-B583F36C1125.root'
]
# options.inputFiles = ['file:/afs/cern.ch/user/d/dryu/BFrag/CMSSW_10_2_15_skim/src/PhysicsTools/BParkingNano/test/BsToPhiJpsi_ToKKMuMu_probefilter_SoftQCDnonD_TuneCP5_13TeV-pythia8-evtgen_E313F59A-3E16-2342-B628-305F2AEEDCC5.root']
#options.inputFiles = ['/store/data/Run2018B/ParkingBPH4/MINIAOD/05May2019-v2/230000/6B5A24B1-0E6E-504B-8331-BD899EB60110.root'] if not options.isMC else \
# ['/store/cmst3/group/bpark/BToKmumu_1000Events_MINIAOD.root']
annotation = '%s nevts:%d' % (outputFileNANO, options.maxEvents)
from Configuration.StandardSequences.Eras import eras
process = cms.Process('BParkNANO', eras.Run2_2018)
# import of standard configurations
process.load('Configuration.StandardSequences.Services_cff')
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.EventContent.EventContent_cff')
import FWCore.ParameterSet.Config as cms
process = cms.Process("TEST")
runOnMC = True
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.parseArguments()
# Modify the defaults
if not options.inputFiles:
options.inputFiles = ['file:rawsim_numEvent100.root']
if options.outputFile == "output.root":
options.outputFile = "test_ntuple.root"
## 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 ),
#SkipEvent = cms.untracked.vstring('ProductNotFound')
)
## 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))
import FWCore.ParameterSet.Config as cms
from ROOT import TMath
###
# Defaults
###
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
options.inputFiles = "file:ME42TagAndProbeTree.root"
options.outputFile = "ME42TagAndProbeTreeAnalysis.root"
options.register ('analysisVar', 1, VarParsing.multiplicity.singleton, VarParsing.varType.int, "Select analyzer to run")
options.register ('inputDirectoryName', 'tagAndProbeTree', VarParsing.multiplicity.singleton, VarParsing.varType.string, "Analyzer input directory")
options.parseArguments()
###
# Analyzer
###
# Prepare Efficiency PSet
Eff = cms.PSet(
# pt
pt_ME42_tight = cms.PSet(
EfficiencyCategoryAndState = cms.vstring("isTightMuon","true"),
UnbinnedVariables = cms.vstring("mass"),
BinnedVariables = cms.PSet(
isME42 = cms.vdouble(0.5,1.5),
eta = cms.vdouble(1.2,1.8),
pt = cms.vdouble(20.0, 40.0, 60.0, 80.0, 100.0, 500.0),
),
BinToPDFmap = cms.vstring("voigtianPlusExponential")
options = VarParsing('analysis')
options.register(
"debug",
False,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"EDM File output"
)
options.parseArguments()
inputFiles = [
'/store/relval/CMSSW_8_1_0_pre9/RelValZMM_13/GEN-SIM-RECO/81X_mcRun2_asymptotic_v2_2023GReco-v1/10000/9041DC3E-5954-E611-8DB7-0025905B8566.root',
]
if len(options.inputFiles) is 0:
options.inputFiles = inputFiles
if len(options.outputFile) is 0:
options.outputFile = 'TnPTree.root'
process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring(options.inputFiles),
)
process.load('FWCore.MessageService.MessageLogger_cfi')
process.MessageLogger.cerr.FwkReport.reportEvery = 1000
process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(options.maxEvents) )
process.tracksSelection = cms.EDFilter("TrackRefSelector",
miniAOD_customizeMC(process)
## add command line parsing for cmsRun
## ussage is:
## cmsRun RecoTauTag/Configuration/test/edmTauVariables_slimmed_cfg.py inputFiles_load==RecoTauTag/Configuration/test/ZTT-validation.py \
## outputFile=MyOutputFile.root maxEvents=100
##
## get more inof here: https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideAboutPythonConfigFile#Passing_Command_Line_Arguments_T
from FWCore.ParameterSet.VarParsing import VarParsing
## define default values
from PhysicsTools.PatAlgos.patInputFiles_cff import filesRelValTTbarPileUpGENSIMRECO
process.source.fileNames = filesRelValTTbarPileUpGENSIMRECO
options = VarParsing('analysis')
options.inputFiles = filesRelValTTbarPileUpGENSIMRECO
options.outputFile = 'patMiniAOD_standard.root'
options.maxEvents = 100
## parse arguments
options.parseArguments()
options.inputFiles.pop(0)
process.source.fileNames = options.inputFiles
process.out.fileName = options.outputFile
process.maxEvents.input = options.maxEvents
## ------------------------------------------------------
# In addition you usually want to change the following
# parameters:
## ------------------------------------------------------
#
# process.GlobalTag.globaltag = ... ## (according to https://twiki.cern.ch/twiki/bin/view/CMS/SWGuideFrontierConditions)
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.maxEvents = 1000
options.parseArguments()
process = cms.Process("PATtest")
## MessageLogger
process.load("FWCore.MessageLogger.MessageLogger_cfi")
process.MessageLogger.cerr.FwkReport.reportEvery = 100
## 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))
## Geometry and Detector Conditions (needed for a few patTuple production steps)
process.load("Configuration.Geometry.GeometryRecoDB_cff")
process.load(
"Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
from Configuration.AlCa.GlobalTag import GlobalTag
import FWCore.ParameterSet.Config as cms
process = cms.Process("L1TMuonEndCap")
runOnMC = False
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.register('verbosity', 1, VarParsing.multiplicity.singleton,
VarParsing.varType.int, 'Verbosity')
options.parseArguments()
# Modify default arguments
if not options.inputFiles:
options.inputFiles = ['file:pippo.root']
if options.outputFile == "output.root":
options.outputFile = "output.root"
# Message logger
process.load("FWCore.MessageService.MessageLogger_cfi")
process.MessageLogger.cerr.FwkReport.reportEvery = 1
# Number of events
process.maxEvents = cms.untracked.PSet(
input=cms.untracked.int32(options.maxEvents))
# Input source
process.source = cms.Source(
"PoolSource",
fileNames=cms.untracked.vstring(options.inputFiles),
)
import FWCore.ParameterSet.Config as cms
###
# Defaults
###
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
options.inputFiles = "file:SingleMu_Run2012C_RECO.root"
options.outputFile = "SingleMu_Run2012C_Skim.root"
options.parseArguments()
###
# Includes
###
process = cms.Process("SingleMuSkim")
process.load("Configuration.StandardSequences.Services_cff")
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.load('FWCore.MessageService.MessageLogger_cfi')
###
# options
###
process.GlobalTag.globaltag = "FT_R_53_V18::All"
process.options = cms.untracked.PSet(
wantSummary = cms.untracked.bool(True),
# SkipEvent = cms.untracked.vstring('ProductNotFound')
)
process.MessageLogger.cerr.FwkReport.reportEvery = 1000
