ZdcNtuple Class Reference

#include <ZdcNtuple.h>

Inheritance diagram for ZdcNtuple:

Public Types
typedef ServiceHandle< StoreGateSvc > &	MetaStorePtr_t
	Type of the metadata store pointer in standalone mode.
typedef const ServiceHandle< StoreGateSvc > &	ConstMetaStorePtr_t
typedef ServiceHandle< StoreGateSvc >	MetaStore_t
	Type of the metadata store variable in Athena.

Public Member Functions
void	setupTriggerHistos ()
	ZdcNtuple (const std::string &name, ISvcLocator *pSvcLocator)
void	processEventInfo ()
void	processVInjInfo ()
bool	processTriggerDecision ()
uint32_t	acceptEvent ()
void	processZdcNtupleFromModules ()
void	processMCEventCollection ()
void	processFCal ()
void	processMBTS ()
void	processInDet ()
void	writeTrack (const xAOD::TrackParticle *, const xAOD::Vertex *vertex, int)
int	trackQuality (const xAOD::TrackParticle *tp, const xAOD::Vertex *vertex)
void	processClusters ()
void	reprocessZdcModule (const xAOD::ZdcModule *zdcMod, bool flipdelay=0)
void	processTriggerTowers ()
void	processGaps ()
void	processProtons ()
double	dR (const double eta1, const double phi1, const double eta2, const double phi2)
virtual StatusCode	initialize () override
virtual StatusCode	execute () override
virtual StatusCode	finalize () override
::StatusCode	requestFileExecute ()
	register this algorithm to have an implementation of fileexecute
::StatusCode	requestBeginInputFile ()
	register this algorithm to have an implementation of beginInputFile
::StatusCode	requestEndInputFile ()
	register this algorithm to have an implementation of endInputFile
void	handle (const Incident &inc)
	receive the given incident
virtual StatusCode	sysInitialize ()
	Initialization method invoked by the framework.
const ServiceHandle< ITHistSvc > &	histSvc () const
	The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a pointer to the THistSvc.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
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
ConstMetaStorePtr_t	inputMetaStore () const
MetaStorePtr_t	inputMetaStore ()
ConstMetaStorePtr_t	outputMetaStore () const
MetaStorePtr_t	outputMetaStore ()

Public Attributes
bool	slimmed
bool	useGRL
std::string	grlFilename
bool	enableOutputTree
bool	enableOutputSamples
bool	enableTrigger
bool	writeOnlyTriggers
bool	enableID
bool	enableCalo
bool	enableClusters
bool	enableTracks
bool	enableMuons
bool	enableElectrons
bool	enablePhotons
bool	enableTT
bool	enableTruth
bool	enableJets
bool	enableTriggerJets
bool	zdcCalib
bool	zdcLaser
bool	zdcInj
bool	zdcOnly
unsigned int	zdcLowGainMode
size_t	trackLimit
bool	trackLimitReject
bool	flipDelay
bool	reprocZdc
std::string	auxSuffix
size_t	nsamplesZdc
bool	lhcf2022
bool	lhcf2022zdc
bool	lhcf2022afp
bool	pbpb2023
bool	enableZDC
bool	enableRPD
bool	enableRPDAmp
bool	enableCentroid
bool	doZdcCalib
std::string	zdcConfig
std::string	outputName
int	outputTreeScaledown
PublicToolHandle< Trig::TrigDecisionTool >	m_trigDecisionTool { this, "TrigDecisionTool", "", "Handle to the TrigDecisionTool" }
asg::AnaToolHandle< IGoodRunsListSelectionTool >	m_grl
asg::AnaToolHandle< ZDC::IZdcAnalysisTool >	m_zdcAnalysisTool
ToolHandle< InDet::IInDetTrackSelectionTool >	m_selTool
SG::ReadHandleKey< xAOD::ZdcModuleContainer >	m_zdcModuleContainerName { this, "ZdcModuleContainerName", "ZdcModules", "" }
SG::ReadHandleKey< xAOD::ZdcModuleContainer >	m_zdcSumContainerName { this, "ZdcSumContainerName", "ZdcSums", "" }
SG::ReadHandleKey< McEventCollection >	m_mcEventCollectionName {this, "MCEventCollectionName", "TruthEvent", ""}
const xAOD::EventInfo *	m_eventInfo
const xAOD::TrigDecision *	m_trigDecision
const xAOD::HIEventShapeContainer *	m_caloSums
const xAOD::HIEventShapeContainer *	m_eventShapes
const xAOD::ForwardEventInfoContainer *	m_mbtsInfo
const xAOD::MBTSModuleContainer *	m_mbtsModules
const xAOD::TrigT2MbtsBitsContainer *	m_trigT2MbtsBits
const xAOD::VertexContainer *	m_primaryVertices
const xAOD::CaloClusterContainer *	m_caloClusters
const xAOD::TrackParticleContainer *	m_trackParticles
const xAOD::EnergySumRoI *	m_lvl1EnergySumRoI
const xAOD::TruthParticleContainer *	m_truthParticleContainer
const xAOD::TriggerTowerContainer *	m_TTcontainer
const xAOD::AFPProtonContainer *	m_afpProtons
std::shared_ptr< ZdcInjPulserAmpMap >	m_zdcInjPulserAmpMap
int	m_nTriggers
int	m_eventCounter
bool	m_isMC
bool	m_setupTrigHist
int	m_scaledownCounter
std::vector< const Trig::ChainGroup * >	m_chainGroups
std::vector< const Trig::ChainGroup * >	m_rerunChainGroups
float	t_L1ET
float	t_L1ET24
TH1 *	h_zdcTriggers
TH1 *	h_zdcTriggersTBP
TTree *	m_outputTree
uint32_t	t_runNumber
uint32_t	t_eventNumber
uint32_t	t_lumiBlock
uint32_t	t_bcid
float	t_vInj
uint8_t	t_bunchGroup
uint32_t	t_passBits
uint32_t	t_extendedLevel1ID
uint32_t	t_timeStamp
uint32_t	t_timeStampNSOffset
float	t_avgIntPerCrossing
float	t_actIntPerCrossing
uint8_t	t_zdcEventInfoError
uint32_t	t_zdcEventInfoErrorWord
uint8_t	t_zdcDecodingError
uint8_t	t_rpdDecodingError
uint64_t	t_trigger
uint32_t	t_trigger_TBP
float	t_prescales [200]
bool	t_decisions [200]
bool	t_rerunDecisions [200]
float	t_mbts_in_e [2][8]
float	t_mbts_out_e [2][4]
float	t_mbts_in_t [2][8]
float	t_mbts_out_t [2][4]
float	t_T2mbts_in_e [2][8]
float	t_T2mbts_out_e [2][4]
float	t_T2mbts_in_t [2][8]
float	t_T2mbts_out_t [2][4]
uint32_t	t_tav [16]
uint32_t	t_tbp [16]
float	t_ZdcAmp [2]
float	t_ZdcAmpErr [2]
float	t_ZdcEnergy [2]
float	t_ZdcEnergyErr [2]
float	t_ZdcNLEnergy [2]
float	t_ZdcNLEnergyErr [2]
float	t_ZdcTime [2]
short	t_ZdcStatus [2]
unsigned int	t_ZdcModuleMask
float	t_ZdcTrigEff [2]
unsigned short	t_ZdcLucrodTriggerSideAmp [2]
unsigned short	t_ZdcLucrodTriggerSideAmpLG [2]
float	t_ZdcTruthTotal [2]
float	t_ZdcTruthInvis [2]
float	t_ZdcTruthEM [2]
float	t_ZdcTruthNonEM [2]
float	t_ZdcTruthEscaped [2]
std::vector< float >	t_ZdcTruthParticlePosx
std::vector< float >	t_ZdcTruthParticlePosy
std::vector< float >	t_ZdcTruthParticlePosz
std::vector< float >	t_ZdcTruthParticleTime
std::vector< float >	t_ZdcTruthParticlePx
std::vector< float >	t_ZdcTruthParticlePy
std::vector< float >	t_ZdcTruthParticlePz
std::vector< float >	t_ZdcTruthParticleEnergy
std::vector< int >	t_ZdcTruthParticlePid
std::vector< int >	t_ZdcTruthParticleStatus
float	t_ZdcModuleAmp [2][4]
float	t_ZdcModuleTime [2][4]
float	t_ZdcModuleFitAmp [2][4]
float	t_ZdcModuleFitT0 [2][4]
float	t_ZdcModuleChisq [2][4]
unsigned int	t_ZdcModuleStatus [2][4]
float	t_ZdcModuleCalibAmp [2][4]
float	t_ZdcModuleCalibTime [2][4]
float	t_ZdcModuleAmpError [2][4]
float	t_ZdcModuleBkgdMaxFraction [2][4]
float	t_ZdcModuleMinDeriv2nd [2][4]
float	t_ZdcModulePresample [2][4]
float	t_ZdcModulePreSampleAmp [2][4]
unsigned short	t_ZdcLucrodTriggerAmp [2][4]
unsigned short	t_ZdcLucrodTriggerAmpLG [2][4]
float	t_ZdcModuleMaxADC [2][4]
float	t_ZdcModuleMaxADCHG [2][4]
float	t_ZdcModuleMaxADCLG [2][4]
float	t_ZdcModulePeakADCHG [2][4]
float	t_ZdcModulePeakADCLG [2][4]
float	t_ZdcModuleAmpLGRefit [2][4]
float	t_ZdcModuleAmpCorrLGRefit [2][4]
float	t_ZdcModuleT0LGRefit [2][4]
float	t_ZdcModuleT0SubLGRefit [2][4]
float	t_ZdcModuleChisqLGRefit [2][4]
float	t_ZdcModuleTruthTotal [2][7]
float	t_ZdcModuleTruthInvis [2][7]
float	t_ZdcModuleTruthEM [2][7]
float	t_ZdcModuleTruthNonEM [2][7]
float	t_ZdcModuleTruthEscaped [2][7]
unsigned int	t_ZdcModuleTruthNphotons [2][7]
float	t_RpdChannelBaseline [2][16]
float	t_RpdChannelPileupExpFitParams [2][16][2]
float	t_RpdChannelPileupStretchedExpFitParams [2][16][3]
float	t_RpdChannelPileupExpFitParamErrs [2][16][2]
float	t_RpdChannelPileupStretchedExpFitParamErrs [2][16][3]
float	t_RpdChannelPileupExpFitMSE [2][16]
float	t_RpdChannelPileupStretchedExpFitMSE [2][16]
float	t_RpdChannelAmplitude [2][16]
float	t_RpdChannelAmplitudeCalib [2][16]
float	t_RpdChannelMaxADC [2][16]
float	t_RpdChannelMaxADCCalib [2][16]
unsigned int	t_RpdChannelMaxSample [2][16]
unsigned int	t_RpdChannelStatus [2][16]
float	t_RpdChannelPileupFrac [2][16]
unsigned int	t_RpdSideStatus [2]
unsigned int	t_RpdModuleTruthNphotons [2][16]
bool	t_centroidDecorationsAvailable
char	t_centroidEventValid
unsigned int	t_centroidStatus [2]
float	t_RPDChannelSubtrAmp [2][16]
float	t_RPDSubtrAmpSum [2]
float	t_xCentroidPreGeomCorPreAvgSubtr [2]
float	t_yCentroidPreGeomCorPreAvgSubtr [2]
float	t_xCentroidPreAvgSubtr [2]
float	t_yCentroidPreAvgSubtr [2]
float	t_xCentroid [2]
float	t_yCentroid [2]
float	t_xRowCentroid [2][4]
float	t_yColCentroid [2][4]
float	t_reactionPlaneAngle [2]
float	t_cosDeltaReactionPlaneAngle
int	t_nvx
float	t_vx [3]
int	t_vxntrk
std::vector< int >	t_vx_trk_index
int	t_vxtype
int	t_pvindex
float	t_vxcov [6]
float	t_vxsumpt2
float	t_puvxz
int	t_puvxntrk
float	t_puvxsumpt
int	t_nstrong
int	t_vxnlooseprimary
int	t_vxnminbias
int	t_vxngoodmuon
int	t_nvtx
std::vector< int8_t >	t_vtx_type
std::vector< float >	t_vtx_x
std::vector< float >	t_vtx_y
std::vector< float >	t_vtx_z
std::vector< int16_t >	t_vtx_ntrk_all
std::vector< float >	t_vtx_sumpt2_all
std::vector< int16_t >	t_vtx_ntrk
std::vector< float >	t_vtx_sumpt2
std::vector< std::vector< int16_t > >	t_vtx_trk_index
float	t_fcalEt
float	t_fcalEtA
float	t_fcalEtC
float	t_fcalEtA_TT
float	t_fcalEtC_TT
float	t_fcalEtA_TTsum
float	t_fcalEtC_TTsum
float	t_totalEt
float	t_totalEt_TTsum
float	t_totalEt24
float	t_totalEt24_TTsum
float	t_edgeGapA
float	t_edgeGapC
float	m_gapPtMin
float	m_gapThresholds [98]
TH1 *	h_TCSigCut
uint16_t	t_mbts_countA
uint16_t	t_mbts_countC
float	t_mbts_timeA
float	t_mbts_timeC
float	t_mbts_timeDiff
uint16_t	t_T2mbts_countAin
uint16_t	t_T2mbts_countCin
uint16_t	t_raw7 [2][4][2][2][7]
uint16_t	t_raw15 [2][4][2][2][15]
uint16_t	t_raw24 [2][4][2][2][24]
uint16_t	t_raw32 [2][4][2][2][32]
uint16_t	t_raw40 [2][4][2][2][40]
uint16_t	t_rpdRaw [2][16][24]
uint16_t	t_rpdRaw32 [2][16][32]
uint16_t	t_rpdRaw40 [2][16][40]
uint32_t	t_ntrk
std::vector< float >	t_trk_pt
std::vector< float >	t_trk_eta
std::vector< float >	t_trk_phi
std::vector< float >	t_trk_e
std::vector< float >	t_trk_theta
std::vector< float >	t_trk_d0
std::vector< float >	t_trk_z0
std::vector< float >	t_trk_vz
std::vector< float >	t_trk_vtxz
std::vector< int8_t >	t_trk_charge
std::vector< int16_t >	t_trk_quality
std::vector< int >	t_trk_index
std::vector< uint8_t >	t_trk_nPixHits
std::vector< uint8_t >	t_trk_nSctHits
std::vector< uint8_t >	t_trk_nPixDead
std::vector< uint8_t >	t_trk_nSctDead
std::vector< uint8_t >	t_trk_nPixHoles
std::vector< uint8_t >	t_trk_nSctHoles
std::vector< uint8_t >	t_trk_nTrtHits
std::vector< uint8_t >	t_trk_nTrtOutliers
std::vector< uint8_t >	t_trk_inPixHits
std::vector< uint8_t >	t_trk_exPixHits
std::vector< uint8_t >	t_trk_ninPixHits
std::vector< uint8_t >	t_trk_nexPixHits
std::vector< float >	t_trk_pixeldEdx
uint32_t	t_nclus
std::vector< float >	t_cc_pt
std::vector< float >	t_cc_eta
std::vector< float >	t_cc_phi
std::vector< float >	t_cc_e
std::vector< float >	t_cc_sig
std::vector< int >	t_cc_layer
std::vector< float >	t_cc_raw_m
std::vector< float >	t_cc_raw_eta
std::vector< float >	t_cc_raw_phi
std::vector< float >	t_cc_raw_e
std::vector< std::vector< float > >	t_cc_raw_samp
float	t_clusEt
float	t_clusEtMax
float	t_clusetaMax
float	t_clusphiMax
int	nProtons
std::vector< double >	proton_pt
std::vector< double >	proton_eta
std::vector< double >	proton_phi
std::vector< double >	proton_e
std::vector< int >	proton_side
std::vector< double >	proton_eLoss
std::vector< double >	proton_t
std::vector< std::vector< int > >	proton_track_stationID
std::vector< std::vector< float > >	proton_track_xLocal
std::vector< std::vector< float > >	proton_track_yLocal
std::vector< std::vector< float > >	proton_track_zLocal
std::vector< std::vector< float > >	proton_track_xSlope
std::vector< std::vector< float > >	proton_track_ySlope
std::vector< std::vector< int > >	proton_track_nClusters
TLorentzVector	p_beam
TLorentzVector	p_scat

Protected Member Functions
virtual void	print () const
	print the state of the algorithm
virtual::StatusCode	fileExecute ()
	perform the action exactly once for each file in the dataset
virtual::StatusCode	beginInputFile ()
	perform the action for the beginning of an input file
virtual::StatusCode	endInputFile ()
	perform the action for the end of an input file
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.
StatusCode	configAthHistogramming (const ServiceHandle< ITHistSvc > &histSvc, const std::string &prefix, const std::string &rootDir, const std::string &histNamePrefix, const std::string &histNamePostfix, const std::string &histTitlePrefix, const std::string &histTitlePostfix)
	To be called by the derived classes to fill the internal configuration.
TH1 *	bookGetPointer (const TH1 &hist, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of histograms.
TH1 *	bookGetPointer (TH1 *hist, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of histograms.
TH1 *	bookGetPointer (TH1 &histRef, std::string tDir="", std::string stream="")
	Simplify the booking and registering (into THistSvc) of histograms.
TTree *	bookGetPointer (const TTree &treeRef, std::string tDir="", std::string stream="")
	Simplify the booking and registering (into THistSvc) of TTrees.
TGraph *	bookGetPointer (const TGraph &graphRef, std::string tDir="", std::string stream="")
	Simplify the booking and registering (into THistSvc) of TGraphs.
TEfficiency *	bookGetPointer (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency *	bookGetPointer (TEfficiency *eff, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TEfficiency.
TEfficiency *	bookGetPointer (TEfficiency &effRef, std::string tDir="", std::string stream="")
	Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode	book (const TH1 &hist, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of histograms.
StatusCode	book (TH1 *hist, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of histograms.
StatusCode	book (TH1 &histRef, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of histograms.
StatusCode	book (const TTree &treeRef, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TTrees.
StatusCode	book (const TGraph &graphRef, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TGraphs.
StatusCode	book (const TEfficiency &eff, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode	book (TEfficiency *eff, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TEfficiency.
StatusCode	book (TEfficiency &effRef, const std::string &tDir="", const std::string &stream="")
	Simplify the booking and registering (into THistSvc) of TEfficiency.
TH1 *	hist (const std::string &histName, const std::string &tDir="", const std::string &stream="")
	Simplify the retrieval of registered histograms of any type.
TH2 *	hist2d (const std::string &histName, const std::string &tDir="", const std::string &stream="")
	Simplify the retrieval of registered 2-d histograms.
TH3 *	hist3d (const std::string &histName, const std::string &tDir="", const std::string &stream="")
	Simplify the retrieval of registered 3-d histograms.
TTree *	tree (const std::string &treeName, const std::string &tDir="", const std::string &stream="")
	Simplify the retrieval of registered TTrees.
TGraph *	graph (const std::string &graphName, const std::string &tDir="", const std::string &stream="")
	Simplify the retrieval of registered TGraphs.
TEfficiency *	efficiency (const std::string &effName, const std::string &tDir="", const std::string &stream="")
	Simplify the retrieval of registered TEfficiency.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t
typedef uint32_t	hash_t
	typedef for the internal hash
typedef std::map< const hash_t, TH1 * >	HistMap_t
	Typedef for convenience.
typedef std::map< const hash_t, TEfficiency * >	EffMap_t
	Typedef for convenience.
typedef std::map< const hash_t, TTree * >	TreeMap_t
	Typedef for convenience.
typedef std::map< const hash_t, TGraph * >	GraphMap_t
	Typedef for convenience.

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>
void	buildBookingString (std::string &bookingString, std::string &histName, std::string &tDir, std::string &stream, bool usePrefixPostfix=false)
	Method to build individual booking string.
void	myReplace (std::string &str, const std::string &oldStr, const std::string &newStr)
	Helper method to replace sub-string.
hash_t	hash (const std::string &histName) const
	Method to calculate a 32-bit hash from a string.

Private Attributes
unsigned int	m_lastRunNumber {0}
ZdcInjPulserAmpMap::Token	m_injMapRunToken {}
MetaStore_t	m_inputMetaStore
	Object accessing the input metadata store.
MetaStore_t	m_outputMetaStore
	Object accessing the output metadata store.
bool	m_hasFileExecute {false}
	the value of hasFileExecute
bool	m_hasBeginInputFile {false}
	the value of hasBeginInputFile
bool	m_hasEndInputFile {false}
	the value of hasEndInputFile
ServiceHandle< ITHistSvc >	m_histSvc
	Default constructor: AthHistogramAlgorithm();.
std::string	m_prefix
	Name of the ROOT output stream (file)
std::string	m_rootDir
	Name of the ROOT directory.
std::string	m_histNamePrefix
	The prefix for the histogram THx name.
std::string	m_histNamePostfix
	The postfix for the histogram THx name.
std::string	m_histTitlePrefix
	The prefix for the histogram THx title.
std::string	m_histTitlePostfix
	The postfix for the histogram THx title.
DataObjIDColl	m_extendedExtraObjects
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
HistMap_t	m_histMap
	The map of histogram names to their pointers.
EffMap_t	m_effMap
	The map of histogram names to their pointers.
TreeMap_t	m_treeMap
	The map of TTree names to their pointers.
GraphMap_t	m_graphMap
	The map of TGraph names to their pointers.
std::string	m_streamName
	Name of the ROOT output stream (file)
std::string	m_name
	Instance name.
MsgStream	m_msg
	Cached Message Stream.

Detailed Description

Definition at line 47 of file ZdcNtuple.h.

Member Typedef Documentation

ConstMetaStorePtr_t

typedef const ServiceHandle< StoreGateSvc >& EL::AnaAlgorithm::ConstMetaStorePtr_t

inherited

Definition at line 111 of file AnaAlgorithm.h.

EffMap_t

typedef std::map< const hash_t, TEfficiency* > AthHistogramming::EffMap_t

privateinherited

Typedef for convenience.

Definition at line 205 of file AthHistogramming.h.

GraphMap_t

typedef std::map< const hash_t, TGraph* > AthHistogramming::GraphMap_t

privateinherited

Typedef for convenience.

Definition at line 219 of file AthHistogramming.h.

hash_t

typedef uint32_t AthHistogramming::hash_t

privateinherited

typedef for the internal hash

Definition at line 169 of file AthHistogramming.h.

HistMap_t

typedef std::map< const hash_t, TH1* > AthHistogramming::HistMap_t

privateinherited

Typedef for convenience.

Definition at line 198 of file AthHistogramming.h.

MetaStore_t

typedef ServiceHandle< StoreGateSvc > EL::AnaAlgorithm::MetaStore_t

inherited

Type of the metadata store variable in Athena.

Definition at line 546 of file AnaAlgorithm.h.

MetaStorePtr_t

typedef ServiceHandle< StoreGateSvc >& EL::AnaAlgorithm::MetaStorePtr_t

inherited

Type of the metadata store pointer in standalone mode.

Definition at line 110 of file AnaAlgorithm.h.

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

TreeMap_t

typedef std::map< const hash_t, TTree* > AthHistogramming::TreeMap_t

privateinherited

Typedef for convenience.

Definition at line 212 of file AthHistogramming.h.

Constructor & Destructor Documentation

ZdcNtuple()

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

Definition at line 25 of file ZdcNtuple.cxx.

  : EL::AnaAlgorithm(name, pSvcLocator),
    m_grl ("GoodRunsListSelectionTool/grl", this),
    //m_zdcAnalysisTool("ZDC::ZdcAnalysisTool/ZdcAnalysisTool", this),
    m_selTool( "InDet::InDetTrackSelectionTool/TrackSelectionTool", this )
{
  // Here you put any code for the base initialization of variables,
  // e.g. initialize all pointers to 0.  Note that you should only put
  // the most basic initialization here, since this method will be
  // called on both the submission and the worker node.  Most of your
  // initialization code will go into histInitialize() and
  // initialize().
  m_setupTrigHist = false;
  m_eventCounter = 0;
 
  declareProperty("isMC",  m_isMC = false, "MC mode");
  declareProperty("enableOutputTree",  enableOutputTree = false, "Enable output tree");
  declareProperty("enableOutputSamples",  enableOutputSamples = false, "Output ZDC and RPD raw data");
  declareProperty("enableTrigger",  enableTrigger = true, "comment");
  declareProperty("enableTracks",  enableTracks = false, "comment");
  declareProperty("trackLimit",  trackLimit = 500, "comment");
  declareProperty("enableID",  enableID = false, "turn on to enable ID tracks & vertices for physics streams");
  declareProperty("enableCalo",  enableCalo = false, "turn on to enable calorimeter energy info for physics streams");
  declareProperty("enableClusters",  enableClusters = false, "turn on to enable calo topo cluster info for physics streams");
  declareProperty("writeOnlyTriggers",  writeOnlyTriggers = false, "comment");
  declareProperty("useGRL",  useGRL = true, "comment");
  declareProperty("grlFilename",  grlFilename = "$ROOTCOREBIN/data/ZdcNtuple/data16_hip8TeV.periodAllYear_DetStatus-v86-pro20-19_DQDefects-00-02-04_PHYS_HeavyIonP_All_Good.xml", "comment");
  declareProperty("slimmed",  slimmed = false, "comment");
  declareProperty("zdcCalib",  zdcCalib = false, "zdc/ZDCCalib file");
  declareProperty("zdcInj",  zdcInj = false, "ZDC injected pulse");
  declareProperty("zdcLaser",  zdcLaser = false, "Run 2 ZDC Laser");
  declareProperty("zdcOnly", zdcOnly = false, "comment");
  declareProperty("zdcLowGainMode",  zdcLowGainMode = 0, "comment");
 
  declareProperty("flipDelay",  flipDelay = 0, "comment");
  declareProperty("reprocZdc",  reprocZdc = 0, "comment");
  declareProperty("auxSuffix",  auxSuffix = "", "comment");
  declareProperty("nsamplesZdc",  nsamplesZdc = 24, "number of samples, 7 = most of Run 2, 24,32,40 = Run 3");
  declareProperty("lhcf2022", lhcf2022 = false,"LHCf2022 general config");
  declareProperty("lhcf2022afp", lhcf2022afp = false,"LHCf2022 AFP-specific config");
  declareProperty("lhcf2022zdc", lhcf2022zdc = false,"LHCf2022 ZDC-specific config");
  declareProperty("pbpb2023", pbpb2023 = true, "PbPb2023 config");
  declareProperty("zdcConfig", zdcConfig = "PbPb2018", "argument to configure ZdcAnalysisTool");
  declareProperty("doZdcCalib", doZdcCalib = false, "perform ZDC energy calibration");
  declareProperty("enableZDC", enableZDC = true);
  declareProperty("enableRPD",enableRPD = false,"enable reading any RPD decorations");
  declareProperty("enableRPDAmp",enableRPDAmp = false,"enable reading RPD amplitudes (also requires enableRPD)");
  declareProperty("enableCentroid",enableCentroid = false,"enable reading centroid decorations (also requires enableRPD)");
 
  declareProperty( "TrackSelectionTool", m_selTool );
 
  m_zdcAnalysisTool.declarePropertyFor (this, "zdcAnalysisTool");
  
}

Member Function Documentation

acceptEvent()

uint32_t ZdcNtuple::acceptEvent

(

)

Definition at line 1986 of file ZdcNtuple.cxx.

{
  uint32_t passbits = 0;
 
  if (!m_isMC)
  {
    if (useGRL)
    {
      if (!m_grl->passRunLB(*m_eventInfo)) {
        passbits += 1; // UPC GRL
      }
    }
 
    /*
      if(!m_grl_mb->passRunLB(*m_eventInfo)){
      passbits += 4; // MB GRL
      }
    */
 
    if (   (m_eventInfo->errorState(xAOD::EventInfo::LAr) == xAOD::EventInfo::Error )
           || (m_eventInfo->errorState(xAOD::EventInfo::Tile) == xAOD::EventInfo::Error )
           || (m_eventInfo->errorState(xAOD::EventInfo::SCT) == xAOD::EventInfo::Error )
           || (m_eventInfo->isEventFlagBitSet(xAOD::EventInfo::Core, 18) ) )
    {
      passbits += 2;
    } // end if event flags check
  } // end if the event is data
 
  return passbits;
}

beginInputFile()

StatusCode EL::AnaAlgorithm::beginInputFile ( )

protectedinherited

perform the action for the beginning of an input file

Ideally you don't use this, but instead rely on meta-data tools instead. However, there are enough people asking for it that I decided to implement it anyways.

\warn To use this you have to call requestBeginInputFile to use this.

\warn If a file is split across multiple jobs this will be called more than once. This only happens for specific batch drivers and/or if it is explicitly configured by the user. With PROOF it could even happen multiple times within the same job, and while PROOF is no longer supported that behavior may come back if support for a similar framework is added in the future. As such, this method should not be used for accounting that relies to be called exactly once per file, take a look at fileExecute if you want something that is guaranteed to be executed exactly once per input file.

\warn The execution order of beginInputFile and fileExecute is currently unspecified.

Definition at line 350 of file AnaAlgorithm.cxx.

  {
    return StatusCode::SUCCESS;
  }

book() [1/8]

StatusCode AthHistogramming::book ( const TEfficiency & eff, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 335 of file AthHistogramming.h.

{
  // We need to create a non-const clone
  TEfficiency* effClone = dynamic_cast< TEfficiency* >( eff.Clone() );
  if ( !effClone ) {
    m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone" << endmsg;
    return StatusCode::FAILURE;
  }
  return this->book( *effClone, tDir, stream );
}

book() [2/8]

StatusCode AthHistogramming::book ( const TGraph & graphRef, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TGraphs.

book() [3/8]

StatusCode AthHistogramming::book ( const TH1 & hist, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 303 of file AthHistogramming.h.

{
  // We need to create a non-const clone
  TH1* histClone = dynamic_cast< TH1* >( hist.Clone() );
  if ( !histClone ) {
    m_msg << MSG::ERROR << "Couldn't create a TH1 clone" << endmsg;
    return StatusCode::FAILURE;
  }
  return this->book( *histClone, tDir, stream );
}

book() [4/8]

StatusCode AthHistogramming::book ( const TTree & treeRef, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TTrees.

Definition at line 403 of file AthHistogramming.h.

{
  // Call the other Book method and see if it returns a valid pointer
  TTree* treePointer = this->bookGetPointer( treeRef, tDir, stream );
  if ( treePointer )
    {
      return StatusCode::SUCCESS;
    }
  else
    {
      return StatusCode::FAILURE;
    }
}

book() [5/8]

StatusCode AthHistogramming::book ( TEfficiency & effRef, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 356 of file AthHistogramming.h.

{
  // Call the other Book method and see if it returns a valid pointer
  TEfficiency* effPointer = this->bookGetPointer( effRef, tDir, stream );
  if ( !effPointer ) {
    m_msg << MSG::ERROR << "Couldn't book a TEfficiency" << endmsg;
    return StatusCode::FAILURE;
  }
  return StatusCode::SUCCESS;
}

book() [6/8]

StatusCode AthHistogramming::book ( TEfficiency * eff, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 346 of file AthHistogramming.h.

{
  if ( !eff ) {
    m_msg << MSG::ERROR << "Got a zero pointer to a TEfficiency" << endmsg;
    return StatusCode::FAILURE;
  }
  return this->book( *eff, tDir, stream );
}

book() [7/8]

StatusCode AthHistogramming::book ( TH1 & histRef, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 324 of file AthHistogramming.h.

{
  // Call the other Book method and see if it returns a valid pointer
  TH1* histPointer = this->bookGetPointer( histRef, tDir, stream );
  if ( !histPointer ) {
    m_msg << MSG::ERROR << "Couldn't book a TH1" << endmsg;
    return StatusCode::FAILURE;
  }
  return StatusCode::SUCCESS;
}

book() [8/8]

StatusCode AthHistogramming::book ( TH1 * hist, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 314 of file AthHistogramming.h.

{
  if ( !hist ) {
    m_msg << MSG::ERROR << "Got a zero pointer to a TH1" << endmsg;
    return StatusCode::FAILURE;
  }
  return this->book( *hist, tDir, stream );
}

bookGetPointer() [1/8]

TEfficiency * AthHistogramming::bookGetPointer ( const TEfficiency & eff, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 281 of file AthHistogramming.h.

{
  // We need to create a non-const clone
  TEfficiency* histClone = dynamic_cast< TEfficiency* >( hist.Clone() );
  if ( !histClone ) {
    m_msg << MSG::ERROR << "Couldn't create a TEfficiency clone in bookGetPointer" << endmsg;
    return 0;
  }
  return this->bookGetPointer( *histClone, tDir, stream );
 
}

bookGetPointer() [2/8]

TGraph * AthHistogramming::bookGetPointer ( const TGraph & graphRef, std::string tDir = "", std::string stream = "" )

protectedinherited

Simplify the booking and registering (into THistSvc) of TGraphs.

Definition at line 427 of file AthHistogramming.cxx.

{
  // Get a pointer
  const TGraph* graphPointer = &graphRef;
 
  // Check that we got a valid pointer
  if ( !graphPointer )
    {
      m_msg << MSG::WARNING
            << "We got an invalid TGraph pointer in the BookGetPointer(TGraph*) method of the class" << m_name
            << "!" << endmsg;
      return NULL;
    }
 
  // Modify the name and title according to the prefixes of this classes instance
  std::string graphName  = graphPointer->GetName();
  const std::string graphTitle = graphPointer->GetTitle();
 
  // Check if the hash for this graphName already exists, i.e., if we have a hash collision
  const hash_t graphHash = this->hash(graphName);
  GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
  if ( it != m_graphMap.end() ) // It does exist!
    {
      m_msg << MSG::WARNING
            << "Detected a hash collision. The hash for the TGraph with name=" << graphName
            << " already exists and points to a TGraph with name=" << it->second->GetName()
            << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
      return NULL;
    }
 
  // Create a clone that has the new name
  TGraph* graphClone = dynamic_cast< TGraph* >( graphPointer->Clone((m_histNamePrefix+graphName+m_histNamePostfix).c_str()) );
  if( !graphClone )
    {
      m_msg << MSG::WARNING
            << "We couldn't clone the TGraph in the BookGetPointer(TGraph&) method of the class" << m_name
            << "!" << endmsg;
      return NULL;
    }
  graphClone->SetTitle ((m_histTitlePrefix+graphTitle+m_histTitlePostfix).c_str());
 
  // Massage the final string to book things
  std::string bookingString("");
  this->buildBookingString( bookingString, graphName, tDir, stream );
 
  // Register the TGraph into the THistSvc
  if ( !((histSvc()->regGraph(bookingString, graphClone)).isSuccess()) )
    {
      m_msg << MSG::WARNING
            << "Problem registering TGraph with name " << graphName
            << ", title " << graphTitle
            << " in " << m_name << "!" << endmsg;
      return NULL;
    }
 
  // Also register it in the local map of string to pointer
  m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphClone ) );
 
  return graphClone;
}

bookGetPointer() [3/8]

TH1 * AthHistogramming::bookGetPointer ( const TH1 & hist, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 260 of file AthHistogramming.h.

{
  // We need to create a non-const clone
  TH1* histClone = dynamic_cast< TH1* >( hist.Clone() );
  if ( !histClone ) {
    m_msg << MSG::ERROR << "Couldn't create a TH1 clone in bookGetPointer" << endmsg;
    return 0;
  }
  return this->bookGetPointer( *histClone, tDir, stream );
 
}

bookGetPointer() [4/8]

TTree * AthHistogramming::bookGetPointer ( const TTree & treeRef, std::string tDir = "", std::string stream = "" )

protectedinherited

Simplify the booking and registering (into THistSvc) of TTrees.

Definition at line 312 of file AthHistogramming.cxx.

{
  // Get a pointer
  const TTree* treePointer = &treeRef;
 
  // Check that we got a valid pointer
  if ( !treePointer )
    {
      m_msg << MSG::WARNING
            << "We got an invalid TTree pointer in the BookGetPointer(TTree*) method of the class" << m_name
            << "!" << endmsg;
      return NULL;
    }
 
  // Modify the name and title according to the prefixes of this classes instance
  std::string treeName  = treePointer->GetName();
  const std::string treeTitle = treePointer->GetTitle();
 
  // Check if the hash for this treeName already exists, i.e., if we have a hash collision
  const hash_t treeHash = this->hash(treeName);
  TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
  if ( it != m_treeMap.end() ) // It does exist!
    {
      m_msg << MSG::WARNING
            << "Detected a hash collision. The hash for the TTree with name=" << treeName
            << " already exists and points to a TTree with name=" << it->second->GetName()
            << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
      return NULL;
    }
 
  // Create a clone that has the new name
  TTree* treeClone = dynamic_cast< TTree* >( treePointer->Clone(treeName.c_str()) );
  if( !treeClone )
    {
      m_msg << MSG::WARNING
            << "We couldn't clone the TTree in the BookGetPointer(TTree&) method of the class" << m_name
            << "!" << endmsg;
      return NULL;
    }
  treeClone->SetTitle (treeTitle.c_str());
 
  // Massage the final string to book things
  std::string bookingString("");
  this->buildBookingString( bookingString, treeName, tDir, stream );
 
  // Register the TTree into the THistSvc
  if ( !((histSvc()->regTree(bookingString, treeClone)).isSuccess()) )
    {
      m_msg << MSG::WARNING
            << "Problem registering TTree with name " << treeName
            << ", title " << treeTitle
            << " in " << m_name << "!" << endmsg;
      return NULL;
    }
 
  // Also register it in the local map of string to pointer
  m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treeClone ) );
 
  return treeClone;
}

bookGetPointer() [5/8]

TEfficiency * AthHistogramming::bookGetPointer ( TEfficiency & effRef, std::string tDir = "", std::string stream = "" )

protectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 146 of file AthHistogramming.cxx.

{
  // Modify the name and title according to the prefixes of this classes instance
  std::string effName(effRef.GetName());
  const std::string effTitle(effRef.GetTitle());
  std::string bookingString("");
 
  this->buildBookingString( bookingString, effName, tDir, stream );
  effRef.SetTitle((m_histTitlePrefix+effTitle+m_histTitlePostfix).c_str() );
  effRef.SetName(effName.c_str());
 
  // Check if the hash for this effName already exists, i.e., if we have a hash collision
  const hash_t effHash = this->hash(effName);
  EffMap_t::const_iterator it = m_effMap.find( effHash );
  if ( it != m_effMap.end() ) // It does exist!
    {
      m_msg << MSG::WARNING
            << "Detected a hash collision. The hash for the TEfficiency with name=" << effName
            << " already exists and points to a TEfficiency with name=" << it->second->GetName()
            << " NOT going to book the new TEfficiency and returning a NULL pointer!" << endmsg;
      return NULL;
    }
 
  // Set the new name and title for the TEfficiency, based on the prefixes that the user set for this class instance
  // Create a clone that has the new name
 
  // Massage the final string to book things
 
  // Register the TEfficiency into the THistSvc
  if ( !((histSvc()->regEfficiency(bookingString, &effRef)).isSuccess()) )
    {
      m_msg << MSG::WARNING
            << "Problem registering TEfficiency with name " << effName
            << ", name prefix " << m_histNamePrefix
            << ", title " << effTitle
            << ", tile prefix " << m_histTitlePrefix
            << ", and tile postfix " << m_histTitlePostfix
            << " in " << m_name << "!" << endmsg;
      return NULL;
    }
 
  // Also register it in the local map of string to pointer
  m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, &effRef ) );
 
  return &effRef;
}

bookGetPointer() [6/8]

TEfficiency * AthHistogramming::bookGetPointer ( TEfficiency * eff, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of TEfficiency.

Definition at line 293 of file AthHistogramming.h.

{
  if ( !hist ) {
    m_msg << MSG::ERROR << "Got a zero pointer to a TEfficiency in bookGetPointer" << endmsg;
    return 0;
  }
  return this->bookGetPointer( *hist, tDir, stream );
}

bookGetPointer() [7/8]

TH1 * AthHistogramming::bookGetPointer ( TH1 & histRef, std::string tDir = "", std::string stream = "" )

protectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 98 of file AthHistogramming.cxx.

{
  // Modify the name and title according to the prefixes of this classes instance
  std::string histName(histRef.GetName());
  const std::string histTitle(histRef.GetTitle());
  std::string bookingString("");
 
  this->buildBookingString( bookingString, histName, tDir, stream );
  histRef.SetTitle((m_histTitlePrefix+histTitle+m_histTitlePostfix).c_str() );
  histRef.SetName(histName.c_str());
 
  // Check if the hash for this histName already exists, i.e., if we have a hash collision
  const hash_t histHash = this->hash(histName);
  HistMap_t::const_iterator it = m_histMap.find( histHash );
  if ( it != m_histMap.end() ) // It does exist!
    {
      m_msg << MSG::WARNING
            << "Detected a hash collision. The hash for the histogram with name=" << histName
            << " already exists and points to a histogram with name=" << it->second->GetName()
            << " NOT going to book the new histogram and returning a NULL pointer!" << endmsg;
      return NULL;
    }
 
  // Set the new name and title for the histogram, based on the prefixes that the user set for this class instance
  // Create a clone that has the new name
 
  // Massage the final string to book things
 
  // Register the histogram into the THistSvc
  if ( !((histSvc()->regHist(bookingString, &histRef)).isSuccess()) )
    {
      m_msg << MSG::WARNING
            << "Problem registering histogram with name " << histName
            << ", name prefix " << m_histNamePrefix
            << ", title " << histTitle
            << ", tile prefix " << m_histTitlePrefix
            << ", and tile postfix " << m_histTitlePostfix
            << " in " << m_name << "!" << endmsg;
      return NULL;
    }
 
  // Also register it in the local map of string to pointer
  m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, &histRef ) );
 
  return &histRef;
}

bookGetPointer() [8/8]

TH1 * AthHistogramming::bookGetPointer ( TH1 * hist, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the booking and registering (into THistSvc) of histograms.

Definition at line 272 of file AthHistogramming.h.

{
  if ( !hist ) {
    m_msg << MSG::ERROR << "Got a zero pointer to a TH1 in bookGetPointer" << endmsg;
    return 0;
  }
  return this->bookGetPointer( *hist, tDir, stream );
}

buildBookingString()

void AthHistogramming::buildBookingString ( std::string & bookingString, std::string & histName, std::string & tDir, std::string & stream, bool usePrefixPostfix = false )

privateinherited

Method to build individual booking string.

Definition at line 560 of file AthHistogramming.cxx.

{
  // Massage the final string to book things
  if(tDir.empty()) tDir = m_rootDir;
  size_t pos = histName.rfind('/');
  if(pos != std::string::npos){
    tDir+='/';
    tDir.append(histName, 0,pos);
    histName.erase(0,pos+1);
  };
  if(stream.empty()) stream = m_streamName;
 
  if(usePrefixPostfix){
    bookingString = "/"+stream+"/"+tDir+"/"+m_histNamePrefix+histName+m_histNamePostfix;
  } else {
    bookingString = "/"+stream+"/"+tDir+"/"+histName;
  }
  while(bookingString.find("//") != std::string::npos){
    this->myReplace(bookingString,"//","/");
  }
 
  return;
}

configAthHistogramming()

StatusCode AthHistogramming::configAthHistogramming ( const ServiceHandle< ITHistSvc > & histSvc, const std::string & prefix, const std::string & rootDir, const std::string & histNamePrefix, const std::string & histNamePostfix, const std::string & histTitlePrefix, const std::string & histTitlePostfix )

protectedinherited

To be called by the derived classes to fill the internal configuration.

Definition at line 66 of file AthHistogramming.cxx.

{
  m_histSvc          = histSvc;
  m_streamName       = prefix;
  m_rootDir          = rootDir;
  m_histNamePrefix   = histNamePrefix;
  m_histNamePostfix  = histNamePostfix;
  m_histTitlePrefix  = histTitlePrefix;
  m_histTitlePostfix = histTitlePostfix;
 
  return StatusCode::SUCCESS;
}

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< 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< 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< 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; }

dR()

double ZdcNtuple::dR	(	const double	eta1,
		const double	phi1,
		const double	eta2,
		const double	phi2 )

Definition at line 2139 of file ZdcNtuple.cxx.

{
  double deta = std::abs(eta1 - eta2);
  double dphi = std::abs(phi1 - phi2) < TMath::Pi() ? std::abs(phi1 - phi2) : 2 * TMath::Pi() - std::abs(phi1 - phi2);
  return std::sqrt(deta * deta + dphi * dphi);
}

efficiency()

TEfficiency * AthHistogramming::efficiency ( const std::string & effName, const std::string & tDir = "", const std::string & stream = "" )

protectedinherited

Simplify the retrieval of registered TEfficiency.

Definition at line 250 of file AthHistogramming.cxx.

{
  // Build a 32 bit hash out of the name
  const hash_t effHash = this->hash(effName);
 
  // See if this entry exists in the map
  EffMap_t::const_iterator it = m_effMap.find( effHash );
  if ( it == m_effMap.end() ) // It doesn't exist!
    { // Let's see into the THistSvc if somebody else has registered the TEfficiency...
 
      // Need to copy the strings as we will massage them from here on
      std::string effNameCopy = effName;
      std::string tDirCopy     = tDir;
      std::string streamCopy   = stream;
 
      // Massage the final string to book things
      std::string bookingString("");
      this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy ,false);
 
      TEfficiency* effPointer(NULL);
      if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
        {
          // Massage the final string to book things
          std::string bookingString("");
          this->buildBookingString( bookingString, effNameCopy, tDirCopy, streamCopy, true );
 
          if ( !((histSvc()->getEfficiency(bookingString, effPointer)).isSuccess()) )
            {
              m_msg << MSG::WARNING
                    << "Problem retrieving the TEfficiency with name (including pre- and post-fixes) "
                    << m_histNamePrefix + effNameCopy + m_histNamePostfix
                    << " or with name " << effNameCopy
                    << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
                    << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
              return NULL;
            }
          // If we get to here, we actually found the TEfficiency in the THistSvc.
          // So let's add it to the local cache map and return its pointer
          m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
          return effPointer;
        }
      // If we get to here, we actually found the TEfficiency in the THistSvc.
      // So let's add it to the local cache map and return its pointer
      m_effMap.insert( m_effMap.end(), std::pair< const hash_t, TEfficiency* >( effHash, effPointer ) );
      return effPointer;
    }
 
  // Return the pointer to the TEfficiency that we got from the local cache map
  return it->second;
}

endInputFile()

StatusCode EL::AnaAlgorithm::endInputFile ( )

protectedinherited

perform the action for the end of an input file

Ideally you don't use this, but instead rely on meta-data tools instead. However, there are enough people asking for it that I decided to implement it anyways.

\warn To use this you have to call requestEndInputFile to use this.

\warn If a file is split across multiple jobs this will be called more than once. This only happens for specific batch drivers and/or if it is explicitly configured by the user. With PROOF it could even happen multiple times within the same job, and while PROOF is no longer supported that behavior may come back if support for a similar framework is added in the future. As such, this method should not be used for accounting that relies to be called exactly once per file, take a look at fileExecute if you want something that is guaranteed to be executed exactly once per input file.

\warn The execution order of endInputFile and fileExecute is currently unspecified.

Definition at line 358 of file AnaAlgorithm.cxx.

  {
    return StatusCode::SUCCESS;
  }

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< 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 ZdcNtuple::execute ( )

overridevirtual

Definition at line 544 of file ZdcNtuple.cxx.

{
  // Here you do everything that needs to be done on every single
  // events, e.g. read input variables, apply cuts, and fill
  // histograms and trees.  This is where most of your actual analysis
  // code will go.
 
 
  // prefilter with a track limit
 
  if (!evtStore())
  {
    ANA_MSG_INFO("*** No event found! ***");
    return StatusCode::SUCCESS;
  }
 
  ANA_CHECK(evtStore()->retrieve( m_eventInfo, "EventInfo"));
  processEventInfo();
  if (zdcInj){
    processVInjInfo();
  }
 
  //tracks used to go here
 
  m_trackParticles = 0;
 
  if ((!(zdcCalib || zdcLaser || zdcOnly || zdcInj)) && enableID)
  {
    ANA_MSG_DEBUG("Trying to extract InDetTrackParticles from evtStore()=" << evtStore());
    ANA_CHECK(evtStore()->retrieve( m_trackParticles, "InDetTrackParticles") );
    size_t n = m_trackParticles->size();
    ANA_MSG_DEBUG("Done w/ extracting InDetTrackParticles with size = " << n);
 
    if (n > trackLimit && trackLimitReject)  return StatusCode::SUCCESS;
  }
 
  bool passTrigger = true;
  m_trigDecision = 0;
  if (enableTrigger)
  {
    ANA_CHECK(evtStore()->retrieve( m_trigDecision, "xTrigDecision"));
    if (!m_setupTrigHist) setupTriggerHistos();
    passTrigger = processTriggerDecision();
  }
 
  if (reprocZdc){
    ANA_MSG_INFO ("Reprocessing ZDC in ZdcNtuple");
    ANA_CHECK(m_zdcAnalysisTool->reprocessZdc());
  }else{
    ANA_MSG_DEBUG ("No ZDC reprocessing");
  }
 
  processZdcNtupleFromModules(); // same model in both cases -- processZdcNtuple() goes straight to the anlaysis tool, which is good for debugging
 
  if(m_isMC){
    processMCEventCollection();
  }
 
  if (!(zdcCalib || zdcLaser || zdcOnly || zdcInj))
  {
 
    // PLEASE NOTE: the commented sections here will be restored once we have a better sense of the Run 3 HI data
 
    // Global E_T quantities for centrality
    if (enableCalo){
        ANA_CHECK(evtStore()->retrieve( m_caloSums, "CaloSums") );
        ANA_CHECK(evtStore()->retrieve( m_eventShapes, "HIEventShape") );
 
        m_lvl1EnergySumRoI = 0;
        ANA_CHECK(evtStore()->retrieve( m_lvl1EnergySumRoI, "LVL1EnergySumRoI") );
 
        processFCal();
    } 
 
    // MBTS quantities, but may require a derivation to be accessible (required STDM6 in pp)
    //ANA_CHECK(evtStore()->retrieve( m_mbtsInfo, "MBTSForwardEventInfo") );
    //ANA_CHECK(evtStore()->retrieve( m_mbtsModules, "MBTSModules") );
    //m_trigT2MbtsBits = 0;
    //ANA_CHECK(evtStore()->retrieve( m_trigT2MbtsBits, "HLT_xAOD__TrigT2MbtsBitsContainer_T2Mbts") );
    //processMBTS();
 
    if (enableID){
        ANA_CHECK(evtStore()->retrieve( m_primaryVertices, "PrimaryVertices") );
        processInDet();     
    }
 
 
    if (enableClusters){
        ANA_CHECK(evtStore()->retrieve( m_caloClusters, "CaloCalTopoClusters"));
        processClusters();       
    }
 
    // Gaps will require some evaluation of Run 3 performance of the clusters
    //processGaps();
 
    if( lhcf2022||lhcf2022zdc||lhcf2022afp ){
      ANA_CHECK(evtStore()->retrieve( m_afpProtons, "AFPProtonContainer"));
      processProtons();
    }
  }
 
  // if trigger enabled, only write out events which pass one of them, unless using MC
 
  if (enableTrigger && !passTrigger && !m_isMC && writeOnlyTriggers) return StatusCode::SUCCESS;
 
 
  if (enableOutputTree)
  {
    tree( "zdcTree" )->Fill();
  }
 
  return StatusCode::SUCCESS;
}

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< 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 & AthAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 50 of file AthAlgorithm.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 Algorithm::extraOutputDeps();
}

fileExecute()

StatusCode EL::AnaAlgorithm::fileExecute ( )

protectedinherited

perform the action exactly once for each file in the dataset

Ideally you don't use this, but instead rely on meta-data tools instead. However, there are enough people asking for it that I decided to implement it anyways.

\warn To use this you have to call requestFileExecute to use this.

\warn The user should not expect this to be called at any particular point in execution. If a file is split between multiple jobs this will be called in only one of these jobs, and not the others. It usually gets called before the first event in a file, but that is not guaranteed and relying on this is a bug.

\warn The execution order of beginInputFile and fileExecute is currently unspecified.

\warn fileExecute does not work with sub-file splitting in Athena, i.e. processing half the events of a file in one job the other half in another job. this should not normally happen, unless you do crazy things like run AthenaMP or explicitly select sub-file splitting in panda. in that case you are on your own.

Definition at line 342 of file AnaAlgorithm.cxx.

  {
    return StatusCode::SUCCESS;
  }

finalize()

StatusCode ZdcNtuple::finalize ( )

overridevirtual

Definition at line 2147 of file ZdcNtuple.cxx.

{
  // This method is the mirror image of initialize(), meaning it gets
  // called after the last event has been processed on the worker node
  // and allows you to finish up any objects you created in
  // initialize() before they are written to disk.  This is actually
  // fairly rare, since this happens separately for each worker node.
  // Most of the time you want to do your post-processing on the
  // submission node after all your histogram outputs have been
  // merged.  This is different from histFinalize() in that it only
  // gets called on worker nodes that processed input events.
 
 
  return StatusCode::SUCCESS;
}

graph()

TGraph * AthHistogramming::graph ( const std::string & graphName, const std::string & tDir = "", const std::string & stream = "" )

protectedinherited

Simplify the retrieval of registered TGraphs.

Definition at line 492 of file AthHistogramming.cxx.

{
  // Build a 32 bit hash out of the name
  const hash_t graphHash = this->hash(graphName);
 
  // See if this entry exists in the map
  GraphMap_t::const_iterator it = m_graphMap.find( graphHash );
  if ( it == m_graphMap.end() ) // It doesn't exist!
    { // Let's see into the THistSvc if somebody else has registered the TGraph...
 
      // Need to copy the strings as we will massage them from here on
      std::string graphNameCopy = graphName;
      std::string tDirCopy      = tDir;
      std::string streamCopy    = stream;
 
      // Massage the final string to book things
      std::string bookingString("");
      this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, true);
 
      TGraph* graphPointer(NULL);
      if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
        {
          // Massage the final string to book things
          std::string bookingString("");
          this->buildBookingString( bookingString, graphNameCopy, tDirCopy, streamCopy, false );
 
          if ( !((histSvc()->getGraph(bookingString, graphPointer)).isSuccess()) )
            {
              m_msg << MSG::WARNING
                    << "Problem retrieving the TGraph with name (including pre- and post-fixes) "
                    << m_histNamePrefix + graphNameCopy + m_histNamePostfix
                    << " or with name " << graphNameCopy
                    << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
                    << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
              return NULL;
            }
          // If we get to here, we actually found the TGraph in the THistSvc.
          // So let's add it to the local cache map and return its pointer
          m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
          return graphPointer;
        }
      // If we get to here, we actually found the TGraph in the THistSvc.
      // So let's add it to the local cache map and return its pointer
      m_graphMap.insert( m_graphMap.end(), std::pair< const hash_t, TGraph* >( graphHash, graphPointer ) );
      return graphPointer;
    }
 
 
  // Return the pointer to the TGraph that we got from the local cache map
  return it->second;
}

handle()

void EL::AnaAlgorithm::handle ( const Incident & inc )

inherited

receive the given incident

Guarantee

basic

Failures

incident handling errors

Definition at line 520 of file AnaAlgorithm.cxx.

  {
    if (inc.type() == IncidentType::BeginInputFile)
    {
      if (m_hasBeginInputFile)
        ANA_CHECK_THROW (beginInputFile ());
      if (m_hasFileExecute)
        ANA_CHECK_THROW (fileExecute ());
    } else if (inc.type() == IncidentType::EndInputFile)
    {
      if (m_hasEndInputFile)
        ANA_CHECK_THROW (endInputFile ());
    } else
    {
      ATH_MSG_WARNING( "Unknown incident type received: " << inc.type() );
    }
  }

hash()

AthHistogramming::hash_t AthHistogramming::hash ( const std::string & histName ) const

inlineprivateinherited

Method to calculate a 32-bit hash from a string.

Definition at line 428 of file AthHistogramming.h.

{
  const uint64_t hash64 = CxxUtils::crc64( histName );
  return (hash_t)(hash64 & 0xFFFFFFFF);
}

hist()

TH1 * AthHistogramming::hist ( const std::string & histName, const std::string & tDir = "", const std::string & stream = "" )

protectedinherited

Simplify the retrieval of registered histograms of any type.

Definition at line 198 of file AthHistogramming.cxx.

{
  // Build a 32 bit hash out of the name
  const hash_t histHash = this->hash(histName);
 
  // See if this entry exists in the map
  HistMap_t::const_iterator it = m_histMap.find( histHash );
  if ( it == m_histMap.end() ) // It doesn't exist!
    { // Let's see into the THistSvc if somebody else has registered the histogram...
 
      // Need to copy the strings as we will massage them from here on
      std::string histNameCopy = histName;
      std::string tDirCopy     = tDir;
      std::string streamCopy   = stream;
 
      // Massage the final string to book things
      std::string bookingString("");
      this->buildBookingString( bookingString, histNameCopy, tDirCopy, streamCopy ,false);
 
      TH1* histPointer(NULL);
      if ( !((histSvc()->getHist(bookingString, histPointer)).isSuccess()) )
        {
          // Massage the final string to book things
          std::string bookingString("");
          this->buildBookingString( bookingString, histNameCopy, tDirCopy, streamCopy, true );
 
          if ( !((histSvc()->getHist(bookingString, histPointer)).isSuccess()) )
            {
              m_msg << MSG::WARNING
                    << "Problem retrieving the histogram with name (including pre- and post-fixes) "
                    << m_histNamePrefix + histNameCopy + m_histNamePostfix
                    << " or with name " << histNameCopy
                    << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
                    << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
              return NULL;
            }
          // If we get to here, we actually found the histogram in the THistSvc.
          // So let's add it to the local cache map and return its pointer
          m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, histPointer ) );
          return histPointer;
        }
      // If we get to here, we actually found the histogram in the THistSvc.
      // So let's add it to the local cache map and return its pointer
      m_histMap.insert( m_histMap.end(), std::pair< const hash_t, TH1* >( histHash, histPointer ) );
      return histPointer;
    }
 
 
  // Return the pointer to the histogram that we got from the local cache map
  return it->second;
}

hist2d()

TH2 * AthHistogramming::hist2d ( const std::string & histName, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the retrieval of registered 2-d histograms.

Definition at line 369 of file AthHistogramming.h.

{
  // Get the TH1 pointer
  TH1* th1Pointer = this->hist(histName, tDir, stream);
  if ( !th1Pointer )
    {
      m_msg << MSG::ERROR
            << "Cannot get a 2-d histogram with name " << histName
            << "... will probably seg-fault!" << endmsg;
      return NULL;
    }
  // If the TH1 pointer is valid, simply return the dynamic_cast
  return dynamic_cast<TH2*>( th1Pointer );
}

hist3d()

TH3 * AthHistogramming::hist3d ( const std::string & histName, const std::string & tDir = "", const std::string & stream = "" )

inlineprotectedinherited

Simplify the retrieval of registered 3-d histograms.

Definition at line 386 of file AthHistogramming.h.

{
  // Get the TH1 pointer
  TH1* th1Pointer = this->hist(histName, tDir, stream);
  if ( !th1Pointer )
    {
      m_msg << MSG::ERROR
            << "Cannot get a 3-d histogram with name " << histName
            << "... will probably seg-fault!" << endmsg;
      return NULL;
    }
  // If the TH1 pointer is valid, simply return the dynamic_cast
  return dynamic_cast<TH3*>( th1Pointer );
}

histSvc()

const ServiceHandle< ITHistSvc > & AthHistogramAlgorithm::histSvc ( ) const

inlineinherited

The standard THistSvc (for writing histograms and TTrees and more to a root file) Returns (kind of) a pointer to the THistSvc.

Definition at line 113 of file AthHistogramAlgorithm.h.

{
  return m_histSvc;
}

initialize()

StatusCode ZdcNtuple::initialize ( )

overridevirtual

Definition at line 99 of file ZdcNtuple.cxx.

{
  // Here you do everything that needs to be done at the very
  // beginning on each worker node, e.g. create histograms and output
  // trees.  This method gets called before any input files are
  // connected.
 
  ANA_MSG_DEBUG("Howdy from Initialize!");
 
  ANA_CHECK(m_zdcModuleContainerName.initialize());
  ANA_CHECK(m_zdcSumContainerName.initialize());
  ATH_CHECK(m_mcEventCollectionName.initialize());
 
  if (enableOutputTree)
  {
 
    ANA_CHECK( book(TTree("zdcTree", "ZDC Tree")));
    m_outputTree = tree( "zdcTree" );
 
    m_outputTree->Branch("runNumber", &t_runNumber, "runNumber/i");
    m_outputTree->Branch("eventNumber", &t_eventNumber, "eventNumber/i");
    m_outputTree->Branch("lumiBlock", &t_lumiBlock, "lumiBlock/i");
    m_outputTree->Branch("bunchGroup", &t_bunchGroup, "bunchGroup/b");
    m_outputTree->Branch("bcid", &t_bcid, "bcid/i");
    if (zdcInj) m_outputTree->Branch("vInj",&t_vInj,"vInj/F");
    m_outputTree->Branch("avgIntPerCrossing", &t_avgIntPerCrossing, "avgIntPerCrossing/F");
    m_outputTree->Branch("actIntPerCrossing", &t_actIntPerCrossing, "actIntPerCrossing/F");
    m_outputTree->Branch("trigger", &t_trigger, "trigger/l");
    m_outputTree->Branch("trigger_TBP", &t_trigger_TBP, "trigger_TBP/i");
    m_outputTree->Branch("tbp", &t_tbp, "tbp[16]/i");
    m_outputTree->Branch("tav", &t_tav, "tav[16]/i");
    m_outputTree->Branch("passBits", &t_passBits, "passBits/i");
    m_outputTree->Branch("extendedLevel1ID",&t_extendedLevel1ID,"extendedLevel1ID/i");
    m_outputTree->Branch("timeStamp",&t_timeStamp,"timeStamp/i");
    m_outputTree->Branch("timeStampNSOffset",&t_timeStampNSOffset,"timeStampNSOffset/i");
    m_outputTree->Branch("zdcEventInfoError",&t_zdcEventInfoError,"zdcEventInfoError/b");
    m_outputTree->Branch("zdcEventInfoErrorWord",&t_zdcEventInfoErrorWord,"zdcEventInfoErrorWord/i");
    // ZDC and RPD decoding errors
    m_outputTree->Branch("zdcDecodingError",&t_zdcDecodingError,"zdcDecodingError/b");
    m_outputTree->Branch("rpdDecodingError",&t_rpdDecodingError,"rpdDecodingError/b");
 
    if (enableZDC)
      {
    if (enableOutputSamples)
      {
        if (nsamplesZdc == 7)
          {
        ANA_MSG_INFO("Setting up for 7 samples");
        m_outputTree->Branch("zdc_raw", &t_raw7, "zdc_raw[2][4][2][2][7]/s"); // 7 samples
          }
        
        if (nsamplesZdc == 15)
          {
        ANA_MSG_INFO("Setting up for 15 samples");
        m_outputTree->Branch("zdc_raw", &t_raw15, "zdc_raw[2][4][2][2][15]/s"); // 15 samples
          }
        
        if (nsamplesZdc == 24)
          {
        ANA_MSG_INFO("Setting up for 24 samples");
        m_outputTree->Branch("zdc_raw", &t_raw24, "zdc_raw[2][4][2][2][24]/s"); // 24 samples
        m_outputTree->Branch("rpd_raw", &t_rpdRaw, "rpd_raw[2][16][24]/s"); // 24 samples
          }
        if (nsamplesZdc == 32)
          {
        ANA_MSG_INFO("Setting up for 32 samples");
        m_outputTree->Branch("zdc_raw", &t_raw32, "zdc_raw[2][4][2][2][32]/s"); // 32 samples
        m_outputTree->Branch("rpd_raw", &t_rpdRaw32, "rpd_raw[2][16][32]/s"); // 32 samples
          }
        if (nsamplesZdc == 40)
          {
        ANA_MSG_INFO("Setting up for 40 samples");
        m_outputTree->Branch("zdc_raw", &t_raw40, "zdc_raw[2][4][2][2][40]/s"); // 40 samples
        m_outputTree->Branch("rpd_raw", &t_rpdRaw40, "rpd_raw[2][16][40]/s"); // 40 samples
          }
      }
    
    m_outputTree->Branch("zdc_ZdcAmp", &t_ZdcAmp, "zdc_ZdcAmp[2]/F");
    m_outputTree->Branch("zdc_ZdcAmpErr", &t_ZdcAmpErr, "zdc_ZdcAmpErr[2]/F");
    m_outputTree->Branch("zdc_ZdcEnergy", &t_ZdcEnergy, "zdc_ZdcEnergy[2]/F");
    m_outputTree->Branch("zdc_ZdcEnergyErr", &t_ZdcEnergyErr, "zdc_ZdcEnergyErr[2]/F");
    m_outputTree->Branch("zdc_ZdcNLEnergy", &t_ZdcNLEnergy, "zdc_ZdcNLEnergy[2]/F");
    m_outputTree->Branch("zdc_ZdcNLEnergyErr", &t_ZdcNLEnergyErr, "zdc_ZdcNLEnergyErr[2]/F");
    m_outputTree->Branch("zdc_ZdcTime", &t_ZdcTime, "zdc_ZdcTime[2]/F");
    m_outputTree->Branch("zdc_ZdcStatus", &t_ZdcStatus, "zdc_ZdcStatus[2]/S");
    m_outputTree->Branch("zdc_ZdcTrigEff", &t_ZdcTrigEff, "zdc_ZdcTrigEff[2]/F");
    m_outputTree->Branch("zdc_ZdcModuleMask", &t_ZdcModuleMask, "zdc_ZdcModuleMask/i");
    m_outputTree->Branch("zdc_ZdcLucrodTriggerSideAmp",&t_ZdcLucrodTriggerSideAmp,"zdc_ZdcLucrodTriggerSideAmp[2]/S");
    m_outputTree->Branch("zdc_ZdcLucrodTriggerSideAmpLG",&t_ZdcLucrodTriggerSideAmpLG,"zdc_ZdcLucrodTriggerSideAmpLG[2]/S");
    
    m_outputTree->Branch("zdc_ZdcModuleAmp", &t_ZdcModuleAmp, "zdc_ZdcModuleAmp[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleTime", &t_ZdcModuleTime, "zdc_ZdcModuleTime[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleFitAmp", &t_ZdcModuleFitAmp, "zdc_ZdcModuleFitAmp[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleFitT0", &t_ZdcModuleFitT0, "zdc_ZdcModuleFitT0[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleStatus", &t_ZdcModuleStatus, "zdc_ZdcModuleStatus[2][4]/i");
    m_outputTree->Branch("zdc_ZdcModuleChisq", &t_ZdcModuleChisq, "zdc_ZdcModuleChisq[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleCalibAmp", &t_ZdcModuleCalibAmp, "zdc_ZdcModuleCalibAmp[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleCalibTime", &t_ZdcModuleCalibTime, "zdc_ZdcModuleCalibTime[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleBkgdMaxFraction", &t_ZdcModuleBkgdMaxFraction, "zdc_ZdcModuleBkgdMaxFraction[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleAmpError", &t_ZdcModuleAmpError, "zdc_ZdcModuleAmpError[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleMinDeriv2nd", &t_ZdcModuleMinDeriv2nd, "zdc_ZdcModuleMinDeriv2nd[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModulePresample", &t_ZdcModulePresample, "zdc_ZdcModulePresample[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModulePreSampleAmp", &t_ZdcModulePreSampleAmp, "zdc_ZdcModulePreSampleAmp[2][4]/F");
    m_outputTree->Branch("zdc_ZdcLucrodTriggerAmp",&t_ZdcLucrodTriggerAmp,"zdc_ZdcLucrodTriggerAmp[2][4]/S");
    m_outputTree->Branch("zdc_ZdcLucrodTriggerAmpLG",&t_ZdcLucrodTriggerAmpLG,"zdc_ZdcLucrodTriggerAmpLG[2][4]/S");
    m_outputTree->Branch("zdc_ZdcModuleMaxADC",&t_ZdcModuleMaxADC,"zdc_ZdcModuleMaxADC[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleMaxADCHG",&t_ZdcModuleMaxADCHG,"zdc_ZdcModuleMaxADCHG[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleMaxADCLG",&t_ZdcModuleMaxADCLG,"zdc_ZdcModuleMaxADCLG[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModulePeakADCHG",&t_ZdcModulePeakADCHG,"zdc_ZdcModulePeakADCHG[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModulePeakADCLG",&t_ZdcModulePeakADCLG,"zdc_ZdcModulePeakADCLG[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleAmpLGRefit", &t_ZdcModuleAmpLGRefit, "zdc_ZdcModuleAmpLGRefit[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleAmpCorrLGRefit", &t_ZdcModuleAmpCorrLGRefit, "zdc_ZdcModuleAmpCorrLGRefit[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleT0LGRefit", &t_ZdcModuleT0LGRefit, "zdc_ZdcModuleT0LGRefit[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleT0SubLGRefit", &t_ZdcModuleT0SubLGRefit, "zdc_ZdcModuleT0SubLGRefit[2][4]/F");
    m_outputTree->Branch("zdc_ZdcModuleChisqLGRefit", &t_ZdcModuleChisqLGRefit, "zdc_ZdcModuleChisqLGRefit[2][4]/F");
    
    if(m_isMC){
      //Modules
      m_outputTree->Branch("zdc_ZdcModuleTruthTotal",&t_ZdcModuleTruthTotal,"zdc_ZdcModuleTruthTotal[2][7]/F");
      m_outputTree->Branch("zdc_ZdcModuleTruthInvisible",&t_ZdcModuleTruthInvis,"zdc_ZdcModuleTruthInvisible[2][7]/F");
      m_outputTree->Branch("zdc_ZdcModuleTruthEM",&t_ZdcModuleTruthEM,"zdc_ZdcModuleTruthEM[2][7]/F");
      m_outputTree->Branch("zdc_ZdcModuleTruthNonEM",&t_ZdcModuleTruthNonEM,"zdc_ZdcModuleTruthNonEM[2][7]/F");
      m_outputTree->Branch("zdc_ZdcModuleTruthEscaped",&t_ZdcModuleTruthEscaped,"zdc_ZdcModuleTruthEscaped[2][7]/F");
      m_outputTree->Branch("zdc_ZdcModuleTruthNphotons",&t_ZdcModuleTruthNphotons,"zdc_ZdcModuleTruthNphotons[2][7]/i");
      m_outputTree->Branch("zdc_RpdModuleTruthNphotons",&t_RpdModuleTruthNphotons,"zdc_RpdModuleTruthNphotons[2][16]/i");
      
      //Sums
      m_outputTree->Branch("zdc_ZdcTruthTotal",&t_ZdcTruthTotal,"zdc_ZdcTruthTotal[2]/F");
      m_outputTree->Branch("zdc_ZdcTruthInvisible",&t_ZdcTruthInvis, "zdc_ZdcTruthInvisible[2]/F");
      m_outputTree->Branch("zdc_ZdcTruthEM",&t_ZdcTruthEM,"zdc_ZdcTruthEM[2]/F");
      m_outputTree->Branch("zdc_ZdcTruthNonEM",&t_ZdcTruthNonEM,"zdc_ZdcTruthNonEM[2]/F");
      m_outputTree->Branch("zdc_ZdcTruthEscaped",&t_ZdcTruthEscaped,"zdc_ZdcTruthEscaped[2]/F");
      
      //Event gen particles
      m_outputTree->Branch("zdc_ZdcTruthParticlePosx",&t_ZdcTruthParticlePosx);
      m_outputTree->Branch("zdc_ZdcTruthParticlePosy",&t_ZdcTruthParticlePosy);
      m_outputTree->Branch("zdc_ZdcTruthParticlePosz",&t_ZdcTruthParticlePosz);
      m_outputTree->Branch("zdc_ZdcTruthParticleTime",&t_ZdcTruthParticleTime);
      m_outputTree->Branch("zdc_ZdcTruthParticlePx",&t_ZdcTruthParticlePx);
      m_outputTree->Branch("zdc_ZdcTruthParticlePy",&t_ZdcTruthParticlePy);
      m_outputTree->Branch("zdc_ZdcTruthParticlePz",&t_ZdcTruthParticlePz);
      m_outputTree->Branch("zdc_ZdcTruthParticleEnergy",&t_ZdcTruthParticleEnergy);
      m_outputTree->Branch("zdc_ZdcTruthParticlePid",&t_ZdcTruthParticlePid);
      m_outputTree->Branch("zdc_ZdcTruthParticleStatus",&t_ZdcTruthParticleStatus);
    }
      }
    if (enableRPD)
      {
      if (enableRPDAmp) {
    m_outputTree->Branch("zdc_RpdChannelBaseline",&t_RpdChannelBaseline,"zdc_RpdChannelBaseline[2][16]/F");
    m_outputTree->Branch("zdc_RpdChannelPileupExpFitParams",&t_RpdChannelPileupExpFitParams,"zdc_RpdChannelPileupExpFitParams[2][16][2]/F");
    m_outputTree->Branch("zdc_RpdChannelPileupStretchedExpFitParams",&t_RpdChannelPileupStretchedExpFitParams,"zdc_RpdChannelPileupStretchedExpFitParams[2][16][3]/F");
    m_outputTree->Branch("zdc_RpdChannelPileupExpFitParamErrs",&t_RpdChannelPileupExpFitParamErrs,"zdc_RpdChannelPileupExpFitParamErrs[2][16][2]/F");
    m_outputTree->Branch("zdc_RpdChannelPileupStretchedExpFitParamErrs",&t_RpdChannelPileupStretchedExpFitParamErrs,"zdc_RpdChannelPileupStretchedExpFitParamErrs[2][16][3]/F");
    m_outputTree->Branch("zdc_RpdChannelPileupExpFitMSE",&t_RpdChannelPileupExpFitMSE,"zdc_RpdChannelPileupExpFitMSE[2][16]/F");
    m_outputTree->Branch("zdc_RpdChannelPileupStretchedExpFitMSE",&t_RpdChannelPileupStretchedExpFitMSE,"zdc_RpdChannelPileupStretchedExpFitMSE[2][16]/F");
    m_outputTree->Branch("zdc_RpdChannelAmplitude",&t_RpdChannelAmplitude,"zdc_RpdChannelAmplitude[2][16]/F");
    m_outputTree->Branch("zdc_RpdChannelAmplitudeCalib",&t_RpdChannelAmplitudeCalib,"zdc_RpdChannelAmplitudeCalib[2][16]/F");
    m_outputTree->Branch("zdc_RpdChannelMaxADC",&t_RpdChannelMaxADC,"zdc_RpdChannelMaxADC[2][16]/F");
    m_outputTree->Branch("zdc_RpdChannelMaxADCCalib",&t_RpdChannelMaxADCCalib,"zdc_RpdChannelMaxADCCalib[2][16]/F");
    m_outputTree->Branch("zdc_RpdChannelMaxSample",&t_RpdChannelMaxSample,"zdc_RpdChannelMaxSample[2][16]/i");
    m_outputTree->Branch("zdc_RpdChannelStatus",&t_RpdChannelStatus,"zdc_RpdChannelStatus[2][16]/i");
    m_outputTree->Branch("zdc_RpdChannelPileupFrac",&t_RpdChannelPileupFrac,"zdc_RpdChannelPileupFrac[2][16]/F");
    m_outputTree->Branch("zdc_RpdSideStatus",&t_RpdSideStatus,"zdc_RpdSideStatus[2]/i");
      }
    if (enableCentroid)
      {
    m_outputTree->Branch("zdc_centroidAvailable", &t_centroidDecorationsAvailable);
    m_outputTree->Branch("zdc_centroidEventValid", &t_centroidEventValid, "zdc_centroidEventValid/B");
    m_outputTree->Branch("zdc_centroidStatus", &t_centroidStatus, "zdc_centroidStatus[2]/i");
    m_outputTree->Branch("zdc_RPDChannelSubtrAmp", &t_RPDChannelSubtrAmp, "zdc_RPDChannelSubtrAmp[2][16]/F");
    m_outputTree->Branch("zdc_RPDSubtrAmpSum", &t_RPDSubtrAmpSum, "zdc_RPDSubtrAmpSum[2]/F");
    m_outputTree->Branch("zdc_xCentroidPreGeomCorPreAvgSubtr", &t_xCentroidPreGeomCorPreAvgSubtr, "zdc_xCentroidPreGeomCorPreAvgSubtr[2]/F");
    m_outputTree->Branch("zdc_yCentroidPreGeomCorPreAvgSubtr", &t_yCentroidPreGeomCorPreAvgSubtr, "zdc_yCentroidPreGeomCorPreAvgSubtr[2]/F");
    m_outputTree->Branch("zdc_xCentroidPreAvgSubtr", &t_xCentroidPreAvgSubtr, "zdc_xCentroidPreAvgSubtr[2]/F");
    m_outputTree->Branch("zdc_yCentroidPreAvgSubtr", &t_yCentroidPreAvgSubtr, "zdc_yCentroidPreAvgSubtr[2]/F");
    m_outputTree->Branch("zdc_xCentroid", &t_xCentroid, "zdc_xCentroid[2]/F");
    m_outputTree->Branch("zdc_yCentroid", &t_yCentroid, "zdc_yCentroid[2]/F");
    m_outputTree->Branch("zdc_xRowCentroid", &t_xRowCentroid, "zdc_xRowCentroid[2][4]/F");
    m_outputTree->Branch("zdc_yColCentroid", &t_yColCentroid, "zdc_yColCentroid[2][4]/F");
    m_outputTree->Branch("zdc_reactionPlaneAngle", &t_reactionPlaneAngle, "zdc_reactionPlaneAngle[2]/F");
    m_outputTree->Branch("zdc_cosDeltaReactionPlaneAngle", &t_cosDeltaReactionPlaneAngle, "zdc_cosDeltaReactionPlaneAngle/F");
      }
    }
    
    if (!(zdcCalib || zdcLaser || zdcOnly || zdcInj))
    {
      m_outputTree->Branch("mbts_in_e", &t_mbts_in_e, "mbts_in_e[2][8]/F");
      m_outputTree->Branch("mbts_in_t", &t_mbts_in_t, "mbts_in_t[2][8]/F");
      m_outputTree->Branch("mbts_out_e", &t_mbts_out_e, "mbts_out_e[2][4]/F");
      m_outputTree->Branch("mbts_out_t", &t_mbts_out_t, "mbts_out_t[2][4]/F");
 
      m_outputTree->Branch("T2mbts_in_e", &t_T2mbts_in_e, "T2mbts_in_e[2][8]/F");
      m_outputTree->Branch("T2mbts_in_t", &t_T2mbts_in_t, "T2mbts_in_t[2][8]/F");
      m_outputTree->Branch("T2mbts_out_e", &t_T2mbts_out_e, "T2mbts_out_e[2][4]/F");
      m_outputTree->Branch("T2mbts_out_t", &t_T2mbts_out_t, "T2mbts_out_t[2][4]/F");
 
      m_outputTree->Branch("L1ET", &t_L1ET, "L1ET/F");
      m_outputTree->Branch("L1ET24", &t_L1ET24, "L1ET24/F");
 
      m_outputTree->Branch("totalEt", &t_totalEt, "totalEt/F");
      m_outputTree->Branch("totalEt_TTsum", &t_totalEt_TTsum, "totalEt_TTsum/F");
 
      m_outputTree->Branch("totalEt24", &t_totalEt24, "totalEt24/F");
      m_outputTree->Branch("totalEt24_TTsum", &t_totalEt24_TTsum, "totalEt24_TTsum/F");
 
      m_outputTree->Branch("fcalEt", &t_fcalEt, "fcalEt/F");
      m_outputTree->Branch("fcalEtA", &t_fcalEtA, "fcalEtA/F");
      m_outputTree->Branch("fcalEtC", &t_fcalEtC, "fcalEtC/F");
      m_outputTree->Branch("fcalEtA_TT", &t_fcalEtA_TT, "fcalEtA_TT/F");
      m_outputTree->Branch("fcalEtC_TT", &t_fcalEtC_TT, "fcalEtC_TT/F");
 
      m_outputTree->Branch("nvx", &t_nvx, "nvx/I");
      m_outputTree->Branch("vx", &t_vx, "vx[3]/F");
      m_outputTree->Branch("pvindex", &t_pvindex, "pvindex/I");
      m_outputTree->Branch("vxntrk", &t_vxntrk, "vxntrk/I");
      m_outputTree->Branch("vx_trk_index", "vector<int>", &t_vx_trk_index);
      m_outputTree->Branch("vxtype", &t_vxtype, "vxtype/I");
      m_outputTree->Branch("vxcov", &t_vxcov, "vxcov[6]/F");
      m_outputTree->Branch("vxsumpt2", &t_vxsumpt2, "vxsumpt2/F");
      m_outputTree->Branch("nstrong", &t_nstrong, "nstrong/I");
      m_outputTree->Branch("puvxntrk", &t_puvxntrk, "puvxntrk/I");
      m_outputTree->Branch("puvxsumpt", &t_puvxsumpt, "puvxsumpt/F");
      m_outputTree->Branch("puvxz", &t_puvxz, "puvxz/F");
      m_outputTree->Branch("vxnlooseprimary", &t_vxnlooseprimary, "vxnlooseprimary/I");
      m_outputTree->Branch("vxnminbias", &t_vxnminbias, "vxnminbias/I");
      m_outputTree->Branch("vxngoodmuon", &t_vxngoodmuon, "vxngoodmuon/I");
 
      m_outputTree->Branch("t_nvtx", &t_nvtx, "nvtx/I");
      m_outputTree->Branch("vtx_type", "vector<int8_t>", &t_vtx_type);
      m_outputTree->Branch("vtx_x", "vector<float>", &t_vtx_x);
      m_outputTree->Branch("vtx_y", "vector<float>", &t_vtx_y);
      m_outputTree->Branch("vtx_z", "vector<float>", &t_vtx_z);
      m_outputTree->Branch("vtx_ntrk_all", "vector<int16_t>", &t_vtx_ntrk_all);
      m_outputTree->Branch("vtx_sumpt2_all", "vector<float>", &t_vtx_sumpt2_all);
      m_outputTree->Branch("vtx_ntrk", "vector<int16_t>", &t_vtx_ntrk);
      m_outputTree->Branch("vtx_sumpt2", "vector<float>", &t_vtx_sumpt2);
      m_outputTree->Branch("vtx_trk_index", "vector< vector<int16_t> >", &t_vtx_trk_index);
 
      m_outputTree->Branch("mbts_countA", &t_mbts_countA, "mbts_countA/s");
      m_outputTree->Branch("mbts_countC", &t_mbts_countC, "mbts_countC/s");
      m_outputTree->Branch("T2mbts_countAin", &t_T2mbts_countAin, "T2mbts_countAin/s");
      m_outputTree->Branch("T2mbts_countCin", &t_T2mbts_countCin, "T2mbts_countCin/s");
      m_outputTree->Branch("mbts_timeA", &t_mbts_timeA, "mbts_timeA/F");
      m_outputTree->Branch("mbts_timeC", &t_mbts_timeC, "mbts_timeC/F");
      m_outputTree->Branch("mbts_timeDiff", &t_mbts_timeDiff, "mbts_timeDiff/F");
 
      t_nclus = 0;
      m_outputTree->Branch("nclus", &t_nclus, "nclus/i");
      m_outputTree->Branch("clusEt", &t_clusEt, "clusEt/F");
      m_outputTree->Branch("clusEtMax", &t_clusEtMax, "clusEtMax/F");
      m_outputTree->Branch("clusetaMax", &t_clusetaMax, "clusetaMax/F");
      m_outputTree->Branch("clusphiMax", &t_clusphiMax, "clusphiMax/F");
 
      m_outputTree->Branch("cc_pt", "vector<float>", &t_cc_pt);
      m_outputTree->Branch("cc_eta", "vector<float>", &t_cc_eta);
      m_outputTree->Branch("cc_phi", "vector<float>", &t_cc_phi);
      m_outputTree->Branch("cc_e", "vector<float>", &t_cc_e);
      m_outputTree->Branch("cc_raw_m", "vector<float>", &t_cc_raw_m);
      m_outputTree->Branch("cc_raw_eta", "vector<float>", &t_cc_raw_eta);
      m_outputTree->Branch("cc_raw_phi", "vector<float>", &t_cc_raw_phi);
      m_outputTree->Branch("cc_raw_e", "vector<float>", &t_cc_raw_e);
      m_outputTree->Branch("cc_raw_samp", "vector<vector<float>>", &t_cc_raw_samp);
      m_outputTree->Branch("cc_sig", "vector<float>", &t_cc_sig);
      m_outputTree->Branch("cc_layer", "vector<int>", &t_cc_layer);
 
      m_outputTree->Branch("edgeGapA", &t_edgeGapA, "edgeGapA/F");
      m_outputTree->Branch("edgeGapC", &t_edgeGapC, "edgeGapC/F");
 
      m_outputTree->Branch("ntrk", &t_ntrk, "ntrk/i");
      if (enableTracks)
        {
          m_outputTree->Branch("trk_pt", "vector<float>", &t_trk_pt);
          m_outputTree->Branch("trk_eta", "vector<float>", &t_trk_eta);
          m_outputTree->Branch("trk_theta", "vector<float>", &t_trk_theta);
          m_outputTree->Branch("trk_phi", "vector<float>", &t_trk_phi);
          m_outputTree->Branch("trk_e", "vector<float>", &t_trk_e);
          m_outputTree->Branch("trk_index", "vector<int>", &t_trk_index);
          m_outputTree->Branch("trk_d0", "vector<float>", &t_trk_d0);
          m_outputTree->Branch("trk_z0", "vector<float>", &t_trk_z0);
          m_outputTree->Branch("trk_vz", "vector<float>", &t_trk_vz);
          m_outputTree->Branch("trk_vtxz", "vector<float>", &t_trk_vtxz);
          m_outputTree->Branch("trk_pixeldEdx", "vector<float>", &t_trk_pixeldEdx);
          m_outputTree->Branch("trk_charge", "vector<int8_t>", &t_trk_charge);
          m_outputTree->Branch("trk_quality", "vector<int16_t>", &t_trk_quality);
          m_outputTree->Branch("trk_nPixHits", "vector<uint8_t>", &t_trk_nPixHits);
          m_outputTree->Branch("trk_nSctHits", "vector<uint8_t>", &t_trk_nSctHits);
          m_outputTree->Branch("trk_nPixDead", "vector<uint8_t>", &t_trk_nPixDead);
          m_outputTree->Branch("trk_nSctDead", "vector<uint8_t>", &t_trk_nSctDead);
          m_outputTree->Branch("trk_nPixHoles", "vector<uint8_t>", &t_trk_nPixHoles);
          m_outputTree->Branch("trk_nSctHoles", "vector<uint8_t>", &t_trk_nSctHoles);
          m_outputTree->Branch("trk_nTrtHits", "vector<uint8_t>", &t_trk_nTrtHits);
          m_outputTree->Branch("trk_nTrtOutliers", "vector<uint8_t>", &t_trk_nTrtOutliers);
          m_outputTree->Branch("trk_inPixHits", "vector<uint8_t>", &t_trk_inPixHits);
          m_outputTree->Branch("trk_exPixHits", "vector<uint8_t>", &t_trk_exPixHits);
          m_outputTree->Branch("trk_ninPixHits", "vector<uint8_t>", &t_trk_ninPixHits);
          m_outputTree->Branch("trk_nexPixHits", "vector<uint8_t>", &t_trk_nexPixHits);
        }
 
      if( lhcf2022 || lhcf2022zdc || lhcf2022afp){
    m_outputTree->Branch("nProtons", &nProtons, "nProtons/i");
    m_outputTree->Branch("proton_pt", "vector<double>", &proton_pt);
    m_outputTree->Branch("proton_eta", "vector<double>", &proton_eta);
    m_outputTree->Branch("proton_phi", "vector<double>", &proton_phi);
    m_outputTree->Branch("proton_e", "vector<double>", &proton_e);
    m_outputTree->Branch("proton_side", "vector<int>", &proton_side);
    m_outputTree->Branch("proton_eLoss", "vector<double>", &proton_eLoss);
    m_outputTree->Branch("proton_t", "vector<double>", &proton_t);
    m_outputTree->Branch("proton_track_stationID", "vector<vector<int>>", &proton_track_stationID);
    m_outputTree->Branch("proton_track_nClusters", "vector<vector<int>>", &proton_track_nClusters);
    m_outputTree->Branch("proton_track_xLocal", "vector<vector<float>>", &proton_track_xLocal);
    m_outputTree->Branch("proton_track_yLocal", "vector<vector<float>>", &proton_track_yLocal);
    m_outputTree->Branch("proton_track_zLocal", "vector<vector<float>>", &proton_track_zLocal);
    m_outputTree->Branch("proton_track_xSlope", "vector<vector<float>>", &proton_track_xSlope);
    m_outputTree->Branch("proton_track_ySlope", "vector<vector<float>>", &proton_track_ySlope);
      }
 
    }
  }
  
  ANA_MSG_DEBUG("Anti-howdy from Initialize!");
  
  
  // Here you do everything that you need to do after the first input
  // file has been connected and before the first event is processed,
  // e.g. create additional histograms based on which variables are
  // available in the input files.  You can also create all of your
  // histograms and trees in here, but be aware that this method
  // doesn't get called if no events are processed.  So any objects
  // you create here won't be available in the output if you have no
  // input events.
 
  ANA_MSG_INFO("enableOutputTree = " << enableOutputTree);
  ANA_MSG_INFO("writeOnlyTriggers = " << writeOnlyTriggers);
  ANA_MSG_INFO("enableOutputSamples = " << enableOutputSamples);
  ANA_MSG_INFO("enableTrigger = " << enableTrigger);
  ANA_MSG_INFO("enableRPD = " << enableRPD);
  ANA_MSG_INFO("enableCentroid = " << enableCentroid);
  ANA_MSG_INFO("zdcCalib = " << zdcCalib);
  ANA_MSG_INFO("zdcInj = " << zdcInj);
  ANA_MSG_INFO("zdcLaser = " << zdcLaser);
  ANA_MSG_INFO("zdcConfig = " << zdcConfig);
  ANA_MSG_INFO("reprocZdc = " << reprocZdc);
  ANA_MSG_INFO("auxSuffix = " << auxSuffix );
  ANA_MSG_INFO("zdcLowGainMode = " << zdcLowGainMode);
  ANA_MSG_INFO("enableID = " << enableID);
  ANA_MSG_INFO("enableCalo = " << enableCalo);
  ANA_MSG_INFO("enableClusters = " << enableClusters);
  ANA_MSG_INFO("trackLimit = " << trackLimit);
  ANA_MSG_INFO("trackLimitReject = " << trackLimitReject);
  ANA_MSG_INFO("isMC = " << m_isMC);
  ANA_MSG_INFO("useGRL = " <<  useGRL);
  ANA_MSG_INFO("grlFilename = " <<  grlFilename);
  ANA_MSG_INFO("slimmed = " <<  slimmed);
  ANA_MSG_INFO("zdcOnly = " << zdcOnly);
  ANA_MSG_INFO("flipDelay = " <<  flipDelay);
  ANA_MSG_INFO("nsamplesZdc = " <<  nsamplesZdc);
  ANA_MSG_INFO("lhcf2022 = " << lhcf2022);
  ANA_MSG_INFO("lhcf2022afp = " << lhcf2022afp);
  ANA_MSG_INFO("lhcf2022zdc = " << lhcf2022zdc);
  ANA_MSG_INFO("doZdcCalib = " << doZdcCalib);
 
  ANA_MSG_DEBUG("initialize: Initialize!");
 
 
  // GRL
 
  if (useGRL)
  {
    const char* fullGRLFilePath = gSystem->ExpandPathName (grlFilename.c_str());
    ANA_MSG_INFO("GRL: " << fullGRLFilePath);
    std::vector<std::string> vecStringGRL;
    vecStringGRL.push_back(fullGRLFilePath);
    ANA_CHECK(m_grl.setProperty( "GoodRunsListVec", vecStringGRL));
    ANA_CHECK(m_grl.setProperty("PassThrough", false)); // if true (default) will ignore result of GRL and will just pass all events
    ANA_CHECK(m_grl.initialize());
 
  }
 
  if (enableTracks)
    {
      ANA_CHECK(m_selTool.retrieve());
    }
 
  if (enableTrigger) // HLT related
  {
    ANA_MSG_INFO("Trying to initialize TDT");
    ANA_CHECK(m_trigDecisionTool.retrieve());
  }
 
  // ZDC re-reco tool
  if (reprocZdc)
  {
    m_zdcAnalysisTool.setTypeAndName("ZDC::ZdcAnalysisTool/ZdcAnalysisTool");
 
    ANA_MSG_INFO("Trying to configure ZDC Analysis Tool!");
 
    ANA_CHECK(m_zdcAnalysisTool.setProperty("FlipEMDelay", flipDelay));
    ANA_CHECK(m_zdcAnalysisTool.setProperty("LowGainMode", zdcLowGainMode));
    ANA_CHECK(m_zdcAnalysisTool.setProperty("DoCalib", doZdcCalib));
    ANA_CHECK(m_zdcAnalysisTool.setProperty("Configuration", zdcConfig));
    ANA_CHECK(m_zdcAnalysisTool.setProperty("AuxSuffix", auxSuffix));
    ANA_CHECK(m_zdcAnalysisTool.setProperty("ForceCalibRun", -1));
    
    ANA_MSG_INFO("Setting up zdcConfig=" << zdcConfig);
    if (zdcConfig=="LHCf2022")
      {
    ANA_CHECK(m_zdcAnalysisTool.setProperty("DoTrigEff", false)); // for now
    ANA_CHECK(m_zdcAnalysisTool.setProperty("DoTimeCalib", false)); // for now
    ANA_CHECK(m_zdcAnalysisTool.setProperty("Configuration", "LHCf2022"));
      }
    else if (zdcConfig == "PbPb2018")
      {
    ANA_CHECK(m_zdcAnalysisTool.setProperty("Configuration", "PbPb2018"));    
      }
    else if (zdcConfig == "pPb2016")
      {
    ANA_CHECK(m_zdcAnalysisTool.setProperty("Configuration", "pPb2016"));
      }
    else if (zdcConfig == "PbPb2015")
      {
    ANA_CHECK(m_zdcAnalysisTool.setProperty("Configuration", "PbPb2015"));
      }
    
    if (flipDelay)
      ANA_MSG_INFO("FLIP ZDC DELAY IN EM MODULES");
    else
      ANA_MSG_INFO("NO FLIP ZDC DELAY IN EM MODULES");
 
 
    ANA_MSG_INFO("Trying to initialize ZDC Analysis Tool!");
    ANA_CHECK(m_zdcAnalysisTool.initialize());
  }
  
  if (zdcInj)
    {
      m_zdcInjPulserAmpMap = std::make_shared<ZdcInjPulserAmpMap>();
      ATH_MSG_INFO( "Using JSON file for injector-pulse voltage at path " << m_zdcInjPulserAmpMap->getFilePath() );
    }
  
  return StatusCode::SUCCESS;
}

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< 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.

inputMetaStore() [1/2]

AnaAlgorithm::MetaStorePtr_t EL::AnaAlgorithm::inputMetaStore ( )

inherited

Definition at line 81 of file AnaAlgorithm.cxx.

  {
#ifdef XAOD_STANDALONE
     return &m_inputMetaStore;
#else
     return m_inputMetaStore;
#endif // XAOD_STANDALONE
  }

inputMetaStore() [2/2]

AnaAlgorithm::ConstMetaStorePtr_t EL::AnaAlgorithm::inputMetaStore ( ) const

inherited

Accessor for the input metadata store

Definition at line 72 of file AnaAlgorithm.cxx.

  {
#ifdef XAOD_STANDALONE
     return &m_inputMetaStore;
#else
     return m_inputMetaStore;
#endif // XAOD_STANDALONE
  }

msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

myReplace()

void AthHistogramming::myReplace ( std::string & str, const std::string & oldStr, const std::string & newStr )

privateinherited

Helper method to replace sub-string.

Definition at line 590 of file AthHistogramming.cxx.

{
  size_t pos = 0;
  while((pos = str.find(oldStr, pos)) != std::string::npos)
    {
      str.replace(pos, oldStr.length(), newStr);
      pos += newStr.length();
    }
}

outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< 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.

outputMetaStore() [1/2]

AnaAlgorithm::MetaStorePtr_t EL::AnaAlgorithm::outputMetaStore ( )

inherited

Definition at line 101 of file AnaAlgorithm.cxx.

  {
#ifdef XAOD_STANDALONE
     return &m_outputMetaStore;
#else
     return m_outputMetaStore;
#endif // XAOD_STANDALONE
  }

outputMetaStore() [2/2]

AnaAlgorithm::ConstMetaStorePtr_t EL::AnaAlgorithm::outputMetaStore ( ) const

inherited

Accessor for the output metadata store

Definition at line 92 of file AnaAlgorithm.cxx.

  {
#ifdef XAOD_STANDALONE
     return &m_outputMetaStore;
#else
     return m_outputMetaStore;
#endif // XAOD_STANDALONE
  }

print()

void EL::AnaAlgorithm::print ( ) const

protectedvirtualinherited

print the state of the algorithm

This is mostly to allow algorithms to add a little debugging information if they feel like it.

Definition at line 336 of file AnaAlgorithm.cxx.

  {}

processClusters()

void ZdcNtuple::processClusters

(

)

Definition at line 1843 of file ZdcNtuple.cxx.

{
  //t_nclus = 0;
 
  t_cc_pt.clear();
  t_cc_eta.clear();
  t_cc_phi.clear();
  t_cc_e.clear();
  t_cc_raw_m.clear();
  t_cc_raw_eta.clear();
  t_cc_raw_phi.clear();
  t_cc_raw_e.clear();
  t_cc_raw_samp.clear();
  t_cc_layer.clear();
  t_cc_sig.clear();
 
  t_nclus = m_caloClusters->size();
 
  t_clusEt = 0;
  t_clusEtMax = -999;
  t_clusetaMax = 0;
  t_clusphiMax = 0;
 
  for (const auto cluster : *m_caloClusters)
  {
    t_cc_pt.push_back(cluster->pt());
    t_cc_eta.push_back(cluster->eta());
    t_cc_phi.push_back(cluster->phi());
    t_cc_e.push_back(cluster->e());
    t_cc_raw_m.push_back(cluster->rawM());
    t_cc_raw_eta.push_back(cluster->rawEta());
    t_cc_raw_phi.push_back(cluster->rawPhi());
    t_cc_raw_e.push_back(cluster->rawE());
 
    std::vector<float> energies;
 
    for (size_t s = CaloSampling::PreSamplerB; s < CaloSampling::Unknown; s++ )
    {
      bool hasSample = cluster->hasSampling( (xAOD::CaloCluster::CaloSample) s );
      float e = 0;
      if (hasSample)
      {
        e = cluster->eSample( (xAOD::CaloCluster::CaloSample) s);
      }
      energies.push_back(e);
    }
    t_cc_raw_samp.push_back(energies);
 
    float et = cluster->e() / TMath::CosH(cluster->eta());
    t_clusEt += et;
    if (et > t_clusEtMax)
    {
      t_clusEtMax = et;
      t_clusetaMax = cluster->eta();
      t_clusphiMax = cluster->phi();
    }
 
    double cell_sig = 0;
    if (!cluster->retrieveMoment(xAOD::CaloCluster::CELL_SIGNIFICANCE, cell_sig)) {ANA_MSG_DEBUG("processClusters() : No CELL_SIGNIFICANCE!");}
    t_cc_sig.push_back(cell_sig);
    double cell_layer = 0;
    if (!cluster->retrieveMoment(xAOD::CaloCluster::CELL_SIG_SAMPLING, cell_layer)) {ANA_MSG_DEBUG("processClusters() : No CELL_SIG_SAMPLING!");}
    t_cc_layer.push_back(static_cast<int>(cell_layer));
    //t_nclus++;
  }
 
  if ( (!enableClusters) || (t_ntrk >= trackLimit) ) // if disabled or if too many tracks
  {
    t_cc_pt.clear();
    t_cc_eta.clear();
    t_cc_phi.clear();
    t_cc_e.clear();
    t_cc_layer.clear();
    t_cc_sig.clear();
  }
  else
  {
    ANA_MSG_DEBUG("processClusters(): keeping clusters");
  }
  return;
}

processEventInfo()

void ZdcNtuple::processEventInfo

(

)

Definition at line 1279 of file ZdcNtuple.cxx.

{
  ANA_MSG_DEBUG( "processing event info");
 
  t_bcid = m_eventInfo->bcid();
  t_runNumber = m_eventInfo->runNumber();
  t_eventNumber = m_eventInfo->eventNumber();
  t_lumiBlock = m_eventInfo->lumiBlock();
  t_bunchGroup = -1;
  t_extendedLevel1ID = m_eventInfo->extendedLevel1ID();
  t_timeStamp = m_eventInfo->timeStamp();
  t_timeStampNSOffset = m_eventInfo->timeStampNSOffset();
  t_passBits = acceptEvent();
  t_avgIntPerCrossing = m_eventInfo->averageInteractionsPerCrossing();
  t_actIntPerCrossing = m_eventInfo->actualInteractionsPerCrossing();
  t_zdcEventInfoError = m_eventInfo->errorState(xAOD::EventInfo::ForwardDet);
  t_zdcEventInfoErrorWord = m_eventInfo->eventFlags(xAOD::EventInfo::ForwardDet);
  
  if ( !(m_eventCounter++ % 1000) || msgLvl(MSG::DEBUG))
  {
    ANA_MSG_INFO("Event# " << m_eventCounter << " Run " << m_eventInfo->runNumber() << " Event " << m_eventInfo->eventNumber() << " LB " << m_eventInfo->lumiBlock() );
  }
 
}

processFCal()

void ZdcNtuple::processFCal

(

)

Definition at line 1617 of file ZdcNtuple.cxx.

{
  ANA_MSG_DEBUG("processFCal: processing FCal");
 
  t_fcalEt = 0.;
  t_fcalEtA = 0.;
  t_fcalEtC = 0.;
  t_fcalEtA_TT = 0.;
  t_fcalEtC_TT = 0.;
 
  if (m_caloSums)
    {
      static const SG::ConstAccessor<std::string> SummaryAcc("Summary");
      for (const auto calosum : *m_caloSums)
    {
      const std::string name = SummaryAcc(*calosum);
      if (name == "FCal")
        {
          t_fcalEt = calosum->et();
          ANA_MSG_DEBUG("processFCal: fcalEt = " << t_fcalEt);
        }
      
      if (name == "All")
        {
          t_totalEt = calosum->et();
          ANA_MSG_DEBUG("processFCal: totalEt = " << t_totalEt);
        }
    }
    }
 
  t_fcalEtA = 0;
  t_fcalEtC = 0;
  t_totalEt24 = 0;
 
  if (m_eventShapes)
    {
      for (const auto eventShape : *m_eventShapes)
    {
      int layer = eventShape->layer();
      float eta = eventShape->etaMin();
      float et = eventShape->et();
      if (layer == 21 || layer == 22 || layer == 23)
        {
          if (eta > 0) t_fcalEtA += et;
          if (eta < 0) t_fcalEtC += et;
        }
      
      if (TMath::Abs(eta) < 2.4)
        {
          t_totalEt24 += et;
        }
    }
    }
 
  t_L1ET = 0;
  t_L1ET24 = 0;
 
  if (m_lvl1EnergySumRoI)
  {
    t_L1ET = m_lvl1EnergySumRoI->energyT();
    //t_L1ET24 = m_lvl1EnergySumRoI->energyTRestricted(); // TBD when limited eta ET available
  }
 
  return;
}

processGaps()

void ZdcNtuple::processGaps

(

)

Definition at line 1683 of file ZdcNtuple.cxx.

{
 
  float eta_min = 5;
  float eta_max = -5;
 
  if (!m_caloClusters) return;
  for (const auto cl : *m_caloClusters)
  {
 
    if (cl->pt() < m_gapPtMin) continue;
 
    int etabin = h_TCSigCut->GetXaxis()->FindBin(cl->eta());
    if (etabin < 1 || etabin > h_TCSigCut->GetNbinsX()) continue;
    float sig_cut = h_TCSigCut->GetBinContent(etabin);
    float sig = cl->getMomentValue(xAOD::CaloCluster::CELL_SIGNIFICANCE);
    int cl_cell_sig_samp = static_cast<int>(cl->getMomentValue(xAOD::CaloCluster::CELL_SIG_SAMPLING));
 
    ANA_MSG_VERBOSE ("gapclus: etabin " << etabin << " sig_cut=" << sig_cut << " sig=" << sig << " samp=" << cl_cell_sig_samp);
 
    if (sig < sig_cut) continue;
 
    if (cl_cell_sig_samp >= CaloSampling::TileBar0 && cl_cell_sig_samp <= CaloSampling::TileExt2) continue;
 
    if (cl->eta() < eta_min) eta_min = cl->eta();
    if (cl->eta() > eta_max) eta_max = cl->eta();
 
  }
 
  t_edgeGapA = 4.9 - eta_max;
  t_edgeGapC = eta_min + 4.9;
  ANA_MSG_DEBUG("processGaps(): egA " << t_edgeGapA << " , egC " << t_edgeGapC);
 
}

processInDet()

void ZdcNtuple::processInDet

(

)

Definition at line 1329 of file ZdcNtuple.cxx.

{
  ANA_MSG_DEBUG("processInDet(): processing tracks & vertices!");
  t_ntrk = 0;
  t_nvx = 0;
  t_vxntrk = 0;
  t_vx_trk_index.clear();
  t_vxsumpt2 = 0;
  t_vxtype = 0;
  t_pvindex = -1;
  t_puvxntrk = 0;
  t_puvxsumpt = 0;
  t_vxnlooseprimary = 0;
  t_vxnminbias = 0;
 
  int i;
  for (i = 0; i < 3; i++) t_vx[i] = 0;
  for (i = 0; i < 6; i++) t_vxcov[i] = 0;
 
  const xAOD::Vertex* primary_vertex = nullptr;
  size_t pv_index = -1;
  size_t vx_index = 0;
  float max_pileup_sumpT = 0.;
  int max_pileup_nTrack = 0;
  float max_pileup_z = 0;
  int nStrongPileup = 0;
 
  t_nvtx = 0;
  t_vtx_type.clear();
  t_vtx_x.clear();
  t_vtx_y.clear();
  t_vtx_z.clear();
  t_vtx_ntrk_all.clear();
  t_vtx_sumpt2_all.clear();
  t_vtx_ntrk.clear();
  t_vtx_sumpt2.clear();
  t_vtx_trk_index.clear();
 
  if (m_primaryVertices)
    {
      ANA_MSG_DEBUG("processInDet: processing vertices");
 
      t_nvx = m_primaryVertices->size();
 
      static const SG::ConstAccessor<float> sumPt2Acc("sumPt2");
      
      // start of new vertex representation
      t_nvtx = m_primaryVertices->size();
      for (const auto vertex : *m_primaryVertices)
    {
      float vtx_sumpt2 = 0;
      int vtx_ntrk = 0;
 
      t_vtx_type.push_back(vertex->vertexType());
      t_vtx_x.push_back(0);
      t_vtx_y.push_back(0);
      t_vtx_z.push_back(vertex->z());
 
      t_vtx_ntrk.push_back(vtx_ntrk);
      t_vtx_sumpt2.push_back(vtx_sumpt2 / 1e6);
      t_vtx_ntrk_all.push_back(vertex->nTrackParticles());
 
      if (sumPt2Acc.isAvailable(*vertex))
        t_vtx_sumpt2_all.push_back(sumPt2Acc(*vertex));
      else
        t_vtx_sumpt2_all.push_back(-1);
 
      std::vector<int16_t> trk_index;
      if ( m_trackParticles && vertex->nTrackParticles() <= trackLimit )
        {
          const std::vector< ElementLink< xAOD::TrackParticleContainer > >& vxTrackParticles = vertex->trackParticleLinks();
          for (size_t itrk = 0; itrk < vxTrackParticles.size(); itrk++)
        {
          ElementLink< xAOD::TrackParticleContainer > trkLink = vxTrackParticles.at(itrk);
          trk_index.push_back(trkLink.index());
        }
        }
      t_vtx_trk_index.push_back(trk_index);
 
      // end of new vertex representation
 
      if (vertex->vertexType() == xAOD::VxType::PriVtx)
        {
          primary_vertex = vertex;
          pv_index = vx_index;
        }
      if (vertex->vertexType() == xAOD::VxType::PileUp)
        {
          float pileup_sumpT = 0;
          int pileup_nTrack = 0;
          for (size_t itr = 0; itr < vertex->nTrackParticles(); itr++)
        {
          int track_quality = trackQuality(vertex->trackParticle(itr), vertex);
          if (track_quality != -1 && (track_quality & 128) != 0)
            {
              pileup_nTrack++;
              pileup_sumpT += vertex->trackParticle(itr)->pt();
            }
        }
          if (pileup_sumpT > max_pileup_sumpT)
        {
          max_pileup_sumpT = pileup_sumpT;
          max_pileup_nTrack = pileup_nTrack;
          max_pileup_z = vertex->z();
        }
          if (pileup_sumpT > 5e3 || pileup_nTrack > 5) nStrongPileup++;
        }
      vx_index++;
    }
    }
 
  t_nstrong = nStrongPileup;
 
  if (primary_vertex != nullptr)
    {
      t_vx[0] = primary_vertex->x();
      t_vx[1] = primary_vertex->y();
      t_vx[2] = primary_vertex->z();
      /*
    const std::vector<float>& cov = primary_vertex->covariance();
      */
      //for (size_t i=0;i<cov.size();i++)
      for (size_t i = 0; i < 6; i++)
    {
      //t_vxcov[i] = cov.at(i);
      t_vxcov[i] = 0;
    }
      t_vxntrk = primary_vertex->nTrackParticles();
      static const SG::ConstAccessor<float> sumPt2Acc("sumPt2");
      if (sumPt2Acc.isAvailable(*primary_vertex))
    t_vxsumpt2 = sumPt2Acc(*primary_vertex);
      else
    t_vxsumpt2 = 0;
      
      t_vxtype = primary_vertex->vertexType();
      t_pvindex = pv_index;
      t_puvxz = max_pileup_z;
      t_puvxsumpt = max_pileup_sumpT;
      t_puvxntrk = max_pileup_nTrack;
 
      const std::vector< ElementLink< xAOD::TrackParticleContainer > >& vxTrackParticles = primary_vertex->trackParticleLinks();
 
      for (size_t itrk = 0; itrk < vxTrackParticles.size(); itrk++)
    {
      ElementLink< xAOD::TrackParticleContainer > trkLink = vxTrackParticles.at(itrk);
      if (!trkLink.isValid()) continue;
      if (m_trackParticles)
        {
          if (m_trackParticles->size() <= trackLimit)
        t_vx_trk_index.push_back(trkLink.index());
        }
    }
 
    }
 
 
  if (m_trackParticles)
    {
      ANA_MSG_DEBUG("processInDet: processing trackss");
 
      t_trk_pt.clear();
      t_trk_eta.clear();
      t_trk_phi.clear();
      t_trk_e.clear();
      t_trk_index.clear();
      t_trk_theta.clear();
      t_trk_charge.clear();
      t_trk_d0.clear();
      t_trk_z0.clear();
      t_trk_vz.clear();
      t_trk_vtxz.clear();
      t_trk_quality.clear();
      t_trk_nPixHits.clear();
      t_trk_nSctHits.clear();
      t_trk_nPixDead.clear();
      t_trk_nSctDead.clear();
      t_trk_nPixHoles.clear();
      t_trk_nSctHoles.clear();
      t_trk_nTrtHits.clear();
      t_trk_nTrtOutliers.clear();
      t_trk_inPixHits.clear();
      t_trk_exPixHits.clear();
      t_trk_ninPixHits.clear();
      t_trk_nexPixHits.clear();
      t_trk_pixeldEdx.clear();
 
      t_ntrk = m_trackParticles->size();
 
      if ( !enableTracks ) return;
 
      if ( t_ntrk <= trackLimit ) // dump all tracks
    {
      int trk_index = 0;
      for (const auto track : *m_trackParticles)
        {
          writeTrack(track, primary_vertex, trk_index++);
        }
    }
      else // write small vertices
    {
      for (const auto vertex : *m_primaryVertices)
        {
          if (vertex->nTrackParticles() <= trackLimit )
        {
          const std::vector< ElementLink< xAOD::TrackParticleContainer > >& vxTrackParticles = vertex->trackParticleLinks();
          for (size_t itrk = 0; itrk < vxTrackParticles.size(); itrk++)
            {
              ElementLink< xAOD::TrackParticleContainer > trkLink = vxTrackParticles.at(itrk);
              writeTrack(*trkLink, vertex, trkLink.index());
            }
        }
        }
    }
    }
 
  return;
}

processMBTS()

void ZdcNtuple::processMBTS

(

)

Definition at line 1718 of file ZdcNtuple.cxx.

{
  ANA_MSG_DEBUG("processMBTS: trying to process!");
  t_mbts_countA = 0;
  t_mbts_countC = 0;
  t_T2mbts_countAin = 0;
  t_T2mbts_countCin = 0;
  t_mbts_timeA = 0.;
  t_mbts_timeC = 0.;
  t_mbts_timeDiff = 0.;
 
  if (m_mbtsInfo->size() > 0)
  {
    t_mbts_countA = m_mbtsInfo->at(0)->countA();
    t_mbts_countC = m_mbtsInfo->at(0)->countC();
    t_mbts_timeA = m_mbtsInfo->at(0)->timeA();
    t_mbts_timeC = m_mbtsInfo->at(0)->timeC();
    t_mbts_timeDiff = m_mbtsInfo->at(0)->timeDiff();
  }
  else
  {
    ANA_MSG_INFO("processMBTS: Warning: MBTS info empty!");
  }
 
  for (int iside = 0; iside < 2; iside++)
  {
    for (int iin = 0; iin < 8; iin++)
    {
      t_mbts_in_e[iside][iin] = 0.;
      t_mbts_in_t[iside][iin] = 0.;
      t_T2mbts_in_e[iside][iin] = 0.;
      t_T2mbts_in_t[iside][iin] = 0.;
    }
    for (int iout = 0; iout < 4; iout++)
    {
      t_mbts_out_e[iside][iout] = 0.;
      t_mbts_out_t[iside][iout] = 0.;
      t_T2mbts_out_e[iside][iout] = 0.;
      t_T2mbts_out_t[iside][iout] = 0.;
    }
  }
 
  ANA_MSG_DEBUG ("filling MBTS");
 
  if (m_mbtsModules == 0)
  {
    ANA_MSG_INFO("processMBTS: no MBTS container?");
    return;
  }
 
  for (const auto mbtsMod : *m_mbtsModules)
  {
    int iside = 1;
    if (mbtsMod->type() < 0) iside = 0.;
    float phibin = 0.;
    int iphibin = -1;
    if (mbtsMod->eta() > 3)
    {
      phibin = mbtsMod->phi() / (2 * TMath::Pi() / 8.) - 0.4;
      iphibin = static_cast<int>(phibin);
      if (iphibin < 0 || iphibin > 7)
      {
        ANA_MSG_INFO("processMBTS: MBTS has bad phi bin");
        continue;
      }
      t_mbts_in_e[iside][iphibin] = mbtsMod->e();
      t_mbts_in_t[iside][iphibin] = mbtsMod->time();
    }
    else
    {
      phibin = mbtsMod->phi() / (2 * TMath::Pi() / 4.) - 0.24;
      iphibin = static_cast<int>(phibin);
      if (iphibin < 0 || iphibin > 3)
      {
        ANA_MSG_INFO("processMBTS: MBTS has bad phi bin");
        continue;
      }
      t_mbts_out_e[iside][iphibin] = mbtsMod->e();
      t_mbts_out_t[iside][iphibin] = mbtsMod->time();
    }
  }
 
  if (!m_trigT2MbtsBits) return;
 
  for (const auto mbtsBits : *m_trigT2MbtsBits)
  {
    const std::vector<float>& energies = mbtsBits->triggerEnergies();
    const std::vector<float>& times = mbtsBits->triggerTimes();
    for (int imbts = 0; imbts < 32; imbts++)
    {
      int side = imbts / 16;
      int ring = (imbts - 16 * side) / 8;
      bool isInner = (ring == 0);
      int index = (imbts - 16 * side - ring * 8);
      if (!isInner)
      {
        if ((index % 2) != 0) continue; // skip odd out ring
        index /= 2;
      }
      int iside = (side == 0) ? 1 : 0; // code maps side 1 into first 16 bits and side -1 into second set
 
      ANA_MSG_VERBOSE ("imbts=" << imbts << " isInner=" << isInner << " iside=" << iside << " index=" << index << " e=" << energies.at(imbts) << " t=" << times.at(imbts));
      if (iside < 2 and index < 8){ //indices in range?
        if (isInner)
        { 
          t_T2mbts_in_e[iside][index] = energies.at(imbts);
          t_T2mbts_in_t[iside][index] = times.at(imbts);
          if (TMath::Abs(times.at(imbts)) < 12.0 && energies.at(imbts) > 40 / 222.)
          {
            if (iside == 0) t_T2mbts_countCin++;
            if (iside == 1) t_T2mbts_countAin++;
          }
        }
        else
        {
          t_T2mbts_out_e[iside][index] = energies.at(imbts);
          t_T2mbts_out_t[iside][index] = times.at(imbts);
        }
      } //indices check
    }
  }
 
  return;
}

processMCEventCollection()

void ZdcNtuple::processMCEventCollection

(

)

Definition at line 1126 of file ZdcNtuple.cxx.

                                        {
  /******************************************
   * Get the McEventCollection (input)
   ******************************************/
  SG::ReadHandle<McEventCollection> mcEventCollection (m_mcEventCollectionName, getContext());
  if (!mcEventCollection.isValid()){
    ANA_MSG_ERROR("Could not retrieve HepMC with key:" << m_mcEventCollectionName.key());
    return;
  }else{
    ANA_MSG_DEBUG("Retrieved HepMC with key: " << m_mcEventCollectionName.key());
  }
 
  /******************************************
   * Clear and resize the output vectors
  ******************************************/
  t_ZdcTruthParticlePosx.clear();
  t_ZdcTruthParticlePosy.clear();
  t_ZdcTruthParticlePosz.clear();
  t_ZdcTruthParticleTime.clear();
  t_ZdcTruthParticlePx.clear();
  t_ZdcTruthParticlePy.clear();
  t_ZdcTruthParticlePz.clear();
  t_ZdcTruthParticleEnergy.clear();
  t_ZdcTruthParticlePid.clear();
  t_ZdcTruthParticleStatus.clear();
 
  /******************************************
   * Sort the particles into sides and add
   * them to the output vectors
  ******************************************/  
  for (unsigned int cntr = 0; cntr < mcEventCollection->size(); ++cntr){
    const HepMC::GenEvent *genEvt = (*mcEventCollection)[cntr];
#ifdef HEPMC3
    for (const auto &vertex : genEvt->vertices()){
      for (const auto &particle : vertex->particles_in()){
#else
    for (const auto &vertex : genEvt->vertex_range()){
      for (auto ip = vertex->particles_in_const_begin();
           ip != vertex->particles_in_const_end();
           ++ip) {
        auto particle = *ip;
#endif
        t_ZdcTruthParticlePosx.push_back(vertex->position().x());
        t_ZdcTruthParticlePosy.push_back(vertex->position().y());
        t_ZdcTruthParticlePosz.push_back(vertex->position().z());
        t_ZdcTruthParticleTime.push_back(vertex->position().t());
        t_ZdcTruthParticlePx.push_back(particle->momentum().x());
        t_ZdcTruthParticlePy.push_back(particle->momentum().y());
        t_ZdcTruthParticlePz.push_back(particle->momentum().z());
        t_ZdcTruthParticleEnergy.push_back(particle->momentum().e());
        t_ZdcTruthParticlePid.push_back(particle->pdg_id());
        t_ZdcTruthParticleStatus.push_back(particle->status());
      } // end loop over particles
    }// end loop over vertices
  }// end loop over HepMC events
}

processProtons()

void ZdcNtuple::processProtons

(

)

Definition at line 1925 of file ZdcNtuple.cxx.

                              {
 
  proton_pt.clear();
  proton_eta.clear();
  proton_phi.clear();
  proton_e.clear();
  proton_side.clear();
  proton_eLoss.clear();
  proton_t.clear();
 
  proton_track_stationID.clear();
  proton_track_nClusters.clear();
  proton_track_xLocal.clear();
  proton_track_yLocal.clear();
  proton_track_zLocal.clear();
  proton_track_xSlope.clear();
  proton_track_ySlope.clear();
 
  proton_track_stationID.resize(m_afpProtons->size(), std::vector<int>(2));
  proton_track_nClusters.resize(m_afpProtons->size(), std::vector<int>(2));
  proton_track_xLocal.resize(m_afpProtons->size(), std::vector<float>(2));
  proton_track_yLocal.resize(m_afpProtons->size(), std::vector<float>(2));
  proton_track_zLocal.resize(m_afpProtons->size(), std::vector<float>(2));
  proton_track_xSlope.resize(m_afpProtons->size(), std::vector<float>(2));
  proton_track_ySlope.resize(m_afpProtons->size(), std::vector<float>(2));
 
  nProtons = 0;
 
  for(const auto * proton: *m_afpProtons){
 
    proton_pt.push_back(proton->pt());
    proton_eta.push_back(proton->eta());
    proton_phi.push_back(proton->phi());
    proton_e.push_back(proton->e());
    proton_side.push_back(proton->side());
 
    proton_eLoss.push_back((6800.-proton->e())/6800.);
    p_scat.SetPtEtaPhiE(proton->pt(), proton->eta(), proton->phi(), proton->e());
    (signbit(proton->eta())) ? p_beam.SetPxPyPzE(0.0, 0.0, -6800.0, 6800.0) : p_beam.SetPxPyPzE(0.0, 0.0, 6800.0,\
                                                6800.0);
 
    proton_t.push_back( (p_beam - p_scat)*(p_beam - p_scat));
 
    for(int i=0; i< int(proton->nTracks()); i++){
 
      proton_track_stationID.at(nProtons).at(i) = proton->track(i)->stationID();
      proton_track_nClusters.at(nProtons).at(i) = proton->track(i)->nClusters();
      proton_track_xLocal.at(nProtons).at(i) = proton->track(i)->xLocal();
      proton_track_yLocal.at(nProtons).at(i) = proton->track(i)->yLocal();
      proton_track_zLocal.at(nProtons).at(i) = proton->track(i)->zLocal();
      proton_track_xSlope.at(nProtons).at(i) = proton->track(i)->xSlope();
      proton_track_ySlope.at(nProtons).at(i) = proton->track(i)->ySlope();
 
    }
 
    nProtons++;
  }
 
  return;
}

processTriggerDecision()

bool ZdcNtuple::processTriggerDecision

(

)

Definition at line 1183 of file ZdcNtuple.cxx.

{
  ANA_MSG_DEBUG ("Processing trigger");
 
  bool passTrigger = false;
 
  t_trigger = 0;
  t_trigger_TBP = 0;
 
  for (int i = 0; i < 16; i++)
  {
    t_tav[i] = 0;
    t_tbp[i] = 0;
  }
 
  if (m_trigDecision)
    {
      for (int i = 0; i < 16; i++)
    {
      t_tbp[i] = m_trigDecision->tbp().at(i);
      t_tav[i] = m_trigDecision->tav().at(i);     
      ANA_MSG_DEBUG( "TD: " << i << " tbp: " << std::hex << t_tbp[i] << "\t" << t_tav[i] );
    }
      t_bunchGroup = m_trigDecision->bgCode();
    }
 
  if (enableTrigger)
  {
 
    int ic = 0;
    for (auto cg : m_chainGroups)
    {
      
      if (zdcCalib)
    {
      std::string name = cg->getListOfTriggers().at(0);
      const unsigned int triggerbits = m_trigDecisionTool->isPassedBits(name);
      // deferred functionality
      //bool tbp = triggerbits&TrigDefs::L1_isPassedBeforePrescale;
      //bool tap = triggerbits&TrigDefs::L1_isPassedAfterPrescale;
      bool tav = triggerbits&TrigDefs::L1_isPassedAfterVeto;
      ANA_MSG_DEBUG("TD: checking trigger name=" << name<< " tav=" << tav);
      if (tav)
        {
          t_trigger += (1 << ic);
          t_decisions[ic] = true;
          t_prescales[ic] = cg->getPrescale();
          passTrigger = true;         
        }
      else
        {
          t_decisions[ic] = 0;
          t_prescales[ic] = 0;
        }
    }
      else
    {
      if (cg->isPassed())
        {
          t_trigger += (1 << ic);
          t_decisions[ic] = true;
          t_prescales[ic] = cg->getPrescale();
          passTrigger = true;
        }
      else
        {
          t_decisions[ic] = 0;
          t_prescales[ic] = 0;
        }
    }
 
      if (cg->isPassedBits()&TrigDefs::EF_passedRaw)
      {
        t_trigger_TBP += (1 << ic);
      }
 
 
      ic++;
    }
 
    int irc = 0;
    for (auto cg : m_rerunChainGroups)
    {
      t_rerunDecisions[irc] = false;
      if (cg->isPassedBits()&TrigDefs::EF_passedRaw)
      {
        t_rerunDecisions[irc] = true;
      }
      irc++;
    }
 
  }
 
  return passTrigger;
}

processTriggerTowers()

void ZdcNtuple::processTriggerTowers

(

)

processVInjInfo()

void ZdcNtuple::processVInjInfo

(

)

Definition at line 1304 of file ZdcNtuple.cxx.

                               {
    // Check for new run number
  //
  if (t_runNumber != m_lastRunNumber) {
    //
    // Get access to the injector pulse steps for this run
    //
    m_injMapRunToken = m_zdcInjPulserAmpMap->lookupRun(t_runNumber, true);
    if (!m_injMapRunToken.isValid()) {
      ANA_MSG_ERROR("Unable to obtain injector pulse steps for run " << t_runNumber);
    }
    else {
      unsigned int startLB = m_zdcInjPulserAmpMap->getFirstLumiBlock(m_injMapRunToken);
      unsigned int nsteps = m_zdcInjPulserAmpMap->getNumSteps(m_injMapRunToken);
      ANA_MSG_DEBUG("Successfully obtained injector pulse steps for run " << t_runNumber
            << ", first LB = " << startLB << ", number of steps = " << nsteps);
    }
 
    // update the last run number to be the current run number
    m_lastRunNumber = t_runNumber;
  }
  
  t_vInj = m_zdcInjPulserAmpMap->getPulserAmplitude(m_injMapRunToken, t_lumiBlock);
}

processZdcNtupleFromModules()

void ZdcNtuple::processZdcNtupleFromModules

(

)

Definition at line 658 of file ZdcNtuple.cxx.

{
  
  SG::ReadHandle<xAOD::ZdcModuleContainer> zdcModules (m_zdcModuleContainerName);
  SG::ReadHandle<xAOD::ZdcModuleContainer> zdcSums (m_zdcSumContainerName);
  
  ANA_MSG_DEBUG ("copying already processed info!");
  
  //Reset the truth separately since it has a different range
  for(int iside : {0,1}){
    for(int imod = 0; imod < 7; ++imod){
      t_ZdcModuleTruthTotal[iside][imod] = 0; 
      t_ZdcModuleTruthInvis[iside][imod] = 0; 
      t_ZdcModuleTruthEM[iside][imod] = 0; 
      t_ZdcModuleTruthNonEM[iside][imod] = 0; 
      t_ZdcModuleTruthEscaped[iside][imod] = 0;
      t_ZdcModuleTruthNphotons[iside][imod] = 0;
    }
    for(int ch = 0; ch < 16; ++ch){
      t_RpdModuleTruthNphotons[iside][ch] = 0;
    }
  }
  
  for (size_t iside = 0; iside < 2; iside++)
    {
      t_ZdcAmp[iside] = 0; t_ZdcEnergy[iside] = 0; t_ZdcEnergyErr[iside] = 0;t_ZdcTime[iside] = 0; t_ZdcStatus[iside] = 0;
      t_ZdcNLEnergy[iside] = 0;t_ZdcNLEnergyErr[iside] = 0;
      t_ZdcTrigEff[iside] = 0;t_ZdcLucrodTriggerSideAmp[iside] = 0; t_ZdcLucrodTriggerSideAmpLG[iside] = 0; t_ZdcTruthTotal[iside] = 0;
      t_ZdcTruthInvis[iside] = 0; t_ZdcTruthEM[iside] = 0; t_ZdcTruthNonEM[iside] = 0;
      t_ZdcTruthEscaped[iside] = 0;
      for (int imod = 0; imod < 4; imod++)
    {
      t_ZdcModuleAmp[iside][imod] = 0; t_ZdcModuleTime[iside][imod] = 0; t_ZdcModuleStatus[iside][imod] = 0;
      
      t_ZdcModuleCalibAmp[iside][imod] = 0; t_ZdcModuleCalibTime[iside][imod] = 0; t_ZdcModuleChisq[iside][imod] = 0; t_ZdcModuleFitAmp[iside][imod] = 0;
      t_ZdcModuleFitT0[iside][imod] = 0; t_ZdcModuleBkgdMaxFraction[iside][imod] = 0; t_ZdcModuleAmpError[iside][imod] = 0;
      t_ZdcModuleMinDeriv2nd[iside][imod] = 0; t_ZdcModulePresample[iside][imod] = 0; t_ZdcModulePreSampleAmp[iside][imod] = 0;
      t_ZdcLucrodTriggerAmp[iside][imod] = 0;t_ZdcLucrodTriggerAmpLG[iside][imod] = 0;
      t_ZdcModuleMaxADC[iside][imod] = 0; t_ZdcModuleMaxADCHG[iside][imod] = 0; t_ZdcModuleMaxADCLG[iside][imod] = 0;
      t_ZdcModulePeakADCHG[iside][imod] = 0; t_ZdcModulePeakADCLG[iside][imod] = 0;
      t_ZdcModuleAmpLGRefit[iside][imod] = 0; t_ZdcModuleAmpCorrLGRefit[iside][imod] = 0; 
      t_ZdcModuleT0LGRefit[iside][imod] = 0; t_ZdcModuleT0SubLGRefit[iside][imod] = 0; t_ZdcModuleChisqLGRefit[iside][imod] = 0;
      
      if (enableOutputSamples)
        {
          for (int ig=0;ig<2;ig++)
        {
          for (int id=0;id<2;id++)
            {
              for (unsigned int isamp=0;isamp<nsamplesZdc;isamp++)
            {
              if (nsamplesZdc==7) t_raw7[iside][imod][ig][id][isamp]=0;
              if (nsamplesZdc==15) t_raw15[iside][imod][ig][id][isamp]=0;
              if (nsamplesZdc==24) t_raw24[iside][imod][ig][id][isamp]=0;
              if (nsamplesZdc==32) t_raw32[iside][imod][ig][id][isamp]=0;
              if (nsamplesZdc==40) t_raw40[iside][imod][ig][id][isamp]=0;
            }
            }
        }
          if (nsamplesZdc==24||nsamplesZdc==32||nsamplesZdc==40)
        {
          for (int ch=0;ch<16;ch++)
            {
              for (unsigned int isamp=0;isamp<nsamplesZdc;isamp++)
            {
              t_rpdRaw[iside][ch][isamp]=0;
            }
            }
        }
        }
      if (enableRPD)
        {
          for (int ch = 0; ch < 16; ch++) {
        t_RpdChannelBaseline[iside][ch] = 0;
        std::fill(t_RpdChannelPileupExpFitParams[iside][ch], t_RpdChannelPileupExpFitParams[iside][ch] + 2, 0);
        std::fill(t_RpdChannelPileupStretchedExpFitParams[iside][ch], t_RpdChannelPileupStretchedExpFitParams[iside][ch] + 3, 0);
        std::fill(t_RpdChannelPileupExpFitParamErrs[iside][ch], t_RpdChannelPileupExpFitParamErrs[iside][ch] + 2, 0);
        std::fill(t_RpdChannelPileupStretchedExpFitParamErrs[iside][ch], t_RpdChannelPileupStretchedExpFitParamErrs[iside][ch] + 3, 0);
        t_RpdChannelPileupExpFitMSE[iside][ch] = 0;
        t_RpdChannelPileupStretchedExpFitMSE[iside][ch] = 0;
        t_RpdChannelAmplitude[iside][ch] = 0;
        t_RpdChannelAmplitudeCalib[iside][ch] = 0;
        t_RpdChannelMaxADC[iside][ch] = 0;
        t_RpdChannelMaxADCCalib[iside][ch] = 0;
        t_RpdChannelMaxSample[iside][ch] = 0;
        t_RpdChannelStatus[iside][ch] = 0;
        t_RpdChannelPileupFrac[iside][ch] = 0;
          }
          t_RpdSideStatus[iside] = 0;
        }
      if (enableCentroid)
        {
          t_centroidStatus[iside] = 0;
          std::fill(t_RPDChannelSubtrAmp[iside], t_RPDChannelSubtrAmp[iside] + 16, 0);
          t_RPDSubtrAmpSum[iside] = 0;
          t_xCentroidPreGeomCorPreAvgSubtr[iside] = 0;
          t_yCentroidPreGeomCorPreAvgSubtr[iside] = 0;
          t_xCentroidPreAvgSubtr[iside] = 0;
          t_yCentroidPreAvgSubtr[iside] = 0;
          t_xCentroid[iside] = 0;
          t_yCentroid[iside] = 0;
          std::fill(t_xRowCentroid[iside], t_xRowCentroid[iside] + 4, 0);
          std::fill(t_yColCentroid[iside], t_yColCentroid[iside] + 4, 0);
          t_reactionPlaneAngle[iside] = 0;
        }
    }
    }
  
  t_ZdcModuleMask = 0;
  if (enableCentroid) {
    t_centroidDecorationsAvailable = false;
    t_centroidEventValid = false;
    t_cosDeltaReactionPlaneAngle = 0;
  }
 
  /*
  if (t_zdcEventInfoError == xAOD::EventInfo::Error)
    {
      ANA_MSG_INFO("ZDC event failed EventInfo error check - aborting!");
      return;
    }
  */
 
  static const SG::ConstAccessor<char> centroidEventValidAcc("centroidEventValid" + auxSuffix);
  static const SG::ConstAccessor<float> cosDeltaReactionPlaneAngleAcc("cosDeltaReactionPlaneAngle" + auxSuffix);
  static const SG::ConstAccessor<float> CalibEnergyAcc("CalibEnergy"+auxSuffix);
  static const SG::ConstAccessor<float> CalibEnergyErrAcc("CalibEnergyErr"+auxSuffix);
  static const SG::ConstAccessor<float> NLCalibEnergyAcc("NLCalibEnergy"+auxSuffix);
  static const SG::ConstAccessor<float> NLCalibEnergyErrAcc("NLCalibEnergyErr"+auxSuffix);
  static const SG::ConstAccessor<float> UncalibSumAcc("UncalibSum"+auxSuffix);
  static const SG::ConstAccessor<float> UncalibSumErrAcc("UncalibSumErr"+auxSuffix);
  static const SG::ConstAccessor<float> AverageTimeAcc("AverageTime"+auxSuffix);
  static const SG::ConstAccessor<unsigned int> StatusAcc("Status"+auxSuffix);
  static const SG::ConstAccessor<unsigned int> ModuleMaskAcc("ModuleMask"+auxSuffix);
  static const SG::ConstAccessor<float> TruthTotalEnergyAcc("TruthTotalEnergy" + auxSuffix);
  static const SG::ConstAccessor<float> TruthInvisibleEnergyAcc("TruthInvisibleEnergy" + auxSuffix);
  static const SG::ConstAccessor<float> TruthEMEnergyAcc("TruthEMEnergy" + auxSuffix);
  static const SG::ConstAccessor<float> TruthNonEMEnergyAcc("TruthNonEMEnergy" + auxSuffix);
  static const SG::ConstAccessor<float> TruthEscapedEnergyAcc("TruthEscapedEnergy" + auxSuffix);
  static const SG::ConstAccessor<unsigned int> centroidStatusAcc("centroidStatus" + auxSuffix);
  static const SG::ConstAccessor<std::vector<float> > RPDChannelSubtrAmpAcc("RPDChannelSubtrAmp" + auxSuffix);
  static const SG::ConstAccessor<float> RPDSubtrAmpSumAcc("RPDSubtrAmpSum" + auxSuffix);
  static const SG::ConstAccessor<float> xCentroidPreGeomCorPreAvgSubtrAcc("xCentroidPreGeomCorPreAvgSubtr" + auxSuffix);
  static const SG::ConstAccessor<float> yCentroidPreGeomCorPreAvgSubtrAcc("yCentroidPreGeomCorPreAvgSubtr" + auxSuffix);
  static const SG::ConstAccessor<float> xCentroidPreAvgSubtrAcc("xCentroidPreAvgSubtr" + auxSuffix);
  static const SG::ConstAccessor<float> yCentroidPreAvgSubtrAcc("yCentroidPreAvgSubtr" + auxSuffix);
  static const SG::ConstAccessor<float> xCentroidAcc("xCentroid" + auxSuffix);
  static const SG::ConstAccessor<float> yCentroidAcc("yCentroid" + auxSuffix);
  static const SG::ConstAccessor<std::vector<float> > xRowCentroidAcc("xRowCentroid" + auxSuffix);
  static const SG::ConstAccessor<std::vector<float> > yColCentroidAcc("yColCentroid" + auxSuffix);
  static const SG::ConstAccessor<float> reactionPlaneAngleAcc("reactionPlaneAngle" + auxSuffix);
  static const SG::ConstAccessor<unsigned int> RPDStatusAcc("RPDStatus" + auxSuffix);
  static const SG::ConstAccessor<unsigned int> nPhotonsAcc("nPhotons" + auxSuffix);
  static const SG::ConstAccessor<float> CalibTimeAcc("CalibTime" + auxSuffix);
  static const SG::ConstAccessor<float> AmplitudeAcc("Amplitude" + auxSuffix);
  static const SG::ConstAccessor<float> TimeAcc("Time" + auxSuffix);
  static const SG::ConstAccessor<float> ChisqAcc("Chisq" + auxSuffix);
  static const SG::ConstAccessor<float> FitAmpAcc("FitAmp" + auxSuffix);
  static const SG::ConstAccessor<float> FitAmpErrorAcc("FitAmpError" + auxSuffix);
  static const SG::ConstAccessor<float> FitT0Acc("FitT0" + auxSuffix);
  static const SG::ConstAccessor<float> BkgdMaxFractionAcc("BkgdMaxFraction" + auxSuffix);
  static const SG::ConstAccessor<float> MinDeriv2ndAcc("MinDeriv2nd" + auxSuffix);
  static const SG::ConstAccessor<float> PresampleAcc("Presample" + auxSuffix);
  static const SG::ConstAccessor<float> PreSampleAmpAcc("PreSampleAmp" + auxSuffix);
  static const SG::ConstAccessor<float> RPDChannelBaselineAcc("RPDChannelBaseline" + auxSuffix);
  static const SG::ConstAccessor<std::vector<float> > RPDChannelPileupExpFitParamsAcc("RPDChannelPileupExpFitParams" + auxSuffix);
  static const SG::ConstAccessor<std::vector<float> > RPDChannelPileupExpFitParamErrsAcc("RPDChannelPileupExpFitParamErrs" + auxSuffix);
  static const SG::ConstAccessor<std::vector<float> > RPDChannelPileupStretchedExpFitParamsAcc("RPDChannelPileupStretchedExpFitParams" + auxSuffix);
  static const SG::ConstAccessor<std::vector<float> > RPDChannelPileupStretchedExpFitParamErrsAcc("RPDChannelPileupStretchedExpFitParamErrs" + auxSuffix);
  static const SG::ConstAccessor<float> RPDChannelPileupExpFitMSEAcc("RPDChannelPileupExpFitMSE" + auxSuffix);
  static const SG::ConstAccessor<float> RPDChannelPileupStretchedExpFitMSEAcc("RPDChannelPileupStretchedExpFitMSE" + auxSuffix);
  static const SG::ConstAccessor<float> RPDChannelAmplitudeAcc("RPDChannelAmplitude" + auxSuffix);
  static const SG::ConstAccessor<float> RPDChannelAmplitudeCalibAcc("RPDChannelAmplitudeCalib" + auxSuffix);
  static const SG::ConstAccessor<float> RPDChannelMaxADCAcc("RPDChannelMaxADC" + auxSuffix);
  static const SG::ConstAccessor<float> RPDChannelMaxADCCalibAcc("RPDChannelMaxADCCalib" + auxSuffix);
  static const SG::ConstAccessor<unsigned int> RPDChannelMaxSampleAcc("RPDChannelMaxSample" + auxSuffix);
  static const SG::ConstAccessor<unsigned int> RPDChannelStatusAcc("RPDChannelStatus" + auxSuffix);
  static const SG::ConstAccessor<float> RPDChannelPileupFracAcc("RPDChannelPileupFrac" + auxSuffix);
  static const SG::ConstAccessor<float> AmpLGRefitAcc("AmpLGRefit" + auxSuffix);
  static const SG::ConstAccessor<float> T0LGRefitAcc("T0LGRefit" + auxSuffix);
  static const SG::ConstAccessor<float> T0SubLGRefitAcc("T0SubLGRefit" + auxSuffix);
  static const SG::ConstAccessor<float> ChisqLGRefitAcc("ChisqLGRefit" + auxSuffix);
 
  static const SG::ConstAccessor<uint16_t> LucrodTriggerAmpAcc("LucrodTriggerAmp");
  static const SG::ConstAccessor<uint16_t> LucrodTriggerAmpLGAcc("LucrodTriggerAmpLG");
  static const SG::ConstAccessor<uint16_t> LucrodTriggerSideAmpAcc("LucrodTriggerSideAmp");
  static const SG::ConstAccessor<uint16_t> LucrodTriggerSideAmpLGAcc("LucrodTriggerSideAmpLG");
  static const SG::ConstAccessor<float> Amplitude0Acc("Amplitude");
  static const SG::ConstAccessor<float> MaxADCAcc("MaxADC");
  static const SG::ConstAccessor<float> MaxADCHGAcc("MaxADCHG");
  static const SG::ConstAccessor<float> MaxADCLGAcc("MaxADCLG");
  static const SG::ConstAccessor<float> PeakADCHGAcc("PeakADCHG");
  static const SG::ConstAccessor<float> PeakADCLGAcc("PeakADCLG");
  static const SG::ConstAccessor<std::vector<uint16_t> > g0d0DataAcc("g0d0Data");
  static const SG::ConstAccessor<std::vector<uint16_t> > g0d1DataAcc("g0d1Data");
  static const SG::ConstAccessor<std::vector<uint16_t> > g1d0DataAcc("g1d0Data");
  static const SG::ConstAccessor<std::vector<uint16_t> > g1d1DataAcc("g1d1Data");
  static const SG::ConstAccessor<std::vector<uint16_t> > g0dataAcc("g0data");
  static const SG::ConstAccessor<std::vector<uint16_t> > g1dataAcc("g1data");
  
  bool rpdErr = m_eventInfo->isEventFlagBitSet(xAOD::EventInfo::ForwardDet, ZdcEventInfo::RPDDECODINGERROR );
  bool zdcErr = m_eventInfo->isEventFlagBitSet(xAOD::EventInfo::ForwardDet, ZdcEventInfo::ZDCDECODINGERROR );
  t_rpdDecodingError = rpdErr;
  t_zdcDecodingError = zdcErr;
  if (enableCentroid) {
    // locate global sum (only if centroid is needed)
    xAOD::ZdcModule const * globalSum = nullptr;
    for (auto const * zdcSum : *zdcSums) {
      if (zdcSum->zdcSide() == 0) {
        globalSum = zdcSum;
        break;
      }
    }
    if (!globalSum) {
      ANA_MSG_ERROR("unable to locate global ZdcSum (side = 0)");
      t_centroidDecorationsAvailable = false;
    }
    else {
      t_centroidDecorationsAvailable = centroidStatusAcc.isAvailable(*globalSum);
    }
  }
 
  if (rpdErr||zdcErr) ANA_MSG_WARNING( "Decoding errors ZDC=" << zdcErr << " RPD=" << rpdErr );
 
  if (zdcSums.ptr())
    {
      ANA_MSG_DEBUG( "accessing ZdcSums" );
      for (const auto zdcSum : *zdcSums)
    {
      if (zdcSum->zdcSide()==0) {
          // trap new global sum
          if (enableCentroid && t_centroidDecorationsAvailable) {
            t_centroidEventValid = centroidEventValidAcc(*zdcSum);
            t_cosDeltaReactionPlaneAngle = cosDeltaReactionPlaneAngleAcc(*zdcSum);
          }
          // no other branches are filled from global sum - skip to the real sides (C=-1 and A=1)
          continue;
        }
      int iside = 0;
      if (zdcSum->zdcSide() > 0) iside = 1;
      
      //static SG::AuxElement::ConstAccessor< float > acc( "CalibEnergy" );
      //t_ZdcEnergy[iside] = acc(*zdcSum);
      
      if (enableZDC && !zdcErr)
        {
          t_ZdcEnergy[iside] = CalibEnergyAcc(*zdcSum);
          t_ZdcEnergyErr[iside] = CalibEnergyErrAcc(*zdcSum);
          t_ZdcNLEnergy[iside] = NLCalibEnergyAcc(*zdcSum);
          t_ZdcNLEnergyErr[iside] = NLCalibEnergyErrAcc(*zdcSum);
          
          t_ZdcAmp[iside] = UncalibSumAcc(*zdcSum);
          t_ZdcAmpErr[iside] = UncalibSumErrAcc(*zdcSum);
          if (LucrodTriggerSideAmpAcc.isAvailable(*zdcSum))
        t_ZdcLucrodTriggerSideAmp[iside] = LucrodTriggerSideAmpAcc(*zdcSum);
          if (LucrodTriggerSideAmpLGAcc.isAvailable(*zdcSum))
        t_ZdcLucrodTriggerSideAmpLG[iside] = LucrodTriggerSideAmpLGAcc(*zdcSum);
          
          ANA_MSG_VERBOSE("processZdcNtupleFromModules: ZdcSum energy = " << t_ZdcEnergy[iside]);
          
          t_ZdcTime[iside] = AverageTimeAcc(*zdcSum);
          t_ZdcStatus[iside] = StatusAcc(*zdcSum);
          t_ZdcModuleMask += ( ModuleMaskAcc(*zdcSum) << 4 * iside);
          
          if(m_isMC){
        ANA_MSG_DEBUG("Filling sum truth");
        t_ZdcTruthTotal  [iside] = TruthTotalEnergyAcc(*zdcSum);
        t_ZdcTruthInvis  [iside] = TruthInvisibleEnergyAcc(*zdcSum);
        t_ZdcTruthEM     [iside] = TruthEMEnergyAcc(*zdcSum);
        t_ZdcTruthNonEM  [iside] = TruthNonEMEnergyAcc(*zdcSum);
        t_ZdcTruthEscaped[iside] = TruthEscapedEnergyAcc(*zdcSum);
          }
          
        }
      
      if (enableRPD)
        {
          if (enableRPDAmp && !rpdErr)
        {
          t_RpdSideStatus[iside] = RPDStatusAcc(*zdcSum);
        }
          if (enableCentroid && t_centroidDecorationsAvailable)
        {
          t_centroidStatus[iside] = centroidStatusAcc(*zdcSum);
          std::vector<float> const& rpdChannelSubtrAmp = RPDChannelSubtrAmpAcc(*zdcSum);
          std::copy(rpdChannelSubtrAmp.begin(), rpdChannelSubtrAmp.end(), t_RPDChannelSubtrAmp[iside]);
          t_RPDSubtrAmpSum[iside] = RPDSubtrAmpSumAcc(*zdcSum);
          t_xCentroidPreGeomCorPreAvgSubtr[iside] = xCentroidPreGeomCorPreAvgSubtrAcc(*zdcSum);
          t_yCentroidPreGeomCorPreAvgSubtr[iside] = yCentroidPreGeomCorPreAvgSubtrAcc(*zdcSum);
          t_xCentroidPreAvgSubtr[iside] = xCentroidPreAvgSubtrAcc(*zdcSum);
          t_yCentroidPreAvgSubtr[iside] = yCentroidPreAvgSubtrAcc(*zdcSum);
          t_xCentroid[iside] = xCentroidAcc(*zdcSum);
          t_yCentroid[iside] = yCentroidAcc(*zdcSum);
          std::vector<float> const& xRowCentroid = xRowCentroidAcc(*zdcSum);
          std::copy(xRowCentroid.begin(), xRowCentroid.end(), t_xRowCentroid[iside]);
          std::vector<float> const& yColCentroid = yColCentroidAcc(*zdcSum);
          std::copy(yColCentroid.begin(), yColCentroid.end(), t_yColCentroid[iside]);
          t_reactionPlaneAngle[iside] = reactionPlaneAngleAcc(*zdcSum);
        }
        }
    }
    }
  
  ANA_MSG_DEBUG(  "accessing ZdcModules" );
  
  if (zdcModules.ptr())
    {
      for (const auto zdcMod : *zdcModules)
    {
      int iside = 0;
      if (zdcMod->zdcSide() > 0) iside = 1;
      int imod = zdcMod->zdcModule();
      
      if (m_isMC){
        //Calib hits are only stored in channel 0 of the RPD
        if(!(imod == 4 && zdcMod->zdcChannel() != 0)){
          t_ZdcModuleTruthTotal  [iside][imod] = TruthTotalEnergyAcc(*zdcMod);
          t_ZdcModuleTruthInvis  [iside][imod] = TruthInvisibleEnergyAcc(*zdcMod);
          t_ZdcModuleTruthEM     [iside][imod] = TruthEMEnergyAcc(*zdcMod);
          t_ZdcModuleTruthNonEM  [iside][imod] = TruthNonEMEnergyAcc(*zdcMod);
          t_ZdcModuleTruthEscaped[iside][imod] = TruthEscapedEnergyAcc(*zdcMod);
          t_ZdcModuleTruthNphotons[iside][imod] = nPhotonsAcc(*zdcMod);
        }
        //Calib hits are stored for all modules
        //Other data is only valid for module 1-4
        if(imod > 4) continue;
      }
      
      ANA_MSG_VERBOSE ("Module " << zdcMod->zdcSide() << " " << zdcMod->zdcModule() << " amp:" << AmplitudeAcc(*zdcMod));
      
      if (zdcMod->zdcType() == 0 && !zdcErr)
        {
          // ZDC energy type modules
          t_ZdcModuleCalibAmp[iside][imod] = CalibEnergyAcc(*zdcMod);
          t_ZdcModuleCalibTime[iside][imod] = CalibTimeAcc(*zdcMod);
          t_ZdcModuleStatus[iside][imod] = StatusAcc(*zdcMod);
          if (t_ZdcModuleAmp[iside][imod] != 0.)
        Warning("processZdcNtupleFromModules", "overwriting side %d module %d!", iside, imod);
          t_ZdcModuleAmp[iside][imod] = AmplitudeAcc(*zdcMod);
          t_ZdcModuleTime[iside][imod] = TimeAcc(*zdcMod);
          
          t_ZdcModuleChisq[iside][imod] = ChisqAcc(*zdcMod);
          t_ZdcModuleFitAmp[iside][imod] = FitAmpAcc(*zdcMod);
          t_ZdcModuleAmpError[iside][imod] = FitAmpErrorAcc(*zdcMod);
          t_ZdcModuleFitT0[iside][imod] = FitT0Acc(*zdcMod);
          t_ZdcModuleBkgdMaxFraction[iside][imod] = BkgdMaxFractionAcc(*zdcMod);
          t_ZdcModuleMinDeriv2nd[iside][imod] = MinDeriv2ndAcc(*zdcMod);
          t_ZdcModulePresample[iside][imod] = PresampleAcc(*zdcMod);
          t_ZdcModulePreSampleAmp[iside][imod] = PreSampleAmpAcc(*zdcMod);
          
          if (AmpLGRefitAcc.isAvailable(*zdcMod)) {
        t_ZdcModuleAmpLGRefit[iside][imod] = AmpLGRefitAcc(*zdcMod);
        t_ZdcModuleT0LGRefit[iside][imod] = T0LGRefitAcc(*zdcMod);
        t_ZdcModuleT0SubLGRefit[iside][imod] = T0SubLGRefitAcc(*zdcMod);
        t_ZdcModuleChisqLGRefit[iside][imod] = ChisqLGRefitAcc(*zdcMod);
          }
          
          if (LucrodTriggerAmpAcc.isAvailable(*zdcMod))
        t_ZdcLucrodTriggerAmp[iside][imod] = LucrodTriggerAmpAcc(*zdcMod);
          if (LucrodTriggerAmpLGAcc.isAvailable(*zdcMod))
        t_ZdcLucrodTriggerAmpLG[iside][imod] = LucrodTriggerAmpLGAcc(*zdcMod);
 
          if (MaxADCAcc.isAvailable(*zdcMod))
        t_ZdcModuleMaxADC[iside][imod] = MaxADCAcc(*zdcMod);
          if (MaxADCHGAcc.isAvailable(*zdcMod))
        t_ZdcModuleMaxADCHG[iside][imod] = MaxADCHGAcc(*zdcMod);
          if (MaxADCLGAcc.isAvailable(*zdcMod))
        t_ZdcModuleMaxADCLG[iside][imod] = MaxADCLGAcc(*zdcMod);
          if (PeakADCHGAcc.isAvailable(*zdcMod))
        t_ZdcModulePeakADCHG[iside][imod] = PeakADCHGAcc(*zdcMod);
          if (PeakADCLGAcc.isAvailable(*zdcMod))
        t_ZdcModulePeakADCLG[iside][imod] = PeakADCLGAcc(*zdcMod);
 
          if (enableOutputSamples && !zdcErr)
        {
          for (unsigned int isamp = 0; isamp < nsamplesZdc; isamp++) // 7 samples
            {
              if (nsamplesZdc == 7)
            {
              t_raw7[iside][imod][0][0][isamp] = g0d0DataAcc(*zdcMod).at(isamp);
              t_raw7[iside][imod][0][1][isamp] = g0d1DataAcc(*zdcMod).at(isamp);
              t_raw7[iside][imod][1][0][isamp] = g1d0DataAcc(*zdcMod).at(isamp);
              t_raw7[iside][imod][1][1][isamp] = g1d1DataAcc(*zdcMod).at(isamp);
            }
              
              if (nsamplesZdc == 15)
            {
              t_raw15[iside][imod][0][0][isamp] = g0d0DataAcc(*zdcMod).at(isamp);
              t_raw15[iside][imod][0][1][isamp] = g0d1DataAcc(*zdcMod).at(isamp);
              t_raw15[iside][imod][1][0][isamp] = g1d0DataAcc(*zdcMod).at(isamp);
              t_raw15[iside][imod][1][1][isamp] = g1d1DataAcc(*zdcMod).at(isamp);
            }
              
              if (nsamplesZdc == 24)
            {
              t_raw24[iside][imod][0][0][isamp] = g0dataAcc(*zdcMod).at(isamp);
              t_raw24[iside][imod][1][0][isamp] = g1dataAcc(*zdcMod).at(isamp);
            }
              if (nsamplesZdc == 32)
            {
              t_raw32[iside][imod][0][0][isamp] = g0dataAcc(*zdcMod).at(isamp);
              t_raw32[iside][imod][1][0][isamp] = g1dataAcc(*zdcMod).at(isamp);
            }
              if (nsamplesZdc == 40)
            {
              t_raw40[iside][imod][0][0][isamp] = g0dataAcc(*zdcMod).at(isamp);
              t_raw40[iside][imod][1][0][isamp] = g1dataAcc(*zdcMod).at(isamp);
            }
            }
        }
        }
      else if (zdcMod->zdcType() == 1 && nsamplesZdc == 24)
        {
          // this is the RPD
          if (enableRPD)
        {
          if (enableRPDAmp && !rpdErr)
            {
              t_RpdChannelBaseline[iside][zdcMod->zdcChannel()] = RPDChannelBaselineAcc(*zdcMod);
              std::vector<float> const &rpdChannelPileupExpFitParams = RPDChannelPileupExpFitParamsAcc(*zdcMod);
              std::copy(rpdChannelPileupExpFitParams.begin(), rpdChannelPileupExpFitParams.end(), t_RpdChannelPileupExpFitParams[iside][zdcMod->zdcChannel()]);
              std::vector<float> const &rpdChannelPileupExpFitParamErrs = RPDChannelPileupExpFitParamErrsAcc(*zdcMod);
              std::copy(rpdChannelPileupExpFitParamErrs.begin(), rpdChannelPileupExpFitParamErrs.end(), t_RpdChannelPileupExpFitParamErrs[iside][zdcMod->zdcChannel()]);
              std::vector<float> const &rpdChannelPileupStretchedExpFitParams = RPDChannelPileupStretchedExpFitParamsAcc(*zdcMod);
              std::copy(rpdChannelPileupStretchedExpFitParams.begin(), rpdChannelPileupStretchedExpFitParams.end(), t_RpdChannelPileupStretchedExpFitParams[iside][zdcMod->zdcChannel()]);
              std::vector<float> const &rpdChannelPileupStretchedExpFitParamErrs = RPDChannelPileupStretchedExpFitParamErrsAcc(*zdcMod);
              std::copy(rpdChannelPileupStretchedExpFitParamErrs.begin(), rpdChannelPileupStretchedExpFitParamErrs.end(), t_RpdChannelPileupStretchedExpFitParamErrs[iside][zdcMod->zdcChannel()]);
              t_RpdChannelPileupExpFitMSE[iside][zdcMod->zdcChannel()] = RPDChannelPileupExpFitMSEAcc(*zdcMod);
              t_RpdChannelPileupStretchedExpFitMSE[iside][zdcMod->zdcChannel()] = RPDChannelPileupStretchedExpFitMSEAcc(*zdcMod);
              t_RpdChannelAmplitude[iside][zdcMod->zdcChannel()] = RPDChannelAmplitudeAcc(*zdcMod);
              t_RpdChannelAmplitudeCalib[iside][zdcMod->zdcChannel()] = RPDChannelAmplitudeCalibAcc(*zdcMod);
              t_RpdChannelMaxADC[iside][zdcMod->zdcChannel()] = RPDChannelMaxADCAcc(*zdcMod);
              t_RpdChannelMaxADCCalib[iside][zdcMod->zdcChannel()] = RPDChannelMaxADCCalibAcc(*zdcMod);
              t_RpdChannelMaxSample[iside][zdcMod->zdcChannel()] = RPDChannelMaxSampleAcc(*zdcMod);
              t_RpdChannelStatus[iside][zdcMod->zdcChannel()] = RPDChannelStatusAcc(*zdcMod);
              t_RpdChannelPileupFrac[iside][zdcMod->zdcChannel()] = RPDChannelPileupFracAcc(*zdcMod);
              if(m_isMC){
            t_RpdModuleTruthNphotons[iside][zdcMod->zdcChannel()] = nPhotonsAcc(*zdcMod);
              }
            }
          if (enableOutputSamples)
            {
              std::vector<uint16_t> const &rpdChannelRaw = g0dataAcc(*zdcMod);
              std::copy(rpdChannelRaw.begin(), rpdChannelRaw.end(), t_rpdRaw[iside][zdcMod->zdcChannel()]);
            }
        }
        }
    }
    }
  else
    {
      ANA_MSG_INFO("No ZdcModules" << auxSuffix << " when expected!");
    }
  
  if (msgLvl (MSG::VERBOSE))
    {
      std::ostringstream message;
      message << "Dump zdc_ZdcModuleAmp: ";
      for (int iside = 0; iside < 2; iside++)
    {
      for (int imod = 0; imod < 4; imod++)
        {
          message << t_ZdcModuleAmp[iside][imod] << " ";
        }
    }
    }
}

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< Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

reprocessZdcModule()

void ZdcNtuple::reprocessZdcModule	(	const xAOD::ZdcModule *	zdcMod,
		bool	flipdelay = 0 )

requestBeginInputFile()

StatusCode EL::AnaAlgorithm::requestBeginInputFile ( )

inherited

register this algorithm to have an implementation of beginInputFile

Guarantee

strong

Failures

beginInputFile not supported

Definition at line 270 of file AnaAlgorithm.cxx.

  {
    m_hasBeginInputFile = true;
 
#ifndef XAOD_STANDALONE
    // Connect to the IncidentSvc:
    ServiceHandle< IIncidentSvc > incSvc( "IncidentSvc", name() );
    ATH_CHECK( incSvc.retrieve() );
 
    // Set up the right callback, but ensure we don't double-register
    // if we are called twice
    incSvc->removeListener( this, IncidentType::BeginInputFile );
    incSvc->addListener( this, IncidentType::BeginInputFile, 0, true );
#endif
 
    return StatusCode::SUCCESS;
  }

requestEndInputFile()

StatusCode EL::AnaAlgorithm::requestEndInputFile ( )

inherited

register this algorithm to have an implementation of endInputFile

Guarantee

strong

Failures

endInputFile not supported

Definition at line 291 of file AnaAlgorithm.cxx.

  {
    m_hasEndInputFile = true;
 
#ifndef XAOD_STANDALONE
    // Connect to the IncidentSvc:
    ServiceHandle< IIncidentSvc > incSvc( "IncidentSvc", name() );
    ATH_CHECK( incSvc.retrieve() );
 
    // Set up the right callback, but ensure we don't double-register
    // if we are called twice
    incSvc->removeListener( this, IncidentType::EndInputFile );
    incSvc->addListener( this, IncidentType::EndInputFile, 0, true );
#endif
 
    return StatusCode::SUCCESS;
  }

requestFileExecute()

StatusCode EL::AnaAlgorithm::requestFileExecute ( )

inherited

register this algorithm to have an implementation of fileexecute

Guarantee

strong

Failures

fileExecute not supported

Definition at line 249 of file AnaAlgorithm.cxx.

  {
    m_hasFileExecute = true;
 
#ifndef XAOD_STANDALONE
    // Connect to the IncidentSvc:
    ServiceHandle< IIncidentSvc > incSvc( "IncidentSvc", name() );
    ATH_CHECK( incSvc.retrieve() );
 
    // Set up the right callback, but ensure we don't double-register
    // if we are called twice
    incSvc->removeListener( this, IncidentType::BeginInputFile );
    incSvc->addListener( this, IncidentType::BeginInputFile, 0, true );
#endif
 
    return StatusCode::SUCCESS;
  }

setupTriggerHistos()

void ZdcNtuple::setupTriggerHistos

(

)

Definition at line 2017 of file ZdcNtuple.cxx.

{
  if (!enableTrigger) return;
 
  m_outputTree = tree( "zdcTree" );
  ANA_MSG_INFO("setupTriggerHistos(): Setting up trigger histos and ntuple = " << m_outputTree);
 
  std::vector<std::string> triggers;
  std::vector<std::string> rerunTriggers;
  bool zdc_triggers = true;
 
  // ZDC triggers
  if (zdc_triggers)
  {
    if (zdcCalib) // lists for calibration data
      {
    if (lhcf2022)
      {
        triggers.push_back("HLT_noalg_ZDCPEB_L1LHCF");
        triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_OR");
      }
 
    if (pbpb2023)
      {
        triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_OR");
        triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_OR_EMPTY");
        triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_OR_UNPAIRED_NONISO");
        triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_A_C");
        triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_A_C_EMPTY");
        triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_A_C_UNPAIRED_NONISO");
        triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_A");
        triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_A_EMPTY");
        triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_A_UNPAIRED_NONISO");
        triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_C");
        triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_C_EMPTY");
        triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_C_UNPAIRED_NONISO");
      }
      }
    else // lists for physics data
      {
    if (lhcf2022)
          {
            triggers.push_back("HLT_noalg_L1LHCF");
          }
    if (lhcf2022afp)
          {
            triggers.push_back("HLT_noalg_AFPPEB_L1AFP_A");
            triggers.push_back("HLT_noalg_AFPPEB_L1AFP_C");
          }
    if (lhcf2022zdc)
          {
            triggers.push_back("HLT_noalg_ZDCPEB_L1ZDC_OR");
            triggers.push_back("HLT_noalg_ZDCPEB_L1LHCF");
            triggers.push_back("HLT_noalg_L1ZDC_OR");
            triggers.push_back("HLT_noalg_L1ZDC_XOR_E2");
            triggers.push_back("HLT_noalg_L1ZDC_XOR_E1_E3");
            triggers.push_back("HLT_noalg_L1ZDC_A_AND_C");
            triggers.push_back("HLT_mb_sptrk_L1ZDC_OR");
            triggers.push_back("HLT_mb_sptrk_L1ZDC_XOR_E2");
            triggers.push_back("HLT_mb_sptrk_L1ZDC_XOR_E1_E3");
            triggers.push_back("HLT_mb_sptrk_L1ZDC_A_AND_C");
            triggers.push_back("HLT_mb_sp100_trk30_hmt_L1ZDC_XOR_E2");
            triggers.push_back("HLT_mb_sp100_trk30_hmt_L1ZDC_XOR_E1_E3");
            triggers.push_back("HLT_mb_sp100_trk30_hmt_L1ZDC_A_AND_C");
          }
      }
  }
 
  //char name[50];
  ANA_MSG_INFO("Adding trigger branches!");
  
  int ic = 0;
  for (auto &trig : triggers)
    {
      const Trig::ChainGroup* cg = m_trigDecisionTool->getChainGroup(trig);
      if (cg->getListOfTriggers().size())
    {
      ANA_MSG_INFO("setupTriggerHistos(): Trigger found = " << trig.c_str() << " bit " << ic);
    }
      else
    {
      ANA_MSG_INFO("setupTriggerHistos(): Trigger NOT found = " << trig.c_str()  << " bit " << ic);
    }
      m_chainGroups.push_back(cg);
      // force all triggers to show up in tree
      TString bname(trig.c_str());
      m_outputTree->Branch(bname, &(t_decisions[ic]), bname + "/O");
      m_outputTree->Branch("ps_" + bname, &(t_prescales[ic]), "ps_" + bname + "/F");
      ic++;
    }
 
  ANA_MSG_INFO( "triggers = " << triggers.size() << " chains = " << m_chainGroups.size() );
 
  int irc = 0;
  ANA_MSG_INFO("Adding rerun trigger branches!");
  for (auto &trig : rerunTriggers)
    {
      const Trig::ChainGroup* cg = m_trigDecisionTool->getChainGroup(trig);
      m_rerunChainGroups.push_back(cg);
      if (cg->getListOfTriggers().size())
    {
      ANA_MSG_INFO("setupTriggerHistos(): Rerun trigger found = " << trig.c_str() << " bit " << irc);
    }
      else
    {
      ANA_MSG_INFO("setupTriggerHistos(): Rerun trigger NOT found = " << trig.c_str()  << " bit " << irc);
    }
      // force all rerun triggers to show up in tree
      TString bname(trig.c_str());
      m_outputTree->Branch(bname, &(t_rerunDecisions[irc]), bname + "/O");
      irc++;
    }
  
  // trigger matching flags for electrons and muons
  
  m_setupTrigHist = true;
 
  ANA_MSG_INFO("setupTriggerHistos(): Finished setting up trigger");
 
}

sysInitialize()

StatusCode AthHistogramAlgorithm::sysInitialize ( )

virtualinherited

Initialization method invoked by the framework.

This method is responsible for any bookkeeping of initialization required by the framework itself. It will in turn invoke the initialize() method of the derived algorithm, and of any sub-algorithms which it creates.

Reimplemented from AthAlgorithm.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 75 of file AthHistogramAlgorithm.cxx.

{
  // ---- stolen from GaudiKernel/Algorithm::sysInitialize -------
  // Bypass the initialization if the algorithm
  // has already been initialized.
  if ( Gaudi::StateMachine::INITIALIZED <= FSMState() ) return StatusCode::SUCCESS;
 
  // Set the Algorithm's properties
  bindPropertiesTo( serviceLocator()->getOptsSvc() );
 
  // Bypass the initialization if the algorithm is disabled.
  // Need to do this after bindPropertiesTo.
  if ( !isEnabled( ) ) return StatusCode::SUCCESS;
  
  // ---- stolen from GaudiKernel/Algorithm::sysInitialize ------- END ---
 
 
  // Get the THistSvc
  ATH_CHECK ( histSvc().retrieve() );
 
  // Configure the underlying AthHistogramming helper
  ATH_CHECK ( AthHistogramming::configAthHistogramming( m_histSvc,
                                                        m_prefix,          m_rootDir, 
                                                        m_histNamePrefix,  m_histNamePostfix,
                                                        m_histTitlePrefix, m_histTitlePostfix ) );
  
  // Print some setup information into the log file
  ATH_MSG_DEBUG ("Initializing " << name() << "...");
  ATH_MSG_DEBUG (" using THistService     = " << m_histSvc );
  ATH_MSG_DEBUG (" using RootStreamName   = " << m_prefix );
  ATH_MSG_DEBUG (" using RootDirName      = " << m_rootDir );
  ATH_MSG_DEBUG (" using HistNamePrefix   = " << m_histNamePrefix );
  ATH_MSG_DEBUG (" using HistNamePostfix  = " << m_histNamePostfix );
  ATH_MSG_DEBUG (" using HistTitlePrefix  = " << m_histTitlePrefix );
  ATH_MSG_DEBUG (" using HistTitlePostfix = " << m_histTitlePostfix );
 
 
  // re-direct to base class...
  return AthAlgorithm::sysInitialize();
}

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< 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.

trackQuality()

int ZdcNtuple::trackQuality	(	const xAOD::TrackParticle *	tp,
		const xAOD::Vertex *	vertex )

Definition at line 1548 of file ZdcNtuple.cxx.

{
 
  if (!track) return -1;
 
  bool pass_looseprimary = false;
  if (m_selTool->accept(*track,vertex))
    {
      pass_looseprimary = true;      
    }
 
  int quality = 0;
  if (pass_looseprimary) quality += 1;
 
  return quality;
 
}

tree()

TTree * AthHistogramming::tree ( const std::string & treeName, const std::string & tDir = "", const std::string & stream = "" )

protectedinherited

Simplify the retrieval of registered TTrees.

Definition at line 378 of file AthHistogramming.cxx.

{
  // Build a 32 bit hash out of the name
  const hash_t treeHash = this->hash(treeName);
 
  // See if this entry exists in the map
  TreeMap_t::const_iterator it = m_treeMap.find( treeHash );
  if ( it == m_treeMap.end() ) // It doesn't exist!
    { // Let's see into the THistSvc if somebody else has registered the TTree...
 
      // Need to copy the strings as we will massage them from here on
      std::string treeNameCopy = treeName;
      std::string tDirCopy     = tDir;
      std::string streamCopy   = stream;
 
      // Massage the final string to book things
      std::string bookingString("");
      this->buildBookingString( bookingString, treeNameCopy, tDirCopy, streamCopy );
 
      TTree* treePointer(NULL);
      if ( !((histSvc()->getTree(bookingString, treePointer)).isSuccess()) )
        {
          m_msg << MSG::WARNING
                << "Problem retrieving the TTree with name " << treeNameCopy
                << " in " << m_name << "... it doesn't exist, neither in the cached map nor in the THistSvc!"
                << " Will return an NULL pointer... you have to handle it correctly!" << endmsg;
          return NULL;
        }
      // If we get to here, we actually found the TTree in the THistSvc.
      // So let's add it to the local cache map and return its pointer
      m_treeMap.insert( m_treeMap.end(), std::pair< const hash_t, TTree* >( treeHash, treePointer ) );
      return treePointer;
    }
 
  // Return the pointer to the TTree that we got from the local cache map
  return it->second;
}

updateVHKA()

void AthCommonDataStore< AthCommonMsg< 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);
      }
    }
  }

writeTrack()

void ZdcNtuple::writeTrack	(	const xAOD::TrackParticle *	track,
		const xAOD::Vertex *	vertex,
		int	trk_index )

Definition at line 1568 of file ZdcNtuple.cxx.

{
  t_trk_pt.push_back(track->pt());
  t_trk_eta.push_back(track->eta());
  t_trk_phi.push_back(track->phi());
  t_trk_e.push_back(track->e());
  t_trk_index.push_back(trk_index);
  t_trk_theta.push_back(track->theta());
  t_trk_charge.push_back(track->charge());
  t_trk_d0.push_back(track->d0());
  t_trk_z0.push_back(track->z0());
  t_trk_vz.push_back(track->vz());
 
  float vtxz = -999.;
  t_trk_vtxz.push_back(vtxz);
 
  static const SG::ConstAccessor<uint8_t> numberOfInnermostPixelLayerHitsAcc("numberOfInnermostPixelLayerHits");
  static const SG::ConstAccessor<uint8_t> expectInnermostPixelLayerHitAcc("expectInnermostPixelLayerHit");
  static const SG::ConstAccessor<uint8_t> numberOfNextToInnermostPixelLayerHitsAcc("numberOfNextToInnermostPixelLayerHits");
  static const SG::ConstAccessor<uint8_t> expectNextToInnermostPixelLayerHitAcc("expectNextToInnermostPixelLayerHit");
  static const SG::ConstAccessor<uint8_t> numberOfSCTHitsAcc("numberOfSCTHits");
  static const SG::ConstAccessor<uint8_t> numberOfPixelHitsAcc("numberOfPixelHits");
  static const SG::ConstAccessor<uint8_t> numberOfSCTDeadSensorsAcc("numberOfSCTDeadSensors");
  static const SG::ConstAccessor<uint8_t> numberOfPixelDeadSensorsAcc("numberOfPixelDeadSensors");
  static const SG::ConstAccessor<uint8_t> numberOfSCTHolesAcc("numberOfSCTHoles");
  static const SG::ConstAccessor<uint8_t> numberOfPixelHolesAcc("numberOfPixelHoles");
  static const SG::ConstAccessor<uint8_t> numberOfTRTHitsAcc("numberOfTRTHits");
  static const SG::ConstAccessor<uint8_t> numberOfTRTOutliersAcc("numberOfTRTOutliers");
 
  t_trk_quality.push_back(trackQuality(track, vertex));
  t_trk_inPixHits.push_back(numberOfInnermostPixelLayerHitsAcc(*track));
  t_trk_exPixHits.push_back(expectInnermostPixelLayerHitAcc(*track));
  t_trk_ninPixHits.push_back(numberOfNextToInnermostPixelLayerHitsAcc(*track));
  t_trk_nexPixHits.push_back(expectNextToInnermostPixelLayerHitAcc(*track));
  t_trk_nSctHits.push_back(numberOfSCTHitsAcc(*track));
  t_trk_nPixHits.push_back(numberOfPixelHitsAcc(*track));
  t_trk_nSctDead.push_back(numberOfSCTDeadSensorsAcc(*track));
  t_trk_nPixDead.push_back(numberOfPixelDeadSensorsAcc(*track));
  t_trk_nSctHoles.push_back(numberOfSCTHolesAcc(*track));
  t_trk_nPixHoles.push_back(numberOfPixelHolesAcc(*track));
  t_trk_nTrtHits.push_back(numberOfTRTHitsAcc(*track));
  t_trk_nTrtOutliers.push_back(numberOfTRTOutliersAcc(*track));
 
  float pixeldEdx = 0;
  track->summaryValue(pixeldEdx, xAOD::SummaryType::pixeldEdx);
  t_trk_pixeldEdx.push_back(pixeldEdx);
}

Member Data Documentation

auxSuffix

std::string ZdcNtuple::auxSuffix

Definition at line 80 of file ZdcNtuple.h.

doZdcCalib

bool ZdcNtuple::doZdcCalib

Definition at line 91 of file ZdcNtuple.h.

enableCalo

bool ZdcNtuple::enableCalo

Definition at line 61 of file ZdcNtuple.h.

enableCentroid

bool ZdcNtuple::enableCentroid

Definition at line 89 of file ZdcNtuple.h.

enableClusters

bool ZdcNtuple::enableClusters

Definition at line 62 of file ZdcNtuple.h.

enableElectrons

bool ZdcNtuple::enableElectrons

Definition at line 65 of file ZdcNtuple.h.

enableID

bool ZdcNtuple::enableID

Definition at line 60 of file ZdcNtuple.h.

enableJets

bool ZdcNtuple::enableJets

Definition at line 69 of file ZdcNtuple.h.

enableMuons

bool ZdcNtuple::enableMuons

Definition at line 64 of file ZdcNtuple.h.

enableOutputSamples

bool ZdcNtuple::enableOutputSamples

Definition at line 57 of file ZdcNtuple.h.

enableOutputTree

bool ZdcNtuple::enableOutputTree

Definition at line 56 of file ZdcNtuple.h.

enablePhotons

bool ZdcNtuple::enablePhotons

Definition at line 66 of file ZdcNtuple.h.

enableRPD

bool ZdcNtuple::enableRPD

Definition at line 87 of file ZdcNtuple.h.

enableRPDAmp

bool ZdcNtuple::enableRPDAmp

Definition at line 88 of file ZdcNtuple.h.

enableTracks

bool ZdcNtuple::enableTracks

Definition at line 63 of file ZdcNtuple.h.

enableTrigger

bool ZdcNtuple::enableTrigger

Definition at line 58 of file ZdcNtuple.h.

enableTriggerJets

bool ZdcNtuple::enableTriggerJets

Definition at line 70 of file ZdcNtuple.h.

enableTruth

bool ZdcNtuple::enableTruth

Definition at line 68 of file ZdcNtuple.h.

enableTT

bool ZdcNtuple::enableTT

Definition at line 67 of file ZdcNtuple.h.

enableZDC

bool ZdcNtuple::enableZDC

Definition at line 86 of file ZdcNtuple.h.

flipDelay

bool ZdcNtuple::flipDelay

Definition at line 78 of file ZdcNtuple.h.

grlFilename

std::string ZdcNtuple::grlFilename

Definition at line 55 of file ZdcNtuple.h.

h_TCSigCut

TH1* ZdcNtuple::h_TCSigCut

Definition at line 336 of file ZdcNtuple.h.

h_zdcTriggers

TH1* ZdcNtuple::h_zdcTriggers

Definition at line 144 of file ZdcNtuple.h.

h_zdcTriggersTBP

TH1* ZdcNtuple::h_zdcTriggersTBP

Definition at line 145 of file ZdcNtuple.h.

lhcf2022

bool ZdcNtuple::lhcf2022

Definition at line 82 of file ZdcNtuple.h.

lhcf2022afp

bool ZdcNtuple::lhcf2022afp

Definition at line 84 of file ZdcNtuple.h.

lhcf2022zdc

bool ZdcNtuple::lhcf2022zdc

Definition at line 83 of file ZdcNtuple.h.

m_afpProtons

const xAOD::AFPProtonContainer* ZdcNtuple::m_afpProtons

Definition at line 127 of file ZdcNtuple.h.

m_caloClusters

const xAOD::CaloClusterContainer* ZdcNtuple::m_caloClusters

Definition at line 122 of file ZdcNtuple.h.

m_caloSums

const xAOD::HIEventShapeContainer* ZdcNtuple::m_caloSums

Definition at line 116 of file ZdcNtuple.h.

m_chainGroups

std::vector<const Trig::ChainGroup*> ZdcNtuple::m_chainGroups

Definition at line 136 of file ZdcNtuple.h.

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_effMap

EffMap_t AthHistogramming::m_effMap

privateinherited

The map of histogram names to their pointers.

Definition at line 208 of file AthHistogramming.h.

m_eventCounter

int ZdcNtuple::m_eventCounter

Definition at line 132 of file ZdcNtuple.h.

m_eventInfo

const xAOD::EventInfo* ZdcNtuple::m_eventInfo

Definition at line 112 of file ZdcNtuple.h.

m_eventShapes

const xAOD::HIEventShapeContainer* ZdcNtuple::m_eventShapes

Definition at line 117 of file ZdcNtuple.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects

privateinherited

Definition at line 79 of file AthAlgorithm.h.

m_gapPtMin

float ZdcNtuple::m_gapPtMin

Definition at line 325 of file ZdcNtuple.h.

m_gapThresholds

float ZdcNtuple::m_gapThresholds[98]

Initial value:

=
    {4.7426,5.11018,5.07498,5.0969,5.10695,5.04098,5.07106,4.98087,5.11647,5.08988,5.16267,
     5.17202,5.23803,5.25314,5.29551,5.35092,5.40863,5.44375,5.38075,5.25022,5.37933,5.25459,5.37719,5.25169,5.73985,
     5.79174,5.79266,5.79588,5.7963,5.81949,5.82273,5.85658,5.85442,5.84779,5.77679,5.83323,5.84524,5.84439,5.84488,
     5.84744,5.84683,5.84524,5.84594,5.84656,5.84639,5.84461,5.84515,5.84206,5.8396,5.84497,5.84801,5.84608,5.84608,
     5.84783,5.84726,5.84844,5.8477,5.84796,5.84757,5.84822,5.84814,5.84617,5.83451,5.77658,5.84309,5.85496,5.85761,
     5.82555,5.82206,5.78982,5.78482,5.7778,5.78327,5.74898,5.25459,5.37503,5.25459,5.37283,5.25169,5.37862,5.44473,
     5.41041,5.34498,5.29551,5.25602,5.2283,5.17428,5.14504,5.09342,5.12256,4.98721,5.07106,5.02642,5.10031,5.11018,
     5.05447,5.10031,4.7426}

Definition at line 327 of file ZdcNtuple.h.

    {4.7426,5.11018,5.07498,5.0969,5.10695,5.04098,5.07106,4.98087,5.11647,5.08988,5.16267,
     5.17202,5.23803,5.25314,5.29551,5.35092,5.40863,5.44375,5.38075,5.25022,5.37933,5.25459,5.37719,5.25169,5.73985,
     5.79174,5.79266,5.79588,5.7963,5.81949,5.82273,5.85658,5.85442,5.84779,5.77679,5.83323,5.84524,5.84439,5.84488,
     5.84744,5.84683,5.84524,5.84594,5.84656,5.84639,5.84461,5.84515,5.84206,5.8396,5.84497,5.84801,5.84608,5.84608,
     5.84783,5.84726,5.84844,5.8477,5.84796,5.84757,5.84822,5.84814,5.84617,5.83451,5.77658,5.84309,5.85496,5.85761,
     5.82555,5.82206,5.78982,5.78482,5.7778,5.78327,5.74898,5.25459,5.37503,5.25459,5.37283,5.25169,5.37862,5.44473,
     5.41041,5.34498,5.29551,5.25602,5.2283,5.17428,5.14504,5.09342,5.12256,4.98721,5.07106,5.02642,5.10031,5.11018,
     5.05447,5.10031,4.7426};

m_graphMap

GraphMap_t AthHistogramming::m_graphMap

privateinherited

The map of TGraph names to their pointers.

Definition at line 222 of file AthHistogramming.h.

m_grl

asg::AnaToolHandle<IGoodRunsListSelectionTool> ZdcNtuple::m_grl

Definition at line 101 of file ZdcNtuple.h.

m_hasBeginInputFile

bool EL::AnaAlgorithm::m_hasBeginInputFile {false}

privateinherited

the value of hasBeginInputFile

Definition at line 590 of file AnaAlgorithm.h.

{false};

m_hasEndInputFile

bool EL::AnaAlgorithm::m_hasEndInputFile {false}

privateinherited

the value of hasEndInputFile

Definition at line 594 of file AnaAlgorithm.h.

{false};

m_hasFileExecute

bool EL::AnaAlgorithm::m_hasFileExecute {false}

privateinherited

the value of hasFileExecute

Definition at line 586 of file AnaAlgorithm.h.

{false};

m_histMap

HistMap_t AthHistogramming::m_histMap

privateinherited

The map of histogram names to their pointers.

Definition at line 201 of file AthHistogramming.h.

m_histNamePostfix

std::string AthHistogramAlgorithm::m_histNamePostfix

privateinherited

The postfix for the histogram THx name.

Definition at line 97 of file AthHistogramAlgorithm.h.

m_histNamePrefix

std::string AthHistogramAlgorithm::m_histNamePrefix

privateinherited

The prefix for the histogram THx name.

Definition at line 94 of file AthHistogramAlgorithm.h.

m_histSvc

ServiceHandle<ITHistSvc> AthHistogramAlgorithm::m_histSvc

privateinherited

Default constructor: AthHistogramAlgorithm();.

a handle on the Hist/TTree registration service

Definition at line 83 of file AthHistogramAlgorithm.h.

m_histTitlePostfix

std::string AthHistogramAlgorithm::m_histTitlePostfix

privateinherited

The postfix for the histogram THx title.

Definition at line 103 of file AthHistogramAlgorithm.h.

m_histTitlePrefix

std::string AthHistogramAlgorithm::m_histTitlePrefix

privateinherited

The prefix for the histogram THx title.

Definition at line 100 of file AthHistogramAlgorithm.h.

m_injMapRunToken

ZdcInjPulserAmpMap::Token ZdcNtuple::m_injMapRunToken {}

private

Definition at line 50 of file ZdcNtuple.h.

{};

m_inputMetaStore

MetaStore_t EL::AnaAlgorithm::m_inputMetaStore

privateinherited

Object accessing the input metadata store.

Definition at line 551 of file AnaAlgorithm.h.

m_isMC

bool ZdcNtuple::m_isMC

Definition at line 133 of file ZdcNtuple.h.

m_lastRunNumber

unsigned int ZdcNtuple::m_lastRunNumber {0}

private

Definition at line 49 of file ZdcNtuple.h.

{0};

m_lvl1EnergySumRoI

const xAOD::EnergySumRoI* ZdcNtuple::m_lvl1EnergySumRoI

Definition at line 124 of file ZdcNtuple.h.

m_mbtsInfo

const xAOD::ForwardEventInfoContainer* ZdcNtuple::m_mbtsInfo

Definition at line 118 of file ZdcNtuple.h.

m_mbtsModules

const xAOD::MBTSModuleContainer* ZdcNtuple::m_mbtsModules

Definition at line 119 of file ZdcNtuple.h.

m_mcEventCollectionName

SG::ReadHandleKey<McEventCollection> ZdcNtuple::m_mcEventCollectionName {this, "MCEventCollectionName", "TruthEvent", ""}

Definition at line 109 of file ZdcNtuple.h.

{this, "MCEventCollectionName", "TruthEvent", ""};

m_msg

MsgStream AthHistogramming::m_msg

privateinherited

Cached Message Stream.

Definition at line 248 of file AthHistogramming.h.

m_name

std::string AthHistogramming::m_name

privateinherited

Instance name.

Definition at line 245 of file AthHistogramming.h.

m_nTriggers

int ZdcNtuple::m_nTriggers

Definition at line 131 of file ZdcNtuple.h.

m_outputMetaStore

MetaStore_t EL::AnaAlgorithm::m_outputMetaStore

privateinherited

Object accessing the output metadata store.

Definition at line 555 of file AnaAlgorithm.h.

m_outputTree

TTree* ZdcNtuple::m_outputTree

Definition at line 147 of file ZdcNtuple.h.

m_prefix

std::string AthHistogramAlgorithm::m_prefix

privateinherited

Name of the ROOT output stream (file)

Definition at line 88 of file AthHistogramAlgorithm.h.

m_primaryVertices

const xAOD::VertexContainer* ZdcNtuple::m_primaryVertices

Definition at line 121 of file ZdcNtuple.h.

m_rerunChainGroups

std::vector<const Trig::ChainGroup*> ZdcNtuple::m_rerunChainGroups

Definition at line 137 of file ZdcNtuple.h.

m_rootDir

std::string AthHistogramAlgorithm::m_rootDir

privateinherited

Name of the ROOT directory.

Definition at line 91 of file AthHistogramAlgorithm.h.

m_scaledownCounter

int ZdcNtuple::m_scaledownCounter

Definition at line 135 of file ZdcNtuple.h.

m_selTool

ToolHandle< InDet::IInDetTrackSelectionTool > ZdcNtuple::m_selTool

Definition at line 103 of file ZdcNtuple.h.

m_setupTrigHist

bool ZdcNtuple::m_setupTrigHist

Definition at line 134 of file ZdcNtuple.h.

m_streamName

std::string AthHistogramming::m_streamName

privateinherited

Name of the ROOT output stream (file)

Definition at line 226 of file AthHistogramming.h.

m_trackParticles

const xAOD::TrackParticleContainer* ZdcNtuple::m_trackParticles

Definition at line 123 of file ZdcNtuple.h.

m_treeMap

TreeMap_t AthHistogramming::m_treeMap

privateinherited

The map of TTree names to their pointers.

Definition at line 215 of file AthHistogramming.h.

m_trigDecision

const xAOD::TrigDecision* ZdcNtuple::m_trigDecision

Definition at line 115 of file ZdcNtuple.h.

m_trigDecisionTool

PublicToolHandle<Trig::TrigDecisionTool> ZdcNtuple::m_trigDecisionTool { this, "TrigDecisionTool", "", "Handle to the TrigDecisionTool" }

Definition at line 99 of file ZdcNtuple.h.

{ this, "TrigDecisionTool", "", "Handle to the TrigDecisionTool" };

m_trigT2MbtsBits

const xAOD::TrigT2MbtsBitsContainer* ZdcNtuple::m_trigT2MbtsBits

Definition at line 120 of file ZdcNtuple.h.

m_truthParticleContainer

const xAOD::TruthParticleContainer* ZdcNtuple::m_truthParticleContainer

Definition at line 125 of file ZdcNtuple.h.

m_TTcontainer

const xAOD::TriggerTowerContainer* ZdcNtuple::m_TTcontainer

Definition at line 126 of file ZdcNtuple.h.

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

m_zdcAnalysisTool

asg::AnaToolHandle<ZDC::IZdcAnalysisTool> ZdcNtuple::m_zdcAnalysisTool

Definition at line 102 of file ZdcNtuple.h.

m_zdcInjPulserAmpMap

std::shared_ptr<ZdcInjPulserAmpMap> ZdcNtuple::m_zdcInjPulserAmpMap

Definition at line 129 of file ZdcNtuple.h.

m_zdcModuleContainerName

SG::ReadHandleKey<xAOD::ZdcModuleContainer> ZdcNtuple::m_zdcModuleContainerName { this, "ZdcModuleContainerName", "ZdcModules", "" }

Definition at line 105 of file ZdcNtuple.h.

{ this, "ZdcModuleContainerName", "ZdcModules", "" };

m_zdcSumContainerName

SG::ReadHandleKey<xAOD::ZdcModuleContainer> ZdcNtuple::m_zdcSumContainerName { this, "ZdcSumContainerName", "ZdcSums", "" }

Definition at line 107 of file ZdcNtuple.h.

{ this, "ZdcSumContainerName", "ZdcSums", "" };

nProtons

int ZdcNtuple::nProtons

Definition at line 407 of file ZdcNtuple.h.

nsamplesZdc

size_t ZdcNtuple::nsamplesZdc

Definition at line 81 of file ZdcNtuple.h.

outputName

std::string ZdcNtuple::outputName

Definition at line 94 of file ZdcNtuple.h.

outputTreeScaledown

int ZdcNtuple::outputTreeScaledown

Definition at line 95 of file ZdcNtuple.h.

p_beam

TLorentzVector ZdcNtuple::p_beam

Definition at line 423 of file ZdcNtuple.h.

p_scat

TLorentzVector ZdcNtuple::p_scat

Definition at line 424 of file ZdcNtuple.h.

pbpb2023

bool ZdcNtuple::pbpb2023

Definition at line 85 of file ZdcNtuple.h.

proton_e

std::vector<double> ZdcNtuple::proton_e

Definition at line 411 of file ZdcNtuple.h.

proton_eLoss

std::vector<double> ZdcNtuple::proton_eLoss

Definition at line 413 of file ZdcNtuple.h.

proton_eta

std::vector<double> ZdcNtuple::proton_eta

Definition at line 409 of file ZdcNtuple.h.

proton_phi

std::vector<double> ZdcNtuple::proton_phi

Definition at line 410 of file ZdcNtuple.h.

proton_pt

std::vector<double> ZdcNtuple::proton_pt

Definition at line 408 of file ZdcNtuple.h.

proton_side

std::vector<int> ZdcNtuple::proton_side

Definition at line 412 of file ZdcNtuple.h.

proton_t

std::vector<double> ZdcNtuple::proton_t

Definition at line 414 of file ZdcNtuple.h.

proton_track_nClusters

std::vector<std::vector<int> > ZdcNtuple::proton_track_nClusters

Definition at line 421 of file ZdcNtuple.h.

proton_track_stationID

std::vector<std::vector<int> > ZdcNtuple::proton_track_stationID

Definition at line 415 of file ZdcNtuple.h.

proton_track_xLocal

std::vector<std::vector<float> > ZdcNtuple::proton_track_xLocal

Definition at line 416 of file ZdcNtuple.h.

proton_track_xSlope

std::vector<std::vector<float> > ZdcNtuple::proton_track_xSlope

Definition at line 419 of file ZdcNtuple.h.

proton_track_yLocal

std::vector<std::vector<float> > ZdcNtuple::proton_track_yLocal

Definition at line 417 of file ZdcNtuple.h.

proton_track_ySlope

std::vector<std::vector<float> > ZdcNtuple::proton_track_ySlope

Definition at line 420 of file ZdcNtuple.h.

proton_track_zLocal

std::vector<std::vector<float> > ZdcNtuple::proton_track_zLocal

Definition at line 418 of file ZdcNtuple.h.

reprocZdc

bool ZdcNtuple::reprocZdc

Definition at line 79 of file ZdcNtuple.h.

slimmed

bool ZdcNtuple::slimmed

Definition at line 53 of file ZdcNtuple.h.

t_actIntPerCrossing

float ZdcNtuple::t_actIntPerCrossing

Definition at line 159 of file ZdcNtuple.h.

t_avgIntPerCrossing

float ZdcNtuple::t_avgIntPerCrossing

Definition at line 158 of file ZdcNtuple.h.

t_bcid

uint32_t ZdcNtuple::t_bcid

Definition at line 151 of file ZdcNtuple.h.

t_bunchGroup

uint8_t ZdcNtuple::t_bunchGroup

Definition at line 153 of file ZdcNtuple.h.

t_cc_e

std::vector<float> ZdcNtuple::t_cc_e

Definition at line 391 of file ZdcNtuple.h.

t_cc_eta

std::vector<float> ZdcNtuple::t_cc_eta

Definition at line 389 of file ZdcNtuple.h.

t_cc_layer

std::vector<int> ZdcNtuple::t_cc_layer

Definition at line 393 of file ZdcNtuple.h.

t_cc_phi

std::vector<float> ZdcNtuple::t_cc_phi

Definition at line 390 of file ZdcNtuple.h.

t_cc_pt

std::vector<float> ZdcNtuple::t_cc_pt

Definition at line 388 of file ZdcNtuple.h.

t_cc_raw_e

std::vector<float> ZdcNtuple::t_cc_raw_e

Definition at line 397 of file ZdcNtuple.h.

t_cc_raw_eta

std::vector<float> ZdcNtuple::t_cc_raw_eta

Definition at line 395 of file ZdcNtuple.h.

t_cc_raw_m

std::vector<float> ZdcNtuple::t_cc_raw_m

Definition at line 394 of file ZdcNtuple.h.

t_cc_raw_phi

std::vector<float> ZdcNtuple::t_cc_raw_phi

Definition at line 396 of file ZdcNtuple.h.

t_cc_raw_samp

std::vector<std::vector<float> > ZdcNtuple::t_cc_raw_samp

Definition at line 398 of file ZdcNtuple.h.

t_cc_sig

std::vector<float> ZdcNtuple::t_cc_sig

Definition at line 392 of file ZdcNtuple.h.

t_centroidDecorationsAvailable

bool ZdcNtuple::t_centroidDecorationsAvailable

Definition at line 265 of file ZdcNtuple.h.

t_centroidEventValid

char ZdcNtuple::t_centroidEventValid

Definition at line 266 of file ZdcNtuple.h.

t_centroidStatus

unsigned int ZdcNtuple::t_centroidStatus[2]

Definition at line 267 of file ZdcNtuple.h.

t_clusEt

float ZdcNtuple::t_clusEt

Definition at line 400 of file ZdcNtuple.h.

t_clusetaMax

float ZdcNtuple::t_clusetaMax

Definition at line 402 of file ZdcNtuple.h.

t_clusEtMax

float ZdcNtuple::t_clusEtMax

Definition at line 401 of file ZdcNtuple.h.

t_clusphiMax

float ZdcNtuple::t_clusphiMax

Definition at line 403 of file ZdcNtuple.h.

t_cosDeltaReactionPlaneAngle

float ZdcNtuple::t_cosDeltaReactionPlaneAngle

Definition at line 279 of file ZdcNtuple.h.

t_decisions

bool ZdcNtuple::t_decisions[200]

Definition at line 171 of file ZdcNtuple.h.

t_edgeGapA

float ZdcNtuple::t_edgeGapA

Definition at line 323 of file ZdcNtuple.h.

t_edgeGapC

float ZdcNtuple::t_edgeGapC

Definition at line 324 of file ZdcNtuple.h.

t_eventNumber

uint32_t ZdcNtuple::t_eventNumber

Definition at line 149 of file ZdcNtuple.h.

t_extendedLevel1ID

uint32_t ZdcNtuple::t_extendedLevel1ID

Definition at line 155 of file ZdcNtuple.h.

t_fcalEt

float ZdcNtuple::t_fcalEt

Definition at line 309 of file ZdcNtuple.h.

t_fcalEtA

float ZdcNtuple::t_fcalEtA

Definition at line 310 of file ZdcNtuple.h.

t_fcalEtA_TT

float ZdcNtuple::t_fcalEtA_TT

Definition at line 312 of file ZdcNtuple.h.

t_fcalEtA_TTsum

float ZdcNtuple::t_fcalEtA_TTsum

Definition at line 314 of file ZdcNtuple.h.

t_fcalEtC

float ZdcNtuple::t_fcalEtC

Definition at line 311 of file ZdcNtuple.h.

t_fcalEtC_TT

float ZdcNtuple::t_fcalEtC_TT

Definition at line 313 of file ZdcNtuple.h.

t_fcalEtC_TTsum

float ZdcNtuple::t_fcalEtC_TTsum

Definition at line 315 of file ZdcNtuple.h.

t_L1ET

float ZdcNtuple::t_L1ET

Definition at line 139 of file ZdcNtuple.h.

t_L1ET24

float ZdcNtuple::t_L1ET24

Definition at line 140 of file ZdcNtuple.h.

t_lumiBlock

uint32_t ZdcNtuple::t_lumiBlock

Definition at line 150 of file ZdcNtuple.h.

t_mbts_countA

uint16_t ZdcNtuple::t_mbts_countA

Definition at line 339 of file ZdcNtuple.h.

t_mbts_countC

uint16_t ZdcNtuple::t_mbts_countC

Definition at line 340 of file ZdcNtuple.h.

t_mbts_in_e

float ZdcNtuple::t_mbts_in_e[2][8]

Definition at line 174 of file ZdcNtuple.h.

t_mbts_in_t

float ZdcNtuple::t_mbts_in_t[2][8]

Definition at line 176 of file ZdcNtuple.h.

t_mbts_out_e

float ZdcNtuple::t_mbts_out_e[2][4]

Definition at line 175 of file ZdcNtuple.h.

t_mbts_out_t

float ZdcNtuple::t_mbts_out_t[2][4]

Definition at line 177 of file ZdcNtuple.h.

t_mbts_timeA

float ZdcNtuple::t_mbts_timeA

Definition at line 341 of file ZdcNtuple.h.

t_mbts_timeC

float ZdcNtuple::t_mbts_timeC

Definition at line 342 of file ZdcNtuple.h.

t_mbts_timeDiff

float ZdcNtuple::t_mbts_timeDiff

Definition at line 343 of file ZdcNtuple.h.

t_nclus

uint32_t ZdcNtuple::t_nclus

Definition at line 387 of file ZdcNtuple.h.

t_nstrong

int ZdcNtuple::t_nstrong

Definition at line 292 of file ZdcNtuple.h.

t_ntrk

uint32_t ZdcNtuple::t_ntrk

Definition at line 359 of file ZdcNtuple.h.

t_nvtx

int ZdcNtuple::t_nvtx

Definition at line 297 of file ZdcNtuple.h.

t_nvx

int ZdcNtuple::t_nvx

Definition at line 281 of file ZdcNtuple.h.

t_passBits

uint32_t ZdcNtuple::t_passBits

Definition at line 154 of file ZdcNtuple.h.

t_prescales

float ZdcNtuple::t_prescales[200]

Definition at line 170 of file ZdcNtuple.h.

t_puvxntrk

int ZdcNtuple::t_puvxntrk

Definition at line 290 of file ZdcNtuple.h.

t_puvxsumpt

float ZdcNtuple::t_puvxsumpt

Definition at line 291 of file ZdcNtuple.h.

t_puvxz

float ZdcNtuple::t_puvxz

Definition at line 289 of file ZdcNtuple.h.

t_pvindex

int ZdcNtuple::t_pvindex

Definition at line 286 of file ZdcNtuple.h.

t_raw15

uint16_t ZdcNtuple::t_raw15[2][4][2][2][15]

Definition at line 349 of file ZdcNtuple.h.

t_raw24

uint16_t ZdcNtuple::t_raw24[2][4][2][2][24]

Definition at line 350 of file ZdcNtuple.h.

t_raw32

uint16_t ZdcNtuple::t_raw32[2][4][2][2][32]

Definition at line 351 of file ZdcNtuple.h.

t_raw40

uint16_t ZdcNtuple::t_raw40[2][4][2][2][40]

Definition at line 352 of file ZdcNtuple.h.

t_raw7

uint16_t ZdcNtuple::t_raw7[2][4][2][2][7]

Definition at line 348 of file ZdcNtuple.h.

t_reactionPlaneAngle

float ZdcNtuple::t_reactionPlaneAngle[2]

Definition at line 278 of file ZdcNtuple.h.

t_rerunDecisions

bool ZdcNtuple::t_rerunDecisions[200]

Definition at line 172 of file ZdcNtuple.h.

t_RpdChannelAmplitude

float ZdcNtuple::t_RpdChannelAmplitude[2][16]

Definition at line 255 of file ZdcNtuple.h.

t_RpdChannelAmplitudeCalib

float ZdcNtuple::t_RpdChannelAmplitudeCalib[2][16]

Definition at line 256 of file ZdcNtuple.h.

t_RpdChannelBaseline

float ZdcNtuple::t_RpdChannelBaseline[2][16]

Definition at line 248 of file ZdcNtuple.h.

t_RpdChannelMaxADC

float ZdcNtuple::t_RpdChannelMaxADC[2][16]

Definition at line 257 of file ZdcNtuple.h.

t_RpdChannelMaxADCCalib

float ZdcNtuple::t_RpdChannelMaxADCCalib[2][16]

Definition at line 258 of file ZdcNtuple.h.

t_RpdChannelMaxSample

unsigned int ZdcNtuple::t_RpdChannelMaxSample[2][16]

Definition at line 259 of file ZdcNtuple.h.

t_RpdChannelPileupExpFitMSE

float ZdcNtuple::t_RpdChannelPileupExpFitMSE[2][16]

Definition at line 253 of file ZdcNtuple.h.

t_RpdChannelPileupExpFitParamErrs

float ZdcNtuple::t_RpdChannelPileupExpFitParamErrs[2][16][2]

Definition at line 251 of file ZdcNtuple.h.

t_RpdChannelPileupExpFitParams

float ZdcNtuple::t_RpdChannelPileupExpFitParams[2][16][2]

Definition at line 249 of file ZdcNtuple.h.

t_RpdChannelPileupFrac

float ZdcNtuple::t_RpdChannelPileupFrac[2][16]

Definition at line 261 of file ZdcNtuple.h.

t_RpdChannelPileupStretchedExpFitMSE

float ZdcNtuple::t_RpdChannelPileupStretchedExpFitMSE[2][16]

Definition at line 254 of file ZdcNtuple.h.

t_RpdChannelPileupStretchedExpFitParamErrs

float ZdcNtuple::t_RpdChannelPileupStretchedExpFitParamErrs[2][16][3]

Definition at line 252 of file ZdcNtuple.h.

t_RpdChannelPileupStretchedExpFitParams

float ZdcNtuple::t_RpdChannelPileupStretchedExpFitParams[2][16][3]

Definition at line 250 of file ZdcNtuple.h.

t_RpdChannelStatus

unsigned int ZdcNtuple::t_RpdChannelStatus[2][16]

Definition at line 260 of file ZdcNtuple.h.

t_RPDChannelSubtrAmp

float ZdcNtuple::t_RPDChannelSubtrAmp[2][16]

Definition at line 268 of file ZdcNtuple.h.

t_rpdDecodingError

uint8_t ZdcNtuple::t_rpdDecodingError

Definition at line 164 of file ZdcNtuple.h.

t_RpdModuleTruthNphotons

unsigned int ZdcNtuple::t_RpdModuleTruthNphotons[2][16]

Definition at line 263 of file ZdcNtuple.h.

t_rpdRaw

uint16_t ZdcNtuple::t_rpdRaw[2][16][24]

Definition at line 354 of file ZdcNtuple.h.

t_rpdRaw32

uint16_t ZdcNtuple::t_rpdRaw32[2][16][32]

Definition at line 355 of file ZdcNtuple.h.

t_rpdRaw40

uint16_t ZdcNtuple::t_rpdRaw40[2][16][40]

Definition at line 356 of file ZdcNtuple.h.

t_RpdSideStatus

unsigned int ZdcNtuple::t_RpdSideStatus[2]

Definition at line 262 of file ZdcNtuple.h.

t_RPDSubtrAmpSum

float ZdcNtuple::t_RPDSubtrAmpSum[2]

Definition at line 269 of file ZdcNtuple.h.

t_runNumber

uint32_t ZdcNtuple::t_runNumber

Definition at line 148 of file ZdcNtuple.h.

t_T2mbts_countAin

uint16_t ZdcNtuple::t_T2mbts_countAin

Definition at line 345 of file ZdcNtuple.h.

t_T2mbts_countCin

uint16_t ZdcNtuple::t_T2mbts_countCin

Definition at line 346 of file ZdcNtuple.h.

t_T2mbts_in_e

float ZdcNtuple::t_T2mbts_in_e[2][8]

Definition at line 179 of file ZdcNtuple.h.

t_T2mbts_in_t

float ZdcNtuple::t_T2mbts_in_t[2][8]

Definition at line 181 of file ZdcNtuple.h.

t_T2mbts_out_e

float ZdcNtuple::t_T2mbts_out_e[2][4]

Definition at line 180 of file ZdcNtuple.h.

t_T2mbts_out_t

float ZdcNtuple::t_T2mbts_out_t[2][4]

Definition at line 182 of file ZdcNtuple.h.

t_tav

uint32_t ZdcNtuple::t_tav[16]

Definition at line 184 of file ZdcNtuple.h.

t_tbp

uint32_t ZdcNtuple::t_tbp[16]

Definition at line 185 of file ZdcNtuple.h.

t_timeStamp

uint32_t ZdcNtuple::t_timeStamp

Definition at line 156 of file ZdcNtuple.h.

t_timeStampNSOffset

uint32_t ZdcNtuple::t_timeStampNSOffset

Definition at line 157 of file ZdcNtuple.h.

t_totalEt

float ZdcNtuple::t_totalEt

Definition at line 317 of file ZdcNtuple.h.

t_totalEt24

float ZdcNtuple::t_totalEt24

Definition at line 319 of file ZdcNtuple.h.

t_totalEt24_TTsum

float ZdcNtuple::t_totalEt24_TTsum

Definition at line 320 of file ZdcNtuple.h.

t_totalEt_TTsum

float ZdcNtuple::t_totalEt_TTsum

Definition at line 318 of file ZdcNtuple.h.

t_trigger

uint64_t ZdcNtuple::t_trigger

Definition at line 166 of file ZdcNtuple.h.

t_trigger_TBP

uint32_t ZdcNtuple::t_trigger_TBP

Definition at line 167 of file ZdcNtuple.h.

t_trk_charge

std::vector<int8_t> ZdcNtuple::t_trk_charge

Definition at line 369 of file ZdcNtuple.h.

t_trk_d0

std::vector<float> ZdcNtuple::t_trk_d0

Definition at line 365 of file ZdcNtuple.h.

t_trk_e

std::vector<float> ZdcNtuple::t_trk_e

Definition at line 363 of file ZdcNtuple.h.

t_trk_eta

std::vector<float> ZdcNtuple::t_trk_eta

Definition at line 361 of file ZdcNtuple.h.

t_trk_exPixHits

std::vector<uint8_t> ZdcNtuple::t_trk_exPixHits

Definition at line 381 of file ZdcNtuple.h.

t_trk_index

std::vector<int> ZdcNtuple::t_trk_index

Definition at line 371 of file ZdcNtuple.h.

t_trk_inPixHits

std::vector<uint8_t> ZdcNtuple::t_trk_inPixHits

Definition at line 380 of file ZdcNtuple.h.

t_trk_nexPixHits

std::vector<uint8_t> ZdcNtuple::t_trk_nexPixHits

Definition at line 383 of file ZdcNtuple.h.

t_trk_ninPixHits

std::vector<uint8_t> ZdcNtuple::t_trk_ninPixHits

Definition at line 382 of file ZdcNtuple.h.

t_trk_nPixDead

std::vector<uint8_t> ZdcNtuple::t_trk_nPixDead

Definition at line 374 of file ZdcNtuple.h.

t_trk_nPixHits

std::vector<uint8_t> ZdcNtuple::t_trk_nPixHits

Definition at line 372 of file ZdcNtuple.h.

t_trk_nPixHoles

std::vector<uint8_t> ZdcNtuple::t_trk_nPixHoles

Definition at line 376 of file ZdcNtuple.h.

t_trk_nSctDead

std::vector<uint8_t> ZdcNtuple::t_trk_nSctDead

Definition at line 375 of file ZdcNtuple.h.

t_trk_nSctHits

std::vector<uint8_t> ZdcNtuple::t_trk_nSctHits

Definition at line 373 of file ZdcNtuple.h.

t_trk_nSctHoles

std::vector<uint8_t> ZdcNtuple::t_trk_nSctHoles

Definition at line 377 of file ZdcNtuple.h.

t_trk_nTrtHits

std::vector<uint8_t> ZdcNtuple::t_trk_nTrtHits

Definition at line 378 of file ZdcNtuple.h.

t_trk_nTrtOutliers

std::vector<uint8_t> ZdcNtuple::t_trk_nTrtOutliers

Definition at line 379 of file ZdcNtuple.h.

t_trk_phi

std::vector<float> ZdcNtuple::t_trk_phi

Definition at line 362 of file ZdcNtuple.h.

t_trk_pixeldEdx

std::vector<float> ZdcNtuple::t_trk_pixeldEdx

Definition at line 384 of file ZdcNtuple.h.

t_trk_pt

std::vector<float> ZdcNtuple::t_trk_pt

Definition at line 360 of file ZdcNtuple.h.

t_trk_quality

std::vector<int16_t> ZdcNtuple::t_trk_quality

Definition at line 370 of file ZdcNtuple.h.

t_trk_theta

std::vector<float> ZdcNtuple::t_trk_theta

Definition at line 364 of file ZdcNtuple.h.

t_trk_vtxz

std::vector<float> ZdcNtuple::t_trk_vtxz

Definition at line 368 of file ZdcNtuple.h.

t_trk_vz

std::vector<float> ZdcNtuple::t_trk_vz

Definition at line 367 of file ZdcNtuple.h.

t_trk_z0

std::vector<float> ZdcNtuple::t_trk_z0

Definition at line 366 of file ZdcNtuple.h.

t_vInj

float ZdcNtuple::t_vInj

Definition at line 152 of file ZdcNtuple.h.

t_vtx_ntrk

std::vector<int16_t> ZdcNtuple::t_vtx_ntrk

Definition at line 304 of file ZdcNtuple.h.

t_vtx_ntrk_all

std::vector<int16_t> ZdcNtuple::t_vtx_ntrk_all

Definition at line 302 of file ZdcNtuple.h.

t_vtx_sumpt2

std::vector<float> ZdcNtuple::t_vtx_sumpt2

Definition at line 305 of file ZdcNtuple.h.

t_vtx_sumpt2_all

std::vector<float> ZdcNtuple::t_vtx_sumpt2_all

Definition at line 303 of file ZdcNtuple.h.

t_vtx_trk_index

std::vector< std::vector<int16_t> > ZdcNtuple::t_vtx_trk_index

Definition at line 306 of file ZdcNtuple.h.

t_vtx_type

std::vector<int8_t> ZdcNtuple::t_vtx_type

Definition at line 298 of file ZdcNtuple.h.

t_vtx_x

std::vector<float> ZdcNtuple::t_vtx_x

Definition at line 299 of file ZdcNtuple.h.

t_vtx_y

std::vector<float> ZdcNtuple::t_vtx_y

Definition at line 300 of file ZdcNtuple.h.

t_vtx_z

std::vector<float> ZdcNtuple::t_vtx_z

Definition at line 301 of file ZdcNtuple.h.

t_vx

float ZdcNtuple::t_vx[3]

Definition at line 282 of file ZdcNtuple.h.

t_vx_trk_index

std::vector<int> ZdcNtuple::t_vx_trk_index

Definition at line 284 of file ZdcNtuple.h.

t_vxcov

float ZdcNtuple::t_vxcov[6]

Definition at line 287 of file ZdcNtuple.h.

t_vxngoodmuon

int ZdcNtuple::t_vxngoodmuon

Definition at line 295 of file ZdcNtuple.h.

t_vxnlooseprimary

int ZdcNtuple::t_vxnlooseprimary

Definition at line 293 of file ZdcNtuple.h.

t_vxnminbias

int ZdcNtuple::t_vxnminbias

Definition at line 294 of file ZdcNtuple.h.

t_vxntrk

int ZdcNtuple::t_vxntrk

Definition at line 283 of file ZdcNtuple.h.

t_vxsumpt2

float ZdcNtuple::t_vxsumpt2

Definition at line 288 of file ZdcNtuple.h.

t_vxtype

int ZdcNtuple::t_vxtype

Definition at line 285 of file ZdcNtuple.h.

t_xCentroid

float ZdcNtuple::t_xCentroid[2]

Definition at line 274 of file ZdcNtuple.h.

t_xCentroidPreAvgSubtr

float ZdcNtuple::t_xCentroidPreAvgSubtr[2]

Definition at line 272 of file ZdcNtuple.h.

t_xCentroidPreGeomCorPreAvgSubtr

float ZdcNtuple::t_xCentroidPreGeomCorPreAvgSubtr[2]

Definition at line 270 of file ZdcNtuple.h.

t_xRowCentroid

float ZdcNtuple::t_xRowCentroid[2][4]

Definition at line 276 of file ZdcNtuple.h.

t_yCentroid

float ZdcNtuple::t_yCentroid[2]

Definition at line 275 of file ZdcNtuple.h.

t_yCentroidPreAvgSubtr

float ZdcNtuple::t_yCentroidPreAvgSubtr[2]

Definition at line 273 of file ZdcNtuple.h.

t_yCentroidPreGeomCorPreAvgSubtr

float ZdcNtuple::t_yCentroidPreGeomCorPreAvgSubtr[2]

Definition at line 271 of file ZdcNtuple.h.

t_yColCentroid

float ZdcNtuple::t_yColCentroid[2][4]

Definition at line 277 of file ZdcNtuple.h.

t_ZdcAmp

float ZdcNtuple::t_ZdcAmp[2]

Definition at line 187 of file ZdcNtuple.h.

t_ZdcAmpErr

float ZdcNtuple::t_ZdcAmpErr[2]

Definition at line 188 of file ZdcNtuple.h.

t_zdcDecodingError

uint8_t ZdcNtuple::t_zdcDecodingError

Definition at line 163 of file ZdcNtuple.h.

t_ZdcEnergy

float ZdcNtuple::t_ZdcEnergy[2]

Definition at line 189 of file ZdcNtuple.h.

t_ZdcEnergyErr

float ZdcNtuple::t_ZdcEnergyErr[2]

Definition at line 190 of file ZdcNtuple.h.

t_zdcEventInfoError

uint8_t ZdcNtuple::t_zdcEventInfoError

Definition at line 160 of file ZdcNtuple.h.

t_zdcEventInfoErrorWord

uint32_t ZdcNtuple::t_zdcEventInfoErrorWord

Definition at line 161 of file ZdcNtuple.h.

t_ZdcLucrodTriggerAmp

unsigned short ZdcNtuple::t_ZdcLucrodTriggerAmp[2][4]

Definition at line 228 of file ZdcNtuple.h.

t_ZdcLucrodTriggerAmpLG

unsigned short ZdcNtuple::t_ZdcLucrodTriggerAmpLG[2][4]

Definition at line 229 of file ZdcNtuple.h.

t_ZdcLucrodTriggerSideAmp

unsigned short ZdcNtuple::t_ZdcLucrodTriggerSideAmp[2]

Definition at line 197 of file ZdcNtuple.h.

t_ZdcLucrodTriggerSideAmpLG

unsigned short ZdcNtuple::t_ZdcLucrodTriggerSideAmpLG[2]

Definition at line 198 of file ZdcNtuple.h.

t_ZdcModuleAmp

float ZdcNtuple::t_ZdcModuleAmp[2][4]

Definition at line 215 of file ZdcNtuple.h.

t_ZdcModuleAmpCorrLGRefit

float ZdcNtuple::t_ZdcModuleAmpCorrLGRefit[2][4]

Definition at line 236 of file ZdcNtuple.h.

t_ZdcModuleAmpError

float ZdcNtuple::t_ZdcModuleAmpError[2][4]

Definition at line 223 of file ZdcNtuple.h.

t_ZdcModuleAmpLGRefit

float ZdcNtuple::t_ZdcModuleAmpLGRefit[2][4]

Definition at line 235 of file ZdcNtuple.h.

t_ZdcModuleBkgdMaxFraction

float ZdcNtuple::t_ZdcModuleBkgdMaxFraction[2][4]

Definition at line 224 of file ZdcNtuple.h.

t_ZdcModuleCalibAmp

float ZdcNtuple::t_ZdcModuleCalibAmp[2][4]

Definition at line 221 of file ZdcNtuple.h.

t_ZdcModuleCalibTime

float ZdcNtuple::t_ZdcModuleCalibTime[2][4]

Definition at line 222 of file ZdcNtuple.h.

t_ZdcModuleChisq

float ZdcNtuple::t_ZdcModuleChisq[2][4]

Definition at line 219 of file ZdcNtuple.h.

t_ZdcModuleChisqLGRefit

float ZdcNtuple::t_ZdcModuleChisqLGRefit[2][4]

Definition at line 239 of file ZdcNtuple.h.

t_ZdcModuleFitAmp

float ZdcNtuple::t_ZdcModuleFitAmp[2][4]

Definition at line 217 of file ZdcNtuple.h.

t_ZdcModuleFitT0

float ZdcNtuple::t_ZdcModuleFitT0[2][4]

Definition at line 218 of file ZdcNtuple.h.

t_ZdcModuleMask

unsigned int ZdcNtuple::t_ZdcModuleMask

Definition at line 195 of file ZdcNtuple.h.

t_ZdcModuleMaxADC

float ZdcNtuple::t_ZdcModuleMaxADC[2][4]

Definition at line 230 of file ZdcNtuple.h.

t_ZdcModuleMaxADCHG

float ZdcNtuple::t_ZdcModuleMaxADCHG[2][4]

Definition at line 231 of file ZdcNtuple.h.

t_ZdcModuleMaxADCLG

float ZdcNtuple::t_ZdcModuleMaxADCLG[2][4]

Definition at line 232 of file ZdcNtuple.h.

t_ZdcModuleMinDeriv2nd

float ZdcNtuple::t_ZdcModuleMinDeriv2nd[2][4]

Definition at line 225 of file ZdcNtuple.h.

t_ZdcModulePeakADCHG

float ZdcNtuple::t_ZdcModulePeakADCHG[2][4]

Definition at line 233 of file ZdcNtuple.h.

t_ZdcModulePeakADCLG

float ZdcNtuple::t_ZdcModulePeakADCLG[2][4]

Definition at line 234 of file ZdcNtuple.h.

t_ZdcModulePresample

float ZdcNtuple::t_ZdcModulePresample[2][4]

Definition at line 226 of file ZdcNtuple.h.

t_ZdcModulePreSampleAmp

float ZdcNtuple::t_ZdcModulePreSampleAmp[2][4]

Definition at line 227 of file ZdcNtuple.h.

t_ZdcModuleStatus

unsigned int ZdcNtuple::t_ZdcModuleStatus[2][4]

Definition at line 220 of file ZdcNtuple.h.

t_ZdcModuleT0LGRefit

float ZdcNtuple::t_ZdcModuleT0LGRefit[2][4]

Definition at line 237 of file ZdcNtuple.h.

t_ZdcModuleT0SubLGRefit

float ZdcNtuple::t_ZdcModuleT0SubLGRefit[2][4]

Definition at line 238 of file ZdcNtuple.h.

t_ZdcModuleTime

float ZdcNtuple::t_ZdcModuleTime[2][4]

Definition at line 216 of file ZdcNtuple.h.

t_ZdcModuleTruthEM

float ZdcNtuple::t_ZdcModuleTruthEM[2][7]

Definition at line 243 of file ZdcNtuple.h.

t_ZdcModuleTruthEscaped

float ZdcNtuple::t_ZdcModuleTruthEscaped[2][7]

Definition at line 245 of file ZdcNtuple.h.

t_ZdcModuleTruthInvis

float ZdcNtuple::t_ZdcModuleTruthInvis[2][7]

Definition at line 242 of file ZdcNtuple.h.

t_ZdcModuleTruthNonEM

float ZdcNtuple::t_ZdcModuleTruthNonEM[2][7]

Definition at line 244 of file ZdcNtuple.h.

t_ZdcModuleTruthNphotons

unsigned int ZdcNtuple::t_ZdcModuleTruthNphotons[2][7]

Definition at line 246 of file ZdcNtuple.h.

t_ZdcModuleTruthTotal

float ZdcNtuple::t_ZdcModuleTruthTotal[2][7]

Definition at line 241 of file ZdcNtuple.h.

t_ZdcNLEnergy

float ZdcNtuple::t_ZdcNLEnergy[2]

Definition at line 191 of file ZdcNtuple.h.

t_ZdcNLEnergyErr

float ZdcNtuple::t_ZdcNLEnergyErr[2]

Definition at line 192 of file ZdcNtuple.h.

t_ZdcStatus

short ZdcNtuple::t_ZdcStatus[2]

Definition at line 194 of file ZdcNtuple.h.

t_ZdcTime

float ZdcNtuple::t_ZdcTime[2]

Definition at line 193 of file ZdcNtuple.h.

t_ZdcTrigEff

float ZdcNtuple::t_ZdcTrigEff[2]

Definition at line 196 of file ZdcNtuple.h.

t_ZdcTruthEM

float ZdcNtuple::t_ZdcTruthEM[2]

Definition at line 201 of file ZdcNtuple.h.

t_ZdcTruthEscaped

float ZdcNtuple::t_ZdcTruthEscaped[2]

Definition at line 203 of file ZdcNtuple.h.

t_ZdcTruthInvis

float ZdcNtuple::t_ZdcTruthInvis[2]

Definition at line 200 of file ZdcNtuple.h.

t_ZdcTruthNonEM

float ZdcNtuple::t_ZdcTruthNonEM[2]

Definition at line 202 of file ZdcNtuple.h.

t_ZdcTruthParticleEnergy

std::vector< float > ZdcNtuple::t_ZdcTruthParticleEnergy

Definition at line 211 of file ZdcNtuple.h.

t_ZdcTruthParticlePid

std::vector< int > ZdcNtuple::t_ZdcTruthParticlePid

Definition at line 212 of file ZdcNtuple.h.

t_ZdcTruthParticlePosx

std::vector< float > ZdcNtuple::t_ZdcTruthParticlePosx

Definition at line 204 of file ZdcNtuple.h.

t_ZdcTruthParticlePosy

std::vector< float > ZdcNtuple::t_ZdcTruthParticlePosy

Definition at line 205 of file ZdcNtuple.h.

t_ZdcTruthParticlePosz

std::vector< float > ZdcNtuple::t_ZdcTruthParticlePosz

Definition at line 206 of file ZdcNtuple.h.

t_ZdcTruthParticlePx

std::vector< float > ZdcNtuple::t_ZdcTruthParticlePx

Definition at line 208 of file ZdcNtuple.h.

t_ZdcTruthParticlePy

std::vector< float > ZdcNtuple::t_ZdcTruthParticlePy

Definition at line 209 of file ZdcNtuple.h.

t_ZdcTruthParticlePz

std::vector< float > ZdcNtuple::t_ZdcTruthParticlePz

Definition at line 210 of file ZdcNtuple.h.

t_ZdcTruthParticleStatus

std::vector< int > ZdcNtuple::t_ZdcTruthParticleStatus

Definition at line 213 of file ZdcNtuple.h.

t_ZdcTruthParticleTime

std::vector< float > ZdcNtuple::t_ZdcTruthParticleTime

Definition at line 207 of file ZdcNtuple.h.

t_ZdcTruthTotal

float ZdcNtuple::t_ZdcTruthTotal[2]

Definition at line 199 of file ZdcNtuple.h.

trackLimit

size_t ZdcNtuple::trackLimit

Definition at line 76 of file ZdcNtuple.h.

trackLimitReject

bool ZdcNtuple::trackLimitReject

Definition at line 77 of file ZdcNtuple.h.

useGRL

bool ZdcNtuple::useGRL

Definition at line 54 of file ZdcNtuple.h.

writeOnlyTriggers

bool ZdcNtuple::writeOnlyTriggers

Definition at line 59 of file ZdcNtuple.h.

zdcCalib

bool ZdcNtuple::zdcCalib

Definition at line 71 of file ZdcNtuple.h.

zdcConfig

std::string ZdcNtuple::zdcConfig

Definition at line 92 of file ZdcNtuple.h.

zdcInj

bool ZdcNtuple::zdcInj

Definition at line 73 of file ZdcNtuple.h.

zdcLaser

bool ZdcNtuple::zdcLaser

Definition at line 72 of file ZdcNtuple.h.

zdcLowGainMode

unsigned int ZdcNtuple::zdcLowGainMode

Definition at line 75 of file ZdcNtuple.h.

zdcOnly

bool ZdcNtuple::zdcOnly

Definition at line 74 of file ZdcNtuple.h.

---

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

• ZdcNtuple.h
• ZdcNtuple.cxx