'EFCalo_isLHVLooseCaloOnly_rel21_20170217&HLT_isLHVLoose_rel21_20170217',
'EFCalo_isLHLooseCaloOnly_rel21_20170217&HLT_isLHLoose_rel21_20170217',
'EFCalo_isLHMediumCaloOnly_rel21_20170217&HLT_isLHMedium_rel21_20170217',
'EFCalo_isLHTightCaloOnly_rel21_20170217&HLT_isLHTight_rel21_20170217',
'EFCalo_isRingerVLoose_v5&HLT_isLHVLoose_rel21_20170217',
'EFCalo_isRingerLoose_v5&HLT_isLHLoose_rel21_20170217',
'EFCalo_isRingerMedium_v5&HLT_isLHMedium_rel21_20170217',
'EFCalo_isRingerTight_v5&HLT_isLHTight_rel21_20170217',
]
############################################################################################
# First event Looper
eventLooper1 = EventLooper(inputFiles=sgnData,
treePath='HLT/Egamma/Expert/support/probes',
dataframe=DataframeEnum.PhysVal,
nov=NOV,
outputFile=histoName,
level=level)
# Second Event Looper
eventLooper2 = EventLooper(inputFiles=bkgData,
treePath='HLT/Egamma/Expert/support/fakes',
dataframe=DataframeEnum.PhysVal,
outputFile=histoName,
nov=NOV,
level=level)
############################################################################################
algZ = EventSelection('EventSelectionZ')
algZ.selectionZ = True
############################################################################################
histoName = 'fudge_MC.root'
NOV = -1
level = LoggingLevel.INFO
mcData = '/home/jodafons/CERN-DATA/Online/data/mc15_13TeV/user.jodafons.mc15_13TeV.361106.Zee.merge.SelectionZee.PhysVal.r0005_GLOBAL/'
ppData = '/home/jodafons/CERN-DATA/Online/data/data16_13TeV/user.jodafons.data16_13TeV.periodG.physicsMain.p3013.SelectionData.PhysVal.r0005_GLOBAL/'
############################################################################################
# First event Looper
eventLooper1 = EventLooper(inputFiles=mcData,
treePath='*/HLT/Egamma/Expert/support/probes',
dataframe=DataframeEnum.PhysVal,
nov=-1,
outputFile=histoName,
level=level)
# Second Event Looper
eventLooper2 = EventLooper(inputFiles=ppData,
treePath='*/HLT/Egamma/Expert/support/probes',
dataframe=DataframeEnum.PhysVal,
outputFile=histoName,
nov=0,
level=level)
############################################################################################
algZ = EventSelection('EventSelectionZ')
algZ.selectionZ = True
help="Use this to switch to Monte Carlo mode")
import sys, os
if len(sys.argv) == 1:
parser.print_help()
sys.exit(1)
args = parser.parse_args()
nov = int(args.nov)
############################################################################################
level = LoggingLevel.INFO
# First event Looper
eventLooper1 = EventLooper(inputFiles=args.inputFiles,
treePath='*/HLT/Egamma/Expert/support/probes',
dataframe=DataframeEnum.PhysVal,
nov=nov,
outputFile=args.outputFile,
level=level)
# Second Event Looper
eventLooper2 = EventLooper(inputFiles=args.inputFiles,
treePath='*/HLT/Egamma/Expert/support/fakes',
dataframe=DataframeEnum.PhysVal,
outputFile=args.outputFile,
nov=nov,
level=level)
############################################################################################
from TrigEgammaDevelopments.tools import EventSelection, QuadrantTool
args = parser.parse_args()
# treat data input
# convert to int
relax = {'Tight':args.relax[0],'Medium':args.relax[1],\
'Loose':args.relax[2],'VLoose':args.relax[3]}
print relax
############################################################################################
level = LoggingLevel.INFO
# First event Looper
eventLooperZ = EventLooper(inputFiles=args.inputFiles,
treePath='*/HLT/Egamma/Expert/support/probes',
dataframe=DataframeEnum.PhysVal,
outputFile=args.outputFile,
nov=args.nov,
level=level)
# Second Event Looper
eventLooperF = EventLooper(inputFiles=args.inputFiles,
treePath='*/HLT/Egamma/Expert/support/fakes',
dataframe=DataframeEnum.PhysVal,
outputFile=args.outputFile,
nov=args.nov,
level=level)
############################################################################################
algZ = EventSelection('EventSelectionZ')
algZ.selectionZ = True
NOV = 5000
#NOV = -1
level = LoggingLevel.DEBUG
#level = LoggingLevel.INFO
#sgnData = '/home/wsfreund/CERN-DATA/Offline/skimmedNtuple/user.wsfreund.mc14_Zee_mcsel_skimntuple_offline_lhcalo_refs_SkimmedNtuple.root'
#bkgData = '/home/wsfreund/CERN-DATA/Offline/skimmedNtuple/user.wsfreund.mc14_JF17_mcsel_skimntuple_offline_lhcalo_refs_SkimmedNtuple.root'
sgnData = '/home/jodafons/CERN-DATA/Online/data/mc15_13TeV/Zee/361106/user.jodafons.mc15_13TeV.361106.Zee.merge.SelectionZee.PhysVal.r0002_GLOBAL/'
bkgData = '/home/jodafons/CERN-DATA/Online/data/mc15_13TeV/JF17/423300/user.jodafons.mc15_13TeV.423300.JF17.SelectionFakes.PhysVal.r0002_GLOBAL/'
############################################################################################
# First event Looper
eventLooper1 = EventLooper(inputFiles=sgnData,
treePath='HLT/Egamma/Expert/support/probes',
dataframe=DataframeEnum.PhysVal,
nov=NOV,
level=level)
# Second Event Looper
eventLooper2 = EventLooper(inputFiles=bkgData,
treePath='HLT/Egamma/Expert/support/fakes',
dataframe=DataframeEnum.PhysVal,
nov=NOV,
level=level)
############################################################################################
algZ = EventSelection('EventSelectionZ')
algZ.selectionZ = True
algZ.selectionFakes = False
help="The offline et cut")
import sys, os
if len(sys.argv) == 1:
parser.print_help()
sys.exit(1)
args = parser.parse_args()
############################################################################################
level = LoggingLevel.INFO
# First event Looper
eventLooper = EventLooper( inputFiles = args.inputFiles,
treePath = '*/HLT/Egamma/Expert/support/fakes' if args.isFakes else\
'*/HLT/Egamma/Expert/support/probes',
dataframe = DataframeEnum.PhysVal,
nov = args.nov,
outputFile = args.outputFile,
level = level)
############################################################################################
from TrigEgammaDevelopments.tools import EventSelection
algE = EventSelection('EventSelection')
algE.selectionZ = False if args.isFakes else True
algE.selectionFakes = True if args.isFakes else False
algE.doTrigger = True
algE.l2EtCut = args.offEtCut
algE.offEtCut = args.offEtCut
algE.selectionMC = True if args.isMC else False