MET Constructor

class METConstructor

Inherits from xAH::Algorithm

Public Functions

METConstructor()
EL::StatusCode setupJob(EL::Job &job)
EL::StatusCode fileExecute()
EL::StatusCode histInitialize()
EL::StatusCode changeInput(bool firstFile)
EL::StatusCode initialize()
EL::StatusCode execute()
EL::StatusCode postExecute()
EL::StatusCode finalize()
EL::StatusCode histFinalize()
StatusCode algInitialize()

Run any initializations commmon to all xAH Algorithms (such as registerInstance). Call this inside histInitialize for best results.

StatusCode algFinalize()

Run any finalizations common to all xAH Algorithms (such as unregisterInstance). Call this inside histFinalize for best results.

StatusCode parseSystValVector()

Parse string of systematic sigma levels in m_systValVectorString into m_systValVector.

Public Members

TString m_referenceMETContainer = "MET_Reference_AntiKt4LCTopo"
TString m_mapName = "METAssoc_AntiKt4LCTopo"
TString m_coreName = "MET_Core_AntiKt4LCTopo"
TString m_outputContainer = "NewRefFinal"
TString m_inputJets = ""
TString m_inputElectrons = ""
TString m_inputPhotons = ""
TString m_inputTaus = ""
TString m_inputMuons = ""
bool m_doElectronCuts = false
bool m_doPhotonCuts = false
bool m_doTauCuts = false
bool m_doMuonCuts = false
bool m_doMuonEloss = false
bool m_doIsolMuonEloss = false
bool m_doJVTCut = false
bool m_dofJVTCut = false
bool m_rebuildUsingTracksInJets = false

Rebuild MET using tracks in calo jets.

bool m_addSoftClusterTerms = false

Include soft cluster terms if rebuilding MET using jet terms (only considered if m_rebuildUsingTracksInJets is false)

bool m_setAFII = false

Force AFII flag in calculation, in case metadata is broken.

bool m_calculateSignificance = false

Enable MET significance calculation.

bool m_significanceTreatPUJets = true

Introduce “resolution” for jets with low JVT, if the analysis is sensitive to pileup jets.

double m_significanceSoftTermReso = 10.0

Set soft term resolution.

bool m_runNominal = true

set to false if you want to run met systematics

std::string m_systName = "All"

do not change it, not useful

float m_systVal = 1.0
bool m_writeSystToMetadata = false

Write systematics names to metadata.

std::string m_SoftTermSystConfigFile = "TrackSoftTerms.config"
std::string m_jetSystematics = ""

Name of jet systematics vector from JetCalibrator.

std::string m_eleSystematics = ""

Name of electron systematics vector from ElectronCalibrator.

std::string m_muonSystematics = ""

Name of muon systematics vector from MuonCalibrator.

std::string m_phoSystematics = ""

Name of photon systematics vector from PhotonCalibrator.

std::string m_outputAlgoSystNames = ""
std::string m_name = "UnnamedAlgorithm"

All algorithms initialized should have a unique name, to differentiate them at the TObject level.

Note, GetName() returns a char* while this returns a std::string.

bool m_debug = false

m_debug is being deprecated

bool m_verbose = false

m_verbose is being deprecated

MSG::Level m_msgLevel = MSG::INFO

debug level

std::string m_systValVectorString = ""

If running systematics, you can run multiple points and store them in here. A comma separated list of working points should be given to m_systValVectorString, and then parsed by calling parseSystValVector.

std::vector<float> m_systValVector
std::string m_eventInfoContainerName = "EventInfo"

If the xAOD has a different EventInfo container name, set it here

int m_isMC = -1

This is an override at the algorithm level to force analyzing MC or not.

Value Meaning
-1 Default, use eventInfo object to determine if data or mc
0 Treat the input as data
1 Treat the input as MC

Protected Functions

bool isMC()

Try to determine if we are running over data or MC. The xAH::Algorithm::m_isMC can be used

to fix the return value. Otherwise the EventInfo object is queried.

An exception is thrown if the type cannot be determined.

Return Value Meaning
0 Data
1 MC

void registerInstance()

Register the given instance under the moniker xAH::Algorithm::m_className

This will increase the reference count by 1.

int numInstances()

Return number of instances registered under the moniker xAH::Algorithm::m_className

This will return the reference count.

Warning

If for some reason the instance wasn’t registered, we spit out a warning.

void unregisterInstance()

Unregister the given instance under the moniker xAH::Algorithm::m_className

This will decrease the reference count by 1.

Warning

If for some reason the instance wasn’t registered, we spit out a warning.

template <typename T>
StatusCode checkToolStore(const std::string &tool_name)

Check whether the input CP tool already exists with this name in the asg::ToolStore

Depending on the outcome, the content of the map xAH::Algorithm::m_toolAlreadyUsed wll be set accordingly.

bool isToolAlreadyUsed(const std::string &tool_name)

Check whether the input CP tool has been already used by any xAH::Algorithm in the current job by scanning xAH::Algorithm::m_toolAlreadyUsed.

template <typename T>
void setToolName(asg::AnaToolHandle<T> &handle, std::string name = "") const

Note: This function does nothing in release 21.1! The native private tool mechanism is used instead.

Sets the name of a tool. If no name is needed, the tool will use the name of the algorithm plus a unique identifier (xAH::Algorithm::getAddress()) appended to ensure the tool is unique and effectively private.

The tool will not be guaranteed unique if two tools of the same type are created without a name passed in. But this is, at this point, up to the user and a more complex scenario than what this function tries to simplify on its own.

std::string getAddress() const

Return a std::string representation of this

Protected Attributes

std::string m_className = "Algorithm"

The moniker by which all instances are tracked in xAH::Algorithm::m_instanceRegistry

xAOD::TEvent *m_event = nullptr

The TEvent object

xAOD::TStore *m_store = nullptr

The TStore object

Private Members

asg::AnaToolHandle<IMETMaker> m_metmaker_handle
asg::AnaToolHandle<IMETSystematicsTool> m_metSyst_handle
asg::AnaToolHandle<IMETSignificance> m_metSignificance_handle
TauAnalysisTools::TauSelectionTool *m_tauSelTool
TString coreMetKey
std::vector<CP::SystematicSet> sysList
int m_numEvent