Python doSystematics Examples

Example #1

algsToRun["postselectTwoLepton"] = ROOT.PostselectTwoLeptonEvents()

#todo move the enums to a separate file since they are shared by multiple algs
algsToRun["calculateRJigsawVariables"] = ROOT.CalculateRJigsawVariables()
algsToRun[
    "calculateRJigsawVariables"].calculatorName = ROOT.CalculateRJigsawVariables.tlsCalculator
algsToRun["calculateRegionVars"] = ROOT.CalculateRegionVars()
algsToRun[
    "calculateRegionVars"].calculatorName = ROOT.CalculateRegionVars.tlsCalculator

for regionName in ["SR"]:
    tmpWriteOutputNtuple = ROOT.WriteOutputNtuple()
    tmpWriteOutputNtuple.outputName = outputFilename
    tmpWriteOutputNtuple.regionName = regionName
    tmpWriteOutputNtuple.systName = ""
    algsToRun["writeOutputNtuple_" + regionName] = tmpWriteOutputNtuple

if options.doSystematics: commonOptions.doSystematics(algsToRun)

job.outputAdd(output)

commonOptions.addAlgsFromDict(job, algsToRun, options.verbosity)

if options.nevents > 0:
    logging.info("Running " + str(options.nevents) + " events")
    job.options().setDouble(ROOT.EL.Job.optMaxEvents, float(options.nevents))

commonOptions.overwriteSubmitDir(options.submitDir, options.doOverwrite)
commonOptions.submitJob(job, options.driver, options.submitDir,
                        options.gridUser, options.gridTag)

Example #2

File: run_lvlv.py Project: JOliver-GH/FactoryTools

algsToRun["preselectDileptonicWW"  ]    = ROOT.PreselectDileptonicWWEvents()
algsToRun["selectDileptonicWW"     ]    = ROOT.SelectDileptonicWWEvents()
algsToRun["postselectDileptonicWW" ]    = ROOT.PostselectDileptonicWWEvents()

#todo move the enums to a separate file since they are shared by multiple algs
algsToRun["calculateRJigsawVariables" ]                = ROOT.CalculateRJigsawVariables()
algsToRun["calculateRJigsawVariables" ].calculatorName = ROOT.CalculateRJigsawVariables.lvlvCalculator
algsToRun["calculateRegionVars" ]                      = ROOT.CalculateRegionVars()
algsToRun["calculateRegionVars" ].calculatorName       = ROOT.CalculateRegionVars.lvlvCalculator

for regionName in ["SR","CR1L","CR0L"]:
    tmpWriteOutputNtuple                     = ROOT.WriteOutputNtuple()
    tmpWriteOutputNtuple.outputName          = outputFilename
    tmpWriteOutputNtuple.regionName          = regionName
    tmpWriteOutputNtuple.systName            = ""
    algsToRun["writeOutputNtuple"+regionName] = tmpWriteOutputNtuple

if options.doSystematics : commonOptions.doSystematics(algsToRun)

job.outputAdd(output);

commonOptions.addAlgsFromDict(job , algsToRun , options.verbosity)

if options.nevents > 0 :
    logging.info("Running " + str(options.nevents) + " events")
    job.options().setDouble (ROOT.EL.Job.optMaxEvents, float(options.nevents));

commonOptions.overwriteSubmitDir(options.submitDir , options.doOverwrite)
commonOptions.submitJob         ( job , options.driver , options.submitDir , options.gridUser , options.gridTag)

Example #3

File: run_zl.py Project: rsmith54/FactoryTools

# These are the calculators that calculate various derived quantities
algsToRun["calculateRJigsawVariables"] = ROOT.CalculateRJigsawVariables()
algsToRun["calculateRJigsawVariables"].calculatorName = ROOT.CalculateRJigsawVariables.zlCalculator
algsToRun["calculateRegionVars"]                      = ROOT.CalculateRegionVars()
algsToRun["calculateRegionVars"].calculatorName       = ROOT.CalculateRegionVars.zlCalculator

# This bit runs filtering on derived variables. e.g. MEff filter.
algsToRun["postselectZeroLepton"]    = ROOT.PostselectZeroLeptonEvents()

# These correspond to writing out the various trees used in the analysis
for regionName in ["SR","CR1L","CR2L","CRY"]:
    tmpWriteOutputNtuple                       = ROOT.WriteOutputNtuple()
    tmpWriteOutputNtuple.outputName            = outputFilename
    tmpWriteOutputNtuple.regionName            = regionName
    # tmpWriteOutputNtuple.systVar            = 0
    algsToRun["writeOutputNtuple_"+regionName] = tmpWriteOutputNtuple

if options.doSystematics : 
	algsToRun = commonOptions.doSystematics(algsToRun,fullChainOnWeightSysts = 0, excludeStrings = ["JET_Rtrk_","TAUS_"])

job.outputAdd(output);
commonOptions.addAlgsFromDict(job , algsToRun , options.verbosity)

if options.nevents > 0 :
    logging.info("Running " + str(options.nevents) + " events")
    job.options().setDouble (ROOT.EL.Job.optMaxEvents, float(options.nevents));

commonOptions.overwriteSubmitDir(options.submitDir , options.doOverwrite)
commonOptions.submitJob         ( job , options.driver , options.submitDir , options.gridUser , options.gridTag)