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)
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)
class METRecalculation(EventFilter):
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)
def get_met(self):
util = self.tool
multisyst = METUtil.MultiSyst()
if Systematics.MET_SCALESOFTTERMS_UP in self.terms:
multisyst.setSyst(Systematics.MET_SCALESOFTTERMS_UP)
return util.getMissingET(METUtil.RefFinal, multisyst)
elif Systematics.MET_SCALESOFTTERMS_DOWN in self.terms:
multisyst.setSyst(Systematics.MET_SCALESOFTTERMS_DOWN)
return util.getMissingET(METUtil.RefFinal, multisyst)
elif Systematics.MET_RESOSOFTTERMS_UP in self.terms:
multisyst.setSyst(Systematics.MET_RESOSOFTTERMS_UP)
return util.getMissingET(METUtil.RefFinal, multisyst)
elif Systematics.MET_RESOSOFTTERMS_DOWN in self.terms:
multisyst.setSyst(Systematics.MET_RESOSOFTTERMS_DOWN)
return util.getMissingET(METUtil.RefFinal, multisyst)
return util.getMissingET(METUtil.RefFinal)
def passes(self, event):
# reset the METUtility
self.tool.reset()
# These set up the systematic "SoftTerms_ptHard"
self.tool.setNvtx(self.tree.nvtxsoftmet)
if self.year == 2012:
return self.passes_12(event)
elif self.year == 2011:
return self.passes_11(event)
else:
raise ValueError("No MET defined for year %d" % self.year)
def passes_11(self, event):
"""
JETS
Always use setJetParameters since they may be recalibrated upstream
https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/MissingETUtilityFAQ#If_I_recalibrate_correct_my_anal
"""
# Due to the jet weights being encoded slightly differently in 2011
# reconstruction, you need to provide one extra argument, i.e. the jet
# pT directly from the D3PD, so that the correct jets can be included
# in the RefJet term, while soft jets (pT<20) will not be.
self.tool.setJetParameters(
event.jet_pt, # recalibrated pT
self.tree.jet_eta_original,
self.tree.jet_phi_original,
event.jet_E,
event.jet_AntiKt4LCTopo_MET_BDTMedium_wet,
event.jet_AntiKt4LCTopo_MET_BDTMedium_wpx,
event.jet_AntiKt4LCTopo_MET_BDTMedium_wpy,
event.jet_AntiKt4LCTopo_MET_BDTMedium_statusWord,
self.tree.jet_pt_original) # <== extra argument for 2011
# Taus
self.tool.setTauParameters(
event.tau_pt,
event.tau_eta,
event.tau_phi,
event.tau_MET_BDTMedium_wet,
event.tau_MET_BDTMedium_wpx,
event.tau_MET_BDTMedium_wpy,
event.tau_MET_BDTMedium_statusWord)
# Electrons
self.tool.setMETTerm(
METUtil.RefEle,
event.MET_RefEle_BDTMedium_etx,
event.MET_RefEle_BDTMedium_ety,
event.MET_RefEle_BDTMedium_sumet)
# Muons
self.tool.setMETTerm(
METUtil.MuonTotal,
event.MET_Muon_Total_Staco_BDTMedium_etx,
event.MET_Muon_Total_Staco_BDTMedium_ety,
event.MET_Muon_Total_Staco_BDTMedium_sumet)
# Muons
# Note that RefMuon is not rebuilt from muons
# -- it is a calorimeter term.
self.tool.setMETTerm(
METUtil.RefMuon,
event.MET_RefMuon_Staco_BDTMedium_etx,
event.MET_RefMuon_Staco_BDTMedium_ety,
event.MET_RefMuon_Staco_BDTMedium_sumet)
# Photons
self.tool.setMETTerm(
METUtil.RefGamma,
event.MET_RefGamma_BDTMedium_etx,
event.MET_RefGamma_BDTMedium_ety,
event.MET_RefGamma_BDTMedium_sumet)
# Soft terms
# In most 2011 samples, the SoftJets and CellOut components of the MET
# are provided as separate pieces, but we since moved to a combined
# treatment of both. Hence, you will need to manually sum the two
# components when providing them to METUtility, like so:
self.tool.setMETTerm(
METUtil.SoftTerms,
event.MET_SoftJets_BDTMedium_etx + event.MET_CellOut_BDTMedium_etx,
event.MET_SoftJets_BDTMedium_ety + event.MET_CellOut_BDTMedium_ety,
event.MET_SoftJets_BDTMedium_sumet + event.MET_CellOut_BDTMedium_sumet)
MET = self.get_met()
if self.verbose:
log.info("Run: {0} Event: {1}".format(
event.RunNumber,
event.EventNumber))
log.info("Recalculated MET: %.3f (original: %.3f)" % (
MET.et(), event.MET_RefFinal_BDTMedium_et))
log.info("Recalculated MET phi: %.3f (original: %.3f)" % (
MET.phi(), event.MET_RefFinal_BDTMedium_phi))
if (abs(MET.et() - event.MET_RefFinal_BDTMedium_et) /
event.MET_RefFinal_BDTMedium_et) > 0.1:
log.warning("Large MET difference!")
# update the MET with the shifted value
self.tree.MET_etx_original = event.MET_RefFinal_BDTMedium_etx
self.tree.MET_ety_original = event.MET_RefFinal_BDTMedium_ety
self.tree.MET_et_original = event.MET_RefFinal_BDTMedium_et
self.tree.MET_sumet_original = event.MET_RefFinal_BDTMedium_sumet
self.tree.MET_phi_original = event.MET_RefFinal_BDTMedium_phi
event.MET_RefFinal_BDTMedium_etx = MET.etx()
event.MET_RefFinal_BDTMedium_ety = MET.ety()
event.MET_RefFinal_BDTMedium_et = MET.et()
event.MET_RefFinal_BDTMedium_sumet = MET.sumet()
event.MET_RefFinal_BDTMedium_phi = MET.phi()
return True
def passes_12(self, event):
if self.refantitau:
# AntiTau MET calculation from Alex Tuna
# If a selected tau matches a JVF jet, clear the corresponding jet
# weights and set the tau weights to 1.0. This must be applied
# after the tau selection but before the jet selection
assert(len(event.taus) == 2)
for tau in event.taus:
# event.taus only contains selected taus at this point
match = False
for jet in event.jets:
# event.jets contains all jets
# Does this jet match the tau?
if dR(tau.eta, tau.phi, jet.eta, jet.phi) < 0.4:
match = True
# Loop through subjets to find the JVF jets used for STVF
for k in xrange(jet.AntiKt4LCTopo_MET_wet.size()):
# If the subjet is a JVF jet, set weight to 0
if jet.AntiKt4LCTopo_MET_statusWord[k] == (0x3300 | 0x0001):
jet.AntiKt4LCTopo_MET_wet[k] = 0.0
jet.AntiKt4LCTopo_MET_wpx[k] = 0.0
jet.AntiKt4LCTopo_MET_wpy[k] = 0.0
# If the tau has a matching jet, set the tau weights to 1
if match:
tau.MET_statusWord[0] = 1
tau.MET_wet[0] = 1.0
tau.MET_wpx[0] = 1.0
tau.MET_wpy[0] = 1.0
# this must be put before setting the jets parameters
self.tool.setJetPUcode(MissingETTags.JPU_JET_JVFCUT)
# JETS
# Always use setJetParameters since they may be recalibrated upstream
# https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/MissingETUtilityFAQ#If_I_recalibrate_correct_my_anal
self.tool.setJetParameters(
event.jet_pt,
self.tree.jet_eta_original,
self.tree.jet_phi_original,
event.jet_E,
event.jet_AntiKt4LCTopo_MET_wet,
event.jet_AntiKt4LCTopo_MET_wpx,
event.jet_AntiKt4LCTopo_MET_wpy,
event.jet_AntiKt4LCTopo_MET_statusWord)
# Taus
self.tool.setTauParameters(
event.tau_pt,
event.tau_eta,
event.tau_phi,
event.tau_MET_wet,
event.tau_MET_wpx,
event.tau_MET_wpy,
event.tau_MET_statusWord)
# Electrons
self.tool.setMETTerm(
METUtil.RefEle,
event.MET_RefEle_etx,
event.MET_RefEle_ety,
event.MET_RefEle_sumet)
# Muons
self.tool.setMETTerm(
METUtil.MuonTotal,
event.MET_Muon_Total_Staco_etx,
event.MET_Muon_Total_Staco_ety,
event.MET_Muon_Total_Staco_sumet)
# Photons
self.tool.setMETTerm(
METUtil.RefGamma,
event.MET_RefGamma_etx,
event.MET_RefGamma_ety,
event.MET_RefGamma_sumet)
# Soft terms
self.tool.setMETTerm(
METUtil.SoftTerms,
event.MET_CellOut_Eflow_STVF_etx,
event.MET_CellOut_Eflow_STVF_ety,
event.MET_CellOut_Eflow_STVF_sumet)
MET = self.get_met()
if self.verbose:
log.info("Run: {0} Event: {1}".format(
event.RunNumber,
event.EventNumber))
log.info("Recalculated MET: %.3f (original: %.3f)" % (
MET.et(), event.MET_RefFinal_STVF_et))
log.info("Recalculated MET phi: %.3f (original: %.3f)" % (
MET.phi(), event.MET_RefFinal_STVF_phi))
if (abs(MET.et() - event.MET_RefFinal_STVF_et) /
event.MET_RefFinal_STVF_et) > 0.1:
log.warning("Large MET difference!")
# update the MET with the shifted value
self.tree.MET_etx_original = event.MET_RefFinal_STVF_etx
self.tree.MET_ety_original = event.MET_RefFinal_STVF_ety
self.tree.MET_et_original = event.MET_RefFinal_STVF_et
self.tree.MET_sumet_original = event.MET_RefFinal_STVF_sumet
self.tree.MET_phi_original = event.MET_RefFinal_STVF_phi
event.MET_RefFinal_STVF_etx = MET.etx()
event.MET_RefFinal_STVF_ety = MET.ety()
event.MET_RefFinal_STVF_et = MET.et()
event.MET_RefFinal_STVF_sumet = MET.sumet()
event.MET_RefFinal_STVF_phi = MET.phi()
return True
class METRecalculation(EventFilter):
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)
def get_met(self):
util = self.tool
multisyst = METUtil.MultiSyst()
if Systematics.MET_SCALESOFTTERMS_UP in self.terms:
multisyst.setSyst(Systematics.MET_SCALESOFTTERMS_UP)
return util.getMissingET(METUtil.RefFinal, multisyst)
elif Systematics.MET_SCALESOFTTERMS_DOWN in self.terms:
multisyst.setSyst(Systematics.MET_SCALESOFTTERMS_DOWN)
return util.getMissingET(METUtil.RefFinal, multisyst)
elif Systematics.MET_RESOSOFTTERMS_UP in self.terms:
multisyst.setSyst(Systematics.MET_RESOSOFTTERMS_UP)
return util.getMissingET(METUtil.RefFinal, multisyst)
elif Systematics.MET_RESOSOFTTERMS_DOWN in self.terms:
multisyst.setSyst(Systematics.MET_RESOSOFTTERMS_DOWN)
return util.getMissingET(METUtil.RefFinal, multisyst)
return util.getMissingET(METUtil.RefFinal)
def passes(self, event):
# reset the METUtility
self.tool.reset()
# These set up the systematic "SoftTerms_ptHard"
self.tool.setNvtx(self.tree.nvtxsoftmet)
if self.year == 2012:
return self.passes_12(event)
elif self.year == 2011:
return self.passes_11(event)
else:
raise ValueError("No MET defined for year %d" % self.year)
def passes_11(self, event):
"""
JETS
Always use setJetParameters since they may be recalibrated upstream
https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/MissingETUtilityFAQ#If_I_recalibrate_correct_my_anal
"""
# Due to the jet weights being encoded slightly differently in 2011
# reconstruction, you need to provide one extra argument, i.e. the jet
# pT directly from the D3PD, so that the correct jets can be included
# in the RefJet term, while soft jets (pT<20) will not be.
self.tool.setJetParameters(
event.jet_pt, # recalibrated pT
self.tree.jet_eta_original,
self.tree.jet_phi_original,
event.jet_E,
event.jet_AntiKt4LCTopo_MET_BDTMedium_wet,
event.jet_AntiKt4LCTopo_MET_BDTMedium_wpx,
event.jet_AntiKt4LCTopo_MET_BDTMedium_wpy,
event.jet_AntiKt4LCTopo_MET_BDTMedium_statusWord,
self.tree.jet_pt_original) # <== extra argument for 2011
# Taus
self.tool.setTauParameters(event.tau_pt, event.tau_eta, event.tau_phi,
event.tau_MET_BDTMedium_wet,
event.tau_MET_BDTMedium_wpx,
event.tau_MET_BDTMedium_wpy,
event.tau_MET_BDTMedium_statusWord)
# Electrons
self.tool.setMETTerm(METUtil.RefEle, event.MET_RefEle_BDTMedium_etx,
event.MET_RefEle_BDTMedium_ety,
event.MET_RefEle_BDTMedium_sumet)
# Muons
self.tool.setMETTerm(METUtil.MuonTotal,
event.MET_Muon_Total_Staco_BDTMedium_etx,
event.MET_Muon_Total_Staco_BDTMedium_ety,
event.MET_Muon_Total_Staco_BDTMedium_sumet)
# Muons
# Note that RefMuon is not rebuilt from muons
# -- it is a calorimeter term.
self.tool.setMETTerm(METUtil.RefMuon,
event.MET_RefMuon_Staco_BDTMedium_etx,
event.MET_RefMuon_Staco_BDTMedium_ety,
event.MET_RefMuon_Staco_BDTMedium_sumet)
# Photons
self.tool.setMETTerm(METUtil.RefGamma,
event.MET_RefGamma_BDTMedium_etx,
event.MET_RefGamma_BDTMedium_ety,
event.MET_RefGamma_BDTMedium_sumet)
# Soft terms
# In most 2011 samples, the SoftJets and CellOut components of the MET
# are provided as separate pieces, but we since moved to a combined
# treatment of both. Hence, you will need to manually sum the two
# components when providing them to METUtility, like so:
self.tool.setMETTerm(
METUtil.SoftTerms,
event.MET_SoftJets_BDTMedium_etx + event.MET_CellOut_BDTMedium_etx,
event.MET_SoftJets_BDTMedium_ety + event.MET_CellOut_BDTMedium_ety,
event.MET_SoftJets_BDTMedium_sumet +
event.MET_CellOut_BDTMedium_sumet)
MET = self.get_met()
if self.verbose:
log.info("Run: {0} Event: {1}".format(event.RunNumber,
event.EventNumber))
log.info("Recalculated MET: %.3f (original: %.3f)" %
(MET.et(), event.MET_RefFinal_BDTMedium_et))
log.info("Recalculated MET phi: %.3f (original: %.3f)" %
(MET.phi(), event.MET_RefFinal_BDTMedium_phi))
if (abs(MET.et() - event.MET_RefFinal_BDTMedium_et) /
event.MET_RefFinal_BDTMedium_et) > 0.1:
log.warning("Large MET difference!")
# update the MET with the shifted value
self.tree.MET_etx_original = event.MET_RefFinal_BDTMedium_etx
self.tree.MET_ety_original = event.MET_RefFinal_BDTMedium_ety
self.tree.MET_et_original = event.MET_RefFinal_BDTMedium_et
self.tree.MET_sumet_original = event.MET_RefFinal_BDTMedium_sumet
self.tree.MET_phi_original = event.MET_RefFinal_BDTMedium_phi
event.MET_RefFinal_BDTMedium_etx = MET.etx()
event.MET_RefFinal_BDTMedium_ety = MET.ety()
event.MET_RefFinal_BDTMedium_et = MET.et()
event.MET_RefFinal_BDTMedium_sumet = MET.sumet()
event.MET_RefFinal_BDTMedium_phi = MET.phi()
return True
def passes_12(self, event):
if self.refantitau:
# AntiTau MET calculation from Alex Tuna
# If a selected tau matches a JVF jet, clear the corresponding jet
# weights and set the tau weights to 1.0. This must be applied
# after the tau selection but before the jet selection
assert (len(event.taus) == 2)
for tau in event.taus:
# event.taus only contains selected taus at this point
match = False
for jet in event.jets:
# event.jets contains all jets
# Does this jet match the tau?
if dR(tau.eta, tau.phi, jet.eta, jet.phi) < 0.4:
match = True
# Loop through subjets to find the JVF jets used for STVF
for k in xrange(jet.AntiKt4LCTopo_MET_wet.size()):
# If the subjet is a JVF jet, set weight to 0
if jet.AntiKt4LCTopo_MET_statusWord[k] == (0x3300 |
0x0001):
jet.AntiKt4LCTopo_MET_wet[k] = 0.0
jet.AntiKt4LCTopo_MET_wpx[k] = 0.0
jet.AntiKt4LCTopo_MET_wpy[k] = 0.0
# If the tau has a matching jet, set the tau weights to 1
if match:
tau.MET_statusWord[0] = 1
tau.MET_wet[0] = 1.0
tau.MET_wpx[0] = 1.0
tau.MET_wpy[0] = 1.0
# this must be put before setting the jets parameters
self.tool.setJetPUcode(MissingETTags.JPU_JET_JVFCUT)
# JETS
# Always use setJetParameters since they may be recalibrated upstream
# https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/MissingETUtilityFAQ#If_I_recalibrate_correct_my_anal
self.tool.setJetParameters(event.jet_pt, self.tree.jet_eta_original,
self.tree.jet_phi_original, event.jet_E,
event.jet_AntiKt4LCTopo_MET_wet,
event.jet_AntiKt4LCTopo_MET_wpx,
event.jet_AntiKt4LCTopo_MET_wpy,
event.jet_AntiKt4LCTopo_MET_statusWord)
# Taus
self.tool.setTauParameters(event.tau_pt, event.tau_eta, event.tau_phi,
event.tau_MET_wet, event.tau_MET_wpx,
event.tau_MET_wpy, event.tau_MET_statusWord)
# Electrons
self.tool.setMETTerm(METUtil.RefEle, event.MET_RefEle_etx,
event.MET_RefEle_ety, event.MET_RefEle_sumet)
# Muons
self.tool.setMETTerm(METUtil.MuonTotal, event.MET_Muon_Total_Staco_etx,
event.MET_Muon_Total_Staco_ety,
event.MET_Muon_Total_Staco_sumet)
# Photons
self.tool.setMETTerm(METUtil.RefGamma, event.MET_RefGamma_etx,
event.MET_RefGamma_ety, event.MET_RefGamma_sumet)
# Soft terms
self.tool.setMETTerm(METUtil.SoftTerms,
event.MET_CellOut_Eflow_STVF_etx,
event.MET_CellOut_Eflow_STVF_ety,
event.MET_CellOut_Eflow_STVF_sumet)
MET = self.get_met()
if self.verbose:
log.info("Run: {0} Event: {1}".format(event.RunNumber,
event.EventNumber))
log.info("Recalculated MET: %.3f (original: %.3f)" %
(MET.et(), event.MET_RefFinal_STVF_et))
log.info("Recalculated MET phi: %.3f (original: %.3f)" %
(MET.phi(), event.MET_RefFinal_STVF_phi))
if (abs(MET.et() - event.MET_RefFinal_STVF_et) /
event.MET_RefFinal_STVF_et) > 0.1:
log.warning("Large MET difference!")
# update the MET with the shifted value
self.tree.MET_etx_original = event.MET_RefFinal_STVF_etx
self.tree.MET_ety_original = event.MET_RefFinal_STVF_ety
self.tree.MET_et_original = event.MET_RefFinal_STVF_et
self.tree.MET_sumet_original = event.MET_RefFinal_STVF_sumet
self.tree.MET_phi_original = event.MET_RefFinal_STVF_phi
event.MET_RefFinal_STVF_etx = MET.etx()
event.MET_RefFinal_STVF_ety = MET.ety()
event.MET_RefFinal_STVF_et = MET.et()
event.MET_RefFinal_STVF_sumet = MET.sumet()
event.MET_RefFinal_STVF_phi = MET.phi()
return True