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 enableCommandLineArguments(process):
options = VarParsing('analysis')
options.parseArguments()
# if options.outputFile=="":
# options.outputFile = "out_%s.root" % process.name_()
# print "output: %s" % options.outputFile
process.source.fileNames = cms.untracked.vstring(options.inputFiles)
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(options.maxEvents)
)
if hasattr(process, "out"):
process.out.fileName = cms.untracked.string(options.outputFile)
return (options.inputFiles, options.outputFile)
# test reco and dump into a tree
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.parseArguments()
from Configuration.StandardSequences.Eras import eras
process = cms.Process('AdvancedMultifit', eras.Run2_25ns)
# 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('Configuration.StandardSequences.GeometryRecoDB_cff')
#process.load('Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
#process.load('Configuration.StandardSequences.EndOfProcess_cff')
#process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
# 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('Configuration.StandardSequences.GeometryRecoDB_cff')
process.load(
'Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
process.load('Configuration.StandardSequences.RawToDigi_Data_cff')
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
# add a list of strings for events to process
options.parseArguments()
process = cms.Process("SimpleGenDumper")
process.load("FWCore.MessageService.MessageLogger_cfi")
process.MessageLogger.cerr.FwkReport.reportEvery = 10000
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10000) )
#process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
#process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(options.maxEvents) )
process.source = cms.Source("PoolSource",
# replace 'myfile.root' with the source file you want to use
#fileNames = cms.untracked.vstring (options.inputFiles),
#fileNames = cms.untracked.vstring('file:/tmp/amassiro/WW1Mevents_TUNE_NoUE_2.root')
fileNames = cms.untracked.vstring('root://cmsxrootd.fnal.gov//store/mc/RunIIWinter15GS/TT_TuneCUETP8M1_13TeV-powheg-pythia8/GEN-SIM/MCRUN2_71_V1-v1/10000/00004E38-D5C7-E411-8F03-0025905B8610.root',
'root://cmsxrootd.fnal.gov//store/mc/RunIIWinter15GS/TT_TuneCUETP8M1_13TeV-powheg-pythia8/GEN-SIM/MCRUN2_71_V1-v1/10000/0061B320-F9C7-E411-ADCC-0025904B2018.root',
'root://cmsxrootd.fnal.gov//store/mc/RunIIWinter15GS/TT_TuneCUETP8M1_13TeV-powheg-pythia8/GEN-SIM/MCRUN2_71_V1-v1/10000/00800509-D7C7-E411-BE5D-00266CFFA204.root',
'root://cmsxrootd.fnal.gov//store/mc/RunIIWinter15GS/TT_TuneCUETP8M1_13TeV-powheg-pythia8/GEN-SIM/MCRUN2_71_V1-v1/10000/00E79439-28C8-E411-BBC3-B083FED429D6.root')
)
process.TFileService = cms.Service("TFileService",
#fileName = cms.string("/tmp/amassiro/WW1Mevents_TUNE_NoUE_2_dump_tree.root"),
fileName = cms.string (options.outputFile),
import FWCore.ParameterSet.Config as cms
import os
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.outputFile = "L1TreeReco.root"
options.register('eActivityCut', 4, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Tau Veto HCAL Activity Threshold")
options.register('hActivityCut', 4, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Tau Veto ECAL Activity Threshold")
options.register('tauThresh', 7, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Tau Seed Threshold")
options.register('tauNThresh', 0, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Tau Neighbor Seed Threshold")
options.register('tauMaxPtTauVeto', 64, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Tau max Pt Tau Veto")
options.register('tauIsoLUTFile',
"L1Trigger/L1TCalorimeter/data/tauIsoLutTest0.12Iso.txt",
VarParsing.multiplicity.singleton, VarParsing.varType.string,
"Tau Isolation Cut")
options.register('tauCalibLUTFile',
"L1Trigger/L1TCalorimeter/data/tauCalibrationLUT_stage1.txt",
VarParsing.multiplicity.singleton, VarParsing.varType.string,
"Tau Calibration LUT")
options.register(
'RCTCalibrationMode',
'L1RCTParametersRcd_L1TDevelCollisions_ExtendedScaleFactorsV1',
VarParsing.multiplicity.singleton, VarParsing.varType.string,
"RCT Calibration Record")
#options.register ('lumilist', "Cert_246908_251883_13TeV_PromptReco_Collisions15_JSON_v2.txt", VarParsing.multiplicity.singleton, VarParsing.varType.string,
# "RCT Calibration Record")
#$Revision: 1.1 $
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('python')
options.register ('filepath',
False,
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"produces pdf weights for uncertainties study")
options.parseArguments()
filepath = options.filepath
process = cms.Process("TheNtupleMaker")
process.load("FWCore.MessageService.MessageLogger_cfi")
# See TheNtupleMaker twiki for a brief explanation
#process.MessageLogger.destinations = cms.untracked.vstring("cerr")
#process.MessageLogger.cerr.FwkReport.reportEvery = 10
#process.MessageLogger.cerr.default.limit = 5
# This is required in order to configure HLTConfigProducer
process.load("L1TriggerConfig.L1GtConfigProducers.L1GtConfig_cff")
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
process.source = cms.Source("PoolSource",
fileNames =
cms.untracked.vstring(
#'file:'+filepath,
#filepath,
#FIXME JSON file here
# New electron SF
# Check muon SF
# TO EXECUTE:
# cmsRun python/HHAnalyzer.py maxEvents=10
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
import os
import sys
options = VarParsing ('analysis')
options.register ('tCut', 0, options.multiplicity.singleton, options.varType.int,
"Trigger cut on/off")
options.register ('Data', False, options.multiplicity.singleton, options.varType.bool,
"Data on/off")
options.parseArguments()
#Un-comment to change output name - default is 'output' (useful for debug)
#options.outputFile = 'test'
# Determine sample name for MC stitching
sample = (options.inputFiles[0]).split('/')[-1].replace('.txt', '') if len(options.inputFiles) > 0 else ''
if sample=='list': sample = (options.inputFiles[0]).split('/')[-3]
process = cms.Process('ALPHA')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.MessageLogger.cerr.threshold = 'ERROR'
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(options.maxEvents) )
import FWCore.ParameterSet.Config as cms
import sys
# To execute test, run
# cmsRun boostedAnalysis_cfg.py isData=False outputFile=testrun maxEvents=100 inputFiles=file:/pnfs/desy.de/cms/tier2/store/user/hmildner/ttHTobb_M125_13TeV_powheg_pythia8/Boostedv2MiniAOD/151017_154254/0000/BoostedTTH_MiniAOD_1.root,file:/pnfs/desy.de/cms/tier2/store/user/hmildner/ttHTobb_M125_13TeV_powheg_pythia8/Boostedv2MiniAOD/151017_154254/0000/BoostedTTH_MiniAOD_2.root
# parse command-line arguments
# https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideCommandLineParsing
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
# The following variables are already defined in VarParsing class:
# 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
import glob, pickle
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
process = cms.Process("count")
process.load("FWCore.MessageService.MessageLogger_cfi")
options = VarParsing ('python')
options.register ('fileList',
'',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"files to cut and count over."
)
options.register ('CMSSW53x',
False,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"Running over a 53x sample"
)
options.register ('saveSourceVariable',
'',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"File name of the process.source variable to be saved in. If nothing given then it won't save."
)
import FWCore.ParameterSet.Config as cms
## removed cleaning from Exo VV package
##
process = cms.Process('NCUANA')
process.options = cms.untracked.PSet(wantSummary = cms.untracked.bool(True))
process.options = cms.untracked.PSet(
allowUnscheduled = cms.untracked.bool(True)
)
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
options.register ('runOnMC',
True,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"runOnMC")
options.register ('isReReco',
True,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"isReReco")
options.register ('useMiniAOD',
True,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"useMiniAOD")
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
import os
options = VarParsing('analysis')
options.register('skipEvents', 0, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Skip this number of events")
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'),
# Example config to produce a merged bad components payload, e.g.
# cmsRun makeMergeBadComponentPayload_example_cfg.py globalTag=auto:run3_data_prompt runNumber=319176 dqmFile=/eos/cms/store/group/comm_dqm/DQMGUI_data/Run2018/ZeroBias/R0003191xx/DQM_V0001_R000319176__ZeroBias__Run2018B-PromptReco-v2__DQMIO.root runStartTime=6574046031825076224
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing as VP
process = cms.Process("BadChannelMerge")
opts = VP("python")
opts.register("globalTag", "DONOTEXIST", VP.multiplicity.singleton,
VP.varType.string, "GlobalTag")
opts.register("dqmFile", "", VP.multiplicity.singleton, VP.varType.string,
"DQM root file")
opts.register("runNumber", 0, VP.multiplicity.singleton, VP.varType.int,
"run number")
opts.register("runStartTime", 0, VP.multiplicity.singleton, VP.varType.int,
"run start time")
opts.register("dbfile", "merged.db", VP.multiplicity.singleton,
VP.varType.string, "SQLite output file")
opts.register("outputTag", "SiStripBadComponents_merged_v0",
VP.multiplicity.singleton, VP.varType.string, "Output tag name")
opts.parseArguments()
notAllSet = False
if opts.globalTag == "DONOTEXIST":
print("ERROR: Global tag must be set")
notAllSet = True
if opts.runStartTime == 0:
print(
"ERROR: Run start time must be set (use the getRunStartTime.py script to get it)"
)
import FWCore.ParameterSet.Config as cms
### CMSSW command line parameter parser
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('python')
## data or MC options
options.register('isData', False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, 'flag to indicate data or MC')
options.register('maxEvts', -1, VarParsing.multiplicity.singleton,
VarParsing.varType.int,
'flag to indicate max events to process')
process = cms.Process("SEXAQ")
process.load("FWCore.MessageService.MessageLogger_cfi")
process.MessageLogger.cerr.FwkReport.reportEvery = cms.untracked.int32(100)
process.options = cms.untracked.PSet(wantSummary=cms.untracked.bool(True))
process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
process.load(
"Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
process.load('Configuration/StandardSequences/MagneticField_38T_cff')
process.load('Configuration/StandardSequences/Reconstruction_cff')
process.load('Configuration/EventContent/EventContent_cff')
from Configuration.AlCa.GlobalTag import GlobalTag
if (options.isData == True):
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
varOptions = VarParsing('analysis')
varOptions.register("isMC", False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "setup MC as in/out put")
varOptions.parseArguments()
process = cms.Process('GETGBR')
process.load('Configuration.StandardSequences.Services_cff')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
process.load(
'Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load(
'Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.MessageLogger.cerr.threshold = ''
process.MessageLogger.cerr.FwkReport.reportEvery = 1000
process.RandomNumberGeneratorService = cms.Service(
"RandomNumberGeneratorService",
calibratedPatElectrons=cms.PSet(
initialSeed=cms.untracked.uint32(1),
engineName=cms.untracked.string('TRandom3')),
calibratedElectrons=cms.PSet(initialSeed=cms.untracked.uint32(1),
engineName=cms.untracked.string('TRandom3')),
)
process.maxEvents = cms.untracked.PSet(
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
#input cmsRun options
options = VarParsing('analysis')
#with open('files_SUSYggH.txt') as f:
# options.inputFiles = f.readlines()
options.inputFiles = 'file:/data/uhussain/signalProduction/steps/CMSSW_8_0_3_patch2/src/Zprime_mchi_5GeV_4.root'
#options.inputFiles ='/store/mc/RunIISummer16MiniAODv2/ZJetsToNuNu_HT-200To400_13TeV-madgraph/MINIAODSIM/PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/120000/029F6BD2-5AC9-E611-BD77-0025907DC9D6.root'
#options.inputFiles = '/store/mc/RunIISpring16MiniAODv2/GluGluHToTauTau_M125_13TeV_powheg_pythia8/MINIAODSIM/PUSpring16RAWAODSIM_reHLT_80X_mcRun2_asymptotic_v14-v1/80000/065837E2-DA38-E611-9157-008CFA50291C.root'
options.outputFile = "Zprime_GenReco_dR02.root"
options.parseArguments()
#name the process
process = cms.Process("TreeProducerFromMiniAOD")
process.load('FWCore/MessageService/MessageLogger_cfi')
process.MessageLogger.cerr.FwkReport.reportEvery = 1
process.MessageLogger.cerr.threshold = cms.untracked.string('INFO')
#process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
#process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff")
#from Configuration.AlCa.GlobalTag_condDBv2 import GlobalTag
#process.GlobalTag = GlobalTag(process.GlobalTag, '80X_mcRun2_asymptotic_2016_TrancheIV_v8')
#process.load("Configuration.StandardSequences.MagneticField_cff")
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
# Auto generated configuration file
# using:
# Revision: 1.19
# Source: /local/reps/CMSSW/CMSSW/Configuration/Applications/python/ConfigBuilder.py,v
# with command line options: step2 --customise SLHCUpgradeSimulations/Configuration/postLS1Customs.customisePostLS1 --conditions auto:run2_mc -s DIGI:pdigi_valid,L1,DIGI2RAW,HLT:@frozen25ns,RAW2DIGI,L1Reco --datatier GEN-SIM-DIGI-RAW-HLTDEBUG -n 10 --magField 38T_PostLS1 --eventcontent FEVTDEBUGHLT --python step2.py --no_exec --filein file:step1.root --fileout file:step2.root
#
# CMSSW_7_4_1
#
import FWCore.ParameterSet.Config as cms
process = cms.Process('HLT')
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.register('mixHistFile', '', VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"histogram file saved for mixing pileups")
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.mix_2015_25ns_Startup_PoissonOOTPU_cfi')
process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
process.load('Configuration.StandardSequences.MagneticField_38T_PostLS1_cff')
process.load('Configuration.StandardSequences.Digi_cff')
process.load('Configuration.StandardSequences.SimL1Emulator_cff')
process.load('Configuration.StandardSequences.DigiToRaw_cff')
process.load('HLTrigger.Configuration.HLT_25ns14e33_v1_cff')
## process.l1GtTriggerMenu = cms.ESSource( "PoolDBESSource"
## , CondCore.DBCommon.CondDBCommon_cfi.CondDBCommon
## , toGet = cms.VPSet(
## cms.PSet(
## record = cms.string( 'L1GtTriggerMenuRcd' ),
## tag = cms.string( 'L1GtStableParameters_CRAFT09_hlt' ) # L1 menu for Fall10 REDIGI (CMSSW_3_8_7)
## )
## )
## )
## process.preferL1GtTriggerMenu = cms.ESPrefer( "PoolDBESSource", "l1GtTriggerMenu" )
import FWCore.PythonUtilities.LumiList as LumiList
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('python')
options.register ('useData',
0,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Run this on real data")
options.register ('doBinPtTrig',
0,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Bin the data in HLT paths by pt")
options.register ('semimuTriggers',
0,
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.maxEvents = -1
options.parseArguments()
muonID = ' userInt("highPtID") == 1'
process = cms.Process("Analyze")
process.load('Configuration.StandardSequences.Services_cff')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.EventContent.EventContent_cff')
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'
)
process.maxEvents = cms.untracked.PSet(
input=cms.untracked.int32(options.maxEvents))
process.source = cms.Source(
"PoolSource",
fileNames=cms.untracked.vstring(
"file:WR_3200_ToLNu_MASSNR_MiniAOD_MuMuJJ.root"),
)
process.options = cms.untracked.PSet(
SkipEvent=cms.untracked.vstring('ProductNotFound'))
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.",
)
# Auto generated configuration file
# using:
# Revision: 1.19
# Source: /local/reps/CMSSW/CMSSW/Configuration/Applications/python/ConfigBuilder.py,v
# with command line options: step3 --conditions 110X_mcRun4_realistic_v3 -n 10 --era Phase2C9 --eventcontent RECOSIM,DQM --runUnscheduled -s RAW2DIGI,RECO:reconstruction_trackingOnly,VALIDATION:@trackingOnlyValidation,DQM:@trackingOnlyDQM --datatier GEN-SIM-RECO,DQMIO --geometry Extended2026D49 --filein file:step2.root --fileout file:step3.root
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.register(
'nThreads', 1,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Number of multiprocessing units",
)
options.register(
'recoFile', 'file:reco.root',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Reco output file",
)
options.parseArguments()
from Configuration.Eras.Era_Phase2C9_cff import Phase2C9
process = cms.Process('RECO',Phase2C9)
# import of standard configurations
process.load('step3_tracking.modules.extras_cff')
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
import sys
process = cms.Process("tnp")
###################################################################
options = dict()
varOptions = VarParsing('analysis')
varOptions.register(
"isMC",
True,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"Compute MC efficiencies"
)
varOptions.parseArguments()
options['HLTProcessName'] = "HLT"
#options['PHOTON_COLL'] = "slimmedPhotons"
options['DIPHOTON_COLL'] = "flashggDiPhotons"
options['PHOTON_CUTS'] = "(abs(superCluster.eta)<2.5) && ((superCluster.energy*sin(superCluster.position.theta))>15.0)"
options['PHOTON_TAG_CUTS'] = "(abs(superCluster.eta)<=2.1) && !(1.4442<=abs(superCluster.eta)<=1.566) && (superCluster.energy*sin(superCluster.position.theta))>30.0"
options['MAXEVENTS'] = cms.untracked.int32(-1)
options['useAOD'] = cms.bool(False)
options['OUTPUTEDMFILENAME'] = 'edmFile.root'
options['DEBUG'] = cms.bool(False)
options['LEADING_PRESELECTION'] = """
(abs(leadingPhoton.superCluster.eta) < 2.5 && abs(subLeadingPhoton.superCluster.eta) < 2.5) &&
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('python')
options.register('runOnData', False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Run this on real data")
options.register('useRawLeptons', False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"Do not correct electrons/muons with smearing/energy scales")
options.register('era', 'era2017', VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"era to use (configurations defined in python/EraConfig.py)")
options.register('outFilename', 'MiniEvents.root',
VarParsing.multiplicity.singleton, VarParsing.varType.string,
"Output file name")
options.register('baseJetCollection', 'slimmedJets',
VarParsing.multiplicity.singleton, VarParsing.varType.string,
"Base jet collection")
options.register('inputFile', None, VarParsing.multiplicity.singleton,
VarParsing.varType.string, "input file to process")
options.register('lumiJson', None, VarParsing.multiplicity.singleton,
VarParsing.varType.string, 'apply this lumi json file')
options.register('saveTree', True, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "save summary tree")
options.register('savePF', False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, 'save PF candidates')
options.parseArguments()
#get the configuration to apply
from TopLJets2015.TopAnalysis.EraConfig import getEraConfiguration
globalTag, jecTag, jecDB = getEraConfiguration(era=options.era,
import FWCore.ParameterSet.Config as cms
import FWCore.Utilities.FileUtils as FileUtils
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
options.setDefault( 'maxEvents', -1)
options.setDefault( 'outputFile',
'ntuple.root' )
options.register( 'runOnMC',
False,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
'True for MC')
options.register( 'globalTag',
'START53_V7A',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
'CMS Global Tag')
options.register( 'wantSummary',
False,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"Print summary at end of job")
options.parseArguments()
import subprocess
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.register('eosdirs', '', VarParsing.multiplicity.list,
VarParsing.varType.string, "files location(s) on EOS")
options.register('localdirs', '', VarParsing.multiplicity.list,
VarParsing.varType.string,
"files location(s) on local filesystem")
options.register('datasets', '', VarParsing.multiplicity.list,
VarParsing.varType.string, "Input dataset(s)")
options.register('outname', 'Electron_iso.root',
VarParsing.multiplicity.singleton, VarParsing.varType.string,
"Output file name")
options.register('usegenpv', True, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Use MC-truth to define PV")
options.register('isTimingSample', False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"Process sample using timing and 4D vertexing")
options.register('saveTracksInfo', False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Save per-track information")
options.register('HGCToySim', False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"Allow to define eta and pt dependent time resolutions")
options.register('targetres', [0.03, 0.05, 0.07, 0.09, 0.15],
VarParsing.multiplicity.list, VarParsing.varType.float,
"Extra time resolution smearings")
options.register('dzCut', 0.1, VarParsing.multiplicity.singleton,
VarParsing.varType.float, "Trk-vtx dz cut")
options.register('debug', False, VarParsing.multiplicity.singleton,
process = cms.Process("EDMNtuple")
## MessageLogger
process.load("FWCore.MessageLogger.MessageLogger_cfi")
import sys
#input stuff for Run/Lumi selection with the "JSON"-formatted files from the PVT group
import FWCore.PythonUtilities.LumiList as LumiList
## Options and Output Report
process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
###############################
####### Parameters ############
###############################
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('python')
options.register('ignoreTrigger',
1,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Ignore trigger in selection")
options.register('ignoreDilepton',
0,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Ignore dilepton selection cut")
options.register ('runData',
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
import sys
process = cms.Process("tnpEGM")
###################################################################
## argument line options
###################################################################
varOptions = VarParsing('analysis')
varOptions.register("isMC", True, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Compute MC efficiencies")
varOptions.register("doEleID", False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Include tree for photon ID SF")
varOptions.register("doPhoID", False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Include tree for photon ID SF")
varOptions.register("doTrigger", False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Include tree for Trigger SF")
varOptions.register("doRECO", False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Include tree for Reco SF")
varOptions.register("calibEn", False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"use EGM smearer to calibrate photon and electron energy")
varOptions.register("isAOD", False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
import FWCore.ParameterSet.Config as cms
from FWCore.PythonUtilities.LumiList import LumiList
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing("analysis")
options.register("runList"
, []
, VarParsing.multiplicity.list
, VarParsing.varType.int
, "Run selection")
options.register("lumiList"
, "/afs/cern.ch/cms/CAF/CMSCOMM/COMM_DQM/certification/Collisions17/13TeV/DCSOnly/json_DCSONLY.txt"
, VarParsing.multiplicity.singleton
, VarParsing.varType.string
, "JSON file")
options.parseArguments()
lumilist = LumiList(filename = options.lumiList)
if len(options.runList) :
runlist = LumiList(runs = options.runList)
lumilist = lumilist & runlist
if not len(lumilist) :
raise RuntimeError, "The resulting LumiList is empty"
process = cms.Process("testRPCCPPF")
process.load("FWCore.MessageService.MessageLogger_cfi")
process.MessageLogger.cerr.FwkReport.reportEvery = 1000
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
# Starting with a skeleton process which gets imported with the following line
from PhysicsTools.PatAlgos.patTemplate_cfg import *
from PhysicsTools.PatAlgos.tools.coreTools import *
###############################
####### Parameters ############
###############################
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('python')
options.register ('useData',
False,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
'Run this on real data')
options.register ('globalTag',
'',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
'Overwrite defaul globalTag')
options.register ('forceCheckClosestZVertex',
False,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Force the check of the closest z vertex")
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")
#readFiles.extend( [
# 'file:/cmshome/gdimperi/GammaJet/JetCorrections/CMSSW_7_3_2/test/test_file_MINIAOD_for_JEC2015.root',
# ])
process.source = cms.Source (
"PoolSource",
fileNames = cms.untracked.vstring(
#'file:/cmshome/gdimperi/GammaJet/JetCorrections/CMSSW_7_3_2/test/test_file_MINIAOD_for_JEC2015.root'
'file:../tuples/GJET_Pythia/GJet_Pythia_80X_file1.root'
# 'file:/cmshome/fpreiato/GammaJet/CMSSW_7_4_14/src/JetMETCorrections/GammaJetFilter/analysis/tuples/QCD_MC/QCD_file1_ReReco.root'
)
)
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing()
options.register ('processedEvents',
'',
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"The number of processed events")
options.register ('crossSection',
'',
VarParsing.multiplicity.singleton,
VarParsing.varType.float,
"Dataset cross section")
options.register ('lowPtHat',
'',
## MessageLogger
process.load("FWCore.MessageLogger.MessageLogger_cfi")
#input stuff for Run/Lumi selection with the "JSON"-formatted files from the PVT group
import FWCore.PythonUtilities.LumiList as LumiList
## Options and Output Report
process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
###############################
####### Parameters ############
###############################
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('python')
options.register ('useFlavorHistory',
0,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Flavor History Mode")
options.register('sampleNameInput',
'Wjets',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Sample name to give histograms")
options.register ('allSys',
1,
Example:
cmsRun DiffractiveZPATTupleMulticrab.py Run=data_MuonA
Optional arguments:
Run = data_MuonP1, data_MuonP2, data_ElectronP1, data_ElectronP2, MC_none, MC_PU, MC_FlatWeight or MC_FlatWeight_and_PU
Authors: D. Figueiredo, R. Arciadiacono and N. Cartiglia
'''
import FWCore.ParameterSet.Config as cms
import os, sys
import atexit
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.register('Run', 'MC_none', VarParsing.multiplicity.singleton,
VarParsing.varType.string, "Option to Run: data or MC.")
options.parseArguments()
process = cms.Process("Analysis")
class config:
pass
config.verbose = True
config.writeEdmOutput = False
config.outputTTreeFile = 'DiffractiveZDataPATTuple.root'
config.runPATSequences = True
import math, sys
from os import environ
from os.path import exists, join, basename
def findFileInPath(theFile):
for s in environ["CMSSW_SEARCH_PATH"].split(":"):
attempt = join(s,theFile)
if exists(attempt):
return attempt
print "BABYMAKER: Could not find file "+theFile+". Not pre-applying JSON"
return None
###### Input parameters parsing
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
options.register('nEventsSample',
100,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Total number of events in dataset for event weight calculation.")
options.register('nEvents',
100,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Number of events to run over.")
options.register('json',
'babymaker/data/json/golden_Cert_271036-284044_13TeV_PromptReco_Collisions16_JSON_NoL1T.txt',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Path to json starting with babymaker/...")
./makeRateTrees_cfg.py
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
#'/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 = "$outputFileName"
options.register(
'eicIsolationThreshold',
3,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"EIC Isolation threshold")
options.register(
'hActivityCut',
0.5,
VarParsing.multiplicity.singleton,
from SUSYPhotonAnalysis.SusyNtuplizer.runOverAOD import configure
from FWCore.ParameterSet.VarParsing import VarParsing
# Call by e.g. "cmsRun runOverAOD_CommandLine.py dataset=53x22Jan2013 inputFiles=file1.root file2.root"
## get command line arguments
opt = VarParsing("analysis")
# 'analysis' contains some default options:
# maxEvents: Number of events to process (singleton integer)
# inputFiles: List of files to process as input (list string)
# secondaryFiles: List of secondary files (if needed; list string)
# outputFile: Name of output file (singleton string)
# secondaryOutput: Name of secondary output (if needed; singleton string)
# register an additional option
opt.register(
"dataset",
"", # default value
VarParsing.multiplicity.singleton,
VarParsing.varType.string, # string, int, or float
"Please choose a dataset token. For examples see config file."
)
opt.register(
"hltPaths",
[],
VarParsing.multiplicity.list,
VarParsing.varType.string,
"Please choose a hltPath token. For examples see config file."
)
opt.register(
"outputName",
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.EventContent.EventContent_cff')
process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
process.load(
'Configuration.StandardSequences.MagneticField_AutoFromDBCurrent_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load(
'Configuration.StandardSequences.FrontierConditions_GlobalTag_condDBv2_cff'
)
process.load("TrackingTools/TransientTrack/TransientTrackBuilder_cfi")
process.MessageLogger.cerr.FwkReport.reportEvery = 1000
# Handle different input options
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.register('targetPID', -1, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "set PID")
options.register('mcChannel', -1, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "set MC channel")
options.parseArguments()
inputFiles = []
if len(options.inputFiles) == 1 and (".root" not in options.inputFiles[0]):
flist = open(options.inputFiles[0])
inputFiles = flist.readlines()
flist.close()
inputFiles = [
'root://se01.grid.nchc.org.tw/{0}'.format(i) for i in inputFiles
]
else:
# SHyFT configuration
from PhysicsTools.PatAlgos.patTemplate_cfg import *
from PhysicsTools.PatAlgos.tools.coreTools import *
###############################
####### Parameters ############
###############################
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('python')
options.register ('useData',
0,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Run this on real data")
options.register ('hltProcess',
'HLT',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"HLT process name to use.")
options.register ('useMuon',
1,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Use muons (1) or electrons (0)")
import os, re
import FWCore.ParameterSet.Config as cms
### CMSSW command line parameter parser
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('python')
## matching cuts
options.register('dRmin',0.3,VarParsing.multiplicity.singleton,VarParsing.varType.float,'dR minimum cut');
options.register('pTmin',20.0,VarParsing.multiplicity.singleton,VarParsing.varType.float,'pT minimum cut');
## useVID
options.register('useGEDVID',False,VarParsing.multiplicity.singleton,VarParsing.varType.bool,'flag to use GED VID in cuts');
options.register('useOOTVID',False,VarParsing.multiplicity.singleton,VarParsing.varType.bool,'flag to use OOT VID in cuts');
## debug it all?
options.register('debug',False,VarParsing.multiplicity.singleton,VarParsing.varType.bool,'flag to dump info');
## data or MC options
options.register('isMC',True,VarParsing.multiplicity.singleton,VarParsing.varType.bool,'flag to indicate data or MC');
## GT to be used
options.register('globalTag','94X_mc2017_realistic_v14',VarParsing.multiplicity.singleton,VarParsing.varType.string,'gloabl tag to be used');
## do a demo run over only 1k events
options.register('demoMode',False,VarParsing.multiplicity.singleton,VarParsing.varType.bool,'flag to run over only 1k events');
## processName
options.register('processName','TREE',VarParsing.multiplicity.singleton,VarParsing.varType.string,'process name to be considered');
## outputFile Name
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
import sys
options = VarParsing('analysis')
options.register(
"isMC",
True,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"Compute MC efficiencies"
)
options.register(
"inputFileName",
"TnPTree_mc.root",
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Input filename"
)
options.register(
"outputFileName",
"test",
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Output filename"
)
options.register(
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
from PhysicsTools.PatAlgos.tools.helpers import getPatAlgosToolsTask
import sys
import os
## Define the process
process = cms.Process("Filter")
patAlgosToolsTask = getPatAlgosToolsTask(process)
process.options = cms.untracked.PSet(wantSummary=cms.untracked.bool(True))
## Set up command line options
options = VarParsing('analysis')
options.register('runOnGenOrAODsim', True, 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 = [
gROOT.ProcessLine(\
"struct MyStruct{\
UInt_t runnr;\
UInt_t luminosityBlock;\
Double_t intgRecLumi;\
Double_t HLTpresc;\
Double_t L1presc;\
};")
from ROOT import MyStruct
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('python')
options.register ('outfilename',
'OutFile.root',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
'name of output root file')
options.register ('infilename',
'TempOut.csv',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
'name of input file from lcr2.py')
options.parseArguments()
./makeRateTrees_cfg.py
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/ETauSkim/skim_12_1_erV.root'
#options.inputFiles = '/store/user/tapas/2012-08-01-CRAB_ZEESkim/skim_10_1_wd2.root'
options.outputFile = "uct_rate_tree_mc.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,
import FWCore.ParameterSet.Config as cms
process = cms.Process('HiForest')
process.options = cms.untracked.PSet(
)
#parse command line arguments
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.register ('isPP',
False,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"Flag if this is a pp simulation")
options.parseArguments()
#####################################################################################
# HiForest labelling info
#####################################################################################
#####################################################################################
# Input source
#####################################################################################
process.source = cms.Source("PoolSource",
duplicateCheckMode = cms.untracked.string("noDuplicateCheck"),
fileNames = cms.untracked.vstring(options.inputFiles[0])
)
# Number of events we want to process, -1 = all events
def SingleTopStep1(process, ):
options = VarParsing('analysis')
options.register('isMC', True, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Run on MC")
options.register('doDebug', False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Run in debugging mode")
options.register('doSkimming', True, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Preselect events")
options.register('doSlimming', True, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Drop unnecessary collections")
options.register('doMuon', True, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Do muon paths")
options.register('doElectron', True, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Do electron paths")
#Tag from https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideFrontierConditions?redirectedfrom=CMS.SWGuideFrontierConditions#2012_MC_production
# Latest for "53Y Releases (MC)"
options.register('globalTag', Config.globalTagMC,
VarParsing.multiplicity.singleton,
VarParsing.varType.string, "Global tag")
options.parseArguments()
process.source.fileNames = cms.untracked.vstring(options.inputFiles)
process.maxEvents = cms.untracked.PSet(
input=cms.untracked.int32(options.maxEvents))
process.out.fileName = cms.untracked.string(options.outputFile)
process.options = cms.untracked.PSet(
wantSummary=cms.untracked.bool(options.doDebug))
if options.doDebug:
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")
postfix = ""
jetCorr = ['L1FastJet', 'L2Relative', 'L3Absolute']
if not options.isMC:
jetCorr += ['L2L3Residual']
usePF2PAT(
process,
runPF2PAT=True,
jetAlgo='AK5',
runOnMC=options.isMC,
postfix=postfix,
jetCorrections=('AK5PFchs', jetCorr),
pvCollection=cms.InputTag('goodOfflinePrimaryVertices'),
#typeIMetCorrections = True
typeIMetCorrections=
False #Type1 MET now applied later using runMETUncertainties
)
# https://twiki.cern.ch/twiki/bin/view/CMSPublic/WorkBookJetEnergyCorrections#JetEnCorPFnoPU2012
process.pfPileUp.Enable = True
process.pfPileUp.checkClosestZVertex = False
#-------------------------------------------------
# selection step 2: vertex filter
#-------------------------------------------------
# https://twiki.cern.ch/twiki/bin/view/CMSPublic/WorkBookJetEnergyCorrections#JetEnCorPFnoPU2012
# https://twiki.cern.ch/twiki/bin/viewauth/CMS/TWikiTopRefEventSel#Cleaning_Filters
process.goodOfflinePrimaryVertices = cms.EDFilter(
"PrimaryVertexObjectFilter",
filterParams=cms.PSet(minNdof=cms.double(4.0),
maxZ=cms.double(24.0),
maxRho=cms.double(2.0)),
filter=cms.bool(True),
src=cms.InputTag('offlinePrimaryVertices'))
#-------------------------------------------------
# Muons
#-------------------------------------------------
#if not maxLeptonIso is None:
# process.pfIsolatedMuons.isolationCut = maxLeptonIso
#Use both isolated and non-isolated muons as a patMuon source
#process.patMuons.pfMuonSource = cms.InputTag("pfMuons")
#process.muonMatch.src = cms.InputTag("pfMuons")
process.patMuons.pfMuonSource = cms.InputTag("pfIsolatedMuons")
process.muonMatch.src = cms.InputTag("pfIsolatedMuons")
process.selectedPatMuons.cut = "pt>20 && abs(eta)<3.0"
# muon ID production (essentially track count embedding) must be here
# because tracks get dropped from the collection after this step, resulting
# in null ptrs.
process.muonsWithID = cms.EDProducer(
'MuonIDProducer',
muonSrc=cms.InputTag("selectedPatMuons"),
primaryVertexSource=cms.InputTag("goodOfflinePrimaryVertices"))
#process.muonClones = cms.EDProducer("MuonShallowCloneProducer",
# src = cms.InputTag("selectedPatMuons")
#)
#-------------------------------------------------
# Electrons
# Implemented as in https://indico.cern.ch/getFile.py/access?contribId=1&resId=0&materialId=slides&confId=208765
#-------------------------------------------------
#if not maxLeptonIso is None:
# process.pfIsolatedElectrons.isolationCut = maxLeptonIso
#Use both isolated and un-isolated electrons as patElectrons.
#NB: no need to change process.electronMatch.src to pfElectrons,
# it's already gsfElectrons, which is a superset of the pfElectrons
process.patElectrons.pfElectronSource = cms.InputTag("pfElectrons")
process.load('EGamma.EGammaAnalysisTools.electronIdMVAProducer_cfi')
process.mvaID = cms.Sequence(process.mvaTrigV0 + process.mvaNonTrigV0)
process.patElectrons.electronIDSources.mvaTrigV0 = cms.InputTag(
"mvaTrigV0")
process.patElectrons.electronIDSources.mvaNonTrigV0 = cms.InputTag(
"mvaNonTrigV0")
process.patPF2PATSequence.replace(process.patElectrons,
process.mvaID * process.patElectrons)
process.selectedPatElectrons.cut = "pt>25 && abs(eta)<3.0"
process.electronsWithID = cms.EDProducer(
'ElectronIDProducer',
electronSrc=cms.InputTag("selectedPatElectrons"),
primaryVertexSource=cms.InputTag("goodOfflinePrimaryVertices"))
#process.electronClones = cms.EDProducer("ElectronShallowCloneProducer",
# src = cms.InputTag("selectedPatElectrons")
#)
#if not maxLeptonIso is None:
# process.pfIsolatedElectrons.isolationCut = maxLeptonIso
#electron dR=0.3
process.pfElectrons.isolationValueMapsCharged = cms.VInputTag(
cms.InputTag("elPFIsoValueCharged03PFId"))
process.pfElectrons.deltaBetaIsolationValueMap = cms.InputTag(
"elPFIsoValuePU03PFId")
process.pfElectrons.isolationValueMapsNeutral = cms.VInputTag(
cms.InputTag("elPFIsoValueNeutral03PFId"),
cms.InputTag("elPFIsoValueGamma03PFId"))
process.pfElectrons.deltaBetaIsolationValueMap = cms.InputTag(
"elPFIsoValuePU03PFId")
process.pfElectrons.isolationValueMapsNeutral = cms.VInputTag(
cms.InputTag("elPFIsoValueNeutral03PFId"),
cms.InputTag("elPFIsoValueGamma03PFId"))
process.patElectrons.isolationValues.pfNeutralHadrons = cms.InputTag(
"elPFIsoValueNeutral03PFId")
process.patElectrons.isolationValues.pfChargedAll = cms.InputTag(
"elPFIsoValueChargedAll03PFId")
process.patElectrons.isolationValues.pfPUChargedHadrons = cms.InputTag(
"elPFIsoValuePU03PFId")
process.patElectrons.isolationValues.pfPhotons = cms.InputTag(
"elPFIsoValueGamma03PFId")
process.patElectrons.isolationValues.pfChargedHadrons = cms.InputTag(
"elPFIsoValueCharged03PFId")
process.pfIsolatedElectrons.isolationValueMapsCharged = cms.VInputTag(
cms.InputTag("elPFIsoValueCharged03PFId"))
process.pfIsolatedElectrons.deltaBetaIsolationValueMap = cms.InputTag(
"elPFIsoValuePU03PFId")
process.pfIsolatedElectrons.isolationValueMapsNeutral = cms.VInputTag(
cms.InputTag("elPFIsoValueNeutral03PFId"),
cms.InputTag("elPFIsoValueGamma03PFId"))
#-------------------------------------------------
# Jets
# MET corrections as https://twiki.cern.ch/twiki/bin/view/CMSPublic/WorkBookMetAnalysis#Type_I_0_with_PAT
#-------------------------------------------------
#pfNoTau == True => remove taus from jets
#process.pfNoTau.enable = noTau
process.selectedPatJets.cut = cms.string("pt>30 && abs(eta)<5.0")
process.load("CMGTools.External.pujetidsequence_cff")
process.patPF2PATSequence += process.puJetIdSqeuence
#process.jetClones = cms.EDProducer("CaloJetShallowCloneProducer",
# src = cms.InputTag("seletedPatJets")
#)
#-------------------------------------------------
# MET uncertainty step
#-------------------------------------------------
#Embed the reference to the original jet in the jets, which is constant during the propagation
process.patJetsWithOwnRef = cms.EDProducer(
'PatObjectOwnRefProducer<pat::Jet>',
src=cms.InputTag("selectedPatJets"))
#Note: this module causes a large memory increase when crossing the file boundary
#Reason - unknown, solution: limit processing to ~1 file.
from PhysicsTools.PatUtils.tools.metUncertaintyTools import runMEtUncertainties
runMEtUncertainties(
process,
electronCollection=cms.InputTag("electronsWithID"),
photonCollection=None,
muonCollection=cms.InputTag("muonsWithID"),
tauCollection="", # "" means emtpy, None means cleanPatTaus
jetCollection=cms.InputTag("patJetsWithOwnRef"),
jetCorrLabel="L3Absolute" if options.isMC else "L2L3Residual",
doSmearJets=options.isMC,
jetCorrPayloadName="AK5PFchs",
addToPatDefaultSequence=False)
process.stpolMetUncertaintySequence = cms.Sequence(
process.metUncertaintySequence)
if not options.doSlimming:
process.out.outputCommands = cms.untracked.vstring('keep *')
else:
process.out.outputCommands = cms.untracked.vstring([
'drop *',
'keep edmMergeableCounter_*_*_*', # Keep the lumi-block counter information
'keep edmTriggerResults_TriggerResults__*', #Keep the trigger results
'keep *_genParticles__*', #keep all the genParticles
#'keep recoVertexs_offlinePrimaryVertices__*', #keep the offline PV-s
'keep recoVertexs_goodOfflinePrimaryVertices__*', #keep the offline PV-s
# Jets
'keep patJets_*__*',
'keep double_*_rho_*', #For rho-corr rel iso
'keep recoGenJets_selectedPatJets_genJets_*', #For Jet MC smearing we need to keep the genJets
"keep *_puJetId_*_*", # input variables
"keep *_puJetMva_*_*", # final MVAs and working point flags
'keep *_jetClones__*',
# Muons
'keep patMuons_*__*',
'keep *_muonClones__*',
# Electrons
'keep patElectrons_*__*',
'keep *_electronClones__*',
# METs
'keep patMETs_*__*',
#ECAL laser corr filter
'keep bool_ecalLaserCorrFilter__*',
#For flavour analyzer
'keep GenEventInfoProduct_generator__*',
#PU info
'keep PileupSummaryInfos_addPileupInfo__*',
##PFCandidates
#'keep recoPFCandidates_*_pfCandidates_PAT',
#'keep recoPFMETs_pfMET__*',
#'keep recoPFMETs_pfMet__*',
#'keep recoGenMETs_genMetTrue__*',
#'keep recoPFCandidates_particleFlow__*',
#'keep recoConversions_allConversions__*',
#'keep recoVertexCompositeCandidates_generalV0Candidates_*_*',
#'keep recoTracks_generalTracks__*',
#'keep recoBeamSpot_offlineBeamSpot__*',
#'keep recoMuons_muons__*',
'keep int_*__PAT',
'keep ints_*__PAT',
'keep double_*__PAT',
'keep doubles_*__PAT',
'keep float_*__PAT',
'keep floats_*__PAT',
])
#FIXME: is this correct?
#Keep events that pass either the muon OR the electron path
process.out.SelectEvents = cms.untracked.PSet(SelectEvents=cms.vstring([]))
#-------------------------------------------------
# Paths
#-------------------------------------------------
process.goodOfflinePVCount = cms.EDProducer(
"CollectionSizeProducer<reco::Vertex>",
src=cms.InputTag("goodOfflinePrimaryVertices"))
process.preCalcSequences = cms.Sequence(process.patJetsWithOwnRef *
process.goodOfflinePVCount)
process.patPF2PATSequence.insert(
process.patPF2PATSequence.index(process.selectedPatMuons) + 1,
process.muonsWithID)
process.patPF2PATSequence.insert(
process.patPF2PATSequence.index(process.selectedPatElectrons) + 1,
process.electronsWithID)
#Need separate paths because of skimming
if options.doMuon:
process.singleTopPathStep1Mu = cms.Path(
process.goodOfflinePrimaryVertices * process.patPF2PATSequence
#* process.muonClones
#* process.electronClones
#* process.jetClones
)
if options.doElectron:
process.singleTopPathStep1Ele = cms.Path(
process.goodOfflinePrimaryVertices * process.patPF2PATSequence
#* process.muonClones
#* process.electronClones
#* process.jetClones
)
if options.doMuon:
process.out.SelectEvents.SelectEvents.append("singleTopPathStep1Mu")
if options.doElectron:
process.out.SelectEvents.SelectEvents.append("singleTopPathStep1Ele")
process.GlobalTag.globaltag = cms.string(options.globalTag)
process.singleTopPathStep1Mu += process.preCalcSequences
process.singleTopPathStep1Ele += process.preCalcSequences
if options.doMuon:
process.singleTopPathStep1Mu += process.stpolMetUncertaintySequence
if options.doElectron:
process.singleTopPathStep1Ele += process.stpolMetUncertaintySequence
if options.isMC:
#https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideBTagJetProbabilityCalibration?redirectedfrom=CMS.SWGuideBTagJetProbabilityCalibration#Calibration_in_53x_Data_and_MC
process.GlobalTag.toGet = cms.VPSet(
cms.PSet(record=cms.string("BTagTrackProbability2DRcd"),
tag=cms.string("TrackProbabilityCalibration_2D_MC53X_v2"),
connect=cms.untracked.string(
"frontier://FrontierPrep/CMS_COND_BTAU")),
cms.PSet(record=cms.string("BTagTrackProbability3DRcd"),
tag=cms.string("TrackProbabilityCalibration_3D_MC53X_v2"),
connect=cms.untracked.string(
"frontier://FrontierPrep/CMS_COND_BTAU")))
else:
#https://twiki.cern.ch/twiki/bin/view/CMSPublic/SWGuideBTagJetProbabilityCalibration?redirectedfrom=CMS.SWGuideBTagJetProbabilityCalibration#Calibration_in_53x_Data_and_MC
process.GlobalTag.toGet = cms.VPSet(
cms.PSet(
record=cms.string("BTagTrackProbability2DRcd"),
tag=cms.string("TrackProbabilityCalibration_2D_Data53X_v2"),
connect=cms.untracked.string(
"frontier://FrontierPrep/CMS_COND_BTAU")),
cms.PSet(
record=cms.string("BTagTrackProbability3DRcd"),
tag=cms.string("TrackProbabilityCalibration_3D_Data53X_v2"),
connect=cms.untracked.string(
"frontier://FrontierPrep/CMS_COND_BTAU")))
process.load('RecoMET.METFilters.ecalLaserCorrFilter_cfi')
process.ecalLaserCorrFilter.taggingMode = True
# https://twiki.cern.ch/twiki/bin/viewauth/CMS/TWikiTopRefEventSel#Cleaning_Filters
process.scrapingFilter = cms.EDFilter(
"FilterOutScraping",
applyfilter=cms.untracked.bool(True),
debugOn=cms.untracked.bool(False),
numtrack=cms.untracked.uint32(10),
thresh=cms.untracked.double(0.25))
#if doElectron:
# process.singleTopPathStep1Ele.insert(0, process.scrapingFilter)
#if doMuon:
# process.singleTopPathStep1Mu.insert(0, process.scrapingFilter)
process.patPF2PATSequence += process.scrapingFilter
process.patPF2PATSequence += process.ecalLaserCorrFilter
#if not onGrid:
# from SingleTopPolarization.Analysis.cmdlineParsing import enableCommandLineArguments
# enableCommandLineArguments(process)
#else:
# process.out.fileName = "step1.root"
if options.doSkimming:
process.out.fileName.setValue(process.out.fileName.value().replace(
".root", "_Skim.root"))
else:
process.out.fileName.setValue(process.out.fileName.value().replace(
".root", "_noSkim.root"))
#-----------------------------------------------
# Skimming
#-----------------------------------------------
#Throw away events before particle flow?
if options.doSkimming:
from SingleTopPolarization.Analysis.eventSkimming_cfg import skimFilters
skimFilters(process)
if options.doMuon:
process.singleTopPathStep1Mu.insert(0, process.muonSkim)
if options.doElectron:
process.singleTopPathStep1Ele.insert(0, process.electronSkim)
#-----------------------------------------------
# Skim efficiency counters
#-----------------------------------------------
#count all processed events
countProcessed(process)
#count events passing mu and ele paths
if options.doMuon:
countInSequence(process, process.singleTopPathStep1Mu)
if options.doElectron:
countInSequence(process, process.singleTopPathStep1Ele)
#-------------------------------------------------
#
#-------------------------------------------------
if not options.doSlimming:
process.out.fileName.setValue(process.out.fileName.value().replace(
".root", "_noSlim.root"))
return process
#from make_file_list import *
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
process = cms.Process("STPOLPDF")
process.load("FWCore.MessageService.MessageLogger_cfi")
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
#filelist = get_file_list()
options = VarParsing('analysis')
options.register ('PowhegTopMassFix', False,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"Apply Powheg scale fix?")
options.parseArguments()
#print "option", options
process.source = cms.Source("PoolSource",
# replace 'myfile.root' with the source file you want to use
fileNames = cms.untracked.vstring(
#filelist
#'file:/hdfs/cms/store/user/joosep/T_t-channel_TuneZ2star_8TeV-powheg-tauola/Jul8_51f69b/e6cee8f160cb239f414a74aa40d88ea9/output_noSkim_1_1_Nq6.root'
)
)
process.pdfWeights = cms.EDProducer('PDFWeightProducer',
PDFSets = cms.vstring([
e.g.: cmsRun runP_Tracker_cfg.py geom="XYZ"
The default component to be monitored is the Tracker. If other
components need to be studied, they must be supplied, one at a time,
at the command line, e.g.: cmsRun runP_Tracker_cfg.py
label="XYZ"
"""
from Validation.Geometry.plot_utils import _LABELS2COMPS
process.load("FWCore.MessageLogger.MessageLogger_cfi")
_ALLOWED_LABELS = _LABELS2COMPS.keys()
options = VarParsing('analysis')
options.register('geom', #name
'Extended2017Plan1', #default value
VarParsing.multiplicity.singleton, # kind of options
VarParsing.varType.string, # type of option
"Select the geometry to be studied" # help message
)
options.register('label', #name
'Tracker', #default value
VarParsing.multiplicity.singleton, # kind of options
VarParsing.varType.string, # type of option
"Select the label to be used to create output files. Default to tracker. If multiple components are selected, it defaults to the join of all components, with '_' as separator." # help message
)
options.setDefault('inputFiles', ['file:single_neutrino_random.root'])
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(
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
# add a list of strings for events to process
options.register ('isMiniAod',
False,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"is miniAod? (default = False). It changes the collection names")
options.register ('debug',
False,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"debug? (default = False).")
options.register ('mcLHERunInfoTag',
'externalLHEProducer',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
'LHE run information')
options.register ('mcLHEEventInfoTag',
'externalLHEProducer',
#'generator',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
'LHE event information')
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
# input cmsRun options
options = VarParsing("analysis")
# with open('files_SUSYggH.txt') as f:
# options.inputFiles = f.readlines()
options.outputFile = "MiniAODeffi_dmf_SUSYggH.root"
options.parseArguments()
# name the process
process = cms.Process("TreeProducerFromMiniAOD")
process.load("FWCore/MessageService/MessageLogger_cfi")
process.MessageLogger.cerr.FwkReport.reportEvery = 5000
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(-1))
process.source = cms.Source("PoolSource", fileNames=cms.untracked.vstring(options.inputFiles))
# 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 dbs:/VBFHToMuMu_M-125_TuneCP5_PSweights_13TeV_powheg_pythia8/RunIISummer16DR80Premix-PUMoriond17_80X_mcRun2_asymptotic_2016_TrancheIV_v6-v1/AODSIM --fileout file:HIG-RunIISummer16MiniAODv3-03813.root --mc --eventcontent MINIAODSIM --runUnscheduled --datatier MINIAODSIM --conditions 94X_mcRun2_asymptotic_v3 --step PAT --nThreads 8 --era Run2_2016,run2_miniAOD_80XLegacy --python_filename HIG-RunIISummer16MiniAODv3-03813_1_cfg.py --no_exec --customise Configuration/DataProcessing/Utils.addMonitoring -n 10000
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing('analysis')
options.register('nEvents', 0, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "nEvents")
options.parseArguments()
from Configuration.StandardSequences.Eras import eras
process = cms.Process('PAT', eras.Run2_2016)
# 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_cff')
process.load('PhysicsTools.PatAlgos.slimming.metFilterPaths_cff')
process.load('Configuration.StandardSequences.PATMC_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load(
'Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
import FWCore.ParameterSet.Config as cms
from DQMServices.Core.DQM_cfg import *
from DQM.L1TMonitor.environment_file_cfi import *
#RunType, and Runkey selection from RCMS
import sys
from FWCore.ParameterSet.VarParsing import VarParsing
#from DQM.Integration.test.dqmPythonTypes import *
runParameters = VarParsing ('analysis')
runParameters.register ('runtype',
'pp_run',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Type of Run in CMS")
runParameters.register ('runkey',
'pp_run',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Run Keys of CMS")
# Fix to allow scram to compile
if len(sys.argv) > 1:
runParameters.parseArguments()
#runType = RunType(['pp_run','cosmic_run','hi_run','hpu_run'])
if not runParameters.runkey.strip():
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.mix_POISSON_average_cfi')
process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
process.load('Configuration.StandardSequences.MagneticField_cff')
process.load('Configuration.StandardSequences.RawToDigi_cff')
process.load('Configuration.StandardSequences.Reconstruction_cff')
process.load('Configuration.StandardSequences.Validation_cff')
process.load('DQMOffline.Configuration.DQMOfflineMC_cff')
process.load(
'Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.load('CNNFiltering.CNNAnalyze.CNNAnalyzePixel_cfi')
options = VarParsing('analysis')
options.register('pileUp', 50, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Pileup")
options.register('skipEvent', 0, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Skip Events")
options.register('numEvents', 100, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "Max Events")
options.register('numFile', 2, VarParsing.multiplicity.singleton,
VarParsing.varType.int, "File Number")
options.register('openDataVM', True, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "On OpenDataVM?")
options.parseArguments()
opendata = [
"file:/lustre/cms/store/user/adiflori/OpenData/opendatas/017E1C89-5F83-D94E-8FE2-93CBE7369C11.root",
"file:/lustre/cms/store/user/adiflori/OpenData/opendatas/18987864-CE50-9341-88A2-47FD4ACD1AE9.root",
process.MessageLogger.cerr.FwkReport.reportEvery = 5000
process.load("RecoTauTag.Configuration.RecoPFTauTag_cff")
from PhysicsTools.PatAlgos.tools.coreTools import *
process.load('Configuration.StandardSequences.Services_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
process.load('JetMETCorrections.Configuration.DefaultJEC_cff')
#--------------------------------------------------------------------------------
################### Take inputs from crab.cfg file ##############
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('python')
options.register ('isMC',
True,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"Sample Type: MC or data")
options.register ('channel',
'ett',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Desired channel: mtt, ett or none")
options.register ('globaltag',
'START53_V21::All',
VarParsing.multiplicity.singleton,
## MessageLogger
process.load("FWCore.MessageLogger.MessageLogger_cfi")
#input stuff for Run/Lumi selection with the "JSON"-formatted files from the PVT group
import FWCore.PythonUtilities.LumiList as LumiList
## Options and Output Report
process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
###############################
####### Parameters ############
###############################
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('python')
options.register ('useFlavorHistory',
0,
VarParsing.multiplicity.singleton,
VarParsing.varType.int,
"Flavor History Mode")
options.register('sampleNameInput',
'Wjets',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Sample name to give histograms")
options.register ('qcdMC',
0,
###############################
####### Parameters ############
###############################
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('python')
options.register('runOnData',
False,
VarParsing.multiplicity.singleton,
VarParsing.varType.bool,
"Run this on real data"
)
options.register('globalTag',
'START42_V14B::All',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Global tag to be used"
)
options.register('outputFilename',
'MyNtupleMaker_output.root',
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Name of the output file"
)
options.register('reportEvery',
10,
VarParsing.multiplicity.singleton,
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
import sys
options = VarParsing('analysis')
options.register("isMC", False, VarParsing.multiplicity.singleton,
VarParsing.varType.bool, "Compute MC efficiencies")
options.register(
"inputFileName",
#"/afs/cern.ch/work/i/ishvetso/public/for_Matteo/TnPTree_mc-powheg.root",
"/data2/pwang/TnP/76X_v2/TnP_Data.root",
#"TnPTree_mc.root",
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"Input filename")
options.register("outputFileName", "Signal", VarParsing.multiplicity.singleton,
VarParsing.varType.string, "Output filename")
options.register(
"idName",
"passingLoose",
#"passingTrigWP90",
VarParsing.multiplicity.singleton,
VarParsing.varType.string,
"ID variable name as in the fitter_tree")
options.register("dirName", "PhotonToRECO", VarParsing.multiplicity.singleton,
VarParsing.varType.string,
# Auto generated configuration file
# using:
# Revision: 1.19
# Source: /local/reps/CMSSW/CMSSW/Configuration/Applications/python/ConfigBuilder.py,v
# with command line options: Configuration/GenProduction/python/HIG-RunIIFall17wmLHEGS-00389-fragment.py --fileout file:step1.root --mc --eventcontent RAWSIM,LHE --datatier GEN-SIM,LHE --conditions 93X_mc2017_realistic_v3 --beamspot Realistic25ns13TeVEarly2017Collision --step LHE,GEN,SIM --nThreads 8 --geometry DB:Extended --era Run2_2017 --python_filename step1.py --no_exec --customise Configuration/DataProcessing/Utils.addMonitoring -n 500
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
options = VarParsing ('analysis')
options.register('jobNum', 0, VarParsing.multiplicity.singleton,VarParsing.varType.int,"jobNum")
options.parseArguments()
from Configuration.StandardSequences.Eras import eras
process = cms.Process('SIM',eras.Run2_2017)
# 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.GeometrySimDB_cff')
process.load('Configuration.StandardSequences.MagneticField_cff')
process.load('Configuration.StandardSequences.Generator_cff')
process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic25ns13TeVEarly2017Collision_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')
import sys
import time
from ROOT import *
from utils import *
from utilsII import *
from glob import glob
from random import shuffle
from FWCore.ParameterSet.VarParsing import VarParsing
import random
datamc = 'mc'
mZ = 91
doGenVersion = True
RelaxGenKin = False
options = VarParsing('python')
options.parseArguments()
weight = 1 #################################WEIGHT
gROOT.SetBatch()
gROOT.SetStyle('Plain')
inputFiles = options.inputFiles
if inputFiles == []:
print 'running on small default DYtoLL sample'
inputFiles = [
"/pnfs/desy.de/cms/tier2/store/user/sbein/NtupleHub/Production2016v1/Summer16.WJetsToLNu_HT-2500ToInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8_ext1_4_RA2AnalysisTree.root"
]
x = len(inputFiles)
c = TChain('TreeMaker2/PreSelection')
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')
)
#FIXME JSON file here
# New electron SF
# Check muon SF
# TO EXECUTE:
# cmsRun python/HHAnalyzer.py maxEvents=10
import FWCore.ParameterSet.Config as cms
from FWCore.ParameterSet.VarParsing import VarParsing
import os
import sys
options = VarParsing('analysis')
options.register('tCut', 0, options.multiplicity.singleton,
options.varType.int, "Trigger cut on/off")
options.parseArguments()
#Un-comment to change output name - default is 'output' (useful for debug)
#options.outputFile = 'test'
# Determine sample name for MC stitching
sample = (options.inputFiles[0]).split('/')[-1].replace(
'.txt', '') if len(options.inputFiles) > 0 else ''
if sample == 'list': sample = (options.inputFiles[0]).split('/')[-3]
process = cms.Process('ALPHA')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.MessageLogger.cerr.threshold = 'ERROR'
process.maxEvents = cms.untracked.PSet(
input=cms.untracked.int32(options.maxEvents))
