IDPerfMonZmumu Class Reference

#include <IDPerfMonZmumu.h>

Inheritance diagram for IDPerfMonZmumu:

Public Member Functions
	IDPerfMonZmumu (const std::string &name, ISvcLocator *pSvcLocator)
	needed for IP resolution studies
	~IDPerfMonZmumu ()
virtual StatusCode	initialize ()
virtual StatusCode	execute (const EventContext &ctx)
	Execute method.
virtual StatusCode	finalize ()
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual bool	isClonable () const override
	Specify if the algorithm is clonable.
virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
virtual bool	filterPassed (const EventContext &ctx) const
	Get filter decision:
virtual void	setFilterPassed (bool state, const EventContext &ctx) const
	Set filter decision:
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
virtual bool	isReEntrant () const override final
	Legacy algorithms are not thread-safe.
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
StatusCode	bookTrees ()
StatusCode	CheckTriggerStatusAndPrescale ()
void	Clear4MuNtupleVariables ()
void	ExtractIDHitsInformation (const xAOD::Muon *muon_pos, const xAOD::Muon *muon_neg)
int	GetMuonQualityValue (std::string qualityname)
void	RegisterHistograms ()
void	ResetCommonNtupleVectors ()
const xAOD::Vertex *	GetDiMuonVertex (const xAOD::TrackParticle *, const xAOD::TrackParticle *)
StatusCode	FillRecParameters (const Trk::Track *track, const xAOD::TrackParticle *trackp_for_unbias, double charge, const xAOD::Vertex *vertex, const EventContext &ctx)
StatusCode	FillRecParametersSimple (const Trk::Track *track, float charge, const xAOD::Vertex *vertex, const EventContext &ctx)
StatusCode	FillRecParametersTP (const xAOD::TrackParticle *trackp, const xAOD::TrackParticle *trackp_for_unbias, double charge, const xAOD::Vertex *vertex, const EventContext &ctx)
StatusCode	FillTruthParameters (const xAOD::TrackParticle *track, const EventContext &ctx)
const xAOD::TruthParticle *	getTruthParticle (const xAOD::IParticle &p)
StatusCode	RunFourLeptonAnalysis ()
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ZmumuEvent	m_xZmm
FourMuonEvent	m_4mu
bool	m_UseTrigger {true}
bool	m_doIsoSelection {true}
bool	m_doIPSelection {true}
bool	m_doMCPSelection {true}
double	m_MassWindowLow {60.}
double	m_MassWindowHigh {120.}
double	m_LeadingMuonPtCut {20.}
double	m_SecondMuonPtCut {15.}
double	m_OpeningAngleCut {0.2}
double	m_Z0GapCut {5.}
bool	m_isMC {}
bool	m_doRefit {}
bool	m_useTrackSelectionTool {}
bool	m_doIP {}
bool	m_doFourMuAnalysis {}
bool	m_storeZmumuNtuple {true}
bool	m_skipMS {}
bool	m_useCustomMuonSelector {}
int	m_minGoodLumiBlock {}
int	m_maxGoodLumiBlock {}
ToolHandle< IegammaTrkRefitterTool >	m_TrackRefitter1
	The track refitter.
ToolHandle< IegammaTrkRefitterTool >	m_TrackRefitter2
	The track refitter.
ToolHandle< Reco::ITrackToVertex >	m_trackToVertexTool
	tool to extrapolate tracks to BL
ToolHandle< Trig::TrigDecisionTool >	m_triggerDecision
ToolHandle< Trig::IMatchingTool >	m_triggerMatching
ToolHandle< InDet::IInDetTrackSelectionTool >	m_selTool
	The track selection Tool.
ToolHandle< Trk::ITrackToVertexIPEstimator >	m_trackToVertexIPEstimator
	Needed for IP resolution studies.
SG::ReadCondHandleKey< InDet::BeamSpotData >	m_beamSpotKey { this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot" }
	used for truth parameters
ToolHandle< Trk::IExtrapolator >	m_extrapolator
ToolHandle< CP::IMuonSelectionTool >	m_muonSelector
	used to pass a custom muon selector
bool	m_validationMode
	< boolean to switch to validation mode
std::string	m_commonTreeName
std::string	m_defaultTreeName
std::string	m_IDTreeName
std::string	m_refit1TreeName
std::string	m_refit2TreeName
std::string	m_truthTreeName
std::string	m_combTreeName
std::string	m_MSTreeName
std::string	m_FourMuTreeName
std::string	m_ValidationTreeDescription
	< validation tree description - second argument in TTree
std::string	m_commonTreeFolder
std::string	m_defaultTreeFolder
std::string	m_IDTreeFolder
std::string	m_refit1TreeFolder
std::string	m_refit2TreeFolder
std::string	m_truthTreeFolder
std::string	m_combTreeFolder
std::string	m_MSTreeFolder
std::string	m_FourMuTreeFolder
	Root Validation Tree.
TTree *	m_commonTree
TTree *	m_defaultTree
TTree *	m_IDTree
TTree *	m_refit1Tree
TTree *	m_refit2Tree
TTree *	m_truthTree
TTree *	m_combTree
TTree *	m_MSTree
TTree *	m_FourMuTree
bool	m_doRemoval {}
bool	m_doDebug {}
ToolHandle< CP::ITrackVertexAssociationTool >	m_Trk2VtxAssociationTool
SG::ReadHandleKey< xAOD::VertexContainer >	m_vertexKey { this, "VertexContainer", "PrimaryVertices", "primary vertex container" }
SG::ReadHandleKey< xAOD::TrackParticleContainer >	m_trackContainerName {this, "trackContainerName", "InDetTrackParticles"}
	IDTtacks.
std::string	m_truthName
std::string	m_trackParticleName
	Track(Particle)TruthCollection input name.
std::string	m_truthLinkVecName
	TrackParticle input name.
TH1F *	m_h_cutflow {}
	link vector to map HepMC onto xAOD truth
unsigned int	m_runNumber {}
unsigned int	m_evtNumber {}
unsigned int	m_lumi_block {}
unsigned int	m_event_mu {}
float	m_event_weight {}
int	m_triggerPrescale {}
std::string	m_triggerName
unsigned int	m_nVertex {}
float	m_beamposX {}
float	m_beamposY {}
double	m_positive_px {}
double	m_positive_py {}
double	m_positive_pt {}
double	m_positive_pz {}
double	m_positive_phi {}
double	m_positive_eta {}
double	m_positive_z0 {}
double	m_positive_d0 {}
double	m_positive_z0_err {}
double	m_positive_d0_err {}
double	m_positive_sigma_pt {}
double	m_positive_z0_PV {}
double	m_positive_d0_PV {}
double	m_positive_z0_PVerr {}
double	m_positive_d0_PVerr {}
double	m_positive_qoverp {}
double	m_positive_sigma_qoverp {}
int	m_positive_1_vtx {}
int	m_positive_parent {}
double	m_positive_2_px {}
double	m_positive_2_py {}
double	m_positive_2_pz {}
double	m_positive_2_z0 {}
double	m_positive_2_d0 {}
double	m_positive_2_z0_err {}
double	m_positive_2_d0_err {}
double	m_positive_2_z0_PV {}
double	m_positive_2_d0_PV {}
double	m_positive_2_z0_PVerr {}
double	m_positive_2_d0_PVerr {}
int	m_positive_2_vtx {}
double	m_negative_px {}
double	m_negative_py {}
double	m_negative_pt {}
double	m_negative_pz {}
double	m_negative_phi {}
double	m_negative_eta {}
double	m_negative_z0 {}
double	m_negative_d0 {}
double	m_negative_z0_err {}
double	m_negative_d0_err {}
double	m_negative_sigma_pt {}
double	m_negative_z0_PV {}
double	m_negative_d0_PV {}
double	m_negative_z0_PVerr {}
double	m_negative_d0_PVerr {}
double	m_negative_qoverp {}
double	m_negative_sigma_qoverp {}
int	m_negative_1_vtx {}
int	m_negative_parent {}
double	m_negative_2_px {}
double	m_negative_2_py {}
double	m_negative_2_pz {}
double	m_negative_2_z0 {}
double	m_negative_2_d0 {}
double	m_negative_2_z0_err {}
double	m_negative_2_d0_err {}
double	m_negative_2_z0_PV {}
double	m_negative_2_d0_PV {}
double	m_negative_2_z0_PVerr {}
double	m_negative_2_d0_PVerr {}
int	m_negative_2_vtx {}
double	m_el_negative1_px {}
double	m_el_negative1_py {}
double	m_el_negative1_pz {}
double	m_el_negative1_z0 {}
double	m_el_negative1_d0 {}
double	m_el_negative1_z0_err {}
double	m_el_negative1_d0_err {}
double	m_el_negative1_z0_PV {}
double	m_el_negative1_d0_PV {}
double	m_el_negative1_z0_PVerr {}
double	m_el_negative1_d0_PVerr {}
int	m_el_negative1_vtx {}
double	m_el_negative2_px {}
double	m_el_negative2_py {}
double	m_el_negative2_pz {}
double	m_el_negative2_z0 {}
double	m_el_negative2_d0 {}
double	m_el_negative2_z0_err {}
double	m_el_negative2_d0_err {}
double	m_el_negative2_z0_PV {}
double	m_el_negative2_d0_PV {}
double	m_el_negative2_z0_PVerr {}
double	m_el_negative2_d0_PVerr {}
int	m_el_negative2_vtx {}
double	m_el_positive1_px {}
double	m_el_positive1_py {}
double	m_el_positive1_pz {}
double	m_el_positive1_z0 {}
double	m_el_positive1_d0 {}
double	m_el_positive1_z0_err {}
double	m_el_positive1_d0_err {}
double	m_el_positive1_z0_PV {}
double	m_el_positive1_d0_PV {}
double	m_el_positive1_z0_PVerr {}
double	m_el_psoitive1_d0_PVerr {}
int	m_el_positive1_vtx {}
double	m_el_positive2_px {}
double	m_el_positive2_py {}
double	m_el_positive2_pz {}
double	m_el_positive2_z0 {}
double	m_el_positive2_d0 {}
double	m_el_positive2_z0_err {}
double	m_el_positive2_d0_err {}
double	m_el_positive2_z0_PV {}
double	m_el_positive2_d0_PV {}
double	m_el_positive2_z0_PVerr {}
double	m_el_psoitive2_d0_PVerr {}
int	m_el_positive2_vtx {}
double	m_4mu_minv {}
double	m_pv_x {}
double	m_pv_y {}
double	m_pv_z {}
unsigned int	m_nTrkInVtx {}
double	m_met {}
double	m_metphi {}
std::vector< float >	m_IDTrack_pt
std::vector< float >	m_CBTrack_pt
std::vector< float >	m_Refit1_pt
std::vector< float >	m_Refit2_pt
std::vector< float >	m_Truth_pt
std::vector< float >	m_IDTrack_eta
std::vector< float >	m_CBTrack_eta
std::vector< float >	m_Refit1_eta
std::vector< float >	m_Refit2_eta
std::vector< float >	m_Truth_eta
std::vector< float >	m_IDTrack_phi
std::vector< float >	m_CBTrack_phi
std::vector< float >	m_Refit1_phi
std::vector< float >	m_Refit2_phi
std::vector< float >	m_Truth_phi
std::vector< float >	m_IDTrack_d0
std::vector< float >	m_CBTrack_d0
std::vector< float >	m_Refit1_d0
std::vector< float >	m_Refit2_d0
std::vector< float >	m_Truth_d0
std::vector< float >	m_IDTrack_z0
std::vector< float >	m_CBTrack_z0
std::vector< float >	m_Refit1_z0
std::vector< float >	m_Refit2_z0
std::vector< float >	m_Truth_z0
std::vector< float >	m_IDTrack_qoverp
std::vector< float >	m_CBTrack_qoverp
std::vector< float >	m_Refit1_qoverp
std::vector< float >	m_Refit2_qoverp
std::vector< float >	m_Truth_qoverp
std::vector< int >	m_Truth_parent
std::vector< float >	m_IDTrack_sigma_pt
std::vector< float >	m_CBTrack_sigma_pt
std::vector< float >	m_Refit1_sigma_pt
std::vector< float >	m_Refit2_sigma_pt
std::vector< float >	m_IDTrack_sigma_d0
std::vector< float >	m_CBTrack_sigma_d0
std::vector< float >	m_Refit1_sigma_d0
std::vector< float >	m_Refit2_sigma_d0
std::vector< float >	m_IDTrack_sigma_z0
std::vector< float >	m_CBTrack_sigma_z0
std::vector< float >	m_Refit1_sigma_z0
std::vector< float >	m_Refit2_sigma_z0
std::vector< float >	m_IDTrack_sigma_qoverp
std::vector< float >	m_CBTrack_sigma_qoverp
std::vector< float >	m_Refit1_sigma_qoverp
std::vector< float >	m_Refit2_sigma_qoverp
std::vector< int >	m_nBLhits
std::vector< int >	m_nPIXhits
std::vector< int >	m_nSCThits
std::vector< int >	m_nTRThits
std::string	m_sTriggerChainName
std::string	m_outputTracksName
std::string	m_MuonQualityName
SG::ReadHandleKey< xAOD::EventInfo >	m_EventInfoKey {this, "EventInfoKey", "EventInfo"}
DataObjIDColl	m_extendedExtraObjects
	Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default).
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default).
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 73 of file IDPerfMonZmumu.h.

Member Typedef Documentation

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

IDPerfMonZmumu()

IDPerfMonZmumu::IDPerfMonZmumu	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

needed for IP resolution studies

missing ET

Definition at line 47 of file IDPerfMonZmumu.cxx.

                           :
  AthAlgorithm(name, pSvcLocator),
  m_isMC(false),
  m_doRefit(false),
  m_TrackRefitter1(""),
  m_TrackRefitter2(""),
  m_trackToVertexTool("Reco::TrackToVertex"),
  m_triggerDecision("Trig::TrigDecisionTool/TrigDecisionTool"),
  m_triggerMatching("Trig::MatchingTool/MatchingTool"),
  m_selTool( "InDet::InDetTrackSelectionTool/TrackSelectionTool"),
  m_trackToVertexIPEstimator("Trk::TrackToVertexIPEstimator"),
  m_extrapolator("Trk::Extrapolator/AtlasExtrapolator"),
  m_muonSelector("CP::MuonSelectionTool/MuonSelectionTool"),
  m_validationMode(true),
 
  m_commonTreeName ("commonTree"),
  m_defaultTreeName("Default_Particle"),
  m_IDTreeName     ("ID_InDetTrackParticle"),
  m_refit1TreeName ("Refit1_SiAndTRT"),
  m_refit2TreeName ("Refit2_SiOnly"),
  m_truthTreeName  ("TruthParams"),
  m_combTreeName   ("CombinedTrackParticle"),
  m_MSTreeName     ("MS_TrackParticle"),
  m_FourMuTreeName ("FourMu"),
 
  m_ValidationTreeDescription("Small Tree for Zmumu fits"),
 
  m_commonTreeFolder ("/ZmumuValidation/common"),
  m_defaultTreeFolder("/ZmumuValidation/default"),
  m_IDTreeFolder     ("/ZmumuValidation/ID"),
  m_refit1TreeFolder ("/ZmumuValidation/refit1"),
  m_refit2TreeFolder ("/ZmumuValidation/refit2"),
  m_truthTreeFolder  ("/ZmumuValidation/truth"),
  m_combTreeFolder   ("/ZmumuValidation/comb"),
  m_MSTreeFolder     ("/ZmumuValidation/ms"),
  m_FourMuTreeFolder ("/ZmumuValidation/fourmu"),
 
  m_commonTree (nullptr),
  m_defaultTree(nullptr),
  m_IDTree(nullptr),
  m_refit1Tree(nullptr),
  m_refit2Tree(nullptr),
  m_truthTree(nullptr),
  m_combTree(nullptr),
  m_MSTree(nullptr),
  m_FourMuTree(nullptr),
  m_doRemoval(true),
  m_doDebug(false),
  m_Trk2VtxAssociationTool("CP::TrackVertexAssociationTool", this)
{
  // Properties that are set from the python scripts.
 
  declareProperty("UseTrigger"    ,    m_UseTrigger     = true);
  declareProperty("doIsoSelection",    m_doIsoSelection = true );
  declareProperty("doIPSelection",     m_doIPSelection = true );
  declareProperty("doMCPSelection",    m_doMCPSelection = true );
  declareProperty("isMC",              m_isMC = false);
  declareProperty("doRefit",           m_doRefit = false);
  declareProperty("doIPextrToPV",      m_doIP = false);
  declareProperty("Extrapolator",      m_extrapolator );
  declareProperty("MassWindowLow",     m_MassWindowLow = 60.0, "Lower cut in mu+mu- invariant mass" );
  declareProperty("MassWindowHigh",    m_MassWindowHigh = 120.0, "Upper cut in mu+mu- invariant mass" );
  declareProperty("OpeningAngle",      m_OpeningAngleCut = 0.2, "Opening angle between the two muons (in radians)");
  declareProperty("OutputTracksName",  m_outputTracksName = "ZmumuTracks");
  declareProperty("PtLeadingMuon",     m_LeadingMuonPtCut = 20., "Pt cut on the leading muon");
  declareProperty("PtSecondMuon",      m_SecondMuonPtCut = 15., "Pt cut on the second muon");
  declareProperty("ReFitterTool1",     m_TrackRefitter1, "ToolHandle for track fitter implementation");
  declareProperty("ReFitterTool2",     m_TrackRefitter2, "ToolHandle for track fitter implementation");
  declareProperty("TrackToVertexTool", m_trackToVertexTool);
  declareProperty("TrackTruthName",    m_truthName="TrackTruthCollection");
  declareProperty("TrackParticleName", m_trackParticleName="CombinedTrackParticle");
  declareProperty("triggerChainName",  m_sTriggerChainName);
  declareProperty("ValidationMode",    m_validationMode);
  declareProperty("xAODTruthLinkVector",      m_truthLinkVecName="xAODTruthLinks");
  declareProperty("Z0Gap",                    m_Z0GapCut = 5.0, "maximum gap between the z0 of both muons (in mm)");
  declareProperty("TrackSelectionTool",       m_selTool );
  declareProperty("UseTrackSelectionTool",    m_useTrackSelectionTool = false);
  declareProperty("TrackToVertexIPEstimator", m_trackToVertexIPEstimator);
 
 
  declareProperty("commonTreeFolder",  m_commonTreeFolder, "/ZmumuValidationUserSel/common" );
  declareProperty("defaultTreeFolder", m_defaultTreeFolder );
  declareProperty("IDTreeFolder",      m_IDTreeFolder );
  declareProperty("refit1TreeFolder",  m_refit1TreeFolder );
  declareProperty("refit2TreeFolder",  m_refit2TreeFolder );
  declareProperty("truthTreeFolder",   m_truthTreeFolder );
  declareProperty("combTreeFolder",    m_combTreeFolder );
  declareProperty("MSTreeFolder",      m_MSTreeFolder, "/ZmumuValidationUserSel/ms" );
 
  declareProperty("UnbiasVertex",      m_doRemoval);
  
  declareProperty("commonTree",  m_commonTreeName, "CommonTree" );
  declareProperty("DefaultTree", m_defaultTreeName );
  declareProperty("IDTree",      m_IDTreeName );
  declareProperty("CBTree",      m_combTreeName );
  declareProperty("Refit1Tree",  m_refit1TreeName );
  declareProperty("Refit2Tree",  m_refit2TreeName );
  declareProperty("MSTree",      m_MSTreeName );
  
  declareProperty("doFourMuAnalysis", m_doFourMuAnalysis = false);
  declareProperty("FourMuTreeFolder", m_FourMuTreeFolder);
 
  declareProperty("StoreZmumuNtuple", m_storeZmumuNtuple = true);
  declareProperty("doZmumuEventDebug", m_doDebug);
  declareProperty("MuonQuality", m_MuonQualityName = "Medium");
  declareProperty("skipMS", m_skipMS = false);
  declareProperty("useCustomMuonSelector", m_useCustomMuonSelector = false );
  declareProperty("MuonSelector", m_muonSelector );
 
  declareProperty( "MinLumiBlock", m_minGoodLumiBlock = 0, "minimum lumiblock number to be accepted");
  declareProperty( "MaxLumiBlock", m_maxGoodLumiBlock = 0, "maximum lumiblock number to be accepted");
 
  return;
}

~IDPerfMonZmumu()

IDPerfMonZmumu::~IDPerfMonZmumu

(

)

Definition at line 165 of file IDPerfMonZmumu.cxx.

{}

Member Function Documentation

bookTrees()

StatusCode IDPerfMonZmumu::bookTrees ( )

private

Definition at line 301 of file IDPerfMonZmumu.cxx.

{
  m_h_cutflow = new TH1F("h_cutflow","cut flow histogram",11, -0.5, 9.5);
 
  ATH_MSG_DEBUG("initialize() ** bookTrees() ** m_commonTree name: " << m_commonTreeName.c_str());    
  ATH_MSG_DEBUG("                              m_defaultTree name: " << m_defaultTreeName.c_str());    
  ATH_MSG_DEBUG("                                   m_IDTree name: " << m_IDTreeName.c_str());    
 
  if (!m_commonTree) {
    ATH_MSG_INFO("initialize() ** defining m_commonTree with name: " << m_commonTreeName.c_str());    
    m_commonTree = new TTree((m_commonTreeName).c_str(), m_ValidationTreeDescription.c_str());
 
    m_commonTree->Branch("runNumber"           , &m_runNumber,  "runNumber/I");
    m_commonTree->Branch("eventNumber"         , &m_evtNumber,  "eventNumber/I");
    m_commonTree->Branch("lumi_block"          , &m_lumi_block, "lumi_block/I");
    m_commonTree->Branch("beamposX"            , &m_beamposX,   "beamposX/F");
    m_commonTree->Branch("beamposY"            , &m_beamposY,   "beamposY/F");
    m_commonTree->Branch("mu"                  , &m_event_mu,   "mu/I");
    m_commonTree->Branch("preScale"            , &m_triggerPrescale, "preScale/I");
    m_commonTree->Branch("mcEventWeight"       , &m_event_weight, "mcEventWeight/F");
    m_commonTree->Branch("IDTrack_pt"          , &m_IDTrack_pt); 
    m_commonTree->Branch("IDTrack_eta"         , &m_IDTrack_eta); 
    m_commonTree->Branch("IDTrack_phi"         , &m_IDTrack_phi);
    m_commonTree->Branch("IDTrack_d0"          , &m_IDTrack_d0);
    m_commonTree->Branch("IDTrack_z0"          , &m_IDTrack_z0);
    m_commonTree->Branch("IDTrack_qoverp"      , &m_IDTrack_qoverp);
    m_commonTree->Branch("IDTrack_sigma_pt"    , &m_IDTrack_sigma_pt);
    m_commonTree->Branch("IDTrack_sigma_d0"    , &m_IDTrack_sigma_d0);
    m_commonTree->Branch("IDTrack_sigma_z0"    , &m_IDTrack_sigma_z0);
    m_commonTree->Branch("IDTrack_sigma_qoverp", &m_IDTrack_sigma_qoverp);
    m_commonTree->Branch("CBTrack_pt"          , &m_CBTrack_pt); 
    m_commonTree->Branch("CBTrack_eta"         , &m_CBTrack_eta); 
    m_commonTree->Branch("CBTrack_phi"         , &m_CBTrack_phi); 
    m_commonTree->Branch("CBTrack_d0"          , &m_CBTrack_d0); 
    m_commonTree->Branch("CBTrack_z0"          , &m_CBTrack_z0); 
    m_commonTree->Branch("CBTrack_qoverp"      , &m_CBTrack_qoverp);
    m_commonTree->Branch("CBTrack_sigma_pt"    , &m_CBTrack_sigma_pt);
    m_commonTree->Branch("CBTrack_sigma_d0"    , &m_CBTrack_sigma_d0);
    m_commonTree->Branch("CBTrack_sigma_z0"    , &m_CBTrack_sigma_z0);
    m_commonTree->Branch("CBTrack_sigma_qoverp", &m_CBTrack_sigma_qoverp);
    m_commonTree->Branch("Refit1_pt"           , &m_Refit1_pt); 
    m_commonTree->Branch("Refit1_eta"          , &m_Refit1_eta); 
    m_commonTree->Branch("Refit1_phi"          , &m_Refit1_phi); 
    m_commonTree->Branch("Refit1_d0"           , &m_Refit1_d0); 
    m_commonTree->Branch("Refit1_z0"           , &m_Refit1_z0); 
    m_commonTree->Branch("Refit1_qoverp"       , &m_Refit1_qoverp); 
    m_commonTree->Branch("Refit1_sigma_pt"     , &m_Refit1_sigma_pt);
    m_commonTree->Branch("Refit1_sigma_d0"     , &m_Refit1_sigma_d0);
    m_commonTree->Branch("Refit1_sigma_z0"     , &m_Refit1_sigma_z0);
    m_commonTree->Branch("Refit1_sigma_qoverp" , &m_Refit1_sigma_qoverp);
    m_commonTree->Branch("Refit2_pt"           , &m_Refit2_pt); 
    m_commonTree->Branch("Refit2_eta"          , &m_Refit2_eta); 
    m_commonTree->Branch("Refit2_phi"          , &m_Refit2_phi); 
    m_commonTree->Branch("Refit2_d0"           , &m_Refit2_d0); 
    m_commonTree->Branch("Refit2_z0"           , &m_Refit2_z0); 
    m_commonTree->Branch("Refit2_qoverp"       , &m_Refit2_qoverp); 
    m_commonTree->Branch("Refit2_sigma_pt"     , &m_Refit2_sigma_pt);
    m_commonTree->Branch("Refit2_sigma_d0"     , &m_Refit2_sigma_d0);
    m_commonTree->Branch("Refit2_sigma_z0"     , &m_Refit2_sigma_z0);
    m_commonTree->Branch("Refit2_sigma_qoverp" , &m_Refit2_sigma_qoverp);
    m_commonTree->Branch("Truth_pt"            , &m_Truth_pt); 
    m_commonTree->Branch("Truth_eta"           , &m_Truth_eta); 
    m_commonTree->Branch("Truth_phi"           , &m_Truth_phi); 
    m_commonTree->Branch("Truth_d0"            , &m_Truth_d0); 
    m_commonTree->Branch("Truth_z0"            , &m_Truth_z0); 
    m_commonTree->Branch("Truth_qoverp"        , &m_Truth_qoverp); 
    m_commonTree->Branch("Truth_parent"        , &m_Truth_parent);
    m_commonTree->Branch("numberOfBLayerHits"  , &m_nBLhits);
    m_commonTree->Branch("numberOfPixelHits"   , &m_nPIXhits);
    m_commonTree->Branch("numberOfSCTHits"     , &m_nSCThits);
    m_commonTree->Branch("numberOfTRTHits"     , &m_nTRThits);
  }
 
  if ( !m_defaultTree){
    ATH_MSG_INFO("initialize() ** defining m_defaultTree with name: " << m_defaultTreeName.c_str());    
    m_defaultTree = new TTree((m_defaultTreeName).c_str(), m_ValidationTreeDescription.c_str());
 
    m_defaultTree->Branch("runNumber"      ,  &m_runNumber,  "runNumber/I");
    m_defaultTree->Branch("eventNumber"    ,  &m_evtNumber,  "eventNumber/I");
    m_defaultTree->Branch("lumi_block"     ,  &m_lumi_block, "lumi_block/I");
    m_defaultTree->Branch("mu"             ,  &m_event_mu,   "mu/I");
 
    m_defaultTree->Branch("Negative_Px",  &m_negative_px,  "Negative_Px/D");
    m_defaultTree->Branch("Negative_Py",  &m_negative_py,  "Negative_Py/D");
    m_defaultTree->Branch("Negative_Pt",  &m_negative_pt,  "Negative_Pt/D");
    m_defaultTree->Branch("Negative_Pz",  &m_negative_pz,  "Negative_Pz/D");
    m_defaultTree->Branch("Negative_Phi", &m_negative_phi, "Negative_Phi/D");
    m_defaultTree->Branch("Negative_eta", &m_negative_eta, "Negative_Eta/D");
    m_defaultTree->Branch("Negative_z0",  &m_negative_z0,  "Negative_z0/D");
    m_defaultTree->Branch("Negative_d0",  &m_negative_d0,  "Negative_d0/D");
    m_defaultTree->Branch("Negative_z0_err",  &m_negative_z0_err,  "Negative_z0_err/D");
    m_defaultTree->Branch("Negative_d0_err",  &m_negative_d0_err,  "Negative_d0_err/D");
    m_defaultTree->Branch("Negative_sigma_pt",  &m_negative_sigma_pt,  "Negative_sigma_pt/D");
 
    m_defaultTree->Branch("Positive_Px",  &m_positive_px,  "Positive_Px/D");
    m_defaultTree->Branch("Positive_Py",  &m_positive_py,  "Positive_Py/D");
    m_defaultTree->Branch("Positive_Pt",  &m_positive_pt,  "Positive_Pt/D");
    m_defaultTree->Branch("Positive_Pz",  &m_positive_pz,  "Positive_Pz/D");
    m_defaultTree->Branch("Positive_Phi", &m_positive_phi, "Positive_Phi/D");
    m_defaultTree->Branch("Positive_eta", &m_positive_eta, "Positive_Eta/D");
    m_defaultTree->Branch("Positive_z0",  &m_positive_z0,  "Positive_z0/D");
    m_defaultTree->Branch("Positive_d0",  &m_positive_d0,  "Positive_d0/D");
    m_defaultTree->Branch("Positive_z0_err",  &m_positive_z0_err,  "Positive_z0_err/D");
    m_defaultTree->Branch("Positive_d0_err",  &m_positive_d0_err,  "Positive_d0_err/D");
    m_defaultTree->Branch("Positive_sigma_pt",  &m_positive_sigma_pt,  "Positive_sigma_pt/D");
    
    if (m_doIP) {
      m_defaultTree->Branch("Negative_d0_PV",      &m_negative_d0_PV   ,  "Negative_d0_PV/D");
      m_defaultTree->Branch("Negative_z0_PV",      &m_negative_z0_PV   ,  "Negative_z0_PV/D");
      m_defaultTree->Branch("Positive_z0_PV",      &m_positive_z0_PV,  "Positive_z0_PV/D");
      m_defaultTree->Branch("Positive_d0_PV",      &m_positive_d0_PV,  "Positive_d0_PV/D");
 
      m_defaultTree->Branch("Negative_d0_PVerr",   &m_negative_d0_PVerr,  "Negative_d0_PVerr/D");
      m_defaultTree->Branch("Negative_z0_PVerr",   &m_negative_z0_PVerr,  "Negative_z0_PVerr/D");
      m_defaultTree->Branch("Positive_z0_PVerr",   &m_positive_z0_PVerr,  "Positive_z0_PVerr/D");
      m_defaultTree->Branch("Positive_d0_PVerr",   &m_positive_d0_PVerr,  "Positive_d0_PVerr/D");
 
      m_defaultTree->Branch("pv_x",      &m_pv_x   ,  "pv_x/D");
      m_defaultTree->Branch("pv_y",      &m_pv_y   ,  "pv_y/D");
      m_defaultTree->Branch("pv_z",      &m_pv_z   ,  "pv_z/D");
      m_defaultTree->Branch("nTrkInVtx", &m_nTrkInVtx,  "nTrkInVtx/I");
    }
  }
  
  bool isTreeNone = false;
  if (m_IDTreeName.find("none") != std::string::npos) isTreeNone = true;
  if ( !m_IDTree and !isTreeNone ){
    ATH_MSG_INFO("initialize() ** defining IDPerfMonZmumu m_IDTree with name: " << m_IDTreeName.c_str());    
    m_IDTree = new TTree((m_IDTreeName).c_str(), m_ValidationTreeDescription.c_str());
    
    m_IDTree->Branch("runNumber"      ,  &m_runNumber,  "runNumber/I");
    m_IDTree->Branch("eventNumber"    ,  &m_evtNumber,  "eventNumber/I");
    m_IDTree->Branch("lumi_bLock"     ,  &m_lumi_block,  "lumi_block/I");
    m_IDTree->Branch("mu"             ,  &m_event_mu,  "mu/I");
 
    m_IDTree->Branch("Negative_Px",  &m_negative_px,  "Negative_Px/D");
    m_IDTree->Branch("Negative_Py",  &m_negative_py,  "Negative_Py/D");
    m_IDTree->Branch("Negative_Pt",  &m_negative_pt,  "Negative_Pt/D");
    m_IDTree->Branch("Negative_Pz",  &m_negative_pz,  "Negative_Pz/D");
    m_IDTree->Branch("Negative_Phi", &m_negative_phi, "Negative_Phi/D");
    m_IDTree->Branch("Negative_Eta", &m_negative_eta, "Negative_Eta/D");
    m_IDTree->Branch("Negative_z0",  &m_negative_z0,  "Negative_z0/D");
    m_IDTree->Branch("Negative_d0",  &m_negative_d0,  "Negative_d0/D");
    m_IDTree->Branch("Negative_z0_err",  &m_negative_z0_err,  "Negative_z0_err/D");
    m_IDTree->Branch("Negative_d0_err",  &m_negative_d0_err,  "Negative_d0_err/D");
    m_IDTree->Branch("Negative_sigma_pt",  &m_negative_sigma_pt,  "Negative_sigma_pt/D");
 
    m_IDTree->Branch("Positive_Px",  &m_positive_px,  "Positive_Px/D");
    m_IDTree->Branch("Positive_Py",  &m_positive_py,  "Positive_Py/D");
    m_IDTree->Branch("Positive_Pt",  &m_positive_pt,  "Positive_Pt/D");
    m_IDTree->Branch("Positive_Pz",  &m_positive_pz,  "Positive_Pz/D");
    m_IDTree->Branch("Positive_Phi", &m_positive_phi, "Positive_Phi/D");
    m_IDTree->Branch("Positive_Eta", &m_positive_eta, "Positive_Eta/D");
    m_IDTree->Branch("Positive_z0",  &m_positive_z0,  "Positive_z0/D");
    m_IDTree->Branch("Positive_d0",  &m_positive_d0,  "Positive_d0/D");
    m_IDTree->Branch("Positive_z0_err",  &m_positive_z0_err,  "Positive_z0_err/D");
    m_IDTree->Branch("Positive_d0_err",  &m_positive_d0_err,  "Positive_d0_err/D");
    m_IDTree->Branch("Positive_sigma_pt",  &m_positive_sigma_pt,  "Positive_sigma_pt/D");
    
    if(m_doIP){
      m_IDTree->Branch("Negative_d0_PV",      &m_negative_d0_PV   ,  "Negative_d0_PV/D");
      m_IDTree->Branch("Negative_z0_PV",      &m_negative_z0_PV   ,  "Negative_z0_PV/D");
      m_IDTree->Branch("Positive_z0_PV",      &m_positive_z0_PV,  "Positive_z0_PV/D");
      m_IDTree->Branch("Positive_d0_PV",      &m_positive_d0_PV,  "Positive_d0_PV/D");
      
      m_IDTree->Branch("Negative_d0_PVerr",   &m_negative_d0_PVerr,  "Negative_d0_PVerr/D");
      m_IDTree->Branch("Negative_z0_PVerr",   &m_negative_z0_PVerr,  "Negative_z0_PVerr/D");
      m_IDTree->Branch("Positive_z0_PVerr",   &m_positive_z0_PVerr,  "Positive_z0_PVerr/D");
      m_IDTree->Branch("Positive_d0_PVerr",   &m_positive_d0_PVerr,  "Positive_d0_PVerr/D");
 
      m_IDTree->Branch("pv_x",      &m_pv_x   ,  "pv_x/D");
      m_IDTree->Branch("pv_y",      &m_pv_y   ,  "pv_y/D");
      m_IDTree->Branch("pv_z",      &m_pv_z   ,  "pv_z/D");
      m_IDTree->Branch("nTrkInVtx", &m_nTrkInVtx,  "nTrkInVtx/I");
    }
  }
 
  // dealing with tree for Refit1 tracks
  isTreeNone = false;
  if (m_refit1TreeName.find("none") != std::string::npos) isTreeNone = true;
  if ( m_doRefit && !m_refit1Tree && !isTreeNone ){
    ATH_MSG_INFO("initialize() ** defining IDPerfMonZmumu m_refit1Tree with name: " << m_refit1TreeName.c_str());    
    m_refit1Tree = new TTree((m_refit1TreeName).c_str(), m_ValidationTreeDescription.c_str()); 
    
    m_refit1Tree->Branch("runNumber"      ,  &m_runNumber,  "runNumber/I");
    m_refit1Tree->Branch("eventNumber"    ,  &m_evtNumber,  "eventNumber/I");
    m_refit1Tree->Branch("lumi_block"     ,  &m_lumi_block,  "lumi_block/I");
    m_refit1Tree->Branch("mu"             ,  &m_event_mu,  "mu/I");
    m_refit1Tree->Branch("preScale"       ,  &m_triggerPrescale, "preScale/I");
 
    m_refit1Tree->Branch("Negative_Px",  &m_negative_px,  "Negative_Px/D");
    m_refit1Tree->Branch("Negative_Py",  &m_negative_py,  "Negative_Py/D");
    m_refit1Tree->Branch("Negative_Pt",  &m_negative_pt,  "Negative_Pt/D");
    m_refit1Tree->Branch("Negative_Pz",  &m_negative_pz,  "Negative_Pz/D");
    m_refit1Tree->Branch("Negative_Phi", &m_negative_phi, "Negative_Phi/D");
    m_refit1Tree->Branch("Negative_Eta", &m_negative_eta, "Negative_Eta/D");
    m_refit1Tree->Branch("Negative_z0",  &m_negative_z0,  "Negative_z0/D");
    m_refit1Tree->Branch("Negative_d0",  &m_negative_d0,  "Negative_d0/D");
    m_refit1Tree->Branch("Negative_z0_err",  &m_negative_z0_err,  "Negative_z0_err/D");
    m_refit1Tree->Branch("Negative_d0_err",  &m_negative_d0_err,  "Negative_d0_err/D");
    m_refit1Tree->Branch("Negative_sigma_pt",  &m_negative_sigma_pt,  "Negative_sigma_pt/D");
 
    m_refit1Tree->Branch("Positive_Px",  &m_positive_px,  "Positive_Px/D");
    m_refit1Tree->Branch("Positive_Py",  &m_positive_py,  "Positive_Py/D");
    m_refit1Tree->Branch("Positive_Pt",  &m_positive_pt,  "Positive_Pt/D");
    m_refit1Tree->Branch("Positive_Pz",  &m_positive_pz,  "Positive_Pz/D");
    m_refit1Tree->Branch("Positive_Phi", &m_positive_phi, "Positive_Phi/D");
    m_refit1Tree->Branch("Positive_eta", &m_positive_eta, "Positive_Eta/D");
    m_refit1Tree->Branch("Positive_z0",  &m_positive_z0,  "Positive_z0/D");
    m_refit1Tree->Branch("Positive_d0",  &m_positive_d0,  "Positive_d0/D");
    m_refit1Tree->Branch("Positive_z0_err",  &m_positive_z0_err,  "Positive_z0_err/D");
    m_refit1Tree->Branch("Positive_d0_err",  &m_positive_d0_err,  "Positive_d0_err/D");
    m_refit1Tree->Branch("Positive_sigma_pt",  &m_positive_sigma_pt,  "Positive_sigma_pt/D");
 
    if(m_doIP){
      m_refit1Tree->Branch("Negative_d0_PV",      &m_negative_d0_PV   ,  "Negative_d0_PV/D");
      m_refit1Tree->Branch("Negative_z0_PV",      &m_negative_z0_PV   ,  "Negative_z0_PV/D");
      m_refit1Tree->Branch("Positive_z0_PV",      &m_positive_z0_PV,     "Positive_z0_PV/D");
      m_refit1Tree->Branch("Positive_d0_PV",      &m_positive_d0_PV,     "Positive_d0_PV/D");
      
      m_refit1Tree->Branch("Negative_d0_PVerr",      &m_negative_d0_PVerr   ,  "Negative_d0_PVerr/D");
      m_refit1Tree->Branch("Negative_z0_PVerr",      &m_negative_z0_PVerr   ,  "Negative_z0_PVerr/D");
      m_refit1Tree->Branch("Positive_z0_PVerr",     &m_positive_z0_PVerr,     "Positive_z0_PVerr/D");
      m_refit1Tree->Branch("Positive_d0_PVerr",     &m_positive_d0_PVerr,     "Positive_d0_PVerr/D");
 
      m_refit1Tree->Branch("pv_x",      &m_pv_x   ,  "pv_x/D");
      m_refit1Tree->Branch("pv_y",      &m_pv_y   ,  "pv_y/D");
      m_refit1Tree->Branch("pv_z",      &m_pv_z   ,  "pv_z/D");
      m_refit1Tree->Branch("nTrkInVtx", &m_nTrkInVtx,  "nTrkInVtx/I");
    }    
  }
 
  // dealing with tree for Refit2 tracks
  isTreeNone = false;
  if (m_refit2TreeName.find("none") != std::string::npos) isTreeNone = true;
  if ( m_doRefit && !m_refit2Tree && !isTreeNone){
    ATH_MSG_INFO("initialize() ** defining m_refit2Tree with name: " << m_refit2TreeName.c_str());
    m_refit2Tree = new TTree((m_refit2TreeName).c_str(), m_ValidationTreeDescription.c_str());
 
    m_refit2Tree->Branch("runNumber"      ,  &m_runNumber,  "runNumber/I");
    m_refit2Tree->Branch("eventNumber"    ,  &m_evtNumber,  "eventNumber/I");
    m_refit2Tree->Branch("lumi_block"     ,  &m_lumi_block,  "lumi_block/I");
    m_refit2Tree->Branch("mu"             ,  &m_event_mu,  "mu/I");
    m_refit2Tree->Branch("preScale"       ,  &m_triggerPrescale, "preScale/I");
 
    m_refit2Tree->Branch("Negative_Px",  &m_negative_px,  "Negative_Px/D");
    m_refit2Tree->Branch("Negative_Py",  &m_negative_py,  "Negative_Py/D");
    m_refit2Tree->Branch("Negative_Pt",  &m_negative_pt,  "Negative_Pt/D");
    m_refit2Tree->Branch("Negative_Pz",  &m_negative_pz,  "Negative_Pz/D");
    m_refit2Tree->Branch("Negative_Phi", &m_negative_phi, "Negative_Phi/D");
    m_refit2Tree->Branch("Negative_Eta", &m_negative_eta, "Negative_Eta/D");
    m_refit2Tree->Branch("Negative_z0",  &m_negative_z0,  "Negative_z0/D");
    m_refit2Tree->Branch("Negative_d0",  &m_negative_d0,  "Negative_d0/D");
    m_refit2Tree->Branch("Negative_z0_err",  &m_negative_z0_err,  "Negative_z0_err/D");
    m_refit2Tree->Branch("Negative_d0_err",  &m_negative_d0_err,  "Negative_d0_err/D");
    m_refit2Tree->Branch("Negative_sigma_pt",  &m_negative_sigma_pt,  "Negative_sigma_pt/D");
 
    m_refit2Tree->Branch("Positive_Px",  &m_positive_px,  "Positive_Px/D");
    m_refit2Tree->Branch("Positive_Py",  &m_positive_py,  "Positive_Py/D");
    m_refit2Tree->Branch("Positive_Pt",  &m_positive_pt,  "Positive_Pt/D");
    m_refit2Tree->Branch("Positive_Pz",  &m_positive_pz,  "Positive_Pz/D");
    m_refit2Tree->Branch("Positive_Phi", &m_positive_phi, "Positive_Phi/D");
    m_refit2Tree->Branch("Positive_Eta", &m_positive_eta, "Positive_Eta/D");
    m_refit2Tree->Branch("Positive_z0",  &m_positive_z0,  "Positive_z0/D");
    m_refit2Tree->Branch("Positive_d0",  &m_positive_d0,  "Positive_d0/D");
    m_refit2Tree->Branch("Positive_z0_err",  &m_positive_z0_err,  "Positive_z0_err/D");
    m_refit2Tree->Branch("Positive_d0_err",  &m_positive_d0_err,  "Positive_d0_err/D");
    m_refit2Tree->Branch("Positive_sigma_pt",  &m_positive_sigma_pt,  "Positive_sigma_pt/D");
 
 
    if(m_doIP){
      m_refit2Tree->Branch("Negative_d0_PV",      &m_negative_d0_PV   ,  "Negative_d0_PV/D");
      m_refit2Tree->Branch("Negative_z0_PV",      &m_negative_z0_PV   ,  "Negative_z0_PV/D");
      m_refit2Tree->Branch("Positive_d0_PV",      &m_positive_d0_PV   ,  "Positive_d0_PV/D");
      m_refit2Tree->Branch("Positive_z0_PV",      &m_positive_z0_PV   ,  "Positive_z0_PV/D");
      
      
      m_refit2Tree->Branch("Negative_z0_PVerr",  &m_negative_z0_PVerr,  "Negative_z0_PVerr/D");
      m_refit2Tree->Branch("Negative_d0_PVerr",  &m_negative_d0_PVerr,  "Negative_d0_PVerr/D");
      m_refit2Tree->Branch("Positive_z0_PVerr",  &m_positive_z0_PVerr,  "Positive_z0_PVerr/D");
      m_refit2Tree->Branch("Positive_d0_PVerr",  &m_positive_d0_PVerr,  "Positive_d0_PVerr/D");
 
      m_refit2Tree->Branch("pv_x",      &m_pv_x   ,  "pv_x/D");
      m_refit2Tree->Branch("pv_y",      &m_pv_y   ,  "pv_y/D");
      m_refit2Tree->Branch("pv_z",      &m_pv_z   ,  "pv_z/D");
 
      m_refit2Tree->Branch("nTrkInVtx", &m_nTrkInVtx,  "nTrkInVtx/I");
    }    
  }
  
  isTreeNone = false;
  if (m_combTreeName.find("none") != std::string::npos) isTreeNone = true;
  if( !m_combTree && !m_skipMS && !isTreeNone){
    ATH_MSG_INFO("initialize() ** defining IDPerfMonZmumu m_combTree with name: " << m_combTreeName.c_str());
    m_combTree = new TTree((m_combTreeName).c_str(), m_ValidationTreeDescription.c_str());
 
    m_combTree->Branch("runNumber",   &m_runNumber,       "runNumber/I");
    m_combTree->Branch("eventNumber", &m_evtNumber,       "eventNumber/I");
    m_combTree->Branch("lumi_block",  &m_lumi_block,      "lumi_block/I");
    m_combTree->Branch("mu",          &m_event_mu,        "mu/I");
    m_combTree->Branch("preScale",    &m_triggerPrescale, "preScale/I");
 
    m_combTree->Branch("Negative_Px",       &m_negative_px,       "Negative_Px/D");
    m_combTree->Branch("Negative_Py",       &m_negative_py,       "Negative_Py/D");
    m_combTree->Branch("Negative_Pt",       &m_negative_pt,       "Negative_Pt/D");
    m_combTree->Branch("Negative_Pz",       &m_negative_pz,       "Negative_Pz/D");
    m_combTree->Branch("Negative_Phi",      &m_negative_phi,      "Negative_Phi/D");
    m_combTree->Branch("Negative_eta",      &m_negative_eta,      "Negative_eta/D");
    m_combTree->Branch("Negative_z0",       &m_negative_z0,       "Negative_z0/D");
    m_combTree->Branch("Negative_d0",       &m_negative_d0,       "Negative_d0/D");
    m_combTree->Branch("Negative_z0_err",   &m_negative_z0_err,   "Negative_z0_err/D");
    m_combTree->Branch("Negative_d0_err",   &m_negative_d0_err,   "Negative_d0_err/D");
    m_combTree->Branch("Negative_sigma_pt", &m_negative_sigma_pt, "Negative_sigma_pt/D");
 
    m_combTree->Branch("Positive_Px",  &m_positive_px,  "Positive_Px/D");
    m_combTree->Branch("Positive_Py",  &m_positive_py,  "Positive_Py/D");
    m_combTree->Branch("Positive_Pt",  &m_positive_pt,  "Positive_Pt/D");
    m_combTree->Branch("Positive_Pz",  &m_positive_pz,  "Positive_Pz/D");
    m_combTree->Branch("Positive_Phi", &m_positive_phi, "Positive_Phi/D");
    m_combTree->Branch("Positive_eta", &m_positive_eta, "Positive_Eta/D");
    m_combTree->Branch("Positive_z0",  &m_positive_z0,  "Positive_z0/D");
    m_combTree->Branch("Positive_d0",  &m_positive_d0,  "Positive_d0/D");
    m_combTree->Branch("Positive_z0_err",  &m_positive_z0_err,  "Positive_z0_err/D");
    m_combTree->Branch("Positive_d0_err",  &m_positive_d0_err,  "Positive_d0_err/D");
    m_combTree->Branch("Positive_sigma_pt",  &m_positive_sigma_pt,  "Positive_sigma_pt/D");
    
    if(m_doIP){
      m_combTree->Branch("Negative_z0_PV",  &m_negative_z0_PV,  "Negative_z0_PV/D");
      m_combTree->Branch("Negative_d0_PV",  &m_negative_d0_PV,  "Negative_d0_PV/D");
      m_combTree->Branch("Positive_z0_PV",  &m_positive_z0_PV,  "Positive_z0_PV/D");
      m_combTree->Branch("Positive_d0_PV",  &m_positive_d0_PV,  "Positive_d0_PV/D");
 
      m_combTree->Branch("Negative_z0_PVerr",  &m_negative_z0_PVerr,  "Negative_z0_PVerr/D");
      m_combTree->Branch("Negative_d0_PVerr",  &m_negative_d0_PVerr,  "Negative_d0_PVerr/D");
      m_combTree->Branch("Positive_z0_PVerr",  &m_positive_z0_PVerr,  "Positive_z0_PVerr/D");
      m_combTree->Branch("Positive_d0_PVerr",  &m_positive_d0_PVerr,  "Positive_d0_PVerr/D");
 
      m_combTree->Branch("pv_x",      &m_pv_x   ,  "pv_x/D");
      m_combTree->Branch("pv_y",      &m_pv_y   ,  "pv_y/D");
      m_combTree->Branch("pv_z",      &m_pv_z   ,  "pv_z/D");
 
      m_combTree->Branch("nTrkInVtx", &m_nTrkInVtx,  "nTrkInVtx/I");
    }
  }
 
  
  isTreeNone = false;
  if (m_MSTreeName.find("none") != std::string::npos) isTreeNone = true;
  if( !m_MSTree && !m_skipMS && !isTreeNone){
    ATH_MSG_INFO("initialize() ** defining IDPerfMon MSTree ");
 
    m_MSTree = new TTree(m_MSTreeName.c_str(), m_ValidationTreeDescription.c_str());
 
    m_MSTree->Branch("runNumber"      ,  &m_runNumber,  "runNumber/I");
    m_MSTree->Branch("eventNumber"      ,  &m_evtNumber,  "eventNumber/I");
    m_MSTree->Branch("lumi_block"      ,  &m_lumi_block,  "lumi_block/I");
    m_MSTree->Branch("preScale"       ,  &m_triggerPrescale, "preScale/I");
    
    m_MSTree->Branch("Negative_Px",  &m_negative_px,  "Negative_Px/D");
    m_MSTree->Branch("Negative_Py",  &m_negative_py,  "Negative_Py/D");
    m_MSTree->Branch("Negative_Pt",  &m_negative_pt,  "Negative_Pt/D");
    m_MSTree->Branch("Negative_Pz",  &m_negative_pz,  "Negative_Pz/D");
    m_MSTree->Branch("Negative_Phi", &m_negative_phi, "Negative_Phi/D");
    m_MSTree->Branch("Negative_Eta", &m_negative_eta, "Negative_Eta/D");
    m_MSTree->Branch("Negative_z0",  &m_negative_z0,  "Negative_z0/D");
    m_MSTree->Branch("Negative_d0",  &m_negative_d0,  "Negative_d0/D");
    m_MSTree->Branch("Negative_z0_err",  &m_negative_z0_err,  "Negative_z0_err/D");
    m_MSTree->Branch("Negative_d0_err",  &m_negative_d0_err,  "Negative_d0_err/D");
    m_MSTree->Branch("Negative_sigma_pt",  &m_negative_sigma_pt,  "Negative_sigma_pt/D");
 
    m_MSTree->Branch("Positive_Px",  &m_positive_px,  "Positive_Px/D");
    m_MSTree->Branch("Positive_Py",  &m_positive_py,  "Positive_Py/D");
    m_MSTree->Branch("Positive_Pt",  &m_positive_pt,  "Positive_Pt/D");
    m_MSTree->Branch("Positive_Pz",  &m_positive_pz,  "Positive_Pz/D");
    m_MSTree->Branch("Positive_Phi", &m_positive_phi, "Positive_Phi/D");
    m_MSTree->Branch("Positive_Eta", &m_positive_eta, "Positive_Eta/D");
    m_MSTree->Branch("Positive_z0",  &m_positive_z0,  "Positive_z0/D");
    m_MSTree->Branch("Positive_d0",  &m_positive_d0,  "Positive_d0/D");
    m_MSTree->Branch("Positive_z0_err",  &m_positive_z0_err,  "Positive_z0_err/D");
    m_MSTree->Branch("Positive_d0_err",  &m_positive_d0_err,  "Positive_d0_err/D");
    m_MSTree->Branch("Positive_sigma_pt",  &m_positive_sigma_pt,  "Positive_sigma_pt/D");
  }
 
  if( m_isMC && !m_truthTree ){
    m_truthTree = new TTree(m_truthTreeName.c_str(), m_ValidationTreeDescription.c_str());
 
    m_truthTree->Branch("runNumber"      ,  &m_runNumber,  "runNumber/I");
    m_truthTree->Branch("eventNumber"      ,  &m_evtNumber,  "eventNumber/I");
    m_truthTree->Branch("lumi_block"      ,  &m_lumi_block,  "lumi_block/I");
    m_truthTree->Branch("preScale"       ,  &m_triggerPrescale, "preScale/I");
 
    m_truthTree->Branch("Negative_Px",  &m_negative_px,  "Negative_Px/D");
    m_truthTree->Branch("Negative_Py",  &m_negative_py,  "Negative_Py/D");
    m_truthTree->Branch("Negative_Pt",  &m_negative_pt,  "Negative_Pt/D");
    m_truthTree->Branch("Negative_Pz",  &m_negative_pz,  "Negative_Pz/D");
    m_truthTree->Branch("Negative_Phi", &m_negative_phi, "Negative_Phi/D");
    m_truthTree->Branch("Negative_Eta", &m_negative_eta, "Negative_Eta/D");
    m_truthTree->Branch("Negative_z0",  &m_negative_z0,  "Negative_z0/D");
    m_truthTree->Branch("Negative_d0",  &m_negative_d0,  "Negative_d0/D");
    m_truthTree->Branch("Negative_z0_err",  &m_negative_z0_err,  "Negative_z0_err/D");
    m_truthTree->Branch("Negative_d0_err",  &m_negative_d0_err,  "Negative_d0_err/D");
    m_truthTree->Branch("Negative_parent",  &m_negative_parent,  "Negative_parent/I");
 
    m_truthTree->Branch("Positive_Px",  &m_positive_px,  "Positive_Px/D");
    m_truthTree->Branch("Positive_Py",  &m_positive_py,  "Positive_Py/D");
    m_truthTree->Branch("Positive_Pt",  &m_positive_pt,  "Positive_Pt/D");
    m_truthTree->Branch("Positive_Pz",  &m_positive_pz,  "Positive_Pz/D");
    m_truthTree->Branch("Positive_Phi", &m_positive_phi, "Positive_Phi/D");
    m_truthTree->Branch("Positive_Eta", &m_positive_eta, "Positive_Eta/D");
    m_truthTree->Branch("Positive_z0",  &m_positive_z0,  "Positive_z0/D");
    m_truthTree->Branch("Positive_d0",  &m_positive_d0,  "Positive_d0/D");
    m_truthTree->Branch("Positive_z0_err",  &m_positive_z0_err,  "Positive_z0_err/D");
    m_truthTree->Branch("Positive_d0_err",  &m_positive_d0_err,  "Positive_d0_err/D");
    m_truthTree->Branch("Positive_parent",  &m_positive_parent,  "Positive_parent/I");
  }
 
  if (m_doFourMuAnalysis && !m_FourMuTree) {
    ATH_MSG_INFO("initialize() ** defining m_FourMuTree ");
    m_FourMuTree = new TTree((m_FourMuTreeName).c_str(), "Four Muon monitoring");
 
    m_FourMuTree->Branch("runNumber"  ,  &m_runNumber,       "runNumber/I");
    m_FourMuTree->Branch("eventNumber",  &m_evtNumber,       "eventNumber/I");
    m_FourMuTree->Branch("lumi_block" ,  &m_lumi_block,      "lumi_block/I");
    m_FourMuTree->Branch("mu"         ,  &m_event_mu,        "mu/I");
    m_FourMuTree->Branch("preScale"   ,  &m_triggerPrescale, "preScale/I");
    m_FourMuTree->Branch("nVertex"    ,  &m_nVertex,         "nVertex/I");
 
    m_FourMuTree->Branch("Negative_1_Px",     &m_negative_px,   "Negative_1_Px/D");
    m_FourMuTree->Branch("Negative_1_Py",     &m_negative_py,   "Negative_1_Py/D");
    m_FourMuTree->Branch("Negative_1_Pz",     &m_negative_pz,   "Negative_1_Pz/D");
    m_FourMuTree->Branch("Negative_1_z0",     &m_negative_z0,   "Negative_1_z0/D");
    m_FourMuTree->Branch("Negative_1_d0",     &m_negative_d0,   "Negative_1_d0/D");
    m_FourMuTree->Branch("Negative_1_z0_err", &m_negative_z0_err,  "Negative_1_z0_err/D");
    m_FourMuTree->Branch("Negative_1_d0_err", &m_negative_d0_err,  "Negative_1_d0_err/D");
    m_FourMuTree->Branch("Negative_1_vtx",    &m_negative_1_vtx,  "Negative_1_vtx/I"); // vertex identifier
 
    m_FourMuTree->Branch("Negative_2_Px",     &m_negative_2_px,  "Negative_2_Px/D");
    m_FourMuTree->Branch("Negative_2_Py",     &m_negative_2_py,  "Negative_2_Py/D");
    m_FourMuTree->Branch("Negative_2_Pz",     &m_negative_2_pz,  "Negative_2_Pz/D");
    m_FourMuTree->Branch("Negative_2_z0",     &m_negative_2_z0,  "Negative_2_z0/D");
    m_FourMuTree->Branch("Negative_2_d0",     &m_negative_2_d0,  "Negative_2_d0/D");
    m_FourMuTree->Branch("Negative_2_z0_err", &m_negative_2_z0_err,  "Negative_2_z0_err/D");
    m_FourMuTree->Branch("Negative_2_d0_err", &m_negative_2_d0_err,  "Negative_2_d0_err/D");
    m_FourMuTree->Branch("Negative_2_vtx",    &m_negative_2_vtx,  "Negative_2_vtx/I");
 
    m_FourMuTree->Branch("Positive_1_Px",     &m_positive_px,     "Positive_1_Px/D");
    m_FourMuTree->Branch("Positive_1_Py",     &m_positive_py,     "Positive_1_Py/D");
    m_FourMuTree->Branch("Positive_1_Pz",     &m_positive_pz,     "Positive_1_Pz/D");
    m_FourMuTree->Branch("Positive_1_z0",     &m_positive_z0,     "Positive_1_z0/D");
    m_FourMuTree->Branch("Positive_1_d0",     &m_positive_d0,     "Positive_1_d0/D");
    m_FourMuTree->Branch("Positive_1_z0_err", &m_positive_z0_err, "Positive_1_z0_err/D");
    m_FourMuTree->Branch("Positive_1_d0_err", &m_positive_d0_err, "Positive_1_d0_err/D");
    m_FourMuTree->Branch("Positive_1_vtx",    &m_positive_1_vtx,  "Positive_1_vtx/I");
 
    m_FourMuTree->Branch("Positive_2_Px",     &m_positive_2_px,    "Positive_1_Px/D");
    m_FourMuTree->Branch("Positive_2_Py",     &m_positive_2_py,    "Positive_1_Py/D");
    m_FourMuTree->Branch("Positive_2_Pz",     &m_positive_2_pz,    "Positive_1_Pz/D");
    m_FourMuTree->Branch("Positive_2_z0",     &m_positive_2_z0,    "Positive_1_z0/D");
    m_FourMuTree->Branch("Positive_2_d0",     &m_positive_2_d0,    "Positive_1_d0/D");
    m_FourMuTree->Branch("Positive_2_z0_err", &m_positive_2_z0_err,"Positive_1_z0_err/D");
    m_FourMuTree->Branch("Positive_2_d0_err", &m_positive_2_d0_err,"Positive_1_d0_err/D");
    m_FourMuTree->Branch("Positive_2_vtx",    &m_positive_2_vtx,   "Positive_2_vtx/I");
 
    // electrons
    m_FourMuTree->Branch("ElNegative_1_Px",     &m_el_negative1_px,     "ElNegative_1_Px/D");
    m_FourMuTree->Branch("ElNegative_1_Py",     &m_el_negative1_py,     "ElNegative_1_Py/D");
    m_FourMuTree->Branch("ElNegative_1_Pz",     &m_el_negative1_pz,     "ElNegative_1_Pz/D");
    m_FourMuTree->Branch("ElNegative_1_z0",     &m_el_negative1_z0,     "ElNegative_1_z0/D");
    m_FourMuTree->Branch("ElNegative_1_d0",     &m_el_negative1_d0,     "ElNegative_1_d0/D");
    m_FourMuTree->Branch("ElNegative_1_z0_err", &m_el_negative1_z0_err, "ElNegative_1_z0_err/D");
    m_FourMuTree->Branch("ElNegative_1_d0_err", &m_el_negative1_d0_err, "ElNegative_1_d0_err/D");
    m_FourMuTree->Branch("ElNegative_1_vtx",    &m_el_negative1_vtx,    "ElNegative_1_vtx/I"); // vertex identifier
 
    m_FourMuTree->Branch("ElNegative_2_Px",     &m_el_negative2_px,     "ElNegative_2_Px/D");
    m_FourMuTree->Branch("ElNegative_2_Py",     &m_el_negative2_py,     "ElNegative_2_Py/D");
    m_FourMuTree->Branch("ElNegative_2_Pz",     &m_el_negative2_pz,     "ElNegative_2_Pz/D");
    m_FourMuTree->Branch("ElNegative_2_z0",     &m_el_negative2_z0,     "ElNegative_2_z0/D");
    m_FourMuTree->Branch("ElNegative_2_d0",     &m_el_negative2_d0,     "ElNegative_2_d0/D");
    m_FourMuTree->Branch("ElNegative_2_z0_err", &m_el_negative2_z0_err, "ElNegative_2_z0_err/D");
    m_FourMuTree->Branch("ElNegative_2_d0_err", &m_el_negative2_d0_err, "ElNegative_2_d0_err/D");
    m_FourMuTree->Branch("ElNegative_2_vtx",    &m_el_negative2_vtx,    "ElNegative_2_vtx/I"); // vertex identifier
 
    m_FourMuTree->Branch("ElPositive_1_Px",     &m_el_positive1_px,     "ElPositive_1_Px/D");
    m_FourMuTree->Branch("ElPositive_1_Py",     &m_el_positive1_py,     "ElPositive_1_Py/D");
    m_FourMuTree->Branch("ElPositive_1_Pz",     &m_el_positive1_pz,     "ElPositive_1_Pz/D");
    m_FourMuTree->Branch("ElPositive_1_z0",     &m_el_positive1_z0,     "ElPositive_1_z0/D");
    m_FourMuTree->Branch("ElPositive_1_d0",     &m_el_positive1_d0,     "ElPositive_1_d0/D");
    m_FourMuTree->Branch("ElPositive_1_z0_err", &m_el_positive1_z0_err, "ElPositive_1_z0_err/D");
    m_FourMuTree->Branch("ElPositive_1_d0_err", &m_el_positive1_d0_err, "ElPositive_1_d0_err/D");
    m_FourMuTree->Branch("ElPositive_1_vtx",    &m_el_positive1_vtx,    "ElPositive_1_vtx/I"); // vertex identifier
 
    m_FourMuTree->Branch("ElPositive_2_Px",     &m_el_positive2_px,     "ElPositive_2_Px/D");
    m_FourMuTree->Branch("ElPositive_2_Py",     &m_el_positive2_py,     "ElPositive_2_Py/D");
    m_FourMuTree->Branch("ElPositive_2_Pz",     &m_el_positive2_pz,     "ElPositive_2_Pz/D");
    m_FourMuTree->Branch("ElPositive_2_z0",     &m_el_positive2_z0,     "ElPositive_2_z0/D");
    m_FourMuTree->Branch("ElPositive_2_d0",     &m_el_positive2_d0,     "ElPositive_2_d0/D");
    m_FourMuTree->Branch("ElPositive_2_z0_err", &m_el_positive2_z0_err, "ElPositive_2_z0_err/D");
    m_FourMuTree->Branch("ElPositive_2_d0_err", &m_el_positive2_d0_err, "ElPositive_2_d0_err/D");
    m_FourMuTree->Branch("ElPositive_2_vtx",    &m_el_positive2_vtx,    "ElPositive_2_vtx/I"); // vertex identifier
 
    // other quantities
    m_FourMuTree->Branch("minv4mu",   &m_4mu_minv,  "minv4mu/D");
    m_FourMuTree->Branch("pv_x",      &m_pv_x   ,   "pv_x/D");
    m_FourMuTree->Branch("pv_y",      &m_pv_y   ,   "pv_y/D");
    m_FourMuTree->Branch("pv_z",      &m_pv_z   ,   "pv_z/D");
    m_FourMuTree->Branch("nTrkInVtx", &m_nTrkInVtx, "nTrkInVtx/I");
    
    m_FourMuTree->Branch("met",       &m_met,     "met/D");
    m_FourMuTree->Branch("metphi",    &m_metphi,  "metphi/D");
  }
  
  // now register the Trees
  ATH_MSG_INFO("initialize() Going to register the mu+mu- trees");
  ServiceHandle<ITHistSvc> tHistSvc("THistSvc", name());
  if (tHistSvc.retrieve().isFailure()){
        ATH_MSG_ERROR("initialize() Could not find Hist Service -> Switching ValidationMode Off !");
    m_validationMode = false;
  }
  
  //m_commonTree cannot be null here
  if ((tHistSvc->regTree(m_commonTreeFolder, m_commonTree)).isSuccess() ) {
    ATH_MSG_INFO("initialize() commonTree succesfully registered!");
  }
  else {
    ATH_MSG_ERROR("initialize() Could not register the validation commonTree -> Switching ValidationMode Off !");
    delete m_commonTree; m_commonTree = nullptr;
    m_validationMode = false;
  }
 
 
  if (m_IDTree) {
    if ((tHistSvc->regTree(m_IDTreeFolder, m_IDTree)).isSuccess() ) {
      ATH_MSG_INFO("initialize() IDTree succesfully registered!");
    }
    else {
      ATH_MSG_ERROR("initialize() Could not register the validation IDTree -> Switching ValidationMode Off !");
      delete m_IDTree; m_IDTree = nullptr;
      m_validationMode = false;
    }
  }
  
  if ( m_skipMS ) {
    ATH_MSG_INFO("CBtree has to be skipped: flag ON");
    delete m_combTree; m_combTree = nullptr;
  }
  else {
    if (m_combTree) {
      if ((tHistSvc->regTree(m_combTreeFolder, m_combTree)).isSuccess() ) {
    ATH_MSG_INFO("initialize() CBTree succesfully registered!");
      }
      else {
    ATH_MSG_ERROR("initialize() Could not register the validation CBTree -> Switching ValidationMode Off !");
    delete m_combTree; m_combTree = nullptr;
    m_validationMode = false;
      }
    }
  }
  
  if ( m_skipMS ) {
    ATH_MSG_INFO("MStree has to be skipped: flag ON");
    delete m_MSTree; m_MSTree = nullptr;
  }
  else {
    if (m_MSTree) {
      if ( (tHistSvc->regTree(m_MSTreeFolder, m_MSTree)).isSuccess()  ){
    ATH_MSG_INFO("initialize() MSTree succesfully registered!");
      }
      else {
    ATH_MSG_ERROR("initialize() Could not register the validation MSTree -> Switching ValidationMode Off !");
    delete m_MSTree; m_MSTree = nullptr;
    m_validationMode = false;
      }
    }
  }
 
  if( m_doRefit ){
    if (m_refit1Tree) {
      if ((tHistSvc->regTree(m_refit1TreeFolder, m_refit1Tree)).isSuccess() ) {
    ATH_MSG_INFO("initialize() Refit1 Tree succesfully registered!");
      }
      else{
    ATH_MSG_ERROR("initialize() Could not register the validation Tree -> Switching ValidationMode Off !");
    delete m_refit1Tree; m_refit1Tree = nullptr;
    m_validationMode = false;
      }
    }
 
    if (m_refit2Tree) {
      if ((tHistSvc->regTree(m_refit2TreeFolder, m_refit2Tree)).isSuccess() ) {
    ATH_MSG_INFO("initialize() Refit2 Tree succesfully registered!");
      }
      else {
    ATH_MSG_ERROR("initialize() Could not register the validation Tree -> Switching ValidationMode Off !");
    delete m_refit2Tree; m_refit2Tree = nullptr;
    m_validationMode = false;
      }
    }
  }
  
  if (m_isMC) {
    if ((tHistSvc->regTree(m_truthTreeFolder, m_truthTree)).isSuccess() ) {
      ATH_MSG_INFO("initialize() truthTree Tree succesfully registered!");
    }
    else {
      ATH_MSG_ERROR("initialize() Could not register the validation truth Tree -> Switching ValidationMode Off !");
      delete m_truthTree; m_truthTree = nullptr;
      m_validationMode = false;
    }
  }
  
  if (m_doFourMuAnalysis) {
    if ((tHistSvc->regTree(m_FourMuTreeFolder, m_FourMuTree)).isSuccess() ) {
      ATH_MSG_INFO("initialize() FourMu Tree succesfully registered!");
    }
    else {
      ATH_MSG_ERROR("initialize() Could not register the validation FourMu Tree -> Switching ValidationMode Off !");
      delete m_FourMuTree; m_FourMuTree = nullptr;
      m_validationMode = false;
    }
  }
 
  return StatusCode::SUCCESS;
}

CheckTriggerStatusAndPrescale()

StatusCode IDPerfMonZmumu::CheckTriggerStatusAndPrescale ( )

private

Definition at line 2304 of file IDPerfMonZmumu.cxx.

{
  m_triggerPrescale = 1; // default value
 
  // check trigger status
  if(m_triggerDecision.retrieve().isFailure()) {
    ATH_MSG_FATAL("CheckTriggerStatusAndPrescale -- Unable to retrieve " << m_triggerDecision << " turn it off");
    return StatusCode::FAILURE;
  }
  else {
    ATH_MSG_DEBUG("CheckTriggerStatusAndPrescale -- retrieved tool: " << m_triggerDecision );
  }
 
  float thisEventTriggerPrescale = 999999.9;
  float thisHLTTriggerPrescale = 999999.9;
  bool thisHLTIsPassed = false;
  std::string thisHLTTriggerName;
  std::string thisEventTriggerName;
  
  // HLT triggers
  const std::vector<std::string> myHLTtriglist = m_triggerDecision->getChainGroup("HLT_mu.*, HLT_2mu.*")->getListOfTriggers();
  for (int i=0; i < (int) myHLTtriglist.size(); i++) {
    if ( m_triggerDecision->isPassed(myHLTtriglist.at(i)) ) {
      ATH_MSG_DEBUG("HLT trigger = " << i << " out of " << myHLTtriglist.size() << " --> " << myHLTtriglist.at(i) <<  "  isPassed? " << m_triggerDecision->isPassed(myHLTtriglist.at(i)) << "  prescale: " << m_triggerDecision->getChainGroup(myHLTtriglist.at(i))->getPrescale() );
    }
    if ( m_triggerDecision->isPassed(myHLTtriglist.at(i)) ) { // trigger is passed
      thisHLTIsPassed = true;
      float thisprescale = m_triggerDecision->getChainGroup(myHLTtriglist.at(i))->getPrescale();
      if (thisprescale < thisHLTTriggerPrescale) {
    thisHLTTriggerPrescale = thisprescale;
    thisHLTTriggerName = myHLTtriglist.at(i);
      }
    }
  }
  if (thisHLTTriggerPrescale < 1) thisHLTTriggerPrescale = 1;
  ATH_MSG_DEBUG("Event HLT trigger prescale = " << thisHLTTriggerPrescale);
 
  //if (thisL1IsPassed && thisHLTIsPassed) {
  if (thisHLTIsPassed) {
    ATH_MSG_DEBUG("Trigger passed -> accept event");
  }
  else{
    ATH_MSG_DEBUG("Trigger Failed -> reject event");
    return StatusCode::FAILURE;
  }
 
  // event prescale
  thisEventTriggerPrescale = std::move(thisHLTTriggerPrescale);
  thisEventTriggerName = thisHLTTriggerName;
  ATH_MSG_DEBUG("CheckTriggerStatusAndPrescale -- Event trigger prescale = " << thisEventTriggerPrescale);
  ATH_MSG_DEBUG("CheckTriggerStatusAndPrescale -- Event trigger name = " << thisEventTriggerName);
  
  m_triggerPrescale = thisEventTriggerPrescale;
  m_triggerName = std::move(thisEventTriggerName);
 
  
  return StatusCode::SUCCESS;
}

Clear4MuNtupleVariables()

void IDPerfMonZmumu::Clear4MuNtupleVariables ( )

private

Definition at line 2454 of file IDPerfMonZmumu.cxx.

{
  m_positive_px = 0.;
  m_positive_py = 0.;
  m_positive_pz = 0.;
  m_positive_d0 = 0.;
  m_positive_z0 = 0.;
  m_positive_d0_err = 0.;
  m_positive_z0_err = 0.;
  m_positive_sigma_pt = 0.;
  m_positive_1_vtx = 0;
  m_positive_parent = 0;
 
  m_negative_px = 0.;
  m_negative_py = 0.;
  m_negative_pz = 0.;
  m_negative_d0 = 0.;
  m_negative_z0 = 0.;
  m_negative_d0_err = 0.;
  m_negative_z0_err = 0.;
  m_negative_sigma_pt = 0.;
  m_negative_1_vtx = 0;
  m_negative_parent = 0;
    
  m_positive_2_px = 0.;
  m_positive_2_py = 0.;
  m_positive_2_pz = 0.;
  m_positive_2_d0 = 0.;
  m_positive_2_z0 = 0.;
  m_positive_2_d0_err = 0.;
  m_positive_2_z0_err = 0.;
  m_positive_2_vtx = 0;
    
  m_negative_2_px = 0.;
  m_negative_2_py = 0.;
  m_negative_2_pz = 0.;
  m_negative_2_d0 = 0.;
  m_negative_2_z0 = 0.;
  m_negative_2_d0_err = 0.;
  m_negative_2_z0_err = 0.;
  m_negative_2_vtx = 0;
 
  // electrons
  m_el_negative1_px = 0.;
  m_el_negative1_py = 0.;
  m_el_negative1_pz = 0.;
  m_el_negative1_d0 = 0.;
  m_el_negative1_z0 = 0.;
  m_el_negative1_d0_err = 0.;
  m_el_negative1_z0_err = 0.;
  m_el_negative1_vtx = 0;
 
  m_el_negative2_px = 0.;
  m_el_negative2_py = 0.;
  m_el_negative2_pz = 0.;
  m_el_negative2_d0 = 0.;
  m_el_negative2_z0 = 0.;
  m_el_negative2_d0_err = 0.;
  m_el_negative2_z0_err = 0.;
  m_el_negative2_vtx = 0;
 
  m_el_positive1_px = 0.;
  m_el_positive1_py = 0.;
  m_el_positive1_pz = 0.;
  m_el_positive1_d0 = 0.;
  m_el_positive1_z0 = 0.;
  m_el_positive1_d0_err = 0.;
  m_el_positive1_z0_err = 0.;
  m_el_positive1_vtx = 0;
  
  m_el_positive2_px = 0.;
  m_el_positive2_py = 0.;
  m_el_positive2_pz = 0.;
  m_el_positive2_d0 = 0.;
  m_el_positive2_z0 = 0.;
  m_el_positive2_d0_err = 0.;
  m_el_positive2_z0_err = 0.;
  m_el_positive2_vtx = 0;
  
  // other quantities
  m_nVertex =        0;
  m_pv_x = 0; m_pv_y = 0; m_pv_z = 0;
  m_4mu_minv = 0.;
 
  return;
}

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl, const SG::VarHandleKeyType &  )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode IDPerfMonZmumu::execute ( const EventContext & ctx )

virtual

Execute method.

– START 4 lepton analysis

– END 4 lepton analysis

Implements AthAlgorithm.

Definition at line 934 of file IDPerfMonZmumu.cxx.

{
  ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** START **");
  
  SG::ReadHandle<xAOD::EventInfo> eventInfo (m_EventInfoKey, ctx);
  SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandleRec { m_beamSpotKey, ctx };
 
  if(eventInfo.isValid()) {
    m_runNumber  = eventInfo->runNumber();
    m_evtNumber  = eventInfo->eventNumber();
    m_lumi_block = eventInfo->lumiBlock();
    m_event_mu   = eventInfo->actualInteractionsPerCrossing();
    m_beamposX   = beamSpotHandleRec->beamPos().x();
    m_beamposY   = beamSpotHandleRec->beamPos().y();
 
    if (eventInfo->mcEventWeights().size()>0) { 
      m_event_weight = eventInfo->mcEventWeights()[0];
    }
    else {
      m_event_weight = 1.; // default
    }      
    ATH_MSG_DEBUG(" Execute() starting on --> Run: " << m_runNumber << "  event: " << m_evtNumber << "   Lumiblock: " << m_lumi_block << "    weight:" << m_event_weight);
  }
  else {
    ATH_MSG_DEBUG(" IDPerfMonZmumu::execute evtStore->retrieve (eventInfo) failed ..  trying another strategy.. ");
    SG::ReadHandle<xAOD::EventInfo> eventInfo2(m_EventInfoKey, ctx);
    if (eventInfo2.isValid()) {
      m_runNumber = eventInfo2->runNumber();
      m_evtNumber = eventInfo2->eventNumber();
      m_lumi_block = eventInfo2->lumiBlock();
      m_event_mu = eventInfo2->actualInteractionsPerCrossing();
      if (eventInfo->mcEventWeights().size()>0) { 
    m_event_weight = eventInfo->mcEventWeights()[0];
      }
      else {
    m_event_weight = 1.; // default
      }      
      ATH_MSG_DEBUG(" Execute() starting on --> Run: " << m_runNumber << "  event: " << m_evtNumber << "   Lumiblock: " << m_lumi_block << "    weight:" << m_event_weight);
    }      
    else{
      ATH_MSG_ERROR("** IDPerfMonZmumu::execute ** Could not retrieve event info.");
    }
  }

  if (m_doFourMuAnalysis) {
    ATH_MSG_DEBUG(" ** IDPerfMonZmumu::execute ** calling FourLeptonAnalysis()...");
    StatusCode fourLeptAnaStatus = this->RunFourLeptonAnalysis();
    if (fourLeptAnaStatus.isSuccess()) ATH_MSG_INFO (" ** IDPerfMonZmumu::execute ** RunFourLeptonAnalysis() SUCCESS -> found a new event");
  }
 
 
  // check if the muon-pair passed the resonance selection cuts:
  ATH_MSG_DEBUG(" ** IDPerfMonZmumu::execute ** calling dimuon analysis m_xZmm.Reco()...");
  if( m_xZmm.Reco( m_lumi_block ) ){
    ATH_MSG_INFO(   "  Run: " << m_runNumber 
            << "  event: " << m_evtNumber 
            << "  Lumiblock: " << m_lumi_block 
            << "  Invariant mass = " << m_xZmm.GetInvMass() << " GeV "
            << "  weight: " << m_event_weight 
            << "  ** SUCCESS **");
  }
  else {
    ATH_MSG_DEBUG(   "  Run: " << m_runNumber
          << "  event: " << m_evtNumber 
          << "  Lumiblock: " << m_lumi_block
          << "  Failed dimuon reconstruction. m_xZmm.Reco() returned FALSE ");
    return StatusCode::SUCCESS;
  }
 
  //
  // If this point is reached -> there is a good mu+mu- pair that stisfies all selection cuts
  //
 
  if (m_UseTrigger) {
    StatusCode isTriggerPassed = CheckTriggerStatusAndPrescale ();
    if (isTriggerPassed == StatusCode::SUCCESS) {
      ATH_MSG_DEBUG("Trigger passed -> accept event");
    }
    else{
      ATH_MSG_DEBUG("Trigger Failed -> reject event --> leave event");
      return StatusCode::SUCCESS;
    }
  }
 
  // std::cout << " ** IDPerfMonZmumu ** extracting muon_pos and muon_neg... " << std::endl;
  const xAOD::Muon* muon_pos = m_xZmm.getCombMuon(m_xZmm.getPosMuon(ZmumuEvent::CB));
  const xAOD::Muon* muon_neg = m_xZmm.getCombMuon(m_xZmm.getNegMuon(ZmumuEvent::CB));
 
  const xAOD::TrackParticle* ppos_comb = nullptr; const xAOD::Vertex* ppos_comb_v = nullptr;
  const xAOD::TrackParticle* pneg_comb = nullptr; const xAOD::Vertex* pneg_comb_v = nullptr;
    
  //To protect against failures of the estimation
  StatusCode success_pos = StatusCode::FAILURE;
  StatusCode success_neg = StatusCode::FAILURE;
 
  if (muon_pos && muon_neg) { // if both combined muons exist and were sucessfully retrieved        
    ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** combined muons exist ** retrieving their m_trackparticleName: " << m_trackParticleName.c_str());
    using enum xAOD::Muon::TrackParticleType;
    if (m_trackParticleName.find("InnerDetectorTrackParticles") != std::string::npos) {
      ATH_MSG_INFO("** IDPerfMonZmumu::execute ** Retrieving InnerDetectorTrackParticles of the accepted muons");
      ppos_comb = muon_pos->trackParticle(InnerDetectorTrackParticle);
      pneg_comb = muon_neg->trackParticle(InnerDetectorTrackParticle);
      
      if (!ppos_comb || !pneg_comb) {
    ATH_MSG_WARNING("** IDPerfMonZmumu::execute ** InnerDetectorTrackParticles are requested but they are not present. Exiting event.");
    return StatusCode::SUCCESS;
      }
    } // InnerDetectorTrackParticles
    if (m_trackParticleName.find("CombinedTrackParticle") != std::string::npos) {
      // 
      ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** Retrieving CombinedTrackParticles of the accepted muons");
      ppos_comb = muon_pos->trackParticle(CombinedTrackParticle);
      pneg_comb = muon_neg->trackParticle(CombinedTrackParticle);
      
      if (!ppos_comb || !pneg_comb){
    ATH_MSG_WARNING( "** IDPerfMonZmumu::execute ** CombinedTrackParticles are requested but they are not present. Exiting event.");
    return StatusCode::SUCCESS;
      }      
    } // combined tracks
    
 
    // double check
    if (ppos_comb && pneg_comb) {
      ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** successfull retrieval of muons as " << m_trackParticleName << ". ppos_comb & pneg_comb both exist");
    }
    else {
      ATH_MSG_WARNING("** IDPerfMonZmumu::execute ** problems retrieving muons as " << m_trackParticleName <<". ppos_comb or pneg_comb are not available");
    }
 
 
    // vertex
    SG::ReadHandle<xAOD::VertexContainer> vertices { m_vertexKey, ctx };
    for (const auto V : *vertices) {
      if (V->vertexType() == xAOD::VxType::VertexType::PriVtx) {
    // primaryVertex =V;
    ppos_comb_v = V;
    pneg_comb_v = V;
    break;
      }
    } // vertex
    
    // before continuing check both particles have vertex
    if (!ppos_comb_v || !pneg_comb_v){
      ATH_MSG_WARNING( "Some or all of the requested particles have no vertex. Exiting event");
      return StatusCode::SUCCESS;
    }      
  } // if (muon_pos && muon_neg) 
  else { // this is the else of --> if (muon_pos && muon_neg)
    ATH_MSG_WARNING("** IDPerfMonZmumu::execute ** Could not find CombinedMuon pos/neg in event"
            << "  Run: " << m_runNumber 
            << "  event: " << m_evtNumber 
            << "  Lumiblock: " << m_lumi_block );
    return StatusCode::SUCCESS;
  }
  
  // 
  // If this point is reached -> there is a good mu+mu- pair and the muons have been associated to a vertex
  // So, ntuple can be filled
  //
  if ( m_xZmm.AcceptEvent() ) { 
    ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** Going to fill ntuples for Run: " << m_runNumber 
          << "  event: " << m_evtNumber 
          << "  Lumiblock: " << m_lumi_block 
          << "  Invariant mass = " << m_xZmm.GetInvMass() << " GeV  ** some preliminaries are due though...");
 
    if (m_useTrackSelectionTool) {
      // The track selection is applied to the InnerDetectorTrackParticle of the muons
      if ( !m_selTool->accept(muon_pos->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle)) ) {
    ATH_MSG_DEBUG("Exiting because the ID segment of muon_pos do not pass the TrackSelection");
    return StatusCode::SUCCESS;
      } 
      if ( !m_selTool->accept(muon_neg->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle)) ) {
    ATH_MSG_DEBUG("Exiting because the ID segment of muon_neg do not pass the TrackSelection");
    return StatusCode::SUCCESS;
      }
      ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** Good news ** pos and neg muon->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle) pass the trackSelection :)");
    } 
    else {
      ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** Track selection is not in use ");
    }
    
    ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** >> before fill rec with default: " << m_trackParticleName << " tracks << ");
    
    if (m_storeZmumuNtuple) {
      // reset vectors
      this->ResetCommonNtupleVectors();
      
      // Fill Inner Detector Tree 
      ATH_MSG_DEBUG("-- >> going to fill ID muons params << --");
      success_pos = FillRecParametersTP (muon_pos->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle), 
                     muon_pos->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle), 
                     muon_pos->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle)->charge(),
                     ppos_comb_v, ctx);
      success_neg = FillRecParametersTP (muon_neg->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle), 
                     muon_neg->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle), 
                     muon_neg->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle)->charge(),
                     pneg_comb_v, ctx);
      
      if (success_pos && success_neg) {
    // before filling the ntuple, extract the hits information. This is filled once and it is used for all track collections
    this->ExtractIDHitsInformation(muon_pos, muon_neg);
    
    // Fill ntuple
    if (m_IDTree != nullptr) {
      ATH_MSG_DEBUG("-- Filling m_IDTree ntuple " << m_IDTree->GetName() << " entry " << m_IDTree->GetEntries()
            << " for run: " << m_runNumber 
            << "  event: " << m_evtNumber 
            << "  Lumiblock: " << m_lumi_block 
            << "  Invariant mass = " << m_xZmm.GetInvMass() << " GeV ");
    }
    ATH_MSG_DEBUG("mu+  --> pxyz " << m_positive_px 
              << ", " << m_positive_py
              << ", " << m_positive_pz
              << "  pt: " << muon_pos->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle)->pt()
              << "  d0: " << m_positive_d0 
              << "  z0: " << m_positive_z0
              << "  d0unc: " << m_positive_d0_err
              << "  sigma_pt: " << m_positive_sigma_pt);
    ATH_MSG_DEBUG("mu-  --> pxyz: " << m_negative_px
              << ", " << m_negative_py
              << ", " << m_negative_pz
              << "  pt: " << muon_neg->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle)->pt()
              << "  d0: " << m_negative_d0
              << "  z0: " << m_negative_z0
              << "  d0unc: " << m_negative_d0_err
              << "  sigma_pt: " << m_negative_sigma_pt);
    // ntuple variables have been filled in FillRecParametersTP
    if (m_IDTree) m_IDTree->Fill();
 
    // Fill vectors for common ntuple
    m_IDTrack_pt.push_back(m_positive_pt);
    m_IDTrack_pt.push_back(m_negative_pt);
    m_IDTrack_eta.push_back(m_positive_eta);
    m_IDTrack_eta.push_back(m_negative_eta);
    m_IDTrack_phi.push_back(m_positive_phi);
    m_IDTrack_phi.push_back(m_negative_phi);
    m_IDTrack_d0.push_back(m_positive_d0);
    m_IDTrack_d0.push_back(m_negative_d0);
    m_IDTrack_z0.push_back(m_positive_z0);
    m_IDTrack_z0.push_back(m_negative_z0);
    m_IDTrack_qoverp.push_back(m_positive_qoverp);
    m_IDTrack_qoverp.push_back(m_negative_qoverp);
    
    m_IDTrack_sigma_pt.push_back(m_positive_sigma_pt);
    m_IDTrack_sigma_pt.push_back(m_negative_sigma_pt);
    m_IDTrack_sigma_d0.push_back(m_positive_d0_err);
    m_IDTrack_sigma_d0.push_back(m_negative_d0_err);
    m_IDTrack_sigma_z0.push_back(m_positive_z0_err);
    m_IDTrack_sigma_z0.push_back(m_negative_z0_err);
    m_IDTrack_sigma_qoverp.push_back(m_positive_sigma_qoverp);
    m_IDTrack_sigma_qoverp.push_back(m_negative_sigma_qoverp);
      }
      else {
    ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** not filling IDTracks in combTree due to problems with muon_xxx->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle)");
      }
      // End of fill ID Tree 
      
      //
      // combined muons ntuple
      if ((!m_skipMS && m_combTree) || m_commonTree ) { // if skipMS -> no combined muons
    success_pos = FillRecParametersTP(muon_pos->trackParticle(xAOD::Muon::TrackParticleType::CombinedTrackParticle), 
                      muon_pos->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle), 
                      ppos_comb->charge(), 
                      ppos_comb_v, ctx);
    ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** success_pos TP pt: " << m_positive_pt << " GeV");
    success_neg = FillRecParametersTP(muon_neg->trackParticle(xAOD::Muon::TrackParticleType::CombinedTrackParticle), 
                      muon_neg->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle), 
                      pneg_comb->charge(), 
                      pneg_comb_v, ctx);
    ATH_MSG_DEBUG("** IDPerfMonZmumu::excute ** success_neg TP pt: " << m_negative_pt << " GeV");
 
    if (success_pos && success_neg) {
      if (m_combTree != nullptr) {
        ATH_MSG_DEBUG("-- Filling m_combTree ntuple " << m_combTree->GetName() << " entry " << m_combTree->GetEntries() 
              << "  run: " << m_runNumber 
              << "  event: " << m_evtNumber 
              << "  Lumiblock: " << m_lumi_block 
              << "  Invariant mass = " << m_xZmm.GetInvMass() << " GeV ");
      }
      ATH_MSG_DEBUG("mu+ --> pxyz " << m_positive_px 
            << ", " << m_positive_py
            << ", " << m_positive_pz
            << "  pt: " << m_positive_pt
            << "  d0: " << m_positive_d0 
            << "  z0: " << m_positive_z0
            << "  d0unc: " << m_positive_d0_err);
      ATH_MSG_DEBUG("mu- --> pxyz: " << m_negative_px
            << ", " << m_negative_py
            << ", " << m_negative_pz 
            << "  pt: " << m_negative_pt
            << "  d0: " << m_negative_d0
            << "  z0: " << m_negative_z0
            << "  d0unc: " << m_negative_d0_err);
      // ntuple variables have been filled in FillRecParameters
      if (m_combTree) m_combTree->Fill();
    }
    // Fill vectors for common ntuple
    m_CBTrack_pt.push_back(m_positive_pt);
    m_CBTrack_pt.push_back(m_negative_pt);
    m_CBTrack_eta.push_back(m_positive_eta);
    m_CBTrack_eta.push_back(m_negative_eta);
    m_CBTrack_phi.push_back(m_positive_phi);
    m_CBTrack_phi.push_back(m_negative_phi);
    m_CBTrack_d0.push_back(m_positive_d0);
    m_CBTrack_d0.push_back(m_negative_d0);
    m_CBTrack_z0.push_back(m_positive_z0);
    m_CBTrack_z0.push_back(m_negative_z0);
    m_CBTrack_qoverp.push_back(m_positive_qoverp);
    m_CBTrack_qoverp.push_back(m_negative_qoverp);
    
    m_CBTrack_sigma_pt.push_back(m_positive_sigma_pt);
    m_CBTrack_sigma_pt.push_back(m_negative_sigma_pt);
    m_CBTrack_sigma_d0.push_back(m_positive_d0_err);
    m_CBTrack_sigma_d0.push_back(m_negative_d0_err);
    m_CBTrack_sigma_z0.push_back(m_positive_z0_err);
    m_CBTrack_sigma_z0.push_back(m_negative_z0_err);
    m_CBTrack_sigma_qoverp.push_back(m_positive_sigma_qoverp);
    m_CBTrack_sigma_qoverp.push_back(m_negative_sigma_qoverp);
      } // if skipMS -> no combined muons
 
      
      
      // MS ntuple
      if (!m_skipMS && m_MSTree) { // if skipMS -> no MS tracks
    ATH_MSG_DEBUG("-- >> going to fill MS muons params << --");
    success_pos = FillRecParametersTP(m_xZmm.getMSTrack(m_xZmm.getPosMuon(ZmumuEvent::CB)),
                      muon_pos->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle),
                      ppos_comb->charge(), 
                      ppos_comb_v, ctx);
 
    success_neg = FillRecParametersTP(m_xZmm.getMSTrack(m_xZmm.getNegMuon(ZmumuEvent::CB)),
                      muon_neg->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle),
                      pneg_comb->charge(),
                      pneg_comb_v, ctx);
    if (success_pos && success_neg) { 
      if (m_MSTree != nullptr) ATH_MSG_DEBUG("-- Filling m_MSTree " << m_MSTree->GetName() << " entry " << m_MSTree->GetEntries() 
                         << "  run: " << m_runNumber 
                         << "  event: " << m_evtNumber 
                         << "  Lumiblock: " << m_lumi_block 
                         << "  Invariant mass = " << m_xZmm.GetInvMass() << " GeV ");
      ATH_MSG_DEBUG("mu+ --> pxyz " << m_positive_px 
            << ", " << m_positive_py
            << "  pt: " << m_xZmm.getMSTrack(m_xZmm.getPosMuon(ZmumuEvent::CB))->pt()
            << "  d0: " << m_positive_d0 
            << "  z0: " << m_positive_z0
            << "  d0unc: " << m_positive_d0_err);
      ATH_MSG_DEBUG("mu- --> pxyz: " << m_negative_px
            << ", " << m_negative_py
            << "  pt: " << m_xZmm.getMSTrack(m_xZmm.getNegMuon(ZmumuEvent::CB))->pt()
            << "  d0: " << m_negative_d0
            << "  z0: " << m_negative_z0
            << "  d0unc: " << m_negative_d0_err);
      // ntuple variables have been filled in FillRecParametersTP
      if (m_MSTree) m_MSTree->Fill();
    }
    else {
      ATH_MSG_INFO("FAILED filling m_MSTree " 
               << "  for run: " << m_runNumber 
               << "  event: " << m_evtNumber 
               << "  Lumiblock: " << m_lumi_block);
    }
    // MS ntuple
      }
    } // if store Zmumu ntuple
    
    // changed refitting to combined particles since run II DESDM_ZMUMU did not store InDetTrackParticles
    if( m_doRefit ) {
      ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** Going to build TrackCollections: muonTrks, muonTrksRefit1 and muonTrksRefit2");
      auto muonTrks        = std::make_unique<TrackCollection>(SG::OWN_ELEMENTS);
      auto muonTrksRefit1  = std::make_unique<TrackCollection>(SG::OWN_ELEMENTS);
      auto muonTrksRefit2  = std::make_unique<TrackCollection>(SG::OWN_ELEMENTS);
          
      ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** Going to refit and store the track parameters "); 
      std::unique_ptr<Trk::Track> defaultMuonTrk1;
      std::unique_ptr<Trk::Track> defaultMuonTrk2;
      Trk::Track* refit1MuonTrk1{};
      Trk::Track* refit2MuonTrk1{};
      Trk::Track* refit1MuonTrk2{};
      Trk::Track* refit2MuonTrk2{};
      
      StatusCode fitStatus;
      //save default and refit track parameters
      if( ppos_comb->track() ) {
        defaultMuonTrk1 = std::make_unique<Trk::Track>(*ppos_comb->track());
        IegammaTrkRefitterTool::Cache fitResult{};
        fitStatus = m_TrackRefitter1->refitTrack( ctx, ppos_comb->track(), fitResult );  
        if (fitStatus.isFailure()) {
          ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** Track Refit1 Failed for ppos_comb->track(). Skipping Event");
          return StatusCode::SUCCESS;
        } else {
          refit1MuonTrk1 = fitResult.refittedTrack.release();
          muonTrksRefit1->push_back(refit1MuonTrk1);
          ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** Track Refit1 Success of ppos_comb->track()." 
            << " Original pt: " << ppos_comb->track()->perigeeParameters()->pT() 
            << "  track refit pt: " << refit1MuonTrk1->perigeeParameters()->pT() );
        }
        fitStatus = m_TrackRefitter2->refitTrack( ctx, ppos_comb->track(), fitResult );
        if (fitStatus.isFailure()) {
          ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** Track Refit2 Failed for ppos_comb->track(). Skipping Event");
          return StatusCode::SUCCESS;
        } else {
          refit2MuonTrk1 = fitResult.refittedTrack.release();
          muonTrksRefit2->push_back(refit2MuonTrk1);
          ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** Track Refit2 Success of ppos_comb->track()." 
            << " Original pt: " << ppos_comb->track()->perigeeParameters()->pT() 
            << "  track refit pt: " << refit1MuonTrk1->perigeeParameters()->pT() );
        }
      }
      
      if( pneg_comb->track() ) {
        defaultMuonTrk2 = std::make_unique<Trk::Track>(*pneg_comb->track());
        
        IegammaTrkRefitterTool::Cache fitResult {};
        fitStatus = m_TrackRefitter1->refitTrack( ctx, pneg_comb->track(), fitResult);
        if (fitStatus.isFailure()) {
          ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** Track Refit1 Failed. Skipping Event");
          return StatusCode::SUCCESS;
        } else {
          refit1MuonTrk2 = fitResult.refittedTrack.release();
          muonTrksRefit1->push_back(refit1MuonTrk2);
          ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** Track Refit1 Success of pneg_comb->track()." 
            << " Original pt: " << pneg_comb->track()->perigeeParameters()->pT() 
            << "  track refit pt: " << refit1MuonTrk2->perigeeParameters()->pT() );
        }
    
    
        fitStatus = m_TrackRefitter2->refitTrack( ctx, pneg_comb->track(), fitResult );
        if (fitStatus.isFailure()) {
          ATH_MSG_DEBUG("Track Refit2 Failed. Skipping Event");
          return StatusCode::SUCCESS;
        } else {
          refit2MuonTrk2 = fitResult.refittedTrack.release();
          muonTrksRefit2->push_back(refit2MuonTrk2);
          ATH_MSG_DEBUG("** IDPerfMonZmumu::execute ** Track Refit2 Success of pneg_comb->track()." 
            << " Original pt: " << pneg_comb->track()->perigeeParameters()->pT() 
            << "  track refit pt: " << refit2MuonTrk2->perigeeParameters()->pT() );
        }
        if (muonTrksRefit1->size() != 2) ATH_MSG_WARNING("** IDPerfMonZmumu::execute ** WARNING ** size of muonTrksRefit1: " << muonTrksRefit1->size());
        if (muonTrksRefit2->size() != 2) ATH_MSG_WARNING("** IDPerfMonZmumu::execute ** WARNING ** size of muonTrksRefit2: " << muonTrksRefit2->size());
      }
      
      //save tracks to storegrate, they can be used by InDetAlignmentMonitoring
      muonTrks->push_back(defaultMuonTrk1.release());
      muonTrks->push_back(defaultMuonTrk2.release());
      auto nMuonTracks = muonTrks->size();
      StatusCode sc = evtStore()->record(std::move(muonTrks), m_outputTracksName);
      if (sc.isSuccess()) {
        ATH_MSG_DEBUG ("Stored "<< nMuonTracks << " " << m_outputTracksName <<" into StoreGate");
      } else{
       ATH_MSG_WARNING ("Failed storing " << m_outputTracksName);
      }
      auto nMuonTracks1 = muonTrksRefit1->size();
      if (nMuonTracks1 > 1) { // keep track collection if at least 2 muons were found
        sc = evtStore()->record(std::move(muonTrksRefit1), m_outputTracksName + "Refit1");
        if (sc.isSuccess()) {
          ATH_MSG_DEBUG ("Stored "<< nMuonTracks1 << " " << m_outputTracksName + "Refit1" << " into StoreGate");
        } else {
          ATH_MSG_WARNING ("Failed storing " << m_outputTracksName + "Refit1");
        }
      }
      auto nMuonTracks2 = muonTrksRefit2->size();
      if (nMuonTracks2 > 1) { // keep track collection if at least 2 muons were found
        sc = evtStore()->record(std::move(muonTrksRefit2), m_outputTracksName + "Refit2");
        if (sc.isSuccess()) {
          ATH_MSG_DEBUG ("Stored "<< nMuonTracks2 << " " << m_outputTracksName + "Refit2" << " into StoreGate");
        } else {
          ATH_MSG_WARNING ("Failed storing " << m_outputTracksName + "Refit2");
        }
      }
      
      //fill refit1 ID parameters
      if (nMuonTracks1<2) {
          ATH_MSG_WARNING("** IDPerfMonZmumu::execute ** Refit1 muon tracks are missing!");
      } else {
        ATH_MSG_DEBUG("** IDPerfMonZmumu::execute **  going to fill refit1tree ");
        success_pos = FillRecParametersSimple (refit1MuonTrk1, ppos_comb->charge(), ppos_comb_v, ctx);
        success_neg = FillRecParametersSimple (refit1MuonTrk2, pneg_comb->charge(), pneg_comb_v, ctx);
        
        if (m_storeZmumuNtuple) {
          if (success_pos && success_neg) {
            if (m_refit1Tree) {
              ATH_MSG_DEBUG("-- Filling ntuple " << m_refit1Tree->GetName() << " entry " << m_refit1Tree->GetEntries() 
                << " for run: " << m_runNumber 
                << "  event: " << m_evtNumber 
                << "  Lumiblock: " << m_lumi_block 
                << "  Invariant mass = " << m_xZmm.GetInvMass() << " GeV ");
            }
            ATH_MSG_DEBUG("Negative p: (" << m_negative_px << ", "
              << m_negative_py << ", "
              << m_negative_pz << ") "
              << "  d0: " << m_negative_d0
              << "  z0: " << m_negative_z0
              << "  sigma_pt: " << m_positive_sigma_pt);
            ATH_MSG_DEBUG("Positive p: (" << m_positive_px << ","
              << m_positive_py << ","
              << m_positive_pz << ") "
              << "  d0: " << m_positive_d0
              << "  z0: " << m_positive_z0    
              << "  sigma_pt: " << m_negative_sigma_pt);
            if (m_refit1Tree) {
              m_refit1Tree->Fill();
            }
            // Fill vectors for common ntuple
            m_Refit1_pt.push_back(m_positive_pt);
            m_Refit1_pt.push_back(m_negative_pt);       
            m_Refit1_eta.push_back(m_positive_eta);
            m_Refit1_eta.push_back(m_negative_eta);
            m_Refit1_phi.push_back(m_positive_phi);
            m_Refit1_phi.push_back(m_negative_phi);
            m_Refit1_d0.push_back(m_positive_d0);
            m_Refit1_d0.push_back(m_negative_d0);
            m_Refit1_z0.push_back(m_positive_z0);
            m_Refit1_z0.push_back(m_negative_z0);
        
            m_Refit1_sigma_pt.push_back(m_positive_sigma_pt);
            m_Refit1_sigma_pt.push_back(m_negative_sigma_pt);
            m_Refit1_sigma_d0.push_back(m_positive_d0_err);
            m_Refit1_sigma_d0.push_back(m_negative_d0_err);
            m_Refit1_sigma_z0.push_back(m_positive_z0_err);
            m_Refit1_sigma_z0.push_back(m_negative_z0_err);
          } // if (success_pos && success_neg) 
        } // if (m_storeZmumuNtuple) 
      } // enough refit1 tracks
      
      //fill refit2 ID parameters
      if (muonTrksRefit2->size()<2) {
          ATH_MSG_WARNING("** IDPerfMonZmumu::execute ** Refit2 muon tracks are missing!");
      } else{
        ATH_MSG_DEBUG("-- >> going to fill refit2params << --");
        success_pos = FillRecParametersSimple (refit2MuonTrk1, ppos_comb->charge(), ppos_comb_v, ctx);
        success_neg = FillRecParametersSimple (refit2MuonTrk2, pneg_comb->charge(), pneg_comb_v, ctx);
        
        if (m_storeZmumuNtuple) {
          if (success_pos && success_neg) {
            if (m_refit2Tree) { 
              ATH_MSG_DEBUG("-- Filling " << m_refit2Tree->GetName() << " entry " << m_refit2Tree->GetEntries() 
                << "  run: " << m_runNumber 
                << "  event: " << m_evtNumber 
                << "  Lumiblock: " << m_lumi_block 
                << "  Invariant mass = " << m_xZmm.GetInvMass() << " GeV ");
            }
            ATH_MSG_DEBUG("Negative p: (" << m_negative_px << ", "
              << m_negative_py << ", "
              << m_negative_pz << ") "
              << "  d0: " << m_negative_d0
              << "  z0: " << m_negative_z0
              << "  sigma_pt: " << m_positive_sigma_pt);
            ATH_MSG_DEBUG("Positive p: (" << m_positive_px << ","
              << m_positive_py << ","
              << m_positive_pz << ") "
              << "  d0: " << m_positive_d0
              << "  z0: " << m_positive_z0    
              << "  sigma_pt: " << m_negative_sigma_pt);
            if (m_refit2Tree) {
              m_refit2Tree->Fill();
            }
            // Fill vectors for common ntuple
            m_Refit2_pt.push_back(m_positive_pt);
            m_Refit2_pt.push_back(m_negative_pt);
            m_Refit2_eta.push_back(m_positive_eta);
            m_Refit2_eta.push_back(m_negative_eta);
            m_Refit2_phi.push_back(m_positive_phi);
            m_Refit2_phi.push_back(m_negative_phi);
            m_Refit2_d0.push_back(m_positive_d0);
            m_Refit2_d0.push_back(m_negative_d0);
            m_Refit2_z0.push_back(m_positive_z0);
            m_Refit2_z0.push_back(m_negative_z0);
        
            m_Refit2_sigma_pt.push_back(m_positive_sigma_pt);
            m_Refit2_sigma_pt.push_back(m_negative_sigma_pt);
            m_Refit2_sigma_d0.push_back(m_positive_d0_err);
            m_Refit2_sigma_d0.push_back(m_negative_d0_err);
            m_Refit2_sigma_z0.push_back(m_positive_z0_err);
            m_Refit2_sigma_z0.push_back(m_negative_z0_err);
          } // if (success_pos && success_neg) 
        } // if (m_storeZmumuNtuple) 
      } // enough refit2 tracks
      
      ATH_MSG_DEBUG("Execute() All NTUPLES filled  Run: " << m_runNumber << "  event: " << m_evtNumber << "  mass: " << m_xZmm.GetInvMass() << " GeV ");
    }
  } // closing -> if ( m_xZmm.AcceptEvent() )
 
  if ( !m_xZmm.AcceptEvent() ) {
    // no good muon pair found
    //failed cuts, continue to next event
    ATH_MSG_DEBUG ("** IDPerfMonZmumu::execute ** No good muon pair found. Leaving Execute(). Run: " << m_runNumber << "  event: " << m_evtNumber);
    return StatusCode::SUCCESS;
  }
 
  //
  // fill truth event iformation even when the reco event has not passed
  //
  if (m_isMC) {
    bool truthStatusIsGood = true;
    if (FillTruthParameters(ppos_comb, ctx).isFailure()){
      truthStatusIsGood = false;
      ATH_MSG_WARNING("Failed to fill truth parameters - skipping event");
    }
    if (FillTruthParameters(pneg_comb, ctx).isFailure()){
      truthStatusIsGood = false;
      ATH_MSG_WARNING("Failed to fill truth parameters - skipping event");
    }
    if (truthStatusIsGood) { 
      ATH_MSG_DEBUG ("-- Filling " << m_truthTree->GetName() << " entry " << m_truthTree->GetEntries() 
             << "  run: " << m_runNumber 
             << "  event: " << m_evtNumber 
             <<  "  Lumiblock: " << m_lumi_block 
             << "  Invariant mass = " << m_xZmm.GetInvMass() << " GeV ");
      ATH_MSG_DEBUG ("Negative p: (" << m_negative_px << ", "
             << m_negative_py << ", "
             << m_negative_pz << ") "
             << "  d0: " << m_negative_d0
             << "  z0: " << m_negative_z0
            << "  parent: " << m_negative_parent);
      ATH_MSG_DEBUG ("Positive p: (" << m_positive_px << ","
             << m_positive_py << ","
             << m_positive_pz << ") "
             << "  d0: " << m_positive_d0
             << "  z0: " << m_positive_z0
             << "  parent: " << m_positive_parent);   
      
      if (m_storeZmumuNtuple && m_truthTree) m_truthTree->Fill();
      m_Truth_pt.push_back(m_positive_pt);
      m_Truth_pt.push_back(m_negative_pt);
      m_Truth_eta.push_back(m_positive_eta);
      m_Truth_eta.push_back(m_negative_eta);
      m_Truth_phi.push_back(m_positive_phi);
      m_Truth_phi.push_back(m_negative_phi);
      m_Truth_d0.push_back(m_positive_d0);
      m_Truth_d0.push_back(m_negative_d0);
      m_Truth_z0.push_back(m_positive_z0);
      m_Truth_z0.push_back(m_negative_z0);
      m_Truth_qoverp.push_back(m_positive_qoverp);
      m_Truth_qoverp.push_back(m_negative_qoverp);
      m_Truth_parent.push_back(m_positive_parent);
      m_Truth_parent.push_back(m_negative_parent);
      
    } // truth info properly filled
    else {
      ATH_MSG_DEBUG("FAILED filling m_truthTree "
            << "  for run: " << m_runNumber 
            << "  event: " << m_evtNumber 
            << "  Lumiblock: " << m_lumi_block);
    }
  } // if (m_isMC)
 
  // fill common tree
  if (m_commonTree ) {
    // fill ntuple if some of the collections is filled
    bool dofill = false;
    if (m_IDTrack_pt.size() >= 2) dofill = true;
    if (m_CBTrack_pt.size() >= 2) dofill = true;
    if (m_Refit1_pt.size() >= 2) dofill = true;
    if (m_Refit2_pt.size() >= 2) dofill = true;
    if (m_Truth_pt.size() >= 2) dofill = true;
    
    if (dofill) {
      ATH_MSG_DEBUG("-- Filling m_commonTree " << m_commonTree->GetName() << " entry " << m_commonTree->GetEntries() 
            << "  run: " << m_runNumber 
            << "  event: " << m_evtNumber 
            << "  Lumiblock: " << m_lumi_block 
            << "  Invariant mass = " << m_xZmm.GetInvMass() << " GeV ");
      m_commonTree->Fill();
    }
    else {
      ATH_MSG_DEBUG("-- Filling m_commonTree " << m_commonTree->GetName() << " FAILED for " 
            << "  run: " << m_runNumber 
            << "  event: " << m_evtNumber 
            << "  Lumiblock: " << m_lumi_block); 
    }
  }
 
  ATH_MSG_DEBUG(" --IDPerfMonZmumu::execute--  event completed -- Run: " << m_runNumber << "  event: " << m_evtNumber);
 
  /* ***********************************************************
   * after all, let's make the collection for ALL MUON ID Tracks
   * ***********************************************************/
 
  ATH_MSG_DEBUG("IDTracks monitoring..");
  int nTracks = 0;
  TrackCollection* IDTracks = new TrackCollection(SG::OWN_ELEMENTS);
 
  // DIFFERENT RETRIEVAL
  const xAOD::MuonContainer* muons = PerfMonServices::getContainer<xAOD::MuonContainer>( PerfMonServices::MUON_COLLECTION );
  if(muons){   
    for (auto muon : *muons) {
      const xAOD::TrackParticle* tp = muon->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle);   
      if(!tp) {
    ATH_MSG_DEBUG(" - tpb = "<< tp << " -> is null? " << (!tp));
    continue;
      }
      const Trk::Track* t = tp->track();
      if ( !t ) {
    ATH_MSG_DEBUG(" - tpb does not have track ");
    continue;
      }
      const Trk::Perigee* tPerigee = t->perigeeParameters() ;
      if ( !tPerigee ) {
    ATH_MSG_DEBUG(" - t does not have perigree ");
    continue;    
      }
      IDTracks->push_back ( new Trk::Track (*t));
      nTracks++; 
    }
  }
  else  ATH_MSG_WARNING ("Could not retrieve Muon container");
 
  ATH_MSG_DEBUG(" --> nTracks = " << nTracks);
 
  StatusCode scid = evtStore()->record(IDTracks, "IDTracks");
  if (scid.isSuccess()) {
    ATH_MSG_DEBUG ("Stored "<< IDTracks->size() << " IDTracks into StoreGate");
  }
  else{
    ATH_MSG_WARNING ("Failed storing IDTracks");
  }
 
  /*********************************************************
   * ..finish
   ********************************************************/
 
  return StatusCode::SUCCESS;
}

ExtractIDHitsInformation()

void IDPerfMonZmumu::ExtractIDHitsInformation ( const xAOD::Muon * muon_pos, const xAOD::Muon * muon_neg )

private

Definition at line 2782 of file IDPerfMonZmumu.cxx.

{
  // hits info
  // positive Muon
  const xAOD::TrackParticle* IDTrkMuPos = muon_pos->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle);
  const xAOD::TrackParticle* IDTrkMuNeg = muon_neg->trackParticle(xAOD::Muon::TrackParticleType::InnerDetectorTrackParticle);
 
  int nBLhits, nhitsPIX, nhitsSCT, nhitsTRT, nContribPixLayers;
  int nPIXholes, nSCTholes;
  if (IDTrkMuPos && IDTrkMuNeg) { // make sure both summaries are available
    uint8_t dummy(-1);
    nBLhits  = IDTrkMuPos->summaryValue( dummy, xAOD::numberOfBLayerHits ) ? dummy :-1;
    nhitsPIX = IDTrkMuPos->summaryValue( dummy, xAOD::numberOfPixelHits ) ? dummy :-1;
    nhitsSCT = IDTrkMuPos->summaryValue( dummy, xAOD::numberOfSCTHits ) ? dummy :-1;
    nhitsTRT = IDTrkMuPos->summaryValue( dummy, xAOD::numberOfTRTHits ) ? dummy :-1;
    nContribPixLayers = IDTrkMuPos->summaryValue( dummy, xAOD::numberOfContribPixelLayers )? dummy :-1;
 
    nPIXholes = IDTrkMuPos->summaryValue( dummy, xAOD::numberOfPixelHoles )? dummy :-1;
    nSCTholes = IDTrkMuPos->summaryValue( dummy, xAOD::numberOfSCTHoles )? dummy :-1;
 
 
    m_nBLhits.push_back(nBLhits);
    m_nPIXhits.push_back(nhitsPIX);
    m_nSCThits.push_back(nhitsSCT);
    m_nTRThits.push_back(nhitsTRT);
 
    ATH_MSG_DEBUG (" -- mu_pos -- HITS --"
           << "  nBLhits:  " << nBLhits
           << "  nhitsPIX: " << nhitsPIX
           << "  nPIXLayers: " << nContribPixLayers
           << "  nhitsSCT: " << nhitsSCT
           << "  Silicon holes: " << nPIXholes + nSCTholes
           << "    nhitsTRT: " << nhitsTRT);
 
    nBLhits  = IDTrkMuNeg->summaryValue( dummy, xAOD::numberOfBLayerHits ) ? dummy :-1;
    nhitsPIX = IDTrkMuNeg->summaryValue( dummy, xAOD::numberOfPixelHits ) ? dummy :-1;
    nhitsSCT = IDTrkMuNeg->summaryValue( dummy, xAOD::numberOfSCTHits ) ? dummy :-1;
    nhitsTRT = IDTrkMuNeg->summaryValue( dummy, xAOD::numberOfTRTHits ) ? dummy :-1;
    nContribPixLayers = IDTrkMuNeg->summaryValue( dummy, xAOD::numberOfContribPixelLayers )? dummy :-1;
    nPIXholes = IDTrkMuNeg->summaryValue( dummy, xAOD::numberOfPixelHoles )? dummy :-1;
    nSCTholes = IDTrkMuNeg->summaryValue( dummy, xAOD::numberOfSCTHoles )? dummy :-1;
 
    m_nBLhits.push_back(nBLhits);
    m_nPIXhits.push_back(nhitsPIX);
    m_nSCThits.push_back(nhitsSCT);
    m_nTRThits.push_back(nhitsTRT);
 
    ATH_MSG_DEBUG (" -- mu_neg -- HITS --"
           << "  nBLhits:  " << nBLhits
           << "  nhitsPIX: " << nhitsPIX
           << "  nPIXLayers: " << nContribPixLayers
           << "  nhitsSCT: " << nhitsSCT
           << "  Silicon holes: " << nPIXholes + nSCTholes
           << "    nhitsTRT: " << nhitsTRT);
 
  }
  return;
}

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

extraOutputDeps()

const DataObjIDColl & AthCommonAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 89 of file AthCommonAlgorithm.cxx.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return BaseAlg::extraOutputDeps();
}

FillRecParameters()

StatusCode IDPerfMonZmumu::FillRecParameters ( const Trk::Track * track, const xAOD::TrackParticle * trackp_for_unbias, double charge, const xAOD::Vertex * vertex, const EventContext & ctx )

private

Definition at line 1803 of file IDPerfMonZmumu.cxx.

{
  if (!track){
    ATH_MSG_DEBUG("* FillRecParameters * Empty Track: track. Skipping.");
    return StatusCode::FAILURE;
  }
 
  if (m_doRemoval && !trackp_for_unbias && !m_skipMS) {
    ATH_MSG_DEBUG("* FillRecParameters * Empty Track: trackp_for_unbias. Skipping.");
    return StatusCode::FAILURE;
  }
 
  const Trk::Perigee* trkPerigee = track->perigeeParameters();
  const Trk::Perigee* trk_for_unbiasPerigee = nullptr;
  if (!m_skipMS)  trk_for_unbiasPerigee = &(trackp_for_unbias->perigeeParameters());
 
  double px = 0;
  double py = 0;
  double pt = 0;
  double pz = 0;
  double phi= 0;
  double eta= 0;
  double d0 = 0;
  double z0 = 0;
  
  double d0_err = 999.;
  double z0_err = 999.;
  
  double PVd0 = 0;
  double PVz0 = 0;
  double PVd0res = 0;
  double PVz0res = 0;
  
 
  if (trkPerigee != nullptr){
    double qOverP   = trkPerigee->parameters()[Trk::qOverP];
    if (qOverP) {
      px = trkPerigee->momentum().x();
      py = trkPerigee->momentum().y();
      pt = std::abs(trkPerigee->pT()); 
      pz = trkPerigee->momentum().z();
      phi= trkPerigee->parameters()[Trk::phi];
      eta= trkPerigee->eta();
      d0 = trkPerigee->parameters()[Trk::d0];
      z0 = trkPerigee->parameters()[Trk::z0];
      d0_err = Amg::error(*trkPerigee->covariance(),Trk::d0);
      z0_err = Amg::error(*trkPerigee->covariance(),Trk::z0);
    }
  }
 
  // access beam spot and extrapolate track till there
  SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandleRec { m_beamSpotKey, ctx };
  auto beamline = m_trackToVertexTool->GetBeamLine(beamSpotHandleRec.cptr());
 
  auto trackBLtemp = m_trackToVertexTool->trackAtBeamline(ctx, *track , beamline.get() );
  const Trk::AtaStraightLine*  atBL = dynamic_cast<const Trk::AtaStraightLine*>(trackBLtemp.get());
 
  if (atBL){
    double qOverP   = atBL->parameters()[Trk::qOverP];
    if(qOverP){
      px = atBL->momentum().x();
      py = atBL->momentum().y();
      pt = std::abs(atBL->pT());
      pz = atBL->momentum().z();
      eta= trkPerigee->eta();
      phi= trkPerigee->parameters()[Trk::phi];
      d0 = atBL->parameters()[Trk::d0];
      z0 = atBL->parameters()[Trk::z0];
    } 
    delete atBL;
  }
  else {
    ATH_MSG_WARNING("FillRecParameters::Failed extrapolation to the BeamLine");
  }
 
  
  // Amg::Vector3D position = beamSpotHandleRec->beamPos();
  TLorentzVector vtrack = TLorentzVector (trkPerigee->momentum().x(),
                      trkPerigee->momentum().y(),
                      trkPerigee->momentum().z(),
                      trkPerigee->momentum().mag());
  
  if(m_doIP && vertex){ //I assume that the vertex is the same of the original track
    std::unique_ptr<Trk::ImpactParametersAndSigma> iPandSigma(nullptr);
    if (!m_skipMS) iPandSigma = m_trackToVertexIPEstimator->estimate(ctx, trk_for_unbiasPerigee,trkPerigee,vertex,m_doRemoval);
    else iPandSigma = m_trackToVertexIPEstimator->estimate(ctx, trkPerigee,vertex);
    if( iPandSigma==0 ){
      ATH_MSG_WARNING ("FillRecParameters::trackToVertexIPEstimator failed !");
      return StatusCode::FAILURE;
    }
    else{
      ATH_MSG_DEBUG("FillRecParameters::trackToVertexIPEstimator success !");
      PVd0 = iPandSigma->IPd0;
      PVd0res = iPandSigma->PVsigmad0;
      //d0res = iPandSigma->sigmad0;  //-> ? 
      PVz0 = iPandSigma->IPz0;
      PVz0res = iPandSigma->PVsigmaz0;
      //z0res = iPandSigma->sigmaz0;  //-> ?
      
      if (vertex) {
    ATH_MSG_INFO(" FillRecParameters::trackToVertexIPEstimator vertex does exist ");
    m_pv_x = vertex->x();
    m_pv_y = vertex->y();
    m_pv_z = vertex->z();
    m_nTrkInVtx = vertex->nTrackParticles();
      }
      
    }
  }
    
  if (charge == 1) {
    m_positive_px = px;
    m_positive_py = py;
    m_positive_pt = pt;
    m_positive_pz = pz;
    m_positive_phi= phi;
    m_positive_eta= eta;
    m_positive_z0 = z0;
    m_positive_d0 = d0;
    m_positive_d0_err = d0_err;
    m_positive_z0_err = z0_err;
    if(m_doIP){
      m_positive_z0_PV    = PVz0;
      m_positive_d0_PV    = PVd0;
      m_positive_z0_PVerr = PVz0res;
      m_positive_d0_PVerr = PVd0res;
    }
    ATH_MSG_DEBUG("(Filled charge == 1 ) (reco)-> px : "<< px <<" py: "<<py <<" pz: "<<pz <<" d0: "<<d0<<" z0: "<< z0);
  }
  if (charge == -1) {
    m_negative_px = px;
    m_negative_py = py;
    m_negative_pt = pt;
    m_negative_pz = pz;
    m_negative_phi= phi;
    m_negative_eta= eta;
    m_negative_z0 = z0;
    m_negative_d0 = d0;
    m_negative_d0_err = d0_err;
    m_negative_z0_err = z0_err;
    if(m_doIP){
      m_negative_z0_PV = PVz0;
      m_negative_d0_PV = PVd0;
      m_negative_z0_PVerr = PVz0res;
      m_negative_d0_PVerr = PVd0res;
    }
    ATH_MSG_DEBUG("(Filled charge == -1 ) (reco)-> px : "<< px <<" py: "<<py <<" pz: "<<pz <<" d0: "<<d0<<" z0: "<<z0 );
  }
 
  return StatusCode::SUCCESS;
}

FillRecParametersSimple()

StatusCode IDPerfMonZmumu::FillRecParametersSimple ( const Trk::Track * track, float charge, const xAOD::Vertex * vertex, const EventContext & ctx )

private

Definition at line 1955 of file IDPerfMonZmumu.cxx.

{
  if (!track){
    ATH_MSG_DEBUG("* FillRecParametersSimple * Empty Track: track. Skipping.");
    return StatusCode::FAILURE;
  }
 
  const Trk::Perigee* trkPerigee = track->perigeeParameters();
  if (not trkPerigee) {
    ATH_MSG_ERROR("trkPerigee pointer is null in IDPerfMonZmumu::FillRecParametersSimple");
    return StatusCode::FAILURE;
  }
  double px = 0;
  double py = 0;
  double pt = 0;
  double pz = 0;
  double phi= 0;
  double eta= 0;
 
  double d0 = 0;
  double z0 = 0;  
  double d0_err = 999.;
  double z0_err = 999.;
  
  //
  double qOverP   = trkPerigee->parameters()[Trk::qOverP];
  if (qOverP) {
    px = trkPerigee->momentum().x();
    py = trkPerigee->momentum().y();
    pt = std::abs(trkPerigee->pT()); 
    pz = trkPerigee->momentum().z();
    phi= trkPerigee->parameters()[Trk::phi];
    eta= trkPerigee->eta();
    d0 = trkPerigee->parameters()[Trk::d0];
    z0 = trkPerigee->parameters()[Trk::z0];
    d0_err = Amg::error(*trkPerigee->covariance(),Trk::d0);
    z0_err = Amg::error(*trkPerigee->covariance(),Trk::z0);
  }
  
  
  SG::ReadHandle<xAOD::EventInfo> eventInfo (m_EventInfoKey, ctx);
  Amg::Vector3D position (eventInfo->beamPosX(), eventInfo->beamPosY(), eventInfo->beamPosZ());
  TLorentzVector vtrack = TLorentzVector (trkPerigee->momentum().x(),
                      trkPerigee->momentum().y(),
                      trkPerigee->momentum().z(),
                      trkPerigee->momentum().mag());
 
  SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandleRec { m_beamSpotKey, ctx }; // This method loads the proper beam spot conditions
 
 
  // access beam spot and extrapolate track till there
  auto beamline = m_trackToVertexTool->GetBeamLine(beamSpotHandleRec.cptr());
 
  auto trackBLtemp = m_trackToVertexTool->trackAtBeamline(ctx, *track , beamline.get() );
  const Trk::AtaStraightLine*  atBL = dynamic_cast<const Trk::AtaStraightLine*>(trackBLtemp.get());
 
  if (atBL){
    double qOverP   = atBL->parameters()[Trk::qOverP];
    if(qOverP){
      px = atBL->momentum().x();
      py = atBL->momentum().y();
      pt = std::abs(atBL->pT());
      pz = atBL->momentum().z();
      eta= trkPerigee->eta();
      phi= trkPerigee->parameters()[Trk::phi];
      d0 = atBL->parameters()[Trk::d0];
      z0 = atBL->parameters()[Trk::z0] - vertex->z(); // correct also for the vertex position along the beam line
    } 
    ATH_MSG_DEBUG("* FillRecParametersSimple *" 
          << " charge " << charge
          << " ** trackToVertexTool->trackAtBeamline ** " 
          << "( " << beamSpotHandleRec->beamPos().x()
          << ", " << beamSpotHandleRec->beamPos().y()
          << ", " << beamSpotHandleRec->beamPos().z()
          << ")  pt: " << pt
          << "  d0: " << d0
          << "  z0: " << z0);
  }
  else {
    ATH_MSG_WARNING("FillRecParametersSimple::Failed extrapolation to the BeamLine --> extrapolate using line approximation");
    float trkd0 =  trkPerigee->parameters()[Trk::d0];
    float trkz0 =  trkPerigee->parameters()[Trk::z0];
    float bsX =    beamSpotHandleRec->beamPos().x();
    float bsY =    beamSpotHandleRec->beamPos().y();
    float bsZ =    beamSpotHandleRec->beamPos().z();
    float btiltX = beamSpotHandleRec->beamTilt(0);
    float btiltY = beamSpotHandleRec->beamTilt(1);
    // compute the corrections for the track parameters considering the beamspot position
    float beamX = bsX + std::tan(btiltX) * (trkz0-bsZ);
    float beamY = bsY + std::tan(btiltY) * (trkz0-bsZ);
    float beamD0 = ( -std::sin(vtrack.Phi())*beamX + std::cos(vtrack.Phi())*beamY );
 
    // correct d0 and z0
    d0 = trkd0 - beamD0;
    z0 = trkz0 - bsZ - vertex->z(); // correct also for the vertex position along the beam line  
 
    ATH_MSG_DEBUG("* FillRecParametersSimple *" 
          << " charge " << charge
          << " ** beamSpotHandleRec " << m_beamSpotKey 
          << "( " << beamSpotHandleRec->beamPos().x()
          << ", " << beamSpotHandleRec->beamPos().y()
          << ", " << beamSpotHandleRec->beamPos().z()
          << ")  pt: " << pt
          << "  d0: " << d0 
          << "  z0: " << z0);
  }
  
  // store the values 
  if (charge == 1) { // positive particle
    m_positive_px = px;
    m_positive_py = py;
    m_positive_pt = pt;
    m_positive_pz = pz;
    m_positive_phi= phi;
    m_positive_eta= eta;
    m_positive_z0 = z0;
    m_positive_d0 = d0;
    m_positive_d0_err = d0_err;
    m_positive_z0_err = z0_err;
  }
  if (charge == -1) { // negative particle
    m_negative_px = px;
    m_negative_py = py;
    m_negative_pt = pt;
    m_negative_pz = pz;
    m_negative_phi= phi;
    m_negative_eta= eta;
    m_negative_z0 = z0;
    m_negative_d0 = d0;
    m_negative_d0_err = d0_err;
    m_negative_z0_err = z0_err;
  }
 
  ATH_MSG_DEBUG("* FillRecParametersSimple * completed * charge " << charge << "  pt: " << pt << "  d0: " << d0 << "  z0: " << z0);
  return StatusCode::SUCCESS;
}

FillRecParametersTP()

StatusCode IDPerfMonZmumu::FillRecParametersTP ( const xAOD::TrackParticle * trackp, const xAOD::TrackParticle * trackp_for_unbias, double charge, const xAOD::Vertex * vertex, const EventContext & ctx )

private

Definition at line 1661 of file IDPerfMonZmumu.cxx.

{
  if (!trackp){
    //ATH_MSG_WARNING("Empty Trackparticle. Skipping.");
    return StatusCode::FAILURE;
  }
 
  if (m_doRemoval && !trackp_for_unbias && !m_skipMS){
    return StatusCode::FAILURE;
  }
    
  double px = 0;
  double py = 0;
  double pt = 0;
  double pz = 0;
  double phi= 0;
  double eta= 0;
  double d0 = 0;
  double z0 = 0;
  double d0res = 0;
  double z0res = 0;
  double PVd0res = 0;
  double PVz0res = 0;
  double PVd0 = 0;
  double PVz0 = 0;
  double sigma_pt = 0;
  double qoverp = 0;
  double sigma_qOverP = 1.; // use a large value by default as 0 is meaningless
 
 
  px = trackp->p4().Px(); 
  py = trackp->p4().Py();
  pt = trackp->p4().Pt(); 
  pz = trackp->p4().Pz();
  phi= trackp->p4().Phi();
  eta= trackp->p4().Eta();
  d0 = trackp->d0();
  z0 = trackp->z0();
  qoverp = trackp->qOverP();
 
  d0res = std::sqrt(trackp->definingParametersCovMatrix()(Trk::d0,Trk::d0));
  z0res = std::sqrt(trackp->definingParametersCovMatrix()(Trk::z0,Trk::z0));
  sigma_qOverP = std::sqrt(trackp->definingParametersCovMatrix()(Trk::qOverP,Trk::qOverP));
  double sigma_theta  = std::sqrt(trackp->definingParametersCovMatrix()(Trk::theta,Trk::theta));
    
  // computing sigma_pt
  // pt = sin(theta) / qOverP --> sigma(pt) = (sin(theta)/ qOverP^2) x sigma(qOverP) // neglecting sigma(sin(theta))
  double sigma_pt_term1 = (sin(trackp->theta()) / pow(qoverp,2)) * sigma_qOverP;
  double sigma_pt_term2 = (1./qoverp) * cos(trackp->theta()) * sigma_theta;
  double sigma_pt_term3 = (cos(trackp->theta()) / pow(qoverp,2)) * sigma_theta * sigma_qOverP;
  sigma_pt = sqrt( pow(sigma_pt_term1,2) + pow(sigma_pt_term2,2) + 2 * sigma_pt_term3 * trackp->definingParametersCovMatrix()(Trk::qOverP,Trk::theta));
 
  //
  if (vertex == nullptr) {
    ATH_MSG_WARNING("in FillRecParametersTP. WARNING: Vertex is NULL");
    return StatusCode::FAILURE;
  }
  if (!m_trackToVertexIPEstimator) {
    ATH_MSG_WARNING("in FillRecParametersTP. WARNING: m_trackToVertexIPEstimator is NULL");
    return StatusCode::FAILURE;
  }
 
  if (m_doIP && vertex != nullptr && m_trackToVertexIPEstimator){
    std::unique_ptr<Trk::ImpactParametersAndSigma> iPandSigma(nullptr);
    ATH_MSG_DEBUG("-- using the trackToVertexIPEstimator --");
    //Calling the estimate(trackp,newtrackp,vertex,doRemoval)
    //The first track is used to unbias the vertex, the second to get the extrapolation
    if (!m_skipMS) iPandSigma = m_trackToVertexIPEstimator->estimate(ctx, trackp_for_unbias, trackp, vertex, m_doRemoval);
    else iPandSigma = m_trackToVertexIPEstimator->estimate(ctx, trackp, vertex);
    ATH_MSG_DEBUG("return from the trackToVertexIPEstimator->estimate()");
    
    if( !iPandSigma ){
      ATH_MSG_WARNING ("FillRecParametersTP::trackToVertexIPEstimator failed !");
      return StatusCode::FAILURE;
    }
    else{
      ATH_MSG_DEBUG("FillRecParametersTP::trackToVertexIPEstimator success !");
      PVd0 = iPandSigma->IPd0;
      PVd0res = iPandSigma->PVsigmad0;
      d0res = iPandSigma->sigmad0;
      PVz0 = iPandSigma->IPz0;
      PVz0res = iPandSigma->PVsigmaz0;
      z0res = iPandSigma->sigmaz0;
    }
  }
 
  if (charge == 1) {
    m_positive_px = px;
    m_positive_py = py;
    m_positive_pt = pt;
    m_positive_pz = pz;
    m_positive_phi= phi;
    m_positive_eta= eta;
    m_positive_z0 = z0;
    m_positive_z0_err = z0res;
    m_positive_d0 = d0;
    m_positive_d0_err = d0res;
    m_positive_sigma_pt = sigma_pt;
    m_positive_qoverp = qoverp;
    m_positive_sigma_qoverp = sigma_qOverP;
    if(m_doIP){
      m_positive_z0_PV = PVz0;
      m_positive_d0_PV = PVd0;
      m_positive_z0_PVerr = PVz0res;
      m_positive_d0_PVerr = PVd0res;
    }
    ATH_MSG_DEBUG("(Filled charge == 1 ) (reco)-> px : "<< px <<" py: "<<py <<" pt: "<<pt <<" pz: "<<pz 
          <<"   d0: "<<m_positive_d0 << " d0res : "<< d0res << " PVd0res : "<< PVd0res <<" z0: "<< m_positive_z0 << " z0res : " << z0res <<  " PVz0res : "<< PVz0res );
 
  } else  if (charge == -1) {
    m_negative_px = px;
    m_negative_py = py;
    m_negative_pt = pt;
    m_negative_pz = pz;
    m_negative_phi= phi;
    m_negative_eta= eta;
    m_negative_z0 = z0;
    m_negative_z0_err = z0res;
    m_negative_d0 = d0;
    m_negative_d0_err = d0res;
    m_negative_sigma_pt = sigma_pt;
    m_negative_qoverp = qoverp;
    m_negative_sigma_qoverp = sigma_qOverP;
    if(m_doIP){
      m_negative_z0_PV = PVz0;
      m_negative_d0_PV = PVd0;
      m_negative_z0_PVerr = PVz0res;
      m_negative_d0_PVerr = PVd0res;
    }
    ATH_MSG_DEBUG("(Filled charge == -1 ) (reco)-> px : "<< px <<" py: "<< py <<" pt: " << pt <<" pz: "<<pz 
          <<"   d0: "<<m_negative_d0 << " d0res : "<< d0res << " PVd0res : "<< PVd0res <<" z0: "<< m_negative_z0 << " z0res : " << z0res <<  " PVz0res : "<< PVz0res );
  }
 
  return StatusCode::SUCCESS;
}

FillTruthParameters()

StatusCode IDPerfMonZmumu::FillTruthParameters ( const xAOD::TrackParticle * track, const EventContext & ctx )

private

Definition at line 2093 of file IDPerfMonZmumu.cxx.

{
 
  if (!trackParticle ){//|| !trackParticle->vertex()){
    ATH_MSG_WARNING("-- FillTruthParameters -- Empty Trackparticle. Skipping.");
    return StatusCode::FAILURE;
  }
 
  const xAOD::TruthParticle* particle = getTruthParticle( *trackParticle );
  double charge = 0;
  if (!particle) {
    ATH_MSG_DEBUG(" -- FillTruthParameters -- Failure while retrieving the truth particle. Returning FAILURE.");
    return StatusCode::FAILURE;
  }
 
  if ( particle->pdgId() ==  13) charge = -1.; // muon-
  if ( particle->pdgId() == -13) charge = 1.; // muon+
  if ( particle->isNeutral() ) {
    ATH_MSG_DEBUG(" -- FillTruthParameters -- reco muon associated to a truth neutral!! Returning FAILURE.");
    return StatusCode::FAILURE;
  }
  if ( charge == 0 ) {
    ATH_MSG_DEBUG(" -- FillTruthParameters -- reco muon associated to a non true muon!! Returning FAILURE.");
    return StatusCode::FAILURE;
  } 
 
  const Amg::Vector3D momentum(particle->px(), particle->py(), particle->pz());  
  const xAOD::TruthVertex * ptruthVertex(0);
  ptruthVertex=particle->prodVtx();
  if (!ptruthVertex){
    ATH_MSG_DEBUG("A production vertex pointer was retrieved, but it is NULL");
    return StatusCode::FAILURE;
  }
 
 
  if (particle->parent(0) != nullptr) {
    static const SG::ConstAccessor<int> truthTypeAcc("truthType");
    static const SG::ConstAccessor<int> truthOriginAcc("truthOrigin");
    if (truthTypeAcc.isAvailable(*particle)) {
      ATH_MSG_DEBUG(" -- FillTruthParameters -- truth particle is good. pdg: " << particle->pdgId() 
            << "  type: " << particle->type()
            << "  nparents " << particle->nParents()
            << "  parent->pdg: " << particle->parent(0)->pdgId()
            << "  truthType: " << truthTypeAcc(*particle)
            << "  truthOrigin: " << truthOriginAcc(*particle)
            );
 
    }
    else {
      ATH_MSG_DEBUG(" -- FillTruthParameters -- truth particle is good. pdg: " << particle->pdgId() 
            << "  type: " << particle->type()
            << "  nparents " << particle->nParents()
            << "  parent->pdg: " << particle->parent(0)->pdgId()
            << "  truthType & truthOrigin: " << " NOT AVAILABLE " 
            );
    }
  }
 
  const auto xPos=ptruthVertex->x();
  const auto yPos=ptruthVertex->y();
  const auto z_truth=ptruthVertex->z();
 
  const Amg::Vector3D position(xPos, yPos, z_truth);
  const Trk::CurvilinearParameters cParameters(position, momentum, charge);
 
  SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandleTruth { m_beamSpotKey, ctx };
  Trk::PerigeeSurface persf( beamSpotHandleTruth->beamPos() );
 
  std::unique_ptr<const Trk::TrackParameters> tP = m_extrapolator->extrapolate(ctx, cParameters,persf, Trk::anyDirection, false);
 
  double px = 0;
  double py = 0;
  double pt = 0;
  double pz = 0;
  double phi= 0;
  double eta= 0;
  double d0 = 0;
  double d0res = 0;
  double PVd0res = 0;
  double z0 = 0;
  double z0res = 0;
  double PVz0res = 0;
 
  double d0recoPos =  m_positive_d0;
  double z0recoPos =  m_positive_z0;
  double d0recoNeg =  m_negative_d0;
  double z0recoNeg =  m_negative_z0;
 
  ATH_MSG_DEBUG("reco IPs (pos): > d0 : "<<d0recoPos << " z0: " << z0recoPos << " trackp z0 : " << trackParticle->z0() << " trackp d0 : " << trackParticle->d0());
  ATH_MSG_DEBUG("reco IPs (neg): > d0 : "<<d0recoNeg << " z0: " << z0recoNeg << " trackp z0 : " << trackParticle->z0() <<" trackp d0 : " << trackParticle->d0() );
 
  double qOverP_truth = 0.;
  if (tP){
    qOverP_truth = tP->parameters()[Trk::qOverP];
    if( qOverP_truth ){
 
      px = tP->momentum().x();
      py = tP->momentum().y();
      pt = tP->pT();
      pz = tP->momentum().z();
      phi= tP->parameters()[Trk::phi];
      eta= tP->eta();
      d0 = tP->parameters()[Trk::d0];
      z0 = tP->parameters()[Trk::z0];
 
      ATH_MSG_DEBUG("cand perig HEP particle (truth) px : "<< tP->momentum().x());
      ATH_MSG_DEBUG("cand perig HEP particle (truth) py : "<< tP->momentum().y());
      ATH_MSG_DEBUG("cand perig HEP particle (truth) pz : "<< tP->momentum().z());
      ATH_MSG_DEBUG("cand perig HEP particle (truth) d0 : "<< tP->parameters()[Trk::d0]);
      ATH_MSG_DEBUG("cand perig HEP particle (truth) z0 : "<< tP->parameters()[Trk::z0]);
      
    }
  }
 
 
  if (charge == 1) {
    m_positive_px = px;
    m_positive_py = py;
    m_positive_pt = pt;
    m_positive_pz = pz;
    m_positive_phi= phi;
    m_positive_eta= eta;
    m_positive_z0 = z0recoPos -z0;
    m_positive_z0_err = z0res;
    m_positive_z0_PVerr = PVz0res;
    m_positive_d0 = d0recoPos -d0;
    m_positive_d0_err = d0res;
    m_positive_d0_PVerr = PVd0res;
    m_positive_qoverp = qOverP_truth;
 
    bool parentfound = false;
    if (particle->nParents()>0) {
      if (particle->parent(0) != nullptr) {
    m_positive_parent = particle->parent(0)->pdgId();
    parentfound = true;
      } 
    }
    if (!parentfound) m_positive_parent = 0;
 
    ATH_MSG_DEBUG(" -- FillTruthParameters -- charge = 1  (truth)-> px : "<< m_positive_px 
          <<"  py: "<<m_positive_py 
          <<"  pz: "<<m_positive_pz 
          <<"  d0: "<<m_positive_d0 
          <<"  z0: "<< m_positive_z0
          <<"  qoverp: " << m_positive_qoverp
         << "  parent: " << m_positive_parent);
 
  } else  if (charge == -1) {
    m_negative_px = px;
    m_negative_py = py;
    m_negative_pt = pt;
    m_negative_pz = pz;
    m_negative_phi= phi;
    m_negative_eta= eta;
    m_negative_z0 = z0recoNeg-z0;
    m_negative_z0_err = z0res;
    m_negative_z0_PVerr = PVz0res;
    m_negative_d0 = d0recoNeg-d0;
    m_negative_d0_err = d0res;
    m_negative_d0_PVerr = PVd0res;
    m_negative_qoverp = qOverP_truth;
 
    bool parentfound = false;
    if (particle->nParents()>0) {
      if (particle->parent(0) != nullptr) {
    m_negative_parent = particle->parent(0)->pdgId();
    parentfound = true;
      }
    }
    if (!parentfound) m_negative_parent = 0;
 
    ATH_MSG_DEBUG(" -- FillTruthParameters-- charge = -1 (truth)-> px : "<< m_negative_px 
          << "  py: "<<m_negative_py 
          << "  pz: "<<m_negative_pz 
          << "  d0: "<<m_negative_d0 
          << "  z0: "<< m_negative_z0
          <<"  qoverp: " << m_negative_qoverp
          << "  parent:" << m_negative_parent);
  }
  return StatusCode::SUCCESS;
}

filterPassed()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Get filter decision:

Definition at line 93 of file AthCommonAlgorithm.h.

                                                           {
    return execState( ctx ).filterPassed();
  }

finalize()

StatusCode IDPerfMonZmumu::finalize ( )

virtual

Definition at line 2297 of file IDPerfMonZmumu.cxx.

{
  m_xZmm.finalize();
  return StatusCode::SUCCESS;
}

GetDiMuonVertex()

const xAOD::Vertex * IDPerfMonZmumu::GetDiMuonVertex ( const xAOD::TrackParticle * muon1, const xAOD::TrackParticle * muon2 )

private

Definition at line 2365 of file IDPerfMonZmumu.cxx.

                                                                                                                {
  
  const xAOD::VertexContainer* vxContainer(0);
  const xAOD::Vertex* myVtx(0);
  vxContainer = PerfMonServices::getContainer<xAOD::VertexContainer> (PerfMonServices::VTX_COLLECTION);
  if (!vxContainer) { 
    return myVtx;
  }
  
  for (int ivtx=0; ivtx< (int) vxContainer->size();ivtx++)  {
    myVtx = (*vxContainer)[ivtx];
    if ((myVtx->vertexType() == xAOD::VxType::PriVtx)) {
      if (m_Trk2VtxAssociationTool->isCompatible(*muon1,*myVtx) && (m_Trk2VtxAssociationTool->isCompatible(*muon2,*myVtx)))
    return myVtx;
    }
  }
  return myVtx;
}

GetMuonQualityValue()

int IDPerfMonZmumu::GetMuonQualityValue ( std::string qualityname )

private

getTruthParticle()

const xAOD::TruthParticle * IDPerfMonZmumu::getTruthParticle ( const xAOD::IParticle & p )

private

A convenience type declaration

A static accessor for the information

Definition at line 2277 of file IDPerfMonZmumu.cxx.

                                                                                  {
  typedef ElementLink< xAOD::TruthParticleContainer > Link_t;
  static const SG::ConstAccessor< Link_t > acc( "truthParticleLink" );
  // Check if such a link exists on the object:
  if( ! acc.isAvailable( p ) ) {
    return nullptr;
  }
  // Get the link:
  const Link_t& link = acc( p );
  // Check if the link is valid:
  if( ! link.isValid() ) {
    return nullptr;
  }
  // Everything has passed, let's return the pointer:
  return *link;
}

initialize()

StatusCode IDPerfMonZmumu::initialize ( )

virtual

Definition at line 170 of file IDPerfMonZmumu.cxx.

{
  ATH_MSG_DEBUG ("** IDPerfMonZmumu::Initialize ** START **");
  // Setup the services
 
  // Retrieve Track fitter and track to vertex
  if (m_doRefit) { // only if track refit is requested
    if (m_TrackRefitter1.retrieve().isSuccess()) {
      ATH_MSG_INFO("Retrieved tool m_TrackRefitter1: " << m_TrackRefitter1 << " SUCCESS ");
    } 
    else {
      ATH_MSG_FATAL("Unable to retrieve m_TrackRefitter1 " << m_TrackRefitter1 << " FAILURE ");
      return StatusCode::FAILURE;
    }
 
    // Retrieve the second fitter
    if (m_TrackRefitter2.retrieve().isSuccess()) {
      ATH_MSG_INFO("Retrieved tool m_TrackRefitter2: " << m_TrackRefitter2 << " SUCCESS ");
    } 
    else {
      ATH_MSG_FATAL("Unable to retrieve m_TrackRefitter2 " << m_TrackRefitter2 << " FAILURE ");
      return StatusCode::FAILURE;
    }
 
    if (m_trackToVertexTool.retrieve().isSuccess()) {
      ATH_MSG_INFO("Retrieved tool m_trackToVertexTool " << m_trackToVertexTool << " SUCCESS ");
    } 
    else {
      ATH_MSG_FATAL("Unable to retrieve m_trackToVertexTool " << m_trackToVertexTool << " FAILURE ");
      return StatusCode::FAILURE;
    }
  }
 
  if(m_useTrackSelectionTool){
    //    ATH_CHECK( m_selTool.retrieve() );
    if (m_selTool.retrieve().isSuccess()){
      ATH_MSG_INFO("Retrieved tool (track selection tool) m_selTool " << m_selTool << " SUCCESS ");
    }
    else {
      ATH_MSG_FATAL("Unable to retrieve (track selection tool) m_selTool " << m_selTool << " FAILURE ");
      return StatusCode::FAILURE;
    }
  }
 
  // initializing the eventInfo "accessor"
  ATH_CHECK(  m_EventInfoKey.initialize() );
 
  if (m_doIP) {
    ATH_MSG_DEBUG("Retrieving tool (trackToVertexIPEstimator)");
    ATH_CHECK (m_trackToVertexIPEstimator.retrieve());
  }
 
 
  ATH_CHECK (m_EventInfoKey.initialize());   // initializing the eventInfo "accessor"
 
  if (m_isMC) ATH_CHECK (m_extrapolator.retrieve()); // this is only used for the truth particles
  
  ATH_CHECK (m_vertexKey.initialize());
  
  ATH_CHECK (m_trackContainerName.initialize());
 
  ATH_MSG_INFO(" -- IDPerfMonZmumu::initialize() -- m_vertexKey: " << m_vertexKey);
 
  ATH_MSG_INFO(" -- IDPerfMonZmumu::initialize() -- init m_beamSpotKey ");
  ATH_CHECK(m_beamSpotKey.initialize());
 
  // START new place for initilization of params
  m_xZmm.setContainer       (PerfMonServices::MUON_COLLECTION);
  m_xZmm.doIsoSelection     (m_doIsoSelection);
  m_xZmm.doIPSelection      (m_doIPSelection);
  m_xZmm.doMCPSelection     (m_doMCPSelection);
  m_xZmm.SetMassWindowLow   (m_MassWindowLow);
  m_xZmm.SetMassWindowHigh  (m_MassWindowHigh);
  m_xZmm.SetLeadingMuonPtCut(m_LeadingMuonPtCut);
  m_xZmm.SetSecondMuonPtCut (m_SecondMuonPtCut);
  m_xZmm.SetOpeningAngleCut (m_OpeningAngleCut);
  m_xZmm.SetZ0GapCut        (m_Z0GapCut);
  m_xZmm.setDebugMode       (m_doDebug);
  m_xZmm.SetMuonQuality     (m_MuonQualityName);
  m_xZmm.SetSkipMSCheck     (m_skipMS);
  m_xZmm.SetMinLumiBlock    (m_minGoodLumiBlock);
  m_xZmm.SetMaxLumiBlock    (m_maxGoodLumiBlock);
  
  if (m_useCustomMuonSelector) {
    m_xZmm.SetMuonSelectionTool (m_muonSelector);    
  }
  // END new place for initialization of params
 
  if (m_doFourMuAnalysis) {
    m_4mu.Init();
    m_4mu.setDebugMode(true);
    if (m_useCustomMuonSelector) m_4mu.SetMuonSelectionTool (m_muonSelector);
    ATH_MSG_DEBUG(" IDPerfMonZmumu FourMuonAnalysis initialization completed " <<  m_Trk2VtxAssociationTool);
  }
 
  // m_Trk2VtxAssociationTool = std::make_unique<CP::TrackVertexAssociationTool>("Loose");
  if ( m_Trk2VtxAssociationTool.retrieve().isSuccess() )  {
    ATH_MSG_DEBUG(" IDPerfMonZmumu Success retrieving tool " <<  m_Trk2VtxAssociationTool);
  }
  else { // is Failure
    ATH_MSG_FATAL("IDPerfMonZmumu Failed to retrieve tool m_Trk2VtxAssociationTool " <<  m_Trk2VtxAssociationTool);
    return StatusCode::FAILURE;
  }
 
  ATH_CHECK ( this->bookTrees() );
 
  if (m_UseTrigger) { // load trigger decission and matching under user request
    if (m_triggerDecision.retrieve().isFailure()) {                                                                                                                                   
      ATH_MSG_FATAL("Unable to retrieve " << m_triggerDecision << " turn it off");                                                                                             
      return StatusCode::FAILURE;                                                                                                                                               
    }                                                                                                                                                                           
    else {                                                                                                                                                                        
      ATH_MSG_DEBUG("retrieved tool: " << m_triggerDecision );                                                                                                                  
    }
 
    if(m_triggerMatching.retrieve().isFailure()) {       
      ATH_MSG_FATAL("Unable to retrieve " << m_triggerMatching << " turn it off");                                                                                                
      return StatusCode::FAILURE;                                                                                                                                                 
    }                                                                                                                                                                             
    else {                                                                                                                                                                         
      ATH_MSG_INFO("retrieved tool: " << m_triggerDecision );                                                                                                                    
    }
  }
 
  ATH_MSG_DEBUG("** IDPerfMonZmumu::Initialize ** Completed **");
  return StatusCode::SUCCESS;
}

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

isClonable()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::isClonable ( ) const

inlineoverridevirtualinherited

Specify if the algorithm is clonable.

Only relevant for non-reentrant algorithms. Actual number of clones needs to be set via the "Cardinality" property.

Reimplemented in AFP_DigiTop, AlgB, AlgT, BCM_Digitization, CscDigitBuilder, CscDigitToCscRDO, G4AtlasAlg, G4RunAlg, HGTD_Digitization, HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ISF::SimKernelMT, ITk::StripDigitization, ITkPixelCablingAlg, ITkStripCablingAlg, LArHitEMapMaker, LArTTL1Maker, LUCID_DigiTop, LVL1::L1TopoSimulation, MergeCalibHits, MergeGenericMuonSimHitColl, MergeHijingPars, MergeMcEventCollection, MergeTrackRecordCollection, MergeTruthJets, MergeTruthParticles, MuonDigitizer, PileUpMTAlg, PixelDigitization, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_Digitization, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, SGInputLoader, Simulation::BeamEffectsAlg, TileHitVecToCnt, TileMuonFitter, TilePulseForTileMuonReceiver, TileRawChannelMaker, TRTDigitization, and ZDC_DigiTop.

Definition at line 68 of file AthCommonAlgorithm.h.

                                           {
    return true;
  }

isReEntrant()

virtual bool AthAlgorithm::isReEntrant ( ) const

inlinefinaloverrideprotectedvirtualinherited

Legacy algorithms are not thread-safe.

Definition at line 47 of file AthAlgorithm.h.

{ return false; }

msg()

MsgStream & AthCommonMsg< Gaudi::Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

RegisterHistograms()

void IDPerfMonZmumu::RegisterHistograms ( )

private

Definition at line 927 of file IDPerfMonZmumu.cxx.

{
  return;
}

renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

renounceArray()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

ResetCommonNtupleVectors()

void IDPerfMonZmumu::ResetCommonNtupleVectors ( )

private

Definition at line 2385 of file IDPerfMonZmumu.cxx.

{
  m_IDTrack_pt.clear();
  m_CBTrack_pt.clear();
  m_Refit1_pt.clear();
  m_Refit2_pt.clear();
  m_Truth_pt.clear();
 
  m_IDTrack_eta.clear();
  m_CBTrack_eta.clear();
  m_Refit1_eta.clear();
  m_Refit2_eta.clear();
  m_Truth_eta.clear();
 
  m_IDTrack_phi.clear();
  m_CBTrack_phi.clear();
  m_Refit1_phi.clear();
  m_Refit2_phi.clear();
  m_Truth_phi.clear();
 
  m_IDTrack_d0.clear();
  m_CBTrack_d0.clear();
  m_Refit1_d0.clear();
  m_Refit2_d0.clear();
  m_Truth_d0.clear();
 
  m_IDTrack_z0.clear();
  m_CBTrack_z0.clear();
  m_Refit1_z0.clear();
  m_Refit2_z0.clear();
  m_Truth_z0.clear();
 
  m_IDTrack_qoverp.clear();
  m_CBTrack_qoverp.clear();
  m_Refit1_qoverp.clear();
  m_Refit2_qoverp.clear();
  m_Truth_qoverp.clear();
 
  m_Truth_parent.clear();
 
  m_IDTrack_sigma_pt.clear();
  m_CBTrack_sigma_pt.clear();
  m_Refit1_sigma_pt.clear();
  m_Refit2_sigma_pt.clear();
 
  m_IDTrack_sigma_d0.clear();
  m_CBTrack_sigma_d0.clear();
  m_Refit1_sigma_d0.clear();
  m_Refit2_sigma_d0.clear();
 
  m_IDTrack_sigma_z0.clear();
  m_CBTrack_sigma_z0.clear();
  m_Refit1_sigma_z0.clear();
  m_Refit2_sigma_z0.clear();
 
  m_IDTrack_sigma_qoverp.clear();
  m_CBTrack_sigma_qoverp.clear();
  m_Refit1_sigma_qoverp.clear();
  m_Refit2_sigma_qoverp.clear();
 
  m_nBLhits.clear();
  m_nPIXhits.clear();
  m_nSCThits.clear();
  m_nTRThits.clear();
 
  return;
}

RunFourLeptonAnalysis()

StatusCode IDPerfMonZmumu::RunFourLeptonAnalysis ( )

private

Definition at line 2541 of file IDPerfMonZmumu.cxx.

{
  StatusCode thisStatus = StatusCode::SUCCESS;
 
  if (m_doFourMuAnalysis) {
    ATH_MSG_DEBUG ("** RunFourLeptonAnalysis ** START **    run: " << m_runNumber << "    event: " << m_evtNumber << "    lumiblock:" << m_lumi_block);  
 
    // Four lepton event
    m_4mu.setContainer(PerfMonServices::MUON_COLLECTION);
    m_4mu.doIsoSelection     (m_doIsoSelection);
    m_4mu.doIPSelection      (m_doIPSelection);
    m_4mu.doMCPSelection     (m_doMCPSelection);
    m_4mu.SetMassWindowLow   (m_MassWindowLow);
    m_4mu.SetMassWindowHigh  (m_MassWindowHigh);
    m_4mu.SetLeadingMuonPtCut(m_LeadingMuonPtCut);
    m_4mu.SetSecondMuonPtCut (m_SecondMuonPtCut);
    m_4mu.SetOpeningAngleCut (m_OpeningAngleCut);
    m_4mu.SetZ0GapCut        (m_Z0GapCut);
    m_4mu.setDebugMode       (m_doDebug);
 
    // This is a report of the number of true electrons in the event
    if (false) {
      const xAOD::TruthParticleContainer* electronTruth = evtStore()->retrieve< const xAOD::TruthParticleContainer >( "egammaTruthParticles" );
      if( electronTruth) {
    ATH_MSG_DEBUG ( "retrieving electron Truth container with key: \"egammaTruthParticles\" SUCCESS ");
    // Check that the auxiliary store association was made successfully:
    if( electronTruth->hasStore() ) {
      ATH_MSG_DEBUG ( "Size of Truth gamma-electrons: " << electronTruth->size() );
      auto tr_itr = electronTruth->begin();
      auto tr_end = electronTruth->end();
      int electronCount = 0;
      for( int i = 0; tr_itr != tr_end; ++tr_itr, ++i ) {
        ATH_MSG_DEBUG( "Truth : Investigating truth electron #" << i << "    pdgID= " << ( *tr_itr )->pdgId());
        if (fabs(( *tr_itr )->pdgId()) == 11) {
          if (( *tr_itr )->pt() > 10000 && fabs(( *tr_itr )->eta()) < 2.47) {
        ATH_MSG_INFO( "        central electron found --> pt: " << ( *tr_itr )->charge()*( *tr_itr )->pt() << " eta  " <<( *tr_itr )->eta()  );
        electronCount++;
          }
        }
      }
      ATH_MSG_INFO( "#central electrons: " << electronCount );
    }
    else {
      ATH_MSG_WARNING ( "electronTruth container has no Store !!! FAILURE");
    }
      }
      else {
    ATH_MSG_WARNING ( "retrieve electron Truth container with key: \"egammaTruthParticles\" FAILURE");
      }
    }
 
    if (m_doDebug) std::cout << " ** RunFourLeptonAnalysis ** calling m_4mu.Reco() .... " << std::endl;  
    bool statusOf4LeptReco = m_4mu.Reco();
 
    // print message 
    if(statusOf4LeptReco) {
      thisStatus = StatusCode::SUCCESS;
      ATH_MSG_INFO ("4-lepton reconstruction SUCCESS. # accepted events " << m_4mu.getAcceptedEvents() << "  Invariant mass = " << m_4mu.GetInvMass() << " GeV ");
      if (m_doDebug) std::cout << " ** RunFourLeptonAnalysis ** 4-lepton reconstruction SUCCESS. # accepted events " << m_4mu.getAcceptedEvents() << "  Invariant mass = " << m_4mu.GetInvMass() << " GeV " << std::endl;
    }
    else {
      thisStatus = StatusCode::FAILURE;
      ATH_MSG_DEBUG ("4-lepton reconstruction FAILURE.  m_4mu.Reco() returned FALSE --> event failed selection");
      if (m_doDebug) std::cout << " ** RunFourLeptonAnalysis ** 4-lepton reconstruction FAILURE.  m_4mu.Reco() returned FALSE --> event failed selection" << std::endl;
    }
 
    if(statusOf4LeptReco) {
      if ( m_4mu.EventPassed() ) {
 
    this->Clear4MuNtupleVariables(); 
 
    const xAOD::TrackParticle* muon1_pos = m_4mu.getIDTrack(m_4mu.getPosMuon(1));
    const xAOD::TrackParticle* muon2_pos = m_4mu.getIDTrack(m_4mu.getPosMuon(2));
    const xAOD::TrackParticle* muon1_neg = m_4mu.getIDTrack(m_4mu.getNegMuon(1));
    const xAOD::TrackParticle* muon2_neg = m_4mu.getIDTrack(m_4mu.getNegMuon(2));
 
    // muons
    if (muon1_pos) {
      m_positive_px = muon1_pos->p4().Px();
      m_positive_py = muon1_pos->p4().Py();
      m_positive_pz = muon1_pos->p4().Pz();
      m_positive_d0 = muon1_pos->d0();
      m_positive_z0 = muon1_pos->z0();
      // m_positive_d0_err = muon1_pos->definingParametersCovMatrix()(0,0);
      m_positive_d0_err = muon1_pos->definingParametersCovMatrixVec()[0];
      m_positive_z0_err = muon1_pos->definingParametersCovMatrix()(1,1);
    }
    if (muon1_neg) {
      m_negative_px = muon1_neg->p4().Px();
      m_negative_py = muon1_neg->p4().Py();
      m_negative_pz = muon1_neg->p4().Pz();
      m_negative_d0 = muon1_neg->d0();
      m_negative_z0 = muon1_neg->z0();
      m_negative_d0_err = muon1_neg->definingParametersCovMatrix()(0,0);
      m_negative_z0_err = muon1_neg->definingParametersCovMatrix()(1,1);
    }
    if (muon2_pos) {
      m_positive_2_px = muon2_pos->p4().Px();
      m_positive_2_py = muon2_pos->p4().Py();
      m_positive_2_pz = muon2_pos->p4().Pz();
      m_positive_2_d0 = muon2_pos->d0();
      m_positive_2_z0 = muon2_pos->z0();
      m_positive_2_d0_err = muon2_pos->definingParametersCovMatrix()(0,0);
      m_positive_2_z0_err = muon2_pos->definingParametersCovMatrix()(1,1);
    }
    if (muon2_neg) {
      m_negative_2_px = muon2_neg->p4().Px();
      m_negative_2_py = muon2_neg->p4().Py();
      m_negative_2_pz = muon2_neg->p4().Pz();
      m_negative_2_d0 = muon2_neg->d0();
      m_negative_2_z0 = muon2_neg->z0();
      m_negative_2_d0_err = muon2_neg->definingParametersCovMatrix()(0,0);
      m_negative_2_z0_err = muon2_neg->definingParametersCovMatrix()(1,1);
    }
    
    // electrons
    const xAOD::TrackParticle* elec1_neg = m_4mu.getELTrack(0);
    const xAOD::TrackParticle* elec2_neg = m_4mu.getELTrack(1);
    const xAOD::TrackParticle* elec1_pos = m_4mu.getELTrack(2);
    const xAOD::TrackParticle* elec2_pos = m_4mu.getELTrack(3);
    if (elec1_neg) {
      m_el_negative1_px = elec1_neg->p4().Px();
      m_el_negative1_py = elec1_neg->p4().Py();
      m_el_negative1_pz = elec1_neg->p4().Pz();
      m_el_negative1_d0 = elec1_neg->d0();
      m_el_negative1_z0 = elec1_neg->z0();
      m_el_negative1_d0_err = elec1_neg->definingParametersCovMatrix()(0,0);
      m_el_negative1_z0_err = elec1_neg->definingParametersCovMatrix()(1,1);
    }
    if (elec2_neg) {
      m_el_negative2_px = elec2_neg->p4().Px();
      m_el_negative2_py = elec2_neg->p4().Py();
      m_el_negative2_pz = elec2_neg->p4().Pz();
      m_el_negative2_d0 = elec2_neg->d0();
      m_el_negative2_z0 = elec2_neg->z0();
      m_el_negative2_d0_err = elec2_neg->definingParametersCovMatrix()(0,0);
      m_el_negative2_z0_err = elec2_neg->definingParametersCovMatrix()(1,1);
    }
    if (elec1_pos) {
      m_el_positive1_px = elec1_pos->p4().Px();
      m_el_positive1_py = elec1_pos->p4().Py();
      m_el_positive1_pz = elec1_pos->p4().Pz();
      m_el_positive1_d0 = elec1_pos->d0();
      m_el_positive1_z0 = elec1_pos->z0();
      m_el_positive1_d0_err = elec1_pos->definingParametersCovMatrix()(0,0);
      m_el_positive1_z0_err = elec1_pos->definingParametersCovMatrix()(1,1);
    }
    if (elec2_pos) {
      m_el_positive2_px = elec2_pos->p4().Px();
      m_el_positive2_py = elec2_pos->p4().Py();
      m_el_positive2_pz = elec2_pos->p4().Pz();
      m_el_positive2_d0 = elec2_pos->d0();
      m_el_positive2_z0 = elec2_pos->z0();
      m_el_positive2_d0_err = elec2_pos->definingParametersCovMatrix()(0,0);
      m_el_positive2_z0_err = elec2_pos->definingParametersCovMatrix()(1,1);
    }
    
    // other quantities
    m_nVertex =        m_4mu.GetNVertex ();
    m_negative_1_vtx = m_4mu.GetVertexMuNeg1();
    m_negative_2_vtx = m_4mu.GetVertexMuNeg2();
    m_positive_1_vtx = m_4mu.GetVertexMuPos1();
    m_positive_2_vtx = m_4mu.GetVertexMuPos2();
 
    m_el_negative1_vtx = m_4mu.GetVertexElec(0);
    m_el_negative2_vtx = m_4mu.GetVertexElec(1);
    m_el_positive1_vtx = m_4mu.GetVertexElec(2);
    m_el_positive2_vtx = m_4mu.GetVertexElec(3);
 
    m_pv_x = 0; m_pv_y = 0; m_pv_z = 0;
 
    m_4mu_minv = m_4mu.GetInvMass();
    
    // Obtain MET 
    std::string metName = "MET_Reference_AntiKt4LCTopo";
    std::string metRefFinalName = "FinalClus";
    const xAOD::MissingETContainer* final_met = nullptr;
 
    m_met = -1; // default value
    m_metphi = -1;
 
    if (!evtStore()->contains<xAOD::MissingETContainer>(metName)) {
      ATH_MSG_WARNING ( "No Collection with name " << metName << " found in StoreGate");
      // return StatusCode::SUCCESS;
    }
    else {
      StatusCode sc = evtStore()->retrieve(final_met, metName);
      if (sc.isFailure()) {
        ATH_MSG_DEBUG ( "Could not retrieve Collection " << metName << " from StoreGate");
        // return StatusCode::SUCCESS;
      }
    }
    const xAOD::MissingET *met = nullptr;
    if (final_met) met = (*final_met)[metRefFinalName];
    if (met) { // load MET values
      m_met = met->met();
      m_metphi = met->phi();
    }
    ATH_MSG_DEBUG (" Zmumu event with MET = " << m_met);
    
    ATH_MSG_INFO (" -- IDPerfMonZmumu::execute -- Accepted event " << m_4mu.getAcceptedEvents() << " with m_4mu.GetInvMass= " << m_4mu_minv);
    ATH_MSG_DEBUG ("**** Filling m_FourMuTree ntuple " << m_FourMuTree->GetName() << " entry " << m_FourMuTree->GetEntries()
               << " for run: " << m_runNumber 
               << "  event: " << m_evtNumber 
               << "  Lumiblock: " << m_lumi_block 
               << "  Invariant mass = " << m_4mu_minv << " GeV ");
    if (muon1_pos != nullptr) ATH_MSG_DEBUG("mu1+  --> pxyz " << muon1_pos->p4().Px() 
                        << ", " << muon1_pos->p4().Py()
                        << ", " << muon1_pos->p4().Pz()
                        << "  pt: " << muon1_pos->pt()
                        << "  d0: " << muon1_pos->d0() 
                        << "  z0: " << muon1_pos->z0()
                        << "  d0unc: " << muon1_pos->definingParametersCovMatrixVec()[0]
                        );
    if (muon1_neg != nullptr) ATH_MSG_DEBUG("mu1-  --> pxyz " << muon1_neg->p4().Px() 
                        << ", " << muon1_neg->p4().Py()
                        << ", " << muon1_neg->p4().Pz()
                        << "  pt: " << muon1_neg->pt()
                        << "  d0: " << muon1_neg->d0() 
                        << "  z0: " << muon1_neg->z0()
                        << "  d0unc: " << muon1_neg->definingParametersCovMatrixVec()[0]
                        );
    if (elec1_pos != nullptr) ATH_MSG_DEBUG("el1+  --> pxyz " << elec1_pos->p4().Px() 
                        << ", " << elec1_pos->p4().Py()
                        << ", " << elec1_pos->p4().Pz()
                        << "  pt: " << elec1_pos->pt()
                        << "  d0: " << elec1_pos->d0() 
                        << "  z0: " << elec1_pos->z0()
                        << "  d0unc: " << elec1_pos->definingParametersCovMatrixVec()[0]
                        );
      // ntuple variables have been filled in FillRecParametersTP
    m_FourMuTree->Fill();
      }
    } // succesful 4mu reco
 
  } // end of fourMuon Analysis
 
  return thisStatus;
}

setFilterPassed()

virtual void AthCommonAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state, const EventContext & ctx ) const

inlinevirtualinherited

Set filter decision:

Reimplemented in AthFilterAlgorithm.

Definition at line 99 of file AthCommonAlgorithm.h.

                                                                            {
    execState( ctx ).setFilterPassed( state );
  }

sysExecute()

StatusCode AthCommonAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx )

overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 80 of file AthCommonAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

sysInitialize()

StatusCode AthCommonAlgorithm< Gaudi::Algorithm >::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, AthHistogramAlgorithm, HypoBase, InputMakerBase, and PyAthena::Alg.

Definition at line 60 of file AthCommonAlgorithm.cxx.

                                                      {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<BaseAlg>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }      
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;  
}

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

updateVHKA()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_4mu

FourMuonEvent IDPerfMonZmumu::m_4mu

private

Definition at line 108 of file IDPerfMonZmumu.h.

m_4mu_minv

double IDPerfMonZmumu::m_4mu_minv {}

private

Definition at line 356 of file IDPerfMonZmumu.h.

{};

m_beamposX

float IDPerfMonZmumu::m_beamposX {}

private

Definition at line 233 of file IDPerfMonZmumu.h.

{};

m_beamposY

float IDPerfMonZmumu::m_beamposY {}

private

Definition at line 234 of file IDPerfMonZmumu.h.

{};

m_beamSpotKey

SG::ReadCondHandleKey<InDet::BeamSpotData> IDPerfMonZmumu::m_beamSpotKey { this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot" }

private

used for truth parameters

Definition at line 155 of file IDPerfMonZmumu.h.

{ this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot" };

m_CBTrack_d0

std::vector<float> IDPerfMonZmumu::m_CBTrack_d0

private

Definition at line 386 of file IDPerfMonZmumu.h.

m_CBTrack_eta

std::vector<float> IDPerfMonZmumu::m_CBTrack_eta

private

Definition at line 374 of file IDPerfMonZmumu.h.

m_CBTrack_phi

std::vector<float> IDPerfMonZmumu::m_CBTrack_phi

private

Definition at line 380 of file IDPerfMonZmumu.h.

m_CBTrack_pt

std::vector<float> IDPerfMonZmumu::m_CBTrack_pt

private

Definition at line 368 of file IDPerfMonZmumu.h.

m_CBTrack_qoverp

std::vector<float> IDPerfMonZmumu::m_CBTrack_qoverp

private

Definition at line 398 of file IDPerfMonZmumu.h.

m_CBTrack_sigma_d0

std::vector<float> IDPerfMonZmumu::m_CBTrack_sigma_d0

private

Definition at line 411 of file IDPerfMonZmumu.h.

m_CBTrack_sigma_pt

std::vector<float> IDPerfMonZmumu::m_CBTrack_sigma_pt

private

Definition at line 406 of file IDPerfMonZmumu.h.

m_CBTrack_sigma_qoverp

std::vector<float> IDPerfMonZmumu::m_CBTrack_sigma_qoverp

private

Definition at line 421 of file IDPerfMonZmumu.h.

m_CBTrack_sigma_z0

std::vector<float> IDPerfMonZmumu::m_CBTrack_sigma_z0

private

Definition at line 416 of file IDPerfMonZmumu.h.

m_CBTrack_z0

std::vector<float> IDPerfMonZmumu::m_CBTrack_z0

private

Definition at line 392 of file IDPerfMonZmumu.h.

m_combTree

TTree* IDPerfMonZmumu::m_combTree

private

Definition at line 199 of file IDPerfMonZmumu.h.

m_combTreeFolder

std::string IDPerfMonZmumu::m_combTreeFolder

private

Definition at line 188 of file IDPerfMonZmumu.h.

m_combTreeName

std::string IDPerfMonZmumu::m_combTreeName

private

Definition at line 172 of file IDPerfMonZmumu.h.

m_commonTree

TTree* IDPerfMonZmumu::m_commonTree

private

Definition at line 193 of file IDPerfMonZmumu.h.

m_commonTreeFolder

std::string IDPerfMonZmumu::m_commonTreeFolder

private

Definition at line 182 of file IDPerfMonZmumu.h.

m_commonTreeName

std::string IDPerfMonZmumu::m_commonTreeName

private

Definition at line 166 of file IDPerfMonZmumu.h.

m_defaultTree

TTree* IDPerfMonZmumu::m_defaultTree

private

Definition at line 194 of file IDPerfMonZmumu.h.

m_defaultTreeFolder

std::string IDPerfMonZmumu::m_defaultTreeFolder

private

Definition at line 183 of file IDPerfMonZmumu.h.

m_defaultTreeName

std::string IDPerfMonZmumu::m_defaultTreeName

private

Definition at line 167 of file IDPerfMonZmumu.h.

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

m_doDebug

bool IDPerfMonZmumu::m_doDebug {}

private

Definition at line 207 of file IDPerfMonZmumu.h.

{};

m_doFourMuAnalysis

bool IDPerfMonZmumu::m_doFourMuAnalysis {}

private

Definition at line 123 of file IDPerfMonZmumu.h.

{};

m_doIP

bool IDPerfMonZmumu::m_doIP {}

private

Definition at line 122 of file IDPerfMonZmumu.h.

{};

m_doIPSelection

bool IDPerfMonZmumu::m_doIPSelection {true}

private

Definition at line 111 of file IDPerfMonZmumu.h.

{true};

m_doIsoSelection

bool IDPerfMonZmumu::m_doIsoSelection {true}

private

Definition at line 110 of file IDPerfMonZmumu.h.

{true};

m_doMCPSelection

bool IDPerfMonZmumu::m_doMCPSelection {true}

private

Definition at line 112 of file IDPerfMonZmumu.h.

{true};

m_doRefit

bool IDPerfMonZmumu::m_doRefit {}

private

Definition at line 120 of file IDPerfMonZmumu.h.

{};

m_doRemoval

bool IDPerfMonZmumu::m_doRemoval {}

private

Definition at line 206 of file IDPerfMonZmumu.h.

{};

m_el_negative1_d0

double IDPerfMonZmumu::m_el_negative1_d0 {}

private

Definition at line 307 of file IDPerfMonZmumu.h.

{};

m_el_negative1_d0_err

double IDPerfMonZmumu::m_el_negative1_d0_err {}

private

Definition at line 309 of file IDPerfMonZmumu.h.

{};

m_el_negative1_d0_PV

double IDPerfMonZmumu::m_el_negative1_d0_PV {}

private

Definition at line 311 of file IDPerfMonZmumu.h.

{};

m_el_negative1_d0_PVerr

double IDPerfMonZmumu::m_el_negative1_d0_PVerr {}

private

Definition at line 313 of file IDPerfMonZmumu.h.

{};

m_el_negative1_px

double IDPerfMonZmumu::m_el_negative1_px {}

private

Definition at line 303 of file IDPerfMonZmumu.h.

{};

m_el_negative1_py

double IDPerfMonZmumu::m_el_negative1_py {}

private

Definition at line 304 of file IDPerfMonZmumu.h.

{};

m_el_negative1_pz

double IDPerfMonZmumu::m_el_negative1_pz {}

private

Definition at line 305 of file IDPerfMonZmumu.h.

{};

m_el_negative1_vtx

int IDPerfMonZmumu::m_el_negative1_vtx {}

private

Definition at line 314 of file IDPerfMonZmumu.h.

{};

m_el_negative1_z0

double IDPerfMonZmumu::m_el_negative1_z0 {}

private

Definition at line 306 of file IDPerfMonZmumu.h.

{};

m_el_negative1_z0_err

double IDPerfMonZmumu::m_el_negative1_z0_err {}

private

Definition at line 308 of file IDPerfMonZmumu.h.

{};

m_el_negative1_z0_PV

double IDPerfMonZmumu::m_el_negative1_z0_PV {}

private

Definition at line 310 of file IDPerfMonZmumu.h.

{};

m_el_negative1_z0_PVerr

double IDPerfMonZmumu::m_el_negative1_z0_PVerr {}

private

Definition at line 312 of file IDPerfMonZmumu.h.

{};

m_el_negative2_d0

double IDPerfMonZmumu::m_el_negative2_d0 {}

private

Definition at line 320 of file IDPerfMonZmumu.h.

{};

m_el_negative2_d0_err

double IDPerfMonZmumu::m_el_negative2_d0_err {}

private

Definition at line 322 of file IDPerfMonZmumu.h.

{};

m_el_negative2_d0_PV

double IDPerfMonZmumu::m_el_negative2_d0_PV {}

private

Definition at line 324 of file IDPerfMonZmumu.h.

{};

m_el_negative2_d0_PVerr

double IDPerfMonZmumu::m_el_negative2_d0_PVerr {}

private

Definition at line 326 of file IDPerfMonZmumu.h.

{};

m_el_negative2_px

double IDPerfMonZmumu::m_el_negative2_px {}

private

Definition at line 316 of file IDPerfMonZmumu.h.

{};

m_el_negative2_py

double IDPerfMonZmumu::m_el_negative2_py {}

private

Definition at line 317 of file IDPerfMonZmumu.h.

{};

m_el_negative2_pz

double IDPerfMonZmumu::m_el_negative2_pz {}

private

Definition at line 318 of file IDPerfMonZmumu.h.

{};

m_el_negative2_vtx

int IDPerfMonZmumu::m_el_negative2_vtx {}

private

Definition at line 327 of file IDPerfMonZmumu.h.

{};

m_el_negative2_z0

double IDPerfMonZmumu::m_el_negative2_z0 {}

private

Definition at line 319 of file IDPerfMonZmumu.h.

{};

m_el_negative2_z0_err

double IDPerfMonZmumu::m_el_negative2_z0_err {}

private

Definition at line 321 of file IDPerfMonZmumu.h.

{};

m_el_negative2_z0_PV

double IDPerfMonZmumu::m_el_negative2_z0_PV {}

private

Definition at line 323 of file IDPerfMonZmumu.h.

{};

m_el_negative2_z0_PVerr

double IDPerfMonZmumu::m_el_negative2_z0_PVerr {}

private

Definition at line 325 of file IDPerfMonZmumu.h.

{};

m_el_positive1_d0

double IDPerfMonZmumu::m_el_positive1_d0 {}

private

Definition at line 333 of file IDPerfMonZmumu.h.

{};

m_el_positive1_d0_err

double IDPerfMonZmumu::m_el_positive1_d0_err {}

private

Definition at line 335 of file IDPerfMonZmumu.h.

{};

m_el_positive1_d0_PV

double IDPerfMonZmumu::m_el_positive1_d0_PV {}

private

Definition at line 337 of file IDPerfMonZmumu.h.

{};

m_el_positive1_px

double IDPerfMonZmumu::m_el_positive1_px {}

private

Definition at line 329 of file IDPerfMonZmumu.h.

{};

m_el_positive1_py

double IDPerfMonZmumu::m_el_positive1_py {}

private

Definition at line 330 of file IDPerfMonZmumu.h.

{};

m_el_positive1_pz

double IDPerfMonZmumu::m_el_positive1_pz {}

private

Definition at line 331 of file IDPerfMonZmumu.h.

{};

m_el_positive1_vtx

int IDPerfMonZmumu::m_el_positive1_vtx {}

private

Definition at line 340 of file IDPerfMonZmumu.h.

{};

m_el_positive1_z0

double IDPerfMonZmumu::m_el_positive1_z0 {}

private

Definition at line 332 of file IDPerfMonZmumu.h.

{};

m_el_positive1_z0_err

double IDPerfMonZmumu::m_el_positive1_z0_err {}

private

Definition at line 334 of file IDPerfMonZmumu.h.

{};

m_el_positive1_z0_PV

double IDPerfMonZmumu::m_el_positive1_z0_PV {}

private

Definition at line 336 of file IDPerfMonZmumu.h.

{};

m_el_positive1_z0_PVerr

double IDPerfMonZmumu::m_el_positive1_z0_PVerr {}

private

Definition at line 338 of file IDPerfMonZmumu.h.

{};

m_el_positive2_d0

double IDPerfMonZmumu::m_el_positive2_d0 {}

private

Definition at line 346 of file IDPerfMonZmumu.h.

{};

m_el_positive2_d0_err

double IDPerfMonZmumu::m_el_positive2_d0_err {}

private

Definition at line 348 of file IDPerfMonZmumu.h.

{};

m_el_positive2_d0_PV

double IDPerfMonZmumu::m_el_positive2_d0_PV {}

private

Definition at line 350 of file IDPerfMonZmumu.h.

{};

m_el_positive2_px

double IDPerfMonZmumu::m_el_positive2_px {}

private

Definition at line 342 of file IDPerfMonZmumu.h.

{};

m_el_positive2_py

double IDPerfMonZmumu::m_el_positive2_py {}

private

Definition at line 343 of file IDPerfMonZmumu.h.

{};

m_el_positive2_pz

double IDPerfMonZmumu::m_el_positive2_pz {}

private

Definition at line 344 of file IDPerfMonZmumu.h.

{};

m_el_positive2_vtx

int IDPerfMonZmumu::m_el_positive2_vtx {}

private

Definition at line 353 of file IDPerfMonZmumu.h.

{};

m_el_positive2_z0

double IDPerfMonZmumu::m_el_positive2_z0 {}

private

Definition at line 345 of file IDPerfMonZmumu.h.

{};

m_el_positive2_z0_err

double IDPerfMonZmumu::m_el_positive2_z0_err {}

private

Definition at line 347 of file IDPerfMonZmumu.h.

{};

m_el_positive2_z0_PV

double IDPerfMonZmumu::m_el_positive2_z0_PV {}

private

Definition at line 349 of file IDPerfMonZmumu.h.

{};

m_el_positive2_z0_PVerr

double IDPerfMonZmumu::m_el_positive2_z0_PVerr {}

private

Definition at line 351 of file IDPerfMonZmumu.h.

{};

m_el_psoitive1_d0_PVerr

double IDPerfMonZmumu::m_el_psoitive1_d0_PVerr {}

private

Definition at line 339 of file IDPerfMonZmumu.h.

{};

m_el_psoitive2_d0_PVerr

double IDPerfMonZmumu::m_el_psoitive2_d0_PVerr {}

private

Definition at line 352 of file IDPerfMonZmumu.h.

{};

m_event_mu

unsigned int IDPerfMonZmumu::m_event_mu {}

private

Definition at line 228 of file IDPerfMonZmumu.h.

{};

m_event_weight

float IDPerfMonZmumu::m_event_weight {}

private

Definition at line 229 of file IDPerfMonZmumu.h.

{};

m_EventInfoKey

SG::ReadHandleKey<xAOD::EventInfo> IDPerfMonZmumu::m_EventInfoKey {this, "EventInfoKey", "EventInfo"}

private

Definition at line 436 of file IDPerfMonZmumu.h.

{this, "EventInfoKey", "EventInfo"};

m_evtNumber

unsigned int IDPerfMonZmumu::m_evtNumber {}

private

Definition at line 226 of file IDPerfMonZmumu.h.

{};

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthCommonAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 108 of file AthCommonAlgorithm.h.

m_extrapolator

ToolHandle<Trk::IExtrapolator> IDPerfMonZmumu::m_extrapolator

private

Definition at line 156 of file IDPerfMonZmumu.h.

m_FourMuTree

TTree* IDPerfMonZmumu::m_FourMuTree

private

Definition at line 201 of file IDPerfMonZmumu.h.

m_FourMuTreeFolder

std::string IDPerfMonZmumu::m_FourMuTreeFolder

private

Root Validation Tree.

Definition at line 190 of file IDPerfMonZmumu.h.

m_FourMuTreeName

std::string IDPerfMonZmumu::m_FourMuTreeName

private

Definition at line 174 of file IDPerfMonZmumu.h.

m_h_cutflow

TH1F* IDPerfMonZmumu::m_h_cutflow {}

private

link vector to map HepMC onto xAOD truth

Definition at line 223 of file IDPerfMonZmumu.h.

{};

m_IDTrack_d0

std::vector<float> IDPerfMonZmumu::m_IDTrack_d0

private

Definition at line 385 of file IDPerfMonZmumu.h.

m_IDTrack_eta

std::vector<float> IDPerfMonZmumu::m_IDTrack_eta

private

Definition at line 373 of file IDPerfMonZmumu.h.

m_IDTrack_phi

std::vector<float> IDPerfMonZmumu::m_IDTrack_phi

private

Definition at line 379 of file IDPerfMonZmumu.h.

m_IDTrack_pt

std::vector<float> IDPerfMonZmumu::m_IDTrack_pt

private

Definition at line 367 of file IDPerfMonZmumu.h.

m_IDTrack_qoverp

std::vector<float> IDPerfMonZmumu::m_IDTrack_qoverp

private

Definition at line 397 of file IDPerfMonZmumu.h.

m_IDTrack_sigma_d0

std::vector<float> IDPerfMonZmumu::m_IDTrack_sigma_d0

private

Definition at line 410 of file IDPerfMonZmumu.h.

m_IDTrack_sigma_pt

std::vector<float> IDPerfMonZmumu::m_IDTrack_sigma_pt

private

Definition at line 405 of file IDPerfMonZmumu.h.

m_IDTrack_sigma_qoverp

std::vector<float> IDPerfMonZmumu::m_IDTrack_sigma_qoverp

private

Definition at line 420 of file IDPerfMonZmumu.h.

m_IDTrack_sigma_z0

std::vector<float> IDPerfMonZmumu::m_IDTrack_sigma_z0

private

Definition at line 415 of file IDPerfMonZmumu.h.

m_IDTrack_z0

std::vector<float> IDPerfMonZmumu::m_IDTrack_z0

private

Definition at line 391 of file IDPerfMonZmumu.h.

m_IDTree

TTree* IDPerfMonZmumu::m_IDTree

private

Definition at line 195 of file IDPerfMonZmumu.h.

m_IDTreeFolder

std::string IDPerfMonZmumu::m_IDTreeFolder

private

Definition at line 184 of file IDPerfMonZmumu.h.

m_IDTreeName

std::string IDPerfMonZmumu::m_IDTreeName

private

Definition at line 168 of file IDPerfMonZmumu.h.

m_isMC

bool IDPerfMonZmumu::m_isMC {}

private

Definition at line 119 of file IDPerfMonZmumu.h.

{};

m_LeadingMuonPtCut

double IDPerfMonZmumu::m_LeadingMuonPtCut {20.}

private

Definition at line 115 of file IDPerfMonZmumu.h.

{20.};

m_lumi_block

unsigned int IDPerfMonZmumu::m_lumi_block {}

private

Definition at line 227 of file IDPerfMonZmumu.h.

{};

m_MassWindowHigh

double IDPerfMonZmumu::m_MassWindowHigh {120.}

private

Definition at line 114 of file IDPerfMonZmumu.h.

{120.};

m_MassWindowLow

double IDPerfMonZmumu::m_MassWindowLow {60.}

private

Definition at line 113 of file IDPerfMonZmumu.h.

{60.};

m_maxGoodLumiBlock

int IDPerfMonZmumu::m_maxGoodLumiBlock {}

private

Definition at line 128 of file IDPerfMonZmumu.h.

{};

m_met

double IDPerfMonZmumu::m_met {}

private

Definition at line 363 of file IDPerfMonZmumu.h.

{};

m_metphi

double IDPerfMonZmumu::m_metphi {}

private

Definition at line 364 of file IDPerfMonZmumu.h.

{};

m_minGoodLumiBlock

int IDPerfMonZmumu::m_minGoodLumiBlock {}

private

Definition at line 127 of file IDPerfMonZmumu.h.

{};

m_MSTree

TTree* IDPerfMonZmumu::m_MSTree

private

Definition at line 200 of file IDPerfMonZmumu.h.

m_MSTreeFolder

std::string IDPerfMonZmumu::m_MSTreeFolder

private

Definition at line 189 of file IDPerfMonZmumu.h.

m_MSTreeName

std::string IDPerfMonZmumu::m_MSTreeName

private

Definition at line 173 of file IDPerfMonZmumu.h.

m_MuonQualityName

std::string IDPerfMonZmumu::m_MuonQualityName

private

Definition at line 433 of file IDPerfMonZmumu.h.

m_muonSelector

ToolHandle<CP::IMuonSelectionTool> IDPerfMonZmumu::m_muonSelector

private

used to pass a custom muon selector

Definition at line 159 of file IDPerfMonZmumu.h.

m_nBLhits

std::vector<int> IDPerfMonZmumu::m_nBLhits

private

Definition at line 425 of file IDPerfMonZmumu.h.

m_negative_1_vtx

int IDPerfMonZmumu::m_negative_1_vtx {}

private

Definition at line 286 of file IDPerfMonZmumu.h.

{};

m_negative_2_d0

double IDPerfMonZmumu::m_negative_2_d0 {}

private

Definition at line 293 of file IDPerfMonZmumu.h.

{};

m_negative_2_d0_err

double IDPerfMonZmumu::m_negative_2_d0_err {}

private

Definition at line 295 of file IDPerfMonZmumu.h.

{};

m_negative_2_d0_PV

double IDPerfMonZmumu::m_negative_2_d0_PV {}

private

Definition at line 297 of file IDPerfMonZmumu.h.

{};

m_negative_2_d0_PVerr

double IDPerfMonZmumu::m_negative_2_d0_PVerr {}

private

Definition at line 299 of file IDPerfMonZmumu.h.

{};

m_negative_2_px

double IDPerfMonZmumu::m_negative_2_px {}

private

Definition at line 289 of file IDPerfMonZmumu.h.

{};

m_negative_2_py

double IDPerfMonZmumu::m_negative_2_py {}

private

Definition at line 290 of file IDPerfMonZmumu.h.

{};

m_negative_2_pz

double IDPerfMonZmumu::m_negative_2_pz {}

private

Definition at line 291 of file IDPerfMonZmumu.h.

{};

m_negative_2_vtx

int IDPerfMonZmumu::m_negative_2_vtx {}

private

Definition at line 300 of file IDPerfMonZmumu.h.

{};

m_negative_2_z0

double IDPerfMonZmumu::m_negative_2_z0 {}

private

Definition at line 292 of file IDPerfMonZmumu.h.

{};

m_negative_2_z0_err

double IDPerfMonZmumu::m_negative_2_z0_err {}

private

Definition at line 294 of file IDPerfMonZmumu.h.

{};

m_negative_2_z0_PV

double IDPerfMonZmumu::m_negative_2_z0_PV {}

private

Definition at line 296 of file IDPerfMonZmumu.h.

{};

m_negative_2_z0_PVerr

double IDPerfMonZmumu::m_negative_2_z0_PVerr {}

private

Definition at line 298 of file IDPerfMonZmumu.h.

{};

m_negative_d0

double IDPerfMonZmumu::m_negative_d0 {}

private

Definition at line 276 of file IDPerfMonZmumu.h.

{};  

m_negative_d0_err

double IDPerfMonZmumu::m_negative_d0_err {}

private

Definition at line 278 of file IDPerfMonZmumu.h.

{};

m_negative_d0_PV

double IDPerfMonZmumu::m_negative_d0_PV {}

private

Definition at line 281 of file IDPerfMonZmumu.h.

{};

m_negative_d0_PVerr

double IDPerfMonZmumu::m_negative_d0_PVerr {}

private

Definition at line 283 of file IDPerfMonZmumu.h.

{};

m_negative_eta

double IDPerfMonZmumu::m_negative_eta {}

private

Definition at line 274 of file IDPerfMonZmumu.h.

{};

m_negative_parent

int IDPerfMonZmumu::m_negative_parent {}

private

Definition at line 287 of file IDPerfMonZmumu.h.

{};

m_negative_phi

double IDPerfMonZmumu::m_negative_phi {}

private

Definition at line 273 of file IDPerfMonZmumu.h.

{};

m_negative_pt

double IDPerfMonZmumu::m_negative_pt {}

private

Definition at line 271 of file IDPerfMonZmumu.h.

{};

m_negative_px

double IDPerfMonZmumu::m_negative_px {}

private

Definition at line 269 of file IDPerfMonZmumu.h.

{};

m_negative_py

double IDPerfMonZmumu::m_negative_py {}

private

Definition at line 270 of file IDPerfMonZmumu.h.

{};

m_negative_pz

double IDPerfMonZmumu::m_negative_pz {}

private

Definition at line 272 of file IDPerfMonZmumu.h.

{};

m_negative_qoverp

double IDPerfMonZmumu::m_negative_qoverp {}

private

Definition at line 284 of file IDPerfMonZmumu.h.

{};

m_negative_sigma_pt

double IDPerfMonZmumu::m_negative_sigma_pt {}

private

Definition at line 279 of file IDPerfMonZmumu.h.

{};

m_negative_sigma_qoverp

double IDPerfMonZmumu::m_negative_sigma_qoverp {}

private

Definition at line 285 of file IDPerfMonZmumu.h.

{};

m_negative_z0

double IDPerfMonZmumu::m_negative_z0 {}

private

Definition at line 275 of file IDPerfMonZmumu.h.

{};

m_negative_z0_err

double IDPerfMonZmumu::m_negative_z0_err {}

private

Definition at line 277 of file IDPerfMonZmumu.h.

{};

m_negative_z0_PV

double IDPerfMonZmumu::m_negative_z0_PV {}

private

Definition at line 280 of file IDPerfMonZmumu.h.

{};

m_negative_z0_PVerr

double IDPerfMonZmumu::m_negative_z0_PVerr {}

private

Definition at line 282 of file IDPerfMonZmumu.h.

{};

m_nPIXhits

std::vector<int> IDPerfMonZmumu::m_nPIXhits

private

Definition at line 426 of file IDPerfMonZmumu.h.

m_nSCThits

std::vector<int> IDPerfMonZmumu::m_nSCThits

private

Definition at line 427 of file IDPerfMonZmumu.h.

m_nTrkInVtx

unsigned int IDPerfMonZmumu::m_nTrkInVtx {}

private

Definition at line 361 of file IDPerfMonZmumu.h.

{};

m_nTRThits

std::vector<int> IDPerfMonZmumu::m_nTRThits

private

Definition at line 428 of file IDPerfMonZmumu.h.

m_nVertex

unsigned int IDPerfMonZmumu::m_nVertex {}

private

Definition at line 232 of file IDPerfMonZmumu.h.

{};

m_OpeningAngleCut

double IDPerfMonZmumu::m_OpeningAngleCut {0.2}

private

Definition at line 117 of file IDPerfMonZmumu.h.

{0.2};

m_outputTracksName

std::string IDPerfMonZmumu::m_outputTracksName

private

Definition at line 432 of file IDPerfMonZmumu.h.

m_positive_1_vtx

int IDPerfMonZmumu::m_positive_1_vtx {}

private

Definition at line 253 of file IDPerfMonZmumu.h.

{};

m_positive_2_d0

double IDPerfMonZmumu::m_positive_2_d0 {}

private

Definition at line 260 of file IDPerfMonZmumu.h.

{};

m_positive_2_d0_err

double IDPerfMonZmumu::m_positive_2_d0_err {}

private

Definition at line 262 of file IDPerfMonZmumu.h.

{};

m_positive_2_d0_PV

double IDPerfMonZmumu::m_positive_2_d0_PV {}

private

Definition at line 264 of file IDPerfMonZmumu.h.

{};

m_positive_2_d0_PVerr

double IDPerfMonZmumu::m_positive_2_d0_PVerr {}

private

Definition at line 266 of file IDPerfMonZmumu.h.

{};

m_positive_2_px

double IDPerfMonZmumu::m_positive_2_px {}

private

Definition at line 256 of file IDPerfMonZmumu.h.

{};

m_positive_2_py

double IDPerfMonZmumu::m_positive_2_py {}

private

Definition at line 257 of file IDPerfMonZmumu.h.

{};

m_positive_2_pz

double IDPerfMonZmumu::m_positive_2_pz {}

private

Definition at line 258 of file IDPerfMonZmumu.h.

{};

m_positive_2_vtx

int IDPerfMonZmumu::m_positive_2_vtx {}

private

Definition at line 267 of file IDPerfMonZmumu.h.

{};

m_positive_2_z0

double IDPerfMonZmumu::m_positive_2_z0 {}

private

Definition at line 259 of file IDPerfMonZmumu.h.

{};

m_positive_2_z0_err

double IDPerfMonZmumu::m_positive_2_z0_err {}

private

Definition at line 261 of file IDPerfMonZmumu.h.

{};

m_positive_2_z0_PV

double IDPerfMonZmumu::m_positive_2_z0_PV {}

private

Definition at line 263 of file IDPerfMonZmumu.h.

{};

m_positive_2_z0_PVerr

double IDPerfMonZmumu::m_positive_2_z0_PVerr {}

private

Definition at line 265 of file IDPerfMonZmumu.h.

{};

m_positive_d0

double IDPerfMonZmumu::m_positive_d0 {}

private

Definition at line 243 of file IDPerfMonZmumu.h.

{};

m_positive_d0_err

double IDPerfMonZmumu::m_positive_d0_err {}

private

Definition at line 245 of file IDPerfMonZmumu.h.

{};

m_positive_d0_PV

double IDPerfMonZmumu::m_positive_d0_PV {}

private

Definition at line 248 of file IDPerfMonZmumu.h.

{};

m_positive_d0_PVerr

double IDPerfMonZmumu::m_positive_d0_PVerr {}

private

Definition at line 250 of file IDPerfMonZmumu.h.

{};

m_positive_eta

double IDPerfMonZmumu::m_positive_eta {}

private

Definition at line 241 of file IDPerfMonZmumu.h.

{};

m_positive_parent

int IDPerfMonZmumu::m_positive_parent {}

private

Definition at line 254 of file IDPerfMonZmumu.h.

{};

m_positive_phi

double IDPerfMonZmumu::m_positive_phi {}

private

Definition at line 240 of file IDPerfMonZmumu.h.

{};

m_positive_pt

double IDPerfMonZmumu::m_positive_pt {}

private

Definition at line 238 of file IDPerfMonZmumu.h.

{};

m_positive_px

double IDPerfMonZmumu::m_positive_px {}

private

Definition at line 236 of file IDPerfMonZmumu.h.

{};

m_positive_py

double IDPerfMonZmumu::m_positive_py {}

private

Definition at line 237 of file IDPerfMonZmumu.h.

{};

m_positive_pz

double IDPerfMonZmumu::m_positive_pz {}

private

Definition at line 239 of file IDPerfMonZmumu.h.

{};

m_positive_qoverp

double IDPerfMonZmumu::m_positive_qoverp {}

private

Definition at line 251 of file IDPerfMonZmumu.h.

{};

m_positive_sigma_pt

double IDPerfMonZmumu::m_positive_sigma_pt {}

private

Definition at line 246 of file IDPerfMonZmumu.h.

{};

m_positive_sigma_qoverp

double IDPerfMonZmumu::m_positive_sigma_qoverp {}

private

Definition at line 252 of file IDPerfMonZmumu.h.

{};

m_positive_z0

double IDPerfMonZmumu::m_positive_z0 {}

private

Definition at line 242 of file IDPerfMonZmumu.h.

{};

m_positive_z0_err

double IDPerfMonZmumu::m_positive_z0_err {}

private

Definition at line 244 of file IDPerfMonZmumu.h.

{};

m_positive_z0_PV

double IDPerfMonZmumu::m_positive_z0_PV {}

private

Definition at line 247 of file IDPerfMonZmumu.h.

{};

m_positive_z0_PVerr

double IDPerfMonZmumu::m_positive_z0_PVerr {}

private

Definition at line 249 of file IDPerfMonZmumu.h.

{};

m_pv_x

double IDPerfMonZmumu::m_pv_x {}

private

Definition at line 358 of file IDPerfMonZmumu.h.

{};

m_pv_y

double IDPerfMonZmumu::m_pv_y {}

private

Definition at line 359 of file IDPerfMonZmumu.h.

{};

m_pv_z

double IDPerfMonZmumu::m_pv_z {}

private

Definition at line 360 of file IDPerfMonZmumu.h.

{};

m_Refit1_d0

std::vector<float> IDPerfMonZmumu::m_Refit1_d0

private

Definition at line 387 of file IDPerfMonZmumu.h.

m_Refit1_eta

std::vector<float> IDPerfMonZmumu::m_Refit1_eta

private

Definition at line 375 of file IDPerfMonZmumu.h.

m_Refit1_phi

std::vector<float> IDPerfMonZmumu::m_Refit1_phi

private

Definition at line 381 of file IDPerfMonZmumu.h.

m_Refit1_pt

std::vector<float> IDPerfMonZmumu::m_Refit1_pt

private

Definition at line 369 of file IDPerfMonZmumu.h.

m_Refit1_qoverp

std::vector<float> IDPerfMonZmumu::m_Refit1_qoverp

private

Definition at line 399 of file IDPerfMonZmumu.h.

m_Refit1_sigma_d0

std::vector<float> IDPerfMonZmumu::m_Refit1_sigma_d0

private

Definition at line 412 of file IDPerfMonZmumu.h.

m_Refit1_sigma_pt

std::vector<float> IDPerfMonZmumu::m_Refit1_sigma_pt

private

Definition at line 407 of file IDPerfMonZmumu.h.

m_Refit1_sigma_qoverp

std::vector<float> IDPerfMonZmumu::m_Refit1_sigma_qoverp

private

Definition at line 422 of file IDPerfMonZmumu.h.

m_Refit1_sigma_z0

std::vector<float> IDPerfMonZmumu::m_Refit1_sigma_z0

private

Definition at line 417 of file IDPerfMonZmumu.h.

m_Refit1_z0

std::vector<float> IDPerfMonZmumu::m_Refit1_z0

private

Definition at line 393 of file IDPerfMonZmumu.h.

m_refit1Tree

TTree* IDPerfMonZmumu::m_refit1Tree

private

Definition at line 196 of file IDPerfMonZmumu.h.

m_refit1TreeFolder

std::string IDPerfMonZmumu::m_refit1TreeFolder

private

Definition at line 185 of file IDPerfMonZmumu.h.

m_refit1TreeName

std::string IDPerfMonZmumu::m_refit1TreeName

private

Definition at line 169 of file IDPerfMonZmumu.h.

m_Refit2_d0

std::vector<float> IDPerfMonZmumu::m_Refit2_d0

private

Definition at line 388 of file IDPerfMonZmumu.h.

m_Refit2_eta

std::vector<float> IDPerfMonZmumu::m_Refit2_eta

private

Definition at line 376 of file IDPerfMonZmumu.h.

m_Refit2_phi

std::vector<float> IDPerfMonZmumu::m_Refit2_phi

private

Definition at line 382 of file IDPerfMonZmumu.h.

m_Refit2_pt

std::vector<float> IDPerfMonZmumu::m_Refit2_pt

private

Definition at line 370 of file IDPerfMonZmumu.h.

m_Refit2_qoverp

std::vector<float> IDPerfMonZmumu::m_Refit2_qoverp

private

Definition at line 400 of file IDPerfMonZmumu.h.

m_Refit2_sigma_d0

std::vector<float> IDPerfMonZmumu::m_Refit2_sigma_d0

private

Definition at line 413 of file IDPerfMonZmumu.h.

m_Refit2_sigma_pt

std::vector<float> IDPerfMonZmumu::m_Refit2_sigma_pt

private

Definition at line 408 of file IDPerfMonZmumu.h.

m_Refit2_sigma_qoverp

std::vector<float> IDPerfMonZmumu::m_Refit2_sigma_qoverp

private

Definition at line 423 of file IDPerfMonZmumu.h.

m_Refit2_sigma_z0

std::vector<float> IDPerfMonZmumu::m_Refit2_sigma_z0

private

Definition at line 418 of file IDPerfMonZmumu.h.

m_Refit2_z0

std::vector<float> IDPerfMonZmumu::m_Refit2_z0

private

Definition at line 394 of file IDPerfMonZmumu.h.

m_refit2Tree

TTree* IDPerfMonZmumu::m_refit2Tree

private

Definition at line 197 of file IDPerfMonZmumu.h.

m_refit2TreeFolder

std::string IDPerfMonZmumu::m_refit2TreeFolder

private

Definition at line 186 of file IDPerfMonZmumu.h.

m_refit2TreeName

std::string IDPerfMonZmumu::m_refit2TreeName

private

Definition at line 170 of file IDPerfMonZmumu.h.

m_runNumber

unsigned int IDPerfMonZmumu::m_runNumber {}

private

Definition at line 225 of file IDPerfMonZmumu.h.

{};

m_SecondMuonPtCut

double IDPerfMonZmumu::m_SecondMuonPtCut {15.}

private

Definition at line 116 of file IDPerfMonZmumu.h.

{15.};

m_selTool

ToolHandle< InDet::IInDetTrackSelectionTool > IDPerfMonZmumu::m_selTool

private

The track selection Tool.

Definition at line 148 of file IDPerfMonZmumu.h.

m_skipMS

bool IDPerfMonZmumu::m_skipMS {}

private

Definition at line 125 of file IDPerfMonZmumu.h.

{};

m_storeZmumuNtuple

bool IDPerfMonZmumu::m_storeZmumuNtuple {true}

private

Definition at line 124 of file IDPerfMonZmumu.h.

{true};

m_sTriggerChainName

std::string IDPerfMonZmumu::m_sTriggerChainName

private

Definition at line 431 of file IDPerfMonZmumu.h.

m_trackContainerName

SG::ReadHandleKey<xAOD::TrackParticleContainer> IDPerfMonZmumu::m_trackContainerName {this, "trackContainerName", "InDetTrackParticles"}

private

IDTtacks.

Definition at line 215 of file IDPerfMonZmumu.h.

{this, "trackContainerName", "InDetTrackParticles"};

m_trackParticleName

std::string IDPerfMonZmumu::m_trackParticleName

private

Track(Particle)TruthCollection input name.

Definition at line 219 of file IDPerfMonZmumu.h.

m_TrackRefitter1

ToolHandle<IegammaTrkRefitterTool> IDPerfMonZmumu::m_TrackRefitter1

private

The track refitter.

Definition at line 132 of file IDPerfMonZmumu.h.

m_TrackRefitter2

ToolHandle<IegammaTrkRefitterTool> IDPerfMonZmumu::m_TrackRefitter2

private

The track refitter.

Definition at line 135 of file IDPerfMonZmumu.h.

m_trackToVertexIPEstimator

ToolHandle< Trk::ITrackToVertexIPEstimator > IDPerfMonZmumu::m_trackToVertexIPEstimator

private

Needed for IP resolution studies.

Definition at line 152 of file IDPerfMonZmumu.h.

m_trackToVertexTool

ToolHandle<Reco::ITrackToVertex> IDPerfMonZmumu::m_trackToVertexTool

private

tool to extrapolate tracks to BL

Definition at line 138 of file IDPerfMonZmumu.h.

m_triggerDecision

ToolHandle<Trig::TrigDecisionTool> IDPerfMonZmumu::m_triggerDecision

private

Definition at line 141 of file IDPerfMonZmumu.h.

m_triggerMatching

ToolHandle<Trig::IMatchingTool> IDPerfMonZmumu::m_triggerMatching

private

Definition at line 145 of file IDPerfMonZmumu.h.

m_triggerName

std::string IDPerfMonZmumu::m_triggerName

private

Definition at line 231 of file IDPerfMonZmumu.h.

m_triggerPrescale

int IDPerfMonZmumu::m_triggerPrescale {}

private

Definition at line 230 of file IDPerfMonZmumu.h.

{};

m_Trk2VtxAssociationTool

ToolHandle< CP::ITrackVertexAssociationTool > IDPerfMonZmumu::m_Trk2VtxAssociationTool

private

Definition at line 209 of file IDPerfMonZmumu.h.

m_Truth_d0

std::vector<float> IDPerfMonZmumu::m_Truth_d0

private

Definition at line 389 of file IDPerfMonZmumu.h.

m_Truth_eta

std::vector<float> IDPerfMonZmumu::m_Truth_eta

private

Definition at line 377 of file IDPerfMonZmumu.h.

m_Truth_parent

std::vector<int> IDPerfMonZmumu::m_Truth_parent

private

Definition at line 403 of file IDPerfMonZmumu.h.

m_Truth_phi

std::vector<float> IDPerfMonZmumu::m_Truth_phi

private

Definition at line 383 of file IDPerfMonZmumu.h.

m_Truth_pt

std::vector<float> IDPerfMonZmumu::m_Truth_pt

private

Definition at line 371 of file IDPerfMonZmumu.h.

m_Truth_qoverp

std::vector<float> IDPerfMonZmumu::m_Truth_qoverp

private

Definition at line 401 of file IDPerfMonZmumu.h.

m_Truth_z0

std::vector<float> IDPerfMonZmumu::m_Truth_z0

private

Definition at line 395 of file IDPerfMonZmumu.h.

m_truthLinkVecName

std::string IDPerfMonZmumu::m_truthLinkVecName

private

TrackParticle input name.

Definition at line 220 of file IDPerfMonZmumu.h.

m_truthName

std::string IDPerfMonZmumu::m_truthName

private

Definition at line 218 of file IDPerfMonZmumu.h.

m_truthTree

TTree* IDPerfMonZmumu::m_truthTree

private

Definition at line 198 of file IDPerfMonZmumu.h.

m_truthTreeFolder

std::string IDPerfMonZmumu::m_truthTreeFolder

private

Definition at line 187 of file IDPerfMonZmumu.h.

m_truthTreeName

std::string IDPerfMonZmumu::m_truthTreeName

private

Definition at line 171 of file IDPerfMonZmumu.h.

m_useCustomMuonSelector

bool IDPerfMonZmumu::m_useCustomMuonSelector {}

private

Definition at line 126 of file IDPerfMonZmumu.h.

{};

m_useTrackSelectionTool

bool IDPerfMonZmumu::m_useTrackSelectionTool {}

private

Definition at line 121 of file IDPerfMonZmumu.h.

{};

m_UseTrigger

bool IDPerfMonZmumu::m_UseTrigger {true}

private

Definition at line 109 of file IDPerfMonZmumu.h.

{true};

m_validationMode

bool IDPerfMonZmumu::m_validationMode

private

< boolean to switch to validation mode

validation tree name - to be acessed by this from root

Definition at line 163 of file IDPerfMonZmumu.h.

m_ValidationTreeDescription

std::string IDPerfMonZmumu::m_ValidationTreeDescription

private

< validation tree description - second argument in TTree

stream/folder to for the TTree to be written out

Definition at line 179 of file IDPerfMonZmumu.h.

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vertexKey

SG::ReadHandleKey<xAOD::VertexContainer> IDPerfMonZmumu::m_vertexKey { this, "VertexContainer", "PrimaryVertices", "primary vertex container" }

private

Definition at line 212 of file IDPerfMonZmumu.h.

{ this, "VertexContainer", "PrimaryVertices", "primary vertex container" };

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

m_xZmm

ZmumuEvent IDPerfMonZmumu::m_xZmm

private

Definition at line 107 of file IDPerfMonZmumu.h.

m_Z0GapCut

double IDPerfMonZmumu::m_Z0GapCut {5.}

private

Definition at line 118 of file IDPerfMonZmumu.h.

{5.};

---

The documentation for this class was generated from the following files:

• IDPerfMonZmumu.h
• IDPerfMonZmumu.cxx