isMC = cms.bool(True),
#AllEvents = copy.deepcopy(allEvents),
# allMuons = copy.deepcopy(allMuons),
#selectedMuons = copy.deepcopy(selectedMuons),
#selectedDimuons = copy.deepcopy(selectedDimuons),
#selectedPhotons = copy.deepcopy(selectedPhotons),
#selectedZgs = copy.deepcopy(selectedZgs),
FullSelection = copy.deepcopy(FullSelection),
)
## Histograms configuration
process.histograms = copy.deepcopy(histograms)
#______________________________________________________________________________
## Event weight
process.eventWeight = cms.PSet(
crossSectionInPicoBarns = cms.double(45.2),
scaleToLumiInInverseFemtoBarns = cms.double(4.9989),
totalProcessedEvents = cms.int64(472380),
)
#______________________________________________________________________________
process.options = cms.PSet(
verbosity = cms.int64(1)
)
if __name__ == '__main__':
import user
onlineDB=cms.string('sqlite_file:dummy_online.db'),
utilDB=cms.string('sqlite_file:dummy_online.db'),
bufferDB=cms.string('sqlite_file:bufferHV.db'),
tag=cms.string('hv_test'),
onlineAuthentication=cms.string('.'),
sinceYear=cms.int32(2009),
sinceMonth=cms.int32(12),
sinceDay=cms.int32(6),
sinceHour=cms.int32(7),
sinceMinute=cms.int32(0),
sinceSecond=cms.int32(1),
untilYear=cms.int32(2009),
untilMonth=cms.int32(12),
untilDay=cms.int32(6),
untilHour=cms.int32(17),
untilMinute=cms.int32(1),
untilSecond=cms.int32(0),
dumpAtStart=cms.bool(True),
dumpAtEnd=cms.bool(False),
bwdTime=cms.int64(3600),
fwdTime=cms.int64(3600),
minTime=cms.int64(1),
mapVersion=cms.string('hvChannels_090326'),
splitVersion=cms.string('hvChanSplit_090818')),
SinceAppendMode=cms.bool(True),
record=cms.string('DTHVStatusRcd'),
loggingOn=cms.untracked.bool(True),
debug=cms.bool(False))
process.p = cms.Path(process.hv_o2o)
options.register('applyWeights', True,
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.bool, "Apply event weights")
# setup any defaults you want
options.outputFile = 'simpleTest.root'
# options.inputFiles= 'file1.root', 'file2.root'
options.maxEvents = -1 # -1 means all events
# get and parse the command line arguments
options.parseArguments()
process = cms.PSet(treeName=cms.string('stage3flat/probe_tree'),
inputFiles=cms.vstring(options.inputFiles),
outputFile=cms.string(options.outputFile),
maxEvents=cms.int64(options.maxEvents),
monitor=cms.bool(options.monitor),
copyObjects=cms.vstring([]),
applyWeights=cms.bool(options.applyWeights),
weight=cms.string("puW*kfW"),
monitored=cms.vstring(['channel']),
filter=cms.string('isTight1 == 1 && isTight2'))
process.channels = cms.VPSet(
cms.PSet(
name=cms.string('mm'),
selection=cms.string('channel == 0'),
),
cms.PSet(
name=cms.string('me'),
selection=cms.string('channel == 3'),
p_int=cms.int32(2147483647), p_int_untracked=cms.untracked.int32(-2147483647), p_int_opt=cms.int32(0), p_int_optuntracked=cms.untracked.int32(7), vint1=cms.vint32(), vint2=cms.vint32(2147483647), vint3=cms.vint32(2147483647, -2147483647), vint4=cms.vint32(2147483647, -2147483647, 0), uint1=cms.uint32(4294967295), uint2=cms.untracked.uint32(0), vuint1=cms.vuint32(), vuint2=cms.vuint32(4294967295), vuint3=cms.vuint32(4294967295, 0), vuint4=cms.vuint32(4294967295, 0, 11), vuint5=cms.vuint32(4294967295, 0, 11, 21, 31, 41), int64v1=cms.int64(9000000000000000000), int64v2=cms.int64(-9000000000000000000), int64v3=cms.int64(0), vint64v1=cms.vint64(), vint64v2=cms.vint64(9000000000000000000), vint64v3=cms.vint64(9000000000000000000, -9000000000000000000), vint64v4=cms.vint64(9000000000000000000, -9000000000000000000, 0), uint64v1=cms.uint64(18000000000000000000), uint64v2=cms.untracked.uint64(0), vuint64v1=cms.vuint64(), vuint64v2=cms.vuint64(18000000000000000000), vuint64v3=cms.vuint64(18000000000000000000, 0), vuint64v4=cms.vuint64(18000000000000000000, 0, 11), doublev1=cms.double(2.2250738585072014e-308), doublev2=cms.untracked.double(0), doublev3=cms.untracked.double(0.3),
# setup any defaults you want
options.outputFile = 'simpleTest.root'
# options.inputFiles= 'file1.root', 'file2.root'
options.maxEvents = -1 # -1 means all events
# get and parse the command line arguments
options.parseArguments()
process = cms.PSet(
treeName = cms.string('stage3flat/probe_tree'),
inputFiles = cms.vstring(options.inputFiles),
outputFile = cms.string(options.outputFile),
maxEvents = cms.int64(options.maxEvents),
monitor = cms.bool(options.monitor),
copyObjects = cms.vstring([]),
applyWeights = cms.bool(options.applyWeights),
weight = cms.string("puW*kfW"),
monitored = cms.vstring(['channel']),
filter = cms.string('isTight1 == 1 && isTight2')
)
process.channels = cms.VPSet(
cms.PSet(
name = cms.string('mm'),
selection = cms.string('channel == 0'),
),
cms.PSet(
name = cms.string('me'),
# Configuration file for EventSetupTest_t
import FWCore.ParameterSet.Config as cms
process = cms.Process("TEST")
process.source = cms.Source("EmptySource",
numberEventsInRun = cms.untracked.uint32(1),
firstRun = cms.untracked.uint32(1)
)
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(1)
)
process.hv = cms.EDAnalyzer("DTTimeUtility",
year = cms.int32(2009),
month = cms.int32( 12),
day = cms.int32( 11),
hour = cms.int32( 0),
minute = cms.int32( 0),
second = cms.int32( 1),
condTime = cms.int64(5445216652344191208),
coralTime = cms.int64(1278584868675000000) )
process.p = cms.Path(process.hv)
'hgg-2013Final8TeV_ID_s12-zllm50-v7n_skim10k.root'
)
process.outputs = cms.PSet(
fileName = cms.string(output_filename),
treeName = cms.string('hPhotonTree'),
variables = copy(process.inputs.variables),
corrections = cms.VPSet(
correction_pset('ph1.sigietaieta', 'ph1.sigietaieta_corr',
'sieieb_qq_2', 'sieiee_qq_2'),
correction_pset('ph2.sigietaieta', 'ph2.sigietaieta_corr',
'sieieb_qq_2', 'sieiee_qq_2'),
*correction_pset_pair('r9')
)
)
#_______________________________________________________________________________
process.maxEntries = cms.PSet(
toProcess = cms.untracked.int64(-1),
reportEvery = cms.untracked.int64(1000)
)
#_______________________________________________________________________________
process.options = cms.PSet(
verbosity = cms.int64(1)
)
#_______________________________________________________________________________
if __name__ == '__main__':
import user
0 ), vuint4 = cms.vuint32( 4294967295, 0, 11 ), vuint5 = cms.vuint32( 4294967295, 0, 11, 21, 31, 41 ), int64v1 = cms.int64(9000000000000000000), int64v2 = cms.int64(-9000000000000000000), int64v3 = cms.int64(0), vint64v1 = cms.vint64(), vint64v2 = cms.vint64(9000000000000000000), vint64v3 = cms.vint64( 9000000000000000000, -9000000000000000000 ), vint64v4 = cms.vint64( 9000000000000000000, -9000000000000000000, 0 ), uint64v1 = cms.uint64(18000000000000000000), uint64v2 = cms.untracked.uint64(0),
# setup any defaults you want
options.outputFile = 'testAnalysis.root'
options.maxEvents = -1 # -1 means all events
# get and parse the command line arguments
options.parseArguments()
process = cms.PSet(
folder = cms.string('treeproducer'),
treeName = cms.string('hwwStep2'),
inputFiles = cms.vstring(options.inputFiles),
outputFile = cms.string(options.outputFile),
maxEvents = cms.int64(options.maxEvents),
firstEvent = cms.int64(0),
copyHistograms = cms.vstring('llCounters','eeCounters','emCounters','meCounters','mmCounters'),
# preselection cuts
analysisTreeName = cms.string('hwwStep3'),
higgsMass = cms.int32(160),
cutFile = cms.string('test/HmassCuts.cfg'),
jetPtMin = cms.double(30.),
jetPtMinDYcut = cms.double(15.),
jetEtaMax = cms.double(5.),
bDiscriminator = cms.string('trackCountingHighEffBJetTags'),
bThreshold = cms.double(2.1),
minMet = cms.double(20.),
minMll = cms.double(12.),
process.logging = cms.PSet(
logLevel = cms.string("trace"),
enableTimeStamp = cms.bool(True),
enableLogging = cms.bool(True)
)
process.tthMEM = cms.PSet(
isMC = cms.bool(True),
is2016 = cms.bool(True),
treeName = cms.string("tree"),
rleSelectionFile = cms.string(""), # run:lumi:evt numbers (one per line)
pdfName = cms.string("MSTW2008lo68cl"),
madgraphFileName = cms.string("tthAnalysis/tthMEM/data/param_card.dat"),
integrationMode = cms.string("VEGAS"),
maxObjFunctionCalls = cms.uint32(150), # just for testing; proper figure: 100k+
startingFromEntry = cms.int64(0),
debugPlots = cms.uint32(1), # use 0 if no debug plots needed
forceGenLevel = cms.bool(True), # True = dumps gen lvl info regardless of clamping
higgsWidth = cms.double(-1.), # use negative number in case of default H width
clampVariables = cms.VPSet(
cms.PSet( var = cms.string("bCosTheta1"), useGen = cms.bool(False), useCfg = cms.bool(False), val = cms.double(0.0)),
cms.PSet( var = cms.string("bPhi1"), useGen = cms.bool(False), useCfg = cms.bool(False), val = cms.double(0.0)),
cms.PSet( var = cms.string("bCosTheta2"), useGen = cms.bool(False), useCfg = cms.bool(False), val = cms.double(0.0)),
cms.PSet( var = cms.string("bPhi2"), useGen = cms.bool(False), useCfg = cms.bool(False), val = cms.double(0.0)),
cms.PSet( var = cms.string("z1"), useGen = cms.bool(False), useCfg = cms.bool(False), val = cms.double(0.0)),
cms.PSet( var = cms.string("tauPhi"), useGen = cms.bool(False), useCfg = cms.bool(False), val = cms.double(0.0)),
cms.PSet( var = cms.string("tauPhiInv"), useGen = cms.bool(False), useCfg = cms.bool(False), val = cms.double(0.0)),
cms.PSet( var = cms.string("tauMinvSquared"), useGen = cms.bool(False), useCfg = cms.bool(False), val = cms.double(0.0))
),
markovChainParams = cms.PSet( # only read if the integrationMode is set to 'markovchain'
mode = cms.string("uniform"),
#This file knows about all the TTH samples as produced by the VHbb group
#LFN to PFN resolution is done with code in samples_base
#Cross-sections also defined in samples_base
from TTH.MEAnalysis.samples_base import *
import FWCore.ParameterSet.Config as cms
import os
samples_dict = {
"sample": cms.PSet(
name = cms.string("sample"),
nickname = cms.string("sample"),
xSec = cms.double("1"),
nGen = cms.int64(0),
skip = cms.bool(False),
isMC = cms.bool(bool(int(os.environ.get("IS_MC", 1)))),
treeName = cms.string(os.environ.get("INPUT_TREE", "vhbb/tree")),
subFiles = cms.vstring(os.environ["INPUT_FILE"]),
)
}
# The following comments couldn't be translated into the new config version:
# Configuration file for EventSetupTest_t
import FWCore.ParameterSet.Config as cms
process = cms.Process("TEST")
process.source = cms.Source("EmptySource",
numberEventsInRun=cms.untracked.uint32(1),
firstRun=cms.untracked.uint32(1))
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(1))
process.hv = cms.EDAnalyzer("DTTimeUtility",
year=cms.int32(2009),
month=cms.int32(12),
day=cms.int32(11),
hour=cms.int32(0),
minute=cms.int32(0),
second=cms.int32(1),
condTime=cms.int64(5445216652344191208),
coralTime=cms.int64(1278584868675000000))
process.p = cms.Path(process.hv)
'skimmed/hgg-2013Final8TeV_ID_s12-zllm50-v7n_skim10k.root',
#'DoubleMu_Run2011A-05Aug2011-v1_mmgSkim_5kevt.root',
]:
process.inputs.fileNames.append(os.path.join(base_path, filename))
#_______________________________________________________________________________
output_filename = os.path.join(
base_path, 'corrected', 'hgg-2013Final8TeV_ID_s12-zllm50-v7n_skim10k.root')
process.outputs = cms.PSet(fileName=cms.string(output_filename),
treeName=cms.string('hPhotonTree'),
variables=copy(process.inputs.variables),
corrections=cms.VPSet(
correction_pset('ph1.sigietaieta',
'ph1.sigietaieta_corr',
'sieieb_qq_2', 'sieiee_qq_2'),
correction_pset('ph2.sigietaieta',
'ph2.sigietaieta_corr',
'sieieb_qq_2', 'sieiee_qq_2'),
*correction_pset_pair('r9')))
#_______________________________________________________________________________
process.maxEntries = cms.PSet(toProcess=cms.untracked.int64(-1),
reportEvery=cms.untracked.int64(1000))
#_______________________________________________________________________________
process.options = cms.PSet(verbosity=cms.int64(1))
#_______________________________________________________________________________
if __name__ == '__main__':
import user
## output label to give it a unique name
output_label=cms.string("SameSign"),
## suffix to print the plots (before the fit)
## blank means do not print
## available options are: eps, png, pdf
suffix=cms.string("png"),
## pile re-weighting histogram
pileup_hist_file=cms.string(analysis_path +
"/data/puWeights_Summer12_53x_Observed.root"),
pileup_hist_name=cms.string("puWeights"),
## max number of events to run on
max_events=cms.int64(-1),
## verbosity (for trouble shooting)
verbose=cms.bool(False),
## mass range for resonance window
mass_low=cms.double(60.0), # GeV
mass_high=cms.double(120.0), # GeV
mass_bin_width=cms.double(2.0), # GeV
# datasets to run on
datasets=cms.VPSet(dy_full, single_el),
## bins for the observables
## supported pt, eta, phi, and # vertices
## note: for eta and phi, no negative bins means use |eta| and |phi|, respectively
activeBranches = copy.deepcopy(activeBranchesMC),
)
process.outputs = cms.PSet(
outputName = cms.string(sample.name + '.root')
)
process.maxEvents = cms.PSet(
toProcess = cms.untracked.int64(int(-1)),
reportEvery = cms.untracked.int64(10000)
)
## Histograms configuration
process.histograms = cms.PSet(
isMC = cms.bool(True),
#NoSelection = copy.deepcopy(noSelection ),
FullSelection = copy.deepcopy(fullSelection),
)
#______________________________________________________________________________
## Event weight
process.eventWeight = cms.PSet(
crossSectionInPicoBarns = cms.double(sample.cross_section_in_pb),
scaleToLumiInInverseFemtoBarns = cms.double(1e-3 * data.lumi_per_pb),
totalProcessedEvents = cms.int64(sample.total_processed_events),
)
if __name__ == '__main__':
import user
lepton_type = cms.string("electron"),
## output label to give it a unique name
output_label = cms.string("SameSign"),
## suffix to print the plots (before the fit)
## blank means do not print
## available options are: eps, png, pdf
suffix = cms.string("png"),
## pile re-weighting histogram
pileup_hist_file = cms.string(analysis_path + "/data/puWeights_Summer12_53x_Observed.root"),
pileup_hist_name = cms.string("puWeights"),
## max number of events to run on
max_events = cms.int64(-1),
## verbosity (for trouble shooting)
verbose = cms.bool(False),
## mass range for resonance window
mass_low = cms.double(60.0), # GeV
mass_high = cms.double(120.0), # GeV
mass_bin_width = cms.double(2.0), # GeV
# datasets to run on
datasets = cms.VPSet(dy_full, single_el),
## bins for the observables
## supported pt, eta, phi, and # vertices
## note: for eta and phi, no negative bins means use |eta| and |phi|, respectively
#This file knows about all the TTH samples as produced by the VHbb group
#LFN to PFN resolution is done with code in samples_base
#Cross-sections also defined in samples_base
from TTH.MEAnalysis.samples_base import *
import FWCore.ParameterSet.Config as cms
import os
samples_dict = {
"sample":
cms.PSet(
name=cms.string("sample"),
nickname=cms.string("sample"),
xSec=cms.double("1"),
nGen=cms.int64(0),
skip=cms.bool(False),
isMC=cms.bool(bool(int(os.environ.get("IS_MC", 1)))),
treeName=cms.string(os.environ.get("INPUT_TREE", "vhbb/tree")),
subFiles=cms.vstring(os.environ["INPUT_FILE"]),
)
}
