def __init__(self, year, tree, terms=None, refantitau=True, verbose=False, very_verbose=False, **kwargs): super(METRecalculation, self).__init__(**kwargs) self.terms = terms or set() self.year = year self.tree = tree self.refantitau = refantitau self.verbose = verbose self.very_verbose = very_verbose # Initialise your METUtility object # https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/MissingETUtility self.tool = METUtility() # configure self.tool.configMissingET( year == 2012, # is 2012 year == 2012) # is STVF if year == 2011: # In 2012, we always rebuild MET from AOD, so the RefMuon term is # not computed. However, for 2011 data, you should activate the # RefMuon term and fill the term. # These are the terms required for MET_RefFinal(_BDTMedium) self.tool.defineMissingET( True, # RefEle True, # RefGamma True, # RefTau True, # RefJet True, # RefMuon True, # MuonTotal True, # Soft ) # The threshold below which jets enter the SoftJets term (JES is not applied) #self.tool.setSoftJetCut(20e3) # Whether to use MUID muons (otherwise STACO). self.tool.setIsMuid(False) # Whether METUtility should scream at you over every little thing self.tool.setVerbosity(self.very_verbose)