The ECL Data Analysis Module. More...

#include <ECLDataAnalysisModule.h>

Inheritance diagram for ECLDataAnalysisModule:

Public Types
enum	EModulePropFlags { c_Input = 1 , c_Output = 2 , c_ParallelProcessingCertified = 4 , c_HistogramManager = 8 , c_InternalSerializer = 16 , c_TerminateInAllProcesses = 32 , c_DontCollectStatistics = 64 }
	Each module can be tagged with property flags, which indicate certain features of the module. More...

typedef ModuleCondition::EAfterConditionPath	EAfterConditionPath
	Forward the EAfterConditionPath definition from the ModuleCondition.

Public Member Functions
	ECLDataAnalysisModule ()
	Constructor of the module.

virtual	~ECLDataAnalysisModule () override
	Destructor of the module.

virtual void	initialize () override
	Initializes the Module.

virtual void	event () override
	event

virtual void	terminate () override
	Termination action.

virtual std::vector< std::string >	getFileNames (bool outputFiles)
	Return a list of output filenames for this modules.

virtual void	beginRun ()
	Called when entering a new run.

virtual void	endRun ()
	This method is called if the current run ends.

const std::string &	getName () const
	Returns the name of the module.

const std::string &	getType () const
	Returns the type of the module (i.e.

const std::string &	getPackage () const
	Returns the package this module is in.

const std::string &	getDescription () const
	Returns the description of the module.

void	setName (const std::string &name)
	Set the name of the module.

void	setPropertyFlags (unsigned int propertyFlags)
	Sets the flags for the module properties.

LogConfig &	getLogConfig ()
	Returns the log system configuration.

void	setLogConfig (const LogConfig &logConfig)
	Set the log system configuration.

void	setLogLevel (int logLevel)
	Configure the log level.

void	setDebugLevel (int debugLevel)
	Configure the debug messaging level.

void	setAbortLevel (int abortLevel)
	Configure the abort log level.

void	setLogInfo (int logLevel, unsigned int logInfo)
	Configure the printed log information for the given level.

void	if_value (const std::string &expression, const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
	Add a condition to the module.

void	if_false (const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
	A simplified version to add a condition to the module.

void	if_true (const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
	A simplified version to set the condition of the module.

bool	hasCondition () const
	Returns true if at least one condition was set for the module.

const ModuleCondition *	getCondition () const
	Return a pointer to the first condition (or nullptr, if none was set)

const std::vector< ModuleCondition > &	getAllConditions () const
	Return all set conditions for this module.

bool	evalCondition () const
	If at least one condition was set, it is evaluated and true returned if at least one condition returns true.

std::shared_ptr< Path >	getConditionPath () const
	Returns the path of the last true condition (if there is at least one, else reaturn a null pointer).

Module::EAfterConditionPath	getAfterConditionPath () const
	What to do after the conditional path is finished.

std::vector< std::shared_ptr< Path > >	getAllConditionPaths () const
	Return all condition paths currently set (no matter if the condition is true or not).

bool	hasProperties (unsigned int propertyFlags) const
	Returns true if all specified property flags are available in this module.

bool	hasUnsetForcedParams () const
	Returns true and prints error message if the module has unset parameters which the user has to set in the steering file.

const ModuleParamList &	getParamList () const
	Return module param list.

template<typename T>
ModuleParam< T > &	getParam (const std::string &name) const
	Returns a reference to a parameter.

bool	hasReturnValue () const
	Return true if this module has a valid return value set.

int	getReturnValue () const
	Return the return value set by this module.

std::shared_ptr< PathElement >	clone () const override
	Create an independent copy of this module.

std::shared_ptr< boost::python::list >	getParamInfoListPython () const
	Returns a python list of all parameters.

Static Public Member Functions
static void	exposePythonAPI ()
	Exposes methods of the Module class to Python.

Protected Member Functions
virtual void	def_initialize ()
	Wrappers to make the methods without "def_" prefix callable from Python.

virtual void	def_beginRun ()
	Wrapper method for the virtual function beginRun() that has the implementation to be used in a call from Python.

virtual void	def_event ()
	Wrapper method for the virtual function event() that has the implementation to be used in a call from Python.

virtual void	def_endRun ()
	This method can receive that the current run ends as a call from the Python side.

virtual void	def_terminate ()
	Wrapper method for the virtual function terminate() that has the implementation to be used in a call from Python.

void	setDescription (const std::string &description)
	Sets the description of the module.

void	setType (const std::string &type)
	Set the module type.

template<typename T>
void	addParam (const std::string &name, T &paramVariable, const std::string &description, const T &defaultValue)
	Adds a new parameter to the module.

template<typename T>
void	addParam (const std::string &name, T &paramVariable, const std::string &description)
	Adds a new enforced parameter to the module.

void	setReturnValue (int value)
	Sets the return value for this module as integer.

void	setReturnValue (bool value)
	Sets the return value for this module as bool.

void	setParamList (const ModuleParamList &params)
	Replace existing parameter list.

Private Member Functions
virtual const char *	eclSimHitArrayName () const
	Default name ECLCalDigits array.

virtual const char *	eclHitArrayName () const
	Default name ECLCalDigits array.

virtual const char *	eclDigitArrayName () const
	Default name ECLDigits array.

virtual const char *	eclCalDigitArrayName () const
	Default name ECLCalDigits array.

virtual const char *	eclShowerArrayName () const
	Default name ECLShower array.

virtual const char *	eclClusterArrayName () const
	Default name ECLClusters array.

virtual const char *	eclConnectedRegionArrayName () const
	Default name ECLConnectedRegions array.

virtual const char *	eclLocalMaximumArrayName () const
	Default name ECLLocalMaxima array.

virtual const char *	eclPureDigitArrayName () const
	Default name ECLPureDigits array.

virtual const char *	eclPureCalDigitArrayName () const
	Default name ECLPureCalDigits array.

virtual const char *	eclPureShowerArrayName () const
	Default name ECLPureShower array.

virtual const char *	eclPureClusterArrayName () const
	Default name ECLPureClusters array.

virtual const char *	eclPureConnectedRegionArrayName () const
	Default name ECLPureConnectedRegions array.

virtual const char *	eclPureLocalMaximumArrayName () const
	Default name ECLPureLocalMaxima array.

std::list< ModulePtr >	getModules () const override
	no submodules, return empty list

std::string	getPathString () const override
	return the module name.

void	setParamPython (const std::string &name, const boost::python::object &pyObj)
	Implements a method for setting boost::python objects.

void	setParamPythonDict (const boost::python::dict &dictionary)
	Implements a method for reading the parameter values from a boost::python dictionary.

Private Attributes
TFile *	m_rootFilePtr
	members of ECLReconstructor Module

std::string	m_rootFileName
	name of the root file

bool	m_writeToRoot
	if true, a rootFile named by m_rootFileName will be filled with info

bool	m_doTracking
	if true, info on tracking will be stored, job will fail if doTracking==1 and the tracking modules are not enabled at python level

bool	m_doPureCsI
	if true, info on pureCsI upgrade is stored

bool	m_doHits
	if true, info on Hits and SimHits is stored

bool	m_doDigits
	if true, info on Digits and CalDigits is stored

StoreArray< Track >	m_tracks
	Tracks storeArray.

StoreArray< TrackFitResult >	m_trackFitResults
	TrackFitResult storeArray.

StoreArray< ECLPidLikelihood >	m_eclPidLikelihoods
	ECLPidLikelihood storeArray.

StoreArray< ECLSimHit >	m_eclSimHits
	Store array: ECLSimHit.

StoreArray< ECLHit >	m_eclHits
	Store array: ECLHit.

StoreArray< ECLDigit >	m_eclDigits
	Store array: ECLDigit.

StoreArray< ECLCalDigit >	m_eclCalDigits
	Store array: ECLCalDigit.

StoreArray< ECLConnectedRegion >	m_eclConnectedRegions
	Store array: ECLConnectedRegion.

StoreArray< ECLShower >	m_eclShowers
	Store array: ECLShower.

StoreArray< ECLCluster >	m_eclClusters
	Store array: ECLCluster.

StoreArray< ECLLocalMaximum >	m_eclLocalMaximums
	Store array: ECLLocalMaximum.

StoreObjPtr< EventLevelClusteringInfo >	m_eventLevelClusteringInfo
	Store object pointer: EventLevelClusteringInfo.

StoreObjPtr< EventMetaData >	m_eventmetadata
	Store object pointer: EventMetaData.

StoreArray< ECLDigit >	m_eclPureDigits
	Store array: ECLPureDigit.

StoreArray< ECLCalDigit >	m_eclPureCalDigits
	Store array: ECLPureCalDigit.

StoreArray< ECLConnectedRegion >	m_eclPureConnectedRegions
	Store array: ECLPureConnectedRegion.

StoreArray< ECLShower >	m_eclPureShowers
	Store array: ECLPureShower.

StoreArray< ECLCluster >	m_eclPureClusters
	Store array: ECLPureCluster.

StoreArray< ECLLocalMaximum >	m_eclPureLocalMaximums
	Store array: ECLPureLocalMaximum.

StoreArray< MCParticle >	m_mcParticles
	MCParticles StoreArray.

TTree *	m_tree
	Root tree and file for saving the output.

int	m_iExperiment
	Experiment number.

int	m_iRun
	Run number.

int	m_iEvent
	Event number.

uint16_t	m_nECLCalDigitsOutOfTimeFWD {0}
	Number of out of time, energetic ECLCalDigits, FWD.

uint16_t	m_nECLCalDigitsOutOfTimeBarrel {0}
	Number of out of time, energetic ECLCalDigits, Barrel.

uint16_t	m_nECLCalDigitsOutOfTimeBWD {0}
	Number of out of time, energetic ECLCalDigits, BWD.

uint8_t	m_nECLShowersRejectedFWD {0}
	Number of photon showers that are rejected before storing to mdst (max.

uint8_t	m_nECLShowersRejectedBarrel {0}
	Number of photon showers that are rejected before storing to mdst (max.

uint8_t	m_nECLShowersRejectedBWD {0}
	Number of photon showers that are rejected before storing to mdst (max.

int	m_eclDigitMultip
	Number of ECLDigits per event.

std::vector< int > *	m_eclDigitIdx
	ECLDigit index.

std::vector< int > *	m_eclDigitToMC
	Index of MCParticle related to that ECLDigit.

std::vector< int > *	m_eclDigitCellId
	Number of ECLDigit CellId.

std::vector< int > *	m_eclDigitAmp
	ECLDigit amplitude.

std::vector< int > *	m_eclDigitTimeFit
	ECLDigit timing.

std::vector< int > *	m_eclDigitFitQuality
	ECLDigit fit quality.

std::vector< int > *	m_eclDigitToCalDigit
	Index of CalDigit related to that ECLDigit.

int	m_eclCalDigitMultip
	Number of ECLCalDigits per event.

std::vector< int > *	m_eclCalDigitIdx
	ECLCalDigit index.

std::vector< int > *	m_eclCalDigitToMC1
	Index of first MCParticle related to ECLCalDigit.

std::vector< int > *	m_eclCalDigitToMC1PDG
	PDG code of first MCParticle related to ECLCalDigit.

std::vector< double > *	m_eclCalDigitToMCWeight1
	Energy contribution of first MCParticle related to ECLCalDigit.

std::vector< int > *	m_eclCalDigitToMC2
	Index of second MCParticle related to ECLCalDigit.

std::vector< int > *	m_eclCalDigitToMC2PDG
	PDG code of second MCParticle related to ECLCalDigit.

std::vector< double > *	m_eclCalDigitToMCWeight2
	Energy contribution of second MCParticle related to ECLCalDigit.

std::vector< int > *	m_eclCalDigitToMC3
	Index of third MCParticle related to ECLCalDigit.

std::vector< int > *	m_eclCalDigitToMC3PDG
	PDG code of third MCParticle related to ECLCalDigit.

std::vector< double > *	m_eclCalDigitToMCWeight3
	Energy contribution of third MCParticle related to ECLCalDigit.

std::vector< int > *	m_eclCalDigitToMC4
	Index of fourth MCParticle related to ECLCalDigit.

std::vector< int > *	m_eclCalDigitToMC4PDG
	PDG code of fourth MCParticle related to ECLCalDigit.

std::vector< double > *	m_eclCalDigitToMCWeight4
	Energy contribution of fourth MCParticle related to ECLCalDigit.

std::vector< int > *	m_eclCalDigitToMC5
	Index of fifth MCParticle related to ECLCalDigit.

std::vector< int > *	m_eclCalDigitToMC5PDG
	PDG code of fifth MCParticle related to ECLCalDigit.

std::vector< double > *	m_eclCalDigitToMCWeight5
	Energy contribution of fifth MCParticle related to ECLCalDigit.

std::vector< double > *	m_eclCalDigitToBkgWeight
	Remaining energy contribution not associated to first five MCParticles related to ECLCalDigit.

std::vector< double > *	m_eclCalDigitSimHitSum
	Full energy contribution related to ECLCalDigit.

std::vector< int > *	m_eclCalDigitToShower
	Index of ECLShower related to that ECLCalDigit.

std::vector< int > *	m_eclCalDigitCellId
	Number of ECLCalDigit CellId.

std::vector< double > *	m_eclCalDigitAmp
	ECLCalDigit amplitude.

std::vector< double > *	m_eclCalDigitTimeFit
	ECLCalDigit timing.

std::vector< int > *	m_eclCalDigitFitQuality
	ECLCalDigit fit quality.

std::vector< int > *	m_eclCalDigitToCR
	Index of CR related to that ECLCalDigit.

std::vector< int > *	m_eclCalDigitToLM
	Index of LM related to that ECLCalDigit.

std::vector< int > *	m_eclCRIdx
	Connected Region ID.

std::vector< int > *	m_eclCRIsTrack
	Int for Connected Region - Track Match.

std::vector< double > *	m_eclCRLikelihoodMIPNGamma
	Connected Region MIP Likelihood.

std::vector< double > *	m_eclCRLikelihoodChargedHadron
	Connected Region Charged Hadron Likelihood.

std::vector< double > *	m_eclCRLikelihoodElectronNGamma
	Connected Region Electron Likelihood.

std::vector< double > *	m_eclCRLikelihoodNGamma
	Connected Region Gamma Likelihood.

std::vector< double > *	m_eclCRLikelihoodNeutralHadron
	Connected Region Neutral Hadron Likelihood.

std::vector< double > *	m_eclCRLikelihoodMergedPi0
	Connected Region Merged Pi0 Likelihood.

int	m_eclLMMultip
	Local Maxima multiplicity.

std::vector< int > *	m_eclLMId
	Local Maximum ID.

std::vector< int > *	m_eclLMType
	Local Maximum type.

std::vector< int > *	m_eclLMCellId
	Local Maximum Cell ID.

int	m_eclSimHitMultip
	Number of ECLSimHits per event.

std::vector< int > *	m_eclSimHitIdx
	Index of ECLSimHit.

std::vector< int > *	m_eclSimHitToMC
	Index of MCParticle related to that ECLSimHit.

std::vector< int > *	m_eclSimHitCellId
	ECLSimHit CellId.

std::vector< int > *	m_eclSimHitPdg
	PDG code of MCParticle associated to that ECLDigit.

std::vector< double > *	m_eclSimHitEnergyDep
	Energy deposition of ECLSimHit.

std::vector< double > *	m_eclSimHitFlightTime
	ECLSimhit Flight Time.

std::vector< double > *	m_eclSimHitX
	ECLSimHit X position.

std::vector< double > *	m_eclSimHitY
	ECLSimHit Y position.

std::vector< double > *	m_eclSimHitZ
	ECLSimHit Z position.

std::vector< double > *	m_eclSimHitPx
	ECLSimHit PX.

std::vector< double > *	m_eclSimHitPy
	ECLSimHit PY.

std::vector< double > *	m_eclSimHitPz
	ECLSimHit PZ.

int	m_eclHitMultip
	Number of ECLHits per event.

std::vector< int > *	m_eclHitIdx
	Index of ECLHits.

std::vector< int > *	m_eclHitToMC
	Index of MCParticle related to ECLHit.

std::vector< int > *	m_eclHitToDigit
	Index of ECLDigit related to ECLHit.

std::vector< int > *	m_eclHitToDigitAmp
	Amplitude of ECLDigit related to ECLHit.

std::vector< int > *	m_eclHitToPureDigit
	Index of ECLDigit related to ECLHit, PureCsI option.

std::vector< int > *	m_eclHitToPureDigitAmp
	Amplitude of ECLDigit related to ECLHit, PureCsI option.

std::vector< int > *	m_eclHitCellId
	ECLHit CellID.

std::vector< double > *	m_eclHitEnergyDep
	ECLHit energy.

std::vector< double > *	m_eclHitTimeAve
	ECLHit time.

int	m_eclClusterMultip
	Number of ECLClusters per event.

int	m_eclClusterTrueMultip
	Number of ECLClusters per event.

int	m_eclClusterGammaMultip
	Number of ECLClusters per event.

std::vector< int > *	m_eclClusterIdx
	ECLCluster index.

std::vector< int > *	m_eclClusterToMC1
	Index of first MCParticle related to ECLCluster.

std::vector< double > *	m_eclClusterToMCWeight1
	Energy contribution of first MCParticle related to ECLCluster.

std::vector< int > *	m_eclClusterToMC1PDG
	PDG code of first MCParticle related to ECLCluster.

std::vector< int > *	m_eclClusterToMC2
	Index of second MCParticle related to ECLCluster.

std::vector< double > *	m_eclClusterToMCWeight2
	Energy contribution of second MCParticle related to ECLCluster.

std::vector< int > *	m_eclClusterToMC2PDG
	PDG code of second MCParticle related to ECLCluster.

std::vector< int > *	m_eclClusterToMC3
	Index of third MCParticle related to ECLCluster.

std::vector< double > *	m_eclClusterToMCWeight3
	Energy contribution of third MCParticle related to ECLCluster.

std::vector< int > *	m_eclClusterToMC3PDG
	PDG code of third MCParticle related to ECLCluster.

std::vector< int > *	m_eclClusterToMC4
	Index of fourth MCParticle related to ECLCluster.

std::vector< double > *	m_eclClusterToMCWeight4
	Energy contribution of fourth MCParticle related to ECLCluster.

std::vector< int > *	m_eclClusterToMC4PDG
	PDG code of fourth MCParticle related to ECLCluster.

std::vector< int > *	m_eclClusterToMC5
	Index of fifth MCParticle related to ECLCluster.

std::vector< double > *	m_eclClusterToMCWeight5
	Energy contribution of 5th MCParticle related to ECLCluster.

std::vector< int > *	m_eclClusterToMC5PDG
	PDG code of fifth MCParticle related to ECLCluster.

std::vector< double > *	m_eclClusterToBkgWeight
	Remaining energy contribution not associated to first five MCParticles related to ECLCluster.

std::vector< double > *	m_eclClusterSimHitSum
	Energy contribution of 1st MCParticle related to ECLCluster.

std::vector< int > *	m_eclClusterToShower
	Index of ECLShower related to ECLCluster.

std::vector< double > *	m_eclClusterEnergy
	ECLCluster energy.

std::vector< double > *	m_eclClusterEnergyError
	ECLCluster energy error.

std::vector< double > *	m_eclClusterTheta
	ECLCluster polar direction.

std::vector< double > *	m_eclClusterThetaError
	ECLCluster error on polar direction.

std::vector< double > *	m_eclClusterPhi
	ECLCluster azimuthal direction.

std::vector< double > *	m_eclClusterPhiError
	ECLCluster error on azimuthal direction.

std::vector< double > *	m_eclClusterR
	ECLCluster distance from IP.

std::vector< double > *	m_eclClusterEnergyDepSum
	ECLCluster simulated energy.

std::vector< double > *	m_eclClusterTiming
	ECLCluster time.

std::vector< double > *	m_eclClusterTimingError
	ECLCluster time error.

std::vector< double > *	m_eclClusterE9oE21
	Ratio of 3x3 over 5x5 crystal matrices energies for ECLCluster.

std::vector< double > *	m_eclClusterHighestE
	Highest energy deposit (per crystal) in ECLCluster.

std::vector< int > *	m_eclClusterCellId
	CellId with highest energy deposit in ECLCluster.

std::vector< int > *	m_eclClusterNofCrystals
	Number of crystals in ECLCluster.

std::vector< int > *	m_eclClusterCrystalHealth
	Crystal health flag.

std::vector< bool > *	m_eclClusterIsTrack
	Flag for charged clusters.

std::vector< double > *	m_eclClusterClosestTrackDist
	Flag for charged clusters.

std::vector< double > *	m_eclClusterDeltaL
	Reconstructed Cluster DeltaL.

std::vector< double > *	m_eclClusterAbsZernike40
	Reconstructed Abs Zernike40.

std::vector< double > *	m_eclClusterAbsZernike51
	Reconstructed Abs Zernike51.

std::vector< double > *	m_eclClusterZernikeMVA
	Zernike MVA.

std::vector< double > *	m_eclClusterE1oE9
	Reconstructed E1 over E9.

std::vector< double > *	m_eclClusterSecondMoment
	Reconstructed Second Moment.

std::vector< double > *	m_eclClusterLAT
	Reconstructed LAT.

std::vector< double > *	m_eclClusterDeltaTime99
	DeltaTime99.

std::vector< int > *	m_eclClusterDetectorRegion
	Cluster Detector Region.

std::vector< int > *	m_eclClusterHasNPhotonHypothesis
	Cluster has n-photon hypothesis.

std::vector< int > *	m_eclClusterHasNeutralHadronHypothesis
	Cluster has neutral hadron hypothesis.

int	m_eclPureDigitMultip
	Number of ECLDigits per event, PureCsI option.

std::vector< int > *	m_eclPureDigitIdx
	ECLDigit index, PureCsI option.

std::vector< int > *	m_eclPureDigitToMC
	Index of MCParticle related to that ECLDigit, PureCsI option.

std::vector< int > *	m_eclPureDigitCellId
	Number of ECLDigit CellId, PureCsI option.

std::vector< int > *	m_eclPureDigitAmp
	ECLDigit amplitude, PureCsI option.

std::vector< int > *	m_eclPureDigitTimeFit
	ECLDigit timing, PureCsI option.

std::vector< int > *	m_eclPureDigitFitQuality
	ECLDigit fit quality, PureCsI option.

std::vector< int > *	m_eclPureDigitToCluster
	ECLDigit To Cluster, PureCsI option.

int	m_eclPureCalDigitMultip
	Number of ECLCalDigits per event, PureCsI option.

std::vector< int > *	m_eclPureCalDigitIdx
	ECLCalDigit index, PureCsI option.

std::vector< int > *	m_eclPureCalDigitToMC1
	Index of first MCParticle related to CalDigit, PureCsI option.

std::vector< int > *	m_eclPureCalDigitToMC1PDG
	PDG code of first MCParticle related to CalDigit, PureCsI option.

std::vector< double > *	m_eclPureCalDigitToMCWeight1
	Energy contribution of first MCParticle related to CalDigit, PureCsI option.

std::vector< int > *	m_eclPureCalDigitToMC2
	Index of second MCParticle related to CalDigit, PureCsI option.

std::vector< int > *	m_eclPureCalDigitToMC2PDG
	PDG code of second MCParticle related to CalDigit, PureCsI option.

std::vector< double > *	m_eclPureCalDigitToMCWeight2
	Energy contribution of second MCParticle related to CalDigit, PureCsI option.

std::vector< int > *	m_eclPureCalDigitToMC3
	Index of third MCParticle related to CalDigit, PureCsI option.

std::vector< int > *	m_eclPureCalDigitToMC3PDG
	PDG code of third MCParticle related to CalDigit, PureCsI option.

std::vector< double > *	m_eclPureCalDigitToMCWeight3
	Energy contribution of third MCParticle related to CalDigit, PureCsI option.

std::vector< int > *	m_eclPureCalDigitToMC4
	Index of fourth MCParticle related to CalDigit, PureCsI option.

std::vector< int > *	m_eclPureCalDigitToMC4PDG
	PDG code of fourth MCParticle related to CalDigit, PureCsI option.

std::vector< double > *	m_eclPureCalDigitToMCWeight4
	Energy contribution of fourth MCParticle related to CalDigit, PureCsI option.

std::vector< int > *	m_eclPureCalDigitToMC5
	Index of fifth MCParticle related to CalDigit, PureCsI option.

std::vector< int > *	m_eclPureCalDigitToMC5PDG
	PDG code of fifth MCParticle related to CalDigit, PureCsI option.

std::vector< double > *	m_eclPureCalDigitToMCWeight5
	Energy contribution of fifth MCParticle related to CalDigit, PureCsI option.

std::vector< double > *	m_eclPureCalDigitToBkgWeight
	Remaining energy contribution not associated to first five MCParticles related to CalDigit, PureCsI option.

std::vector< double > *	m_eclPureCalDigitSimHitSum
	Full energy contribution related to CalDigit, PureCsI option.

std::vector< int > *	m_eclPureCalDigitToShower
	Index of ECLShower related to that CalDigit, PureCsI option.

std::vector< int > *	m_eclPureCalDigitCellId
	Number of CalDigit CellId, PureCsI option.

std::vector< double > *	m_eclPureCalDigitAmp
	CalDigit amplitude, PureCsI option.

std::vector< double > *	m_eclPureCalDigitTimeFit
	CalDigit timing, PureCsI option.

std::vector< int > *	m_eclPureCalDigitFitQuality
	CalDigit fit quality, PureCsI option.

std::vector< int > *	m_eclPureCalDigitToCR
	Index of CR related to that CalDigit, PureCsI option.

std::vector< int > *	m_eclPureCalDigitToLM
	Index of LM related to that CalDigit, PureCsI option.

std::vector< int > *	m_eclPureCRIdx
	Connected Region ID, PureCsI option.

std::vector< int > *	m_eclPureCRIsTrack
	Int for Connected Region - Track Match, PureCsI option.

std::vector< double > *	m_eclPureCRLikelihoodMIPNGamma
	Connected Region MIP Likelihood, PureCsI option.

std::vector< double > *	m_eclPureCRLikelihoodChargedHadron
	Connected Region Charged Hadron Likelihood, PureCsI option.

std::vector< double > *	m_eclPureCRLikelihoodElectronNGamma
	Connected Region Electron Likelihood, PureCsI option.

std::vector< double > *	m_eclPureCRLikelihoodNGamma
	Connected Region Gamma Likelihood, PureCsI option.

std::vector< double > *	m_eclPureCRLikelihoodNeutralHadron
	Connected Region Neutral Hadron Likelihood, PureCsI option.

std::vector< double > *	m_eclPureCRLikelihoodMergedPi0
	Connected Region Merged Pi0 Likelihood, PureCsI option.

int	m_eclPureLMMultip
	Local Maxima multiplicity, PureCsI option.

std::vector< int > *	m_eclPureLMId
	Local Maximum ID, PureCsI option.

std::vector< int > *	m_eclPureLMType
	Local Maximum type, PureCsI option.

std::vector< int > *	m_eclPureLMCellId
	Local Maximum Cell ID, PureCsI option.

int	m_eclPureClusterMultip
	Number of ECLClusterss per event, PureCsI option.

std::vector< int > *	m_eclPureClusterIdx
	ECLCluster index, PureCsI option.

std::vector< int > *	m_eclPureClusterToMC1
	Index of first MCParticle related to ECLCluster, PureCsI option.

std::vector< double > *	m_eclPureClusterToMCWeight1
	Energy contribution of first MCParticle related to ECLCluster, PureCsI option.

std::vector< int > *	m_eclPureClusterToMC1PDG
	PDG code of first MCParticle related to ECLCluster, PureCsI option.

std::vector< int > *	m_eclPureClusterToMC2
	Index of second MCParticle related to ECLCluster, PureCsI option.

std::vector< double > *	m_eclPureClusterToMCWeight2
	Energy contribution of second MCParticle related to ECLCluster, PureCsI option.

std::vector< int > *	m_eclPureClusterToMC2PDG
	PDG code of second MCParticle related to ECLCluster, PureCsI option.

std::vector< int > *	m_eclPureClusterToMC3
	Index of third MCParticle related to ECLCluster, PureCsI option.

std::vector< double > *	m_eclPureClusterToMCWeight3
	Energy contribution of third MCParticle related to ECLCluster, PureCsI option.

std::vector< int > *	m_eclPureClusterToMC3PDG
	PDG code of third MCParticle related to ECLCluster, PureCsI option.

std::vector< int > *	m_eclPureClusterToMC4
	Index of fourth MCParticle related to ECLCluster, PureCsI option.

std::vector< double > *	m_eclPureClusterToMCWeight4
	Energy contribution of fourth MCParticle related to ECLCluster, PureCsI option.

std::vector< int > *	m_eclPureClusterToMC4PDG
	PDG code of fourth MCParticle related to ECLCluster, PureCsI option.

std::vector< int > *	m_eclPureClusterToMC5
	Index of fifth MCParticle related to ECLCluster, PureCsI option.

std::vector< double > *	m_eclPureClusterToMCWeight5
	Energy contribution of 5th MCParticle related to ECLCluster, PureCsI option.

std::vector< int > *	m_eclPureClusterToMC5PDG
	PDG code of fifth MCParticle related to ECLCluster, PureCsI option.

std::vector< double > *	m_eclPureClusterToBkgWeight
	Remaining energy contribution not associated to first five MCParticles related to ECLCluster, PureCsI option.

std::vector< double > *	m_eclPureClusterEnergy
	Cluster energy, PureCsI option.

std::vector< double > *	m_eclPureClusterEnergyError
	Cluster energy error, PureCsI option.

std::vector< double > *	m_eclPureClusterTheta
	Cluster polar direction, PureCsI option.

std::vector< double > *	m_eclPureClusterThetaError
	Cluster error on polar direction, PureCsI option.

std::vector< double > *	m_eclPureClusterPhi
	Cluster azimuthal direction, PureCsI option.

std::vector< double > *	m_eclPureClusterPhiError
	Cluster error on azimuthal direction, PureCsI option.

std::vector< double > *	m_eclPureClusterR
	Cluster distance from IP, PureCsI option.

std::vector< double > *	m_eclPureClusterEnergyDepSum
	Cluster simulated energy, PureCsI option.

std::vector< double > *	m_eclPureClusterTiming
	Cluster time, PureCsI option.

std::vector< double > *	m_eclPureClusterTimingError
	Cluster time error, PureCsI option.

std::vector< double > *	m_eclPureClusterE9oE21
	Ratio of 3x3 over 5x5 crystal matrices energies for Cluster, PureCsI option.

std::vector< double > *	m_eclPureClusterHighestE
	Highest energy deposit (per crystal) in Cluster, PureCsI option.

std::vector< int > *	m_eclPureClusterCellId
	CellId with highest energy deposit in Cluster, PureCsI option.

std::vector< double > *	m_eclPureClusterLat
	Cluster shape parameter LAT, PureCsI option.

std::vector< int > *	m_eclPureClusterNofCrystals
	Number of crystals in Cluster, PureCsI option.

std::vector< int > *	m_eclPureClusterCrystalHealth
	Crystal health flag, PureCsI option.

std::vector< bool > *	m_eclPureClusterIsTrack
	Flag for charged clusters, PureCsI option.

std::vector< double > *	m_eclPureClusterDeltaL
	Reconstructed Cluster DeltaL, PureCsI option.

std::vector< double > *	m_eclPureClusterClosestTrackDist
	Reconstructed Distance to Closest Track, PureCsI option.

std::vector< double > *	m_eclPureClusterAbsZernike40
	Reconstructed Zernike40, PureCsI option.

std::vector< double > *	m_eclPureClusterAbsZernike51
	Reconstructed Zernike51, PureCsI option.

std::vector< double > *	m_eclPureClusterZernikeMVA
	Output of MVA classifier based on Zernike Momenta, PureCsI option.

std::vector< double > *	m_eclPureClusterSecondMoment
	Reconstructed Cluster Second Moment, PureCsI option.

std::vector< double > *	m_eclPureClusterE1oE9
	Reconstructed E1oE9, PureCsI option.

std::vector< double > *	m_eclPureClusterDeltaTime99
	Reconstructed DeltaT99, PureCsI option.

std::vector< int > *	m_eclPureClusterDetectorRegion
	Clusters detector region, PureCsI option.

std::vector< int > *	m_eclPureClusterHasNPhotonHypothesis
	Cluster has n-photon hypothesis, PureCsI option.

std::vector< int > *	m_eclPureClusterHasNeutralHadronHypothesis
	Cluster has neutral hadron hypothesis, PureCsI option.

int	m_eclShowerMultip
	Number of ECLShowers per event.

std::vector< int > *	m_eclShowerIdx
	Shower Index.

std::vector< int > *	m_eclShowerToMC1
	Index of first MCParticle related to ECLShower.

std::vector< double > *	m_eclShowerToMCWeight1
	Energy contribution of first MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC1PDG
	PDG code of first MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC1Moth
	Mother index of first MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC1MothPDG
	PDG code of parent of first MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC1GMoth
	GMother index of first MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC1GMothPDG
	PDG code of Gparent of first MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC2
	Index of second MCParticle related to ECLShower.

std::vector< double > *	m_eclShowerToMCWeight2
	Energy contribution of second MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC2PDG
	PDG code of second MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC2Moth
	Mother index of first MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC2MothPDG
	PDG code of parent of first MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC2GMoth
	GMother index of second MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC2GMothPDG
	PDG code of Gparent of second MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC3
	Index of third MCParticle related to ECLShower.

std::vector< double > *	m_eclShowerToMCWeight3
	Energy contribution of third MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC3PDG
	PDG code of third MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC3Moth
	Mother index of first MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC3MothPDG
	PDG code of parent of first MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC3GMoth
	GMother index of third MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC3GMothPDG
	PDG code of Gparent of third MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC4
	Index of fourth MCParticle related to ECLShower.

std::vector< double > *	m_eclShowerToMCWeight4
	Energy contribution of fourth MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC4PDG
	PDG code of fourth MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC4Moth
	Mother index of fourth MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC4MothPDG
	PDG code of parent of fourth MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC4GMoth
	GMother index of fourth MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC4GMothPDG
	PDG code of Gparent of fourth MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC5
	Index of fifth MCParticle related to ECLShower.

std::vector< double > *	m_eclShowerToMCWeight5
	Energy contribution of fifth MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC5PDG
	PDG code of fifth MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC5Moth
	Mother index of fifth MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC5MothPDG
	PDG code of parent of fifth MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC5GMoth
	GMother index of fifth MCParticle related to ECLShower.

std::vector< int > *	m_eclShowerToMC5GMothPDG
	PDG code of Gparent of fifth MCParticle related to ECLShower.

std::vector< double > *	m_eclShowerToBkgWeight
	Remaining energy contribution not associated to first five MCParticles related to ECLShower.

std::vector< int > *	m_eclShowerToLM1
	Index of first maximum related to ECLShower.

std::vector< int > *	m_eclShowerToLM2
	Index of 2nd maximum related to ECLShower.

std::vector< int > *	m_eclShowerToLM3
	Index of 3rd maximum related to ECLShower.

std::vector< int > *	m_eclShowerToLM4
	Index of 4th maximum related to ECLShower.

std::vector< int > *	m_eclShowerToLM5
	Index of 5th maximum related to ECLShower.

std::vector< double > *	m_eclShowerSimHitSum
	Full energy contribution related to ECLShower.

std::vector< double > *	m_eclShowerUncEnergy
	Shower bare energy.

std::vector< double > *	m_eclShowerEnergy
	Shower Energy.

std::vector< double > *	m_eclShowerTheta
	Shower Theta.

std::vector< double > *	m_eclShowerPhi
	Shower Phi.

std::vector< double > *	m_eclShowerR
	Shower R.

std::vector< double > *	m_eclShowerNHits
	Shower NHits.

std::vector< double > *	m_eclShowerE9oE21
	Shower E9oE21.

std::vector< double > *	m_eclShowerTime
	Shower Timing.

std::vector< double > *	m_eclShowerT99
	Shower T99.

std::vector< int > *	m_eclShowerConnectedRegionId
	Matched Connetcted Region Idx.

std::vector< int > *	m_eclShowerHypothesisId
	Shower Particle Hypothesis ID.

std::vector< int > *	m_eclShowerCentralCellId
	Cell ID for most energetic crystal.

std::vector< double > *	m_eclShowerEnergyError
	Shower Energy Error.

std::vector< double > *	m_eclShowerThetaError
	Shower Theta Error.

std::vector< double > *	m_eclShowerPhiError
	Shower Phi Error.

std::vector< double > *	m_eclShowerTimeResolution
	Shower Time Resolution.

std::vector< double > *	m_eclShowerHighestEnergy
	Shower Highest Energy Crystal Energy.

std::vector< double > *	m_eclShowerLateralEnergy
	Shower Lateral Energy.

std::vector< double > *	m_eclShowerMinTrkDistance
	Shower Min Dist to Track.

std::vector< double > *	m_eclShowerTrkDepth
	Shower Track Depth.

std::vector< double > *	m_eclShowerShowerDepth
	Shower Depth.

std::vector< double > *	m_eclShowerAbsZernike40
	Shower Zernike40 Moment.

std::vector< double > *	m_eclShowerAbsZernike51
	Shower Zernike51 Moment.

std::vector< double > *	m_eclShowerZernikeMVA
	Shower ZernikeMVA.

std::vector< double > *	m_eclShowerSecondMoment
	Shower Second Moment.

std::vector< double > *	m_eclShowerE1oE9
	Shower E1/E9.

std::vector< int > *	m_eclShowerIsTrack
	Shower Track Match.

std::vector< bool > *	m_eclShowerIsCluster
	Shower Cluster Match.

std::vector< int > *	m_eclShowerMCVtxInEcl
	Int, 1 if particle decays (interacts) in ECL, 0 otherwise.

std::vector< int > *	m_eclShowerMCFlightMatch
	Int, 1 if particle flight direction is "well" reconstructed in ECL, 0 otherwise.

std::vector< int > *	m_eclShowerMCFFlightMatch
	Int, 1 if primary particle flight direction is "well" reconstructed in ECL, 0 otherwise, DEBUG PURPOSE.

std::vector< double > *	m_eclShowerHighestE1mE2
	Energy difference for 2 highest energy deposits in shower.

std::vector< double > *	m_eclShowerNumberOfCrystalsForEnergy
	Number of crystals used for energy calculation.

int	m_mcMultip
	Multiplicity of MCParticles.

std::vector< int > *	m_mcIdx
	MCParticle index.

std::vector< int > *	m_mcPdg
	MCParticle PDG code.

std::vector< int > *	m_mcMothPdg
	MCParticle mother particle PDG code.

std::vector< int > *	m_mcGMothPdg
	MCParticle grandmother particle PDG code.

std::vector< int > *	m_mcGGMothPdg
	MCParticle greand-grandmother particle PDG code.

std::vector< double > *	m_mcEnergy
	MCParticle energyx.

std::vector< double > *	m_mcPx
	MCParticle momentum X direction.

std::vector< double > *	m_mcPy
	MCParticle momentum Y direction.

std::vector< double > *	m_mcPz
	MCParticle momentum Z direction.

std::vector< double > *	m_mcDecayVtxX
	MCParticle decay vertex X.

std::vector< double > *	m_mcDecayVtxY
	MCParticle decay vertex Y.

std::vector< double > *	m_mcDecayVtxZ
	MCParticle decay vertex Z.

std::vector< double > *	m_mcProdVtxX
	MCParticle production vertex X.

std::vector< double > *	m_mcProdVtxY
	MCParticle production vertex Y.

std::vector< double > *	m_mcProdVtxZ
	MCParticle production vertex Z.

std::vector< int > *	m_mcSecondaryPhysProc
	Flag for secondary physics process.

int	m_trkMultip
	Track Multiplicity.

std::vector< int > *	m_trkPdg
	Track PDG code.

std::vector< int > *	m_trkIdx
	Track index.

std::vector< int > *	m_trkCharge
	Track charge.

std::vector< double > *	m_trkPx
	Track momentum along X direction.

std::vector< double > *	m_trkPy
	Track momentum along Y direction.

std::vector< double > *	m_trkPz
	Track momentum along Z direction.

std::vector< double > *	m_trkP
	Track momentum.

std::vector< double > *	m_trkTheta
	Track polar direction.

std::vector< double > *	m_trkPhi
	Track azimuthal direction.

std::vector< double > *	m_trkX
	Track DOCA X (?)

std::vector< double > *	m_trkY
	Track DOCA Y (?)

std::vector< double > *	m_trkZ
	Track DOCA Z (?)

std::vector< int > *	m_eclpidtrkIdx
	PID track index.

std::vector< double > *	m_eclpidEnergy
	PID track energy.

std::vector< double > *	m_eclpidEop
	PID track E/p.

std::vector< double > *	m_eclpidE9E21
	PID track ration of 3x3 over 5x5 crystal matrices energies.

std::vector< int > *	m_eclpidNCrystals
	PID track number of crystals.

std::vector< int > *	m_eclpidNClusters
	PID track number of clusters.

std::vector< double > *	m_eclLogLikeEl
	PID track electron likelihood.

std::vector< double > *	m_eclLogLikeMu
	PID track muon likelihood.

std::vector< double > *	m_eclLogLikePi
	PID track pion likelihood.

std::string	m_name
	The name of the module, saved as a string (user-modifiable)

std::string	m_type
	The type of the module, saved as a string.

std::string	m_package
	Package this module is found in (may be empty).

std::string	m_description
	The description of the module.

unsigned int	m_propertyFlags
	The properties of the module as bitwise or (with \|) of EModulePropFlags.

LogConfig	m_logConfig
	The log system configuration of the module.

ModuleParamList	m_moduleParamList
	List storing and managing all parameter of the module.

bool	m_hasReturnValue
	True, if the return value is set.

int	m_returnValue
	The return value.

std::vector< ModuleCondition >	m_conditions
	Module condition, only non-null if set.

Detailed Description

The ECL Data Analysis Module.

this module dump an ntuple containing ECL-related infos starting from mdst

Definition at line 51 of file ECLDataAnalysisModule.h.

Member Typedef Documentation

◆ EAfterConditionPath

typedef ModuleCondition::EAfterConditionPath EAfterConditionPath

inherited

Forward the EAfterConditionPath definition from the ModuleCondition.

Definition at line 88 of file Module.h.

Member Enumeration Documentation

◆ EModulePropFlags

enum EModulePropFlags

inherited

Each module can be tagged with property flags, which indicate certain features of the module.

Enumerator
c_Input	This module is an input module (reads data).
c_Output	This module is an output module (writes data).
c_ParallelProcessingCertified	This module can be run in parallel processing mode safely (All I/O must be done through the data store, in particular, the module must not write any files.)
c_HistogramManager	This module is used to manage histograms accumulated by other modules.
c_InternalSerializer	This module is an internal serializer/deserializer for parallel processing.
c_TerminateInAllProcesses	When using parallel processing, call this module's terminate() function in all processes(). This will also ensure that there is exactly one process (single-core if no parallel modules found) or at least one input, one main and one output process.
c_DontCollectStatistics	No statistics is collected for this module.

Definition at line 77 of file Module.h.

                          {
      c_Input                       = 1,  
      c_Output                      = 2,  
      c_ParallelProcessingCertified = 4,  
      c_HistogramManager            = 8, 
      c_InternalSerializer          = 16,  
      c_TerminateInAllProcesses     = 32,  
      c_DontCollectStatistics       = 64,  
    };

Constructor & Destructor Documentation

◆ ECLDataAnalysisModule()

ECLDataAnalysisModule ( )

Constructor of the module.

Definition at line 50 of file ECLDataAnalysisModule.cc.

  : Module(),
    m_rootFilePtr(0),
    m_writeToRoot(1),
    m_doTracking(0),
    m_eclSimHits(eclSimHitArrayName()),
    m_eclHits(eclHitArrayName()),
    m_eclDigits(eclDigitArrayName()),
    m_eclCalDigits(eclCalDigitArrayName()),
    m_eclConnectedRegions(eclConnectedRegionArrayName()),
    m_eclShowers(eclShowerArrayName()),
    m_eclClusters(eclClusterArrayName()),
    m_eclLocalMaximums(eclLocalMaximumArrayName()),
    m_eclPureDigits(eclPureDigitArrayName()),
    m_eclPureCalDigits(eclPureCalDigitArrayName()),
    m_eclPureConnectedRegions(eclPureConnectedRegionArrayName()),
    m_eclPureShowers(eclPureShowerArrayName()),
    m_eclPureClusters(eclPureClusterArrayName()),
    m_eclPureLocalMaximums(eclPureLocalMaximumArrayName()),
    m_tree(0),
    m_iExperiment(0),
    m_iRun(0),
    m_iEvent(0),
//Digit
    m_eclDigitMultip(0),
    m_eclDigitIdx(0),
    m_eclDigitToMC(0),
    m_eclDigitCellId(0),
    m_eclDigitAmp(0),
    m_eclDigitTimeFit(0),
    m_eclDigitFitQuality(0),
    m_eclDigitToCalDigit(0),
//CalDigit
    m_eclCalDigitMultip(0),
    m_eclCalDigitIdx(0),
    m_eclCalDigitToMC1(0),
    m_eclCalDigitToMC1PDG(0),
    m_eclCalDigitToMCWeight1(0),
    m_eclCalDigitToMC2(0),
    m_eclCalDigitToMC2PDG(0),
    m_eclCalDigitToMCWeight2(0),
    m_eclCalDigitToMC3(0),
    m_eclCalDigitToMC3PDG(0),
    m_eclCalDigitToMCWeight3(0),
    m_eclCalDigitToMC4(0),
    m_eclCalDigitToMC4PDG(0),
    m_eclCalDigitToMCWeight4(0),
    m_eclCalDigitToMC5(0),
    m_eclCalDigitToMC5PDG(0),
    m_eclCalDigitToMCWeight5(0),
    m_eclCalDigitToBkgWeight(0),
    m_eclCalDigitSimHitSum(0),
    m_eclCalDigitToShower(0),
    m_eclCalDigitCellId(0),
    m_eclCalDigitAmp(0),
    m_eclCalDigitTimeFit(0),
    m_eclCalDigitFitQuality(0),
    m_eclCalDigitToCR(0),
    m_eclCalDigitToLM(0),
//Connected Region
    m_eclCRIdx(0),
    m_eclCRIsTrack(0),
    m_eclCRLikelihoodMIPNGamma(0),
    m_eclCRLikelihoodChargedHadron(0),
    m_eclCRLikelihoodElectronNGamma(0),
    m_eclCRLikelihoodNGamma(0),
    m_eclCRLikelihoodNeutralHadron(0),
    m_eclCRLikelihoodMergedPi0(0),
//Local Maxima
    m_eclLMMultip(0),
    m_eclLMId(0),
    m_eclLMType(0),
    m_eclLMCellId(0),
//SimHit
    m_eclSimHitMultip(0),
    m_eclSimHitIdx(0),
    m_eclSimHitToMC(0),
    m_eclSimHitCellId(0),
    m_eclSimHitPdg(0),
    m_eclSimHitEnergyDep(0),
    m_eclSimHitFlightTime(0),
    m_eclSimHitX(0),
    m_eclSimHitY(0),
    m_eclSimHitZ(0),
    m_eclSimHitPx(0),
    m_eclSimHitPy(0),
    m_eclSimHitPz(0),
//Hit
    m_eclHitMultip(0),
    m_eclHitIdx(0),
    m_eclHitToMC(0),
    m_eclHitToDigit(0),
    m_eclHitToDigitAmp(0),
    m_eclHitToPureDigit(0),
    m_eclHitToPureDigitAmp(0),
    m_eclHitCellId(0),
    m_eclHitEnergyDep(0),
    m_eclHitTimeAve(0),
//Cluster
    m_eclClusterMultip(0),
    m_eclClusterTrueMultip(0),
    m_eclClusterGammaMultip(0),
    m_eclClusterIdx(0),
    m_eclClusterToMC1(0),
    m_eclClusterToMCWeight1(0),
    m_eclClusterToMC1PDG(0),
    m_eclClusterToMC2(0),
    m_eclClusterToMCWeight2(0),
    m_eclClusterToMC2PDG(0),
    m_eclClusterToMC3(0),
    m_eclClusterToMCWeight3(0),
    m_eclClusterToMC3PDG(0),
    m_eclClusterToMC4(0),
    m_eclClusterToMCWeight4(0),
    m_eclClusterToMC4PDG(0),
    m_eclClusterToMC5(0),
    m_eclClusterToMCWeight5(0),
    m_eclClusterToMC5PDG(0),
    m_eclClusterToBkgWeight(0),
    m_eclClusterSimHitSum(0),
    m_eclClusterToShower(0),
    m_eclClusterEnergy(0),
    m_eclClusterEnergyError(0),
    m_eclClusterTheta(0),
    m_eclClusterThetaError(0),
    m_eclClusterPhi(0),
    m_eclClusterPhiError(0),
    m_eclClusterR(0),
    m_eclClusterEnergyDepSum(0),
    m_eclClusterTiming(0),
    m_eclClusterTimingError(0),
    m_eclClusterE9oE21(0),
    m_eclClusterHighestE(0),
    m_eclClusterCellId(0),
    m_eclClusterNofCrystals(0),
    m_eclClusterCrystalHealth(0),
    m_eclClusterIsTrack(0),
    m_eclClusterClosestTrackDist(0),
    m_eclClusterDeltaL(0),
    m_eclClusterAbsZernike40(0),
    m_eclClusterAbsZernike51(0),
    m_eclClusterZernikeMVA(0),
    m_eclClusterE1oE9(0),
    m_eclClusterSecondMoment(0),
    m_eclClusterLAT(0),
    m_eclClusterDeltaTime99(0),
    m_eclClusterDetectorRegion(0),
    m_eclClusterHasNPhotonHypothesis(0),
    m_eclClusterHasNeutralHadronHypothesis(0),
 
//PureDigits
    m_eclPureDigitMultip(0),
    m_eclPureDigitIdx(0),
    m_eclPureDigitToMC(0),
    m_eclPureDigitCellId(0),
    m_eclPureDigitAmp(0),
    m_eclPureDigitTimeFit(0),
    m_eclPureDigitFitQuality(0),
    m_eclPureDigitToCluster(0),
//PureCalDigit
    m_eclPureCalDigitMultip(0),
    m_eclPureCalDigitIdx(0),
    m_eclPureCalDigitToMC1(0),
    m_eclPureCalDigitToMC1PDG(0),
    m_eclPureCalDigitToMCWeight1(0),
    m_eclPureCalDigitToMC2(0),
    m_eclPureCalDigitToMC2PDG(0),
    m_eclPureCalDigitToMCWeight2(0),
    m_eclPureCalDigitToMC3(0),
    m_eclPureCalDigitToMC3PDG(0),
    m_eclPureCalDigitToMCWeight3(0),
    m_eclPureCalDigitToMC4(0),
    m_eclPureCalDigitToMC4PDG(0),
    m_eclPureCalDigitToMCWeight4(0),
    m_eclPureCalDigitToMC5(0),
    m_eclPureCalDigitToMC5PDG(0),
    m_eclPureCalDigitToMCWeight5(0),
    m_eclPureCalDigitToBkgWeight(0),
    m_eclPureCalDigitSimHitSum(0),
    m_eclPureCalDigitToShower(0),
    m_eclPureCalDigitCellId(0),
    m_eclPureCalDigitAmp(0),
    m_eclPureCalDigitTimeFit(0),
    m_eclPureCalDigitFitQuality(0),
    m_eclPureCalDigitToCR(0),
    m_eclPureCalDigitToLM(0),
//Connected Region
    m_eclPureCRIdx(0),
    m_eclPureCRIsTrack(0),
    m_eclPureCRLikelihoodMIPNGamma(0),
    m_eclPureCRLikelihoodChargedHadron(0),
    m_eclPureCRLikelihoodElectronNGamma(0),
    m_eclPureCRLikelihoodNGamma(0),
    m_eclPureCRLikelihoodNeutralHadron(0),
    m_eclPureCRLikelihoodMergedPi0(0),
//Local Maxima
    m_eclPureLMMultip(0),
    m_eclPureLMId(0),
    m_eclPureLMType(0),
    m_eclPureLMCellId(0),
//PureCluster
    m_eclPureClusterMultip(0),
    m_eclPureClusterIdx(0),
    m_eclPureClusterToMC1(0),
    m_eclPureClusterToMCWeight1(0),
    m_eclPureClusterToMC1PDG(0),
    m_eclPureClusterToMC2(0),
    m_eclPureClusterToMCWeight2(0),
    m_eclPureClusterToMC2PDG(0),
    m_eclPureClusterToMC3(0),
    m_eclPureClusterToMCWeight3(0),
    m_eclPureClusterToMC3PDG(0),
    m_eclPureClusterToMC4(0),
    m_eclPureClusterToMCWeight4(0),
    m_eclPureClusterToMC4PDG(0),
    m_eclPureClusterToMC5(0),
    m_eclPureClusterToMCWeight5(0),
    m_eclPureClusterToMC5PDG(0),
    m_eclPureClusterToBkgWeight(0),
    m_eclPureClusterEnergy(0),
    m_eclPureClusterEnergyError(0),
    m_eclPureClusterTheta(0),
    m_eclPureClusterThetaError(0),
    m_eclPureClusterPhi(0),
    m_eclPureClusterPhiError(0),
    m_eclPureClusterR(0),
    m_eclPureClusterEnergyDepSum(0),
    m_eclPureClusterTiming(0),
    m_eclPureClusterTimingError(0),
    m_eclPureClusterE9oE21(0),
    m_eclPureClusterHighestE(0),
    m_eclPureClusterCellId(0),
    m_eclPureClusterLat(0),
    m_eclPureClusterNofCrystals(0),
    m_eclPureClusterCrystalHealth(0),
    m_eclPureClusterIsTrack(0),
    m_eclPureClusterDeltaL(0),
    m_eclPureClusterClosestTrackDist(0),
    m_eclPureClusterAbsZernike40(0),
    m_eclPureClusterAbsZernike51(0),
    m_eclPureClusterZernikeMVA(0),
    m_eclPureClusterSecondMoment(0),
    m_eclPureClusterE1oE9(0),
    m_eclPureClusterDeltaTime99(0),
    m_eclPureClusterDetectorRegion(0),
    m_eclPureClusterHasNPhotonHypothesis(0),
    m_eclPureClusterHasNeutralHadronHypothesis(0),
 
//Shower
    m_eclShowerMultip(0),
    m_eclShowerIdx(0),
    m_eclShowerToMC1(0),
    m_eclShowerToMCWeight1(0),
    m_eclShowerToMC1PDG(0),
    m_eclShowerToMC1Moth(0),
    m_eclShowerToMC1MothPDG(0),
    m_eclShowerToMC1GMoth(0),
    m_eclShowerToMC1GMothPDG(0),
    m_eclShowerToMC2(0),
    m_eclShowerToMCWeight2(0),
    m_eclShowerToMC2PDG(0),
    m_eclShowerToMC2Moth(0),
    m_eclShowerToMC2MothPDG(0),
    m_eclShowerToMC2GMoth(0),
    m_eclShowerToMC2GMothPDG(0),
    m_eclShowerToMC3(0),
    m_eclShowerToMCWeight3(0),
    m_eclShowerToMC3PDG(0),
    m_eclShowerToMC3Moth(0),
    m_eclShowerToMC3MothPDG(0),
    m_eclShowerToMC3GMoth(0),
    m_eclShowerToMC3GMothPDG(0),
    m_eclShowerToMC4(0),
    m_eclShowerToMCWeight4(0),
    m_eclShowerToMC4PDG(0),
    m_eclShowerToMC4Moth(0),
    m_eclShowerToMC4MothPDG(0),
    m_eclShowerToMC4GMoth(0),
    m_eclShowerToMC4GMothPDG(0),
    m_eclShowerToMC5(0),
    m_eclShowerToMCWeight5(0),
    m_eclShowerToMC5PDG(0),
    m_eclShowerToMC5Moth(0),
    m_eclShowerToMC5MothPDG(0),
    m_eclShowerToMC5GMoth(0),
    m_eclShowerToMC5GMothPDG(0),
    m_eclShowerToBkgWeight(0),
    m_eclShowerToLM1(0),
    m_eclShowerToLM2(0),
    m_eclShowerToLM3(0),
    m_eclShowerToLM4(0),
    m_eclShowerToLM5(0),
    m_eclShowerSimHitSum(0),
    m_eclShowerUncEnergy(0),
    m_eclShowerEnergy(0),
    m_eclShowerTheta(0),
    m_eclShowerPhi(0),
    m_eclShowerR(0),
    m_eclShowerNHits(0),
    m_eclShowerE9oE21(0),
    m_eclShowerTime(0),
    m_eclShowerT99(0),
    m_eclShowerConnectedRegionId(0),
    m_eclShowerHypothesisId(0),
    m_eclShowerCentralCellId(0),
    m_eclShowerEnergyError(0),
    m_eclShowerThetaError(0),
    m_eclShowerPhiError(0),
    m_eclShowerTimeResolution(0),
    m_eclShowerHighestEnergy(0),
    m_eclShowerLateralEnergy(0),
    m_eclShowerMinTrkDistance(0),
    m_eclShowerTrkDepth(0),
    m_eclShowerShowerDepth(0),
    m_eclShowerAbsZernike40(0),
    m_eclShowerAbsZernike51(0),
    m_eclShowerZernikeMVA(0),
    m_eclShowerSecondMoment(0),
    m_eclShowerE1oE9(0),
    m_eclShowerIsTrack(0),
    m_eclShowerIsCluster(0),
    m_eclShowerMCVtxInEcl(0),
    m_eclShowerMCFlightMatch(0),
    m_eclShowerMCFFlightMatch(0),
    m_eclShowerHighestE1mE2(0),
    m_eclShowerNumberOfCrystalsForEnergy(0),
//MC
    m_mcMultip(0),
    m_mcIdx(0),
    m_mcPdg(0),
    m_mcMothPdg(0),
    m_mcGMothPdg(0),
    m_mcGGMothPdg(0),
    m_mcEnergy(0),
    m_mcPx(0),
    m_mcPy(0),
    m_mcPz(0),
    m_mcDecayVtxX(0),
    m_mcDecayVtxY(0),
    m_mcDecayVtxZ(0),
    m_mcProdVtxX(0),
    m_mcProdVtxY(0),
    m_mcProdVtxZ(0),
    m_mcSecondaryPhysProc(0),
//Tracks
    m_trkMultip(0),
    m_trkPdg(0),
    m_trkIdx(0),
    m_trkCharge(0),
    m_trkPx(0),
    m_trkPy(0),
    m_trkPz(0),
    m_trkP(0),
    m_trkTheta(0),
    m_trkPhi(0),
    m_trkX(0),
    m_trkY(0),
    m_trkZ(0),
//PID
    m_eclpidtrkIdx(0),
    m_eclpidEnergy(0),
    m_eclpidEop(0),
    m_eclpidE9E21(0),
    m_eclpidNCrystals(0),
    m_eclpidNClusters(0),
    m_eclLogLikeEl(0),
    m_eclLogLikeMu(0),
    m_eclLogLikePi(0)
 
{
  //Set module properties
  setDescription("This module produces an ntuple with ECL-related quantities starting from mdst");
  addParam("writeToRoot", m_writeToRoot,
           "set true if you want to save the information in a root file named by parameter 'rootFileName'", bool(true));
  addParam("rootFileName", m_rootFileName,
           "fileName used for root file where info are saved. Will be ignored if parameter 'writeToRoot' is false (standard)",
           string("eclDataAnalysis"));
  addParam("doTracking", m_doTracking, "set true if you want to save the information from TrackFitResults'rootFileName'",
           bool(false));
  addParam("doHits", m_doHits, "set true if you want to save the Hit and SimHit informations'", bool(false));
  addParam("doDigits", m_doDigits, "set true if you want to save the Digits and CalDigits informations'", bool(false));
  addParam("doPureCsI", m_doPureCsI, "set true if you want to save the information for upgrade option'", bool(false));
 
}

◆ ~ECLDataAnalysisModule()

~ECLDataAnalysisModule ( )

overridevirtual

Destructor of the module.

Definition at line 435 of file ECLDataAnalysisModule.cc.

436{

437}

Member Function Documentation

◆ beginRun()

virtual void beginRun ( void )

inlinevirtualinherited

Called when entering a new run.

Called at the beginning of each run, the method gives you the chance to change run dependent constants like alignment parameters, etc.

This method can be implemented by subclasses.

Definition at line 146 of file Module.h.

146{};

◆ clone()

std::shared_ptr< PathElement > clone ( ) const

overridevirtualinherited

Create an independent copy of this module.

Note that parameters are shared, so changing them on a cloned module will also affect the original module.

Implements PathElement.

Definition at line 179 of file Module.cc.

{
  ModulePtr newModule = ModuleManager::Instance().registerModule(getType());
  newModule->m_moduleParamList.setParameters(getParamList());
  newModule->setName(getName());
  newModule->m_package = m_package;
  newModule->m_propertyFlags = m_propertyFlags;
  newModule->m_logConfig = m_logConfig;
  newModule->m_conditions = m_conditions;
 
  return newModule;
}

◆ def_beginRun()

virtual void def_beginRun ( )

inlineprotectedvirtualinherited

Wrapper method for the virtual function beginRun() that has the implementation to be used in a call from Python.

Reimplemented in PyModule.

Definition at line 425 of file Module.h.

425{ beginRun(); }

◆ def_endRun()

virtual void def_endRun ( )

inlineprotectedvirtualinherited

This method can receive that the current run ends as a call from the Python side.

For regular C++-Modules that forwards the call to the regular endRun() method.

Reimplemented in PyModule.

Definition at line 438 of file Module.h.

438{ endRun(); }

◆ def_event()

virtual void def_event ( )

inlineprotectedvirtualinherited

Wrapper method for the virtual function event() that has the implementation to be used in a call from Python.

Reimplemented in PyModule.

Definition at line 431 of file Module.h.

431{ event(); }

◆ def_initialize()

virtual void def_initialize ( )

inlineprotectedvirtualinherited

Wrappers to make the methods without "def_" prefix callable from Python.

Overridden in PyModule. Wrapper method for the virtual function initialize() that has the implementation to be used in a call from Python.

Reimplemented in PyModule.

Definition at line 419 of file Module.h.

419{ initialize(); }

◆ def_terminate()

virtual void def_terminate ( )

inlineprotectedvirtualinherited

Wrapper method for the virtual function terminate() that has the implementation to be used in a call from Python.

Reimplemented in PyModule.

Definition at line 444 of file Module.h.

444{ terminate(); }

◆ eclCalDigitArrayName()

virtual const char * eclCalDigitArrayName ( ) const

inlineprivatevirtual

Default name ECLCalDigits array.

Definition at line 113 of file ECLDataAnalysisModule.h.

114 { return "ECLCalDigits" ; }

◆ eclClusterArrayName()

virtual const char * eclClusterArrayName ( ) const

inlineprivatevirtual

Default name ECLClusters array.

Definition at line 119 of file ECLDataAnalysisModule.h.

120 { return "ECLClusters" ; }

◆ eclConnectedRegionArrayName()

virtual const char * eclConnectedRegionArrayName ( ) const

inlineprivatevirtual

Default name ECLConnectedRegions array.

Definition at line 122 of file ECLDataAnalysisModule.h.

123 { return "ECLConnectedRegions" ; }

◆ eclDigitArrayName()

virtual const char * eclDigitArrayName ( ) const

inlineprivatevirtual

Default name ECLDigits array.

Definition at line 110 of file ECLDataAnalysisModule.h.

111 { return "ECLDigits" ; }

◆ eclHitArrayName()

virtual const char * eclHitArrayName ( ) const

inlineprivatevirtual

Default name ECLCalDigits array.

Definition at line 107 of file ECLDataAnalysisModule.h.

108 { return "ECLHits" ; }

◆ eclLocalMaximumArrayName()

virtual const char * eclLocalMaximumArrayName ( ) const

inlineprivatevirtual

Default name ECLLocalMaxima array.

Definition at line 125 of file ECLDataAnalysisModule.h.

126 { return "ECLLocalMaximums" ; }

◆ eclPureCalDigitArrayName()

virtual const char * eclPureCalDigitArrayName ( ) const

inlineprivatevirtual

Default name ECLPureCalDigits array.

Definition at line 146 of file ECLDataAnalysisModule.h.

147 { return "ECLCalDigitsPureCsI" ; }

◆ eclPureClusterArrayName()

virtual const char * eclPureClusterArrayName ( ) const

inlineprivatevirtual

Default name ECLPureClusters array.

Definition at line 152 of file ECLDataAnalysisModule.h.

153 { return "ECLClustersPureCsI" ; }

◆ eclPureConnectedRegionArrayName()

virtual const char * eclPureConnectedRegionArrayName ( ) const

inlineprivatevirtual

Default name ECLPureConnectedRegions array.

Definition at line 155 of file ECLDataAnalysisModule.h.

156 { return "ECLConnectedRegionsPureCsI" ; }

◆ eclPureDigitArrayName()

virtual const char * eclPureDigitArrayName ( ) const

inlineprivatevirtual

Default name ECLPureDigits array.

Definition at line 143 of file ECLDataAnalysisModule.h.

144 { return "ECLDigitsPureCsI" ; }

◆ eclPureLocalMaximumArrayName()

virtual const char * eclPureLocalMaximumArrayName ( ) const

inlineprivatevirtual

Default name ECLPureLocalMaxima array.

Definition at line 158 of file ECLDataAnalysisModule.h.

159 { return "ECLLocalMaximumsPureCsI" ; }

◆ eclPureShowerArrayName()

virtual const char * eclPureShowerArrayName ( ) const

inlineprivatevirtual

Default name ECLPureShower array.

Definition at line 149 of file ECLDataAnalysisModule.h.

150 { return "ECLShowersPureCsI" ; }

◆ eclShowerArrayName()

virtual const char * eclShowerArrayName ( ) const

inlineprivatevirtual

Default name ECLShower array.

Definition at line 116 of file ECLDataAnalysisModule.h.

117 { return "ECLShowers" ; }

◆ eclSimHitArrayName()

virtual const char * eclSimHitArrayName ( ) const

inlineprivatevirtual

Default name ECLCalDigits array.

Definition at line 104 of file ECLDataAnalysisModule.h.

105 { return "ECLSimHits" ; }

◆ endRun()

virtual void endRun ( void )

inlinevirtualinherited

This method is called if the current run ends.

Use this method to store information, which should be aggregated over one run.

This method can be implemented by subclasses.

Definition at line 165 of file Module.h.

165{};

◆ evalCondition()

bool evalCondition ( ) const

inherited

If at least one condition was set, it is evaluated and true returned if at least one condition returns true.

If no condition or result value was defined, the method returns false. Otherwise, the condition is evaluated and true returned, if at least one condition returns true. To speed up the evaluation, the condition strings were already parsed in the method if_value().

Returns: True if at least one condition and return value exists and at least one condition expression was evaluated to true.

Definition at line 96 of file Module.cc.

{
  if (m_conditions.empty()) return false;
 
  //okay, a condition was set for this Module...
  if (!m_hasReturnValue) {
    B2FATAL("A condition was set for '" << getName() << "', but the module did not set a return value!");
  }
 
  for (const auto& condition : m_conditions) {
    if (condition.evaluate(m_returnValue)) {
      return true;
    }
  }
  return false;
}

◆ event()

void event ( void )

overridevirtual

event

Digits

CalDigits

Hit

Clusters

Showers

Pure Digits

Pure Clusters

MC

Tracks

PID

Reimplemented from Module.

Definition at line 861 of file ECLDataAnalysisModule.cc.

{
 
  B2DEBUG(1, "  ++++++++++++++ ECLDataAnalysisModule");
 
  //EventLevelClusteringInfo
  m_nECLCalDigitsOutOfTimeFWD = 0;
  m_nECLCalDigitsOutOfTimeBarrel = 0;
  m_nECLCalDigitsOutOfTimeBWD = 0;
  m_nECLShowersRejectedFWD = 0;
  m_nECLShowersRejectedBarrel = 0;
  m_nECLShowersRejectedBWD = 0;
 
  if (m_doDigits == 1) {
    m_eclDigitMultip = 0;
    m_eclDigitIdx->clear();
    m_eclDigitToMC->clear();
    m_eclDigitCellId->clear();
    m_eclDigitAmp->clear();
    m_eclDigitTimeFit->clear();
    m_eclDigitFitQuality->clear();
    m_eclDigitToCalDigit->clear();
    m_eclCalDigitMultip = 0;
    m_eclCalDigitCellId->clear();
    m_eclCalDigitAmp->clear();
    m_eclCalDigitTimeFit->clear();
    m_eclCalDigitFitQuality->clear();
    m_eclCalDigitIdx->clear();
    m_eclCalDigitToMC1->clear();
    m_eclCalDigitToMCWeight1->clear();
    m_eclCalDigitToMC2->clear();
    m_eclCalDigitToMCWeight2->clear();
    m_eclCalDigitToMC3->clear();
    m_eclCalDigitToMCWeight3->clear();
    m_eclCalDigitToMC4->clear();
    m_eclCalDigitToMCWeight4->clear();
    m_eclCalDigitToMC5->clear();
    m_eclCalDigitToMCWeight5->clear();
    m_eclCalDigitToMC1PDG->clear();
    m_eclCalDigitToMC2PDG->clear();
    m_eclCalDigitToMC3PDG->clear();
    m_eclCalDigitToMC4PDG->clear();
    m_eclCalDigitToMC5PDG->clear();
    m_eclCalDigitToBkgWeight->clear();
    m_eclCalDigitSimHitSum->clear();
    m_eclCalDigitToShower->clear();
    m_eclCalDigitToCR->clear();
    m_eclCalDigitToLM->clear();
  }
 
  m_eclCRIdx->clear();
  m_eclCRIsTrack->clear();
  m_eclCRLikelihoodMIPNGamma->clear();
  m_eclCRLikelihoodChargedHadron->clear();
  m_eclCRLikelihoodElectronNGamma->clear();
  m_eclCRLikelihoodNGamma->clear();
  m_eclCRLikelihoodNeutralHadron->clear();
  m_eclCRLikelihoodMergedPi0->clear();
 
  if (m_doHits == 1) {
    //SimHit
    m_eclSimHitMultip = 0;
    m_eclSimHitCellId->clear();
    m_eclSimHitPdg->clear();
    m_eclSimHitEnergyDep->clear();
    m_eclSimHitFlightTime->clear();
    m_eclSimHitIdx->clear();
    m_eclSimHitToMC->clear();
    m_eclSimHitX->clear();
    m_eclSimHitY->clear();
    m_eclSimHitZ->clear();
    m_eclSimHitPx->clear();
    m_eclSimHitPy->clear();
    m_eclSimHitPz->clear();
    m_eclHitMultip = 0;
    m_eclHitCellId->clear();
    m_eclHitEnergyDep->clear();
    m_eclHitTimeAve->clear();
    m_eclHitIdx->clear();
    m_eclHitToMC->clear();
    m_eclHitToDigit->clear();
    m_eclHitToDigitAmp->clear();
  }
 
  m_eclLMMultip = 0;
  m_eclLMId->clear();
  m_eclLMType->clear();
  m_eclLMCellId->clear();
  m_eclClusterMultip = 0;
  m_eclClusterTrueMultip = 0;
  m_eclClusterGammaMultip = 0;
  m_eclClusterEnergy->clear();
  m_eclClusterEnergyError->clear();
  m_eclClusterTheta->clear();
  m_eclClusterThetaError->clear();
  m_eclClusterPhi->clear();
  m_eclClusterPhiError->clear();
  m_eclClusterR->clear();
  m_eclClusterIdx->clear();
  m_eclClusterToMC1->clear();
  m_eclClusterToMCWeight1->clear();
  m_eclClusterToMC1PDG->clear();
  m_eclClusterToMC2->clear();
  m_eclClusterToMCWeight2->clear();
  m_eclClusterToMC2PDG->clear();
  m_eclClusterToMC3->clear();
  m_eclClusterToMCWeight3->clear();
  m_eclClusterToMC3PDG->clear();
  m_eclClusterToMC4->clear();
  m_eclClusterToMCWeight4->clear();
  m_eclClusterToMC4PDG->clear();
  m_eclClusterToMC5->clear();
  m_eclClusterToMCWeight5->clear();
  m_eclClusterToMC5PDG->clear();
  m_eclClusterToBkgWeight->clear();
  m_eclClusterSimHitSum->clear();
  m_eclClusterEnergyDepSum->clear();
  m_eclClusterTiming->clear();
  m_eclClusterTimingError->clear();
  m_eclClusterE9oE21->clear();
  m_eclClusterHighestE->clear();
  m_eclClusterCellId->clear();
  m_eclClusterNofCrystals->clear();
  m_eclClusterCrystalHealth->clear();
  m_eclClusterIsTrack->clear();
  m_eclClusterClosestTrackDist->clear();
  m_eclClusterDeltaL->clear();
  m_eclClusterToShower->clear();
  m_eclClusterAbsZernike40->clear();
  m_eclClusterAbsZernike51->clear();
  m_eclClusterZernikeMVA->clear();
  m_eclClusterE1oE9->clear();
  m_eclClusterSecondMoment->clear();
  m_eclClusterLAT->clear();
  m_eclClusterDeltaTime99->clear();
  m_eclClusterDetectorRegion->clear();
  m_eclClusterHasNPhotonHypothesis->clear();
  m_eclClusterHasNeutralHadronHypothesis->clear();
  m_eclShowerMultip = 0;
  m_eclShowerIdx->clear();
  m_eclShowerToMC1->clear();
  m_eclShowerToMCWeight1->clear();
  m_eclShowerToMC1PDG->clear();
  m_eclShowerToMC1Moth->clear();
  m_eclShowerToMC1MothPDG->clear();
  m_eclShowerToMC1GMoth->clear();
  m_eclShowerToMC1GMothPDG->clear();
  m_eclShowerToMC2->clear();
  m_eclShowerToMCWeight2->clear();
  m_eclShowerToMC2PDG->clear();
  m_eclShowerToMC2Moth->clear();
  m_eclShowerToMC2MothPDG->clear();
  m_eclShowerToMC2GMoth->clear();
  m_eclShowerToMC2GMothPDG->clear();
  m_eclShowerToMC3->clear();
  m_eclShowerToMCWeight3->clear();
  m_eclShowerToMC3PDG->clear();
  m_eclShowerToMC3Moth->clear();
  m_eclShowerToMC3MothPDG->clear();
  m_eclShowerToMC3GMoth->clear();
  m_eclShowerToMC3GMothPDG->clear();
  m_eclShowerToMC4->clear();
  m_eclShowerToMCWeight4->clear();
  m_eclShowerToMC4PDG->clear();
  m_eclShowerToMC4Moth->clear();
  m_eclShowerToMC4MothPDG->clear();
  m_eclShowerToMC4GMoth->clear();
  m_eclShowerToMC4GMothPDG->clear();
  m_eclShowerToMC5->clear();
  m_eclShowerToMCWeight5->clear();
  m_eclShowerToMC5PDG->clear();
  m_eclShowerToMC5Moth->clear();
  m_eclShowerToMC5MothPDG->clear();
  m_eclShowerToMC5GMoth->clear();
  m_eclShowerToMC5GMothPDG->clear();
  m_eclShowerToBkgWeight->clear();
  m_eclShowerToLM1->clear();
  m_eclShowerToLM2->clear();
  m_eclShowerToLM3->clear();
  m_eclShowerToLM4->clear();
  m_eclShowerToLM5->clear();
  m_eclShowerSimHitSum->clear();
  m_eclShowerUncEnergy->clear();
  m_eclShowerEnergy->clear();
  m_eclShowerTheta->clear();
  m_eclShowerPhi->clear();
  m_eclShowerR->clear();
  m_eclShowerNHits->clear();
  m_eclShowerE9oE21->clear();
  m_eclShowerTime->clear();
  m_eclShowerT99->clear();
  m_eclShowerConnectedRegionId->clear();
  m_eclShowerHypothesisId->clear();
  m_eclShowerCentralCellId->clear();
  m_eclShowerEnergyError->clear();
  m_eclShowerThetaError->clear();
  m_eclShowerPhiError->clear();
  m_eclShowerTimeResolution->clear();
  m_eclShowerHighestEnergy->clear();
  m_eclShowerLateralEnergy->clear();
  m_eclShowerMinTrkDistance->clear();
  m_eclShowerTrkDepth->clear();
  m_eclShowerShowerDepth->clear();
  m_eclShowerAbsZernike40->clear();
  m_eclShowerAbsZernike51->clear();
  m_eclShowerZernikeMVA->clear();
  m_eclShowerSecondMoment->clear();
  m_eclShowerE1oE9->clear();
  m_eclShowerIsTrack->clear();
  m_eclShowerIsCluster->clear();
  m_eclShowerMCVtxInEcl->clear();
  m_eclShowerMCFlightMatch->clear();
  m_eclShowerMCFFlightMatch->clear();
  m_eclShowerHighestE1mE2->clear();
  m_eclShowerNumberOfCrystalsForEnergy->clear();
  if (m_doPureCsI == true) {
    m_eclPureDigitMultip = 0;
    m_eclPureDigitIdx->clear();
    m_eclPureDigitToMC->clear();
    m_eclPureDigitCellId->clear();
    m_eclPureDigitAmp->clear();
    m_eclPureDigitTimeFit->clear();
    m_eclPureDigitFitQuality->clear();
    m_eclPureDigitToCluster->clear();
    m_eclHitToPureDigit->clear();
    m_eclHitToPureDigitAmp->clear();
 
    //Pure CalDigits
    m_eclPureCalDigitMultip = 0;
    m_eclPureCalDigitCellId->clear();
    m_eclPureCalDigitAmp->clear();
    m_eclPureCalDigitTimeFit->clear();
    m_eclPureCalDigitFitQuality->clear();
    m_eclPureCalDigitIdx->clear();
    m_eclPureCalDigitToMC1->clear();
    m_eclPureCalDigitToMCWeight1->clear();
    m_eclPureCalDigitToMC2->clear();
    m_eclPureCalDigitToMCWeight2->clear();
    m_eclPureCalDigitToMC3->clear();
    m_eclPureCalDigitToMCWeight3->clear();
    m_eclPureCalDigitToMC4->clear();
    m_eclPureCalDigitToMCWeight4->clear();
    m_eclPureCalDigitToMC5->clear();
    m_eclPureCalDigitToMCWeight5->clear();
    m_eclPureCalDigitToMC1PDG->clear();
    m_eclPureCalDigitToMC2PDG->clear();
    m_eclPureCalDigitToMC3PDG->clear();
    m_eclPureCalDigitToMC4PDG->clear();
    m_eclPureCalDigitToMC5PDG->clear();
    m_eclPureCalDigitToBkgWeight->clear();
    m_eclPureCalDigitSimHitSum->clear();
    m_eclPureCalDigitToShower->clear();
    m_eclPureCalDigitToCR->clear();
    m_eclPureCalDigitToLM->clear();
 
    m_eclPureCRIdx->clear();
    m_eclPureCRIsTrack->clear();
    m_eclPureCRLikelihoodMIPNGamma->clear();
    m_eclPureCRLikelihoodChargedHadron->clear();
    m_eclPureCRLikelihoodElectronNGamma->clear();
    m_eclPureCRLikelihoodNGamma->clear();
    m_eclPureCRLikelihoodNeutralHadron->clear();
    m_eclPureCRLikelihoodMergedPi0->clear();
 
    m_eclPureLMMultip = 0;
    m_eclPureLMId->clear();
    m_eclPureLMType->clear();
    m_eclPureLMCellId->clear();
    m_eclPureClusterMultip = 0;
    m_eclPureClusterToMC1->clear();
    m_eclPureClusterToMCWeight1->clear();
    m_eclPureClusterToMC1PDG->clear();
    m_eclPureClusterToMC2->clear();
    m_eclPureClusterToMCWeight2->clear();
    m_eclPureClusterToMC2PDG->clear();
    m_eclPureClusterToMC3->clear();
    m_eclPureClusterToMCWeight3->clear();
    m_eclPureClusterToMC3PDG->clear();
    m_eclPureClusterToMC4->clear();
    m_eclPureClusterToMCWeight4->clear();
    m_eclPureClusterToMC4PDG->clear();
    m_eclPureClusterToMC5->clear();
    m_eclPureClusterToMCWeight5->clear();
    m_eclPureClusterToMC5PDG->clear();
    m_eclPureClusterToBkgWeight->clear();
    m_eclPureClusterEnergy->clear();
    m_eclPureClusterEnergyError->clear();
    m_eclPureClusterTheta->clear();
    m_eclPureClusterEnergyDepSum->clear();
    m_eclPureClusterThetaError->clear();
    m_eclPureClusterPhi->clear();
    m_eclPureClusterPhiError->clear();
    m_eclPureClusterR->clear();
    m_eclPureClusterIdx->clear();
    m_eclPureClusterTiming->clear();
    m_eclPureClusterTimingError->clear();
    m_eclPureClusterE9oE21->clear();
    m_eclPureClusterHighestE->clear();
    m_eclPureClusterCellId->clear();
    m_eclPureClusterLat->clear();
    m_eclPureClusterNofCrystals->clear();
    m_eclPureClusterCrystalHealth->clear();
    m_eclPureClusterClosestTrackDist->clear();
    m_eclPureClusterAbsZernike40->clear();
    m_eclPureClusterAbsZernike51->clear();
    m_eclPureClusterZernikeMVA->clear();
    m_eclPureClusterSecondMoment->clear();
    m_eclPureClusterIsTrack->clear();
    m_eclPureClusterDeltaL->clear();
    m_eclPureClusterE1oE9->clear();
    m_eclPureClusterDeltaTime99->clear();
    m_eclPureClusterDetectorRegion->clear();
    m_eclPureClusterHasNPhotonHypothesis->clear();
    m_eclPureClusterHasNeutralHadronHypothesis->clear();
  }
  m_mcMultip = 0;
  m_mcIdx->clear();
  m_mcPdg->clear();
  m_mcMothPdg->clear();
  m_mcGMothPdg->clear();
  m_mcGGMothPdg->clear();
  m_mcEnergy->clear();
  m_mcPx->clear();
  m_mcPy->clear();
  m_mcPz->clear();
  m_mcDecayVtxX->clear();
  m_mcDecayVtxY->clear();
  m_mcDecayVtxZ->clear();
  m_mcProdVtxX->clear();
  m_mcProdVtxY->clear();
  m_mcProdVtxZ->clear();
  m_mcSecondaryPhysProc->clear();
  if (m_doTracking == true) {
    m_trkMultip = 0;
    m_trkIdx->clear();
    m_trkPdg->clear();
    m_trkCharge->clear();
    m_trkPx->clear();
    m_trkPy->clear();
    m_trkPz->clear();
    m_trkP->clear();
    m_trkPhi->clear();
    m_trkTheta->clear();
    m_trkPhi->clear();
    m_trkX->clear();
    m_trkY->clear();
    m_trkZ->clear();
    m_eclpidtrkIdx->clear();
    m_eclpidEnergy->clear();
    m_eclpidEop->clear();
    m_eclpidE9E21->clear();
    m_eclpidNCrystals->clear();
    m_eclpidNClusters->clear();
    m_eclLogLikeEl->clear();
    m_eclLogLikeMu->clear();
    m_eclLogLikePi->clear();
  }
 
 
  if (m_eventmetadata) {
    m_iExperiment = m_eventmetadata->getExperiment();
    m_iRun = m_eventmetadata->getRun();
    m_iEvent = m_eventmetadata->getEvent();
  } else {
    m_iExperiment = -1;
    m_iRun = -1;
    m_iEvent = -1;
  }
 
  //EventLevelClusteringInfo
  m_nECLCalDigitsOutOfTimeFWD = m_eventLevelClusteringInfo->getNECLCalDigitsOutOfTimeFWD();
  m_nECLCalDigitsOutOfTimeBarrel = m_eventLevelClusteringInfo->getNECLCalDigitsOutOfTimeBarrel();
  m_nECLCalDigitsOutOfTimeBWD  = m_eventLevelClusteringInfo->getNECLCalDigitsOutOfTimeBWD();
  m_nECLShowersRejectedFWD = m_eventLevelClusteringInfo->getNECLShowersRejectedFWD();
  m_nECLShowersRejectedBarrel = m_eventLevelClusteringInfo->getNECLShowersRejectedBarrel();
  m_nECLShowersRejectedBWD = m_eventLevelClusteringInfo->getNECLShowersRejectedBWD();
 
  if (m_doDigits == 1) {
    //DIGITS
    m_eclDigitMultip = m_eclDigits.getEntries();
    for (int idigits = 0; idigits < m_eclDigits.getEntries() ; idigits++) {
      ECLDigit* aECLDigits = m_eclDigits[idigits];
      m_eclDigitIdx->push_back(idigits);
      m_eclDigitCellId->push_back(aECLDigits->getCellId());
      m_eclDigitAmp->push_back(aECLDigits->getAmp());
      m_eclDigitTimeFit->push_back(aECLDigits->getTimeFit());
      m_eclDigitFitQuality->push_back(aECLDigits->getQuality());
 
      if (aECLDigits->getRelated<MCParticle>() != (nullptr)) {
        const MCParticle* mc_digit = aECLDigits->getRelated<MCParticle>();
        m_eclDigitToMC->push_back(mc_digit->getArrayIndex());
      } else
        m_eclDigitToMC->push_back(-1);
 
      if (aECLDigits->getRelated<ECLCalDigit>() != (nullptr)) {
        const ECLCalDigit* cal_digit = aECLDigits->getRelated<ECLCalDigit>();
        m_eclDigitToCalDigit->push_back(cal_digit->getArrayIndex());
      } else
        m_eclDigitToCalDigit->push_back(-1);
    }
 
    //CAL DIGITS
    m_eclCalDigitMultip = m_eclCalDigits.getEntries();
    for (uint icaldigits = 0; icaldigits < (uint)m_eclCalDigits.getEntries() ; icaldigits++) {
      ECLCalDigit* aECLCalDigits = m_eclCalDigits[icaldigits];
 
      m_eclCalDigitIdx->push_back(icaldigits);
      m_eclCalDigitCellId->push_back(aECLCalDigits->getCellId());
      m_eclCalDigitAmp->push_back(aECLCalDigits->getEnergy());
      m_eclCalDigitTimeFit->push_back(aECLCalDigits->getTime());
      m_eclCalDigitFitQuality->push_back(aECLCalDigits->isFailedFit());
 
      double sumHit = 0;
      int idx[10];
      for (int i = 0; i < 10; i++)
        idx[i] = -1;
 
      double wi[10];
      for (int i = 0; i < 10; i++)
        wi[i] = -1;
 
      int ii = 0;
      sumHit = 0;
 
      auto digitMCRelations = aECLCalDigits->getRelationsTo<MCParticle>();
      for (unsigned int i = 0; i < digitMCRelations.size(); ++i) {
        if (ii < 10) {
          const auto mcParticle = digitMCRelations.object(i);
          idx[ii] = mcParticle->getIndex() - 1;
          wi[ii] = digitMCRelations.weight(i);
          sumHit = sumHit + digitMCRelations.weight(i);
          ii++;
        }
      }
 
      //Re-ordering based on contribution
      int y = 0;
      while (y < 10) {
        for (int i = 0; i < 9; i++) {
          if (((idx[i]) > -1) && ((idx[i + 1]) > -1)) {
            if (wi[i] < wi[i + 1]) {
              int temp = idx[i];
              idx[i] = idx[i + 1];
              idx[i + 1] = temp;
 
              double wtemp = wi[i];
              wi[i] = wi[i + 1];
              wi[i + 1] = wtemp;
            }
          }
        }
        y++;
      }
 
      m_eclCalDigitToBkgWeight->push_back(aECLCalDigits->getEnergy() - sumHit);
      m_eclCalDigitSimHitSum->push_back(sumHit);
      if (idx[0] > -1) {
        m_eclCalDigitToMCWeight1->push_back(wi[0]);
        m_eclCalDigitToMC1->push_back(idx[0]);
        const MCParticle* amcParticle = m_mcParticles[idx[0]];
        m_eclCalDigitToMC1PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclCalDigitToMCWeight1->push_back(-1);
        m_eclCalDigitToMC1->push_back(-1);
        m_eclCalDigitToMC1PDG->push_back(-1);
      }
      if (idx[1] > -1) {
        m_eclCalDigitToMCWeight2->push_back(wi[1]);
        m_eclCalDigitToMC2->push_back(idx[1]);
        const MCParticle* amcParticle = m_mcParticles[idx[1]];
        m_eclCalDigitToMC2PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclCalDigitToMCWeight2->push_back(-1);
        m_eclCalDigitToMC2->push_back(-1);
        m_eclCalDigitToMC2PDG->push_back(-1);
      }
      if (idx[2] > -1) {
        m_eclCalDigitToMCWeight3->push_back(wi[2]);
        m_eclCalDigitToMC3->push_back(idx[2]);
        const MCParticle* amcParticle = m_mcParticles[idx[2]];
        m_eclCalDigitToMC3PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclCalDigitToMCWeight3->push_back(-1);
        m_eclCalDigitToMC3->push_back(-1);
        m_eclCalDigitToMC3PDG->push_back(-1);
      }
      if (idx[3] > -1) {
        m_eclCalDigitToMCWeight4->push_back(wi[3]);
        m_eclCalDigitToMC4->push_back(idx[3]);
        const MCParticle* amcParticle = m_mcParticles[idx[3]];
        m_eclCalDigitToMC4PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclCalDigitToMCWeight4->push_back(-1);
        m_eclCalDigitToMC4->push_back(-1);
        m_eclCalDigitToMC4PDG->push_back(-1);
      }
      if (idx[4] > -1) {
        m_eclCalDigitToMCWeight5->push_back(wi[4]);
        m_eclCalDigitToMC5->push_back(idx[4]);
        const MCParticle* amcParticle = m_mcParticles[idx[4]];
        m_eclCalDigitToMC5PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclCalDigitToMCWeight5->push_back(-1);
        m_eclCalDigitToMC5->push_back(-1);
        m_eclCalDigitToMC5PDG->push_back(-1);
      }
 
      if (aECLCalDigits->getRelated<ECLShower>() != (nullptr)) {
        const ECLShower* shower_caldigit = aECLCalDigits->getRelated<ECLShower>();
        m_eclCalDigitToShower->push_back(shower_caldigit->getArrayIndex());
      } else
        m_eclCalDigitToShower->push_back(-1);
 
      if (aECLCalDigits->getRelated<ECLConnectedRegion>() != (nullptr)) {
        const ECLConnectedRegion* cr_caldigit = aECLCalDigits->getRelated<ECLConnectedRegion>();
        m_eclCalDigitToCR->push_back(cr_caldigit->getCRId());
      } else
        m_eclCalDigitToCR->push_back(-1);
 
      if (aECLCalDigits->getRelated<ECLLocalMaximum>() != (nullptr)) {
        const ECLLocalMaximum* lm_caldigit = aECLCalDigits->getRelated<ECLLocalMaximum>();
        m_eclCalDigitToLM->push_back(lm_caldigit->getLMId());
      } else
        m_eclCalDigitToLM->push_back(-1);
 
    }
  }
 
  //CR
  int CRmultip = m_eclConnectedRegions.getEntries();
  for (int i = 0; i < CRmultip; i++) {
    const ECLConnectedRegion* aECLCR = m_eclConnectedRegions[i];
    m_eclCRIdx->push_back(aECLCR->getCRId());
    m_eclCRIsTrack->push_back(aECLCR->isTrack());
    m_eclCRLikelihoodMIPNGamma->push_back(aECLCR->getLikelihoodMIPNGamma());
    m_eclCRLikelihoodChargedHadron->push_back(aECLCR->getLikelihoodChargedHadron());
    m_eclCRLikelihoodElectronNGamma->push_back(aECLCR->getLikelihoodElectronNGamma());
    m_eclCRLikelihoodNGamma->push_back(aECLCR->getLikelihoodNGamma());
    m_eclCRLikelihoodNeutralHadron->push_back(aECLCR->getLikelihoodNeutralHadron());
    m_eclCRLikelihoodMergedPi0->push_back(aECLCR->getLikelihoodMergedPi0());
  }
 
  if (m_doHits == 1) {
 
    //SIM HITS
    m_eclSimHitMultip = m_eclSimHits.getEntries();
    for (int isimhits = 0; isimhits < m_eclSimHits.getEntries() ; isimhits++) {
      ECLSimHit* aECLSimHits = m_eclSimHits[isimhits];
 
      m_eclSimHitIdx->push_back(isimhits);
      m_eclSimHitCellId->push_back(aECLSimHits->getCellId());
      m_eclSimHitPdg->push_back(aECLSimHits->getPDGCode());
      m_eclSimHitEnergyDep->push_back(aECLSimHits->getEnergyDep());
      m_eclSimHitFlightTime->push_back(aECLSimHits->getFlightTime());
      m_eclSimHitX->push_back(aECLSimHits->getPosition().x());
      m_eclSimHitY->push_back(aECLSimHits->getPosition().y());
      m_eclSimHitZ->push_back(aECLSimHits->getPosition().z());
      m_eclSimHitPx->push_back(aECLSimHits->getMomentum().x());
      m_eclSimHitPy->push_back(aECLSimHits->getMomentum().y());
      m_eclSimHitPz->push_back(aECLSimHits->getMomentum().z());
 
      if (aECLSimHits->getRelated<MCParticle>() != (nullptr)) {
        const MCParticle* mc_simhit = aECLSimHits->getRelated<MCParticle>();
        m_eclSimHitToMC->push_back(mc_simhit->getArrayIndex());
      } else
        m_eclSimHitToMC->push_back(-1);
    }
 
    //HITS
    m_eclHitMultip = m_eclHits.getEntries();
    for (int ihits = 0; ihits < m_eclHits.getEntries() ; ihits++) {
      ECLHit* aECLHits = m_eclHits[ihits];
      m_eclHitIdx->push_back(ihits);
      m_eclHitCellId->push_back(aECLHits->getCellId());
      m_eclHitEnergyDep->push_back(aECLHits->getEnergyDep());
      m_eclHitTimeAve->push_back(aECLHits->getTimeAve());
 
      if (aECLHits->getRelated<ECLDigit>() != (nullptr)) {
        const ECLDigit* hit_digit = aECLHits->getRelated<ECLDigit>();
        m_eclHitToDigit->push_back(hit_digit->getArrayIndex());
        m_eclHitToDigitAmp->push_back(hit_digit->getAmp());
      } else {
        m_eclHitToDigit->push_back(-1);
        m_eclHitToDigitAmp->push_back(-1);
      }
 
      if (m_doPureCsI == true) {
        if (aECLHits->getRelated<ECLDigit>(eclPureDigitArrayName()) != (nullptr)) {
          const ECLDigit* hit_pdigit = aECLHits->getRelated<ECLDigit>(eclPureDigitArrayName());
          m_eclHitToPureDigit->push_back(hit_pdigit->getArrayIndex());
          m_eclHitToPureDigitAmp->push_back(hit_pdigit->getAmp());
        } else {
          m_eclHitToPureDigit->push_back(-1);
          m_eclHitToPureDigitAmp->push_back(-1);
        }
      }
 
      if (aECLHits->getRelated<MCParticle>() != (nullptr)) {
        const MCParticle* mc_hit = aECLHits->getRelated<MCParticle>();
        m_eclHitToMC->push_back(mc_hit->getArrayIndex());
      } else
        m_eclHitToMC->push_back(-1);
    }
  }
 
  //LM
  m_eclLMMultip = m_eclLocalMaximums.getEntries();
  for (unsigned int ilms = 0; ilms < (unsigned int)m_eclLocalMaximums.getEntries() ; ilms++) {
    const ECLLocalMaximum* aECLLMs = m_eclLocalMaximums[ilms];
    m_eclLMId->push_back(aECLLMs->getLMId());
    m_eclLMType->push_back(aECLLMs->getType());
    m_eclLMCellId->push_back(aECLLMs->getCellId());
  }
 
  //CLUSTERS
  m_eclClusterMultip = m_eclClusters.getEntries();
  for (unsigned int iclusters = 0; iclusters < (unsigned int)m_eclClusters.getEntries() ; iclusters++) {
    ECLCluster* aECLClusters = m_eclClusters[iclusters];
    m_eclClusterIdx->push_back(iclusters);
 
    double clusterE = 0.0;
    if (aECLClusters->hasHypothesis(ECLCluster::EHypothesisBit::c_nPhotons)) clusterE = aECLClusters->getEnergy(
            ECLCluster::EHypothesisBit::c_nPhotons);
    else clusterE = aECLClusters->getEnergy(ECLCluster::EHypothesisBit::c_neutralHadron);
    m_eclClusterEnergy->push_back(clusterE);
 
    m_eclClusterEnergyError->push_back(aECLClusters->getUncertaintyEnergy());
    m_eclClusterTheta->push_back(aECLClusters->getTheta());
    m_eclClusterThetaError->push_back(aECLClusters->getUncertaintyTheta());
    m_eclClusterPhi->push_back(aECLClusters->getPhi());
    m_eclClusterPhiError->push_back(aECLClusters->getUncertaintyPhi());
    m_eclClusterR->push_back(aECLClusters->getR());
    m_eclClusterEnergyDepSum->push_back(aECLClusters->getEnergyRaw());
    m_eclClusterTiming->push_back(aECLClusters->getTime());
    m_eclClusterTimingError->push_back(aECLClusters->getDeltaTime99());
    m_eclClusterE9oE21->push_back(aECLClusters->getE9oE21());
    m_eclClusterHighestE->push_back(aECLClusters->getEnergyHighestCrystal());
    m_eclClusterCellId->push_back(aECLClusters->getMaxECellId());
    m_eclClusterNofCrystals->push_back(aECLClusters->getNumberOfCrystals());
    m_eclClusterCrystalHealth->push_back(aECLClusters->getStatus());
 
    m_eclClusterIsTrack->push_back(aECLClusters->isTrack());
    m_eclClusterClosestTrackDist->push_back(aECLClusters->getMinTrkDistance());
    m_eclClusterDeltaL->push_back(aECLClusters->getDeltaL());
 
    m_eclClusterAbsZernike40->push_back(aECLClusters->getAbsZernike40());
    m_eclClusterAbsZernike51->push_back(aECLClusters->getAbsZernike51());
    m_eclClusterZernikeMVA->push_back(aECLClusters->getZernikeMVA());
    m_eclClusterE1oE9->push_back(aECLClusters->getE1oE9());
    m_eclClusterSecondMoment->push_back(aECLClusters->getSecondMoment());
    m_eclClusterLAT->push_back(aECLClusters->getLAT());
    m_eclClusterDeltaTime99->push_back(aECLClusters->getDeltaTime99());
    m_eclClusterDetectorRegion->push_back(aECLClusters->getDetectorRegion());
    m_eclClusterHasNPhotonHypothesis->push_back(aECLClusters->hasHypothesis(ECLCluster::EHypothesisBit::c_nPhotons));
    m_eclClusterHasNeutralHadronHypothesis->push_back(aECLClusters->hasHypothesis(ECLCluster::EHypothesisBit::c_neutralHadron));
 
    if (aECLClusters->getRelated<ECLShower>() != (nullptr)) {
      const ECLShower* shower_cluster = aECLClusters->getRelated<ECLShower>();
      m_eclClusterToShower->push_back(shower_cluster->getArrayIndex());
    } else
      m_eclClusterToShower->push_back(-1);
 
    //Dump MC Info - Multiple Matching
    double sumHit = 0;
    int idx[10];
    for (int i = 0; i < 10; i++)
      idx[i] = -1;
 
    double wi[10];
    for (int i = 0; i < 10; i++)
      wi[i] = -1;
 
    int ii = 0;
    sumHit = 0;
 
    auto clusterMCRelations = aECLClusters->getRelationsTo<MCParticle>();
    for (unsigned int i = 0; i < clusterMCRelations.size(); ++i) {
      if (ii < 10) {
        const auto mcParticle = clusterMCRelations.object(i);
        idx[ii] = mcParticle->getIndex() - 1;
        wi[ii] = clusterMCRelations.weight(i);
        sumHit = sumHit + clusterMCRelations.weight(i);
        ii++;
      }
    }
 
    //Re-ordering based on contribution
    int y = 0;
    while (y < 10) {
      for (int i = 0; i < 9; i++) {
        if (((idx[i]) > -1) && ((idx[i + 1]) > -1)) {
          if (wi[i] < wi[i + 1]) {
            int temp = idx[i];
            idx[i] = idx[i + 1];
            idx[i + 1] = temp;
 
            double wtemp = wi[i];
            wi[i] = wi[i + 1];
            wi[i + 1] = wtemp;
          }
        }
      }
      y++;
    }
 
    m_eclClusterToBkgWeight->push_back(clusterE - sumHit);
    m_eclClusterSimHitSum->push_back(sumHit);
    if (idx[0] > -1) {
      m_eclClusterToMCWeight1->push_back(wi[0]);
      m_eclClusterToMC1->push_back(idx[0]);
      const MCParticle* amcParticle = m_mcParticles[idx[0]];
      m_eclClusterToMC1PDG->push_back(amcParticle->getPDG());
    } else {
      m_eclClusterToMCWeight1->push_back(-1);
      m_eclClusterToMC1->push_back(-1);
      m_eclClusterToMC1PDG->push_back(-1);
    }
    if (idx[1] > -1) {
      m_eclClusterToMCWeight2->push_back(wi[1]);
      m_eclClusterToMC2->push_back(idx[1]);
      const MCParticle* amcParticle = m_mcParticles[idx[1]];
      m_eclClusterToMC2PDG->push_back(amcParticle->getPDG());
    } else {
      m_eclClusterToMCWeight2->push_back(-1);
      m_eclClusterToMC2->push_back(-1);
      m_eclClusterToMC2PDG->push_back(-1);
    }
    if (idx[2] > -1) {
      m_eclClusterToMCWeight3->push_back(wi[2]);
      m_eclClusterToMC3->push_back(idx[2]);
      const MCParticle* amcParticle = m_mcParticles[idx[2]];
      m_eclClusterToMC3PDG->push_back(amcParticle->getPDG());
    } else {
      m_eclClusterToMCWeight3->push_back(-1);
      m_eclClusterToMC3->push_back(-1);
      m_eclClusterToMC3PDG->push_back(-1);
    }
    if (idx[3] > -1) {
      m_eclClusterToMCWeight4->push_back(wi[3]);
      m_eclClusterToMC4->push_back(idx[3]);
      const MCParticle* amcParticle = m_mcParticles[idx[3]];
      m_eclClusterToMC4PDG->push_back(amcParticle->getPDG());
    } else {
      m_eclClusterToMCWeight4->push_back(-1);
      m_eclClusterToMC4->push_back(-1);
      m_eclClusterToMC4PDG->push_back(-1);
    }
    if (idx[4] > -1) {
      m_eclClusterToMCWeight5->push_back(wi[4]);
      m_eclClusterToMC5->push_back(idx[4]);
      const MCParticle* amcParticle = m_mcParticles[idx[4]];
      m_eclClusterToMC5PDG->push_back(amcParticle->getPDG());
    } else {
      m_eclClusterToMCWeight5->push_back(-1);
      m_eclClusterToMC5->push_back(-1);
      m_eclClusterToMC5PDG->push_back(-1);
    }
 
  }
 
  if (m_doPureCsI == true) {
 
    m_eclPureDigitMultip = m_eclPureDigits.getEntries();
    for (int idigits = 0; idigits < m_eclPureDigits.getEntries() ; idigits++) {
      ECLDigit* aECLPureDigits = m_eclPureDigits[idigits];
 
      m_eclPureDigitIdx->push_back(idigits);
      m_eclPureDigitCellId->push_back(aECLPureDigits->getCellId());
      m_eclPureDigitAmp->push_back(aECLPureDigits->getAmp());
      m_eclPureDigitTimeFit->push_back(aECLPureDigits->getTimeFit());
      m_eclPureDigitFitQuality->push_back(aECLPureDigits->getQuality());
 
      if (aECLPureDigits->getRelated<MCParticle>() != (nullptr)) {
        const MCParticle* mc_digit = aECLPureDigits->getRelated<MCParticle>();
        m_eclPureDigitToMC->push_back(mc_digit->getArrayIndex());
      } else
        m_eclPureDigitToMC->push_back(-1);
    }
 
    //PURE CAL DIGITS
    m_eclPureCalDigitMultip = m_eclPureCalDigits.getEntries();
    for (uint icaldigits = 0; icaldigits < (uint)m_eclPureCalDigits.getEntries() ; icaldigits++) {
      ECLCalDigit* aECLPureCalDigits = m_eclPureCalDigits[icaldigits];
 
      m_eclPureCalDigitIdx->push_back(icaldigits);
      m_eclPureCalDigitCellId->push_back(aECLPureCalDigits->getCellId());
      m_eclPureCalDigitAmp->push_back(aECLPureCalDigits->getEnergy());
      m_eclPureCalDigitTimeFit->push_back(aECLPureCalDigits->getTime());
      m_eclPureCalDigitFitQuality->push_back(aECLPureCalDigits->isFailedFit());
 
      double sumHit = 0;
      int idx[10];
      for (int i = 0; i < 10; i++)
        idx[i] = -1;
 
      double wi[10];
      for (int i = 0; i < 10; i++)
        wi[i] = -1;
 
      int ii = 0;
      sumHit = 0;
 
      auto digitMCRelations = aECLPureCalDigits->getRelationsTo<MCParticle>();
      for (unsigned int i = 0; i < digitMCRelations.size(); ++i) {
        if (ii < 10) {
          const auto mcParticle = digitMCRelations.object(i);
          idx[ii] = mcParticle->getIndex() - 1;
          wi[ii] = digitMCRelations.weight(i);
          sumHit = sumHit + digitMCRelations.weight(i);
          ii++;
        }
      }
 
      //Re-ordering based on contribution
      int y = 0;
      while (y < 10) {
        for (int i = 0; i < 9; i++) {
          if (((idx[i]) > -1) && ((idx[i + 1]) > -1)) {
            if (wi[i] < wi[i + 1]) {
              int temp = idx[i];
              idx[i] = idx[i + 1];
              idx[i + 1] = temp;
 
              double wtemp = wi[i];
              wi[i] = wi[i + 1];
              wi[i + 1] = wtemp;
            }
          }
        }
        y++;
      }
 
      m_eclPureCalDigitToBkgWeight->push_back(aECLPureCalDigits->getEnergy() - sumHit);
      m_eclPureCalDigitSimHitSum->push_back(sumHit);
      if (idx[0] > -1) {
        m_eclPureCalDigitToMCWeight1->push_back(wi[0]);
        m_eclPureCalDigitToMC1->push_back(idx[0]);
        const MCParticle* amcParticle = m_mcParticles[idx[0]];
        m_eclPureCalDigitToMC1PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclPureCalDigitToMCWeight1->push_back(-1);
        m_eclPureCalDigitToMC1->push_back(-1);
        m_eclPureCalDigitToMC1PDG->push_back(-1);
      }
      if (idx[1] > -1) {
        m_eclPureCalDigitToMCWeight2->push_back(wi[1]);
        m_eclPureCalDigitToMC2->push_back(idx[1]);
        const MCParticle* amcParticle = m_mcParticles[idx[1]];
        m_eclPureCalDigitToMC2PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclPureCalDigitToMCWeight2->push_back(-1);
        m_eclPureCalDigitToMC2->push_back(-1);
        m_eclPureCalDigitToMC2PDG->push_back(-1);
      }
      if (idx[2] > -1) {
        m_eclPureCalDigitToMCWeight3->push_back(wi[2]);
        m_eclPureCalDigitToMC3->push_back(idx[2]);
        const MCParticle* amcParticle = m_mcParticles[idx[2]];
        m_eclPureCalDigitToMC3PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclPureCalDigitToMCWeight3->push_back(-1);
        m_eclPureCalDigitToMC3->push_back(-1);
        m_eclPureCalDigitToMC3PDG->push_back(-1);
      }
      if (idx[3] > -1) {
        m_eclPureCalDigitToMCWeight4->push_back(wi[3]);
        m_eclPureCalDigitToMC4->push_back(idx[3]);
        const MCParticle* amcParticle = m_mcParticles[idx[3]];
        m_eclPureCalDigitToMC4PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclPureCalDigitToMCWeight4->push_back(-1);
        m_eclPureCalDigitToMC4->push_back(-1);
        m_eclPureCalDigitToMC4PDG->push_back(-1);
      }
      if (idx[4] > -1) {
        m_eclPureCalDigitToMCWeight5->push_back(wi[4]);
        m_eclPureCalDigitToMC5->push_back(idx[4]);
        const MCParticle* amcParticle = m_mcParticles[idx[4]];
        m_eclPureCalDigitToMC5PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclPureCalDigitToMCWeight5->push_back(-1);
        m_eclPureCalDigitToMC5->push_back(-1);
        m_eclPureCalDigitToMC5PDG->push_back(-1);
      }
 
      if (aECLPureCalDigits->getRelated<ECLShower>() != (nullptr)) {
        const ECLShower* shower_caldigit = aECLPureCalDigits->getRelated<ECLShower>();
        m_eclPureCalDigitToShower->push_back(shower_caldigit->getArrayIndex());
      } else
        m_eclPureCalDigitToShower->push_back(-1);
 
      if (aECLPureCalDigits->getRelated<ECLConnectedRegion>(eclPureConnectedRegionArrayName()) != (nullptr)) {
        const ECLConnectedRegion* cr_caldigit = aECLPureCalDigits->getRelated<ECLConnectedRegion>(eclPureConnectedRegionArrayName());
        m_eclPureCalDigitToCR->push_back(cr_caldigit->getCRId());
      } else
        m_eclPureCalDigitToCR->push_back(-1);
 
      if (aECLPureCalDigits->getRelated<ECLLocalMaximum>(eclPureLocalMaximumArrayName()) != (nullptr)) {
        const ECLLocalMaximum* lm_caldigit = aECLPureCalDigits->getRelated<ECLLocalMaximum>(eclPureLocalMaximumArrayName());
        m_eclPureCalDigitToLM->push_back(lm_caldigit->getLMId());
      } else
        m_eclPureCalDigitToLM->push_back(-1);
 
    }
 
    //CR
    int PureCRmultip = m_eclPureConnectedRegions.getEntries();
    for (int i = 0; i < PureCRmultip; i++) {
      const ECLConnectedRegion* aECLPureCR = m_eclPureConnectedRegions[i];
      m_eclPureCRIdx->push_back(aECLPureCR->getCRId());
      m_eclPureCRIsTrack->push_back(aECLPureCR->isTrack());
      m_eclPureCRLikelihoodMIPNGamma->push_back(aECLPureCR->getLikelihoodMIPNGamma());
      m_eclPureCRLikelihoodChargedHadron->push_back(aECLPureCR->getLikelihoodChargedHadron());
      m_eclPureCRLikelihoodElectronNGamma->push_back(aECLPureCR->getLikelihoodElectronNGamma());
      m_eclPureCRLikelihoodNGamma->push_back(aECLPureCR->getLikelihoodNGamma());
      m_eclPureCRLikelihoodNeutralHadron->push_back(aECLPureCR->getLikelihoodNeutralHadron());
      m_eclPureCRLikelihoodMergedPi0->push_back(aECLPureCR->getLikelihoodMergedPi0());
    }
 
    //LM
    m_eclPureLMMultip = m_eclPureLocalMaximums.getEntries();
    for (unsigned int pure_ilms = 0; pure_ilms < (unsigned int)m_eclPureLocalMaximums.getEntries() ; pure_ilms++) {
      const ECLLocalMaximum* aECLLMs = m_eclPureLocalMaximums[pure_ilms];
      m_eclPureLMId->push_back(aECLLMs->getLMId());
      m_eclPureLMType->push_back(aECLLMs->getType());
      m_eclPureLMCellId->push_back(aECLLMs->getCellId());
    }
 
    m_eclPureClusterMultip = m_eclPureClusters.getEntries();
    for (unsigned int iclusters = 0; iclusters < (unsigned int)m_eclPureClusters.getEntries() ; iclusters++) {
      ECLCluster* aECLClusters = m_eclPureClusters[iclusters];
      m_eclPureClusterIdx->push_back(iclusters);
 
      double clusterE = 0.0;
      if (aECLClusters->hasHypothesis(ECLCluster::EHypothesisBit::c_nPhotons)) clusterE = aECLClusters->getEnergy(
              ECLCluster::EHypothesisBit::c_nPhotons);
      else clusterE = aECLClusters->getEnergy(ECLCluster::EHypothesisBit::c_neutralHadron);
      m_eclPureClusterEnergy->push_back(clusterE);
 
      m_eclPureClusterEnergyError->push_back(aECLClusters->getUncertaintyEnergy());
      m_eclPureClusterTheta->push_back(aECLClusters->getTheta());
      m_eclPureClusterThetaError->push_back(aECLClusters->getUncertaintyTheta());
      m_eclPureClusterPhi->push_back(aECLClusters->getPhi());
      m_eclPureClusterPhiError->push_back(aECLClusters->getUncertaintyPhi());
      m_eclPureClusterR->push_back(aECLClusters->getR());
      m_eclPureClusterEnergyDepSum->push_back(aECLClusters->getEnergyRaw());
      m_eclPureClusterTiming->push_back(aECLClusters->getTime());
      m_eclPureClusterTimingError->push_back(aECLClusters->getDeltaTime99());
      m_eclPureClusterE9oE21->push_back(aECLClusters->getE9oE21());
      m_eclPureClusterHighestE->push_back(aECLClusters->getEnergyHighestCrystal());
      m_eclPureClusterCellId->push_back(aECLClusters->getMaxECellId());
      m_eclPureClusterLat->push_back(aECLClusters->getLAT());
      m_eclPureClusterNofCrystals->push_back(aECLClusters->getNumberOfCrystals());
      m_eclPureClusterCrystalHealth->push_back(aECLClusters->getStatus());
 
      m_eclPureClusterClosestTrackDist->push_back(aECLClusters->getMinTrkDistance());
      m_eclPureClusterAbsZernike40->push_back(aECLClusters->getAbsZernike40());
      m_eclPureClusterAbsZernike51->push_back(aECLClusters->getAbsZernike51());
      m_eclPureClusterZernikeMVA->push_back(aECLClusters->getZernikeMVA());
      m_eclPureClusterSecondMoment->push_back(aECLClusters->getSecondMoment());
 
      m_eclPureClusterIsTrack->push_back(aECLClusters->isTrack());
      m_eclPureClusterDeltaL->push_back(aECLClusters->getDeltaL());
 
      m_eclPureClusterE1oE9->push_back(aECLClusters->getE1oE9());
      m_eclPureClusterDeltaTime99->push_back(aECLClusters->getDeltaTime99());
      m_eclPureClusterDetectorRegion->push_back(aECLClusters->getDetectorRegion());
      m_eclPureClusterHasNPhotonHypothesis->push_back(aECLClusters->hasHypothesis(ECLCluster::EHypothesisBit::c_nPhotons));
      m_eclPureClusterHasNeutralHadronHypothesis->push_back(aECLClusters->hasHypothesis(ECLCluster::EHypothesisBit::c_neutralHadron));
 
      //Dump MC Info - Multiple Matching
      double sumHit = 0;
      int idx[10];
      for (int i = 0; i < 10; i++)
        idx[i] = -1;
 
      double wi[10];
      for (int i = 0; i < 10; i++)
        wi[i] = -1;
 
      int ii = 0;
      sumHit = 0;
 
      auto clusterMCRelations = aECLClusters->getRelationsTo<MCParticle>();
      for (unsigned int i = 0; i < clusterMCRelations.size(); ++i) {
        if (ii < 10) {
          const auto mcParticle = clusterMCRelations.object(i);
          idx[ii] = mcParticle->getIndex() - 1;
          wi[ii] = clusterMCRelations.weight(i);
          sumHit = sumHit + clusterMCRelations.weight(i);
          ii++;
        }
      }
 
      //Re-ordering based on contribution
      int y = 0;
      while (y < 10) {
        for (int i = 0; i < 9; i++) {
          if (((idx[i]) > -1) && ((idx[i + 1]) > -1)) {
            if (wi[i] < wi[i + 1]) {
              int temp = idx[i];
              idx[i] = idx[i + 1];
              idx[i + 1] = temp;
 
              double wtemp = wi[i];
              wi[i] = wi[i + 1];
              wi[i + 1] = wtemp;
            }
          }
        }
        y++;
      }
 
      m_eclPureClusterToBkgWeight->push_back(clusterE - sumHit);
      if (idx[0] > -1) {
        m_eclPureClusterToMCWeight1->push_back(wi[0]);
        m_eclPureClusterToMC1->push_back(idx[0]);
        const MCParticle* amcParticle = m_mcParticles[idx[0]];
        m_eclPureClusterToMC1PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclPureClusterToMCWeight1->push_back(-1);
        m_eclPureClusterToMC1->push_back(-1);
        m_eclPureClusterToMC1PDG->push_back(-1);
      }
      if (idx[1] > -1) {
        m_eclPureClusterToMCWeight2->push_back(wi[1]);
        m_eclPureClusterToMC2->push_back(idx[1]);
        const MCParticle* amcParticle = m_mcParticles[idx[1]];
        m_eclPureClusterToMC2PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclPureClusterToMCWeight2->push_back(-1);
        m_eclPureClusterToMC2->push_back(-1);
        m_eclPureClusterToMC2PDG->push_back(-1);
      }
      if (idx[2] > -1) {
        m_eclPureClusterToMCWeight3->push_back(wi[2]);
        m_eclPureClusterToMC3->push_back(idx[2]);
        const MCParticle* amcParticle = m_mcParticles[idx[2]];
        m_eclPureClusterToMC3PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclPureClusterToMCWeight3->push_back(-1);
        m_eclPureClusterToMC3->push_back(-1);
        m_eclPureClusterToMC3PDG->push_back(-1);
      }
      if (idx[3] > -1) {
        m_eclPureClusterToMCWeight4->push_back(wi[3]);
        m_eclPureClusterToMC4->push_back(idx[3]);
        const MCParticle* amcParticle = m_mcParticles[idx[3]];
        m_eclPureClusterToMC4PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclPureClusterToMCWeight4->push_back(-1);
        m_eclPureClusterToMC4->push_back(-1);
        m_eclPureClusterToMC4PDG->push_back(-1);
      }
      if (idx[4] > -1) {
        m_eclPureClusterToMCWeight5->push_back(wi[4]);
        m_eclPureClusterToMC5->push_back(idx[4]);
        const MCParticle* amcParticle = m_mcParticles[idx[4]];
        m_eclPureClusterToMC5PDG->push_back(amcParticle->getPDG());
      } else {
        m_eclPureClusterToMCWeight5->push_back(-1);
        m_eclPureClusterToMC5->push_back(-1);
        m_eclPureClusterToMC5PDG->push_back(-1);
      }
 
    }
  }
 
  m_eclShowerMultip = m_eclShowers.getEntries();
  for (uint ishowers = 0; ishowers < (uint)m_eclShowers.getEntries() ; ishowers++) {
    ECLShower* aECLShowers = m_eclShowers[ishowers];
    m_eclShowerIdx->push_back(ishowers);
    m_eclShowerEnergy->push_back(aECLShowers->getEnergy());
    m_eclShowerTheta->push_back(aECLShowers->getTheta());
    m_eclShowerPhi->push_back(aECLShowers->getPhi());
    m_eclShowerR->push_back(aECLShowers->getR());
    m_eclShowerNHits->push_back(aECLShowers->getNumberOfCrystals());
    m_eclShowerE9oE21->push_back(aECLShowers->getE9oE21());
    m_eclShowerUncEnergy->push_back(aECLShowers->getEnergyRaw());
    m_eclShowerTime->push_back(aECLShowers->getTime());
    m_eclShowerT99->push_back(aECLShowers->getDeltaTime99());
    m_eclShowerConnectedRegionId->push_back(aECLShowers->getConnectedRegionId());
    m_eclShowerHypothesisId->push_back(aECLShowers->getHypothesisId());
    m_eclShowerCentralCellId->push_back(aECLShowers->getCentralCellId());
    m_eclShowerEnergyError->push_back(aECLShowers->getUncertaintyEnergy());
    m_eclShowerThetaError->push_back(aECLShowers->getUncertaintyTheta());
    m_eclShowerPhiError->push_back(aECLShowers->getUncertaintyPhi());
    m_eclShowerTimeResolution->push_back(aECLShowers->getDeltaTime99());
    m_eclShowerHighestEnergy->push_back(aECLShowers->getEnergyHighestCrystal());
    m_eclShowerLateralEnergy->push_back(aECLShowers->getLateralEnergy());
    m_eclShowerMinTrkDistance->push_back(aECLShowers->getMinTrkDistance());
    m_eclShowerTrkDepth->push_back(aECLShowers->getTrkDepth());
    m_eclShowerShowerDepth->push_back(aECLShowers->getShowerDepth());
    m_eclShowerAbsZernike40->push_back(aECLShowers->getAbsZernikeMoment(4, 0));
    m_eclShowerAbsZernike51->push_back(aECLShowers->getAbsZernikeMoment(5, 1));
    m_eclShowerZernikeMVA->push_back(aECLShowers->getZernikeMVA());
    m_eclShowerSecondMoment->push_back(aECLShowers->getSecondMoment());
    m_eclShowerE1oE9->push_back(aECLShowers->getE1oE9());
    m_eclShowerIsTrack->push_back(aECLShowers->getIsTrack());
    m_eclShowerNumberOfCrystalsForEnergy->push_back(aECLShowers->getNumberOfCrystalsForEnergy());
 
    double fe = 0.;
    double se = 0.;
 
    auto showerDigitRelations = aECLShowers->getRelationsTo<ECLCalDigit>();
    for (unsigned int i = 0; i < showerDigitRelations.size(); ++i) {
      const auto aECLCalDigits = showerDigitRelations.object(i);
      if (aECLCalDigits->getEnergy() > fe) {
        se = fe;
        fe = aECLCalDigits->getEnergy();
      }
    }
    if (fe > 0 && se > 0)
      m_eclShowerHighestE1mE2->push_back(fe - se);
    else
      m_eclShowerHighestE1mE2->push_back(-1);
 
    int lm1[5] = { -1, -1, -1, -1, -1};
 
    auto showerLMRelations = aECLShowers->getRelationsTo<ECLLocalMaximum>();
    for (unsigned int i = 0; i < showerLMRelations.size(); ++i) {
      const auto aECLLM = showerLMRelations.object(i);
      lm1[i] = aECLLM->getLMId();
    }
    m_eclShowerToLM1->push_back(lm1[0]);
    m_eclShowerToLM2->push_back(lm1[1]);
    m_eclShowerToLM3->push_back(lm1[2]);
    m_eclShowerToLM4->push_back(lm1[3]);
    m_eclShowerToLM5->push_back(lm1[4]);
 
    double sumHit = 0;
 
    int idx[10];
    for (int i = 0; i < 10; i++)
      idx[i] = -1;
 
    double wi[10];
    for (int i = 0; i < 10; i++)
      wi[i] = -1;
 
    int ii = 0;
    sumHit = 0;
 
    auto showerMCRelations = aECLShowers->getRelationsTo<MCParticle>();
    for (unsigned int i = 0; i < showerMCRelations.size(); ++i) {
      if (ii < 10) {
        const auto mcParticle = showerMCRelations.object(i);
        idx[ii] = mcParticle->getIndex() - 1;
        wi[ii] = showerMCRelations.weight(i);
        sumHit = sumHit + showerMCRelations.weight(i);
        ii++;
      }
    }
 
    //Re-ordering based on contribution
    int y = 0;
    while (y < 10) {
      for (int i = 0; i < 9; i++) {
        if (((idx[i]) > -1) && ((idx[i + 1]) > -1)) {
          if (wi[i] < wi[i + 1]) {
            int temp = idx[i];
            idx[i] = idx[i + 1];
            idx[i + 1] = temp;
 
            double wtemp = wi[i];
            wi[i] = wi[i + 1];
            wi[i + 1] = wtemp;
          }
        }
      }
      y++;
    }
 
    int no_Primary = 1;
 
    for (unsigned int i = 0; i < showerMCRelations.size(); ++i) {
      const auto mcParticle = showerMCRelations.object(i);
      if (mcParticle->getSecondaryPhysicsProcess() == 0 && mcParticle->getPDG() == Const::Klong.getPDGCode()) {
        double vtxx = mcParticle->getDecayVertex().X();
        double vtxy = mcParticle->getDecayVertex().Y();
        double vtxz = mcParticle->getDecayVertex().Z();
        if ((TMath::Sqrt(vtxx * vtxx + vtxy * vtxy) > 118) || (vtxz > 196.16) || (vtxz < -102.16))
          no_Primary = 0;
      } else if (mcParticle->getSecondaryPhysicsProcess() != 0 && mcParticle->getMother()->getPDG() == Const::Klong.getPDGCode()) {
        double vtxx = mcParticle->getProductionVertex().X();
        double vtxy = mcParticle->getProductionVertex().Y();
        double vtxz = mcParticle->getProductionVertex().Z();
        if ((TMath::Sqrt(vtxx * vtxx + vtxy * vtxy) > 118) || (vtxz > 196.16) || (vtxz < -102.16))
          no_Primary = 0;
      }
    }
 
    if (no_Primary == 0)
      m_eclShowerMCVtxInEcl->push_back(1);
    else
      m_eclShowerMCVtxInEcl->push_back(0);
 
    double no_fMatch = 0;
    double no_fFMatch = 0;
 
    for (unsigned int i = 0; i < showerMCRelations.size(); ++i) {
      const auto mcParticle = showerMCRelations.object(i);
      if (mcParticle->getSecondaryPhysicsProcess() == 0) {
        double vtxx = mcParticle->getDecayVertex().X();
        double vtxy = mcParticle->getDecayVertex().Y();
        double vtxz = mcParticle->getDecayVertex().Z();
        double px = mcParticle->getMomentum().X();
        double py = mcParticle->getMomentum().Y();
        double pz = mcParticle->getMomentum().Z();
        double p = TMath::Sqrt(px * px + py * py + pz * pz);
        double pTheta = TMath::ACos(pz / p);
        double pPhi = 10.;
        if (py > 0) {
          if (px > 0)
            pPhi = TMath::ATan(py / px);
          if (px < 0)
            pPhi = TMath::ATan(py / px) + M_PI;
        } else {
          if (px > 0)
            pPhi = TMath::ATan(py / px) ;
          if (px < 0)
            pPhi = TMath::ATan(py / px) - M_PI;
        }
        if ((TMath::Sqrt(vtxx * vtxx + vtxy * vtxy) > 118) || (vtxz > 196.16) || (vtxz < -102.16)) {
          if (TMath::Abs(aECLShowers->getTheta() - pTheta) < 0.05 && TMath::Abs(aECLShowers->getPhi() - pPhi) < 0.05)
            no_fMatch = 1;
          no_fFMatch = 1;
        }
      } else if (mcParticle->getMother()->getPDG() == Const::Klong.getPDGCode()) {
        double vtxx = mcParticle->getProductionVertex().X();
        double vtxy = mcParticle->getProductionVertex().Y();
        double vtxz = mcParticle->getProductionVertex().Z();
        const MCParticle* amcParticle = mcParticle->getMother();
        double px = amcParticle->getMomentum().X();
        double py = amcParticle->getMomentum().Y();
        double pz = amcParticle->getMomentum().Z();
        double p = TMath::Sqrt(px * px + py * py + pz * pz);
        double pTheta = TMath::ACos(pz / p);
        double pPhi = 10.;
        if (py > 0) {
          if (px > 0)
            pPhi = TMath::ATan(py / px);
          if (px < 0)
            pPhi = TMath::ATan(py / px) + M_PI;
        } else {
          if (px > 0)
            pPhi = TMath::ATan(py / px) ;
          if (px < 0)
            pPhi = TMath::ATan(py / px) - M_PI;
        }
        if ((TMath::Sqrt(vtxx * vtxx + vtxy * vtxy) > 118) || (vtxz > 196.16) || (vtxz < -102.16)) {
          if (TMath::Abs(aECLShowers->getTheta() - pTheta) < 0.05 && TMath::Abs(aECLShowers->getPhi() - pPhi) < 0.05)
            no_fMatch = 1;
        }
      }
    }
 
    if (no_fMatch == 1)
      m_eclShowerMCFlightMatch->push_back(1);
    else
      m_eclShowerMCFlightMatch->push_back(0);
 
    if (no_fFMatch == 1)
      m_eclShowerMCFFlightMatch->push_back(1);
    else
      m_eclShowerMCFFlightMatch->push_back(0);
 
    if (idx[0] > -1) {
      m_eclShowerToMCWeight1->push_back(wi[0]);
      m_eclShowerToMC1->push_back(idx[0]);
      const MCParticle* amcParticle = m_mcParticles[idx[0]];
      m_eclShowerToMC1PDG->push_back(amcParticle->getPDG());
      if (amcParticle->getMother() != nullptr) {
        m_eclShowerToMC1Moth->push_back(amcParticle->getMother()->getIndex());
        m_eclShowerToMC1MothPDG->push_back(amcParticle->getMother()->getPDG());
        if (amcParticle->getMother()->getMother() != nullptr) {
          m_eclShowerToMC1GMoth->push_back(amcParticle->getMother()->getMother()->getIndex());
          m_eclShowerToMC1GMothPDG->push_back(amcParticle->getMother()->getMother()->getPDG());
        } else {
          m_eclShowerToMC1GMoth->push_back(-999);
          m_eclShowerToMC1GMothPDG->push_back(-999);
        }
      } else {
        m_eclShowerToMC1Moth->push_back(-999);
        m_eclShowerToMC1MothPDG->push_back(-999);
        m_eclShowerToMC1GMoth->push_back(-999);
        m_eclShowerToMC1GMothPDG->push_back(-999);
      }
    } else {
      m_eclShowerToMCWeight1->push_back(-1);
      m_eclShowerToMC1->push_back(-1);
      m_eclShowerToMC1PDG->push_back(-1);
      m_eclShowerToMC1Moth->push_back(-1);
      m_eclShowerToMC1MothPDG->push_back(-1);
      m_eclShowerToMC1GMoth->push_back(-1);
      m_eclShowerToMC1GMothPDG->push_back(-1);
    }
    if (idx[1] > -1) {
      m_eclShowerToMCWeight2->push_back(wi[1]);
      m_eclShowerToMC2->push_back(idx[1]);
      const MCParticle* amcParticle = m_mcParticles[idx[1]];
      m_eclShowerToMC2PDG->push_back(amcParticle->getPDG());
      if (amcParticle->getMother() != nullptr) {
        m_eclShowerToMC2Moth->push_back(amcParticle->getMother()->getIndex());
        m_eclShowerToMC2MothPDG->push_back(amcParticle->getMother()->getPDG());
        if (amcParticle->getMother()->getMother() != nullptr) {
          m_eclShowerToMC2GMoth->push_back(amcParticle->getMother()->getMother()->getIndex());
          m_eclShowerToMC2GMothPDG->push_back(amcParticle->getMother()->getMother()->getPDG());
        } else {
          m_eclShowerToMC2GMoth->push_back(-999);
          m_eclShowerToMC2GMothPDG->push_back(-999);
        }
      } else {
        m_eclShowerToMC2Moth->push_back(-999);
        m_eclShowerToMC2MothPDG->push_back(-999);
        m_eclShowerToMC2GMoth->push_back(-999);
        m_eclShowerToMC2GMothPDG->push_back(-999);
      };
    } else {
      m_eclShowerToMCWeight2->push_back(-1);
      m_eclShowerToMC2->push_back(-1);
      m_eclShowerToMC2PDG->push_back(-1);
      m_eclShowerToMC2Moth->push_back(-1);
      m_eclShowerToMC2MothPDG->push_back(-1);
      m_eclShowerToMC2GMoth->push_back(-1);
      m_eclShowerToMC2GMothPDG->push_back(-1);
    }
    if (idx[2] > -1) {
      m_eclShowerToMCWeight3->push_back(wi[2]);
      m_eclShowerToMC3->push_back(idx[2]);
      const MCParticle* amcParticle = m_mcParticles[idx[2]];
      m_eclShowerToMC3PDG->push_back(amcParticle->getPDG());
      if (amcParticle->getMother() != nullptr) {
        m_eclShowerToMC3Moth->push_back(amcParticle->getMother()->getIndex());
        m_eclShowerToMC3MothPDG->push_back(amcParticle->getMother()->getPDG());
        if (amcParticle->getMother()->getMother() != nullptr) {
          m_eclShowerToMC3GMoth->push_back(amcParticle->getMother()->getMother()->getIndex());
          m_eclShowerToMC3GMothPDG->push_back(amcParticle->getMother()->getMother()->getPDG());
        } else {
          m_eclShowerToMC3GMoth->push_back(-999);
          m_eclShowerToMC3GMothPDG->push_back(-999);
        }
      } else {
        m_eclShowerToMC3Moth->push_back(-999);
        m_eclShowerToMC3MothPDG->push_back(-999);
        m_eclShowerToMC3GMoth->push_back(-999);
        m_eclShowerToMC3GMothPDG->push_back(-999);
      }
    } else {
      m_eclShowerToMCWeight3->push_back(-1);
      m_eclShowerToMC3->push_back(-1);
      m_eclShowerToMC3PDG->push_back(-1);
      m_eclShowerToMC3Moth->push_back(-1);
      m_eclShowerToMC3MothPDG->push_back(-1);
      m_eclShowerToMC3GMoth->push_back(-1);
      m_eclShowerToMC3GMothPDG->push_back(-1);
    }
    if (idx[3] > -1) {
      m_eclShowerToMCWeight4->push_back(wi[3]);
      m_eclShowerToMC4->push_back(idx[3]);
      const MCParticle* amcParticle = m_mcParticles[idx[3]];
      m_eclShowerToMC4PDG->push_back(amcParticle->getPDG());
      if (amcParticle->getMother() != nullptr) {
        m_eclShowerToMC4Moth->push_back(amcParticle->getMother()->getIndex());
        m_eclShowerToMC4MothPDG->push_back(amcParticle->getMother()->getPDG());
        if (amcParticle->getMother()->getMother() != nullptr) {
          m_eclShowerToMC4GMoth->push_back(amcParticle->getMother()->getMother()->getIndex());
          m_eclShowerToMC4GMothPDG->push_back(amcParticle->getMother()->getMother()->getPDG());
        } else {
          m_eclShowerToMC4GMoth->push_back(-999);
          m_eclShowerToMC4GMothPDG->push_back(-999);
        }
      } else {
        m_eclShowerToMC4Moth->push_back(-999);
        m_eclShowerToMC4MothPDG->push_back(-999);
        m_eclShowerToMC4GMoth->push_back(-999);
        m_eclShowerToMC4GMothPDG->push_back(-999);
      }
    } else {
      m_eclShowerToMCWeight4->push_back(-1);
      m_eclShowerToMC4->push_back(-1);
      m_eclShowerToMC4PDG->push_back(-1);
      m_eclShowerToMC4Moth->push_back(-1);
      m_eclShowerToMC4MothPDG->push_back(-1);
      m_eclShowerToMC4GMoth->push_back(-1);
      m_eclShowerToMC4GMothPDG->push_back(-1);
    }
    if (idx[4] > -1) {
      m_eclShowerToMCWeight5->push_back(wi[4]);
      m_eclShowerToMC5->push_back(idx[4]);
      const MCParticle* amcParticle = m_mcParticles[idx[4]];
      m_eclShowerToMC5PDG->push_back(amcParticle->getPDG());
      if (amcParticle->getMother() != nullptr) {
        m_eclShowerToMC5Moth->push_back(amcParticle->getMother()->getIndex());
        m_eclShowerToMC5MothPDG->push_back(amcParticle->getMother()->getPDG());
        if (amcParticle->getMother()->getMother() != nullptr) {
          m_eclShowerToMC5GMoth->push_back(amcParticle->getMother()->getMother()->getIndex());
          m_eclShowerToMC5GMothPDG->push_back(amcParticle->getMother()->getMother()->getPDG());
        } else {
          m_eclShowerToMC5GMoth->push_back(-999);
          m_eclShowerToMC5GMothPDG->push_back(-999);
        }
      } else {
        m_eclShowerToMC5Moth->push_back(-999);
        m_eclShowerToMC5MothPDG->push_back(-999);
        m_eclShowerToMC5GMoth->push_back(-999);
        m_eclShowerToMC5GMothPDG->push_back(-999);
      }
    } else {
      m_eclShowerToMCWeight5->push_back(-1);
      m_eclShowerToMC5->push_back(-1);
      m_eclShowerToMC5PDG->push_back(-1);
      m_eclShowerToMC5Moth->push_back(-1);
      m_eclShowerToMC5MothPDG->push_back(-1);
      m_eclShowerToMC5GMoth->push_back(-1);
      m_eclShowerToMC5GMothPDG->push_back(-1);
    }
    m_eclShowerToBkgWeight->push_back(aECLShowers->getEnergy() - sumHit);
  }
 
  m_mcMultip = m_mcParticles.getEntries();
  for (int imcpart = 0; imcpart < m_mcParticles.getEntries(); imcpart++) {
    const MCParticle* amcParticle = m_mcParticles[imcpart];
    m_mcIdx->push_back(amcParticle->getArrayIndex());
    m_mcPdg->push_back(amcParticle->getPDG());
    if (amcParticle->getMother() != nullptr) m_mcMothPdg->push_back(amcParticle->getMother()->getPDG());
    else m_mcMothPdg->push_back(-999);
    if (amcParticle->getMother() != nullptr
        && amcParticle->getMother()->getMother() != nullptr) m_mcGMothPdg->push_back(amcParticle->getMother()->getMother()->getPDG());
    else m_mcGMothPdg->push_back(-999);
    if (amcParticle->getMother() != nullptr && amcParticle->getMother()->getMother() != nullptr
        && amcParticle->getMother()->getMother()->getMother() != nullptr)
      m_mcGGMothPdg->push_back(amcParticle->getMother()->getMother()->getMother()->getPDG());
    else m_mcGGMothPdg->push_back(-999);
    m_mcEnergy->push_back(amcParticle->getEnergy());
    m_mcSecondaryPhysProc->push_back(amcParticle->getSecondaryPhysicsProcess());
 
    m_mcPx->push_back(amcParticle->getMomentum().X());
    m_mcPy->push_back(amcParticle->getMomentum().Y());
    m_mcPz->push_back(amcParticle->getMomentum().Z());
 
    m_mcDecayVtxX->push_back(amcParticle->getDecayVertex().X());
    m_mcDecayVtxY->push_back(amcParticle->getDecayVertex().Y());
    m_mcDecayVtxZ->push_back(amcParticle->getDecayVertex().Z());
 
    m_mcProdVtxX->push_back(amcParticle->getProductionVertex().X());
    m_mcProdVtxY->push_back(amcParticle->getProductionVertex().Y());
    m_mcProdVtxZ->push_back(amcParticle->getProductionVertex().Z());
 
  }
 
  if (m_doTracking == true) {
    m_trkMultip = 0;
    for (const Track& itrk : m_tracks) {
      const TrackFitResult* atrk = itrk.getTrackFitResult(Const::pion);
      if (atrk == nullptr) continue;
 
      m_trkIdx->push_back(m_trkMultip);
      m_trkPdg->push_back(atrk->getParticleType().getPDGCode());
      m_trkCharge->push_back(atrk->getChargeSign());
      m_trkPx->push_back(atrk->getMomentum().X());
      m_trkPy->push_back(atrk->getMomentum().Y());
      m_trkPz->push_back(atrk->getMomentum().Z());
      m_trkP->push_back(atrk->getMomentum().R());
      m_trkTheta->push_back(atrk->getMomentum().Theta());
      m_trkPhi->push_back(atrk->getMomentum().Phi());
      m_trkX->push_back(atrk->getPosition().X());
      m_trkY->push_back(atrk->getPosition().Y());
      m_trkZ->push_back(atrk->getPosition().Z());
 
      const ECLPidLikelihood* eclpid = itrk.getRelatedTo<ECLPidLikelihood>() ;
 
      if (eclpid != nullptr) {
        m_eclpidtrkIdx -> push_back(m_trkMultip);
        m_eclpidEnergy -> push_back(eclpid-> energy());
        m_eclpidEop    -> push_back(eclpid-> eop());
        m_eclpidE9E21  -> push_back(eclpid-> e9e25());
        m_eclpidNCrystals -> push_back(eclpid-> nCrystals());
        m_eclpidNClusters -> push_back(eclpid-> nClusters());
        m_eclLogLikeEl -> push_back(eclpid-> getLogLikelihood(Const::electron));
        m_eclLogLikeMu -> push_back(eclpid-> getLogLikelihood(Const::muon));
        m_eclLogLikePi -> push_back(eclpid-> getLogLikelihood(Const::pion));
      } else {
        m_eclpidtrkIdx -> push_back(m_trkMultip);
        m_eclpidEnergy -> push_back(0);
        m_eclpidEop    -> push_back(0);
        m_eclpidE9E21  -> push_back(0);
        m_eclpidNCrystals -> push_back(0);
        m_eclpidNClusters -> push_back(0);
        m_eclLogLikeEl -> push_back(0);
        m_eclLogLikeMu -> push_back(0);
        m_eclLogLikePi -> push_back(0);
      }
      m_trkMultip++;
    }
  }
 
  m_tree->Fill();
}

◆ exposePythonAPI()

void exposePythonAPI ( )

staticinherited

Exposes methods of the Module class to Python.

Definition at line 325 of file Module.cc.

{
  // to avoid confusion between std::arg and boost::python::arg we want a shorthand namespace as well
  namespace bp = boost::python;
 
  docstring_options options(true, true, false); //userdef, py sigs, c++ sigs
 
  void (Module::*setReturnValueInt)(int) = &Module::setReturnValue;
 
  enum_<Module::EAfterConditionPath>("AfterConditionPath",
                                     R"(Determines execution behaviour after a conditional path has been executed:
 
.. attribute:: END
 
  End processing of this path after the conditional path. (this is the default for if_value() etc.)
 
.. attribute:: CONTINUE
 
  After the conditional path, resume execution after this module.)")
  .value("END", Module::EAfterConditionPath::c_End)
  .value("CONTINUE", Module::EAfterConditionPath::c_Continue)
  ;
 
  /* Do not change the names of >, <, ... we use them to serialize conditional paths */
  enum_<Belle2::ModuleCondition::EConditionOperators>("ConditionOperator")
  .value(">", Belle2::ModuleCondition::EConditionOperators::c_GT)
  .value("<", Belle2::ModuleCondition::EConditionOperators::c_ST)
  .value(">=", Belle2::ModuleCondition::EConditionOperators::c_GE)
  .value("<=", Belle2::ModuleCondition::EConditionOperators::c_SE)
  .value("==", Belle2::ModuleCondition::EConditionOperators::c_EQ)
  .value("!=", Belle2::ModuleCondition::EConditionOperators::c_NE)
  ;
 
  enum_<Module::EModulePropFlags>("ModulePropFlags",
                                  R"(Flags to indicate certain low-level features of modules, see :func:`Module.set_property_flags()`, :func:`Module.has_properties()`. Most useful flags are:
 
.. attribute:: PARALLELPROCESSINGCERTIFIED
 
    This module can be run in parallel processing mode safely (All I/O must be done through the data store, in particular, the module must not write any files.)
 
.. attribute:: HISTOGRAMMANAGER
 
  This module is used to manage histograms accumulated by other modules
 
.. attribute:: TERMINATEINALLPROCESSES
 
  When using parallel processing, call this module's terminate() function in all processes. This will also ensure that there is exactly one process (single-core if no parallel modules found) or at least one input, one main and one output process.
)")
  .value("INPUT", Module::EModulePropFlags::c_Input)
  .value("OUTPUT", Module::EModulePropFlags::c_Output)
  .value("PARALLELPROCESSINGCERTIFIED", Module::EModulePropFlags::c_ParallelProcessingCertified)
  .value("HISTOGRAMMANAGER", Module::EModulePropFlags::c_HistogramManager)
  .value("INTERNALSERIALIZER", Module::EModulePropFlags::c_InternalSerializer)
  .value("TERMINATEINALLPROCESSES", Module::EModulePropFlags::c_TerminateInAllProcesses)
  ;
 
  //Python class definition
  class_<Module, PyModule> module("Module", R"(
Base class for Modules.
 
A module is the smallest building block of the framework.
A typical event processing chain consists of a Path containing
modules. By inheriting from this base class, various types of
modules can be created. To use a module, please refer to
:func:`Path.add_module()`.  A list of modules is available by running
``basf2 -m`` or ``basf2 -m package``, detailed information on parameters is
given by e.g. ``basf2 -m RootInput``.
 
The 'Module Development' section in the manual provides detailed information
on how to create modules, setting parameters, or using return values/conditions:
https://xwiki.desy.de/xwiki/rest/p/f4fa4/#HModuleDevelopment
 
)");
  module
  .def("__str__", &Module::getPathString)
  .def("name", &Module::getName, return_value_policy<copy_const_reference>(),
       "Returns the name of the module. Can be changed via :func:`set_name() <Module.set_name()>`, use :func:`type() <Module.type()>` for identifying a particular module class.")
  .def("type", &Module::getType, return_value_policy<copy_const_reference>(),
       "Returns the type of the module (i.e. class name minus 'Module')")
  .def("set_name", &Module::setName, args("name"), R"(
Set custom name, e.g. to distinguish multiple modules of the same type.
 
>>> path.add_module('EventInfoSetter')
>>> ro = path.add_module('RootOutput', branchNames=['EventMetaData'])
>>> ro.set_name('RootOutput_metadata_only')
>>> print(path)
[EventInfoSetter -> RootOutput_metadata_only]
 
)")
  .def("description", &Module::getDescription, return_value_policy<copy_const_reference>(),
       "Returns the description of this module.")
  .def("package", &Module::getPackage, return_value_policy<copy_const_reference>(),
       "Returns the package this module belongs to.")
  .def("available_params", &_getParamInfoListPython,
       "Return list of all module parameters as `ModuleParamInfo` instances")
  .def("has_properties", &Module::hasProperties, (bp::arg("properties")),
       R"DOCSTRING(Allows to check if the module has the given properties out of `ModulePropFlags` set.
 
>>> if module.has_properties(ModulePropFlags.PARALLELPROCESSINGCERTIFIED):
>>>     ...
 
Parameters:
  properties (int): bitmask of `ModulePropFlags` to check for.
)DOCSTRING")
  .def("set_property_flags", &Module::setPropertyFlags, args("property_mask"),
       "Set module properties in the form of an OR combination of `ModulePropFlags`.");
  {
    // python signature is too crowded, make ourselves
    docstring_options subOptions(true, false, false); //userdef, py sigs, c++ sigs
    module
    .def("if_value", &Module::if_value,
         (bp::arg("expression"), bp::arg("condition_path"), bp::arg("after_condition_path")= Module::EAfterConditionPath::c_End),
         R"DOCSTRING(if_value(expression, condition_path, after_condition_path=AfterConditionPath.END)
 
Sets a conditional sub path which will be executed after this
module if the return value set in the module passes the given ``expression``.
 
Modules can define a return value (int or bool) using ``setReturnValue()``,
which can be used in the steering file to split the Path based on this value, for example
 
>>> module_with_condition.if_value("<1", another_path)
 
In case the return value of the ``module_with_condition`` for a given event is
less than 1, the execution will be diverted into ``another_path`` for this event.
 
You could for example set a special return value if an error occurs, and divert
the execution into a path containing :b2:mod:`RootOutput` if it is found;
saving only the data producing/produced by the error.
 
After a conditional path has executed, basf2 will by default stop processing
the path for this event. This behaviour can be changed by setting the
``after_condition_path`` argument.
 
Parameters:
  expression (str): Expression to determine if the conditional path should be executed.
      This should be one of the comparison operators ``<``, ``>``, ``<=``,
      ``>=``, ``==``, or ``!=`` followed by a numerical value for the return value
  condition_path (Path): path to execute in case the expression is fulfilled
  after_condition_path (AfterConditionPath): What to do once the ``condition_path`` has been executed.
)DOCSTRING")
    .def("if_false", &Module::if_false,
         (bp::arg("condition_path"), bp::arg("after_condition_path")= Module::EAfterConditionPath::c_End),
         R"DOC(if_false(condition_path, after_condition_path=AfterConditionPath.END)
 
Sets a conditional sub path which will be executed after this module if
the return value of the module evaluates to False. This is equivalent to
calling `if_value` with ``expression=\"<1\"``)DOC")
    .def("if_true", &Module::if_true,
         (bp::arg("condition_path"), bp::arg("after_condition_path")= Module::EAfterConditionPath::c_End),
         R"DOC(if_true(condition_path, after_condition_path=AfterConditionPath.END)
 
Sets a conditional sub path which will be executed after this module if
the return value of the module evaluates to True. It is equivalent to
calling `if_value` with ``expression=\">=1\"``)DOC");
  }
  module
  .def("has_condition", &Module::hasCondition,
       "Return true if a conditional path has been set for this module "
       "using `if_value`, `if_true` or `if_false`")
  .def("get_all_condition_paths", &_getAllConditionPathsPython,
       "Return a list of all conditional paths set for this module using "
       "`if_value`, `if_true` or `if_false`")
  .def("get_all_conditions", &_getAllConditionsPython,
       "Return a list of all conditional path expressions set for this module using "
       "`if_value`, `if_true` or `if_false`")
  .add_property("logging", make_function(&Module::getLogConfig, return_value_policy<reference_existing_object>()),
                &Module::setLogConfig)

◆ getAfterConditionPath()

Module::EAfterConditionPath getAfterConditionPath ( ) const

inherited

What to do after the conditional path is finished.

(defaults to c_End if no condition is set)

Definition at line 133 of file Module.cc.

{
  if (m_conditions.empty()) return EAfterConditionPath::c_End;
 
  //okay, a condition was set for this Module...
  if (!m_hasReturnValue) {
    B2FATAL("A condition was set for '" << getName() << "', but the module did not set a return value!");
  }
 
  for (const auto& condition : m_conditions) {
    if (condition.evaluate(m_returnValue)) {
      return condition.getAfterConditionPath();
    }
  }
 
  return EAfterConditionPath::c_End;
}

◆ getAllConditionPaths()

std::vector< std::shared_ptr< Path > > getAllConditionPaths ( ) const

inherited

Return all condition paths currently set (no matter if the condition is true or not).

Definition at line 150 of file Module.cc.

{
  std::vector<std::shared_ptr<Path>> allConditionPaths;
  for (const auto& condition : m_conditions) {
    allConditionPaths.push_back(condition.getPath());
  }
 
  return allConditionPaths;
}

◆ getAllConditions()

const std::vector< ModuleCondition > & getAllConditions ( ) const

inlineinherited

Return all set conditions for this module.

Definition at line 323 of file Module.h.

    {
      return m_conditions;
    }

◆ getCondition()

const ModuleCondition * getCondition ( ) const

inlineinherited

Return a pointer to the first condition (or nullptr, if none was set)

Definition at line 313 of file Module.h.

    {
      if (m_conditions.empty()) {
        return nullptr;
      } else {
        return &m_conditions.front();
      }
    }

◆ getConditionPath()

std::shared_ptr< Path > getConditionPath ( ) const

inherited

Returns the path of the last true condition (if there is at least one, else reaturn a null pointer).

Definition at line 113 of file Module.cc.

{
  PathPtr p;
  if (m_conditions.empty()) return p;
 
  //okay, a condition was set for this Module...
  if (!m_hasReturnValue) {
    B2FATAL("A condition was set for '" << getName() << "', but the module did not set a return value!");
  }
 
  for (const auto& condition : m_conditions) {
    if (condition.evaluate(m_returnValue)) {
      return condition.getPath();
    }
  }
 
  // if none of the conditions were true, return a null pointer.
  return p;
}

◆ getDescription()

const std::string & getDescription ( ) const

inlineinherited

Returns the description of the module.

Definition at line 201 of file Module.h.

201{return m_description;}

◆ getFileNames()

virtual std::vector< std::string > getFileNames ( bool outputFiles )

inlinevirtualinherited

Return a list of output filenames for this modules.

This will be called when basf2 is run with "--dry-run" if the module has set either the c_Input or c_Output properties.

If the parameter outputFiles is false (for modules with c_Input) the list of input filenames should be returned (if any). If outputFiles is true (for modules with c_Output) the list of output files should be returned (if any).

If a module has sat both properties this member is called twice, once for each property.

The module should return the actual list of requested input or produced output filenames (including handling of input/output overrides) so that the grid system can handle input/output files correctly.

This function should return the same value when called multiple times. This is especially important when taking the input/output overrides from Environment as they get consumed when obtained so the finalized list of output files should be stored for subsequent calls.

Reimplemented in RootInputModule, RootOutputModule, and StorageRootOutputModule.

Definition at line 133 of file Module.h.

    {
      return std::vector<std::string>();
    }

◆ getLogConfig()

LogConfig & getLogConfig ( )

inlineinherited

Returns the log system configuration.

Definition at line 224 of file Module.h.

224{return m_logConfig;}

◆ getModules()

std::list< ModulePtr > getModules ( ) const

inlineoverrideprivatevirtualinherited

no submodules, return empty list

Implements PathElement.

Definition at line 505 of file Module.h.

505{ return std::list<ModulePtr>(); }

◆ getName()

const std::string & getName ( ) const

inlineinherited

Returns the name of the module.

This can be changed via e.g. set_name() in the steering file to give more useful names if there is more than one module of the same type.

For identifying the type of a module, using getType() (or type() in Python) is recommended.

Definition at line 186 of file Module.h.

186{return m_name;}

◆ getPackage()

const std::string & getPackage ( ) const

inlineinherited

Returns the package this module is in.

Definition at line 196 of file Module.h.

196{return m_package;}

◆ getParamInfoListPython()

std::shared_ptr< boost::python::list > getParamInfoListPython ( ) const

inherited

Returns a python list of all parameters.

Each item in the list consists of the name of the parameter, a string describing its type, a python list of all default values and the description of the parameter.

Returns: A python list containing the parameters of this parameter list.

Definition at line 279 of file Module.cc.

{
  return m_moduleParamList.getParamInfoListPython();
}

◆ getParamList()

const ModuleParamList & getParamList ( ) const

inlineinherited

Return module param list.

Definition at line 362 of file Module.h.

362{ return m_moduleParamList; }

◆ getPathString()

std::string getPathString ( ) const

overrideprivatevirtualinherited

return the module name.

Implements PathElement.

Definition at line 192 of file Module.cc.

{
 
  std::string output = getName();
 
  for (const auto& condition : m_conditions) {
    output += condition.getString();
  }
 
  return output;
}

◆ getReturnValue()

int getReturnValue ( ) const

inlineinherited

Return the return value set by this module.

This value is only meaningful if hasReturnValue() is true

Definition at line 380 of file Module.h.

380{ return m_returnValue; }

◆ getType()

const std::string & getType ( ) const

inherited

Returns the type of the module (i.e.

class name minus 'Module')

Definition at line 41 of file Module.cc.

{
  if (m_type.empty())
    B2FATAL("Module type not set for " << getName());
  return m_type;
}

◆ hasCondition()

bool hasCondition ( ) const

inlineinherited

Returns true if at least one condition was set for the module.

Definition at line 310 of file Module.h.

310{ return not m_conditions.empty(); };

◆ hasProperties()

bool hasProperties ( unsigned int propertyFlags ) const

inherited

Returns true if all specified property flags are available in this module.

Parameters

propertyFlags Ored EModulePropFlags which should be compared with the module flags.

Definition at line 160 of file Module.cc.

{
  return (propertyFlags & m_propertyFlags) == propertyFlags;
}

◆ hasReturnValue()

bool hasReturnValue ( ) const

inlineinherited

Return true if this module has a valid return value set.

Definition at line 377 of file Module.h.

377{ return m_hasReturnValue; }

◆ hasUnsetForcedParams()

bool hasUnsetForcedParams ( ) const

inherited

Returns true and prints error message if the module has unset parameters which the user has to set in the steering file.

Definition at line 166 of file Module.cc.

{
  auto missing = m_moduleParamList.getUnsetForcedParams();
  std::string allMissing = "";
  for (const auto& s : missing)
    allMissing += s + " ";
  if (!missing.empty())
    B2ERROR("The following required parameters of Module '" << getName() << "' were not specified: " << allMissing <<
            "\nPlease add them to your steering file.");
  return !missing.empty();
}

◆ if_false()

void if_false	(	const std::shared_ptr< Path > &	path,
		EAfterConditionPath	afterConditionPath = EAfterConditionPath::c_End )

inherited

A simplified version to add a condition to the module.

Please note that successive calls of this function will add more than one condition to the module. If more than one condition results in true, only the last of them will be used.

Please be careful: Avoid creating cyclic paths, e.g. by linking a condition to a path which is processed before the path where this module is located in.

It is equivalent to the if_value() method, using the expression "<1". This method is meant to be used together with the setReturnValue(bool value) method.

Parameters

path	Shared pointer to the Path which will be executed if the return value is false.
afterConditionPath	What to do after executing 'path'.

Definition at line 85 of file Module.cc.

{
  if_value("<1", path, afterConditionPath);
}

◆ if_true()

void if_true	(	const std::shared_ptr< Path > &	path,
		EAfterConditionPath	afterConditionPath = EAfterConditionPath::c_End )

inherited

A simplified version to set the condition of the module.

Please note that successive calls of this function will add more than one condition to the module. If more than one condition results in true, only the last of them will be used.

Please be careful: Avoid creating cyclic paths, e.g. by linking a condition to a path which is processed before the path where this module is located in.

It is equivalent to the if_value() method, using the expression ">=1". This method is meant to be used together with the setReturnValue(bool value) method.

Parameters

path	Shared pointer to the Path which will be executed if the return value is true.
afterConditionPath	What to do after executing 'path'.

Definition at line 90 of file Module.cc.

{
  if_value(">=1", path, afterConditionPath);
}

◆ if_value()

void if_value	(	const std::string &	expression,
		const std::shared_ptr< Path > &	path,
		EAfterConditionPath	afterConditionPath = EAfterConditionPath::c_End )

inherited

Add a condition to the module.

Please note that successive calls of this function will add more than one condition to the module. If more than one condition results in true, only the last of them will be used.

See https://xwiki.desy.de/xwiki/rest/p/a94f2 or ModuleCondition for a description of the syntax.

Please be careful: Avoid creating cyclic paths, e.g. by linking a condition to a path which is processed before the path where this module is located in.

Parameters

expression	The expression of the condition.
path	Shared pointer to the Path which will be executed if the condition is evaluated to true.
afterConditionPath	What to do after executing 'path'.

Definition at line 79 of file Module.cc.

{
  m_conditions.emplace_back(expression, path, afterConditionPath);
}

◆ initialize()

void initialize ( void )

overridevirtual

Initializes the Module.

SHOWERS

Reimplemented from Module.

Definition at line 440 of file ECLDataAnalysisModule.cc.

{
 
  B2INFO("[ECLDataAnalysis Module]: Starting initialization of ECLDataAnalysis Module.");
 
  m_eventLevelClusteringInfo.registerInDataStore();
 
  m_eclSimHits.registerInDataStore(eclSimHitArrayName());
  m_eclHits.registerInDataStore(eclHitArrayName());
 
  m_eclDigits.registerInDataStore(eclDigitArrayName());
  m_eclCalDigits.registerInDataStore(eclCalDigitArrayName());
  m_eclConnectedRegions.registerInDataStore(eclConnectedRegionArrayName());
  m_eclShowers.registerInDataStore(eclShowerArrayName());
  m_eclClusters.registerInDataStore(eclClusterArrayName());
  m_eclLocalMaximums.registerInDataStore(eclLocalMaximumArrayName());
 
  m_eclCalDigits.registerRelationTo(m_mcParticles);
  m_eclDigits.registerRelationTo(m_mcParticles);
  m_eclShowers.registerRelationTo(m_mcParticles);
  m_eclClusters.registerRelationTo(m_mcParticles);
 
  if (m_doPureCsI == 1) {
    m_eclPureDigits.registerInDataStore(eclPureDigitArrayName());
    m_eclPureCalDigits.registerInDataStore(eclPureCalDigitArrayName());
    m_eclPureConnectedRegions.registerInDataStore(eclPureConnectedRegionArrayName());
    m_eclPureShowers.registerInDataStore(eclPureShowerArrayName());
    m_eclPureClusters.registerInDataStore(eclPureClusterArrayName());
    m_eclPureLocalMaximums.registerInDataStore(eclPureLocalMaximumArrayName());
 
    m_eclPureCalDigits.registerRelationTo(m_mcParticles);
    m_eclPureDigits.registerRelationTo(m_mcParticles);
    m_eclPureShowers.registerRelationTo(m_mcParticles);
    m_eclPureClusters.registerRelationTo(m_mcParticles);
  }
 
  if (m_doTracking == true) {
    m_tracks.isRequired();
    m_trackFitResults.isRequired();
    m_eclPidLikelihoods.isRequired();
  }
 
  if (m_writeToRoot == true) {
    m_rootFilePtr = new TFile(m_rootFileName.c_str(), "RECREATE");
  } else
    m_rootFilePtr = nullptr;
 
  // initialize tree
  m_tree     = new TTree("m_tree", "ECL Analysis tree");
 
  m_tree->Branch("expNo", &m_iExperiment, "expNo/I");
  m_tree->Branch("runNo", &m_iRun, "runNo/I");
  m_tree->Branch("evtNo", &m_iEvent, "evtNo/I");
 
  //EventLevelClusteringInfo
  m_tree->Branch("eclNumOutOfTimeDigitsFwd",     &m_nECLCalDigitsOutOfTimeFWD,         "eclNumOutOfTimeDigitsFwd/s");
  m_tree->Branch("eclNumOutOfTimeDigitsBrl",     &m_nECLCalDigitsOutOfTimeBarrel,         "eclNumOutOfTimeDigitsBrl/s");
  m_tree->Branch("eclNumOutOfTimeDigitsBwd",     &m_nECLCalDigitsOutOfTimeBWD,         "eclNumOutOfTimeDigitsBwd/s");
  m_tree->Branch("eclNumRejectedShowersFwd",     &m_nECLShowersRejectedFWD,         "eclNumRejectedShowersFwd/b");
  m_tree->Branch("eclNumRejectedShowersBrl",     &m_nECLShowersRejectedBarrel,         "eclNumRejectedShowersBrl/b");
  m_tree->Branch("eclNumRejectedShowersBwd",     &m_nECLShowersRejectedBWD,         "eclNumRejectedShowersBwd/b");
 
  if (m_doDigits == 1) {
    m_tree->Branch("eclDigitMultip",     &m_eclDigitMultip,         "ecldigit_Multip/I");
    m_tree->Branch("eclDigitIdx",        "std::vector<int>",         &m_eclDigitIdx);
    m_tree->Branch("eclDigitToMC",      "std::vector<int>",          &m_eclDigitToMC);
    m_tree->Branch("eclDigitCellId",     "std::vector<int>",         &m_eclDigitCellId);
    m_tree->Branch("eclDigitAmp",        "std::vector<int>",         &m_eclDigitAmp);
    m_tree->Branch("eclDigitTimeFit",    "std::vector<int>",         &m_eclDigitTimeFit);
    m_tree->Branch("eclDigitFitQuality",    "std::vector<int>",         &m_eclDigitFitQuality);
    m_tree->Branch("eclDigitToCalDigit",      "std::vector<int>",          &m_eclDigitToCalDigit);
 
    m_tree->Branch("eclCalDigitMultip",     &m_eclCalDigitMultip,         "eclCaldigit_Multip/I");
    m_tree->Branch("eclCalDigitIdx",        "std::vector<int>",         &m_eclCalDigitIdx);
    m_tree->Branch("eclCalDigitToMC1",      "std::vector<int>",       &m_eclCalDigitToMC1);
    m_tree->Branch("eclCalDigitToMC1PDG",      "std::vector<int>",       &m_eclCalDigitToMC1PDG);
    m_tree->Branch("eclCalDigitToMCWeight1",      "std::vector<double>",       &m_eclCalDigitToMCWeight1);
    m_tree->Branch("eclCalDigitToMC2",      "std::vector<int>",       &m_eclCalDigitToMC2);
    m_tree->Branch("eclCalDigitToMC2PDG",      "std::vector<int>",       &m_eclCalDigitToMC2PDG);
    m_tree->Branch("eclCalDigitToMCWeight2",      "std::vector<double>",       &m_eclCalDigitToMCWeight2);
    m_tree->Branch("eclCalDigitToMC3",      "std::vector<int>",       &m_eclCalDigitToMC3);
    m_tree->Branch("eclCalDigitToMC3PDG",      "std::vector<int>",       &m_eclCalDigitToMC3PDG);
    m_tree->Branch("eclCalDigitToMCWeight3",      "std::vector<double>",       &m_eclCalDigitToMCWeight3);
    m_tree->Branch("eclCalDigitToMC4",      "std::vector<int>",       &m_eclCalDigitToMC4);
    m_tree->Branch("eclCalDigitToMC4PDG",      "std::vector<int>",       &m_eclCalDigitToMC4PDG);
    m_tree->Branch("eclCalDigitToMCWeight4",      "std::vector<double>",       &m_eclCalDigitToMCWeight4);
    m_tree->Branch("eclCalDigitToMC5",      "std::vector<int>",       &m_eclCalDigitToMC5);
    m_tree->Branch("eclCalDigitToMC5PDG",      "std::vector<int>",       &m_eclCalDigitToMC5PDG);
    m_tree->Branch("eclCalDigitToMCWeight5",      "std::vector<double>",       &m_eclCalDigitToMCWeight5);
    m_tree->Branch("eclCalDigitToBkgWeight",      "std::vector<double>",       &m_eclCalDigitToBkgWeight);
    m_tree->Branch("eclCalDigitSimHitSum",      "std::vector<double>",       &m_eclCalDigitSimHitSum);
    m_tree->Branch("eclCalDigitToShower",      "std::vector<int>",          &m_eclCalDigitToShower);
    m_tree->Branch("eclCalDigitCellId",     "std::vector<int>",         &m_eclCalDigitCellId);
    m_tree->Branch("eclCalDigitAmp",        "std::vector<double>",         &m_eclCalDigitAmp);
    m_tree->Branch("eclCalDigitTimeFit",    "std::vector<double>",         &m_eclCalDigitTimeFit);
    m_tree->Branch("eclCalDigitFitQuality",    "std::vector<int>",         &m_eclCalDigitFitQuality);
    m_tree->Branch("eclCalDigitToCR",      "std::vector<int>",          &m_eclCalDigitToCR);
    m_tree->Branch("eclCalDigitToLM",      "std::vector<int>",          &m_eclCalDigitToLM);
  }
 
  m_tree->Branch("eclCRIdx", "std::vector<int>", &m_eclCRIdx);
  m_tree->Branch("eclCRIsTrack", "std::vector<int>", &m_eclCRIsTrack);
  m_tree->Branch("eclCRLikelihoodMIPNGamma", "std::vector<double>", &m_eclCRLikelihoodMIPNGamma);
  m_tree->Branch("eclCRLikelihoodChargedHadron", "std::vector<double>", &m_eclCRLikelihoodChargedHadron);
  m_tree->Branch("eclCRLikelihoodElectronNGamma", "std::vector<double>", &m_eclCRLikelihoodElectronNGamma);
  m_tree->Branch("eclCRLikelihoodNGamma", "std::vector<double>", &m_eclCRLikelihoodNGamma);
  m_tree->Branch("eclCRLikelihoodNeutralHadron", "std::vector<double>", &m_eclCRLikelihoodNeutralHadron);
  m_tree->Branch("eclCRLikelihoodMergedPi0", "std::vector<double>", &m_eclCRLikelihoodMergedPi0);
 
  if (m_doHits == 1) {
    m_tree->Branch("eclSimHitMultip",     &m_eclSimHitMultip,      "eclSimHitMultip/I");
    m_tree->Branch("eclSimHitIdx",     "std::vector<int>",       &m_eclSimHitIdx);
    m_tree->Branch("eclSimHitToMC",      "std::vector<int>",       &m_eclSimHitToMC);
    m_tree->Branch("eclSimHitCellId",     "std::vector<int>",       &m_eclSimHitCellId);
    m_tree->Branch("eclSimHitPdg",        "std::vector<int>",       &m_eclSimHitPdg);
    m_tree->Branch("eclSimHitEnergyDep",      "std::vector<double>",    &m_eclSimHitEnergyDep);
    m_tree->Branch("eclSimHitFlightTime",      "std::vector<double>",    &m_eclSimHitFlightTime);
    m_tree->Branch("eclSimHitX",          "std::vector<double>",    &m_eclSimHitX);
    m_tree->Branch("eclSimHitY",          "std::vector<double>",    &m_eclSimHitY);
    m_tree->Branch("eclSimHitZ",          "std::vector<double>",    &m_eclSimHitZ);
    m_tree->Branch("eclSimHitPx",         "std::vector<double>",    &m_eclSimHitPx);
    m_tree->Branch("eclSimHitPy",         "std::vector<double>",    &m_eclSimHitPy);
    m_tree->Branch("eclSimHitPz",         "std::vector<double>",    &m_eclSimHitPz);
 
    m_tree->Branch("eclHitMultip",     &m_eclHitMultip,      "eclHitMultip/I");
    m_tree->Branch("eclHitIdx",     "std::vector<int>",       &m_eclHitIdx);
    m_tree->Branch("eclHitToMC",      "std::vector<int>",       &m_eclHitToMC);
    m_tree->Branch("eclHitToDigit",      "std::vector<int>",       &m_eclHitToDigit);
    m_tree->Branch("eclHitToDigitAmp",      "std::vector<int>",       &m_eclHitToDigitAmp);
    m_tree->Branch("eclHitCellId",     "std::vector<int>",    &m_eclHitCellId);
    m_tree->Branch("eclHitEnergyDep",      "std::vector<double>", &m_eclHitEnergyDep);
    m_tree->Branch("eclHitTimeAve",       "std::vector<double>", &m_eclHitTimeAve);
  }
 
  m_tree->Branch("eclLMMultip",     &m_eclLMMultip,     "eclLMMultip/I");
  m_tree->Branch("eclLMId",    "std::vector<int>",     &m_eclLMId);
  m_tree->Branch("eclLMType",    "std::vector<int>",     &m_eclLMType);
  m_tree->Branch("eclLMCellId",     "std::vector<int>",     &m_eclLMCellId);
 
  m_tree->Branch("eclClusterMultip",     &m_eclClusterMultip,     "eclClusterMultip/I");
  m_tree->Branch("eclClusterTrueMultip",     &m_eclClusterTrueMultip,     "eclClusterTrueMultip/I");
  m_tree->Branch("eclClusterGammaMultip",     &m_eclClusterGammaMultip,     "eclClusterGammaMultip/I");
  m_tree->Branch("eclClusterIdx",     "std::vector<int>",       &m_eclClusterIdx);
  m_tree->Branch("eclClusterToMC1",      "std::vector<int>",       &m_eclClusterToMC1);
  m_tree->Branch("eclClusterToMCWeight1",      "std::vector<double>",       &m_eclClusterToMCWeight1);
  m_tree->Branch("eclClusterToMC1PDG",      "std::vector<int>",       &m_eclClusterToMC1PDG);
  m_tree->Branch("eclClusterToMC2",      "std::vector<int>",       &m_eclClusterToMC2);
  m_tree->Branch("eclClusterToMCWeight2",      "std::vector<double>",       &m_eclClusterToMCWeight2);
  m_tree->Branch("eclClusterToMC2PDG",      "std::vector<int>",       &m_eclClusterToMC2PDG);
  m_tree->Branch("eclClusterToMC3",      "std::vector<int>",       &m_eclClusterToMC3);
  m_tree->Branch("eclClusterToMCWeight3",      "std::vector<double>",       &m_eclClusterToMCWeight3);
  m_tree->Branch("eclClusterToMC3PDG",      "std::vector<int>",       &m_eclClusterToMC3PDG);
  m_tree->Branch("eclClusterToMC4",      "std::vector<int>",       &m_eclClusterToMC4);
  m_tree->Branch("eclClusterToMCWeight4",      "std::vector<double>",       &m_eclClusterToMCWeight4);
  m_tree->Branch("eclClusterToMC4PDG",      "std::vector<int>",       &m_eclClusterToMC4PDG);
  m_tree->Branch("eclClusterToMC5",      "std::vector<int>",       &m_eclClusterToMC5);
  m_tree->Branch("eclClusterToMCWeight5",      "std::vector<double>",       &m_eclClusterToMCWeight5);
  m_tree->Branch("eclClusterToMC5PDG",      "std::vector<int>",       &m_eclClusterToMC5PDG);
  m_tree->Branch("eclClusterToBkgWeight",      "std::vector<double>",       &m_eclClusterToBkgWeight);
  m_tree->Branch("eclClusterSimHitSum",      "std::vector<double>",       &m_eclClusterSimHitSum);
  m_tree->Branch("eclClusterToShower",      "std::vector<int>",       &m_eclClusterToShower);
  m_tree->Branch("eclClusterEnergy",     "std::vector<double>",    &m_eclClusterEnergy);
  m_tree->Branch("eclClusterEnergyError",  "std::vector<double>",    &m_eclClusterEnergyError);
  m_tree->Branch("eclClusterTheta",      "std::vector<double>",    &m_eclClusterTheta);
  m_tree->Branch("eclClusterThetaError",   "std::vector<double>",    &m_eclClusterThetaError);
  m_tree->Branch("eclClusterPhi",        "std::vector<double>",    &m_eclClusterPhi);
  m_tree->Branch("eclClusterPhiError",     "std::vector<double>",    &m_eclClusterPhiError);
  m_tree->Branch("eclClusterR",          "std::vector<double>",    &m_eclClusterR);
  m_tree->Branch("eclClusterEnergyDepSum",   "std::vector<double>",    &m_eclClusterEnergyDepSum);
  m_tree->Branch("eclClusterTiming",     "std::vector<double>",    &m_eclClusterTiming);
  m_tree->Branch("eclClusterTimingError",  "std::vector<double>",    &m_eclClusterTimingError);
  m_tree->Branch("eclClusterE9oE21",     "std::vector<double>",    &m_eclClusterE9oE21);
  m_tree->Branch("eclClusterHighestE",   "std::vector<double>",    &m_eclClusterHighestE);
  m_tree->Branch("eclClusterCellId",   "std::vector<int>",    &m_eclClusterCellId);
  m_tree->Branch("eclClusterNofCrystals",   "std::vector<int>",       &m_eclClusterNofCrystals);
  m_tree->Branch("eclClusterCrystalHealth", "std::vector<int>",       &m_eclClusterCrystalHealth);
  m_tree->Branch("eclClusterIsTrack",    "std::vector<bool>",       &m_eclClusterIsTrack);
  m_tree->Branch("eclClusterClosestTrackDist",    "std::vector<double>",       &m_eclClusterClosestTrackDist);
  m_tree->Branch("eclClusterDeltaL",     "std::vector<double>",    &m_eclClusterDeltaL);
  m_tree->Branch("eclClusterAbsZernike40",     "std::vector<double>",    &m_eclClusterAbsZernike40);
  m_tree->Branch("eclClusterAbsZernike51",     "std::vector<double>",    &m_eclClusterAbsZernike51);
  m_tree->Branch("eclClusterZernikeMVA",     "std::vector<double>",    &m_eclClusterZernikeMVA);
  m_tree->Branch("eclClusterE1oE9",     "std::vector<double>",    &m_eclClusterE1oE9);
  m_tree->Branch("eclClusterSecondMoment",     "std::vector<double>",    &m_eclClusterSecondMoment);
  m_tree->Branch("eclClusterLAT",     "std::vector<double>",    &m_eclClusterLAT);
  m_tree->Branch("eclClusterDeltaTime99",     "std::vector<double>",    &m_eclClusterDeltaTime99);
  m_tree->Branch("eclClusterDetectorRegion",     "std::vector<int>",    &m_eclClusterDetectorRegion);
  m_tree->Branch("eclClusterHasNPhotonHypothesis",     "std::vector<int>",    &m_eclClusterHasNPhotonHypothesis);
  m_tree->Branch("eclClusterHasNeutralHadronHypothesis",     "std::vector<int>",    &m_eclClusterHasNeutralHadronHypothesis);
 
  if (m_doPureCsI == true) {
    m_tree->Branch("eclHitToPureDigit",      "std::vector<int>",       &m_eclHitToPureDigit);
    m_tree->Branch("eclHitToPureDigitAmp",      "std::vector<int>",       &m_eclHitToPureDigitAmp);
    //Pure Digit
    m_tree->Branch("eclPureDigitMultip",     &m_eclPureDigitMultip,         "eclPureDigit_Multip/I");
    m_tree->Branch("eclPureDigitIdx",        "std::vector<int>",         &m_eclPureDigitIdx);
    m_tree->Branch("eclPureDigitToMC",      "std::vector<int>",          &m_eclPureDigitToMC);
    m_tree->Branch("eclPureDigitCellId",     "std::vector<int>",         &m_eclPureDigitCellId);
    m_tree->Branch("eclPureDigitAmp",        "std::vector<int>",         &m_eclPureDigitAmp);
    m_tree->Branch("eclPureDigitTimeFit",    "std::vector<int>",         &m_eclPureDigitTimeFit);
    m_tree->Branch("eclPureDigitFitQuality",    "std::vector<int>",         &m_eclPureDigitFitQuality);
    m_tree->Branch("eclPureDigitToCluster",        "std::vector<int>",         &m_eclPureDigitToCluster);
    //Pure CalDigit
    m_tree->Branch("eclPureCalDigitMultip",     &m_eclPureCalDigitMultip, "eclPureCalDigit_Multip/I");
    m_tree->Branch("eclPureCalDigitIdx",      "std::vector<int>",       &m_eclPureCalDigitIdx);
    m_tree->Branch("eclPureCalDigitToMC1",      "std::vector<int>",       &m_eclPureCalDigitToMC1);
    m_tree->Branch("eclPureCalDigitToMC1PDG",      "std::vector<int>",       &m_eclPureCalDigitToMC1PDG);
    m_tree->Branch("eclPureCalDigitToMCWeight1",      "std::vector<double>",       &m_eclPureCalDigitToMCWeight1);
    m_tree->Branch("eclPureCalDigitToMC2",      "std::vector<int>",       &m_eclPureCalDigitToMC2);
    m_tree->Branch("eclPureCalDigitToMC2PDG",      "std::vector<int>",       &m_eclPureCalDigitToMC2PDG);
    m_tree->Branch("eclPureCalDigitToMCWeight2",      "std::vector<double>",       &m_eclPureCalDigitToMCWeight2);
    m_tree->Branch("eclPureCalDigitToMC3",      "std::vector<int>",       &m_eclPureCalDigitToMC3);
    m_tree->Branch("eclPureCalDigitToMC3PDG",      "std::vector<int>",       &m_eclPureCalDigitToMC3PDG);
    m_tree->Branch("eclPureCalDigitToMCWeight3",      "std::vector<double>",       &m_eclPureCalDigitToMCWeight3);
    m_tree->Branch("eclPureCalDigitToMC4",      "std::vector<int>",       &m_eclPureCalDigitToMC4);
    m_tree->Branch("eclPureCalDigitToMC4PDG",      "std::vector<int>",       &m_eclPureCalDigitToMC4PDG);
    m_tree->Branch("eclPureCalDigitToMCWeight4",      "std::vector<double>",       &m_eclPureCalDigitToMCWeight4);
    m_tree->Branch("eclPureCalDigitToMC5",      "std::vector<int>",       &m_eclPureCalDigitToMC5);
    m_tree->Branch("eclPureCalDigitToMC5PDG",      "std::vector<int>",       &m_eclPureCalDigitToMC5PDG);
    m_tree->Branch("eclPureCalDigitToMCWeight5",      "std::vector<double>",       &m_eclPureCalDigitToMCWeight5);
    m_tree->Branch("eclPureCalDigitToBkgWeight",      "std::vector<double>",       &m_eclPureCalDigitToBkgWeight);
    m_tree->Branch("eclPureCalDigitSimHitSum",      "std::vector<double>",       &m_eclPureCalDigitSimHitSum);
    m_tree->Branch("eclPureCalDigitToShower",      "std::vector<int>",       &m_eclPureCalDigitToShower);
    m_tree->Branch("eclPureCalDigitCellId",      "std::vector<int>",       &m_eclPureCalDigitCellId);
    m_tree->Branch("eclPureCalDigitAmp",      "std::vector<double>",       &m_eclPureCalDigitAmp);
    m_tree->Branch("eclPureCalDigitTimeFit",      "std::vector<double>",       &m_eclPureCalDigitTimeFit);
    m_tree->Branch("eclPureCalDigitFitQuality",      "std::vector<int>",       &m_eclPureCalDigitFitQuality);
    m_tree->Branch("eclPureCalDigitToCR",      "std::vector<int>",       &m_eclPureCalDigitToCR);
    m_tree->Branch("eclPureCalDigitToLM",      "std::vector<int>",       &m_eclPureCalDigitToLM);
    m_tree->Branch("eclPureLMMultip",     &m_eclPureLMMultip,     "eclPureLMMultip/I");
 
    m_tree->Branch("eclPureCRIdx", "std::vector<int>", &m_eclPureCRIdx);
    m_tree->Branch("eclPureCRIsTrack", "std::vector<int>", &m_eclPureCRIsTrack);
    m_tree->Branch("eclPureCRLikelihoodMIPNGamma", "std::vector<double>", &m_eclPureCRLikelihoodMIPNGamma);
    m_tree->Branch("eclPureCRLikelihoodChargedHadron", "std::vector<double>", &m_eclPureCRLikelihoodChargedHadron);
    m_tree->Branch("eclPureCRLikelihoodElectronNGamma", "std::vector<double>", &m_eclPureCRLikelihoodElectronNGamma);
    m_tree->Branch("eclPureCRLikelihoodNGamma", "std::vector<double>", &m_eclPureCRLikelihoodNGamma);
    m_tree->Branch("eclPureCRLikelihoodNeutralHadron", "std::vector<double>", &m_eclPureCRLikelihoodNeutralHadron);
    m_tree->Branch("eclPureCRLikelihoodMergedPi0", "std::vector<double>", &m_eclPureCRLikelihoodMergedPi0);
 
    m_tree->Branch("eclPureLMId",    "std::vector<int>",     &m_eclPureLMId);
    m_tree->Branch("eclPureLMType",    "std::vector<int>",     &m_eclPureLMType);
    m_tree->Branch("eclPureLMCellId",     "std::vector<int>",     &m_eclPureLMCellId);
    //Pure Clusters
    m_tree->Branch("eclPureClusterMultip",     &m_eclPureClusterMultip,     "eclPureClusterMultip/I");
    m_tree->Branch("eclPureClusterIdx",     "std::vector<int>",       &m_eclPureClusterIdx);
    m_tree->Branch("eclPureClusterToMC1",      "std::vector<int>",       &m_eclPureClusterToMC1);
    m_tree->Branch("eclPureClusterToMCWeight1",      "std::vector<double>",       &m_eclPureClusterToMCWeight1);
    m_tree->Branch("eclPureClusterToMC1PDG",      "std::vector<int>",       &m_eclPureClusterToMC1PDG);
    m_tree->Branch("eclPureClusterToMC2",      "std::vector<int>",       &m_eclPureClusterToMC2);
    m_tree->Branch("eclPureClusterToMCWeight2",      "std::vector<double>",       &m_eclPureClusterToMCWeight2);
    m_tree->Branch("eclPureClusterToMC2PDG",      "std::vector<int>",       &m_eclPureClusterToMC2PDG);
    m_tree->Branch("eclPureClusterToMC3",      "std::vector<int>",       &m_eclPureClusterToMC3);
    m_tree->Branch("eclPureClusterToMCWeight3",      "std::vector<double>",       &m_eclPureClusterToMCWeight3);
    m_tree->Branch("eclPureClusterToMC3PDG",      "std::vector<int>",       &m_eclPureClusterToMC3PDG);
    m_tree->Branch("eclPureClusterToMC4",      "std::vector<int>",       &m_eclPureClusterToMC4);
    m_tree->Branch("eclPureClusterToMCWeight4",      "std::vector<double>",       &m_eclPureClusterToMCWeight4);
    m_tree->Branch("eclPureClusterToMC4PDG",      "std::vector<int>",       &m_eclPureClusterToMC4PDG);
    m_tree->Branch("eclPureClusterToMC5",      "std::vector<int>",       &m_eclPureClusterToMC5);
    m_tree->Branch("eclPureClusterToMCWeight5",      "std::vector<double>",       &m_eclPureClusterToMCWeight5);
    m_tree->Branch("eclPureClusterToMC5PDG",      "std::vector<int>",       &m_eclPureClusterToMC5PDG);
    m_tree->Branch("eclPureClusterToBkgWeight",      "std::vector<double>",       &m_eclPureClusterToBkgWeight);
    m_tree->Branch("eclPureClusterEnergy",     "std::vector<double>",    &m_eclPureClusterEnergy);
    m_tree->Branch("eclPureClusterEnergyError",  "std::vector<double>",    &m_eclPureClusterEnergyError);
    m_tree->Branch("eclPureClusterTheta",      "std::vector<double>",    &m_eclPureClusterTheta);
    m_tree->Branch("eclPureClusterThetaError",   "std::vector<double>",    &m_eclPureClusterThetaError);
    m_tree->Branch("eclPureClusterPhi",        "std::vector<double>",    &m_eclPureClusterPhi);
    m_tree->Branch("eclPureClusterPhiError",     "std::vector<double>",    &m_eclPureClusterPhiError);
    m_tree->Branch("eclPureClusterR",          "std::vector<double>",    &m_eclPureClusterR);
    m_tree->Branch("eclPureClusterEnergyDepSum",   "std::vector<double>",    &m_eclPureClusterEnergyDepSum);
    m_tree->Branch("eclPureClusterTiming",     "std::vector<double>",    &m_eclPureClusterTiming);
    m_tree->Branch("eclPureClusterTimingError",  "std::vector<double>",    &m_eclPureClusterTimingError);
    m_tree->Branch("eclPureClusterE9oE21",     "std::vector<double>",    &m_eclPureClusterE9oE21);
    m_tree->Branch("eclPureClusterHighestE",   "std::vector<double>",    &m_eclPureClusterHighestE);
    m_tree->Branch("eclPureClusterCellId",   "std::vector<int>",    &m_eclPureClusterCellId);
    m_tree->Branch("eclPureClusterLat",        "std::vector<double>",    &m_eclPureClusterLat);
    m_tree->Branch("eclPureClusterNofCrystals",   "std::vector<int>",       &m_eclPureClusterNofCrystals);
    m_tree->Branch("eclPureClusterCrystalHealth", "std::vector<int>",       &m_eclPureClusterCrystalHealth);
    m_tree->Branch("eclPureClusterClosestTrackDist",    "std::vector<double>",       &m_eclPureClusterClosestTrackDist);
    m_tree->Branch("eclPureClusterAbsZernike40",    "std::vector<double>",       &m_eclPureClusterAbsZernike40);
    m_tree->Branch("eclPureClusterAbsZernike51",    "std::vector<double>",       &m_eclPureClusterAbsZernike51);
    m_tree->Branch("eclPureClusterZernikeMVA",    "std::vector<double>",       &m_eclPureClusterZernikeMVA);
    m_tree->Branch("eclPureClusterSecondMoment",    "std::vector<double>",       &m_eclPureClusterSecondMoment);
    m_tree->Branch("eclPureClusterIsTrack",    "std::vector<bool>",       &m_eclPureClusterIsTrack);
    m_tree->Branch("eclPureClusterDeltaL",     "std::vector<double>",    &m_eclPureClusterDeltaL);
    m_tree->Branch("eclPureClusterE1oE9",         "std::vector<double>",    &m_eclPureClusterE1oE9);
    m_tree->Branch("eclPureClusterDeltaTime99",         "std::vector<double>",    &m_eclPureClusterDeltaTime99);
    m_tree->Branch("eclPureClusterDetectorRegion",    "std::vector<int>",       &m_eclPureClusterDetectorRegion);
    m_tree->Branch("eclPureClusterHasNPhotonHypothesis",     "std::vector<int>",    &m_eclPureClusterHasNPhotonHypothesis);
    m_tree->Branch("eclPureClusterHasNeutralHadronHypothesis",     "std::vector<int>",    &m_eclPureClusterHasNeutralHadronHypothesis);
  }
  m_tree->Branch("eclShowerMultip",     &m_eclShowerMultip,     "eclShowerMultip/I");
  m_tree->Branch("eclShowerIdx",     "std::vector<int>",       &m_eclShowerIdx);
  m_tree->Branch("eclShowerToMC1",      "std::vector<int>",       &m_eclShowerToMC1);
  m_tree->Branch("eclShowerToMCWeight1",      "std::vector<double>",       &m_eclShowerToMCWeight1);
  m_tree->Branch("eclShowerToMC1PDG",      "std::vector<int>",       &m_eclShowerToMC1PDG);
  m_tree->Branch("eclShowerToMC1Moth",      "std::vector<int>",       &m_eclShowerToMC1Moth);
  m_tree->Branch("eclShowerToMC1MothPDG",      "std::vector<int>",       &m_eclShowerToMC1MothPDG);
  m_tree->Branch("eclShowerToMC1GMoth",      "std::vector<int>",       &m_eclShowerToMC1GMoth);
  m_tree->Branch("eclShowerToMC1GMothPDG",      "std::vector<int>",       &m_eclShowerToMC1GMothPDG);
  m_tree->Branch("eclShowerToMC2",      "std::vector<int>",       &m_eclShowerToMC2);
  m_tree->Branch("eclShowerToMCWeight2",      "std::vector<double>",       &m_eclShowerToMCWeight2);
  m_tree->Branch("eclShowerToMC2PDG",      "std::vector<int>",       &m_eclShowerToMC2PDG);
  m_tree->Branch("eclShowerToMC2Moth",      "std::vector<int>",       &m_eclShowerToMC2Moth);
  m_tree->Branch("eclShowerToMC2MothPDG",      "std::vector<int>",       &m_eclShowerToMC2MothPDG);
  m_tree->Branch("eclShowerToMC2GMoth",      "std::vector<int>",       &m_eclShowerToMC2GMoth);
  m_tree->Branch("eclShowerToMC2GMothPDG",      "std::vector<int>",       &m_eclShowerToMC2GMothPDG);
  m_tree->Branch("eclShowerToMC3",      "std::vector<int>",       &m_eclShowerToMC3);
  m_tree->Branch("eclShowerToMCWeight3",      "std::vector<double>",       &m_eclShowerToMCWeight3);
  m_tree->Branch("eclShowerToMC3PDG",      "std::vector<int>",       &m_eclShowerToMC3PDG);
  m_tree->Branch("eclShowerToMC3Moth",      "std::vector<int>",       &m_eclShowerToMC3Moth);
  m_tree->Branch("eclShowerToMC3MothPDG",      "std::vector<int>",       &m_eclShowerToMC3MothPDG);
  m_tree->Branch("eclShowerToMC3GMoth",      "std::vector<int>",       &m_eclShowerToMC3GMoth);
  m_tree->Branch("eclShowerToMC3GMothPDG",      "std::vector<int>",       &m_eclShowerToMC3GMothPDG);
  m_tree->Branch("eclShowerToMC4",      "std::vector<int>",       &m_eclShowerToMC4);
  m_tree->Branch("eclShowerToMCWeight4",      "std::vector<double>",       &m_eclShowerToMCWeight4);
  m_tree->Branch("eclShowerToMC4PDG",      "std::vector<int>",       &m_eclShowerToMC4PDG);
  m_tree->Branch("eclShowerToMC4Moth",      "std::vector<int>",       &m_eclShowerToMC4Moth);
  m_tree->Branch("eclShowerToMC4MothPDG",      "std::vector<int>",       &m_eclShowerToMC4MothPDG);
  m_tree->Branch("eclShowerToMC4GMoth",      "std::vector<int>",       &m_eclShowerToMC4GMoth);
  m_tree->Branch("eclShowerToMC4GMothPDG",      "std::vector<int>",       &m_eclShowerToMC4GMothPDG);
  m_tree->Branch("eclShowerToMC5",      "std::vector<int>",       &m_eclShowerToMC5);
  m_tree->Branch("eclShowerToMCWeight5",      "std::vector<double>",       &m_eclShowerToMCWeight5);
  m_tree->Branch("eclShowerToMC5PDG",      "std::vector<int>",       &m_eclShowerToMC5PDG);
  m_tree->Branch("eclShowerToMC5Moth",      "std::vector<int>",       &m_eclShowerToMC5Moth);
  m_tree->Branch("eclShowerToMC5MothPDG",      "std::vector<int>",       &m_eclShowerToMC5MothPDG);
  m_tree->Branch("eclShowerToMC5GMoth",      "std::vector<int>",       &m_eclShowerToMC5GMoth);
  m_tree->Branch("eclShowerToMC5GMothPDG",      "std::vector<int>",       &m_eclShowerToMC5GMothPDG);
  m_tree->Branch("eclShowerToBkgWeight",      "std::vector<double>",       &m_eclShowerToBkgWeight);
  m_tree->Branch("eclShowerToLM1",      "std::vector<int>",       &m_eclShowerToLM1);
  m_tree->Branch("eclShowerToLM2",      "std::vector<int>",       &m_eclShowerToLM2);
  m_tree->Branch("eclShowerToLM3",      "std::vector<int>",       &m_eclShowerToLM3);
  m_tree->Branch("eclShowerToLM4",      "std::vector<int>",       &m_eclShowerToLM4);
  m_tree->Branch("eclShowerToLM5",      "std::vector<int>",       &m_eclShowerToLM5);
  m_tree->Branch("eclShowerSimHitSum",      "std::vector<double>",       &m_eclShowerSimHitSum);
  m_tree->Branch("eclShowerEnergy",     "std::vector<double>",    &m_eclShowerEnergy);
  m_tree->Branch("eclShowerUncEnergy",  "std::vector<double>",    &m_eclShowerUncEnergy);
  m_tree->Branch("eclShowerTheta",      "std::vector<double>",    &m_eclShowerTheta);
  m_tree->Branch("eclShowerPhi",        "std::vector<double>",    &m_eclShowerPhi);
  m_tree->Branch("eclShowerR",          "std::vector<double>",    &m_eclShowerR);
  m_tree->Branch("eclShowerNHits",   "std::vector<double>",       &m_eclShowerNHits);
  m_tree->Branch("eclShowerE9oE21",     "std::vector<double>",    &m_eclShowerE9oE21);
  m_tree->Branch("eclShowerTime",     "std::vector<double>",    &m_eclShowerTime);
  m_tree->Branch("eclShowerT99",     "std::vector<double>",    &m_eclShowerT99);
  m_tree->Branch("eclShowerConnectedRegionId",     "std::vector<int>",      &m_eclShowerConnectedRegionId);
  m_tree->Branch("eclShowerHypothesisId",     "std::vector<int>",    &m_eclShowerHypothesisId);
  m_tree->Branch("eclShowerCentralCellId",     "std::vector<int>",      &m_eclShowerCentralCellId);
  m_tree->Branch("eclShowerEnergyError",     "std::vector<double>",       &m_eclShowerEnergyError);
  m_tree->Branch("eclShowerThetaError",     "std::vector<double>",      &m_eclShowerThetaError);
  m_tree->Branch("eclShowerPhiError",     "std::vector<double>",      &m_eclShowerPhiError);
  m_tree->Branch("eclShowerTimeResolution",     "std::vector<double>",      &m_eclShowerTimeResolution);
  m_tree->Branch("eclShowerHighestEnergy",     "std::vector<double>",      &m_eclShowerHighestEnergy);
  m_tree->Branch("eclShowerLateralEnergy",     "std::vector<double>",      &m_eclShowerLateralEnergy);
  m_tree->Branch("eclShowerMinTrkDistance",     "std::vector<double>",    &m_eclShowerMinTrkDistance);
  m_tree->Branch("eclShowerTrkDepth",     "std::vector<double>",    &m_eclShowerTrkDepth);
  m_tree->Branch("eclShowerShowerDepth",     "std::vector<double>",      &m_eclShowerShowerDepth);
  m_tree->Branch("eclShowerAbsZernike40",     "std::vector<double>",       &m_eclShowerAbsZernike40);
  m_tree->Branch("eclShowerAbsZernike51",     "std::vector<double>",       &m_eclShowerAbsZernike51);
  m_tree->Branch("eclShowerZernikeMVA",     "std::vector<double>",       &m_eclShowerZernikeMVA);
  m_tree->Branch("eclShowerSecondMoment",     "std::vector<double>",        &m_eclShowerSecondMoment);
  m_tree->Branch("eclShowerE1oE9",     "std::vector<double>",      &m_eclShowerE1oE9);
  m_tree->Branch("eclShowerIsTrack",     "std::vector<int>",    &m_eclShowerIsTrack);
  m_tree->Branch("eclShowerIsCluster",     "std::vector<bool>",        &m_eclShowerIsCluster);
  m_tree->Branch("eclShowerMCVtxInEcl",     "std::vector<int>",    &m_eclShowerMCVtxInEcl);
  m_tree->Branch("eclShowerMCFlightMatch",     "std::vector<int>",    &m_eclShowerMCFlightMatch);
  m_tree->Branch("eclShowerMCFFlightMatch",     "std::vector<int>",    &m_eclShowerMCFFlightMatch);
  m_tree->Branch("eclShowerHighestE1mE2",     "std::vector<double>",    &m_eclShowerHighestE1mE2);
  m_tree->Branch("eclShowerNumberOfCrystalsForEnergy",     "std::vector<double>",    &m_eclShowerNumberOfCrystalsForEnergy);
 
  //MC Truth
  m_tree->Branch("mcMultip",     &m_mcMultip,           "mcMultip/I");
  m_tree->Branch("mcIdx",        "std::vector<int>",    &m_mcIdx);
  m_tree->Branch("mcPdg",        "std::vector<int>",    &m_mcPdg);
  m_tree->Branch("mcMothPdg",    "std::vector<int>",    &m_mcMothPdg);
  m_tree->Branch("mcGMothPdg",   "std::vector<int>",    &m_mcGMothPdg);
  m_tree->Branch("mcGGMothPdg",  "std::vector<int>",    &m_mcGGMothPdg);
  m_tree->Branch("mcEnergy",     "std::vector<double>", &m_mcEnergy);
  m_tree->Branch("mcPx",         "std::vector<double>", &m_mcPx);
  m_tree->Branch("mcPy",         "std::vector<double>", &m_mcPy);
  m_tree->Branch("mcPz",         "std::vector<double>", &m_mcPz);
  m_tree->Branch("mcDecVtxx",    "std::vector<double>", &m_mcDecayVtxX);
  m_tree->Branch("mcDecVtxy",    "std::vector<double>", &m_mcDecayVtxY);
  m_tree->Branch("mcDecVtxz",    "std::vector<double>", &m_mcDecayVtxZ);
  m_tree->Branch("mcProdVtxx",   "std::vector<double>", &m_mcProdVtxX);
  m_tree->Branch("mcProdVtxy",   "std::vector<double>", &m_mcProdVtxY);
  m_tree->Branch("mcProdVtxz",   "std::vector<double>", &m_mcProdVtxZ);
  m_tree->Branch("mcSecProc",    "std::vector<int>",    &m_mcSecondaryPhysProc);
 
  if (m_doTracking == true) {
    m_tree->Branch("trkMultip",     &m_trkMultip,          "trkMulti/I");
    m_tree->Branch("trkIdx",     "std::vector<int>",       &m_trkIdx);
    m_tree->Branch("trkPdg",        "std::vector<int>",    &m_trkPdg);
    m_tree->Branch("trkCharge",        "std::vector<int>",    &m_trkCharge);
    m_tree->Branch("trkPx",         "std::vector<double>", &m_trkPx);
    m_tree->Branch("trkPy",         "std::vector<double>", &m_trkPy);
    m_tree->Branch("trkPz",         "std::vector<double>", &m_trkPz);
    m_tree->Branch("trkP",         "std::vector<double>", &m_trkP);
    m_tree->Branch("trkTheta",         "std::vector<double>", &m_trkTheta);
    m_tree->Branch("trkPhi",         "std::vector<double>", &m_trkPhi);
    m_tree->Branch("trkPosx",       "std::vector<double>", &m_trkX);
    m_tree->Branch("trkPosy",       "std::vector<double>", &m_trkY);
    m_tree->Branch("trkPosz",      "std::vector<double>",  &m_trkZ);
 
    m_tree->Branch("eclpidtrkIdx",     "std::vector<int>",  &m_eclpidtrkIdx);
    m_tree->Branch("eclpidEnergy",     "std::vector<double>",  &m_eclpidEnergy);
    m_tree->Branch("eclpidEop",        "std::vector<double>",  &m_eclpidEop);
    m_tree->Branch("eclpidE9E21",      "std::vector<double>",  &m_eclpidE9E21);
    m_tree->Branch("eclpidNCrystals",  "std::vector<int>",  &m_eclpidNCrystals);
    m_tree->Branch("eclpidNClusters",  "std::vector<int>",  &m_eclpidNClusters);
    m_tree->Branch("eclLogLikeEl",      "std::vector<double>",  &m_eclLogLikeEl);
    m_tree->Branch("eclLogLikeMu",      "std::vector<double>",  &m_eclLogLikeMu);
    m_tree->Branch("eclLogLikePi",      "std::vector<double>",  &m_eclLogLikePi);
  }
 
  B2INFO("[ECLDataAnalysis Module]: Initialization of ECLDataAnalysis Module completed.");
 
}

◆ setAbortLevel()

void setAbortLevel ( int abortLevel )

inherited

Configure the abort log level.

Definition at line 67 of file Module.cc.

{
  m_logConfig.setAbortLevel(static_cast<LogConfig::ELogLevel>(abortLevel));
}

◆ setDebugLevel()

void setDebugLevel ( int debugLevel )

inherited

Configure the debug messaging level.

Definition at line 61 of file Module.cc.

{
  m_logConfig.setDebugLevel(debugLevel);
}

◆ setDescription()

void setDescription ( const std::string & description )

protectedinherited

Sets the description of the module.

Parameters

description A description of the module.

Definition at line 214 of file Module.cc.

{
  m_description = description;
}

◆ setLogConfig()

void setLogConfig ( const LogConfig & logConfig )

inlineinherited

Set the log system configuration.

Definition at line 229 of file Module.h.

229{m_logConfig = logConfig;}

◆ setLogInfo()

void setLogInfo	(	int	logLevel,
		unsigned int	logInfo )

inherited

Configure the printed log information for the given level.

Parameters

logLevel	The log level (one of LogConfig::ELogLevel)
logInfo	What kind of info should be printed? ORed combination of LogConfig::ELogInfo flags.

Definition at line 73 of file Module.cc.

{
  m_logConfig.setLogInfo(static_cast<LogConfig::ELogLevel>(logLevel), logInfo);
}

◆ setLogLevel()

void setLogLevel ( int logLevel )

inherited

Configure the log level.

Definition at line 55 of file Module.cc.

{
  m_logConfig.setLogLevel(static_cast<LogConfig::ELogLevel>(logLevel));
}

◆ setName()

void setName ( const std::string & name )

inlineinherited

Set the name of the module.

Note: The module name is set when using the REG_MODULE macro, but the module can be renamed before calling process() using the set_name() function in your steering file.

Parameters

name	The name of the module

Definition at line 213 of file Module.h.

213{ m_name = name; };

◆ setParamList()

void setParamList ( const ModuleParamList & params )

inlineprotectedinherited

Replace existing parameter list.

Definition at line 500 of file Module.h.

500{ m_moduleParamList = params; }

◆ setParamPython()

void setParamPython	(	const std::string &	name,
		const boost::python::object &	pyObj )

privateinherited

Implements a method for setting boost::python objects.

The method supports the following types: list, dict, int, double, string, bool The conversion of the python object to the C++ type and the final storage of the parameter value is done in the ModuleParam class.

Parameters

name	The unique name of the parameter.
pyObj	The object which should be converted and stored as the parameter value.

Definition at line 234 of file Module.cc.

{
  LogSystem& logSystem = LogSystem::Instance();
  logSystem.updateModule(&(getLogConfig()), getName());
  try {
    m_moduleParamList.setParamPython(name, pyObj);
  } catch (std::runtime_error& e) {
    throw std::runtime_error("Cannot set parameter '" + name + "' for module '"
                             + m_name + "': " + e.what());
  }
 
  logSystem.updateModule(nullptr);
}

◆ setParamPythonDict()

void setParamPythonDict ( const boost::python::dict & dictionary )

privateinherited

Implements a method for reading the parameter values from a boost::python dictionary.

The key of the dictionary has to be the name of the parameter and the value has to be of one of the supported parameter types.

Parameters

dictionary The python dictionary from which the parameter values are read.

Definition at line 249 of file Module.cc.

{
 
  LogSystem& logSystem = LogSystem::Instance();
  logSystem.updateModule(&(getLogConfig()), getName());
 
  boost::python::list dictKeys = dictionary.keys();
  int nKey = boost::python::len(dictKeys);
 
  //Loop over all keys in the dictionary
  for (int iKey = 0; iKey < nKey; ++iKey) {
    boost::python::object currKey = dictKeys[iKey];
    boost::python::extract<std::string> keyProxy(currKey);
 
    if (keyProxy.check()) {
      const boost::python::object& currValue = dictionary[currKey];
      setParamPython(keyProxy, currValue);
    } else {
      B2ERROR("Setting the module parameters from a python dictionary: invalid key in dictionary!");
    }
  }
 
  logSystem.updateModule(nullptr);
}

◆ setPropertyFlags()

void setPropertyFlags ( unsigned int propertyFlags )

inherited

Sets the flags for the module properties.

Parameters

propertyFlags bitwise OR of EModulePropFlags

Definition at line 208 of file Module.cc.

{
  m_propertyFlags = propertyFlags;
}

◆ setReturnValue() [1/2]

void setReturnValue ( bool value )

protectedinherited

Sets the return value for this module as bool.

The bool value is saved as an integer with the convention 1 meaning true and 0 meaning false. The value can be used in the steering file to divide the analysis chain into several paths.

Parameters

value The value of the return value.

Definition at line 227 of file Module.cc.

{
  m_hasReturnValue = true;
  m_returnValue = value;
}

◆ setReturnValue() [2/2]

void setReturnValue ( int value )

protectedinherited

Sets the return value for this module as integer.

The value can be used in the steering file to divide the analysis chain into several paths.

Parameters

value The value of the return value.

Definition at line 220 of file Module.cc.

{
  m_hasReturnValue = true;
  m_returnValue = value;
}

◆ setType()

void setType ( const std::string & type )

protectedinherited

Set the module type.

Only for use by internal modules (which don't use the normal REG_MODULE mechanism).

Definition at line 48 of file Module.cc.

{
  if (!m_type.empty())
    B2FATAL("Trying to change module type from " << m_type << " is not allowed, the value is assumed to be fixed.");
  m_type = type;
}

◆ terminate()

void terminate ( void )

overridevirtual

Termination action.

Reimplemented from Module.

Definition at line 2380 of file ECLDataAnalysisModule.cc.

{
 
  if (m_rootFilePtr != nullptr) {
    m_rootFilePtr->cd();
    m_tree->Write();
  }
 
}

Member Data Documentation

◆ m_conditions

std::vector<ModuleCondition> m_conditions

privateinherited

Module condition, only non-null if set.

Definition at line 520 of file Module.h.

◆ m_description

std::string m_description

privateinherited

The description of the module.

Definition at line 510 of file Module.h.

◆ m_doDigits

bool m_doDigits

private

if true, info on Digits and CalDigits is stored

Definition at line 86 of file ECLDataAnalysisModule.h.

◆ m_doHits

bool m_doHits

private

if true, info on Hits and SimHits is stored

Definition at line 85 of file ECLDataAnalysisModule.h.

◆ m_doPureCsI

bool m_doPureCsI

private

if true, info on pureCsI upgrade is stored

Definition at line 84 of file ECLDataAnalysisModule.h.

◆ m_doTracking

bool m_doTracking

private

if true, info on tracking will be stored, job will fail if doTracking==1 and the tracking modules are not enabled at python level

Definition at line 83 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitAmp

std::vector<double>* m_eclCalDigitAmp

private

ECLCalDigit amplitude.

Definition at line 215 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitCellId

std::vector<int>* m_eclCalDigitCellId

private

Number of ECLCalDigit CellId.

Definition at line 214 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitFitQuality

std::vector<int>* m_eclCalDigitFitQuality

private

ECLCalDigit fit quality.

Definition at line 217 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitIdx

std::vector<int>* m_eclCalDigitIdx

private

ECLCalDigit index.

Definition at line 194 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitMultip

int m_eclCalDigitMultip

private

Number of ECLCalDigits per event.

Definition at line 193 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigits

StoreArray<ECLCalDigit> m_eclCalDigits

private

Store array: ECLCalDigit.

Definition at line 95 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitSimHitSum

std::vector<double>* m_eclCalDigitSimHitSum

private

Full energy contribution related to ECLCalDigit.

Definition at line 212 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitTimeFit

std::vector<double>* m_eclCalDigitTimeFit

private

ECLCalDigit timing.

Definition at line 216 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToBkgWeight

std::vector<double>* m_eclCalDigitToBkgWeight

private

Remaining energy contribution not associated to first five MCParticles related to ECLCalDigit.

Definition at line 211 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToCR

std::vector<int>* m_eclCalDigitToCR

private

Index of CR related to that ECLCalDigit.

Definition at line 218 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToLM

std::vector<int>* m_eclCalDigitToLM

private

Index of LM related to that ECLCalDigit.

Definition at line 219 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMC1

std::vector<int>* m_eclCalDigitToMC1

private

Index of first MCParticle related to ECLCalDigit.

Definition at line 195 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMC1PDG

std::vector<int>* m_eclCalDigitToMC1PDG

private

PDG code of first MCParticle related to ECLCalDigit.

Definition at line 196 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMC2

std::vector<int>* m_eclCalDigitToMC2

private

Index of second MCParticle related to ECLCalDigit.

Definition at line 198 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMC2PDG

std::vector<int>* m_eclCalDigitToMC2PDG

private

PDG code of second MCParticle related to ECLCalDigit.

Definition at line 199 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMC3

std::vector<int>* m_eclCalDigitToMC3

private

Index of third MCParticle related to ECLCalDigit.

Definition at line 201 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMC3PDG

std::vector<int>* m_eclCalDigitToMC3PDG

private

PDG code of third MCParticle related to ECLCalDigit.

Definition at line 202 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMC4

std::vector<int>* m_eclCalDigitToMC4

private

Index of fourth MCParticle related to ECLCalDigit.

Definition at line 204 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMC4PDG

std::vector<int>* m_eclCalDigitToMC4PDG

private

PDG code of fourth MCParticle related to ECLCalDigit.

Definition at line 205 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMC5

std::vector<int>* m_eclCalDigitToMC5

private

Index of fifth MCParticle related to ECLCalDigit.

Definition at line 207 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMC5PDG

std::vector<int>* m_eclCalDigitToMC5PDG

private

PDG code of fifth MCParticle related to ECLCalDigit.

Definition at line 208 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMCWeight1

std::vector<double>* m_eclCalDigitToMCWeight1

private

Energy contribution of first MCParticle related to ECLCalDigit.

Definition at line 197 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMCWeight2

std::vector<double>* m_eclCalDigitToMCWeight2

private

Energy contribution of second MCParticle related to ECLCalDigit.

Definition at line 200 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMCWeight3

std::vector<double>* m_eclCalDigitToMCWeight3

private

Energy contribution of third MCParticle related to ECLCalDigit.

Definition at line 203 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMCWeight4

std::vector<double>* m_eclCalDigitToMCWeight4

private

Energy contribution of fourth MCParticle related to ECLCalDigit.

Definition at line 206 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToMCWeight5

std::vector<double>* m_eclCalDigitToMCWeight5

private

Energy contribution of fifth MCParticle related to ECLCalDigit.

Definition at line 209 of file ECLDataAnalysisModule.h.

◆ m_eclCalDigitToShower

std::vector<int>* m_eclCalDigitToShower

private

Index of ECLShower related to that ECLCalDigit.

Definition at line 213 of file ECLDataAnalysisModule.h.

◆ m_eclClusterAbsZernike40

std::vector<double>* m_eclClusterAbsZernike40

private

Reconstructed Abs Zernike40.

Definition at line 301 of file ECLDataAnalysisModule.h.

◆ m_eclClusterAbsZernike51

std::vector<double>* m_eclClusterAbsZernike51

private

Reconstructed Abs Zernike51.

Definition at line 302 of file ECLDataAnalysisModule.h.

◆ m_eclClusterCellId

std::vector<int>* m_eclClusterCellId

private

CellId with highest energy deposit in ECLCluster.

Definition at line 295 of file ECLDataAnalysisModule.h.

◆ m_eclClusterClosestTrackDist

std::vector<double>* m_eclClusterClosestTrackDist

private

Flag for charged clusters.

Definition at line 299 of file ECLDataAnalysisModule.h.

◆ m_eclClusterCrystalHealth

std::vector<int>* m_eclClusterCrystalHealth

private

Crystal health flag.

Definition at line 297 of file ECLDataAnalysisModule.h.

◆ m_eclClusterDeltaL

std::vector<double>* m_eclClusterDeltaL

private

Reconstructed Cluster DeltaL.

Definition at line 300 of file ECLDataAnalysisModule.h.

◆ m_eclClusterDeltaTime99

std::vector<double>* m_eclClusterDeltaTime99

private

DeltaTime99.

Definition at line 307 of file ECLDataAnalysisModule.h.

◆ m_eclClusterDetectorRegion

std::vector<int>* m_eclClusterDetectorRegion

private

Cluster Detector Region.

Definition at line 308 of file ECLDataAnalysisModule.h.

◆ m_eclClusterE1oE9

std::vector<double>* m_eclClusterE1oE9

private

Reconstructed E1 over E9.

Definition at line 304 of file ECLDataAnalysisModule.h.

◆ m_eclClusterE9oE21

std::vector<double>* m_eclClusterE9oE21

private

Ratio of 3x3 over 5x5 crystal matrices energies for ECLCluster.

Definition at line 293 of file ECLDataAnalysisModule.h.

◆ m_eclClusterEnergy

std::vector<double>* m_eclClusterEnergy

private

ECLCluster energy.

Definition at line 283 of file ECLDataAnalysisModule.h.

◆ m_eclClusterEnergyDepSum

std::vector<double>* m_eclClusterEnergyDepSum

private

ECLCluster simulated energy.

Definition at line 290 of file ECLDataAnalysisModule.h.

◆ m_eclClusterEnergyError

std::vector<double>* m_eclClusterEnergyError

private

ECLCluster energy error.

Definition at line 284 of file ECLDataAnalysisModule.h.

◆ m_eclClusterGammaMultip

int m_eclClusterGammaMultip

private

Number of ECLClusters per event.

Definition at line 262 of file ECLDataAnalysisModule.h.

◆ m_eclClusterHasNeutralHadronHypothesis

std::vector<int>* m_eclClusterHasNeutralHadronHypothesis

private

Cluster has neutral hadron hypothesis.

Definition at line 310 of file ECLDataAnalysisModule.h.

◆ m_eclClusterHasNPhotonHypothesis

std::vector<int>* m_eclClusterHasNPhotonHypothesis

private

Cluster has n-photon hypothesis.

Definition at line 309 of file ECLDataAnalysisModule.h.

◆ m_eclClusterHighestE

std::vector<double>* m_eclClusterHighestE

private

Highest energy deposit (per crystal) in ECLCluster.

Definition at line 294 of file ECLDataAnalysisModule.h.

◆ m_eclClusterIdx

std::vector<int>* m_eclClusterIdx

private

ECLCluster index.

Definition at line 263 of file ECLDataAnalysisModule.h.

◆ m_eclClusterIsTrack

std::vector<bool>* m_eclClusterIsTrack

private

Flag for charged clusters.

Definition at line 298 of file ECLDataAnalysisModule.h.

◆ m_eclClusterLAT

std::vector<double>* m_eclClusterLAT

private

Reconstructed LAT.

Definition at line 306 of file ECLDataAnalysisModule.h.

◆ m_eclClusterMultip

int m_eclClusterMultip

private

Number of ECLClusters per event.

Definition at line 260 of file ECLDataAnalysisModule.h.

◆ m_eclClusterNofCrystals

std::vector<int>* m_eclClusterNofCrystals

private

Number of crystals in ECLCluster.

Definition at line 296 of file ECLDataAnalysisModule.h.

◆ m_eclClusterPhi

std::vector<double>* m_eclClusterPhi

private

ECLCluster azimuthal direction.

Definition at line 287 of file ECLDataAnalysisModule.h.

◆ m_eclClusterPhiError

std::vector<double>* m_eclClusterPhiError

private

ECLCluster error on azimuthal direction.

Definition at line 288 of file ECLDataAnalysisModule.h.

◆ m_eclClusterR

std::vector<double>* m_eclClusterR

private

ECLCluster distance from IP.

Definition at line 289 of file ECLDataAnalysisModule.h.

◆ m_eclClusters

StoreArray<ECLCluster> m_eclClusters

private

Store array: ECLCluster.

Definition at line 98 of file ECLDataAnalysisModule.h.

◆ m_eclClusterSecondMoment

std::vector<double>* m_eclClusterSecondMoment

private

Reconstructed Second Moment.

Definition at line 305 of file ECLDataAnalysisModule.h.

◆ m_eclClusterSimHitSum

std::vector<double>* m_eclClusterSimHitSum

private

Energy contribution of 1st MCParticle related to ECLCluster.

Definition at line 281 of file ECLDataAnalysisModule.h.

◆ m_eclClusterTheta

std::vector<double>* m_eclClusterTheta

private

ECLCluster polar direction.

Definition at line 285 of file ECLDataAnalysisModule.h.

◆ m_eclClusterThetaError

std::vector<double>* m_eclClusterThetaError

private

ECLCluster error on polar direction.

Definition at line 286 of file ECLDataAnalysisModule.h.

◆ m_eclClusterTiming

std::vector<double>* m_eclClusterTiming

private

ECLCluster time.

Definition at line 291 of file ECLDataAnalysisModule.h.

◆ m_eclClusterTimingError

std::vector<double>* m_eclClusterTimingError

private

ECLCluster time error.

Definition at line 292 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToBkgWeight

std::vector<double>* m_eclClusterToBkgWeight

private

Remaining energy contribution not associated to first five MCParticles related to ECLCluster.

Definition at line 280 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMC1

std::vector<int>* m_eclClusterToMC1

private

Index of first MCParticle related to ECLCluster.

Definition at line 264 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMC1PDG

std::vector<int>* m_eclClusterToMC1PDG

private

PDG code of first MCParticle related to ECLCluster.

Definition at line 266 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMC2

std::vector<int>* m_eclClusterToMC2

private

Index of second MCParticle related to ECLCluster.

Definition at line 267 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMC2PDG

std::vector<int>* m_eclClusterToMC2PDG

private

PDG code of second MCParticle related to ECLCluster.

Definition at line 269 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMC3

std::vector<int>* m_eclClusterToMC3

private

Index of third MCParticle related to ECLCluster.

Definition at line 270 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMC3PDG

std::vector<int>* m_eclClusterToMC3PDG

private

PDG code of third MCParticle related to ECLCluster.

Definition at line 272 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMC4

std::vector<int>* m_eclClusterToMC4

private

Index of fourth MCParticle related to ECLCluster.

Definition at line 273 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMC4PDG

std::vector<int>* m_eclClusterToMC4PDG

private

PDG code of fourth MCParticle related to ECLCluster.

Definition at line 275 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMC5

std::vector<int>* m_eclClusterToMC5

private

Index of fifth MCParticle related to ECLCluster.

Definition at line 276 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMC5PDG

std::vector<int>* m_eclClusterToMC5PDG

private

PDG code of fifth MCParticle related to ECLCluster.

Definition at line 278 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMCWeight1

std::vector<double>* m_eclClusterToMCWeight1

private

Energy contribution of first MCParticle related to ECLCluster.

Definition at line 265 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMCWeight2

std::vector<double>* m_eclClusterToMCWeight2

private

Energy contribution of second MCParticle related to ECLCluster.

Definition at line 268 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMCWeight3

std::vector<double>* m_eclClusterToMCWeight3

private

Energy contribution of third MCParticle related to ECLCluster.

Definition at line 271 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMCWeight4

std::vector<double>* m_eclClusterToMCWeight4

private

Energy contribution of fourth MCParticle related to ECLCluster.

Definition at line 274 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToMCWeight5

std::vector<double>* m_eclClusterToMCWeight5

private

Energy contribution of 5th MCParticle related to ECLCluster.

Definition at line 277 of file ECLDataAnalysisModule.h.

◆ m_eclClusterToShower

std::vector<int>* m_eclClusterToShower

private

Index of ECLShower related to ECLCluster.

Definition at line 282 of file ECLDataAnalysisModule.h.

◆ m_eclClusterTrueMultip

int m_eclClusterTrueMultip

private

Number of ECLClusters per event.

Definition at line 261 of file ECLDataAnalysisModule.h.

◆ m_eclClusterZernikeMVA

std::vector<double>* m_eclClusterZernikeMVA

private

Zernike MVA.

Definition at line 303 of file ECLDataAnalysisModule.h.

◆ m_eclConnectedRegions

StoreArray<ECLConnectedRegion> m_eclConnectedRegions

private

Store array: ECLConnectedRegion.

Definition at line 96 of file ECLDataAnalysisModule.h.

◆ m_eclCRIdx

std::vector<int>* m_eclCRIdx

private

Connected Region ID.

Definition at line 221 of file ECLDataAnalysisModule.h.

◆ m_eclCRIsTrack

std::vector<int>* m_eclCRIsTrack

private

Int for Connected Region - Track Match.

Definition at line 222 of file ECLDataAnalysisModule.h.

◆ m_eclCRLikelihoodChargedHadron

std::vector<double>* m_eclCRLikelihoodChargedHadron

private

Connected Region Charged Hadron Likelihood.

Definition at line 224 of file ECLDataAnalysisModule.h.

◆ m_eclCRLikelihoodElectronNGamma

std::vector<double>* m_eclCRLikelihoodElectronNGamma

private

Connected Region Electron Likelihood.

Definition at line 225 of file ECLDataAnalysisModule.h.

◆ m_eclCRLikelihoodMergedPi0

std::vector<double>* m_eclCRLikelihoodMergedPi0

private

Connected Region Merged Pi0 Likelihood.

Definition at line 228 of file ECLDataAnalysisModule.h.

◆ m_eclCRLikelihoodMIPNGamma

std::vector<double>* m_eclCRLikelihoodMIPNGamma

private

Connected Region MIP Likelihood.

Definition at line 223 of file ECLDataAnalysisModule.h.

◆ m_eclCRLikelihoodNeutralHadron

std::vector<double>* m_eclCRLikelihoodNeutralHadron

private

Connected Region Neutral Hadron Likelihood.

Definition at line 227 of file ECLDataAnalysisModule.h.

◆ m_eclCRLikelihoodNGamma

std::vector<double>* m_eclCRLikelihoodNGamma

private

Connected Region Gamma Likelihood.

Definition at line 226 of file ECLDataAnalysisModule.h.

◆ m_eclDigitAmp

std::vector<int>* m_eclDigitAmp

private

ECLDigit amplitude.

Definition at line 188 of file ECLDataAnalysisModule.h.

◆ m_eclDigitCellId

std::vector<int>* m_eclDigitCellId

private

Number of ECLDigit CellId.

Definition at line 187 of file ECLDataAnalysisModule.h.

◆ m_eclDigitFitQuality

std::vector<int>* m_eclDigitFitQuality

private

ECLDigit fit quality.

Definition at line 190 of file ECLDataAnalysisModule.h.

◆ m_eclDigitIdx

std::vector<int>* m_eclDigitIdx

private

ECLDigit index.

Definition at line 185 of file ECLDataAnalysisModule.h.

◆ m_eclDigitMultip

int m_eclDigitMultip

private

Number of ECLDigits per event.

Definition at line 184 of file ECLDataAnalysisModule.h.

◆ m_eclDigits

StoreArray<ECLDigit> m_eclDigits

private

Store array: ECLDigit.

Definition at line 94 of file ECLDataAnalysisModule.h.

◆ m_eclDigitTimeFit

std::vector<int>* m_eclDigitTimeFit

private

ECLDigit timing.

Definition at line 189 of file ECLDataAnalysisModule.h.

◆ m_eclDigitToCalDigit

std::vector<int>* m_eclDigitToCalDigit

private

Index of CalDigit related to that ECLDigit.

Definition at line 191 of file ECLDataAnalysisModule.h.

◆ m_eclDigitToMC

std::vector<int>* m_eclDigitToMC

private

Index of MCParticle related to that ECLDigit.

Definition at line 186 of file ECLDataAnalysisModule.h.

◆ m_eclHitCellId

std::vector<int>* m_eclHitCellId

private

ECLHit CellID.

Definition at line 256 of file ECLDataAnalysisModule.h.

◆ m_eclHitEnergyDep

std::vector<double>* m_eclHitEnergyDep

private

ECLHit energy.

Definition at line 257 of file ECLDataAnalysisModule.h.

◆ m_eclHitIdx

std::vector<int>* m_eclHitIdx

private

Index of ECLHits.

Definition at line 250 of file ECLDataAnalysisModule.h.

◆ m_eclHitMultip

int m_eclHitMultip

private

Number of ECLHits per event.

Definition at line 249 of file ECLDataAnalysisModule.h.

◆ m_eclHits

StoreArray<ECLHit> m_eclHits

private

Store array: ECLHit.

Definition at line 93 of file ECLDataAnalysisModule.h.

◆ m_eclHitTimeAve

std::vector<double>* m_eclHitTimeAve

private

ECLHit time.

Definition at line 258 of file ECLDataAnalysisModule.h.

◆ m_eclHitToDigit

std::vector<int>* m_eclHitToDigit

private

Index of ECLDigit related to ECLHit.

Definition at line 252 of file ECLDataAnalysisModule.h.

◆ m_eclHitToDigitAmp

std::vector<int>* m_eclHitToDigitAmp

private

Amplitude of ECLDigit related to ECLHit.

Definition at line 253 of file ECLDataAnalysisModule.h.

◆ m_eclHitToMC

std::vector<int>* m_eclHitToMC

private

Index of MCParticle related to ECLHit.

Definition at line 251 of file ECLDataAnalysisModule.h.

◆ m_eclHitToPureDigit

std::vector<int>* m_eclHitToPureDigit

private

Index of ECLDigit related to ECLHit, PureCsI option.

Definition at line 254 of file ECLDataAnalysisModule.h.

◆ m_eclHitToPureDigitAmp

std::vector<int>* m_eclHitToPureDigitAmp

private

Amplitude of ECLDigit related to ECLHit, PureCsI option.

Definition at line 255 of file ECLDataAnalysisModule.h.

◆ m_eclLMCellId

std::vector<int>* m_eclLMCellId

private

Local Maximum Cell ID.

Definition at line 233 of file ECLDataAnalysisModule.h.

◆ m_eclLMId

std::vector<int>* m_eclLMId

private

Local Maximum ID.

Definition at line 231 of file ECLDataAnalysisModule.h.

◆ m_eclLMMultip

int m_eclLMMultip

private

Local Maxima multiplicity.

Definition at line 230 of file ECLDataAnalysisModule.h.

◆ m_eclLMType

std::vector<int>* m_eclLMType

private

Local Maximum type.

Definition at line 232 of file ECLDataAnalysisModule.h.

◆ m_eclLocalMaximums

StoreArray<ECLLocalMaximum> m_eclLocalMaximums

private

Store array: ECLLocalMaximum.

Definition at line 99 of file ECLDataAnalysisModule.h.

◆ m_eclLogLikeEl

std::vector<double>* m_eclLogLikeEl

private

PID track electron likelihood.

Definition at line 539 of file ECLDataAnalysisModule.h.

◆ m_eclLogLikeMu

std::vector<double>* m_eclLogLikeMu

private

PID track muon likelihood.

Definition at line 540 of file ECLDataAnalysisModule.h.

◆ m_eclLogLikePi

std::vector<double>* m_eclLogLikePi

private

PID track pion likelihood.

Definition at line 541 of file ECLDataAnalysisModule.h.

◆ m_eclpidE9E21

std::vector<double>* m_eclpidE9E21

private

PID track ration of 3x3 over 5x5 crystal matrices energies.

Definition at line 536 of file ECLDataAnalysisModule.h.

◆ m_eclpidEnergy

std::vector<double>* m_eclpidEnergy

private

PID track energy.

Definition at line 534 of file ECLDataAnalysisModule.h.

◆ m_eclpidEop

std::vector<double>* m_eclpidEop

private

PID track E/p.

Definition at line 535 of file ECLDataAnalysisModule.h.

◆ m_eclPidLikelihoods

StoreArray<ECLPidLikelihood> m_eclPidLikelihoods

private

ECLPidLikelihood storeArray.

Definition at line 91 of file ECLDataAnalysisModule.h.

◆ m_eclpidNClusters

std::vector<int>* m_eclpidNClusters

private

PID track number of clusters.

Definition at line 538 of file ECLDataAnalysisModule.h.

◆ m_eclpidNCrystals

std::vector<int>* m_eclpidNCrystals

private

PID track number of crystals.

Definition at line 537 of file ECLDataAnalysisModule.h.

◆ m_eclpidtrkIdx

std::vector<int>* m_eclpidtrkIdx

private

PID track index.

Definition at line 533 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitAmp

std::vector<double>* m_eclPureCalDigitAmp

private

CalDigit amplitude, PureCsI option.

Definition at line 348 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitCellId

std::vector<int>* m_eclPureCalDigitCellId

private

Number of CalDigit CellId, PureCsI option.

Definition at line 347 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitFitQuality

std::vector<int>* m_eclPureCalDigitFitQuality

private

CalDigit fit quality, PureCsI option.

Definition at line 350 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitIdx

std::vector<int>* m_eclPureCalDigitIdx

private

ECLCalDigit index, PureCsI option.

Definition at line 322 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitMultip

int m_eclPureCalDigitMultip

private

Number of ECLCalDigits per event, PureCsI option.

Definition at line 321 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigits

StoreArray<ECLCalDigit> m_eclPureCalDigits

private

Store array: ECLPureCalDigit.

Definition at line 132 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitSimHitSum

std::vector<double>* m_eclPureCalDigitSimHitSum

private

Full energy contribution related to CalDigit, PureCsI option.

Definition at line 345 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitTimeFit

std::vector<double>* m_eclPureCalDigitTimeFit

private

CalDigit timing, PureCsI option.

Definition at line 349 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToBkgWeight

std::vector<double>* m_eclPureCalDigitToBkgWeight

private

Remaining energy contribution not associated to first five MCParticles related to CalDigit, PureCsI option.

Definition at line 344 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToCR

std::vector<int>* m_eclPureCalDigitToCR

private

Index of CR related to that CalDigit, PureCsI option.

Definition at line 351 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToLM

std::vector<int>* m_eclPureCalDigitToLM

private

Index of LM related to that CalDigit, PureCsI option.

Definition at line 352 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMC1

std::vector<int>* m_eclPureCalDigitToMC1

private

Index of first MCParticle related to CalDigit, PureCsI option.

Definition at line 323 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMC1PDG

std::vector<int>* m_eclPureCalDigitToMC1PDG

private

PDG code of first MCParticle related to CalDigit, PureCsI option.

Definition at line 324 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMC2

std::vector<int>* m_eclPureCalDigitToMC2

private

Index of second MCParticle related to CalDigit, PureCsI option.

Definition at line 327 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMC2PDG

std::vector<int>* m_eclPureCalDigitToMC2PDG

private

PDG code of second MCParticle related to CalDigit, PureCsI option.

Definition at line 328 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMC3

std::vector<int>* m_eclPureCalDigitToMC3

private

Index of third MCParticle related to CalDigit, PureCsI option.

Definition at line 331 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMC3PDG

std::vector<int>* m_eclPureCalDigitToMC3PDG

private

PDG code of third MCParticle related to CalDigit, PureCsI option.

Definition at line 332 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMC4

std::vector<int>* m_eclPureCalDigitToMC4

private

Index of fourth MCParticle related to CalDigit, PureCsI option.

Definition at line 335 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMC4PDG

std::vector<int>* m_eclPureCalDigitToMC4PDG

private

PDG code of fourth MCParticle related to CalDigit, PureCsI option.

Definition at line 336 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMC5

std::vector<int>* m_eclPureCalDigitToMC5

private

Index of fifth MCParticle related to CalDigit, PureCsI option.

Definition at line 339 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMC5PDG

std::vector<int>* m_eclPureCalDigitToMC5PDG

private

PDG code of fifth MCParticle related to CalDigit, PureCsI option.

Definition at line 340 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMCWeight1

std::vector<double>* m_eclPureCalDigitToMCWeight1

private

Energy contribution of first MCParticle related to CalDigit, PureCsI option.

Definition at line 326 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMCWeight2

std::vector<double>* m_eclPureCalDigitToMCWeight2

private

Energy contribution of second MCParticle related to CalDigit, PureCsI option.

Definition at line 330 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMCWeight3

std::vector<double>* m_eclPureCalDigitToMCWeight3

private

Energy contribution of third MCParticle related to CalDigit, PureCsI option.

Definition at line 334 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMCWeight4

std::vector<double>* m_eclPureCalDigitToMCWeight4

private

Energy contribution of fourth MCParticle related to CalDigit, PureCsI option.

Definition at line 338 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToMCWeight5

std::vector<double>* m_eclPureCalDigitToMCWeight5

private

Energy contribution of fifth MCParticle related to CalDigit, PureCsI option.

Definition at line 342 of file ECLDataAnalysisModule.h.

◆ m_eclPureCalDigitToShower

std::vector<int>* m_eclPureCalDigitToShower

private

Index of ECLShower related to that CalDigit, PureCsI option.

Definition at line 346 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterAbsZernike40

std::vector<double>* m_eclPureClusterAbsZernike40

private

Reconstructed Zernike40, PureCsI option.

Definition at line 411 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterAbsZernike51

std::vector<double>* m_eclPureClusterAbsZernike51

private

Reconstructed Zernike51, PureCsI option.

Definition at line 412 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterCellId

std::vector<int>* m_eclPureClusterCellId

private

CellId with highest energy deposit in Cluster, PureCsI option.

Definition at line 404 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterClosestTrackDist

std::vector<double>* m_eclPureClusterClosestTrackDist

private

Reconstructed Distance to Closest Track, PureCsI option.

Definition at line 410 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterCrystalHealth

std::vector<int>* m_eclPureClusterCrystalHealth

private

Crystal health flag, PureCsI option.

Definition at line 407 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterDeltaL

std::vector<double>* m_eclPureClusterDeltaL

private

Reconstructed Cluster DeltaL, PureCsI option.

Definition at line 409 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterDeltaTime99

std::vector<double>* m_eclPureClusterDeltaTime99

private

Reconstructed DeltaT99, PureCsI option.

Definition at line 416 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterDetectorRegion

std::vector<int>* m_eclPureClusterDetectorRegion

private

Clusters detector region, PureCsI option.

Definition at line 417 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterE1oE9

std::vector<double>* m_eclPureClusterE1oE9

private

Reconstructed E1oE9, PureCsI option.

Definition at line 415 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterE9oE21

std::vector<double>* m_eclPureClusterE9oE21

private

Ratio of 3x3 over 5x5 crystal matrices energies for Cluster, PureCsI option.

Definition at line 402 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterEnergy

std::vector<double>* m_eclPureClusterEnergy

private

Cluster energy, PureCsI option.

Definition at line 392 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterEnergyDepSum

std::vector<double>* m_eclPureClusterEnergyDepSum

private

Cluster simulated energy, PureCsI option.

Definition at line 399 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterEnergyError

std::vector<double>* m_eclPureClusterEnergyError

private

Cluster energy error, PureCsI option.

Definition at line 393 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterHasNeutralHadronHypothesis

std::vector<int>* m_eclPureClusterHasNeutralHadronHypothesis

private

Cluster has neutral hadron hypothesis, PureCsI option.

Definition at line 419 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterHasNPhotonHypothesis

std::vector<int>* m_eclPureClusterHasNPhotonHypothesis

private

Cluster has n-photon hypothesis, PureCsI option.

Definition at line 418 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterHighestE

std::vector<double>* m_eclPureClusterHighestE

private

Highest energy deposit (per crystal) in Cluster, PureCsI option.

Definition at line 403 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterIdx

std::vector<int>* m_eclPureClusterIdx

private

ECLCluster index, PureCsI option.

Definition at line 369 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterIsTrack

std::vector<bool>* m_eclPureClusterIsTrack

private

Flag for charged clusters, PureCsI option.

Definition at line 408 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterLat

std::vector<double>* m_eclPureClusterLat

private

Cluster shape parameter LAT, PureCsI option.

Definition at line 405 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterMultip

int m_eclPureClusterMultip

private

Number of ECLClusterss per event, PureCsI option.

Definition at line 368 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterNofCrystals

std::vector<int>* m_eclPureClusterNofCrystals

private

Number of crystals in Cluster, PureCsI option.

Definition at line 406 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterPhi

std::vector<double>* m_eclPureClusterPhi

private

Cluster azimuthal direction, PureCsI option.

Definition at line 396 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterPhiError

std::vector<double>* m_eclPureClusterPhiError

private

Cluster error on azimuthal direction, PureCsI option.

Definition at line 397 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterR

std::vector<double>* m_eclPureClusterR

private

Cluster distance from IP, PureCsI option.

Definition at line 398 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusters

StoreArray<ECLCluster> m_eclPureClusters

private

Store array: ECLPureCluster.

Definition at line 138 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterSecondMoment

std::vector<double>* m_eclPureClusterSecondMoment

private

Reconstructed Cluster Second Moment, PureCsI option.

Definition at line 414 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterTheta

std::vector<double>* m_eclPureClusterTheta

private

Cluster polar direction, PureCsI option.

Definition at line 394 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterThetaError

std::vector<double>* m_eclPureClusterThetaError

private

Cluster error on polar direction, PureCsI option.

Definition at line 395 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterTiming

std::vector<double>* m_eclPureClusterTiming

private

Cluster time, PureCsI option.

Definition at line 400 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterTimingError

std::vector<double>* m_eclPureClusterTimingError

private

Cluster time error, PureCsI option.

Definition at line 401 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToBkgWeight

std::vector<double>* m_eclPureClusterToBkgWeight

private

Remaining energy contribution not associated to first five MCParticles related to ECLCluster, PureCsI option.

Definition at line 391 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMC1

std::vector<int>* m_eclPureClusterToMC1

private

Index of first MCParticle related to ECLCluster, PureCsI option.

Definition at line 370 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMC1PDG

std::vector<int>* m_eclPureClusterToMC1PDG

private

PDG code of first MCParticle related to ECLCluster, PureCsI option.

Definition at line 373 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMC2

std::vector<int>* m_eclPureClusterToMC2

private

Index of second MCParticle related to ECLCluster, PureCsI option.

Definition at line 374 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMC2PDG

std::vector<int>* m_eclPureClusterToMC2PDG

private

PDG code of second MCParticle related to ECLCluster, PureCsI option.

Definition at line 377 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMC3

std::vector<int>* m_eclPureClusterToMC3

private

Index of third MCParticle related to ECLCluster, PureCsI option.

Definition at line 378 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMC3PDG

std::vector<int>* m_eclPureClusterToMC3PDG

private

PDG code of third MCParticle related to ECLCluster, PureCsI option.

Definition at line 381 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMC4

std::vector<int>* m_eclPureClusterToMC4

private

Index of fourth MCParticle related to ECLCluster, PureCsI option.

Definition at line 382 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMC4PDG

std::vector<int>* m_eclPureClusterToMC4PDG

private

PDG code of fourth MCParticle related to ECLCluster, PureCsI option.

Definition at line 385 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMC5

std::vector<int>* m_eclPureClusterToMC5

private

Index of fifth MCParticle related to ECLCluster, PureCsI option.

Definition at line 386 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMC5PDG

std::vector<int>* m_eclPureClusterToMC5PDG

private

PDG code of fifth MCParticle related to ECLCluster, PureCsI option.

Definition at line 389 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMCWeight1

std::vector<double>* m_eclPureClusterToMCWeight1

private

Energy contribution of first MCParticle related to ECLCluster, PureCsI option.

Definition at line 372 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMCWeight2

std::vector<double>* m_eclPureClusterToMCWeight2

private

Energy contribution of second MCParticle related to ECLCluster, PureCsI option.

Definition at line 376 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMCWeight3

std::vector<double>* m_eclPureClusterToMCWeight3

private

Energy contribution of third MCParticle related to ECLCluster, PureCsI option.

Definition at line 380 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMCWeight4

std::vector<double>* m_eclPureClusterToMCWeight4

private

Energy contribution of fourth MCParticle related to ECLCluster, PureCsI option.

Definition at line 384 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterToMCWeight5

std::vector<double>* m_eclPureClusterToMCWeight5

private

Energy contribution of 5th MCParticle related to ECLCluster, PureCsI option.

Definition at line 388 of file ECLDataAnalysisModule.h.

◆ m_eclPureClusterZernikeMVA

std::vector<double>* m_eclPureClusterZernikeMVA

private

Output of MVA classifier based on Zernike Momenta, PureCsI option.

Definition at line 413 of file ECLDataAnalysisModule.h.

◆ m_eclPureConnectedRegions

StoreArray<ECLConnectedRegion> m_eclPureConnectedRegions

private

Store array: ECLPureConnectedRegion.

Definition at line 134 of file ECLDataAnalysisModule.h.

◆ m_eclPureCRIdx

std::vector<int>* m_eclPureCRIdx

private

Connected Region ID, PureCsI option.

Definition at line 354 of file ECLDataAnalysisModule.h.

◆ m_eclPureCRIsTrack

std::vector<int>* m_eclPureCRIsTrack

private

Int for Connected Region - Track Match, PureCsI option.

Definition at line 355 of file ECLDataAnalysisModule.h.

◆ m_eclPureCRLikelihoodChargedHadron

std::vector<double>* m_eclPureCRLikelihoodChargedHadron

private

Connected Region Charged Hadron Likelihood, PureCsI option.

Definition at line 357 of file ECLDataAnalysisModule.h.

◆ m_eclPureCRLikelihoodElectronNGamma

std::vector<double>* m_eclPureCRLikelihoodElectronNGamma

private

Connected Region Electron Likelihood, PureCsI option.

Definition at line 358 of file ECLDataAnalysisModule.h.

◆ m_eclPureCRLikelihoodMergedPi0

std::vector<double>* m_eclPureCRLikelihoodMergedPi0

private

Connected Region Merged Pi0 Likelihood, PureCsI option.

Definition at line 361 of file ECLDataAnalysisModule.h.

◆ m_eclPureCRLikelihoodMIPNGamma

std::vector<double>* m_eclPureCRLikelihoodMIPNGamma

private

Connected Region MIP Likelihood, PureCsI option.

Definition at line 356 of file ECLDataAnalysisModule.h.

◆ m_eclPureCRLikelihoodNeutralHadron

std::vector<double>* m_eclPureCRLikelihoodNeutralHadron

private

Connected Region Neutral Hadron Likelihood, PureCsI option.

Definition at line 360 of file ECLDataAnalysisModule.h.

◆ m_eclPureCRLikelihoodNGamma

std::vector<double>* m_eclPureCRLikelihoodNGamma

private

Connected Region Gamma Likelihood, PureCsI option.

Definition at line 359 of file ECLDataAnalysisModule.h.

◆ m_eclPureDigitAmp

std::vector<int>* m_eclPureDigitAmp

private

ECLDigit amplitude, PureCsI option.

Definition at line 316 of file ECLDataAnalysisModule.h.

◆ m_eclPureDigitCellId

std::vector<int>* m_eclPureDigitCellId

private

Number of ECLDigit CellId, PureCsI option.

Definition at line 315 of file ECLDataAnalysisModule.h.

◆ m_eclPureDigitFitQuality

std::vector<int>* m_eclPureDigitFitQuality

private

ECLDigit fit quality, PureCsI option.

Definition at line 318 of file ECLDataAnalysisModule.h.

◆ m_eclPureDigitIdx

std::vector<int>* m_eclPureDigitIdx

private

ECLDigit index, PureCsI option.

Definition at line 313 of file ECLDataAnalysisModule.h.

◆ m_eclPureDigitMultip

int m_eclPureDigitMultip

private

Number of ECLDigits per event, PureCsI option.

Definition at line 312 of file ECLDataAnalysisModule.h.

◆ m_eclPureDigits

StoreArray<ECLDigit> m_eclPureDigits

private

Store array: ECLPureDigit.

Definition at line 130 of file ECLDataAnalysisModule.h.

◆ m_eclPureDigitTimeFit

std::vector<int>* m_eclPureDigitTimeFit

private

ECLDigit timing, PureCsI option.

Definition at line 317 of file ECLDataAnalysisModule.h.

◆ m_eclPureDigitToCluster

std::vector<int>* m_eclPureDigitToCluster

private

ECLDigit To Cluster, PureCsI option.

Definition at line 319 of file ECLDataAnalysisModule.h.

◆ m_eclPureDigitToMC

std::vector<int>* m_eclPureDigitToMC

private

Index of MCParticle related to that ECLDigit, PureCsI option.

Definition at line 314 of file ECLDataAnalysisModule.h.

◆ m_eclPureLMCellId

std::vector<int>* m_eclPureLMCellId

private

Local Maximum Cell ID, PureCsI option.

Definition at line 366 of file ECLDataAnalysisModule.h.

◆ m_eclPureLMId

std::vector<int>* m_eclPureLMId

private

Local Maximum ID, PureCsI option.

Definition at line 364 of file ECLDataAnalysisModule.h.

◆ m_eclPureLMMultip

int m_eclPureLMMultip

private

Local Maxima multiplicity, PureCsI option.

Definition at line 363 of file ECLDataAnalysisModule.h.

◆ m_eclPureLMType

std::vector<int>* m_eclPureLMType

private

Local Maximum type, PureCsI option.

Definition at line 365 of file ECLDataAnalysisModule.h.

◆ m_eclPureLocalMaximums

StoreArray<ECLLocalMaximum> m_eclPureLocalMaximums

private

Store array: ECLPureLocalMaximum.

Definition at line 140 of file ECLDataAnalysisModule.h.

◆ m_eclPureShowers

StoreArray<ECLShower> m_eclPureShowers

private

Store array: ECLPureShower.

Definition at line 136 of file ECLDataAnalysisModule.h.

◆ m_eclShowerAbsZernike40

std::vector<double>* m_eclShowerAbsZernike40

private

Shower Zernike40 Moment.

Definition at line 487 of file ECLDataAnalysisModule.h.

◆ m_eclShowerAbsZernike51

std::vector<double>* m_eclShowerAbsZernike51

private

Shower Zernike51 Moment.

Definition at line 488 of file ECLDataAnalysisModule.h.

◆ m_eclShowerCentralCellId

std::vector<int>* m_eclShowerCentralCellId

private

Cell ID for most energetic crystal.

Definition at line 477 of file ECLDataAnalysisModule.h.

◆ m_eclShowerConnectedRegionId

std::vector<int>* m_eclShowerConnectedRegionId

private

Matched Connetcted Region Idx.

Definition at line 475 of file ECLDataAnalysisModule.h.

◆ m_eclShowerE1oE9

std::vector<double>* m_eclShowerE1oE9

private

Shower E1/E9.

Definition at line 491 of file ECLDataAnalysisModule.h.

◆ m_eclShowerE9oE21

std::vector<double>* m_eclShowerE9oE21

private

Shower E9oE21.

Definition at line 472 of file ECLDataAnalysisModule.h.

◆ m_eclShowerEnergy

std::vector<double>* m_eclShowerEnergy

private

Shower Energy.

Definition at line 467 of file ECLDataAnalysisModule.h.

◆ m_eclShowerEnergyError

std::vector<double>* m_eclShowerEnergyError

private

Shower Energy Error.

Definition at line 478 of file ECLDataAnalysisModule.h.

◆ m_eclShowerHighestE1mE2

std::vector<double>* m_eclShowerHighestE1mE2

private

Energy difference for 2 highest energy deposits in shower.

Definition at line 498 of file ECLDataAnalysisModule.h.

◆ m_eclShowerHighestEnergy

std::vector<double>* m_eclShowerHighestEnergy

private

Shower Highest Energy Crystal Energy.

Definition at line 482 of file ECLDataAnalysisModule.h.

◆ m_eclShowerHypothesisId

std::vector<int>* m_eclShowerHypothesisId

private

Shower Particle Hypothesis ID.

Definition at line 476 of file ECLDataAnalysisModule.h.

◆ m_eclShowerIdx

std::vector<int>* m_eclShowerIdx

private

Shower Index.

Definition at line 422 of file ECLDataAnalysisModule.h.

◆ m_eclShowerIsCluster

std::vector<bool>* m_eclShowerIsCluster

private

Shower Cluster Match.

Definition at line 493 of file ECLDataAnalysisModule.h.

◆ m_eclShowerIsTrack

std::vector<int>* m_eclShowerIsTrack

private

Shower Track Match.

Definition at line 492 of file ECLDataAnalysisModule.h.

◆ m_eclShowerLateralEnergy

std::vector<double>* m_eclShowerLateralEnergy

private

Shower Lateral Energy.

Definition at line 483 of file ECLDataAnalysisModule.h.

◆ m_eclShowerMCFFlightMatch

std::vector<int>* m_eclShowerMCFFlightMatch

private

Int, 1 if primary particle flight direction is "well" reconstructed in ECL, 0 otherwise, DEBUG PURPOSE.

Definition at line 497 of file ECLDataAnalysisModule.h.

◆ m_eclShowerMCFlightMatch

std::vector<int>* m_eclShowerMCFlightMatch

private

Int, 1 if particle flight direction is "well" reconstructed in ECL, 0 otherwise.

Definition at line 495 of file ECLDataAnalysisModule.h.

◆ m_eclShowerMCVtxInEcl

std::vector<int>* m_eclShowerMCVtxInEcl

private

Int, 1 if particle decays (interacts) in ECL, 0 otherwise.

Definition at line 494 of file ECLDataAnalysisModule.h.

◆ m_eclShowerMinTrkDistance

std::vector<double>* m_eclShowerMinTrkDistance

private

Shower Min Dist to Track.

Definition at line 484 of file ECLDataAnalysisModule.h.

◆ m_eclShowerMultip

int m_eclShowerMultip

private

Number of ECLShowers per event.

Definition at line 421 of file ECLDataAnalysisModule.h.

◆ m_eclShowerNHits

std::vector<double>* m_eclShowerNHits

private

Shower NHits.

Definition at line 471 of file ECLDataAnalysisModule.h.

◆ m_eclShowerNumberOfCrystalsForEnergy

std::vector<double>* m_eclShowerNumberOfCrystalsForEnergy

private

Number of crystals used for energy calculation.

Definition at line 499 of file ECLDataAnalysisModule.h.

◆ m_eclShowerPhi

std::vector<double>* m_eclShowerPhi

private

Shower Phi.

Definition at line 469 of file ECLDataAnalysisModule.h.

◆ m_eclShowerPhiError

std::vector<double>* m_eclShowerPhiError

private

Shower Phi Error.

Definition at line 480 of file ECLDataAnalysisModule.h.

◆ m_eclShowerR

std::vector<double>* m_eclShowerR

private

Shower R.

Definition at line 470 of file ECLDataAnalysisModule.h.

◆ m_eclShowers

StoreArray<ECLShower> m_eclShowers

private

Store array: ECLShower.

Definition at line 97 of file ECLDataAnalysisModule.h.

◆ m_eclShowerSecondMoment

std::vector<double>* m_eclShowerSecondMoment

private

Shower Second Moment.

Definition at line 490 of file ECLDataAnalysisModule.h.

◆ m_eclShowerShowerDepth

std::vector<double>* m_eclShowerShowerDepth

private

Shower Depth.

Definition at line 486 of file ECLDataAnalysisModule.h.

◆ m_eclShowerSimHitSum

std::vector<double>* m_eclShowerSimHitSum

private

Full energy contribution related to ECLShower.

Definition at line 465 of file ECLDataAnalysisModule.h.

◆ m_eclShowerT99

std::vector<double>* m_eclShowerT99

private

Shower T99.

Definition at line 474 of file ECLDataAnalysisModule.h.

◆ m_eclShowerTheta

std::vector<double>* m_eclShowerTheta

private

Shower Theta.

Definition at line 468 of file ECLDataAnalysisModule.h.

◆ m_eclShowerThetaError

std::vector<double>* m_eclShowerThetaError

private

Shower Theta Error.

Definition at line 479 of file ECLDataAnalysisModule.h.

◆ m_eclShowerTime

std::vector<double>* m_eclShowerTime

private

Shower Timing.

Definition at line 473 of file ECLDataAnalysisModule.h.

◆ m_eclShowerTimeResolution

std::vector<double>* m_eclShowerTimeResolution

private

Shower Time Resolution.

Definition at line 481 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToBkgWeight

std::vector<double>* m_eclShowerToBkgWeight

private

Remaining energy contribution not associated to first five MCParticles related to ECLShower.

Definition at line 459 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToLM1

std::vector<int>* m_eclShowerToLM1

private

Index of first maximum related to ECLShower.

Definition at line 460 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToLM2

std::vector<int>* m_eclShowerToLM2

private

Index of 2nd maximum related to ECLShower.

Definition at line 461 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToLM3

std::vector<int>* m_eclShowerToLM3

private

Index of 3rd maximum related to ECLShower.

Definition at line 462 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToLM4

std::vector<int>* m_eclShowerToLM4

private

Index of 4th maximum related to ECLShower.

Definition at line 463 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToLM5

std::vector<int>* m_eclShowerToLM5

private

Index of 5th maximum related to ECLShower.

Definition at line 464 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC1

std::vector<int>* m_eclShowerToMC1

private

Index of first MCParticle related to ECLShower.

Definition at line 423 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC1GMoth

std::vector<int>* m_eclShowerToMC1GMoth

private

GMother index of first MCParticle related to ECLShower.

Definition at line 428 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC1GMothPDG

std::vector<int>* m_eclShowerToMC1GMothPDG

private

PDG code of Gparent of first MCParticle related to ECLShower.

Definition at line 429 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC1Moth

std::vector<int>* m_eclShowerToMC1Moth

private

Mother index of first MCParticle related to ECLShower.

Definition at line 426 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC1MothPDG

std::vector<int>* m_eclShowerToMC1MothPDG

private

PDG code of parent of first MCParticle related to ECLShower.

Definition at line 427 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC1PDG

std::vector<int>* m_eclShowerToMC1PDG

private

PDG code of first MCParticle related to ECLShower.

Definition at line 425 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC2

std::vector<int>* m_eclShowerToMC2

private

Index of second MCParticle related to ECLShower.

Definition at line 430 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC2GMoth

std::vector<int>* m_eclShowerToMC2GMoth

private

GMother index of second MCParticle related to ECLShower.

Definition at line 435 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC2GMothPDG

std::vector<int>* m_eclShowerToMC2GMothPDG

private

PDG code of Gparent of second MCParticle related to ECLShower.

Definition at line 436 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC2Moth

std::vector<int>* m_eclShowerToMC2Moth

private

Mother index of first MCParticle related to ECLShower.

Definition at line 433 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC2MothPDG

std::vector<int>* m_eclShowerToMC2MothPDG

private

PDG code of parent of first MCParticle related to ECLShower.

Definition at line 434 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC2PDG

std::vector<int>* m_eclShowerToMC2PDG

private

PDG code of second MCParticle related to ECLShower.

Definition at line 432 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC3

std::vector<int>* m_eclShowerToMC3

private

Index of third MCParticle related to ECLShower.

Definition at line 437 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC3GMoth

std::vector<int>* m_eclShowerToMC3GMoth

private

GMother index of third MCParticle related to ECLShower.

Definition at line 442 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC3GMothPDG

std::vector<int>* m_eclShowerToMC3GMothPDG

private

PDG code of Gparent of third MCParticle related to ECLShower.

Definition at line 443 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC3Moth

std::vector<int>* m_eclShowerToMC3Moth

private

Mother index of first MCParticle related to ECLShower.

Definition at line 440 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC3MothPDG

std::vector<int>* m_eclShowerToMC3MothPDG

private

PDG code of parent of first MCParticle related to ECLShower.

Definition at line 441 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC3PDG

std::vector<int>* m_eclShowerToMC3PDG

private

PDG code of third MCParticle related to ECLShower.

Definition at line 439 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC4

std::vector<int>* m_eclShowerToMC4

private

Index of fourth MCParticle related to ECLShower.

Definition at line 444 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC4GMoth

std::vector<int>* m_eclShowerToMC4GMoth

private

GMother index of fourth MCParticle related to ECLShower.

Definition at line 449 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC4GMothPDG

std::vector<int>* m_eclShowerToMC4GMothPDG

private

PDG code of Gparent of fourth MCParticle related to ECLShower.

Definition at line 450 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC4Moth

std::vector<int>* m_eclShowerToMC4Moth

private

Mother index of fourth MCParticle related to ECLShower.

Definition at line 447 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC4MothPDG

std::vector<int>* m_eclShowerToMC4MothPDG

private

PDG code of parent of fourth MCParticle related to ECLShower.

Definition at line 448 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC4PDG

std::vector<int>* m_eclShowerToMC4PDG

private

PDG code of fourth MCParticle related to ECLShower.

Definition at line 446 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC5

std::vector<int>* m_eclShowerToMC5

private

Index of fifth MCParticle related to ECLShower.

Definition at line 451 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC5GMoth

std::vector<int>* m_eclShowerToMC5GMoth

private

GMother index of fifth MCParticle related to ECLShower.

Definition at line 456 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC5GMothPDG

std::vector<int>* m_eclShowerToMC5GMothPDG

private

PDG code of Gparent of fifth MCParticle related to ECLShower.

Definition at line 457 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC5Moth

std::vector<int>* m_eclShowerToMC5Moth

private

Mother index of fifth MCParticle related to ECLShower.

Definition at line 454 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC5MothPDG

std::vector<int>* m_eclShowerToMC5MothPDG

private

PDG code of parent of fifth MCParticle related to ECLShower.

Definition at line 455 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMC5PDG

std::vector<int>* m_eclShowerToMC5PDG

private

PDG code of fifth MCParticle related to ECLShower.

Definition at line 453 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMCWeight1

std::vector<double>* m_eclShowerToMCWeight1

private

Energy contribution of first MCParticle related to ECLShower.

Definition at line 424 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMCWeight2

std::vector<double>* m_eclShowerToMCWeight2

private

Energy contribution of second MCParticle related to ECLShower.

Definition at line 431 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMCWeight3

std::vector<double>* m_eclShowerToMCWeight3

private

Energy contribution of third MCParticle related to ECLShower.

Definition at line 438 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMCWeight4

std::vector<double>* m_eclShowerToMCWeight4

private

Energy contribution of fourth MCParticle related to ECLShower.

Definition at line 445 of file ECLDataAnalysisModule.h.

◆ m_eclShowerToMCWeight5

std::vector<double>* m_eclShowerToMCWeight5

private

Energy contribution of fifth MCParticle related to ECLShower.

Definition at line 452 of file ECLDataAnalysisModule.h.

◆ m_eclShowerTrkDepth

std::vector<double>* m_eclShowerTrkDepth

private

Shower Track Depth.

Definition at line 485 of file ECLDataAnalysisModule.h.

◆ m_eclShowerUncEnergy

std::vector<double>* m_eclShowerUncEnergy

private

Shower bare energy.

Definition at line 466 of file ECLDataAnalysisModule.h.

◆ m_eclShowerZernikeMVA

std::vector<double>* m_eclShowerZernikeMVA

private

Shower ZernikeMVA.

Definition at line 489 of file ECLDataAnalysisModule.h.

◆ m_eclSimHitCellId

std::vector<int>* m_eclSimHitCellId

private

ECLSimHit CellId.

Definition at line 238 of file ECLDataAnalysisModule.h.

◆ m_eclSimHitEnergyDep

std::vector<double>* m_eclSimHitEnergyDep

private

Energy deposition of ECLSimHit.

Definition at line 240 of file ECLDataAnalysisModule.h.

◆ m_eclSimHitFlightTime

std::vector<double>* m_eclSimHitFlightTime

private

ECLSimhit Flight Time.

Definition at line 241 of file ECLDataAnalysisModule.h.

◆ m_eclSimHitIdx

std::vector<int>* m_eclSimHitIdx

private

Index of ECLSimHit.

Definition at line 236 of file ECLDataAnalysisModule.h.

◆ m_eclSimHitMultip

int m_eclSimHitMultip

private

Number of ECLSimHits per event.

Definition at line 235 of file ECLDataAnalysisModule.h.

◆ m_eclSimHitPdg

std::vector<int>* m_eclSimHitPdg

private

PDG code of MCParticle associated to that ECLDigit.

Definition at line 239 of file ECLDataAnalysisModule.h.

◆ m_eclSimHitPx

std::vector<double>* m_eclSimHitPx

private

ECLSimHit PX.

Definition at line 245 of file ECLDataAnalysisModule.h.

◆ m_eclSimHitPy

std::vector<double>* m_eclSimHitPy

private

ECLSimHit PY.

Definition at line 246 of file ECLDataAnalysisModule.h.

◆ m_eclSimHitPz

std::vector<double>* m_eclSimHitPz

private

ECLSimHit PZ.

Definition at line 247 of file ECLDataAnalysisModule.h.

◆ m_eclSimHits

StoreArray<ECLSimHit> m_eclSimHits

private

Store array: ECLSimHit.

Definition at line 92 of file ECLDataAnalysisModule.h.

◆ m_eclSimHitToMC

std::vector<int>* m_eclSimHitToMC

private

Index of MCParticle related to that ECLSimHit.

Definition at line 237 of file ECLDataAnalysisModule.h.

◆ m_eclSimHitX

std::vector<double>* m_eclSimHitX

private

ECLSimHit X position.

Definition at line 242 of file ECLDataAnalysisModule.h.

◆ m_eclSimHitY

std::vector<double>* m_eclSimHitY

private

ECLSimHit Y position.

Definition at line 243 of file ECLDataAnalysisModule.h.

◆ m_eclSimHitZ

std::vector<double>* m_eclSimHitZ

private

ECLSimHit Z position.

Definition at line 244 of file ECLDataAnalysisModule.h.

◆ m_eventLevelClusteringInfo

StoreObjPtr<EventLevelClusteringInfo> m_eventLevelClusteringInfo

private

Store object pointer: EventLevelClusteringInfo.

Definition at line 100 of file ECLDataAnalysisModule.h.

◆ m_eventmetadata

StoreObjPtr<EventMetaData> m_eventmetadata

private

Store object pointer: EventMetaData.

Definition at line 101 of file ECLDataAnalysisModule.h.

◆ m_hasReturnValue

bool m_hasReturnValue

privateinherited

True, if the return value is set.

Definition at line 517 of file Module.h.

◆ m_iEvent

int m_iEvent

private

Event number.

Definition at line 168 of file ECLDataAnalysisModule.h.

◆ m_iExperiment

int m_iExperiment

private

Experiment number.

Definition at line 166 of file ECLDataAnalysisModule.h.

◆ m_iRun

int m_iRun

private

Run number.

Definition at line 167 of file ECLDataAnalysisModule.h.

◆ m_logConfig

LogConfig m_logConfig

privateinherited

The log system configuration of the module.

Definition at line 513 of file Module.h.

◆ m_mcDecayVtxX

std::vector<double>* m_mcDecayVtxX

private

MCParticle decay vertex X.

Definition at line 511 of file ECLDataAnalysisModule.h.

◆ m_mcDecayVtxY

std::vector<double>* m_mcDecayVtxY

private

MCParticle decay vertex Y.

Definition at line 512 of file ECLDataAnalysisModule.h.

◆ m_mcDecayVtxZ

std::vector<double>* m_mcDecayVtxZ

private

MCParticle decay vertex Z.

Definition at line 513 of file ECLDataAnalysisModule.h.

◆ m_mcEnergy

std::vector<double>* m_mcEnergy

private

MCParticle energyx.

Definition at line 507 of file ECLDataAnalysisModule.h.

◆ m_mcGGMothPdg

std::vector<int>* m_mcGGMothPdg

private

MCParticle greand-grandmother particle PDG code.

Definition at line 506 of file ECLDataAnalysisModule.h.

◆ m_mcGMothPdg

std::vector<int>* m_mcGMothPdg

private

MCParticle grandmother particle PDG code.

Definition at line 505 of file ECLDataAnalysisModule.h.

◆ m_mcIdx

std::vector<int>* m_mcIdx

private

MCParticle index.

Definition at line 502 of file ECLDataAnalysisModule.h.

◆ m_mcMothPdg

std::vector<int>* m_mcMothPdg

private

MCParticle mother particle PDG code.

Definition at line 504 of file ECLDataAnalysisModule.h.

◆ m_mcMultip

int m_mcMultip

private

Multiplicity of MCParticles.

Definition at line 501 of file ECLDataAnalysisModule.h.

◆ m_mcParticles

StoreArray<MCParticle> m_mcParticles

private

MCParticles StoreArray.

Definition at line 161 of file ECLDataAnalysisModule.h.

◆ m_mcPdg

std::vector<int>* m_mcPdg

private

MCParticle PDG code.

Definition at line 503 of file ECLDataAnalysisModule.h.

◆ m_mcProdVtxX

std::vector<double>* m_mcProdVtxX

private

MCParticle production vertex X.

Definition at line 514 of file ECLDataAnalysisModule.h.

◆ m_mcProdVtxY

std::vector<double>* m_mcProdVtxY

private

MCParticle production vertex Y.

Definition at line 515 of file ECLDataAnalysisModule.h.

◆ m_mcProdVtxZ

std::vector<double>* m_mcProdVtxZ

private

MCParticle production vertex Z.

Definition at line 516 of file ECLDataAnalysisModule.h.

◆ m_mcPx

std::vector<double>* m_mcPx

private

MCParticle momentum X direction.

Definition at line 508 of file ECLDataAnalysisModule.h.

◆ m_mcPy

std::vector<double>* m_mcPy

private

MCParticle momentum Y direction.

Definition at line 509 of file ECLDataAnalysisModule.h.

◆ m_mcPz

std::vector<double>* m_mcPz

private

MCParticle momentum Z direction.

Definition at line 510 of file ECLDataAnalysisModule.h.

◆ m_mcSecondaryPhysProc

std::vector<int>* m_mcSecondaryPhysProc

private

Flag for secondary physics process.

Definition at line 517 of file ECLDataAnalysisModule.h.

◆ m_moduleParamList

ModuleParamList m_moduleParamList

privateinherited

List storing and managing all parameter of the module.

Definition at line 515 of file Module.h.

◆ m_name

std::string m_name

privateinherited

The name of the module, saved as a string (user-modifiable)

Definition at line 507 of file Module.h.

◆ m_nECLCalDigitsOutOfTimeBarrel

uint16_t m_nECLCalDigitsOutOfTimeBarrel {0}

private

Number of out of time, energetic ECLCalDigits, Barrel.

Definition at line 174 of file ECLDataAnalysisModule.h.

174{0};

◆ m_nECLCalDigitsOutOfTimeBWD

uint16_t m_nECLCalDigitsOutOfTimeBWD {0}

private

Number of out of time, energetic ECLCalDigits, BWD.

Definition at line 176 of file ECLDataAnalysisModule.h.

176{0};

◆ m_nECLCalDigitsOutOfTimeFWD

uint16_t m_nECLCalDigitsOutOfTimeFWD {0}

private

Number of out of time, energetic ECLCalDigits, FWD.

Definition at line 172 of file ECLDataAnalysisModule.h.

172{0};

◆ m_nECLShowersRejectedBarrel

uint8_t m_nECLShowersRejectedBarrel {0}

private

Number of photon showers that are rejected before storing to mdst (max.

255), Barrel.

Definition at line 180 of file ECLDataAnalysisModule.h.

180{0};

◆ m_nECLShowersRejectedBWD

uint8_t m_nECLShowersRejectedBWD {0}

private

Number of photon showers that are rejected before storing to mdst (max.

255), BWD.

Definition at line 182 of file ECLDataAnalysisModule.h.

182{0};

◆ m_nECLShowersRejectedFWD

uint8_t m_nECLShowersRejectedFWD {0}

private

Number of photon showers that are rejected before storing to mdst (max.

255), FWD.

Definition at line 178 of file ECLDataAnalysisModule.h.

178{0};

◆ m_package

std::string m_package

privateinherited

Package this module is found in (may be empty).

Definition at line 509 of file Module.h.

◆ m_propertyFlags

unsigned int m_propertyFlags

privateinherited

The properties of the module as bitwise or (with |) of EModulePropFlags.

Definition at line 511 of file Module.h.

◆ m_returnValue

int m_returnValue

privateinherited

The return value.

Definition at line 518 of file Module.h.

◆ m_rootFileName

std::string m_rootFileName

private

name of the root file

Definition at line 81 of file ECLDataAnalysisModule.h.

◆ m_rootFilePtr

TFile* m_rootFilePtr

private

members of ECLReconstructor Module

root file used for storing info

Definition at line 80 of file ECLDataAnalysisModule.h.

◆ m_trackFitResults

StoreArray<TrackFitResult> m_trackFitResults

private

TrackFitResult storeArray.

Definition at line 90 of file ECLDataAnalysisModule.h.

◆ m_tracks

StoreArray<Track> m_tracks

private

Tracks storeArray.

Definition at line 89 of file ECLDataAnalysisModule.h.

◆ m_tree

TTree* m_tree

private

Root tree and file for saving the output.

Definition at line 163 of file ECLDataAnalysisModule.h.

◆ m_trkCharge

std::vector<int>* m_trkCharge

private

Track charge.

Definition at line 522 of file ECLDataAnalysisModule.h.

◆ m_trkIdx

std::vector<int>* m_trkIdx

private

Track index.

Definition at line 521 of file ECLDataAnalysisModule.h.

◆ m_trkMultip

int m_trkMultip

private

Track Multiplicity.

Definition at line 519 of file ECLDataAnalysisModule.h.

◆ m_trkP

std::vector<double>* m_trkP

private

Track momentum.

Definition at line 526 of file ECLDataAnalysisModule.h.

◆ m_trkPdg

std::vector<int>* m_trkPdg

private

Track PDG code.

Definition at line 520 of file ECLDataAnalysisModule.h.

◆ m_trkPhi

std::vector<double>* m_trkPhi

private

Track azimuthal direction.

Definition at line 528 of file ECLDataAnalysisModule.h.

◆ m_trkPx

std::vector<double>* m_trkPx

private

Track momentum along X direction.

Definition at line 523 of file ECLDataAnalysisModule.h.

◆ m_trkPy

std::vector<double>* m_trkPy

private

Track momentum along Y direction.

Definition at line 524 of file ECLDataAnalysisModule.h.

◆ m_trkPz

std::vector<double>* m_trkPz

private

Track momentum along Z direction.

Definition at line 525 of file ECLDataAnalysisModule.h.

◆ m_trkTheta

std::vector<double>* m_trkTheta

private

Track polar direction.

Definition at line 527 of file ECLDataAnalysisModule.h.

◆ m_trkX

std::vector<double>* m_trkX

private

Track DOCA X (?)

Definition at line 529 of file ECLDataAnalysisModule.h.

◆ m_trkY

std::vector<double>* m_trkY

private

Track DOCA Y (?)

Definition at line 530 of file ECLDataAnalysisModule.h.

◆ m_trkZ

std::vector<double>* m_trkZ

private

Track DOCA Z (?)

Definition at line 531 of file ECLDataAnalysisModule.h.

◆ m_type

std::string m_type

privateinherited

The type of the module, saved as a string.

Definition at line 508 of file Module.h.

◆ m_writeToRoot

bool m_writeToRoot

private

if true, a rootFile named by m_rootFileName will be filled with info

Definition at line 82 of file ECLDataAnalysisModule.h.

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

ecl/modules/eclDataAnalysis/include/ECLDataAnalysisModule.h
ecl/modules/eclDataAnalysis/src/ECLDataAnalysisModule.cc

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ EAfterConditionPath

Member Enumeration Documentation

◆ EModulePropFlags

Constructor & Destructor Documentation

◆ ECLDataAnalysisModule()

◆ ~ECLDataAnalysisModule()

Member Function Documentation

◆ beginRun()

◆ clone()

◆ def_beginRun()

◆ def_endRun()

◆ def_event()

◆ def_initialize()

◆ def_terminate()

◆ eclCalDigitArrayName()

◆ eclClusterArrayName()

◆ eclConnectedRegionArrayName()

◆ eclDigitArrayName()

◆ eclHitArrayName()

◆ eclLocalMaximumArrayName()

◆ eclPureCalDigitArrayName()

◆ eclPureClusterArrayName()

◆ eclPureConnectedRegionArrayName()

◆ eclPureDigitArrayName()

◆ eclPureLocalMaximumArrayName()

◆ eclPureShowerArrayName()

◆ eclShowerArrayName()

◆ eclSimHitArrayName()

◆ endRun()

◆ evalCondition()

◆ event()

◆ exposePythonAPI()

◆ getAfterConditionPath()

◆ getAllConditionPaths()

◆ getAllConditions()

◆ getCondition()

◆ getConditionPath()

◆ getDescription()

◆ getFileNames()

◆ getLogConfig()

◆ getModules()

◆ getName()

◆ getPackage()

◆ getParamInfoListPython()

◆ getParamList()

◆ getPathString()

◆ getReturnValue()

◆ getType()

◆ hasCondition()

◆ hasProperties()

◆ hasReturnValue()

◆ hasUnsetForcedParams()

◆ if_false()

◆ if_true()

◆ if_value()

◆ initialize()

◆ setAbortLevel()

◆ setDebugLevel()

◆ setDescription()

◆ setLogConfig()

◆ setLogInfo()

◆ setLogLevel()

◆ setName()

◆ setParamList()

◆ setParamPython()

◆ setParamPythonDict()

◆ setPropertyFlags()

◆ setReturnValue() [1/2]

◆ setReturnValue() [2/2]

◆ setType()

◆ terminate()

Member Data Documentation