MergeFilesAndCalculateEfficiencies_cfg.py

Go to the documentation of this file.

#!/usr/bin/env python
"""
 
     MergeFilesAndCalculateEfficiencies.py
 
        Merges multiple root files containing the numerator and denominator histograms produced by the
        TauTagValidation package.  The efficiency (num/denominator) is then computed 
        as defined Validation/RecoTau/python/RecoTauValidation_cff.py and stored in OutputFile_Eff.root
 
     Usage: cmsRun MergeFilesAndCalculateEfficiencies.py OutputFile InputFiles
 
     Example: ./MergeFilesAndCalculateEfficiencies.py CMSSW_3_1_0_Signal.root CMSSW_3_1_0_ZTT_*.root
 
"""
from __future__ import print_function
 
import os
import sys
import FWCore.ParameterSet.Config as cms
from Validation.RecoTau.ValidationOptions_cff import allowedOptions
 
if len(sys.argv) < 5:
   print("Error. Expected at least 3 arguments\n\nUsage: MergeFilesAndCalculateEfficiencies.py dataType OutputFile InputFileGlob")
   sys.exit()
 
dataType   = sys.argv[2]
OutputFile = sys.argv[3]
Inputs     = sys.argv[4:]
 
if not dataType in allowedOptions['eventType']:
   print("Error. The first argument must be the dataType. Types availables are:")
   print(allowedOptions['eventType'])
   sys.exit()
 
for aFile in Inputs:
   if not os.path.exists(aFile):
      print("Input file %s does not exist!" % aFile)
      sys.exit()
 
if os.path.exists(OutputFile):
   GotGoodValue = False
   userInput = ""
   while not GotGoodValue:
      userInput = raw_input("Output file %s exists; replace it? [yn] " % OutputFile).strip()
      if userInput != 'y' and userInput != 'n':
         print("Please enter y or n")
      else:
         GotGoodValue = True
   if userInput == 'n':
      sys.exit()
 
# Merge files using hadd utility
commandString  = "hadd -f %s " % OutputFile
for aFile in Inputs:
   commandString += aFile
   commandString += " "
 
os.system(commandString)
 
print("Running cmsRun command to generate efficiencies")
 
process = cms.Process("TEST")
 
process.source = cms.Source("EmptySource")
 
process.maxEvents = cms.untracked.PSet(
    input = cms.untracked.int32(1)
)
 
process.DQMStore = cms.Service("DQMStore")
process.load("Validation.RecoTau.dataTypes.ValidateTausOn%s_cff"%dataType)
 
process.loadFile   = cms.EDAnalyzer("TauDQMFileLoader",
      myFiles = cms.PSet(
         inputFileNames = cms.vstring(OutputFile),
         scaleFactor = cms.double(1.),
         )
)
 
process.saveTauEff = cms.EDAnalyzer("TauDQMSimpleFileSaver",
#  outputFileName = cms.string(OutputFile.replace('.root', '_Eff.root'))
  outputFileName = cms.string(OutputFile)
)
 
process.p = cms.Path(
      process.loadFile*
      getattr(process,'TauEfficiencies%s'%dataType)*
      process.saveTauEff
      )