TileTBAANtuple Class Reference

class to produce TileCal commissioning ntuples More...

#include <TileTBAANtuple.h>

Inheritance diagram for TileTBAANtuple:

Public Member Functions
	TileTBAANtuple (const std::string &name, ISvcLocator *pSvcLocator)
virtual	~TileTBAANtuple ()=default
StatusCode	ntuple_initialize (const EventContext &ctx)
	Alg standard interface function LF TODO: We could have a problem with the new feature introduced by Sasha that quit empty fragments and therefore it can increase or decrease the number of RawChannels fragments dynamically.
StatusCode	ntuple_clear ()
virtual StatusCode	initialize () override
virtual StatusCode	execute (const EventContext &ctx) override
	Execute method.
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
bool	filterPassed (const EventContext &ctx) const
void	setFilterPassed (bool state, const EventContext &ctx) const
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const
const EventContext &	getContext () const
	Deprecated methods (use the ones with EventContext).
bool	filterPassed () const
void	setFilterPassed (bool state) const

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

Private Types
enum	{ NOT_SETUP = -9999 }
enum	{ MAX_MINIDRAWER = 4 , MAX_CHAN = 48 , MAX_DMU = 16 }
typedef std::map< unsignedint, unsignedint, std::less< unsignedint > >::iterator	drawerMap_iterator
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
StatusCode	storeRawChannels (const EventContext &ctx, const SG::ReadHandleKey< TileRawChannelContainer > &containerKey, bool calibMode, std::vector< std::array< float, MAX_CHAN > > *eneVec, std::vector< std::array< float, MAX_CHAN > > *timeVec, std::vector< std::array< float, MAX_CHAN > > *chi2Vec, std::vector< std::array< float, MAX_CHAN > > *pedVec, bool saveDQstatus=false)
	/ Fill ntuple with data from TRC.
StatusCode	storeDigits (const EventContext &ctx, const SG::ReadHandleKey< TileDigitsContainer > &containerKey, const TileDQstatus *dqStatus=nullptr)
	/ Fill Ntuple with info from TileDigits / Return true if the collection is empty, / which means that there are no RawChanels either.
StatusCode	storeDigitsFlx (const EventContext &ctx, const SG::ReadHandleKey< TileDigitsContainer > &containerKey)
StatusCode	storeBeamElements (const EventContext &ctx)
StatusCode	storeCells (const EventContext &ctx)
StatusCode	storeLaser (const EventContext &ctx)
StatusCode	storeHitVector (const EventContext &ctx)
StatusCode	storeHitContainer (const EventContext &ctx)
	/ Fill Ntuple with MC truth info from simulation / Namely, hit energies corrected by photoelectron statistics and Birks' law
void	storeHit (const TileHit *hit, int fragType, int fragId, std::array< float, MAX_CHAN > &ehitVec, std::array< float, MAX_CHAN > &thitVec, const TileSamplingFraction *samplingFraction)
StatusCode	initList (const EventContext &ctx)
StatusCode	initListFlx (const EventContext &ctx)
StatusCode	initNTuple (void)
void	getEta (void)
void	TRIGGER_addBranch (void)
	///////////////////////////////////////////////////////////////////////////
void	MUON_addBranch (void)
	//////////////////////////////////////////////////////////////////////////// /Add MUON variables to the Tree
void	ECAL_addBranch (void)
	//////////////////////////////////////////////////////////////////////////// /Add ECAL variables to the Tree
void	QDC_addBranch (void)
	//////////////////////////////////////////////////////////////////////////// /Add QDC variables to the Tree
void	LASER_addBranch (void)
	//////////////////////////////////////////////////////////////////////////// /Add Tree LASER variables Tree
void	ADDER_addBranch (void)
	//////////////////////////////////////////////////////////////////////////// /Add Tree ADDER variables Tree
void	CISPAR_addBranch (void)
	//////////////////////////////////////////////////////////////////////////// /Add Tree CISPAR variables Tree
void	BEAM_addBranch (void)
	//////////////////////////////////////////////////////////////////////////// /Add Tree BEAM variables Tree
void	DIGI_addBranch (void)
	//////////////////////////////////////////////////////////////////////////// /Add Tree DIGI variables Tree
void	FELIX_addBranch (void)
void	HIT_addBranch (void)
	//////////////////////////////////////////////////////////////////////////// /Add Tree HIT variables Tree
void	ENETOTAL_addBranch (void)
	//////////////////////////////////////////////////////////////////////////// /Add Tree ENETOTAL variables Tree
void	COINCBOARD_addBranch (void)
	//////////////////////////////////////////////////////////////////////////// /Add Tree COINCBOARD variables Tree
void	LASEROBJ_addBranch (void)
void	TRIGGER_clearBranch (void)
	//////////////////////////////////////////////////////////////////////////// Clear Tree TRIGGER variables
void	MUON_clearBranch (void)
	//////////////////////////////////////////////////////////////////////////// Clear Tree MUON variables
void	ECAL_clearBranch (void)
	//////////////////////////////////////////////////////////////////////////// Clear Tree ECAL variables
void	QDC_clearBranch (void)
	//////////////////////////////////////////////////////////////////////////// Clear Tree QDC variables
void	LASER_clearBranch (void)
	//////////////////////////////////////////////////////////////////////////// Clear Tree LASER variables
void	ADDER_clearBranch (void)
	//////////////////////////////////////////////////////////////////////////// Clear Tree ADDER variables
void	CISPAR_clearBranch (void)
	//////////////////////////////////////////////////////////////////////////// Clear Tree CISPAR variables
void	BEAM_clearBranch (void)
	//////////////////////////////////////////////////////////////////////////// Clear Tree BEAM variables
void	DIGI_clearBranch (void)
	//////////////////////////////////////////////////////////////////////////// /Clear Tree DIGI variables ////////////////////////////////////////////////////////////////////////////
void	FELIX_clearBranch (void)
void	HIT_clearBranch (void)
	//////////////////////////////////////////////////////////////////////////// /Clear Tree HIT variables ////////////////////////////////////////////////////////////////////////////
void	ENETOTAL_clearBranch (void)
	//////////////////////////////////////////////////////////////////////////// Clear Tree ENETOTAL variables
void	COINCBOARD_clearBranch (void)
	//////////////////////////////////////////////////////////////////////////// Clear Tree COINCBOARD variables
void	LASEROBJ_clearBranch (void)
template<typename T>
void	clear_init_minus1 (std::vector< T > &vec)
template<typename T, size_t N>
void	clear_init_minus1 (std::vector< std::array< T, N > > &vec)
template<typename T, size_t N>
void	clear_init_zero (std::vector< std::array< T, N > > &vec)
void	clear_samples (std::vector< std::unique_ptr< int[]> > &vec, const std::vector< int > &nsamples, int nchan=MAX_CHAN)
int	digiChannel2PMT (int fragType, int chan)
short	CheckDMUFormat (uint32_t header)
	bit_31 of the DMU header must be 1 and bit_17 of the DMU header must be 0
short	CheckDMUParity (uint32_t header)
	Parity of the DMU header should be odd.
void	checkIsPropertySetup (float property, const std::string &name)
void	setupPropertyDefaultValue (float property, float defaultValue, const std::string &name)
void	setupBeamChambersBeforeTB2015 (void)
void	setupBeamChambersTB2015 (void)
void	setupBeamChambersTB2016_2020 (void)
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
SG::ReadCondHandleKey< TileSamplingFraction >	m_samplingFractionKey
	Name of TileSamplingFraction in condition store.
SG::ReadHandleKey< TileDigitsContainer >	m_digitsContainerKey
SG::ReadHandleKey< TileDigitsContainer >	m_digitsContainerFlxKey
SG::ReadHandleKey< TileBeamElemContainer >	m_beamElemContainerKey
SG::ReadHandleKey< TileRawChannelContainer >	m_flatRawChannelContainerKey
SG::ReadHandleKey< TileRawChannelContainer >	m_fitRawChannelContainerKey
SG::ReadHandleKey< TileRawChannelContainer >	m_optRawChannelContainerKey
SG::ReadHandleKey< TileRawChannelContainer >	m_dspRawChannelContainerKey
SG::ReadHandleKey< TileRawChannelContainer >	m_fitcRawChannelContainerKey
SG::ReadHandleKey< TileRawChannelContainer >	m_flxFitRawChannelContainerKey
SG::ReadHandleKey< TileRawChannelContainer >	m_flxOptRawChannelContainerKey
SG::ReadHandleKey< TileLaserObject >	m_laserObjectKey
SG::ReadHandleKey< TileHitContainer >	m_hitContainerKey
SG::ReadHandleKey< TileHitVector >	m_hitVectorKey
SG::ReadHandleKey< CaloCellContainer >	m_cellContainerKey
ToolHandle< TileCondToolEmscale >	m_tileToolEmscale
ToolHandle< TileRawChannelBuilderFlatFilter >	m_adderFilterAlgTool
SG::ReadHandleKey< TileDQstatus >	m_dqStatusKey
Gaudi::Property< bool >	m_calibrateEnergy {this, "CalibrateEnergy", true, "Calibrate energy"}
Gaudi::Property< bool >	m_useDspUnits {this, "UseDspUnits", false, "Use DSP untis"}
Gaudi::Property< int >	m_finalUnit {this, "OfflineUnits", TileRawChannelUnit::MegaElectronVolts, "Calibrate everything to this level"}
Gaudi::Property< bool >	m_calibMode {this, "CalibMode", false, "If data should be put in calib mode"}
Gaudi::Property< bool >	m_unpackAdder {this, "UnpackAdder", false, "Unpack adder"}
Gaudi::Property< bool >	m_completeNtuple {this, "CompleteNtuple", true, "Complete the ntuple"}
Gaudi::Property< bool >	m_bsInput {this, "BSInput", true, "Bytestream input"}
Gaudi::Property< bool >	m_pmtOrder {this, "PMTOrder", true, "Change channel ordering to pmt ordering in the ntuple"}
Gaudi::Property< int >	m_nSamples {this, "NSamples", NOT_SETUP, "Number of samples"}
Gaudi::Property< int >	m_nSamplesFlx {this, "NSamplesFelix", NOT_SETUP, "Number of samples for FELIX"}
Gaudi::Property< unsigned int >	m_nDrawers {this, "NDrawers", 6, "Number of drawers"}
Gaudi::Property< unsigned int >	m_nDrawersFlx {this, "NDrawersFelix", 0, "Number of drawers for FELIX"}
Gaudi::Property< int >	m_TBperiod {this, "TBperiod", 2016, "Tuned for 2016 testbeam by default"}
Gaudi::Property< int >	m_eventsPerFile {this, "EventsPerFile", 200000, "Number of events per file"}
Gaudi::Property< Long64_t >	m_treeSize {this, "TreeSize", 16000000000LL, "Size of tree"}
Gaudi::Property< std::string >	m_streamName {this, "StreamName", "AANT", "Name of the output stream"}
Gaudi::Property< std::string >	m_ntupleID {this, "NTupleID", "h1000", "Name of the ntuple ID"}
Gaudi::Property< std::vector< std::string > >	m_rosName {this, "rosName", {"B", "A", "C", "D", "E"}, "Name of arrays in ntuple for different ROSes"}
Gaudi::Property< std::vector< std::string > >	m_drawerList {this, "drawerList", {"-1"}, "List of frag IDs in correct order; Setup drawer list from data"}
Gaudi::Property< std::vector< int > >	m_drawerType {this, "drawerType", {}, "Type of every drawer 1-4: B+, B-, EB+, EB-; Take drawer type from Frag ID (doesn't work for 2003)"}
Gaudi::Property< std::vector< std::string > >	m_beamFragList {this, "beamFragList", {}, "List of beam frag IDs to store in the ntuple"}
Gaudi::Property< float >	m_beamBN2X1 {this, "BN2X1", 0.0, "Params for Beam TDC: Beam chamber: -2"}
Gaudi::Property< float >	m_beamBN2X2 {this, "BN2X2", 0.2, "Params for Beam TDC: Beam chamber: -2"}
Gaudi::Property< float >	m_beamBN2Y1 {this, "BN2Y1", 0.0, "Params for Beam TDC: Beam chamber: -2"}
Gaudi::Property< float >	m_beamBN2Y2 {this, "BN2Y2", 0.2, "Params for Beam TDC: Beam chamber: -2"}
Gaudi::Property< float >	m_beamBN1X1 {this, "BN1X1", 0.0, "Params for Beam TDC: Beam chamber: -1"}
Gaudi::Property< float >	m_beamBN1X2 {this, "BN1X2", 0.2, "Params for Beam TDC: Beam chamber: -1"}
Gaudi::Property< float >	m_beamBN1Y1 {this, "BN1Y1", 0.0, "Params for Beam TDC: Beam chamber: -1"}
Gaudi::Property< float >	m_beamBN1Y2 {this, "BN1Y2", 0.2, "Params for Beam TDC: Beam chamber: -1"}
Gaudi::Property< float >	m_beamBC0X1 {this, "BC0X1", 0.0, "Params for Beam TDC: Beam chamber: 0"}
Gaudi::Property< float >	m_beamBC0X2 {this, "BC0X2", 0.2, "Params for Beam TDC: Beam chamber: 0"}
Gaudi::Property< float >	m_beamBC0Y1 {this, "BC0Y1", 0.0, "Params for Beam TDC: Beam chamber: 0"}
Gaudi::Property< float >	m_beamBC0Y2 {this, "BC0Y2", 0.2, "Params for Beam TDC: Beam chamber: 0"}
Gaudi::Property< float >	m_beamBC0Z {this, "BC0Z", 17138.0, "Params for Beam TDC: Beam chamber: 0"}
Gaudi::Property< float >	m_beamBC1X1 {this, "BC1X1", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}
Gaudi::Property< float >	m_beamBC1X2 {this, "BC1X2", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}
Gaudi::Property< float >	m_beamBC1Y1 {this, "BC1Y1", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}
Gaudi::Property< float >	m_beamBC1Y2 {this, "BC1Y2", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}
Gaudi::Property< float >	m_beamBC1Z {this, "BC1Z", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}
Gaudi::Property< float >	m_beamBC1Z_0 {this, "BC1Z_0", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}
Gaudi::Property< float >	m_beamBC1Z_90 {this, "BC1Z_90", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}
Gaudi::Property< float >	m_beamBC1Z_min90 {this, "BC1Z_min90", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}
Gaudi::Property< float >	m_beamBC2X1 {this, "BC2X1", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}
Gaudi::Property< float >	m_beamBC2X2 {this, "BC2X2", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}
Gaudi::Property< float >	m_beamBC2Y1 {this, "BC2Y1", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}
Gaudi::Property< float >	m_beamBC2Y2 {this, "BC2Y2", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}
Gaudi::Property< float >	m_beamBC2Z {this, "BC2Z", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}
Gaudi::Property< float >	m_beamBC2Z_0 {this, "BC2Z_0", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}
Gaudi::Property< float >	m_beamBC2Z_90 {this, "BC2Z_90", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}
Gaudi::Property< float >	m_beamBC2Z_min90 {this, "BC2Z_min90", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}
Gaudi::Property< float >	m_radius {this, "Radius", 2280.0, "Inner radius of calo, for CTB 2004 only"}
Gaudi::Property< std::string >	m_etaFileName {this, "EtaFileName", "TileEtaCTB.txt", "File name with ETA, for CTB 2004 only"}
ServiceHandle< ITHistSvc >	m_thistSvc
TTree *	m_ntuplePtr {nullptr}
bool	m_ntupleCreated {false}
int	m_evtNr {-1}
int	m_evTime {}
int	m_run {}
int	m_evt {}
int	m_trigType {}
int	m_dspFlags {}
std::vector< int >	m_l1ID
std::vector< int >	m_l1Type
std::vector< int >	m_evBCID
std::vector< int >	m_evType
std::vector< int >	m_frBCID
float	m_muBackHit {}
float	m_muBackSum {}
std::array< float, 14 >	m_muBack {}
std::array< float, 2 >	m_muCalib {}
std::array< float, 8 >	m_ecal {}
std::array< uint32_t, 33 >	m_qdc {}
int	m_las_BCID {}
int	m_las_Filt {}
double	m_las_ReqAmp {}
double	m_las_MeasAmp {}
int	m_las_D1_ADC {}
int	m_las_D2_ADC {}
int	m_las_D3_ADC {}
int	m_las_D4_ADC {}
double	m_las_D1_Ped {}
double	m_las_D2_Ped {}
double	m_las_D3_Ped {}
double	m_las_D4_Ped {}
double	m_las_D1_Ped_RMS {}
double	m_las_D2_Ped_RMS {}
double	m_las_D3_Ped_RMS {}
double	m_las_D4_Ped_RMS {}
double	m_las_D1_Alpha {}
double	m_las_D2_Alpha {}
double	m_las_D3_Alpha {}
double	m_las_D4_Alpha {}
double	m_las_D1_Alpha_RMS {}
double	m_las_D2_Alpha_RMS {}
double	m_las_D3_Alpha_RMS {}
double	m_las_D4_Alpha_RMS {}
double	m_las_D1_AlphaPed {}
double	m_las_D2_AlphaPed {}
double	m_las_D3_AlphaPed {}
double	m_las_D4_AlphaPed {}
double	m_las_D1_AlphaPed_RMS {}
double	m_las_D2_AlphaPed_RMS {}
double	m_las_D3_AlphaPed_RMS {}
double	m_las_D4_AlphaPed_RMS {}
int	m_las_PMT1_ADC {}
int	m_las_PMT2_ADC {}
int	m_las_PMT1_TDC {}
int	m_las_PMT2_TDC {}
double	m_las_PMT1_Ped {}
double	m_las_PMT2_Ped {}
double	m_las_PMT1_Ped_RMS {}
double	m_las_PMT2_Ped_RMS {}
double	m_las_Temperature {}
int	m_lasFlag {}
float	m_las0 {}
float	m_las1 {}
float	m_las2 {}
float	m_las3 {}
std::array< float, 4 >	m_lasExtra {}
int	m_commonPU {}
std::vector< int * >	m_adder
std::vector< int >	m_adderPayload
std::array< float, 16 >	m_eneAdd {}
std::array< float, 16 >	m_timeAdd {}
int	m_cispar [16] {}
uint32_t	m_s1cou {}
uint32_t	m_s2cou {}
uint32_t	m_s3cou {}
uint32_t	m_cher1 {}
uint32_t	m_cher2 {}
uint32_t	m_cher3 {}
uint32_t	m_muTag {}
uint32_t	m_muHalo {}
uint32_t	m_muVeto {}
int	m_s2extra {}
int	m_s3extra {}
int	m_sc1 {}
int	m_sc2 {}
std::array< int, 16 >	m_tof {}
std::array< int, 16 >	m_btdc1 {}
std::array< int, 16 >	m_btdc2 {}
std::array< int, 16 >	m_scaler {}
std::vector< std::vector< int > > *	m_btdc {}
int	m_tjitter {}
int	m_tscTOF {}
int	m_btdcNhit [16] {}
int	m_btdcNchMultiHit [2] {}
float	m_xChN2 {}
float	m_yChN2 {}
float	m_xChN1 {}
float	m_yChN1 {}
float	m_xCha0 {}
float	m_yCha0 {}
float	m_xCha1 {}
float	m_yCha1 {}
float	m_xCha2 {}
float	m_yCha2 {}
float	m_xCha1_0 {}
float	m_yCha1_0 {}
float	m_xCha2_0 {}
float	m_yCha2_0 {}
float	m_xImp {}
float	m_yImp {}
float	m_xImp_0 {}
float	m_yImp_0 {}
float	m_xImp_90 {}
float	m_yImp_90 {}
float	m_xImp_min90 {}
float	m_yImp_min90 {}
std::vector< int >	m_evtVec
std::vector< short >	m_rodBCIDVec
std::vector< short >	m_sizeVec
std::vector< int >	m_evtflxVec
std::vector< short >	m_rodBCIDflxVec
std::vector< short >	m_sizeflxVec
std::vector< std::array< int, MAX_CHAN > >	m_gainflxVec
std::vector< std::unique_ptr< int[]> >	m_sampleflxVec
std::vector< std::array< int, MAX_DMU > >	m_bcidVec
std::vector< std::array< uint32_t, MAX_DMU > >	m_DMUheaderVec
std::vector< std::array< short, MAX_DMU > >	m_DMUformatErrVec
std::vector< std::array< short, MAX_DMU > >	m_DMUparityErrVec
std::vector< std::array< short, MAX_DMU > >	m_DMUmemoryErrVec
std::vector< std::array< short, MAX_DMU > >	m_DMUDstrobeErrVec
std::vector< std::array< short, MAX_DMU > >	m_DMUSstrobeErrVec
std::vector< std::array< short, MAX_DMU > >	m_DMUHeadparityErrVec
std::vector< std::array< short, MAX_DMU > >	m_DMUDataparityErrVec
std::vector< std::array< int, 2 > >	m_dmuMaskVec
std::vector< std::array< int, 2 > >	m_slinkCRCVec
std::vector< std::array< int, MAX_CHAN > >	m_gainVec
std::vector< std::unique_ptr< int[]> >	m_sampleVec
std::vector< std::array< int, MAX_DMU > >	m_feCRCVec
std::vector< std::array< int, MAX_DMU > >	m_rodCRCVec
std::vector< std::array< float, MAX_CHAN > >	m_eneVec
std::vector< std::array< float, MAX_CHAN > >	m_timeVec
std::vector< std::array< float, MAX_CHAN > >	m_pedFlatVec
std::vector< std::array< float, MAX_CHAN > >	m_chi2FlatVec
std::vector< std::array< float, MAX_CHAN > >	m_efitVec
std::vector< std::array< float, MAX_CHAN > >	m_tfitVec
std::vector< std::array< float, MAX_CHAN > >	m_pedfitVec
std::vector< std::array< float, MAX_CHAN > >	m_chi2Vec
std::vector< std::array< float, MAX_CHAN > >	m_efitcVec
std::vector< std::array< float, MAX_CHAN > >	m_tfitcVec
std::vector< std::array< float, MAX_CHAN > >	m_pedfitcVec
std::vector< std::array< float, MAX_CHAN > >	m_chi2cVec
std::vector< std::array< float, MAX_CHAN > >	m_eOptVec
std::vector< std::array< float, MAX_CHAN > >	m_tOptVec
std::vector< std::array< float, MAX_CHAN > >	m_pedOptVec
std::vector< std::array< float, MAX_CHAN > >	m_chi2OptVec
std::vector< std::array< float, MAX_CHAN > >	m_eDspVec
std::vector< std::array< float, MAX_CHAN > >	m_tDspVec
std::vector< std::array< float, MAX_CHAN > >	m_chi2DspVec
std::vector< std::array< float, MAX_CHAN > >	m_eflxfitVec
std::vector< std::array< float, MAX_CHAN > >	m_tflxfitVec
std::vector< std::array< float, MAX_CHAN > >	m_chi2flxfitVec
std::vector< std::array< float, MAX_CHAN > >	m_pedflxfitVec
std::vector< std::array< float, MAX_CHAN > >	m_eflxoptVec
std::vector< std::array< float, MAX_CHAN > >	m_tflxoptVec
std::vector< std::array< float, MAX_CHAN > >	m_chi2flxoptVec
std::vector< std::array< float, MAX_CHAN > >	m_pedflxoptVec
std::vector< short >	m_ROD_GlobalCRCVec
std::vector< std::array< short, MAX_DMU > >	m_ROD_DMUBCIDVec
std::vector< std::array< short, MAX_DMU > >	m_ROD_DMUmemoryErrVec
std::vector< std::array< short, MAX_DMU > >	m_ROD_DMUSstrobeErrVec
std::vector< std::array< short, MAX_DMU > >	m_ROD_DMUDstrobeErrVec
std::vector< std::array< short, MAX_DMU > >	m_ROD_DMUHeadformatErrVec
std::vector< std::array< short, MAX_DMU > >	m_ROD_DMUDataformatErrVec
std::vector< std::array< short, 2 > >	m_ROD_DMUMaskVec
std::vector< std::array< int, MAX_MINIDRAWER > >	m_mdL1idflxVec
std::vector< std::array< int, MAX_MINIDRAWER > >	m_mdBcidflxVec
std::vector< std::array< int, MAX_MINIDRAWER > >	m_mdModuleflxVec
std::vector< std::array< int, MAX_MINIDRAWER > >	m_mdRunTypeflxVec
std::vector< std::array< int, MAX_MINIDRAWER > >	m_mdPedLoflxVec
std::vector< std::array< int, MAX_MINIDRAWER > >	m_mdPedHiflxVec
std::vector< std::array< int, MAX_MINIDRAWER > >	m_mdRunflxVec
std::vector< std::array< int, MAX_MINIDRAWER > >	m_mdChargeflxVec
std::vector< std::array< int, MAX_MINIDRAWER > >	m_mdChargeTimeflxVec
std::vector< std::array< int, MAX_MINIDRAWER > >	m_mdCapacitorflxVec
std::array< float, 4 >	m_LarEne {}
std::array< float, 3 >	m_BarEne {}
std::array< float, 3 >	m_ExtEne {}
std::array< float, 3 >	m_GapEne {}
std::array< unsigned int, 96 >	m_coincTrig1 {}
std::array< unsigned int, 96 >	m_coincTrig2 {}
std::array< unsigned int, 96 >	m_coincTrig3 {}
std::array< unsigned int, 96 >	m_coincTrig4 {}
std::array< unsigned int, 96 >	m_coincTrig5 {}
std::array< unsigned int, 96 >	m_coincTrig6 {}
std::array< unsigned int, 96 >	m_coincTrig7 {}
std::array< unsigned int, 96 >	m_coincTrig8 {}
int	m_coincFlag1 {}
int	m_coincFlag2 {}
int	m_coincFlag3 {}
int	m_coincFlag4 {}
int	m_coincFlag5 {}
int	m_coincFlag6 {}
int	m_coincFlag7 {}
int	m_coincFlag8 {}
std::map< unsigned int, unsigned int, std::less< unsigned int > >	m_drawerMap
std::map< unsigned int, unsigned int, std::less< unsigned int > >	m_drawerFlxMap
bool	m_beamIdList [32]
int	m_runNumber {}
float	m_eta {}
float	m_theta {}
std::vector< std::array< float, MAX_CHAN > >	m_ehitVec
std::vector< std::array< float, MAX_CHAN > >	m_thitVec
std::vector< std::array< float, MAX_CHAN > >	m_ehitCnt
std::vector< std::array< float, MAX_CHAN > >	m_thitCnt
bool	m_calibrateEnergyThisEvent {false}
TileRawChannelUnit::UNIT	m_rchUnit
	Unit for TileRawChannels (ADC, pCb, MeV).
TileRawChannelUnit::UNIT	m_dspUnit
	Unit for TileRawChannels in DSP.
const TileID *	m_tileID {nullptr}
const TileHWID *	m_tileHWID {nullptr}
const TileCablingService *	m_cabling {nullptr}
std::map< int, int >	m_nSamplesInDrawerMap
std::map< int, int >	m_nSamplesFlxInDrawerMap
std::vector< int >	m_nSamplesInDrawer
std::vector< int >	m_nSamplesFlxInDrawer
bool	m_saveFelixData {false}
DataObjIDColl	m_extendedExtraObjects
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default).
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default).
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

class to produce TileCal commissioning ntuples

Definition at line 89 of file TileTBAANtuple.h.

Member Typedef Documentation

drawerMap_iterator

typedef std::map<unsignedint,unsignedint,std::less<unsignedint>>::iterator TileTBAANtuple::drawerMap_iterator

private

Definition at line 605 of file TileTBAANtuple.h.

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

private

Enumerator
NOT_SETUP

Definition at line 105 of file TileTBAANtuple.h.

{NOT_SETUP = -9999};

anonymous enum

anonymous enum

private

Enumerator
MAX_MINIDRAWER
MAX_CHAN
MAX_DMU

Definition at line 106 of file TileTBAANtuple.h.

{MAX_MINIDRAWER = 4, MAX_CHAN = 48, MAX_DMU = 16};

Constructor & Destructor Documentation

TileTBAANtuple()

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

Definition at line 82 of file TileTBAANtuple.cxx.

  : AthAlgorithm(name, pSvcLocator)
  , m_thistSvc("THistSvc", name)
  , m_rchUnit(TileRawChannelUnit::MegaElectronVolts)
  , m_dspUnit(TileRawChannelUnit::ADCcounts)
{
 
  char frg[6] = "0x000";
  m_beamFragList.clear();
  for (unsigned int i=0; i<sizeof(m_beamIdList)/sizeof(bool); ++i) {
    m_beamIdList[i] = false;
    // no coins trig by default (this is the case for 2004)
    if (i <= ECAL_ADC_FRAG || i >= COMMON_ADC1_FRAG ) {
      sprintf(frg,"0x%3.3x",i);
      m_beamFragList.value().push_back((std::string)frg);
    }
  }
 
}

~TileTBAANtuple()

virtual TileTBAANtuple::~TileTBAANtuple ( )

virtualdefault

Member Function Documentation

ADDER_addBranch()

void TileTBAANtuple::ADDER_addBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// /Add Tree ADDER variables Tree

////////////////////////////////////////////////////////////////////////////

Definition at line 2762 of file TileTBAANtuple.cxx.

                                         {
 
  if (m_unpackAdder) {
    if (m_beamIdList[ADD_FADC_FRAG]) {
      m_adder.resize(16);
      m_adderPayload.resize(16*16);
      for (int j = 0; j < 16; j++) {
        m_adder[j] = m_adderPayload.data() + j * 16;
      }
 
      m_ntuplePtr->Branch("Adder", *m_adder.data(), "m_adder[16][16]/I");
      m_ntuplePtr->Branch("EneAdd", &m_eneAdd, "m_eneAdd[16]/F");
      m_ntuplePtr->Branch("TimeAdd", &m_timeAdd, "m_timeAdd[16]/F");
    }
  }
}

ADDER_clearBranch()

void TileTBAANtuple::ADDER_clearBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// Clear Tree ADDER variables

////////////////////////////////////////////////////////////////////////////

Definition at line 2785 of file TileTBAANtuple.cxx.

{
 
 
}

BEAM_addBranch()

void TileTBAANtuple::BEAM_addBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// /Add Tree BEAM variables Tree

////////////////////////////////////////////////////////////////////////////

Definition at line 2823 of file TileTBAANtuple.cxx.

                                        {
 
    if (((m_TBperiod >= 2015 || m_unpackAdder) && m_beamIdList[BEAM_ADC_FRAG])
    || ((!m_unpackAdder || m_TBperiod >= 2022) && m_beamIdList[COMMON_ADC1_FRAG])) {
 
    m_ntuplePtr->Branch("S1cou", &m_s1cou, "S1cou/S");
    m_ntuplePtr->Branch("S2cou", &m_s2cou, "S2cou/S");
    if (m_TBperiod >= 2016 && m_TBperiod < 2022) {
      m_ntuplePtr->Branch("Cher3", &m_s3cou, "Cher3/S"); // dropped S3 in favor of Cher3 in September 2016 TB
     // m_ntuplePtr->Branch("SpmFinger", &m_muVeto, "SmpFinger/S"); //  muVeto replaced by SPM for testing
      m_ntuplePtr->Branch("SiPM1", &m_muTag, "SiPM1/S"); //  muTag replaced by SiPM1 for testing
      m_ntuplePtr->Branch("SiPM2", &m_muHalo, "SiPM2/S");// muHalo repalce by SiPM1 for testing
    } else {
      m_ntuplePtr->Branch("S3cou", &m_s3cou, "S3cou/S");
    }
 
    m_ntuplePtr->Branch("Cher1", &m_cher1, "Cher1/S");
    m_ntuplePtr->Branch("Cher2", &m_cher2, "Cher2/S");
    if (m_TBperiod >= 2022) {
      m_ntuplePtr->Branch("Cher3", &m_cher3, "Cher3/S");
    }
 
    if (m_TBperiod < 2015) {
      m_ntuplePtr->Branch("MuTag", &m_muTag, "MuTag/S");
      m_ntuplePtr->Branch("MuHalo", &m_muHalo, "MuHalo/S");
      m_ntuplePtr->Branch("MuVeto", &m_muVeto, "MuVeto/S");
    //} else if (m_TBperiod == 2015) {
      // nothing
    } else if (m_TBperiod >= 2016) {
      m_ntuplePtr->Branch("SCalo1", &m_muTag, "SCalo1/S");
      m_ntuplePtr->Branch("SCalo2", &m_muHalo, "SCalo2/S");
      //m_ntuplePtr->Branch("SCalo3", &m_muVeto, "SCalo3/S");
    }
  }
 
  if (m_TBperiod >= 2015) {
    if (m_beamIdList[COMMON_TDC1_FRAG] || (m_TBperiod >= 2022 && m_beamIdList[COMMON_TOF_FRAG])) {
      m_btdc = new std::vector<std::vector<int> >(16);
      m_ntuplePtr->Branch("btdc1", &m_btdc1, "m_btdc1[16]/I");
      m_ntuplePtr->Branch("btdc2", &m_btdc2, "m_btdc2[16]/I");
      m_ntuplePtr->Branch("btdc", &m_btdc);
      m_ntuplePtr->Branch("tjitter", &m_tjitter, "tjitter/I");
      m_ntuplePtr->Branch("tscTOF", &m_tscTOF, "tscTOF/I");
      m_ntuplePtr->Branch("btdcNhit", m_btdcNhit, "btdcNhit[16]/I");
      m_ntuplePtr->Branch("btdcNchMultiHit", m_btdcNchMultiHit, "btdcNchMultiHit[2]/I");
    }
    if (m_beamIdList[COMMON_TOF_FRAG]) {
      m_ntuplePtr->Branch("tof", &m_tof, "m_tof[16]/I");
    }
    if (m_beamIdList[COMMON_PTN_FRAG]) {
      m_ntuplePtr->Branch("scaler", &m_scaler, "m_scaler[16]/I");
    }
  } else if (!m_unpackAdder) {
    if (m_beamIdList[COMMON_ADC2_FRAG]) {
      m_ntuplePtr->Branch("S2extra", &m_s2extra, "S2extra/S");
      m_ntuplePtr->Branch("S2cou", &m_s3extra, "S3extra/S");
    }
    if (m_beamIdList[BEAM_ADC_FRAG]) {
      m_ntuplePtr->Branch("SC1", &m_sc1, "SC1/S");
      m_ntuplePtr->Branch("SC2", &m_sc2, "SC2/S");
    }
    if (m_beamIdList[COMMON_PTN_FRAG]) {
      m_ntuplePtr->Branch("pu", &m_commonPU, "pu/S");
    }
    if (m_beamIdList[COMMON_TOF_FRAG]) {
      m_ntuplePtr->Branch("tof", &m_tof, "m_tof[8]/I");
    }
    if (m_beamIdList[COMMON_TDC1_FRAG]) {
      m_ntuplePtr->Branch("btdc1", &m_btdc1, "m_btdc1/I");
    }
    if (m_beamIdList[COMMON_TDC2_FRAG]) {
      m_ntuplePtr->Branch("btdc2", &m_btdc2, "m_btdc2/I");
    }
 
    if (m_beamIdList[COMMON_TDC1_FRAG]) {
      m_ntuplePtr->Branch("XchN2", &m_xChN2, "XchN2/F");
      m_ntuplePtr->Branch("YchN2", &m_yChN2, "YchN2/F");
      m_ntuplePtr->Branch("XchN1", &m_xChN1, "XchN1/F");
      m_ntuplePtr->Branch("YchN1", &m_yChN1, "YchN1/F");
    }
    if (m_beamIdList[COMMON_TDC2_FRAG]) {
      m_ntuplePtr->Branch("Xcha0", &m_xCha0, "Xcha0/F");
      m_ntuplePtr->Branch("Ycha0", &m_yCha0, "Ycha0/F");
    }
  } else {
    if (m_beamIdList[BEAM_TDC_FRAG]) {
      m_ntuplePtr->Branch("btdc1", &m_btdc1, "m_btdc1[16]/I");
      m_ntuplePtr->Branch("btdc2", &m_btdc2, "m_btdc2[16]/I");
    }
  }
 
  if ((m_TBperiod >= 2015 && m_beamIdList[COMMON_TDC1_FRAG])
      || (m_unpackAdder && m_beamIdList[BEAM_TDC_FRAG])
      || (!m_unpackAdder && m_beamIdList[COMMON_TDC2_FRAG])
      || (m_TBperiod >= 2022 && m_beamIdList[COMMON_TOF_FRAG])) {
 
    m_ntuplePtr->Branch("Xcha1", &m_xCha1, "Xcha1/F");
    m_ntuplePtr->Branch("Ycha1", &m_yCha1, "Ycha1/F");
    m_ntuplePtr->Branch("Xcha2", &m_xCha2, "Xcha2/F");
    m_ntuplePtr->Branch("Ycha2", &m_yCha2, "Ycha2/F");
    m_ntuplePtr->Branch("Xcha1_0", &m_xCha1_0, "Xcha1_0/F");
    m_ntuplePtr->Branch("Ycha1_0", &m_yCha1_0, "Ycha1_0/F");
    m_ntuplePtr->Branch("Xcha2_0", &m_xCha2_0, "Xcha2_0/F");
    m_ntuplePtr->Branch("Ycha2_0", &m_yCha2_0, "Ycha2_0/F");
    m_ntuplePtr->Branch("Ximp", &m_xImp, "Ximp/F");
    m_ntuplePtr->Branch("Yimp", &m_yImp, "Yimp/F");
    m_ntuplePtr->Branch("Ximp_0", &m_xImp_0, "Ximp_0/F");
    m_ntuplePtr->Branch("Yimp_0", &m_yImp_0, "Yimp_0/F");
    m_ntuplePtr->Branch("Ximp_90", &m_xImp_90, "Ximp_90/F");
    m_ntuplePtr->Branch("Yimp_90", &m_yImp_90, "Yimp_90/F");
    m_ntuplePtr->Branch("Ximp_min90", &m_xImp_min90, "Ximp_min90/F");
    m_ntuplePtr->Branch("Yimp_min90", &m_yImp_min90, "Yimp_min90/F");
  }
 
}

BEAM_clearBranch()

void TileTBAANtuple::BEAM_clearBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// Clear Tree BEAM variables

////////////////////////////////////////////////////////////////////////////

Definition at line 2945 of file TileTBAANtuple.cxx.

                                          {
 
  if ((m_unpackAdder && m_beamIdList[BEAM_ADC_FRAG])
      || (!m_unpackAdder && m_beamIdList[COMMON_ADC1_FRAG])) {
 
    m_s1cou = 0;
    m_s2cou = 0;
    m_s3cou = 0;
    m_cher1 = 0;
    m_cher2 = 0;
    m_cher3 = 0;
    m_muTag = 0;
    m_muHalo = 0;
    m_muVeto = 0;
  }
  if (!m_unpackAdder) {
    if (m_beamIdList[COMMON_ADC2_FRAG]) {
 
      m_s2extra = 0;
      m_s3extra = 0;
    }
    if (m_beamIdList[BEAM_ADC_FRAG]) {
      m_sc1 = 0;
      m_sc2 = 0;
    }
    if (m_beamIdList[COMMON_PTN_FRAG]) {
      m_commonPU = 0;
    }
    if (m_beamIdList[COMMON_TDC1_FRAG]) {
      m_xChN2 = 0.;
      m_yChN2 = 0.;
      m_xChN1 = 0.;
      m_yChN1 = 0.;
    }
    if (m_beamIdList[COMMON_TDC2_FRAG]) {
      m_xCha0 = 0.;
      m_yCha0 = 0.;
    }
  }
 
  if ((m_TBperiod >= 2015 && m_beamIdList[COMMON_TDC1_FRAG])
      || (m_unpackAdder && m_beamIdList[BEAM_TDC_FRAG])
      || (!m_unpackAdder && m_beamIdList[COMMON_TDC2_FRAG])
      || (m_TBperiod >= 2022 && m_beamIdList[COMMON_TOF_FRAG])) {
 
    m_xCha1 = 0.;
    m_yCha1 = 0.;
    m_xCha2 = 0.;
    m_yCha2 = 0.;
    m_xCha1_0 = 0.;
    m_yCha1_0 = 0.;
    m_xCha2_0 = 0.;
    m_yCha2_0 = 0.;
    m_xImp = 0.;
    m_yImp = 0.;
  }
 
  m_muBack.fill(0.0F);
 
  for (int i=0; i<16; i+=2) {
    m_tof[i] = +0xFFFF;
    m_tof[i+1] = -0xFFFF;
  }
 
  for (int i=0; i<16; i+=2) {
    m_btdc1[i] = +0xFFFF;
    m_btdc1[i+1] = -0xFFFF;
  }
 
  for (int i=0; i<16; i+=2) {
    m_btdc2[i] = +0xFFFF;
    m_btdc2[i+1] = -0xFFFF;
  }
 
  if (m_btdc) {
    for (std::vector<int>& btdc_amplitudes : *m_btdc) {
      btdc_amplitudes.clear();
    }
  }
 
  memset(&m_scaler, 0, sizeof(m_scaler));
 
  memset(m_btdcNhit,0,sizeof(m_btdcNhit));
  memset(m_btdcNchMultiHit,0,sizeof(m_btdcNchMultiHit));
}

CheckDMUFormat()

short TileTBAANtuple::CheckDMUFormat ( uint32_t header )

inlineprivate

bit_31 of the DMU header must be 1 and bit_17 of the DMU header must be 0

Definition at line 302 of file TileTBAANtuple.h.

                                                 {
      if (((header >> 31 & 0x1) == 1) && ((header >> 17 & 0x1) == 0))
        return 0; // no error
      else
        return 1; //error
    }

CheckDMUParity()

short TileTBAANtuple::CheckDMUParity ( uint32_t header )

inlineprivate

Parity of the DMU header should be odd.

Definition at line 311 of file TileTBAANtuple.h.

                                                 {
      uint32_t parity(0);
      for (int i = 0; i < 32; ++i)
        parity += header >> i & 0x1;
 
      if ((parity % 2) == 1)
        return 0; //no error
      else
        return 1; //error
    }

checkIsPropertySetup()

void TileTBAANtuple::checkIsPropertySetup ( float property, const std::string & name )

inlineprivate

Definition at line 322 of file TileTBAANtuple.h.

                                                                     {
      if (property < NOT_SETUP) {
        ATH_MSG_ERROR("The following property should be set up via JO: " << name);
      }
    }

CISPAR_addBranch()

void TileTBAANtuple::CISPAR_addBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// /Add Tree CISPAR variables Tree

////////////////////////////////////////////////////////////////////////////

Definition at line 2797 of file TileTBAANtuple.cxx.

                                          {
 
  if (m_beamIdList[DIGI_PAR_FRAG & 0x1F]) {
 
    m_ntuplePtr->Branch("cispar", m_cispar, "cispar[16]/I");
  }
}

CISPAR_clearBranch()

void TileTBAANtuple::CISPAR_clearBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// Clear Tree CISPAR variables

////////////////////////////////////////////////////////////////////////////

Definition at line 2811 of file TileTBAANtuple.cxx.

{
  memset(m_cispar,-1,sizeof(m_cispar));
}

clear_init_minus1() [1/2]

template<typename T, size_t N>

void TileTBAANtuple::clear_init_minus1 ( std::vector< std::array< T, N > > & vec )

private

Definition at line 3752 of file TileTBAANtuple.cxx.

                                                                    {
  for (std::array<T,N>& arr : vec) {
    std::fill(arr.begin(), arr.end(), static_cast<T>(-1));
  }
}

clear_init_minus1() [2/2]

template<typename T>

void TileTBAANtuple::clear_init_minus1 ( std::vector< T > & vec )

private

Definition at line 3747 of file TileTBAANtuple.cxx.

                                                        {
  std::fill(vec.begin(), vec.end(), static_cast<T>(-1));
}

clear_init_zero()

template<typename T, size_t N>

void TileTBAANtuple::clear_init_zero ( std::vector< std::array< T, N > > & vec )

private

Definition at line 3759 of file TileTBAANtuple.cxx.

                                                                  {
  for (std::array<T,N>& arr : vec) {
    std::fill(arr.begin(), arr.end(), static_cast<T>(0));
  }
}

clear_samples()

void TileTBAANtuple::clear_samples ( std::vector< std::unique_ptr< int[]> > & vec, const std::vector< int > & nsamples, int nchan = MAX_CHAN )

private

Definition at line 3765 of file TileTBAANtuple.cxx.

{
  for (unsigned int i = 0; i < vec.size(); ++i) {
    std::fill(vec[i].get(), vec[i].get() + nsamples.at(i) * nchan, -1);
  }
}

COINCBOARD_addBranch()

void TileTBAANtuple::COINCBOARD_addBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// /Add Tree COINCBOARD variables Tree

////////////////////////////////////////////////////////////////////////////

Definition at line 3065 of file TileTBAANtuple.cxx.

                                              {
  if (m_beamIdList[COIN_TRIG1_FRAG]) {
    m_ntuplePtr->Branch("CoincTrig1", &m_coincTrig1, "m_coincTrig1[96]/I");
    m_ntuplePtr->Branch("CoincFlag1", &m_coincFlag1, "CoincFlag1/I");
  }
 
  if (m_beamIdList[COIN_TRIG2_FRAG]) {
    m_ntuplePtr->Branch("CoincTrig2", &m_coincTrig2, "m_coincTrig2[96]/I");
    m_ntuplePtr->Branch("CoincFlag2", &m_coincFlag2, "CoincFlag2/I");
  }
 
  if (m_beamIdList[COIN_TRIG3_FRAG]) {
    m_ntuplePtr->Branch("CoincTrig3", &m_coincTrig3, "m_coincTrig3[96]/I");
    m_ntuplePtr->Branch("CoincFlag3", &m_coincFlag3, "CoincFlag3/I");
  }
 
  if (m_beamIdList[COIN_TRIG4_FRAG]) {
    m_ntuplePtr->Branch("CoincTrig4", &m_coincTrig4, "m_coincTrig4[96]/I");
    m_ntuplePtr->Branch("CoincFlag4", &m_coincFlag4, "CoincFlag4/I");
  }
  if (m_beamIdList[COIN_TRIG5_FRAG]) {
    m_ntuplePtr->Branch("CoincTrig5", &m_coincTrig5, "m_coincTrig5[96]/I");
    m_ntuplePtr->Branch("CoincFlag5", &m_coincFlag5, "CoincFlag5/I");
  }
 
  if (m_beamIdList[COIN_TRIG6_FRAG]) {
    m_ntuplePtr->Branch("CoincTrig6", &m_coincTrig6, "m_coincTrig6[96]/I");
    m_ntuplePtr->Branch("CoincFlag6", &m_coincFlag6, "CoincFlag6/I");
  }
 
  if (m_beamIdList[COIN_TRIG7_FRAG]) {
    m_ntuplePtr->Branch("CoincTrig7", &m_coincTrig7, "m_coincTrig7[96]/I");
    m_ntuplePtr->Branch("CoincFlag7", &m_coincFlag7, "CoincFlag7/I");
  }
 
  if (m_beamIdList[COIN_TRIG8_FRAG]) {
    m_ntuplePtr->Branch("CoincTrig8", &m_coincTrig8, "m_coincTrig8[96]/I");
    m_ntuplePtr->Branch("CoincFlag8", &m_coincFlag8, "CoincFlag8/I");
  }
 
}

COINCBOARD_clearBranch()

void TileTBAANtuple::COINCBOARD_clearBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// Clear Tree COINCBOARD variables

////////////////////////////////////////////////////////////////////////////

Definition at line 3113 of file TileTBAANtuple.cxx.

{
 
 
 
}

declareGaudiProperty()

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

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

declareProperty()

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

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

detStore()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

DIGI_addBranch()

void TileTBAANtuple::DIGI_addBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// /Add Tree DIGI variables Tree

////////////////////////////////////////////////////////////////////////////

Reserve is needed, because, if the vector is displaced, then the position declared in the Branch i not longer valid. Can be done a cleaner way?

Definition at line 3126 of file TileTBAANtuple.cxx.

{
 

  m_bcidVec.reserve(MAX_DRAWERS);
  m_DMUheaderVec.reserve(MAX_DRAWERS);
  m_DMUformatErrVec.reserve(MAX_DRAWERS);
  m_DMUparityErrVec.reserve(MAX_DRAWERS);
  m_DMUmemoryErrVec.reserve(MAX_DRAWERS);
  m_DMUDstrobeErrVec.reserve(MAX_DRAWERS);
  m_DMUHeadparityErrVec.reserve(MAX_DRAWERS);
  m_DMUDataparityErrVec.reserve(MAX_DRAWERS);
  m_DMUSstrobeErrVec.reserve(MAX_DRAWERS);
  m_dmuMaskVec.reserve(MAX_DRAWERS);
  m_slinkCRCVec.reserve(MAX_DRAWERS);
  m_gainVec.reserve(MAX_DRAWERS);
  m_sampleVec.reserve(MAX_DRAWERS);
  m_feCRCVec.reserve(MAX_DRAWERS);
  m_rodCRCVec.reserve(MAX_DRAWERS);
  m_eneVec.reserve(MAX_DRAWERS);
  m_timeVec.reserve(MAX_DRAWERS);
  m_pedFlatVec.reserve(MAX_DRAWERS);
  m_chi2FlatVec.reserve(MAX_DRAWERS);
  m_efitVec.reserve(MAX_DRAWERS);
  m_tfitVec.reserve(MAX_DRAWERS);
  m_pedfitVec.reserve(MAX_DRAWERS);
  m_chi2Vec.reserve(MAX_DRAWERS);
  m_efitcVec.reserve(MAX_DRAWERS);
  m_tfitcVec.reserve(MAX_DRAWERS);
  m_pedfitcVec.reserve(MAX_DRAWERS);
  m_chi2cVec.reserve(MAX_DRAWERS);
  m_eOptVec.reserve(MAX_DRAWERS);
  m_tOptVec.reserve(MAX_DRAWERS);
  m_pedOptVec.reserve(MAX_DRAWERS);
  m_chi2OptVec.reserve(MAX_DRAWERS);
  m_eDspVec.reserve(MAX_DRAWERS);
  m_tDspVec.reserve(MAX_DRAWERS);
  m_chi2DspVec.reserve(MAX_DRAWERS);
  m_ROD_GlobalCRCVec.reserve(MAX_DRAWERS);
  m_ROD_DMUBCIDVec.reserve(MAX_DRAWERS);
  m_ROD_DMUmemoryErrVec.reserve(MAX_DRAWERS);
  m_ROD_DMUSstrobeErrVec.reserve(MAX_DRAWERS);
  m_ROD_DMUDstrobeErrVec.reserve(MAX_DRAWERS);
  m_ROD_DMUHeadformatErrVec.reserve(MAX_DRAWERS);
  m_ROD_DMUDataformatErrVec.reserve(MAX_DRAWERS);
  m_ROD_DMUMaskVec.reserve(MAX_DRAWERS);
 
  std::ostringstream oss;
  oss << m_nSamples;
  std::string nSampStr=oss.str();
 
  unsigned int listSize = std::min(m_nDrawers.value(), static_cast<unsigned int>(m_drawerMap.size()));
 
  if (listSize > 0) {
 
    std::string digit[10] = { "0", "1", "2", "3", "4", "5", "6", "7", "8", "9" };
    std::vector<std::string> suffixArr;
    unsigned int length;
    bool testbeam = TileCablingService::getInstance()->getTestBeam();
 
    if (m_calibMode) {
 
      length = 2 * m_nDrawers;
      suffixArr.resize(length);
 
      for (unsigned int i = 0; i < listSize; ++i) {
        unsigned int ros = m_drawerType[i];
        unsigned int drawer = strtol(m_drawerList[i].data(), NULL, 0) & 0x3F;
        std::string digits;
        if (m_TBperiod >= 2010) {
          ++drawer; // count modules from 1
          digits = digit[drawer / 10] + digit[drawer % 10];
        } else if (testbeam) {
          digits = digit[drawer & 7];
        } else {
          ++drawer; // count modules from 1
          digits = digit[drawer / 10] + digit[drawer % 10];
        }
 
        if (ros == 0) {
          std::string suff = m_drawerList[i];
          suff.replace(suff.find("0x"), 2, "");
          suffixArr[i] = suff + "lo";
          suffixArr[i + m_nDrawers] = suff + "hi";
        } else {
          suffixArr[i] = m_rosName[ros] + digits + "lo";
          suffixArr[i + m_nDrawers] = m_rosName[ros] + digits + "hi";
        }
      }
    } else {
 
      length = m_nDrawers;
      suffixArr.resize(length);
 
      for (unsigned int i = 0; i < listSize; ++i) {
        unsigned int ros = m_drawerType[i];
        unsigned int drawer = strtol(m_drawerList[i].data(), NULL, 0) & 0x3F;
        std::string digits;
        if (m_TBperiod >= 2010) {
          ++drawer; // count modules from 1
          digits = digit[drawer / 10] + digit[drawer % 10];
        } else if (testbeam) {
          digits = digit[drawer & 7];
        } else {
          ++drawer; // count modules from 1
          digits = digit[drawer / 10] + digit[drawer % 10];
        }
 
        if (ros == 0) {
          std::string suff = m_drawerList[i];
          suff.replace(suff.find("0x"), 2, "");
          suffixArr[i] = std::move(suff);
        } else {
          suffixArr[i] = m_rosName[ros] + digits;
        }
      }
    }
 
    m_nSamplesInDrawer.reserve(length);
 
    m_evtVec.resize(length);
    m_rodBCIDVec.resize(length);
    m_sizeVec.resize(length);
    m_ROD_GlobalCRCVec.resize(length);
 
    for (unsigned int i = 0; i < length; i++) {
 
      int nSamplesInDrawer(m_nSamples);
      int frag = std::stoi(m_drawerList[i%m_nDrawers], nullptr, 0);
 
      auto it = m_nSamplesInDrawerMap.find(frag);
      if (it != m_nSamplesInDrawerMap.end()) {
        nSamplesInDrawer = it->second;
      }
      m_nSamplesInDrawer.push_back(nSamplesInDrawer);
 
      m_bcidVec.push_back(std::array<int, MAX_DMU>()); // U
      m_DMUheaderVec.push_back(std::array<uint32_t, MAX_DMU>()); // U32
      m_DMUformatErrVec.push_back(std::array<short, MAX_DMU>()); // U32
      m_DMUparityErrVec.push_back(std::array<short, MAX_DMU>()); // U32
      m_DMUmemoryErrVec.push_back(std::array<short, MAX_DMU>()); // U32
      m_DMUDstrobeErrVec.push_back(std::array<short, MAX_DMU>()); // U32
      m_DMUSstrobeErrVec.push_back(std::array<short, MAX_DMU>()); // U32
      m_DMUHeadparityErrVec.push_back(std::array<short, MAX_DMU>());
      m_DMUDataparityErrVec.push_back(std::array<short, MAX_DMU>());
 
      m_dmuMaskVec.push_back(std::array<int, 2>()); // U(2)
      m_slinkCRCVec.push_back(std::array<int, 2>()); // U(2)
      m_gainVec.push_back(std::array<int, MAX_CHAN>()); // U(48/96)
      m_sampleVec.push_back(std::make_unique<int[]>(MAX_CHAN * nSamplesInDrawer)); // U(48/96,9)
      m_feCRCVec.push_back(std::array<int, MAX_DMU>()); //U
      m_rodCRCVec.push_back(std::array<int, MAX_DMU>()); //U
 
      m_eneVec.push_back(std::array<float, MAX_CHAN>());
      m_timeVec.push_back(std::array<float, MAX_CHAN>());
      m_pedFlatVec.push_back(std::array<float, MAX_CHAN>());
      m_chi2FlatVec.push_back(std::array<float, MAX_CHAN>());
 
      m_efitVec.push_back(std::array<float, MAX_CHAN>());
      m_tfitVec.push_back(std::array<float, MAX_CHAN>());
      m_pedfitVec.push_back(std::array<float, MAX_CHAN>());
      m_chi2Vec.push_back(std::array<float, MAX_CHAN>());
 
      m_efitcVec.push_back(std::array<float, MAX_CHAN>());
      m_tfitcVec.push_back(std::array<float, MAX_CHAN>());
      m_pedfitcVec.push_back(std::array<float, MAX_CHAN>());
      m_chi2cVec.push_back(std::array<float, MAX_CHAN>());
 
      m_eOptVec.push_back(std::array<float, MAX_CHAN>());
      m_tOptVec.push_back(std::array<float, MAX_CHAN>());
      m_pedOptVec.push_back(std::array<float, MAX_CHAN>());
      m_chi2OptVec.push_back(std::array<float, MAX_CHAN>());
 
      m_eDspVec.push_back(std::array<float, MAX_CHAN>());
      m_tDspVec.push_back(std::array<float, MAX_CHAN>());
      m_chi2DspVec.push_back(std::array<float, MAX_CHAN>());
 
      m_ROD_DMUBCIDVec.push_back(std::array<short, MAX_DMU>());
      m_ROD_DMUmemoryErrVec.push_back(std::array<short, MAX_DMU>());
      m_ROD_DMUSstrobeErrVec.push_back(std::array<short, MAX_DMU>());
      m_ROD_DMUDstrobeErrVec.push_back(std::array<short, MAX_DMU>());
      m_ROD_DMUHeadformatErrVec.push_back(std::array<short, MAX_DMU>());
      m_ROD_DMUDataformatErrVec.push_back(std::array<short, MAX_DMU>());
      m_ROD_DMUMaskVec.push_back(std::array<short, 2>());
 
      if (i % m_nDrawers < listSize) {
 
        ATH_MSG_DEBUG( "Adding items for " << suffixArr[i] );
       // create ntuple layout
        if (m_bsInput) {
          m_ntuplePtr->Branch(("Evt"+suffixArr[i]).c_str(), &m_evtVec.data()[i], ("Evt"+suffixArr[i]+"/I").c_str()); // int
          m_ntuplePtr->Branch(("rodBCID"+suffixArr[i]).c_str(), &m_rodBCIDVec.data()[i], ("rodBCID"+suffixArr[i]+"/S").c_str());//int
          m_ntuplePtr->Branch(("Size"+suffixArr[i]).c_str(), &m_sizeVec.data()[i], ("Size"+suffixArr[i]+"/S").c_str()); // short
          m_ntuplePtr->Branch(("BCID"+suffixArr[i]).c_str(), &m_bcidVec.back(), ("bcid"+suffixArr[i]+"[16]/I").c_str()); // int
          m_ntuplePtr->Branch(("DMUheader"+suffixArr[i]).c_str(), &m_DMUheaderVec.back(), ("DMUheader"+suffixArr[i]+"[16]/i").c_str()); // uint32
          m_ntuplePtr->Branch(("DMUformatErr"+suffixArr[i]).c_str(), &m_DMUformatErrVec.back(), ("DMUformatErr"+suffixArr[i]+"[16]/S").c_str()); // short
          m_ntuplePtr->Branch(("DMUparityErr"+suffixArr[i]).c_str(), &m_DMUparityErrVec.back(), ("DMUparityErr"+suffixArr[i]+"[16]/S").c_str()); // short
          m_ntuplePtr->Branch(("DMUmemoryErr"+suffixArr[i]).c_str(), &m_DMUmemoryErrVec.back(), ("DMUmemoryErr"+suffixArr[i]+"[16]/S").c_str()); // short
          m_ntuplePtr->Branch(("DMUSstrobeErr"+suffixArr[i]).c_str(), &m_DMUSstrobeErrVec.back(), ("DMUSstrobeErr"+suffixArr[i]+"[16]/S").c_str()); // short
          m_ntuplePtr->Branch(("DMUDstrobeErr"+suffixArr[i]).c_str(), &m_DMUDstrobeErrVec.back(), ("DMUDstrobeErr"+suffixArr[i]+"[16]/S").c_str()); // short
          m_ntuplePtr->Branch(("DMUMask"+suffixArr[i]).c_str(), &m_dmuMaskVec.back(), ("dmumask"+suffixArr[i]+"[2]/I").c_str()); // int
          m_ntuplePtr->Branch(("SlinkCRC"+suffixArr[i]).c_str(), &m_slinkCRCVec.back(), ("crc"+suffixArr[i]+"[2]/I").c_str()); // int
          m_ntuplePtr->Branch(("DMUHeadparityErr"+suffixArr[i]).c_str(), &m_DMUHeadparityErrVec.back(), ("DMUHeadparityErr"+suffixArr[i]+"[16]/s").c_str()); // unsigned short
          m_ntuplePtr->Branch(("DMUDataparityErr"+suffixArr[i]).c_str(), &m_DMUDataparityErrVec.back(), ("DMUDataparityErr"+suffixArr[i]+"[16]/s").c_str()); // unsigned short
        }
 
        m_ntuplePtr->Branch(("Gain"+suffixArr[i]).c_str(),&m_gainVec.back(), ("gain"+suffixArr[i]+"[48]/I").c_str()); // int
 
        if (nSamplesInDrawer > 0) {
          nSampStr = std::to_string(nSamplesInDrawer);
          m_ntuplePtr->Branch(("Sample" + suffixArr[i]).c_str(), m_sampleVec.back().get(),
                              ("sample" + suffixArr[i] + "[48]["+nSampStr+"]/I").c_str()); // size m_nsample and type int
        }
 
        if (m_bsInput) {
          m_ntuplePtr->Branch(("feCRC" + suffixArr[i]).c_str(), &m_feCRCVec.back(), ("fe_crc" + suffixArr[i] + "[16]/I").c_str()); // int
          m_ntuplePtr->Branch(("rodCRC" + suffixArr[i]).c_str(), &m_rodCRCVec.back(), ("rod_crc" + suffixArr[i] + "[16]/I").c_str()); // int
        }
 
        if (!m_flatRawChannelContainerKey.empty()) {
 
          m_ntuplePtr->Branch(("Ene" + suffixArr[i]).c_str(), &m_eneVec.back(), ("ene" + suffixArr[i] + "[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Time" + suffixArr[i]).c_str(), &m_timeVec.back(), ("time" + suffixArr[i] + "[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Ped" + suffixArr[i]).c_str(), &m_pedFlatVec.back(), ("pedflat" + suffixArr[i] + "[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Chi2ene" + suffixArr[i]).c_str(), &m_chi2FlatVec.back(), ("chiflat" + suffixArr[i] + "[48]/F").c_str()); // float
 
        }
 
        if (!m_fitRawChannelContainerKey.empty()) {
 
          m_ntuplePtr->Branch(("Efit" + suffixArr[i]).c_str(), &m_efitVec.back(), ("efit" + suffixArr[i] + "[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Tfit" + suffixArr[i]).c_str(), &m_tfitVec.back(), ("tfit" + suffixArr[i] + "[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Pedfit" + suffixArr[i]).c_str(), &m_pedfitVec.back(), ("pedfit" + suffixArr[i] + "[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Chi2fit" + suffixArr[i]).c_str(), &m_chi2Vec.back(), ("chifit" + suffixArr[i] + "[48]/F").c_str()); // float
 
        }
 
        if (!m_fitcRawChannelContainerKey.empty()) {
 
          m_ntuplePtr->Branch(("Efitc" + suffixArr[i]).c_str(), &m_efitcVec.back(), ("efitc" + suffixArr[i] + "[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Tfitc" + suffixArr[i]).c_str(), &m_tfitcVec.back(), ("tfitc" + suffixArr[i] + "[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Pedfitc" + suffixArr[i]).c_str(), &m_pedfitcVec.back(), ("pedfitc" + suffixArr[i] + "[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Chi2fitc" + suffixArr[i]).c_str(), &m_chi2cVec.back(), ("chifitc" + suffixArr[i] + "[48]/F").c_str()); // float
 
        }
 
        if (!m_optRawChannelContainerKey.empty()) {
 
          m_ntuplePtr->Branch(("Eopt"+suffixArr[i]).c_str(), &m_eOptVec.back(), ("eOpt"+suffixArr[i]+"[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Topt"+suffixArr[i]).c_str(), &m_tOptVec.back(), ("tOpt"+suffixArr[i]+"[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Pedopt"+suffixArr[i]).c_str(), &m_pedOptVec.back(), ("pedOpt"+suffixArr[i]+"[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Chi2opt"+suffixArr[i]).c_str(), &m_chi2OptVec.back(), ("chiOpt"+suffixArr[i]+"[48]/F").c_str()); // float
 
        }
 
        if (!m_dspRawChannelContainerKey.empty()) {
 
          m_ntuplePtr->Branch(("Edsp"+suffixArr[i]).c_str(), &m_eDspVec.back(), ("eDsp"+suffixArr[i]+"[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Tdsp"+suffixArr[i]).c_str(), &m_tDspVec.back(), ("tDsp"+suffixArr[i]+"[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("Chi2dsp"+suffixArr[i]).c_str(), &m_chi2DspVec.back(), ("chiDsp"+suffixArr[i]+"[48]/F").c_str()); // float
          m_ntuplePtr->Branch(("ROD_GlobalCRC"+suffixArr[i]).c_str(), &m_ROD_GlobalCRCVec.data()[i], ("ROD_GlobalCRC"+suffixArr[i]+"/s").c_str()); // unsigned short
          m_ntuplePtr->Branch(("ROD_DMUBCID"+suffixArr[i]).c_str(), &m_ROD_DMUBCIDVec.back(), ("ROD_DMUBCID"+suffixArr[i]+"[16]/s").c_str()); // unsigned short
          m_ntuplePtr->Branch(("ROD_DMUmemoryErr"+suffixArr[i]).c_str(), &m_ROD_DMUmemoryErrVec.back(), ("ROD_DMUmemoryErr"+suffixArr[i]+"[16]/s").c_str()); // unsigned short
          m_ntuplePtr->Branch(("ROD_DMUSstrobeErr"+suffixArr[i]).c_str(), &m_ROD_DMUSstrobeErrVec.back(), ("ROD_DMUSstrobeErr"+suffixArr[i]+"[16]/s").c_str()); // unsigned short
          m_ntuplePtr->Branch(("ROD_DMUDstrobeErr"+suffixArr[i]).c_str(), &m_ROD_DMUDstrobeErrVec.back(), ("ROD_DMUDstrobeErr"+suffixArr[i]+"[16]/s").c_str()); // unsigned short
          m_ntuplePtr->Branch(("ROD_DMUHeadformatErr"+suffixArr[i]).c_str(), &m_ROD_DMUHeadformatErrVec.back(), ("ROD_DMUHeadformatErr"+suffixArr[i]+"[16]/s").c_str()); // unsigned short
          m_ntuplePtr->Branch(("ROD_DMUDataformatErr"+suffixArr[i]).c_str(), &m_ROD_DMUDataformatErrVec.back(), ("ROD_DMUDataformatErr"+suffixArr[i]+"[16]/s").c_str()); // unsigned short
          m_ntuplePtr->Branch(("ROD_DMUMask"+suffixArr[i]).c_str(), &m_ROD_DMUMaskVec.back(), ("ROD_DMUMask"+suffixArr[i]+"[2]/s").c_str()); // unsigned short
 
        }
      }
    }
  }
}

DIGI_clearBranch()

void TileTBAANtuple::DIGI_clearBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// /Clear Tree DIGI variables ////////////////////////////////////////////////////////////////////////////

Definition at line 3408 of file TileTBAANtuple.cxx.

{
  clear_init_minus1(m_evtVec);
  clear_init_minus1(m_rodBCIDVec);
  clear_init_minus1(m_sizeVec);
 
  clear_init_minus1(m_bcidVec);
  clear_init_zero(m_DMUheaderVec);
  clear_init_minus1(m_DMUformatErrVec);
  clear_init_minus1(m_DMUparityErrVec);
  clear_init_minus1(m_DMUmemoryErrVec);
  clear_init_minus1(m_DMUDstrobeErrVec);
  clear_init_minus1(m_DMUSstrobeErrVec);
  clear_init_minus1(m_DMUDataparityErrVec);
  clear_init_minus1(m_DMUHeadparityErrVec);
 
  clear_init_minus1(m_ROD_GlobalCRCVec);
  clear_init_minus1(m_ROD_DMUBCIDVec);
  clear_init_minus1(m_ROD_DMUmemoryErrVec);
  clear_init_minus1(m_ROD_DMUSstrobeErrVec);
  clear_init_minus1(m_ROD_DMUDstrobeErrVec);
  clear_init_minus1(m_ROD_DMUHeadformatErrVec);
  clear_init_minus1(m_ROD_DMUDataformatErrVec);
  clear_init_minus1(m_ROD_DMUMaskVec);
 
  clear_init_minus1(m_dmuMaskVec);
  clear_init_minus1(m_slinkCRCVec);
  clear_init_minus1(m_gainVec);
  clear_samples(m_sampleVec, m_nSamplesInDrawer);
  clear_init_minus1(m_feCRCVec);
  clear_init_minus1(m_rodCRCVec);
 
  clear_init_zero(m_eneVec);
  clear_init_zero(m_timeVec);
  clear_init_zero(m_pedFlatVec);
  clear_init_zero(m_chi2FlatVec);
  clear_init_zero(m_efitVec);
  clear_init_zero(m_tfitVec);
  clear_init_zero(m_pedfitVec);
  clear_init_zero(m_chi2Vec);
  clear_init_zero(m_efitcVec);
  clear_init_zero(m_tfitcVec);
  clear_init_zero(m_pedfitcVec);
  clear_init_zero(m_chi2cVec);
  clear_init_zero(m_eOptVec);
  clear_init_zero(m_tOptVec);
  clear_init_zero(m_pedOptVec);
  clear_init_zero(m_chi2OptVec);
  clear_init_zero(m_eDspVec);
  clear_init_zero(m_tDspVec);
  clear_init_zero(m_chi2DspVec);
 
}

digiChannel2PMT()

int TileTBAANtuple::digiChannel2PMT ( int fragType, int chan )

inlineprivate

Definition at line 296 of file TileTBAANtuple.h.

                                                       {
      return (abs(m_cabling->channel2hole(fragType, chan)) - 1);
    }

ECAL_addBranch()

void TileTBAANtuple::ECAL_addBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// /Add ECAL variables to the Tree

////////////////////////////////////////////////////////////////////////////

Definition at line 2540 of file TileTBAANtuple.cxx.

{
 
  if (m_beamIdList[ECAL_ADC_FRAG]) {
    m_ntuplePtr->Branch("Ecal",&m_ecal,"m_ecal[8]/F");
  }
 
}

ECAL_clearBranch()

void TileTBAANtuple::ECAL_clearBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// Clear Tree ECAL variables

////////////////////////////////////////////////////////////////////////////

Definition at line 2606 of file TileTBAANtuple.cxx.

{
 
}

ENETOTAL_addBranch()

void TileTBAANtuple::ENETOTAL_addBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// /Add Tree ENETOTAL variables Tree

////////////////////////////////////////////////////////////////////////////

Definition at line 3037 of file TileTBAANtuple.cxx.

{
  if (m_completeNtuple) {
    m_ntuplePtr->Branch("LarSmp", &m_LarEne, "m_LarEne[4]/F");
    m_ntuplePtr->Branch("BarSmp", &m_BarEne, "m_BarEne[3]/F");
    m_ntuplePtr->Branch("ExtSmp", &m_ExtEne, "m_ExtEne[3]/F");
    m_ntuplePtr->Branch("GapSmp", &m_GapEne, "m_GapEne[3]/F");
  }
}

ENETOTAL_clearBranch()

void TileTBAANtuple::ENETOTAL_clearBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// Clear Tree ENETOTAL variables

////////////////////////////////////////////////////////////////////////////

Definition at line 3053 of file TileTBAANtuple.cxx.

{
 
 
}

evtStore()

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode TileTBAANtuple::execute ( const EventContext & ctx )

overridevirtual

Execute method.

Provides access to the EventContext if needed but is non-const as opposed to AthReentrantAlgorithm.

Implements AthAlgorithm.

Definition at line 362 of file TileTBAANtuple.cxx.

                                                          {
 
  const TileDQstatus* dqStatus = (!m_dqStatusKey.empty()) ? SG::makeHandle(m_dqStatusKey, ctx).get() : nullptr;
 
  if (m_evtNr < 0) {
 
    //bool calibMode  = (m_beamInfo->calibMode() == 1);
    //if ( calibMode != m_calibMode ) {
    //  ATH_MSG_INFO( "Calib mode from data is " << calibMode );
    //  ATH_MSG_INFO( "  Overwriting calib mode " );
    //  m_calibMode = calibMode;
    //}
 
    if (ntuple_initialize(ctx).isFailure()) {
      ATH_MSG_ERROR( "ntuple_initialize failed" );
    }
 
  }
 
  m_dspFlags = 0;
  if (ntuple_clear().isFailure()) {
    ATH_MSG_ERROR( "ntuple_clear failed" );
  }
 
  if (m_evtNr % 1000 == 0)
    ATH_MSG_INFO( m_evtNr << " events processed so far" );
 
  m_run = ctx.eventID().run_number();
  m_evt = ctx.eventID().event_number();
 
  //Get timestamp of the event
  if (ctx.eventID().time_stamp() > 0) {
    m_evTime = ctx.eventID().time_stamp();
  }
 
  // store BeamElements
  bool empty = storeBeamElements(ctx).isFailure();
 
  //store Laser Object
  empty &= storeLaser(ctx).isFailure();
 
  if (m_drawerMap.size() > 0) {
 
    // store TileDigits
    if (m_nSamples > 0) {
      empty &= (storeDigits(ctx, m_digitsContainerKey, dqStatus).isFailure());
    }
    if (m_nSamplesFlx > 0) {
      empty &= (storeDigitsFlx(ctx, m_digitsContainerFlxKey).isFailure());
    }
 
    // store TileRawChannels
    // start from DSP channels - so we can find out what is the DSP units
    empty &= (storeRawChannels(ctx, m_dspRawChannelContainerKey, m_calibMode, &m_eDspVec, &m_tDspVec, &m_chi2DspVec, nullptr, true).isFailure());
    empty &= (storeRawChannels(ctx, m_fitRawChannelContainerKey, m_calibMode, &m_efitVec, &m_tfitVec, &m_chi2Vec, &m_pedfitVec).isFailure());
    empty &= (storeRawChannels(ctx, m_flxFitRawChannelContainerKey, true, &m_eflxfitVec, &m_tflxfitVec, &m_chi2flxfitVec, &m_pedflxfitVec).isFailure());
    empty &= (storeRawChannels(ctx, m_optRawChannelContainerKey, m_calibMode, &m_eOptVec, &m_tOptVec, &m_chi2OptVec, &m_pedOptVec).isFailure());
    empty &= (storeRawChannels(ctx, m_flxOptRawChannelContainerKey, true, &m_eflxoptVec, &m_tflxoptVec, &m_chi2flxoptVec, &m_pedflxoptVec).isFailure());
    empty &= (storeRawChannels(ctx, m_fitcRawChannelContainerKey, m_calibMode, &m_efitcVec, &m_tfitcVec, &m_chi2cVec, &m_pedfitcVec).isFailure());
    empty &= (storeRawChannels(ctx, m_flatRawChannelContainerKey, m_calibMode, &m_eneVec, &m_timeVec, &m_chi2FlatVec, &m_pedFlatVec).isFailure());
 
    empty &= (storeHitVector(ctx).isFailure());
    empty &= (storeHitContainer(ctx).isFailure());
  }
 
 
  if (m_completeNtuple && m_TBperiod < 2015) {
    // store energy per sampling from all calorimeters
    empty &= (storeCells(ctx).isFailure());
  }
 
  // increase event nr
  // this number can be different from real event number if we skip events
 
  if (empty) {
    ATH_MSG_WARNING( "Error in execute " );
  }
 
  m_ntuplePtr->Fill();
 
  ++m_evtNr;
 
  // Execution completed.
  ATH_MSG_DEBUG( "execute() completed successfully" );
 
  return StatusCode::SUCCESS;
}

extraDeps_update_handler()

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

protectedinherited

Add StoreName to extra input/output deps as needed.

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

extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 50 of file AthAlgorithm.cxx.

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

FELIX_addBranch()

void TileTBAANtuple::FELIX_addBranch ( void )

private

Definition at line 3463 of file TileTBAANtuple.cxx.

{
 
  m_rodBCIDflxVec.reserve(MAX_DRAWERS);
  m_sizeflxVec.reserve(MAX_DRAWERS);
  m_evtflxVec.reserve(MAX_DRAWERS);
 
  m_eflxfitVec.reserve(MAX_DRAWERS);
  m_tflxfitVec.reserve(MAX_DRAWERS);
  m_chi2flxfitVec.reserve(MAX_DRAWERS);
  m_pedflxfitVec.reserve(MAX_DRAWERS);
  m_eflxoptVec.reserve(MAX_DRAWERS);
  m_tflxoptVec.reserve(MAX_DRAWERS);
  m_chi2flxoptVec.reserve(MAX_DRAWERS);
  m_pedflxoptVec.reserve(MAX_DRAWERS);
  m_gainflxVec.reserve(MAX_DRAWERS);
  m_sampleflxVec.reserve(MAX_DRAWERS);
 
  m_mdL1idflxVec.reserve(MAX_DRAWERS);
  m_mdBcidflxVec.reserve(MAX_DRAWERS);
  m_mdModuleflxVec.reserve(MAX_DRAWERS);
  m_mdRunTypeflxVec.reserve(MAX_DRAWERS);
  m_mdPedLoflxVec.reserve(MAX_DRAWERS);
  m_mdPedHiflxVec.reserve(MAX_DRAWERS);
  m_mdRunflxVec.reserve(MAX_DRAWERS);
  m_mdChargeflxVec.reserve(MAX_DRAWERS);
  m_mdChargeTimeflxVec.reserve(MAX_DRAWERS);
  m_mdCapacitorflxVec.reserve(MAX_DRAWERS);
 
  std::string nSampStrFlx = std::to_string(m_nSamplesFlx);
  unsigned int listSize = std::min(m_nDrawersFlx.value(), static_cast<unsigned int>(m_drawerFlxMap.size()));
 
  bool testbeam = TileCablingService::getInstance()->getTestBeam();
  std::vector<std::string> moduleNames(m_nDrawersFlx, "");
  unsigned int length =  2 * m_nDrawersFlx;
  std::vector<std::string> suffixArr(length, "");
  m_nSamplesFlxInDrawer.resize(length);
 
  for (const std::pair<const unsigned int, unsigned int>& fragAndDrawer : m_drawerFlxMap) {
    unsigned int frag = fragAndDrawer.first;
    unsigned int ros = frag >> 8;
    unsigned int drawer = frag & 0x3F;
    unsigned int drawerIndex = fragAndDrawer.second;
 
    std::ostringstream drawerName;
    drawerName << m_rosName[ros];
    if (testbeam) {
      drawerName << (drawer & 7);
    } else {
      ++drawer; // count modules from 1
      drawerName << std::setw(2) << std::setfill('0') << drawer;
    }
 
    moduleNames.at(drawerIndex) = drawerName.str();
    suffixArr.at(drawerIndex) = drawerName.str() + "lo";
    suffixArr.at(drawerIndex + m_nDrawersFlx) = drawerName.str() + "hi";
 
    int nSamples(m_nSamplesFlx);
    auto it = m_nSamplesFlxInDrawerMap.find(frag);
    if (it != m_nSamplesFlxInDrawerMap.end()) {
      nSamples = it->second;
    }
 
    m_nSamplesFlxInDrawer[drawerIndex] = nSamples;
    m_nSamplesFlxInDrawer[drawerIndex + m_nDrawersFlx] = nSamples;
  }
 
  m_rodBCIDflxVec.resize(m_nDrawersFlx);
  m_sizeflxVec.resize(m_nDrawersFlx);
  m_evtflxVec.resize(m_nDrawersFlx);
 
  for (unsigned int i = 0; i < length; ++i) {
 
    int nSamplesInDrawer = m_nSamplesFlxInDrawer[i];
 
    m_eflxfitVec.push_back(std::array<float, MAX_CHAN>());
    m_tflxfitVec.push_back(std::array<float, MAX_CHAN>());
    m_pedflxfitVec.push_back(std::array<float, MAX_CHAN>());
    m_chi2flxfitVec.push_back(std::array<float, MAX_CHAN>());
 
    m_eflxoptVec.push_back(std::array<float, MAX_CHAN>());
    m_tflxoptVec.push_back(std::array<float, MAX_CHAN>());
    m_pedflxoptVec.push_back(std::array<float, MAX_CHAN>());
    m_chi2flxoptVec.push_back(std::array<float, MAX_CHAN>());
 
    if (i < m_nDrawersFlx) {
      m_mdL1idflxVec.push_back(std::array<int, MAX_MINIDRAWER>());
      m_mdBcidflxVec.push_back(std::array<int, MAX_MINIDRAWER>());
      m_mdModuleflxVec.push_back(std::array<int, MAX_MINIDRAWER>());
      m_mdRunTypeflxVec.push_back(std::array<int, MAX_MINIDRAWER>());
      m_mdPedLoflxVec.push_back(std::array<int, MAX_MINIDRAWER>());
      m_mdPedHiflxVec.push_back(std::array<int, MAX_MINIDRAWER>());
      m_mdRunflxVec.push_back(std::array<int, MAX_MINIDRAWER>());
      m_mdChargeflxVec.push_back(std::array<int, MAX_MINIDRAWER>());
      m_mdChargeTimeflxVec.push_back(std::array<int, MAX_MINIDRAWER>());
      m_mdCapacitorflxVec.push_back(std::array<int, MAX_MINIDRAWER>());
 
      std::string suffix = moduleNames[i];
      ATH_MSG_DEBUG( "Adding items for " << suffix );
 
      if (m_bsInput) {
        m_ntuplePtr->Branch(("FlxEvt"+suffix).c_str(), &m_evtflxVec.data()[i]); // int
        m_ntuplePtr->Branch(("FlxRodBCID"+suffix).c_str(), &m_rodBCIDflxVec.data()[i]);//int
        m_ntuplePtr->Branch(("FlxSize"+suffix).c_str(), &m_sizeflxVec.data()[i]); // short
      }
 
      m_ntuplePtr->Branch(("FlxL1ID"+suffix).c_str(), &m_mdL1idflxVec.back()); // int
      m_ntuplePtr->Branch(("FlxBCID"+suffix).c_str(), &m_mdBcidflxVec.back()); // int
      m_ntuplePtr->Branch(("FlxModule"+suffix).c_str(), &m_mdModuleflxVec.back()); // int
      m_ntuplePtr->Branch(("FlxRun"+suffix).c_str(), &m_mdRunflxVec.back()); // int
      m_ntuplePtr->Branch(("FlxRunType"+suffix).c_str(), &m_mdRunTypeflxVec.back()); // int
      m_ntuplePtr->Branch(("FlxPedLo"+suffix).c_str(), &m_mdPedLoflxVec.back()); // int
      m_ntuplePtr->Branch(("FlxPedHi"+suffix).c_str(), &m_mdPedHiflxVec.back()); // int
      m_ntuplePtr->Branch(("FlxCharge"+suffix).c_str(), &m_mdChargeflxVec.back()); // int
      m_ntuplePtr->Branch(("FlxChargeTime"+suffix).c_str(), &m_mdChargeTimeflxVec.back()); // int
      m_ntuplePtr->Branch(("FlxCapacitor"+suffix).c_str(), &m_mdCapacitorflxVec.back()); // int
 
    }
 
    m_gainflxVec.push_back(std::array<int, MAX_CHAN>());
    m_sampleflxVec.push_back(std::make_unique<int[]>(MAX_CHAN * nSamplesInDrawer)); // U(48/96,9)
 
    if (i % m_nDrawersFlx < listSize) {
      ATH_MSG_DEBUG( "Adding items for " << suffixArr[i] );
 
      m_ntuplePtr->Branch(("FlxGain"+suffixArr[i]).c_str(), &m_gainflxVec.back()); // int
 
      if (nSamplesInDrawer > 0) {
        nSampStrFlx = std::to_string(nSamplesInDrawer);
        m_ntuplePtr->Branch(("FlxSample" + suffixArr[i]).c_str(), m_sampleflxVec.back().get(),
                            ("Flxsample" + suffixArr[i] + "[48]["+nSampStrFlx+"]/I").c_str()); // size m_nsample and type int
      }
 
      if (!m_flxFitRawChannelContainerKey.empty()) {
        m_ntuplePtr->Branch(("FlxEfit" + suffixArr[i]).c_str(), &m_eflxfitVec.back()); // float
        m_ntuplePtr->Branch(("FlxTfit" + suffixArr[i]).c_str(), &m_tflxfitVec.back()); // float
        m_ntuplePtr->Branch(("FlxPedfit" + suffixArr[i]).c_str(), &m_pedflxfitVec.back()); // float
        m_ntuplePtr->Branch(("FlxChi2fit" + suffixArr[i]).c_str(), &m_chi2flxfitVec.back()); // float
      }
 
      if (!m_flxOptRawChannelContainerKey.empty()) {
        m_ntuplePtr->Branch(("FlxEOpt" + suffixArr[i]).c_str(), &m_eflxoptVec.back()); // float
        m_ntuplePtr->Branch(("FlxTOpt" + suffixArr[i]).c_str(), &m_tflxoptVec.back()); // float
        m_ntuplePtr->Branch(("FlxPedOpt" + suffixArr[i]).c_str(), &m_pedflxoptVec.back()); // float
        m_ntuplePtr->Branch(("FlxChi2Opt" + suffixArr[i]).c_str(), &m_chi2flxoptVec.back()); // float
      }
    }
  }
}

FELIX_clearBranch()

void TileTBAANtuple::FELIX_clearBranch ( void )

private

Definition at line 3615 of file TileTBAANtuple.cxx.

{
  clear_init_minus1(m_evtflxVec);
  clear_init_minus1(m_rodBCIDflxVec);
  clear_init_minus1(m_sizeflxVec);
 
  clear_init_minus1(m_gainflxVec);
  clear_samples(m_sampleflxVec, m_nSamplesFlxInDrawer);
 
  clear_init_zero(m_eflxfitVec);
  clear_init_zero(m_tflxfitVec );
  clear_init_zero(m_chi2flxfitVec);
  clear_init_zero(m_pedflxfitVec);
 
  clear_init_zero(m_eflxoptVec);
  clear_init_zero(m_tflxoptVec);
  clear_init_zero(m_chi2flxoptVec);
  clear_init_zero(m_pedflxoptVec);
 
  clear_init_minus1(m_mdL1idflxVec);
  clear_init_minus1(m_mdBcidflxVec);
  clear_init_minus1(m_mdModuleflxVec);
  clear_init_minus1(m_mdRunTypeflxVec);
  clear_init_minus1(m_mdRunflxVec);
  clear_init_minus1(m_mdChargeflxVec);
  clear_init_minus1(m_mdChargeTimeflxVec);
  clear_init_minus1(m_mdCapacitorflxVec);
  clear_init_minus1(m_mdPedLoflxVec);
  clear_init_minus1(m_mdPedHiflxVec);
}

filterPassed() [1/2]

bool AthAlgorithm::filterPassed ( ) const

inherited

Definition at line 94 of file AthAlgorithm.cxx.

                                      {
  return filterPassed( Gaudi::Hive::currentContext() );
}

filterPassed() [2/2]

bool AthAlgorithm::filterPassed ( const EventContext & ctx ) const

inherited

Definition at line 98 of file AthAlgorithm.cxx.

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

getContext()

const EventContext & AthAlgorithm::getContext ( ) const

inherited

Deprecated methods (use the ones with EventContext).

Definition at line 90 of file AthAlgorithm.cxx.

                                                   {
  return Gaudi::Hive::currentContext();
}

getEta()

void TileTBAANtuple::getEta ( void )

private

Definition at line 2396 of file TileTBAANtuple.cxx.

                            {
 
  // Get eta from an ASCII file with the following structure :
  // runnumber eta
 
  ATH_MSG_INFO( "Get eta for run " << m_runNumber );
 
  // Find the full path to filename:
  std::string fileName = PathResolver::find_file(m_etaFileName, "DATAPATH");
  ATH_MSG_INFO( "Reading file  " << fileName );
 
  if (fileName.size() == 0) {
 
    ATH_MSG_WARNING( "Could not find input file " << m_etaFileName );
    ATH_MSG_WARNING( "Skip reading of eta value " );
 
  } else {
 
    std::ifstream etafile;
    etafile.open(fileName.c_str());
 
    if (etafile.good()) {
 
      int runNumber = 0;
      float eta = 0;
      //coverity[TAINTED_SCALAR]
      while ((runNumber != m_runNumber) && (!etafile.eof())) {
        etafile >> runNumber >> eta;
      }
 
      if (runNumber != m_runNumber) {
        ATH_MSG_INFO( "Run " << m_runNumber << " has not been found, keep eta and theta at zero" );
 
        m_eta = m_theta = 0.0;
      } else {
        m_eta = eta;
        m_theta = (M_PI_2 - 2 * atan(exp(m_eta)));
        ATH_MSG_INFO( "Run " << m_runNumber << " has been found with eta=" << m_eta << ", theta =" << m_theta );
 
      }
 
    } else {
 
      ATH_MSG_WARNING( "Problem with file " << fileName );
      ATH_MSG_WARNING( "Skip reading of eta value " );
 
    }
 
    etafile.close();
  }
}

HIT_addBranch()

void TileTBAANtuple::HIT_addBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// /Add Tree HIT variables Tree

////////////////////////////////////////////////////////////////////////////

Definition at line 3652 of file TileTBAANtuple.cxx.

{
  m_ehitVec.clear();
  m_thitVec.clear();
  m_ehitCnt.clear();
  m_thitCnt.clear();
 
  m_ehitVec.reserve(MAX_DRAWERS);
  m_thitVec.reserve(MAX_DRAWERS);
  m_ehitCnt.reserve(MAX_DRAWERS);
  m_thitCnt.reserve(MAX_DRAWERS);
 
  unsigned int listSize = std::min(m_nDrawers.value(), static_cast<unsigned int>(m_drawerMap.size()));
 
  if (listSize > 0) {
 
    std::string digit[10] = { "0", "1", "2", "3", "4", "5", "6", "7", "8", "9" };
    std::vector<std::string> suffixArr;
    unsigned int length;
    bool testbeam = TileCablingService::getInstance()->getTestBeam();
 
    length = m_nDrawers;
    suffixArr.resize(length);
 
    for (unsigned int i = 0; i < listSize; ++i) {
      unsigned int ros = m_drawerType[i];
      unsigned int drawer = strtol(m_drawerList[i].data(), NULL, 0) & 0x3F;
      std::string digits;
      if (m_TBperiod >= 2010) {
        ++drawer; // count modules from 1
        digits = digit[drawer / 10] + digit[drawer % 10];
      } else if (testbeam) {
        digits = digit[drawer & 7];
      } else {
        ++drawer; // count modules from 1
        digits = digit[drawer / 10] + digit[drawer % 10];
      }
 
      if (ros == 0) {
        std::string suff = m_drawerList[i];
        suff.replace(suff.find("0x"), 2, "");
        suffixArr[i] = std::move(suff);
      } else {
        suffixArr[i] = m_rosName[ros] + digits;
      }
    }
 
    for (unsigned int i = 0; i < length; i++) {
      
      if (i % m_nDrawers < listSize)
        ATH_MSG_DEBUG( "Adding items for " << suffixArr[i] );
      
      m_ehitVec.push_back(std::array<float, MAX_CHAN>());
      m_thitVec.push_back(std::array<float, MAX_CHAN>());
      m_ehitCnt.push_back(std::array<float, MAX_CHAN>());
      m_thitCnt.push_back(std::array<float, MAX_CHAN>());
 
      if (i%m_nDrawers < listSize) {
 
        if (!m_hitVectorKey.empty()) {
 
          if (i % m_nDrawers < listSize)
            ATH_MSG_DEBUG( "Adding G4 hit info for " << suffixArr[i] );
      
          m_ntuplePtr->Branch(("EhitG4"+suffixArr[i]).c_str(),&m_ehitVec.back(),("eHitG4"+suffixArr[i]+"[48]/F").c_str());
          m_ntuplePtr->Branch(("ThitG4"+suffixArr[i]).c_str(),&m_thitVec.back(),("tHitG4"+suffixArr[i]+"[48]/F").c_str());
        }
 
        if (!m_hitContainerKey.empty()) {
 
          if (i % m_nDrawers < listSize)
            ATH_MSG_DEBUG( "Adding G4 corrected hit info for " << suffixArr[i] );
      
          m_ntuplePtr->Branch(("EhitSim"+suffixArr[i]).c_str(),&m_ehitCnt.back(),("eHitSim"+suffixArr[i]+"[48]/F").c_str());
          m_ntuplePtr->Branch(("ThitSim"+suffixArr[i]).c_str(),&m_thitCnt.back(),("tHitSim"+suffixArr[i]+"[48]/F").c_str());
        }
      }
    }
  }
}

HIT_clearBranch()

void TileTBAANtuple::HIT_clearBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// /Clear Tree HIT variables ////////////////////////////////////////////////////////////////////////////

Definition at line 3738 of file TileTBAANtuple.cxx.

{
  clear_init_zero(m_ehitVec);
  clear_init_zero(m_thitVec);
  clear_init_zero(m_ehitCnt);
  clear_init_zero(m_thitCnt);
}

initialize()

StatusCode TileTBAANtuple::initialize ( )

overridevirtual

Definition at line 102 of file TileTBAANtuple.cxx.

                                      {
  ATH_CHECK( m_samplingFractionKey.initialize(SG::AllowEmpty) );
  m_saveFelixData = !(m_digitsContainerFlxKey.empty() && m_flxOptRawChannelContainerKey.empty() && m_flxFitRawChannelContainerKey.empty());
 
  ATH_CHECK( m_digitsContainerKey.initialize(SG::AllowEmpty) );
  ATH_CHECK( m_digitsContainerFlxKey.initialize(SG::AllowEmpty) );
  ATH_CHECK( m_beamElemContainerKey.initialize(SG::AllowEmpty) );
  ATH_CHECK( m_flatRawChannelContainerKey.initialize(SG::AllowEmpty) );
  ATH_CHECK( m_fitRawChannelContainerKey.initialize(SG::AllowEmpty) );
  ATH_CHECK( m_optRawChannelContainerKey.initialize(SG::AllowEmpty) );
  ATH_CHECK( m_dspRawChannelContainerKey.initialize(SG::AllowEmpty) );
  ATH_CHECK( m_fitcRawChannelContainerKey.initialize(SG::AllowEmpty) );
  ATH_CHECK( m_flxFitRawChannelContainerKey.initialize(SG::AllowEmpty) );
  ATH_CHECK( m_flxOptRawChannelContainerKey.initialize(SG::AllowEmpty) );
  ATH_CHECK( m_laserObjectKey.initialize(SG::AllowEmpty) );
  ATH_CHECK( m_hitContainerKey.initialize(SG::AllowEmpty) );
  ATH_CHECK( m_hitVectorKey.initialize(SG::AllowEmpty) );
  ATH_CHECK( m_cellContainerKey.initialize(!m_cellContainerKey.empty() && m_completeNtuple && m_TBperiod < 2015) );
 
  ATH_CHECK( m_tileToolEmscale.retrieve() );
  ATH_CHECK( m_adderFilterAlgTool.retrieve(EnableTool{m_unpackAdder}) );
 
  ATH_CHECK( m_dqStatusKey.initialize(SG::AllowEmpty) );
 
  return StatusCode::SUCCESS;
}

initList()

StatusCode TileTBAANtuple::initList ( const EventContext & ctx )

private

Definition at line 2188 of file TileTBAANtuple.cxx.

                                                           {
 
  unsigned int size = m_drawerList.size();
 
  if (size > 0) {
 
    if (m_digitsContainerKey.empty()) { // empty name, nothing to do
 
      ATH_MSG_WARNING( "can't retrieve Digits from TDS" );
      ATH_MSG_WARNING( "can't set up fragment list for ntuple" );
 
      if (m_nSamples != 0) {
        ATH_MSG_WARNING( "Disable digit samples in ntuple" );
        m_nSamples = 0;
      }
 
      return StatusCode::SUCCESS;
    }
 
    // Read Digits from TES
    const TileDigitsContainer* digitsCnt = SG::makeHandle (m_digitsContainerKey, ctx).get();
 
    int frag = strtol(m_drawerList[0].data(), NULL, 0);
    if (frag < 0) { // setup frags IDs from the data
 
      std::vector<unsigned int> frags;
      // Go through all TileDigitsCollections
      for (const TileDigitsCollection* digitsCollection : *digitsCnt) {
        if (!digitsCollection->empty()) {
          // determine type of frag
          frags.push_back(digitsCollection->identify());
        }
      }
      size = frags.size();
 
      if (size > 0) {
 
        if (size < m_nDrawers) {
          ATH_MSG_INFO( "decreasing m_nDrawers from " << m_nDrawers << " to " << size );
          m_nDrawers = size;
        }
 
        unsigned int rosOrder[5] = { 2, 1, 3, 4, 0 };
        unsigned int dr = 0;
        char frg[6] = "0x000";
 
        m_drawerList.clear();
        // m_drawerType.clear();
        m_drawerMap.clear();
 
        msg(MSG::INFO) << "setting drawerList from data " << MSG::hex;
        for (unsigned int ir = 0; ir < 5; ++ir) {
          for (unsigned int i = 0; i < size; ++i) {
            unsigned int frag = frags[i];
            if (frag >> 8 == rosOrder[ir]) {
              sprintf(frg, "0x%3.3x", frag);
              m_drawerList.value().push_back((std::string) frg);
              if (dr == m_drawerType.size()) m_drawerType.value().push_back(frag >> 8);
              m_drawerMap[frag] = dr;
              msg(MSG::INFO) << " 0x" << frag;
              ++dr;
            }
          }
        }
 
        msg(MSG::INFO) << MSG::dec << endmsg;
 
        size = m_drawerType.size();
        if (size < m_nDrawers) {
          m_drawerType.value().resize(m_nDrawers);
          for (; size < m_nDrawers; ++size)
            m_drawerType[size] = 0;
        }
 
        msg(MSG::INFO) << MSG::INFO << "drawerType ";
        for (unsigned int dr = 0; dr < size; ++dr)
          msg(MSG::INFO) << " " << m_drawerType[dr];
        msg(MSG::INFO) << endmsg;
 
        if (size > m_nDrawers) {
          ATH_MSG_INFO( "increasing m_nDrawers from " << m_nDrawers << " to " << size );
          m_nDrawers = size;
        }
 
        if (size < 1) size = 1;
        if (m_eventsPerFile == 0) {
          m_eventsPerFile = static_cast<int>(200 / size) * 1000;
          ATH_MSG_INFO( "Number of events per file was 0, set it to 200k/" << size << " = " << m_eventsPerFile );
        }
 
      } else {
 
        ATH_MSG_ERROR( "can't find any TileDigits collections" );
        ATH_MSG_ERROR( "can't set up fragment list for ntuple" );
      }
    }
 
    // once again - check number of samples in the data
    // but do not print any ERRORs now
 
    // Go through all TileDigitsCollections
    for (const TileDigitsCollection* digitsCollection : *digitsCnt) {
      if (!digitsCollection->empty()) {
        int siz = digitsCollection->front()->samples().size();
        m_nSamplesInDrawerMap[digitsCollection->identify()] = siz;
        if (siz > m_nSamples && m_nSamples != 0) {
          ATH_MSG_WARNING( "Increasing number of digit samples in ntuple from " << m_nSamples << " to " << siz );
          m_nSamples = siz;
        }
      }
    }
  }
 
  return StatusCode::SUCCESS;
}

initListFlx()

StatusCode TileTBAANtuple::initListFlx ( const EventContext & ctx )

private

Definition at line 2305 of file TileTBAANtuple.cxx.

                                                              {
 
  unsigned int listSize = m_drawerList.size();
 
  if (listSize > 0) {
 
    if (m_digitsContainerFlxKey.empty()) { // empty name, nothing to do
 
      ATH_MSG_WARNING( "can't retrieve FELIX Digits from TDS" );
      ATH_MSG_WARNING( "can't set up FELIX fragment list for ntuple" );
 
      if (m_nSamples != 0) {
        ATH_MSG_WARNING( "Disable FELIX digit samples in ntuple" );
        m_nSamplesFlx = 0;
      }
 
      return StatusCode::SUCCESS;
    }
 
    // Read Digits from TES
    const TileDigitsContainer* digitsCntFlx = SG::makeHandle (m_digitsContainerFlxKey, ctx).get();
 
    if (listSize == m_nDrawers) {
      // Only legacy drawers in the list, setup FELIX frags IDs from the data
 
      std::vector<unsigned int> frags;
      // Go through all TileDigitsCollections
      for (const TileDigitsCollection* digitsCollection : *digitsCntFlx) {
        if (!digitsCollection->empty()) {
          // Determine type of frag
          frags.push_back(digitsCollection->identify());
        }
      }
 
      unsigned int nFrags = frags.size();
 
      if (nFrags > 0) {
 
        if (nFrags != m_nDrawersFlx) {
          ATH_MSG_INFO( "changing m_nDrawersFlx from " << m_nDrawersFlx.value() << " to " << nFrags );
          m_nDrawersFlx = nFrags;
        }
 
        m_drawerFlxMap.clear();
 
        std::ostringstream os;
        os << "setting FELIX drawers from data " << std::hex;
        unsigned int drawerIndex = 0;
        for (unsigned int frag : frags) {
          m_drawerFlxMap[frag] = drawerIndex;
          os << " 0x" << frag;
          ++drawerIndex;
        }
        os << std::dec;
 
        ATH_MSG_INFO(os.str());
 
        if (m_eventsPerFile == 0) {
          int nDrawersAll = m_nDrawers.value() + m_nDrawersFlx.value();
          m_eventsPerFile = static_cast<int>(200 / nDrawersAll) * 1000;
          ATH_MSG_INFO( "Number of events per file was 0, set it to 200k/" << nDrawersAll << " = " << m_eventsPerFile );
        }
 
      } else {
        ATH_MSG_ERROR( "can't find any FELIX TileDigits collections" );
        ATH_MSG_ERROR( "can't set up FELIX fragment list for ntuple" );
      }
    }
 
 
    // once again - check number of samples in the data
    // but do not print any ERRORs now
 
    for (const TileDigitsCollection* digitsCollection : *digitsCntFlx) {
      if (!digitsCollection->empty()) {
        int siz = digitsCollection->front()->samples().size();
        m_nSamplesFlxInDrawerMap[digitsCollection->identify()] = siz;
        if (siz > m_nSamplesFlx && m_nSamplesFlx != 0) {
          ATH_MSG_WARNING( "Increasing number of FELIX  digit samples in ntuple from " << m_nSamplesFlx << " to " << siz );
          m_nSamplesFlx = siz;
        }
      }
    }
  }
 
  return StatusCode::SUCCESS;
}

initNTuple()

StatusCode TileTBAANtuple::initNTuple ( void )

private

Definition at line 2019 of file TileTBAANtuple.cxx.

                                          {
  MsgStream log(msgSvc(), name());
 
  m_evtVec.clear();
  m_bcidVec.clear();
  m_DMUheaderVec.clear();
  m_DMUformatErrVec.clear();
  m_DMUparityErrVec.clear();
  m_DMUmemoryErrVec.clear();
  m_DMUDstrobeErrVec.clear();
  m_DMUSstrobeErrVec.clear();
  m_DMUHeadparityErrVec.clear();
  m_DMUDataparityErrVec.clear();
  m_rodBCIDVec.clear();
  m_sizeVec.clear();
  m_dmuMaskVec.clear();
  m_slinkCRCVec.clear();
  m_gainVec.clear();
  m_sampleVec.clear();
  m_feCRCVec.clear();
  m_rodCRCVec.clear();
 
  m_evtflxVec.clear();
  m_rodBCIDflxVec.clear();
  m_sizeflxVec.clear();
  m_gainflxVec.clear();
  m_sampleflxVec.clear();
 
  m_eneVec.clear();
  m_timeVec.clear();
  m_pedFlatVec.clear();
  m_chi2FlatVec.clear();
 
  m_efitVec.clear();
  m_tfitVec.clear();
  m_pedfitVec.clear();
  m_chi2Vec.clear();
 
  m_efitcVec.clear();
  m_tfitcVec.clear();
  m_pedfitcVec.clear();
  m_chi2cVec.clear();
 
  m_eOptVec.clear();
  m_tOptVec.clear();
  m_pedOptVec.clear();
  m_chi2OptVec.clear();
 
  m_eflxfitVec.clear();
  m_tflxfitVec.clear();
  m_pedflxfitVec.clear();
  m_chi2flxfitVec.clear();
 
  m_eflxoptVec.clear();
  m_tflxoptVec.clear();
  m_pedflxoptVec.clear();
  m_chi2flxoptVec.clear();
 
  m_eDspVec.clear();
  m_tDspVec.clear();
  m_chi2DspVec.clear();
 
  m_ROD_GlobalCRCVec.clear();
  m_ROD_DMUBCIDVec.clear();
  m_ROD_DMUmemoryErrVec.clear();
  m_ROD_DMUSstrobeErrVec.clear();
  m_ROD_DMUDstrobeErrVec.clear();
  m_ROD_DMUHeadformatErrVec.clear();
  m_ROD_DMUDataformatErrVec.clear();
  m_ROD_DMUMaskVec.clear();
 
  //Ntuple creation
 
  auto tree = std::make_unique<TTree>(m_ntupleID.value().c_str(), "TileBEAM-Ntuple");
  tree->SetMaxTreeSize(m_treeSize);
  m_ntuplePtr = tree.get();
  if (m_thistSvc->regTree("/" + m_streamName + "/" + m_ntupleID, std::move(tree)).isFailure()) {
    ATH_MSG_ERROR( "Problem registering TileRec Tree" );
    m_ntupleCreated = false;
  } else {
    m_ntupleCreated = true;
  }
 
  TRIGGER_addBranch();
  MUON_addBranch();
 
  if (m_TBperiod < 2015) {
    ECAL_addBranch();
    LASER_addBranch();
    ADDER_addBranch();
    ENETOTAL_addBranch();
    COINCBOARD_addBranch();
  }
 
  if (m_TBperiod > 2015) {
    QDC_addBranch();
  }
 
  CISPAR_addBranch();
  BEAM_addBranch();
  DIGI_addBranch(); //working now
  if (m_saveFelixData) {
    FELIX_addBranch();
  }
 
  HIT_addBranch();
 
  return StatusCode::SUCCESS;
}

inputHandles()

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

overridevirtualinherited

Return this algorithm's input handles.

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

isReEntrant()

virtual bool AthAlgorithm::isReEntrant ( ) const

inlinefinaloverrideprotectedvirtualinherited

Legacy algorithms are not thread-safe.

Definition at line 111 of file AthAlgorithm.h.

{ return false; }

LASER_addBranch()

void TileTBAANtuple::LASER_addBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// /Add Tree LASER variables Tree

////////////////////////////////////////////////////////////////////////////

Definition at line 2629 of file TileTBAANtuple.cxx.

                                         {
 
  if (!m_laserObjectKey.empty() > 0) {
 
    m_ntuplePtr->Branch("LASER_BCID", &m_las_BCID, "LASER_BCID/I");
 
    m_ntuplePtr->Branch("LASER_FILTER", &m_las_Filt, "LASER_FILTER/I");
    m_ntuplePtr->Branch("LASER_REQAMP", &m_las_ReqAmp, "LASER_REQAMP/D");
    m_ntuplePtr->Branch("LASER_MEASAMP", &m_las_MeasAmp, "LASER_MEASAMP/D");
 
    m_ntuplePtr->Branch("LASER_Diode_1_ADC", &m_las_D1_ADC, "LASER_Diode_1_ADC/I");
    m_ntuplePtr->Branch("LASER_Diode_2_ADC", &m_las_D2_ADC, "LASER_Diode_2_ADC/I");
    m_ntuplePtr->Branch("LASER_Diode_3_ADC", &m_las_D3_ADC, "LASER_Diode_3_ADC/I");
    m_ntuplePtr->Branch("LASER_Diode_4_ADC", &m_las_D4_ADC, "LASER_Diode_4_ADC/I");
 
    m_ntuplePtr->Branch("LASER_Diode_1_Ped", &m_las_D1_Ped, "LASER_Diode_1_Ped/D");
    m_ntuplePtr->Branch("LASER_Diode_2_Ped", &m_las_D2_Ped, "LASER_Diode_2_Ped/D");
    m_ntuplePtr->Branch("LASER_Diode_3_Ped", &m_las_D3_Ped, "LASER_Diode_3_Ped/D");
    m_ntuplePtr->Branch("LASER_Diode_4_Ped", &m_las_D4_Ped, "LASER_Diode_4_Ped/D");
 
    m_ntuplePtr->Branch("LASER_Diode_1_Ped_RMS", &m_las_D1_Ped_RMS, "LASER_Diode_1_Ped_RMS/D");
    m_ntuplePtr->Branch("LASER_Diode_2_Ped_RMS", &m_las_D2_Ped_RMS, "LASER_Diode_2_Ped_RMS/D");
    m_ntuplePtr->Branch("LASER_Diode_3_Ped_RMS", &m_las_D1_Ped_RMS, "LASER_Diode_3_Ped_RMS/D");
    m_ntuplePtr->Branch("LASER_Diode_4_Ped_RMS", &m_las_D1_Ped_RMS, "LASER_Diode_4_Ped_RMS/D");
 
    m_ntuplePtr->Branch("LASER_Diode_1_Alpha", &m_las_D1_Alpha, "LASER_Diode_1_Alpha/D");
    m_ntuplePtr->Branch("LASER_Diode_2_Alpha", &m_las_D2_Alpha, "LASER_Diode_2_Alpha/D");
    m_ntuplePtr->Branch("LASER_Diode_3_Alpha", &m_las_D3_Alpha, "LASER_Diode_3_Alpha/D");
    m_ntuplePtr->Branch("LASER_Diode_4_Alpha", &m_las_D4_Alpha, "LASER_Diode_4_Alpha/D");
 
    m_ntuplePtr->Branch("LASER_Diode_1_Alpha_RMS", &m_las_D1_Alpha_RMS, "LASER_Diode_1_Alpha_RMS/D");
    m_ntuplePtr->Branch("LASER_Diode_2_Alpha_RMS", &m_las_D2_Alpha_RMS, "LASER_Diode_2_Alpha_RMS/D");
    m_ntuplePtr->Branch("LASER_Diode_3_Alpha_RMS", &m_las_D3_Alpha_RMS, "LASER_Diode_3_Alpha_RMS/D");
    m_ntuplePtr->Branch("LASER_Diode_4_Alpha_RMS", &m_las_D4_Alpha_RMS, "LASER_Diode_4_Alpha_RMS/D");
 
    m_ntuplePtr->Branch("LASER_Diode_1_AlphaPed", &m_las_D1_AlphaPed, "LASER_Diode_1_AlphaPed/D");
    m_ntuplePtr->Branch("LASER_Diode_2_AlphaPed", &m_las_D2_AlphaPed, "LASER_Diode_2_AlphaPed/D");
    m_ntuplePtr->Branch("LASER_Diode_3_AlphaPed", &m_las_D3_AlphaPed, "LASER_Diode_3_AlphaPed/D");
    m_ntuplePtr->Branch("LASER_Diode_4_AlphaPed", &m_las_D4_AlphaPed, "LASER_Diode_4_AlphaPed/D");
 
    m_ntuplePtr->Branch("LASER_Diode_1_AlphaPed_RMS", &m_las_D1_AlphaPed_RMS, "LASER_Diode_1_AlphaPed_RMS/D");
    m_ntuplePtr->Branch("LASER_Diode_2_AlphaPed_RMS", &m_las_D2_AlphaPed_RMS, "LASER_Diode_2_AlphaPed_RMS/D");
    m_ntuplePtr->Branch("LASER_Diode_3_AlphaPed_RMS", &m_las_D3_AlphaPed_RMS, "LASER_Diode_3_AlphaPed_RMS/D");
    m_ntuplePtr->Branch("LASER_Diode_4_AlphaPed_RMS", &m_las_D4_AlphaPed_RMS, "LASER_Diode_4_AlphaPed_RMS/D");
 
    m_ntuplePtr->Branch("LASER_PMT_1_ADC", &m_las_PMT1_ADC, "LASER_PMT_1_ADC/I");
    m_ntuplePtr->Branch("LASER_PMT_2_ADC", &m_las_PMT2_ADC, "LASER_PMT_2_ADC/I");
 
    m_ntuplePtr->Branch("LASER_PMT_1_TDC", &m_las_PMT1_TDC, "LASER_PMT_1_TDC/I");
    m_ntuplePtr->Branch("LASER_PMT_2_TDC", &m_las_PMT2_TDC, "LASER_PMT_2_TDC/I");
 
    m_ntuplePtr->Branch("LASER_PMT_1_Ped", &m_las_PMT1_Ped, "LASER_PMT_1_Ped/D");
    m_ntuplePtr->Branch("LASER_PMT_2_Ped", &m_las_PMT2_Ped, "LASER_PMT_2_Ped/D");
 
    m_ntuplePtr->Branch("LASER_PMT_1_Ped_RMS", &m_las_PMT1_Ped_RMS, "LASER_PMT_1_Ped_RMS/D");
    m_ntuplePtr->Branch("LASER_PMT_2_Ped_RMS", &m_las_PMT2_Ped_RMS, "LASER_PMT_2_Ped_RMS/D");
 
    m_ntuplePtr->Branch("LASER_HEAD_Temp", &m_las_Temperature, "LASER_HEAD_Temp/D");
  }
 
  if (m_beamIdList[LASE_PTN_FRAG]) {
    m_ntuplePtr->Branch("LasFlag", &m_lasFlag, "LasFlag/s");
  }
 
  if ((m_unpackAdder && m_beamIdList[BEAM_ADC_FRAG])
      || (!m_unpackAdder && m_beamIdList[LASE_ADC_FRAG])) {
    m_ntuplePtr->Branch("Las0", &m_las0, "Las0/F");
    m_ntuplePtr->Branch("Las1", &m_las1, "Las1/F");
    m_ntuplePtr->Branch("Las2", &m_las2, "Las2/F");
    m_ntuplePtr->Branch("Las3", &m_las3, "Las3/F");
    m_ntuplePtr->Branch("LasExtra", &m_lasExtra, "m_lasExtra[4]/F");
  }
 
}

LASER_clearBranch()

void TileTBAANtuple::LASER_clearBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// Clear Tree LASER variables

////////////////////////////////////////////////////////////////////////////

Definition at line 2710 of file TileTBAANtuple.cxx.

{
 
  if (!m_laserObjectKey.empty()) {
 
    m_las_BCID = 0;
    m_las_D1_ADC = 0;
    m_las_D2_ADC = 0;
    m_las_D3_ADC = 0;
    m_las_D4_ADC = 0;
 
    m_las_D1_Ped = 0;
    m_las_D2_Ped = 0;
    m_las_D3_Ped = 0;
    m_las_D4_Ped = 0;
 
    m_las_D1_Alpha = 0;
    m_las_D2_Alpha = 0;
    m_las_D3_Alpha = 0;
    m_las_D4_Alpha = 0;
 
    m_las_PMT1_ADC = 0;
    m_las_PMT2_ADC = 0;
 
    m_las_PMT1_TDC = 0;
    m_las_PMT2_TDC = 0;
 
    m_las_PMT1_Ped = 0;
    m_las_PMT2_Ped = 0;
  }
 
  if (m_beamIdList[LASE_PTN_FRAG]) {
    m_lasFlag=0;
  }
 
  if ((m_unpackAdder && m_beamIdList[BEAM_ADC_FRAG])
      || (!m_unpackAdder && m_beamIdList[LASE_ADC_FRAG])) {
 
    m_las0 = 0.;
    m_las1 = 0.;
    m_las2 = 0.;
    m_las3 = 0.;
 
  }
}

LASEROBJ_addBranch()

void TileTBAANtuple::LASEROBJ_addBranch ( void )

private

LASEROBJ_clearBranch()

void TileTBAANtuple::LASEROBJ_clearBranch ( void )

private

msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

MUON_addBranch()

void TileTBAANtuple::MUON_addBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// /Add MUON variables to the Tree

////////////////////////////////////////////////////////////////////////////

Definition at line 2512 of file TileTBAANtuple.cxx.

{
 
  if ( m_beamIdList[MUON_ADC_FRAG] ||
       m_beamIdList[ADDR_ADC_FRAG] )
    {
      m_ntuplePtr->Branch("MuBackHit",&m_muBackHit,"MuBackHit/F");
      m_ntuplePtr->Branch("MuBackSum",&m_muBackSum,"MuBackSum/F");
 
      if (m_TBperiod < 2015) {
        m_ntuplePtr->Branch("MuBack",&m_muBack,"m_muBack[14]");
        m_ntuplePtr->Branch("MuCalib",&m_muCalib,"m_muCalib[2]");
      } else if (m_TBperiod == 2015) {
        m_ntuplePtr->Branch("MuBack",&m_muBack,"MuBack[8]/F");
      } else {
        m_ntuplePtr->Branch("MuBack",&m_muBack,"MuBack[12]/F");
      }
    }
 
 
}

MUON_clearBranch()

void TileTBAANtuple::MUON_clearBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// Clear Tree MUON variables

////////////////////////////////////////////////////////////////////////////

Definition at line 2594 of file TileTBAANtuple.cxx.

                                          {
  m_muBackHit = 0.;
  m_muBackSum = 0.;
 
}

ntuple_clear()

StatusCode TileTBAANtuple::ntuple_clear

(

)

Definition at line 1969 of file TileTBAANtuple.cxx.

                                        {
 
  TRIGGER_clearBranch();
 
  if (m_beamIdList[MUON_ADC_FRAG] || m_beamIdList[ADDR_ADC_FRAG]) {
    MUON_clearBranch();
  }
 
  if (m_beamIdList[ECAL_ADC_FRAG]) {
    ECAL_clearBranch();
    if (m_TBperiod > 2015) QDC_clearBranch();
  }
 
  LASER_clearBranch();
 
  if (m_unpackAdder) {
    if (m_beamIdList[ADD_FADC_FRAG]) {
      ADDER_clearBranch();
    }
  }
 
  if (m_beamIdList[DIGI_PAR_FRAG & 0x1F]) {
    CISPAR_clearBranch();
  }
 
  if (((m_TBperiod >= 2015 || m_unpackAdder) && m_beamIdList[BEAM_ADC_FRAG])
      || (!m_unpackAdder && m_beamIdList[COMMON_ADC1_FRAG])
      || (m_TBperiod >= 2022 && m_beamIdList[COMMON_TOF_FRAG])) {
    BEAM_clearBranch();
    ATH_MSG_VERBOSE( "clear branch");
  }
 
  if (m_completeNtuple && m_TBperiod < 2015) {
    ENETOTAL_clearBranch();
  }
 
  COINCBOARD_clearBranch();
 
  DIGI_clearBranch(); // working now
  if (m_saveFelixData) {
    FELIX_clearBranch();
  }
 
  HIT_clearBranch();
  ATH_MSG_DEBUG( "clear() successfully" );
 
  return StatusCode::SUCCESS;
}

ntuple_initialize()

StatusCode TileTBAANtuple::ntuple_initialize

(

const EventContext &

ctx

)

Alg standard interface function LF TODO: We could have a problem with the new feature introduced by Sasha that quit empty fragments and therefore it can increase or decrease the number of RawChannels fragments dynamically.

The code seems to be capable of react if a module is missing during a run, but if rawchans for a module appear after the beginning of a run (ntuple varialbes are branched af the beginnning of run), its quantities will not be in the ntuple.

if first event: determine mode (normal or calib) by examining the first TileDigitCollection initialize NTuple accordingly

Definition at line 135 of file TileTBAANtuple.cxx.

                                                                    {
 
  CHECK( m_thistSvc.retrieve() );
 
  // find TileCablingService
  m_cabling = TileCablingService::getInstance();
 
  // retrieve TileID helper from det store
 
  CHECK( detStore()->retrieve(m_tileID) );
 
  CHECK( detStore()->retrieve(m_tileHWID) );
 
  if (m_nSamples < 0) {
    m_nSamples = 7;
  }
  if (m_nSamplesFlx < 0) {
    m_nSamplesFlx = m_saveFelixData ? 32 : 0;
  }
 
  if (m_TBperiod >= 2015)  {
    m_unpackAdder = false;
 
    if (m_TBperiod == 2015) {
      m_nDrawers = 2;
      m_drawerList.value().resize(m_nDrawers);
      m_drawerType.value().resize(m_nDrawers);
      m_drawerList[0] = "0x200"; m_drawerType[0] = 2; // barrel neg
      m_drawerList[1] = "0x401"; m_drawerType[1] = 4; // ext.barrel neg
    } else if (m_TBperiod == 2016 || m_TBperiod == 2018 || m_TBperiod == 2021 || m_TBperiod == 2022) {
      m_nDrawers = 5;
      m_drawerList.value().resize(m_nDrawers); m_drawerType.value().resize(m_nDrawers);
      m_drawerList[0] = "0x100"; m_drawerType[0] = 1; // M0 pos
      m_drawerList[1] = "0x101"; m_drawerType[1] = 1; // barrel pos
      m_drawerList[2] = "0x200"; m_drawerType[2] = 2; // M0 neg
      m_drawerList[3] = "0x201"; m_drawerType[3] = 2; // barrel neg
      m_drawerList[4] = "0x402"; m_drawerType[4] = 4; // ext.barrel neg
    } else if (m_TBperiod == 2017) {
      m_nDrawers = 6;
      m_drawerList.value().resize(m_nDrawers); m_drawerType.value().resize(m_nDrawers);
      m_drawerList[0] = "0x100"; m_drawerType[0] = 1; // M0 pos
      m_drawerList[1] = "0x101"; m_drawerType[1] = 1; // barrel pos
      m_drawerList[2] = "0x200"; m_drawerType[2] = 2; // M0 neg
      m_drawerList[3] = "0x201"; m_drawerType[3] = 2; // barrel neg
      m_drawerList[4] = "0x203"; m_drawerType[4] = 2; // barrel neg
      m_drawerList[5] = "0x402"; m_drawerType[5] = 4; // ext.barrel neg
    } else if (m_TBperiod == 2019) {
      m_nDrawers = 7;
      m_drawerList.value().resize(m_nDrawers); m_drawerType.value().resize(m_nDrawers);
      m_drawerList[0] = "0x100"; m_drawerType[0] = 1; // M0 pos
      m_drawerList[1] = "0x101"; m_drawerType[1] = 1; // barrel pos
      m_drawerList[2] = "0x200"; m_drawerType[2] = 2; // M0 neg
      m_drawerList[3] = "0x201"; m_drawerType[3] = 2; // barrel neg
      m_drawerList[4] = "0x203"; m_drawerType[4] = 2; // barrel neg
      m_drawerList[5] = "0x402"; m_drawerType[5] = 4; // ext.barrel neg
      m_drawerList[6] = "0x405"; m_drawerType[6] = 4; // ext.barrel neg
    }
 
    if (m_TBperiod < 2016) {
      setupBeamChambersTB2015();
    } else if (m_TBperiod < 2021) {
      setupBeamChambersTB2016_2020();
    } else {
      // Strating from 2021 the following properties should be setup via JO
      checkIsPropertySetup(m_beamBC1X1, "BC1X1");
      checkIsPropertySetup(m_beamBC1X2, "BC1X2");
      checkIsPropertySetup(m_beamBC1Y1, "BC1Y1");
      checkIsPropertySetup(m_beamBC1Y2, "BC1Y2");
      checkIsPropertySetup(m_beamBC1Z, "BC1Z");
      checkIsPropertySetup(m_beamBC1Z_0, "BC1Z_0");
      checkIsPropertySetup(m_beamBC1Z_90, "BC1Z_90");
      checkIsPropertySetup(m_beamBC1Z_min90, "BC1Z_min90");
 
      checkIsPropertySetup(m_beamBC2X1, "BC2X1");
      checkIsPropertySetup(m_beamBC2X2, "BC2X2");
      checkIsPropertySetup(m_beamBC2Y1, "BC2Y1");
      checkIsPropertySetup(m_beamBC2Y2, "BC2Y2");
      checkIsPropertySetup(m_beamBC2Z, "BC2Z");
      checkIsPropertySetup(m_beamBC2Z_0, "BC2Z_0");
      checkIsPropertySetup(m_beamBC2Z_90, "BC2Z_90");
      checkIsPropertySetup(m_beamBC2Z_min90, "BC2Z_min90");
    }
 
  } else {
    setupBeamChambersBeforeTB2015();
  }
 
 
  if (m_unpackAdder) {
    // get TileRawChannelBuilderFlatFilter for adder energy calculation
 
    StatusCode sc;
    sc &= m_adderFilterAlgTool->setProperty("TileRawChannelContainer", "TileAdderFlat");
    sc &= m_adderFilterAlgTool->setProperty("calibrateEnergy", "true");
 
    sc &= m_adderFilterAlgTool->setProperty("PedStart", "0");
    sc &= m_adderFilterAlgTool->setProperty("PedLength", "1");
    sc &= m_adderFilterAlgTool->setProperty("PedOffset", "0");
    sc &= m_adderFilterAlgTool->setProperty("SignalStart", "1");
    sc &= m_adderFilterAlgTool->setProperty("SignalLength", "15");
    sc &= m_adderFilterAlgTool->setProperty("FilterLength", "5");
    sc &= m_adderFilterAlgTool->setProperty("FrameLength", "16");
    sc &= m_adderFilterAlgTool->setProperty("DeltaCutLo", "9.5");
    sc &= m_adderFilterAlgTool->setProperty("DeltaCutHi", "9.5");
    sc &= m_adderFilterAlgTool->setProperty("RMSCutLo", "1.0");
    sc &= m_adderFilterAlgTool->setProperty("RMSCutHi", "1.0");
    if (sc.isFailure()) {
      ATH_MSG_ERROR("Failure setting properties of " << m_adderFilterAlgTool);
      return StatusCode::FAILURE;
    }
  }
 
  if (m_finalUnit < TileRawChannelUnit::ADCcounts
      || m_finalUnit > TileRawChannelUnit::OnlineMegaElectronVolts) {
 
    m_finalUnit = -1;
    if (!m_useDspUnits && m_calibrateEnergy) {
      m_useDspUnits = true;
      ATH_MSG_INFO( "Final offline units are not set, will use DSP units" );
    }
  }
 
  if (!m_calibrateEnergy) {
    if (m_useDspUnits) {
      ATH_MSG_INFO( "calibrateEnergy is disabled, don't want to use DSP units"  );
      m_useDspUnits = false;
    }
  }
 
  ATH_MSG_INFO( "TestBeam period " << m_TBperiod.value() );
  ATH_MSG_INFO( "calibMode " << m_calibMode.value() );
  ATH_MSG_INFO( "calibrateEnergy " << m_calibrateEnergy.value() );
  ATH_MSG_INFO( "offlineUnits " << m_finalUnit.value() );
  ATH_MSG_INFO( "useDspUnits " << m_useDspUnits.value() );
  ATH_MSG_INFO( "number of samples " << m_nSamples.value() );
 
 
  msg(MSG::INFO) << "drawerList " << MSG::hex;
  unsigned int size = m_drawerList.size();
  for (unsigned int dr = 0; dr < size; ++dr) {
    int frag = strtol(m_drawerList[dr].data(), NULL, 0);
    if (frag >= 0) {
      m_drawerMap[frag] = dr;
      if (dr == m_drawerType.size()) m_drawerType.value().push_back(frag >> 8);
      msg(MSG::INFO) << " 0x" << frag;
    } else {
      msg(MSG::INFO) << " " << m_drawerList[dr];
      // put negative number in first element (flag to read frag ID from data)
      m_drawerList[0] = m_drawerList[dr];
      m_drawerList.value().resize(1);
      // m_drawerType.clear();
      m_drawerMap.clear();
      size = 0;
      break;
    }
  }
  if (size == 0) {
    if (m_drawerList.size() > 0)
      msg(MSG::INFO) << " - negative number, will read frag IDs from the data" << MSG::dec << endmsg;
    else
      msg(MSG::INFO) << "is empty, no drawer fragments in ntuple" << MSG::dec << endmsg;
  } else {
    msg(MSG::INFO) << MSG::dec << endmsg;
 
    size = m_drawerType.size();
    if (size < m_nDrawers) {
      m_drawerType.value().resize(m_nDrawers);
      for (; size < m_nDrawers; ++size)
        m_drawerType[size] = 0;
    }
 
    msg(MSG::INFO) << MSG::INFO << "drawerType ";
    for (unsigned int dr = 0; dr < size; ++dr)
      msg(MSG::INFO) << " " << m_drawerType[dr];
    msg(MSG::INFO) << endmsg;
 
    if (size > m_nDrawers) {
      ATH_MSG_INFO( "increasing m_nDrawers from " << m_nDrawers << " to " << size );
      m_nDrawers = size;
    }
  }
 
  msg(MSG::INFO) << MSG::INFO << "Beam Frag List " << MSG::hex;
  size = m_beamFragList.size();
  for (unsigned int dr = 0; dr < size; ++dr) {
    int frag = strtol(m_beamFragList[dr].data(), NULL, 0);
    if (frag >= 0) {
      m_beamIdList[frag & 0x1F] = true;
      msg(MSG::INFO) << " 0x" << frag;
    }
  }
  if (size == 0) {
    msg(MSG::INFO) << "is empty, no beam fragments in ntuple" << MSG::dec << endmsg;
  } else {
    msg(MSG::INFO) << MSG::dec << endmsg;
  }
 
  // set event number to 0 before first event
  m_evtNr = 0;
  m_calibrateEnergyThisEvent = m_calibrateEnergy;
 
  // find event and beam ROD header, calib mode for digitizers

  if (m_evtNr == 0) {
    if (m_unpackAdder) // in 2003 event number starts from 1
      ++m_evtNr;
  }
 
  if (initList(ctx).isFailure()) {
    ATH_MSG_ERROR( " Error during drawer list initialization"  );
  }
 
  if (m_saveFelixData && initListFlx(ctx).isFailure()) {
    ATH_MSG_ERROR( " Error during drawer list initialization"  );
  }
 
  if (initNTuple().isFailure()) {
    ATH_MSG_ERROR( " Error during ntuple initialization" );
  }
 
  ATH_MSG_INFO( "initialization completed" );
  return StatusCode::SUCCESS;
}

outputHandles()

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

overridevirtualinherited

Return this algorithm's output handles.

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

QDC_addBranch()

void TileTBAANtuple::QDC_addBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// /Add QDC variables to the Tree

////////////////////////////////////////////////////////////////////////////

Definition at line 2555 of file TileTBAANtuple.cxx.

{
 
  if (m_beamIdList[ECAL_ADC_FRAG]) {
    for(unsigned i=0; i<33; ++i){ m_qdc[i]=0.0; }
    m_ntuplePtr->Branch("qdc", &m_qdc, "qdc[33]/i");
  }
 
}

QDC_clearBranch()

void TileTBAANtuple::QDC_clearBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// Clear Tree QDC variables

////////////////////////////////////////////////////////////////////////////

Definition at line 2617 of file TileTBAANtuple.cxx.

{
  for(unsigned i=0; i<33; ++i){ m_qdc[i]=0.0; }
}

renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

setFilterPassed() [1/2]

void AthAlgorithm::setFilterPassed ( bool state ) const

inherited

Definition at line 102 of file AthAlgorithm.cxx.

                                                     {
  setFilterPassed( state, Gaudi::Hive::currentContext() );
}

setFilterPassed() [2/2]

void AthAlgorithm::setFilterPassed ( bool state, const EventContext & ctx ) const

inherited

Definition at line 106 of file AthAlgorithm.cxx.

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

setupBeamChambersBeforeTB2015()

void TileTBAANtuple::setupBeamChambersBeforeTB2015 ( void )

private

Definition at line 3772 of file TileTBAANtuple.cxx.

                                                       {
  // Setup default value for the following properties, if they are not setup via JO
  setupPropertyDefaultValue(m_beamBC1X1, -0.938, "BC1X1");
  setupPropertyDefaultValue(m_beamBC1X2, 0.1747, "BC1X2");
  setupPropertyDefaultValue(m_beamBC1Y1, 0.125, "BC1Y1");
  setupPropertyDefaultValue(m_beamBC1Y2, 0.1765, "BC1Y2");
  setupPropertyDefaultValue(m_beamBC1Z, 13788.0, "BC1Z");
  setupPropertyDefaultValue(m_beamBC1Z_0, 13788.0, "BC1Z_0");
  setupPropertyDefaultValue(m_beamBC1Z_90, 13788.0, "BC1Z_90");
  setupPropertyDefaultValue(m_beamBC1Z_min90, 13788.0, "BC1Z_min90");
 
  setupPropertyDefaultValue(m_beamBC2X1, -0.9369, "BC2X1");
  setupPropertyDefaultValue(m_beamBC2X2, 0.191, "BC2X2");
  setupPropertyDefaultValue(m_beamBC2Y1, -1.29, "BC2Y1");
  setupPropertyDefaultValue(m_beamBC2Y2, 0.187, "BC2Y2");
  setupPropertyDefaultValue(m_beamBC2Z, 9411.0, "BC2Z");
  setupPropertyDefaultValue(m_beamBC2Z_0, 9411.0, "BC2Z_0");
  setupPropertyDefaultValue(m_beamBC2Z_90, 9411.0, "BC2Z_90");
  setupPropertyDefaultValue(m_beamBC2Z_min90, 9411.0, "BC2Z_min90");
}

setupBeamChambersTB2015()

void TileTBAANtuple::setupBeamChambersTB2015 ( void )

private

Definition at line 3793 of file TileTBAANtuple.cxx.

                                                 {
  setupPropertyDefaultValue(m_beamBC1X1, -0.0462586, "BC1X1");
  setupPropertyDefaultValue(m_beamBC1X2, -0.175666, "BC1X2");
  setupPropertyDefaultValue(m_beamBC1Y1, -0.051923, "BC1Y1");
  setupPropertyDefaultValue(m_beamBC1Y2, -0.176809, "BC1Y2");
  setupPropertyDefaultValue(m_beamBC1Z, 13000. + 2760, "BC1Z");
  setupPropertyDefaultValue(m_beamBC1Z_0, 13788.0, "BC1Z_0");
  setupPropertyDefaultValue(m_beamBC1Z_90, 13788.0, "BC1Z_90");
  setupPropertyDefaultValue(m_beamBC1Z_min90, 13788.0, "BC1Z_min90");
 
  setupPropertyDefaultValue(m_beamBC2X1, 0.25202, "BC2X1");
  setupPropertyDefaultValue(m_beamBC2X2, -0.18053, "BC2X2");
  setupPropertyDefaultValue(m_beamBC2Y1, 0.0431688, "BC2Y1");
  setupPropertyDefaultValue(m_beamBC2Y2, -0.181128, "BC2Y2");
  setupPropertyDefaultValue(m_beamBC2Z, 2760, "BC2Z");
  setupPropertyDefaultValue(m_beamBC2Z_0, 9411.0, "BC2Z_0");
  setupPropertyDefaultValue(m_beamBC2Z_90, 9411.0, "BC2Z_90");
  setupPropertyDefaultValue(m_beamBC2Z_min90, 9411.0, "BC2Z_min90");
}

setupBeamChambersTB2016_2020()

void TileTBAANtuple::setupBeamChambersTB2016_2020 ( void )

private

Definition at line 3813 of file TileTBAANtuple.cxx.

                                                      {
  // 2016 settings:
  //
  // https://pcata007.cern.ch/elog/TB2016/88
  //
  //
  // The calibration has been done with the following runs :
  // BC1:
  // Center : 610212
  // Left/up : 610210
  // Right/down : 610209
  //
  // BC2:
  // Center : 610256
  // Left/up : 610321
  // Right/down : 610320
  //
  // Here are the new constants :
  //
  // BC1
  // horizontal slope = -0.171928
  // horizontal offset = -0.047624
  //
  // vertical slope = -0.172942
  // vertical offset = -0.0958677
  //
  // BC2
  // horizontal slope = -0.175698
  // horizontal offset = -1.04599
  //
  // vertical slope = -0.174535
  // vertical offset = -3.10666
 
  // June 2016 calibration
  //m_beamBC1X1 = -0.047624;
  //m_beamBC1X2 = -0.171928;
  //m_beamBC1Y1 = -0.0958677;
  //m_beamBC1Y2 = -0.172942;
  //m_beamBC1Z  = 13000. + 2760 /* 2600. */;
 
  //m_beamBC2X1 = -1.04599;
  //m_beamBC2X2 = -0.175698;
  //m_beamBC2Y1 = -3.10666;
  //m_beamBC2Y2 = -0.174535;
  //m_beamBC2Z  = 2760 /* 2600. */;
 
  // September 2016 calibration, https://pcata007.cern.ch/elog/TB2016/300 (joakim.olsson@cern.ch)
  //m_beamBC1X1 = 0.100857923042;
  //m_beamBC1X2 = -0.172098;
  //m_beamBC1Y1 = -0.133045996607;
  //m_beamBC1Y2 = -0.172855178323;
  //m_beamBC1Z  = 13000. + 2760 /* 2600. */;
  //
  //m_beamBC2X1 = 0.271555258578 ;
  //m_beamBC2X2 = -0.173463 ;
  //m_beamBC2Y1 = 0.305483228502;
  //m_beamBC2Y2 = -0.173805131744 ;
  //m_beamBC2Z  = 2760 /* 2600. */;
 
  // June 2017 calibration, https://pcata007.cern.ch/elog/TB2017/550 (schae@cern.ch)
  //m_beamBC1X1 =  0.153584934082;
  //m_beamBC1X2 = -0.175220;
  //m_beamBC1Y1 = -0.493246053303;
  //m_beamBC1Y2 = -0.176567356723;
  //m_beamBC1Z  = 13000. + 2760 /* 2600. */;
  //
  //m_beamBC2X1 = 0.414611893278;
  //m_beamBC2X2 = -0.176122;
  //m_beamBC2Y1 = 0.150807740888;
  //m_beamBC2Y2 = -0.173472808704;
  //m_beamBC2Z  = 2760 /* 2600. */;
 
  // August 2017 calibration, https://pcata007.cern.ch/elog/TB2017/550 (schae@cern.ch)
  //m_beamBC1X1 =  0.181797;
  //m_beamBC1X2 = -0.175657;
  //m_beamBC1Y1 = -0.128910;
  //m_beamBC1Y2 = -0.175965;
  //m_beamBC1Z  = 13000. + 2760 /* 2600. */;
  //
  //m_beamBC2X1 = 0.611502;
  //m_beamBC2X2 = -0.183116;
  //m_beamBC2Y1 = 0.541212;
  //m_beamBC2Y2 = -0.183115;
  //m_beamBC2Z  = 2760 /* 2600. */;
 
  // September 2017 calibration, https://pcata007.cern.ch/elog/TB2017/550 (schae@cern.ch)
  //m_beamBC1X1 =  0.181797;
  //m_beamBC1X2 = -0.175657;
  //m_beamBC1Y1 = -0.128910;
  //m_beamBC1Y2 = -0.175965;
  //m_beamBC1Z  = 13000. + 2760 /* 2600. */;
 
  //m_beamBC2X1 = 0.622896039922;
  //m_beamBC2X2 = -0.176735;
  //m_beamBC2Y1 = 0.195954125116;
  //m_beamBC2Y2 = -0.176182117624;
  //m_beamBC2Z  = 2760 /* 2600. */;
 
    // September 2017 calibration with additional precision from Survey
  setupPropertyDefaultValue(m_beamBC1X1, 0.681797, "BC1X1");
  setupPropertyDefaultValue(m_beamBC1X2, -0.175657, "BC1X2");
  setupPropertyDefaultValue(m_beamBC1Y1, -2.02891, "BC1Y1");
  setupPropertyDefaultValue(m_beamBC1Y2, -0.175965, "BC1Y2");
  setupPropertyDefaultValue(m_beamBC1Z, 17348.8, "BC1Z");
  setupPropertyDefaultValue(m_beamBC1Z_0, 17348.8, "BC1Z_0");
  setupPropertyDefaultValue(m_beamBC1Z_90, 15594.05, "BC1Z_90");
  setupPropertyDefaultValue(m_beamBC1Z_min90, 15571.8, "BC1Z_min90");
 
  setupPropertyDefaultValue(m_beamBC2X1, -24.377104, "BC2X1");
  setupPropertyDefaultValue(m_beamBC2X2, -0.176735, "BC2X2");
  setupPropertyDefaultValue(m_beamBC2Y1, 17.895954, "BC2Y1");
  setupPropertyDefaultValue(m_beamBC2Y2, -0.176182117624, "BC2Y2");
  setupPropertyDefaultValue(m_beamBC2Z, 4404.2, "BC2Z");
  setupPropertyDefaultValue(m_beamBC2Z_0, 4420.7, "BC2Z_0");
  setupPropertyDefaultValue(m_beamBC2Z_90, 2649.45, "BC2Z_90");
  setupPropertyDefaultValue(m_beamBC2Z_min90, 2627.2, "BC2Z_min90");
 
}

setupPropertyDefaultValue()

void TileTBAANtuple::setupPropertyDefaultValue ( float property, float defaultValue, const std::string & name )

inlineprivate

Definition at line 328 of file TileTBAANtuple.h.

                                                                                              {
      if (property < NOT_SETUP) {
        property = defaultValue;
        ATH_MSG_INFO("The following property is not set up via JO, using default value: " << name << "=" << defaultValue);
      }
    }

storeBeamElements()

StatusCode TileTBAANtuple::storeBeamElements ( const EventContext & ctx )

private

Definition at line 525 of file TileTBAANtuple.cxx.

                                                                    {
 
  if ( m_beamElemContainerKey.empty()) {
    return StatusCode::SUCCESS;
  }
 
  // Read Beam Elements from TES
  const TileBeamElemContainer* beamElemCnt = SG::makeHandle (m_beamElemContainerKey, ctx).get();
 
  TileBeamElemContainer::const_iterator collItr = beamElemCnt->begin();
  TileBeamElemContainer::const_iterator lastColl = beamElemCnt->end();
 
  if ( m_completeNtuple ) {
    // Store ROD header info from collection (just from first one)
      int nDrawersAll = m_nDrawers.value() + m_nDrawersFlx.value();
    if ( collItr!=lastColl ) {
      m_l1ID.at(nDrawersAll) = (*collItr)->getLvl1Id();
      m_l1Type.at(nDrawersAll) = (*collItr)->getLvl1Type();
      m_evType.at(nDrawersAll) = (*collItr)->getDetEvType();
      m_evBCID.at(nDrawersAll) = (*collItr)->getRODBCID();
    } else {
      m_l1ID.at(nDrawersAll) = 0xFFFFFFFF;
      m_l1Type.at(nDrawersAll) = 0xFFFFFFFF;
      m_evType.at(nDrawersAll) = 0xFFFFFFFF;
      m_evBCID.at(nDrawersAll) = 0xFFFFFFFF;
    }
  }
 
 
  m_trigType = 0;
 
  for(; collItr != lastColl; ++collItr) {
 
    TileBeamElemCollection::const_iterator beamItr=(*collItr)->begin();
    TileBeamElemCollection::const_iterator lastBeam=(*collItr)->end();
 
    if (msgLvl(MSG::VERBOSE)) {
 
      for (; beamItr != lastBeam; ++beamItr) {
        HWIdentifier id = (*beamItr)->adc_HWID();
        std::vector<uint32_t> digits = (*beamItr)->get_digits();
        msg(MSG::VERBOSE) << " --- TileBeamElem -- Identifier " << m_tileHWID->to_string(id)
                          << MSG::hex << " frag: 0x" << (*collItr)->identify()
                          << MSG::dec << " channel " << m_tileHWID->channel(id)
                          << " digits size " << digits.size() << endmsg;
        msg(MSG::VERBOSE) << " --- TileBeamElem -- BeamElem : ";
        for (unsigned int k = 0; k < digits.size(); k++)
          msg(MSG::VERBOSE) << digits[k] << " ";
        msg(MSG::VERBOSE) << endmsg;
      }
      //restore iterator
      beamItr = (*collItr)->begin();
    }
 
    int frag = (*collItr)->identify();
    ATH_MSG_DEBUG( " frag: " << frag );
    ATH_MSG_DEBUG( " trigType " << (*collItr)->getLvl1Type() );
 
    if ( m_trigType == 0 && (*collItr)->getLvl1Type() != 0 ) // take it from the ROD header
      m_trigType = (*collItr)->getLvl1Type();
 
    // unpack only fragments which we want to store in ntuple
    if ( m_beamIdList[frag&0x1F] ) {
 
      for (; beamItr != lastBeam; ++beamItr) {
 
        HWIdentifier id = (*beamItr)->adc_HWID();
        std::vector<uint32_t> digits = (*beamItr)->get_digits();
        int cha = m_tileHWID->channel(id);
        int dsize = digits.size();
 
        if ( dsize <= 0 ) {
 
          WRONG_SAMPLE(frag,cha,dsize);
 
        } else if ( dsize != 16 && frag == ADD_FADC_FRAG ) {
 
          WRONG_SAMPLE(frag,cha,dsize);
 
        } else if ( dsize != 1              && frag != ADD_FADC_FRAG    &&
                    frag != BEAM_TDC_FRAG   && frag != COMMON_TDC1_FRAG &&
                    frag != COMMON_TOF_FRAG && frag != COMMON_TDC2_FRAG &&
                    !(frag == ECAL_ADC_FRAG)) {
 
          WRONG_SAMPLE(frag,cha,dsize);
 
        } else {
 
          uint32_t amplitude = digits[0];
          SIGNAL_FOUND(frag, cha, amplitude);
 
          switch (frag) {
 
          case BEAM_TDC_FRAG:
 
            FRAG_FOUND(frag,cha,dsize);
            if(cha < 8) m_btdc1[cha] = amplitude;
            else if(cha < 16) m_btdc2[cha-8] = amplitude;
            else WRONG_CHANNEL(frag,cha);
            break;
 
          case BEAM_ADC_FRAG:
        
            if ( m_TBperiod >= 2015 ) {
              switch(cha) {
                  // BEAM
                case 0: m_s1cou = amplitude; break;
                case 1: m_s2cou = amplitude; break;
                case 2: m_s3cou = amplitude; break;
                case 3: m_cher1 = amplitude; break; // ATH_MSG_VERBOSE("load beam adc " << m_cher1); break;
                case 4: m_cher2 = amplitude; break;
                case 5: m_muTag = amplitude; break;
                case 6: m_muHalo= amplitude; break;
                case 7: m_muVeto= amplitude; break;
                default: WRONG_CHANNEL(frag, cha);
              }
            } else if ( m_unpackAdder ) {
              switch(cha) {
                  // BEAM
                case 0: m_s1cou = amplitude; break;
                case 1: m_s2cou = amplitude; break;
                case 2: m_s3cou = amplitude; break;
                case 3: m_cher1 = amplitude; break; // swap of Cher1
                case 4: m_muTag = amplitude; break; // and S4 in 2003 data
                case 5: m_cher2 = amplitude; break;
                case 6: m_muHalo= amplitude; break;
                case 7: m_muVeto= amplitude; break;
                // LASER
                case 8:  m_las0 = amplitude; break;
                case 9:  m_las1 = amplitude; break;
                case 10: m_las2 = amplitude; break;
                case 11: m_las3 = amplitude; break;
                case 12: m_lasExtra[0] = amplitude; break;
                case 13: m_lasExtra[1] = amplitude; break;
                case 14: m_lasExtra[2] = amplitude; break;
                case 15: m_lasExtra[3] = amplitude; break;
                default: WRONG_CHANNEL(frag, cha);
              }
            } else { // 2004 data
              switch(cha) {
                // BEAM
                case 0: m_sc1 = amplitude; break;
                case 1: m_sc2 = amplitude; break;
                case 2:
                case 3:
                case 4:
                case 5:
                case 6:
                case 7: break;
                default: WRONG_CHANNEL(frag, cha);
              }
            }
            break;
 
          case MUON_ADC_FRAG:
 
            // first 14 are m_muBack
            if(cha < 14) m_muBack[cha] = amplitude;
            // last 2 are m_muCalib
            else if (cha < 16) m_muCalib[cha - 14] = amplitude;
            else WRONG_CHANNEL(frag, cha);
            break;
 
          case ADDR_ADC_FRAG:
 
            // second half of Muon Wall in 2004
            if(cha < 6) m_muBack[cha + 8] = amplitude;
            // last 2 are m_muCalib
            else if (cha < 8) m_muCalib[cha - 6] = amplitude;
            else WRONG_CHANNEL(frag,cha);
            break;
 
          case LASE_PTN_FRAG:
 
            if (cha == 0) {
              // laser pattern unit
              m_lasFlag  = amplitude;
              if (amplitude & 0xFF00) m_trigType = amplitude >> 8;
            } else {
              WRONG_CHANNEL(frag, cha);
            }
            break;
 
          case LASE_ADC_FRAG:
 
            // laser in 2004
            switch(cha) {
              case 0: m_las0 = amplitude; break;
              case 1: m_las1 = amplitude; break;
              case 2: m_las2 = amplitude; break;
              case 3: m_las3 = amplitude; break;
              case 4: m_lasExtra[0] = amplitude; break;
              case 5: m_lasExtra[1] = amplitude; break;
              case 6: m_lasExtra[2] = amplitude; break;
              case 7: m_lasExtra[3] = amplitude; break;
              default: WRONG_CHANNEL(frag,cha);
            }
            break;
 
          case ADD_FADC_FRAG:
 
            if (m_unpackAdder) {
 
              for (int k = 0; k < dsize; k++) {
                //m_addx[k]=k;
                //check how the matrix is filled
                m_adder[cha][k] = digits[k];
              }
                // FlatFiler adders
              double ene, tim;
              m_adderFilterAlgTool->flatFilter(digits, 0, ene, tim);
              m_eneAdd[cha] = ene;
              m_timeAdd[cha] = tim;
            }
 
            break;
 
          case ECAL_ADC_FRAG:
 
            if (m_TBperiod > 2015) {
 
              if(cha < 15) {
                m_qdc[cha] = amplitude;
                ATH_MSG_VERBOSE( "QDC: " << cha << " amp: " << amplitude);
              } else if (cha == 15) {
                for (int idx = 0; idx < dsize && idx < 18; ++idx) {
                  m_qdc[idx + 15] = digits[idx];
                  ATH_MSG_VERBOSE("QDC2: " << cha << " amp: " << amplitude);
                }
              } else {
                WRONG_CHANNEL(frag, cha);
              }
 
            } else {
              if(cha < 8) m_ecal[cha] = amplitude;
              else WRONG_CHANNEL(frag, cha);
            }
 
            break;
 
          case DIGI_PAR_FRAG:
 
            if(cha < 16) m_cispar[cha] = amplitude; //m_cispar->at(cha)=amplitude;
            else WRONG_CHANNEL(frag,cha);
            break;
 
          case COMMON_ADC1_FRAG:
            if (m_TBperiod > 2015) {
              if (cha < 16) {
                if (m_run > 2211444) {
                  switch(cha) {
                    // BEAM
                    case 0: m_s1cou = amplitude; break;
                    case 1: m_s2cou = amplitude; break;
                    case 2: { 
                      if (m_run < 2310000) {
                        m_muBack[10] = amplitude;
                      } else {
                        m_s3cou = amplitude;
                      }
                    }
                      break;
                    case 3: m_cher1 = amplitude; break;
                    case 4: m_cher2 = amplitude; break;
                    case 5: m_cher3 = amplitude; break;
                    default: m_muBack[cha - 6] = amplitude;
                  }
                } else {
                  switch(cha) {
                    // BEAM
                    case 0: m_s1cou = amplitude; break;
                    case 1: m_s2cou = amplitude; break;
                    case 2: m_s3cou = amplitude; break;
                    case 3: m_cher1 = amplitude; break;
                    case 4: m_cher2 = amplitude; break;
                    case 5: m_cher3 = amplitude; break;
                  }
                }
              } else {
                WRONG_CHANNEL(frag, cha);
              }
            } else {
              switch(cha) {
                // BEAM
                case 0: m_s1cou = amplitude; break;
                case 1: m_s2cou = amplitude; break;
                case 2: m_s3cou = amplitude; break;
                case 3: m_muTag = amplitude; break;
                case 4: m_cher1 = amplitude; break;
                case 5: m_cher2 = amplitude; break;
                case 6: m_muHalo= amplitude; break;
                case 7: m_muVeto= amplitude; break;
                default: WRONG_CHANNEL(frag, cha);
              }
            }
            break;
 
          case COMMON_ADC2_FRAG:
 
            if (m_TBperiod > 2015) {
              if(cha < 14) {
                m_muBack[cha] = amplitude;
              } else {
                WRONG_CHANNEL(frag, cha);
              }
            } else {
              if ( ! m_unpackAdder ) {
                switch(cha) {
                  // BEAM
                  case 0: break;
                  case 1: m_s2extra = amplitude; break;
                  case 2: m_s3extra = amplitude; break;
                  case 3:
                  case 4:
                  case 5:
                  case 6:
                  case 7: break;
                  default: WRONG_CHANNEL(frag, cha);
                }
              }
            }
            break;
 
          case COMMON_PTN_FRAG:
            if (m_run > 2310000 && cha < 16) {
              m_scaler[cha] = amplitude;
            } else if (cha == 0) {
              m_commonPU = amplitude;
            } else {
              WRONG_CHANNEL(frag, cha);
            }
            break;
 
          case COMMON_TOF_FRAG:
 
            if (m_TBperiod >= 2022) {
              if (cha > 11) { // The first 12 channels are (can be) connected to BC1 and BC2, the last 4 channels are supposed to be TOF
                if(cha < 16) {
                  m_tof[cha] = amplitude;
                  ATH_MSG_VERBOSE( "TOF: " << cha << " amp: " << amplitude);
                } else {
                  WRONG_CHANNEL(frag, cha);
                }
                break;
              }
              // Fall through to case COMMON_TDC1_FRAG to unpack the first 12 channels of BC1 and BC2
              [[fallthrough]]; // silent the warning on fall through
            } else if (m_TBperiod > 2015) {
 
              if(cha < 16) {
                m_tof[cha] = amplitude;
                ATH_MSG_VERBOSE( "TOF: " << cha << " amp: " << amplitude);
              } else {
                WRONG_CHANNEL(frag, cha);
              }
              break;
            } else {
              if(cha < 8) m_tof[cha] = amplitude;
              else WRONG_CHANNEL(frag, cha);
              break;
            }
 
          case COMMON_TDC1_FRAG:
 
            FRAG_FOUND(frag,cha,dsize);
            if ((cha > 11) && (cha < 16) && (m_run > 2211136)) {
              m_tof[cha] = amplitude;
              ATH_MSG_VERBOSE( "TOF: " << cha << " amp: " << amplitude);
            } if(cha < 16) {
              m_btdc1[cha] = amplitude;
              ATH_MSG_VERBOSE( "TDC: " << cha << " amp: " << amplitude);
              if (m_btdcNhit[cha]==0) {
                m_btdc2[cha] = amplitude;
              }
              (*m_btdc)[cha].push_back(amplitude);
              ++m_btdcNhit[cha];
            } else WRONG_CHANNEL(frag, cha);
            break;
 
          case COMMON_TDC2_FRAG:
 
            FRAG_FOUND(frag,cha,dsize);
            if(cha < 16) m_btdc2[cha] = amplitude;
            else WRONG_CHANNEL(frag, cha);
            break;
 
          case COIN_TRIG1_FRAG:
 
            if(cha < 3) {
              int idx = cha * 32;
              for (int ibit = 0; ibit < 32; ++ibit){
                m_coincTrig1[idx++] = (amplitude >> ibit) & 1;
              }
            } else if (cha == 3) {
              m_coincFlag1 = amplitude;
            } else WRONG_CHANNEL(frag, cha);
 
            break;
 
          case COIN_TRIG2_FRAG:
 
            if(cha < 3) {
              int idx = cha * 32;
              for (int ibit=0; ibit < 32; ++ibit){
                m_coincTrig2[idx++] = (amplitude >> ibit) & 1;
              }
            } else if (cha == 3) {
              m_coincFlag2 = amplitude;
            } else WRONG_CHANNEL(frag, cha);
 
            break;
 
          case COIN_TRIG3_FRAG:
 
            if(cha < 3) {
              int idx = cha * 32;
              for (int ibit = 0; ibit < 32; ++ibit){
                m_coincTrig3[idx++] = (amplitude >> ibit) & 1;
              }
            } else if (cha == 3) {
              m_coincFlag3 = amplitude;
            } else WRONG_CHANNEL(frag,cha);
 
            break;
 
          case COIN_TRIG4_FRAG:
 
            if(cha < 3) {
              int idx = cha * 32;
              for (int ibit = 0; ibit < 32; ++ibit){
                m_coincTrig4[idx++] = (amplitude >> ibit) & 1;
              }
            } else if (cha == 3) {
              m_coincFlag4 = amplitude;
            } else WRONG_CHANNEL(frag, cha);
 
            break;
 
          case COIN_TRIG5_FRAG:
 
            if(cha < 3) {
              int idx = cha * 32;
              for (int ibit = 0; ibit < 32; ++ibit){
                m_coincTrig5[idx++] = (amplitude >> ibit) & 1;
              }
            } else if (cha == 3) {
              m_coincFlag5 = amplitude;
            } else WRONG_CHANNEL(frag, cha);
 
            break;
 
          case COIN_TRIG6_FRAG:
 
            if(cha < 3) {
              int idx = cha * 32;
              for (int ibit = 0; ibit < 32; ++ibit){
                m_coincTrig6[idx++] = (amplitude >> ibit) & 1;
              }
            } else if (cha == 3) {
              m_coincFlag6 = amplitude;
            } else WRONG_CHANNEL(frag, cha);
 
            break;
 
          case COIN_TRIG7_FRAG:
 
            if(cha < 3) {
              int idx = cha * 32;
              for (int ibit = 0; ibit < 32; ++ibit){
                m_coincTrig7[idx++] = (amplitude >> ibit) & 1;
              }
            } else if (cha == 3) {
              m_coincFlag7 = amplitude;
            } else WRONG_CHANNEL(frag,cha);
 
            break;
 
          case COIN_TRIG8_FRAG:
 
            if(cha < 3) {
              int idx = cha * 32;
              for (int ibit = 0; ibit < 32; ++ibit){
                m_coincTrig8[idx++] = (amplitude >> ibit) & 1;
              }
            } else if (cha == 3) {
              m_coincFlag8 = amplitude;
            } else WRONG_CHANNEL(frag, cha);
 
            break;
 
          default:
            break;
          }
        }
      }
    }
  }
 
  for (int i=0; i<8; ++i) {
    if (m_btdcNhit[i] > 1) ++m_btdcNchMultiHit[i>>2];
  }
  // calculate beam coords in Beam Chambers
  if ( m_TBperiod >= 2015 ) {
 
//      For BC1
//      -------
//      m_xCha1 = -0.0462586 + (-0.175666)*(m_btdc1[1] - m_btdc1[0]);
//      m_yCha1 = -0.051923 + (-0.176809)*(m_btdc1[2] - m_btdc1[3]);
//
//      For BC2
//      -------
//      m_xCha2 = 0.25202 + (-0.18053)*(m_btdc1[5] - m_btdc1[4]);
//      m_yCha2 = 0.0431688 + (-0.181128)*(m_btdc1[6] - m_btdc1[7]);
 
      if (m_run > 2211444) {
        m_xCha1 = m_beamBC1X1 + m_beamBC1X2*(m_btdc1[8] - m_btdc1[0]);
        m_yCha1 = m_beamBC1Y1 + m_beamBC1Y2*(m_btdc1[9] - m_btdc1[3]);
      } else {
        m_xCha1 = m_beamBC1X1 + m_beamBC1X2*(m_btdc1[1] - m_btdc1[0]);
        m_yCha1 = m_beamBC1Y1 + m_beamBC1Y2*(m_btdc1[2] - m_btdc1[3]);
      }
      if (m_run > 612543 && m_run< 614141) {
        m_xCha2 = m_beamBC2X1 + m_beamBC2X2*(m_btdc1[5] - m_btdc1[6]);
        m_yCha2 = m_beamBC2Y1 + m_beamBC2Y2*(m_btdc1[4] - m_btdc1[7]);
      } else {
        m_xCha2 = m_beamBC2X1 + m_beamBC2X2*(m_btdc1[5] - m_btdc1[4]);
        m_yCha2 = m_beamBC2Y1 + m_beamBC2Y2*(m_btdc1[6] - m_btdc1[7]);
      }
 
      // Using the first value from the TDC
      if (m_run > 2211444) {
        m_xCha1_0 = m_beamBC1X1 + m_beamBC1X2*(m_btdc2[8] - m_btdc2[0]);
        m_yCha1_0 = m_beamBC1Y1 + m_beamBC1Y2*(m_btdc2[9] - m_btdc2[3]);
      } else {
        m_xCha1_0 = m_beamBC1X1 + m_beamBC1X2*(m_btdc2[1] - m_btdc2[0]);
        m_yCha1_0 = m_beamBC1Y1 + m_beamBC1Y2*(m_btdc2[2] - m_btdc2[3]);
      }
      m_xCha2_0 = m_beamBC2X1 + m_beamBC2X2*(m_btdc2[5] - m_btdc2[4]);
      m_yCha2_0 = m_beamBC2Y1 + m_beamBC2Y2*(m_btdc2[6] - m_btdc2[7]);
 
      m_tjitter = m_btdc1[8];
      m_tscTOF  = m_btdc1[14];
 
      m_xImp  = m_xCha2 + (m_xCha2 - m_xCha1)*m_beamBC2Z/(m_beamBC1Z - m_beamBC2Z);
      m_yImp  = m_yCha2 + (m_yCha2 - m_yCha1)*m_beamBC2Z/(m_beamBC1Z - m_beamBC2Z);
 
// Work in progress
 
      m_xImp_0  = m_xCha2_0 + (m_xCha2_0 - m_xCha1_0)*m_beamBC2Z_0/(m_beamBC1Z_0 - m_beamBC2Z_0);
      m_yImp_0  = m_yCha2_0 + (m_yCha2_0 - m_yCha1_0)*m_beamBC2Z_0/(m_beamBC1Z_0 - m_beamBC2Z_0);
      m_xImp_90  = m_xCha2_0 + (m_xCha2_0 - m_xCha1_0)*m_beamBC2Z_90/(m_beamBC1Z_90 - m_beamBC2Z_90);
      m_yImp_90  = m_yCha2_0 + (m_yCha2_0 - m_yCha1_0)*m_beamBC2Z_90/(m_beamBC1Z_90 - m_beamBC2Z_90);
      m_xImp_min90  = m_xCha2_0 + (m_xCha2_0 - m_xCha1_0)*m_beamBC2Z_min90/(m_beamBC1Z_min90 - m_beamBC2Z_min90);
      m_yImp_min90  = m_yCha2_0 + (m_yCha2_0 - m_yCha1_0)*m_beamBC2Z_min90/(m_beamBC1Z_min90 - m_beamBC2Z_min90);
 
      if (m_run > 2211444) {
        ATH_MSG_DEBUG( "BC1x  : ( "<< m_btdc1[0] <<" - "<< m_btdc1[8] <<" )\t" << m_xCha1 );
        ATH_MSG_DEBUG( "BC1y  : ( "<< m_btdc1[2] <<" - "<< m_btdc1[9] <<" )\t" << m_yCha1 );
      } else {
        ATH_MSG_DEBUG( "BC1x  : ( "<< m_btdc1[0] <<" - "<< m_btdc1[1] <<" )\t" << m_xCha1 );
        ATH_MSG_DEBUG( "BC1y  : ( "<< m_btdc1[2] <<" - "<< m_btdc1[3] <<" )\t" << m_yCha1 );
      }
      ATH_MSG_DEBUG( "BC2x  : ( "<< m_btdc1[4] <<" - "<< m_btdc1[5] <<" )\t" << m_xCha2 );
      ATH_MSG_DEBUG( "BC2y  : ( "<< m_btdc1[6] <<" - "<< m_btdc1[7] <<" )\t" << m_yCha2 );
 
  } else if ( m_unpackAdder ) { // this is 2003 data
 
    if ( m_beamIdList[BEAM_TDC_FRAG] ) {
      m_xCha1 = m_beamBC1X1 + m_beamBC1X2*(m_btdc1[6] - m_btdc1[7]); // last two channels of TDC !!!
      m_yCha1 = m_beamBC1Y1 + m_beamBC1Y2*(m_btdc1[2] - m_btdc1[3]);
      m_xCha2 = m_beamBC2X1 + m_beamBC2X2*(m_btdc2[0] - m_btdc2[1]);
      m_yCha2 = m_beamBC2Y1 + m_beamBC2Y2*(m_btdc2[2] - m_btdc2[3]);
 
      m_xImp  = m_xCha2 + (m_xCha2 - m_xCha1)*m_beamBC2Z/(m_beamBC1Z - m_beamBC2Z);
      m_yImp  = m_yCha2 + (m_yCha2 - m_yCha1)*m_beamBC2Z/(m_beamBC1Z - m_beamBC2Z);
    }
 
  } else { // this is 2004 data
 
    if ( m_beamIdList[COMMON_TDC1_FRAG] ) {
      m_yChN2 = m_beamBN2Y1 + m_beamBN2Y2*(m_btdc1[0] - m_btdc1[1]);
      m_xChN2 = m_beamBN2X1 + m_beamBN2X2*(m_btdc1[2] - m_btdc1[3]);
      m_yChN1 = m_beamBN1Y1 + m_beamBN1Y2*(m_btdc1[4] - m_btdc1[5]);
      m_xChN1 = m_beamBN1X1 + m_beamBN1X2*(m_btdc1[6] - m_btdc1[7]);
 
      m_yCha0 = m_beamBC0Y1 + m_beamBC0Y2*(m_btdc1[8] - m_btdc1[9]);
      m_xCha0 = m_beamBC0X1 + m_beamBC0X2*(m_btdc1[10]- m_btdc1[11]);
 
      ATH_MSG_DEBUG( "BC-2x : ( "<< m_btdc1[2] <<" - "<< m_btdc1[3] <<" )\t" <<m_xChN2  );
      ATH_MSG_DEBUG( "BC-2y : ( "<< m_btdc1[0] <<" - "<< m_btdc1[1] <<" )\t" <<m_yChN2 );
      ATH_MSG_DEBUG( "BC-1x : ( "<< m_btdc1[6] <<" - "<< m_btdc1[7] <<" )\t" <<m_xChN1 );
      ATH_MSG_DEBUG( "BC-1y : ( "<< m_btdc1[4] <<" - "<< m_btdc1[5] <<" )\t" <<m_yChN1 );
      ATH_MSG_DEBUG( "BC0x  : ( "<< m_btdc1[10] <<" - "<< m_btdc1[11] <<" )\t" <<m_xCha0 );
      ATH_MSG_DEBUG( "BC0y  : ( "<< m_btdc1[8] <<" - "<< m_btdc1[9] <<" )\t"   <<m_yCha0 );
 
    }
 
    if ( m_beamIdList[COMMON_TDC2_FRAG] ) {
      m_yCha1 = m_beamBC1Y1 + m_beamBC1Y2*(m_btdc2[0] - m_btdc2[1]);
      m_xCha1 = m_beamBC1X1 + m_beamBC1X2*(m_btdc2[2] - m_btdc2[3]);
      m_yCha2 = m_beamBC2Y1 + m_beamBC2Y2*(m_btdc2[4] - m_btdc2[5]);
      m_xCha2 = m_beamBC2X1 + m_beamBC2X2*(m_btdc2[6] - m_btdc2[7]);
 
      m_xImp = 0.0;
      m_yImp = 0.0;
 
      ATH_MSG_DEBUG( "BC1x  : ( "<< m_btdc2[2] <<" - "<< m_btdc2[3] <<" )\t" <<m_xCha1 );
      ATH_MSG_DEBUG( "BC1y  : ( "<< m_btdc2[0] <<" - "<< m_btdc2[1] <<" )\t" <<m_yCha1 );
      ATH_MSG_DEBUG( "BC2x  : ( "<< m_btdc2[6] <<" - "<< m_btdc2[7] <<" )\t" <<m_xCha2 );
      ATH_MSG_DEBUG( "BC2y  : ( "<< m_btdc2[4] <<" - "<< m_btdc2[5] <<" )\t" <<m_yCha2 );
 
    }
 
    if ( m_beamIdList[COMMON_TDC1_FRAG] && m_beamIdList[COMMON_TDC2_FRAG] ) {

      // Get run number and eta
      if ( 0==m_runNumber ){
        const xAOD::EventInfo* eventInfo(0);
        if (evtStore()->retrieve(eventInfo).isFailure()){
          ATH_MSG_ERROR( "No EventInfo object found! Can't read run number!" );
          m_runNumber = -1;
        } else {
          m_runNumber = eventInfo->runNumber();
          getEta();
        }
      }
 
      // Computation of X,Y imp on TileCal/LAr front
      float tanBx = (m_xCha1-m_xCha0) / (m_beamBC0Z - m_beamBC1Z);
      float tanBy = (m_yCha1-m_yCha0) / (m_beamBC0Z - m_beamBC1Z);
 
      float Ximp = m_xCha1 + m_beamBC1Z * tanBx;
      float Yimp = m_yCha1 + m_beamBC1Z * tanBy;
 
      if ( 0.0 != cos(m_theta) * (1 + tanBx * tan(m_theta)) ){
        m_xImp = (Ximp + m_radius *(tanBx * (cos(m_theta) - 1) -sin(m_theta))) / (cos(m_theta) * (1 + tanBx * tan(m_theta)));
        m_yImp = Yimp + tanBy * (m_radius * (1 - cos(m_theta)) -Ximp * sin(m_theta)) / (cos(m_theta) * (1+tanBx * tan(m_theta)));
      }
    }
  }
 
 
  // do not apply Cesium and Laser calibration for CIS events
  // m_calibrateEnergyThisEvent = m_calibrateEnergy && (m_trigType != 8);
 
  return StatusCode::SUCCESS;
}

storeCells()

StatusCode TileTBAANtuple::storeCells ( const EventContext & ctx )

private

Definition at line 2129 of file TileTBAANtuple.cxx.

                                                             {
 
  if (m_cellContainerKey.empty()) { // empty name, nothing to do
    return StatusCode::FAILURE;
  }
 
  //Retrieve Cell collection from SG
  const CaloCellContainer* cellContainer = SG::makeHandle (m_cellContainerKey, ctx).get();
 
  //Loop over all cells in container. Sum up the Energy and fill 2DHistograms
  ATH_MSG_DEBUG( "succeeded retrieving cellContainer from SG" );
 
  ATH_MSG_DEBUG( "TileTBAANtuple : about to iterate over CaloCells" );
 
  m_LarEne[0] = m_LarEne[1] = m_LarEne[2] = m_LarEne[3] = 0.0;
  m_BarEne[0] = m_BarEne[1] = m_BarEne[2] = 0.0;
  m_ExtEne[0] = m_ExtEne[1] = m_ExtEne[2] = 0.0;
  m_GapEne[0] = m_GapEne[1] = m_GapEne[2] = 0.0;
 
  for (const CaloCell* cell : *cellContainer) {
    //Decode cell information
    const double energy = cell->energy();
    const CaloDetDescrElement* caloDDE = cell->caloDDE(); //pointer to the DetectorDescriptionElement
    const CaloCell_ID::CaloSample sampl = caloDDE->getSampling(); //To which sampling belongs this cell?
 
    if (sampl == CaloCell_ID::PreSamplerB) {
      m_LarEne[0] += energy;
    } else if (sampl == CaloCell_ID::EMB1) {
      m_LarEne[1] += energy;
    } else if (sampl == CaloCell_ID::EMB2) {
      m_LarEne[2] += energy;
    } else if (sampl == CaloCell_ID::EMB3) {
      m_LarEne[3] += energy;
    } else if (sampl == CaloCell_ID::TileBar0) {
      m_BarEne[0] += energy;
    } else if (sampl == CaloCell_ID::TileBar1) {
      m_BarEne[1] += energy;
    } else if (sampl == CaloCell_ID::TileBar2) {
      m_BarEne[2] += energy;
    } else if (sampl == CaloCell_ID::TileExt0) {
      m_ExtEne[0] += energy;
    } else if (sampl == CaloCell_ID::TileExt1) {
      m_ExtEne[1] += energy;
    } else if (sampl == CaloCell_ID::TileExt2) {
      m_ExtEne[2] += energy;
    } else if (sampl == CaloCell_ID::TileGap1) {
      m_GapEne[1] += energy;
    } else if (sampl == CaloCell_ID::TileGap2) {
      m_GapEne[2] += energy;
    } else if (sampl == CaloCell_ID::TileGap3) {
      m_GapEne[0] += energy;
    }
  }
 
  return StatusCode::SUCCESS;
}

storeDigits()

StatusCode TileTBAANtuple::storeDigits ( const EventContext & ctx, const SG::ReadHandleKey< TileDigitsContainer > & containerKey, const TileDQstatus * dqStatus = nullptr )

private

/ Fill Ntuple with info from TileDigits / Return true if the collection is empty, / which means that there are no RawChanels either.

Store ROD header info from collection (should be just one per ROD, but we don't know how many RODs we have, so store it for every collection)

Full DMU header, stored for debugging

bit_31==1 and bit_17==0

parity must be an odd number

memory parity error bit_25

single strobe error bit_24 (it is recovered)

double strobe error bit_23 (cannot be recovered)

Full DMU header, stored for debugging

bit_31==1 and bit_17==0

parity must be an odd number

memory parity error bit_25

single strobe error bit_24 (it is recovered)

double strobe error bit_23 (cannot be recovered)

Full DMU header, stored for debugging

bit_31==1 and bit_17==0

parity must be an odd number

memory parity error bit_25

single strobe error bit_24 (it is recovered)

double strobe error bit_23 (cannot be recovered)

Definition at line 1344 of file TileTBAANtuple.cxx.

                                                                                                                                                      {
 
  if (containerKey.empty()) { // empty name, nothing to do
    return StatusCode::FAILURE;
  }
 
  // Read Digits from TES
  const TileDigitsContainer* digitsCnt = SG::makeHandle (containerKey, ctx).get();
 
  bool emptyColl = true;
 
  // drawerIndex is 0 - m_nDrawers-1, fragType is 1-4 B+/B-/EB+/EB-
  int drawerIndex, fragType, channel;
 
  std::vector<float> sampleVec;
  std::vector<uint32_t> headerVec;
  std::vector<uint32_t> headerVecHi;
  uint32_t CRCmask;
  uint32_t fe_crc;
  uint32_t rod_crc;
 
  // Go through all TileDigitsCollections
  for (const TileDigitsCollection* digitsCollection : * digitsCnt) {
    // determine type of frag
    int fragId = digitsCollection->identify();
    int ros = (fragId >> 8);
    int drawer = fragId & 0x3F;
 
    drawerMap_iterator itr = m_drawerMap.find(fragId);
    if ( itr != m_drawerMap.end() ) {
      drawerIndex = (*itr).second;
    } else {
      drawerIndex= -1;
    }
 
    if (drawerIndex < 0) {
      if ( !digitsCollection->empty() )
        ATH_MSG_DEBUG( "frag id 0x" << MSG::hex << fragId << MSG::dec <<" was not found among valid frag IDs when storing TRC!" );
 
    } else {
      fragType = m_drawerType[drawerIndex];
 
      ATH_MSG_DEBUG( "Event# " << m_evtNr
                     << " Frag id 0x" << MSG::hex << fragId << MSG::dec
                     << " index " << drawerIndex
                     << " Calib " << m_calibMode );
 
      ATH_MSG_DEBUG( "       Size=" << digitsCollection->getFragSize()
                    << " BCID=" << digitsCollection->getFragBCID()<<MSG::hex
                    << " CRC=0x" << (digitsCollection->getFragCRC()&0xffff)
                    << " DMUMask=0x" << (digitsCollection->getFragDMUMask()&0xffff)<<MSG::dec );
 
      ATH_MSG_DEBUG( "       Lvl1ID=" << digitsCollection->getLvl1Id()
                     << " Lvl1Type=" << digitsCollection->getLvl1Type()
                     << " EvBCID=" << digitsCollection->getRODBCID()
                     << " EvType=" << digitsCollection->getDetEvType() );
 
      ATH_MSG_DEBUG( "       Header=" << digitsCollection->getFragChipHeaderWords() );
 
      if (m_completeNtuple) {
        m_l1ID.at(drawerIndex) = digitsCollection->getLvl1Id();
        m_l1Type.at(drawerIndex) = digitsCollection->getLvl1Type();
        m_evType.at(drawerIndex) = digitsCollection->getDetEvType();
        m_evBCID.at(drawerIndex) = digitsCollection->getRODBCID();
        // store FrBCID
        m_frBCID.at(drawerIndex) = digitsCollection->getFragBCID();
      }
 
      if(m_calibMode) {
        // Digits in calib mode
        // check gain for first digits in collection
 
        int dcnt=0;
        int drawerIndexHi = drawerIndex + m_nDrawers;
        // non empty collection
        if (!digitsCollection->empty()) {
          // store evtnr, bcid,crc, size
          // Same for lo and hi, because they come from the same fragment
          m_rodBCIDVec.at(drawerIndex) = digitsCollection->getRODBCID();
          m_sizeVec.at(drawerIndex) = digitsCollection->getFragSize();
          m_sizeVec.at(drawerIndexHi) = digitsCollection->getFragSize();
          m_evtVec.at(drawerIndex) = m_evtNr;
          m_evtVec.at(drawerIndexHi) = m_evtNr;
 
          headerVec = digitsCollection->getFragChipHeaderWords();
          headerVecHi = digitsCollection->getFragChipHeaderWordsHigh();
          CRCmask = digitsCollection->getFragDMUMask(); //mask of FE+ROD DMU crc check (16bit+16bit) 0xffffffff == All ok
          fe_crc = CRCmask & 0xFFFF;
          rod_crc = CRCmask >> 16;
 
          unsigned int headsize = std::min(16U, static_cast<unsigned int>(headerVec.size()));
          unsigned int headsizehi = std::min(16U, static_cast<unsigned int>(headerVecHi.size()));
 
          for (unsigned int ih = 0; ih < headsize; ++ih) {
 
            m_bcidVec.at(drawerIndex)[ih] = (headerVec[ih] & 0xFFF);
            m_DMUheaderVec.at(drawerIndex)[ih] = headerVec[ih]; 
            m_DMUformatErrVec.at(drawerIndex)[ih] = CheckDMUFormat(headerVec[ih]); 
            m_DMUparityErrVec.at(drawerIndex)[ih] = CheckDMUParity(headerVec[ih]); 
            m_DMUmemoryErrVec.at(drawerIndex)[ih] = (headerVec[ih] >> 25 & 0x1); 
            m_DMUSstrobeErrVec.at(drawerIndex)[ih] = (headerVec[ih] >> 24 & 0x1); 
            m_DMUDstrobeErrVec.at(drawerIndex)[ih] = (headerVec[ih] >> 23 & 0x1); 
            m_feCRCVec.at(drawerIndex)[ih] = (fe_crc >> ih & 0x1);
            m_rodCRCVec.at(drawerIndex)[ih] = (rod_crc >> ih & 0x1);
 
            if (dqStatus) {
              m_DMUHeadparityErrVec.at(drawerIndex)[ih] = dqStatus->checkHeaderParityErr(ros, drawer, ih, 0);
              m_DMUDataparityErrVec.at(drawerIndex)[ih] = dqStatus->checkSampleParityErr(ros, drawer, ih, 0);
            }
          }
 
          for (unsigned int ihhi = 0; ihhi < headsizehi; ++ihhi) {
            m_bcidVec.at(drawerIndexHi)[ihhi] = (headerVecHi[ihhi] & 0xFFF);
            m_DMUheaderVec.at(drawerIndexHi)[ihhi] = headerVecHi[ihhi]; 
            m_DMUformatErrVec.at(drawerIndexHi)[ihhi] = CheckDMUFormat(headerVecHi[ihhi]); 
            m_DMUparityErrVec.at(drawerIndexHi)[ihhi] = CheckDMUParity(headerVecHi[ihhi]); 
            m_DMUmemoryErrVec.at(drawerIndexHi)[ihhi] = (headerVecHi[ihhi] >> 25 & 0x1); 
            m_DMUSstrobeErrVec.at(drawerIndexHi)[ihhi] = (headerVecHi[ihhi] >> 24 & 0x1); 
            m_DMUDstrobeErrVec.at(drawerIndexHi)[ihhi] = (headerVecHi[ihhi] >> 23 & 0x1); 
            m_feCRCVec.at(drawerIndex)[ihhi] = -1; //Variables must be filled anyway, empty variables are not allowed
            m_rodCRCVec.at(drawerIndex)[ihhi] = -1; //Variables must be filled anyway, empty variables are not allowed
 
            if (dqStatus) {
              m_DMUHeadparityErrVec.at(drawerIndex)[ihhi] = dqStatus->checkHeaderParityErr(ros, drawer, ihhi, 1);
              m_DMUDataparityErrVec.at(drawerIndex)[ihhi] = dqStatus->checkSampleParityErr(ros, drawer, ihhi, 1);
            }
          }
 
          m_slinkCRCVec.at(drawerIndex)[0] = (digitsCollection->getFragCRC() >> 16) & 0xffff;
          m_dmuMaskVec.at(drawerIndex)[0] = (digitsCollection->getFragDMUMask() >> 16) & 0xffff;
          m_slinkCRCVec.at(drawerIndex)[1] = digitsCollection->getFragCRC() & 0xffff;
          m_dmuMaskVec.at(drawerIndex)[1] = digitsCollection->getFragDMUMask() & 0xffff;
 
          m_slinkCRCVec.at(drawerIndexHi)[0] = (digitsCollection->getFragCRC() >> 16) & 0xffff;
          m_dmuMaskVec.at(drawerIndexHi)[0] = (digitsCollection->getFragDMUMask() >> 16) & 0xffff;
          m_slinkCRCVec.at(drawerIndexHi)[1] = digitsCollection->getFragCRC() & 0xffff;
          m_dmuMaskVec.at(drawerIndexHi)[1] = digitsCollection->getFragDMUMask() & 0xffff;
 
          // go through all TileDigits in collection
          for (const TileDigits* tile_digits : *digitsCollection) {
 
            emptyColl = false;
 
            HWIdentifier hwid = tile_digits->adc_HWID();
            // determine gain
            int gain = m_tileHWID->adc(hwid);
            // add m_nDrawers to index if hi gain
            int index = (gain == 1) ? drawerIndexHi : drawerIndex;
            int nSamplesInDrawer = m_nSamplesInDrawer[index];
 
            // determine channel
            channel = m_tileHWID->channel(hwid);
            // cabling for testbeam (convert to pmt#-1)
 
            if ((m_TBperiod < 2015 ||
                 (m_TBperiod==2015 && fragType<3) ||
                 ((m_TBperiod==2016 || m_TBperiod==2021) && ((fragId&0xFF)<4 && fragId != 0x201)) ||
                 (m_TBperiod==2017 && ((fragId&0xFF)<4 && !(fragId == 0x201 || fragId == 0x203))) ||
                 (m_TBperiod==2018 && ((fragId&0xFF)<4 && !(fragId == 0x201 || fragId == 0x402))) ||
                 (m_TBperiod==2019 && ((fragId&0xFF)<5 && !(fragId == 0x201 || fragId == 0x203 || fragId >= 0x402))) ||
                 (m_TBperiod==2022 && ((fragId&0xFF)<4 && !(fragId == 0x201 || (m_run >= 2210456 && fragId == 0x402)))))
                && fragType > 0 && m_pmtOrder)
 
              channel = digiChannel2PMT(fragType, channel);
 
 
            // gain determined for all digits in collection
            m_gainVec.at(index)[channel] = gain;
            ATH_MSG_DEBUG( "Storing TD for channel: " << channel
                           << " with gain " << m_tileHWID->adc(hwid)
                           << " index " << index );
 
            // get digits
            sampleVec = tile_digits->samples();
            int siz = sampleVec.size();
 
            if (msgLvl(MSG::DEBUG)) {
              msg(MSG::DEBUG) << "Digits(" << siz << ")." << (dcnt++) << " {";
              for (int i = 0; i < siz; i++) {
                msg(MSG::DEBUG) << static_cast<int>(sampleVec[i]) << " ";
              }
 
              if (siz > nSamplesInDrawer) {
                msg(MSG::DEBUG) << "} ONLY " << nSamplesInDrawer << " digits saved to ntuple" << endmsg;
              } else {
                msg(MSG::DEBUG) << "}" << endmsg;
              }
            }
 
            if (siz > nSamplesInDrawer) siz = nSamplesInDrawer;
            std::transform(sampleVec.begin(), sampleVec.begin() + siz, &m_sampleVec.at(index).get()[0] + nSamplesInDrawer * channel, [] (float v) {return static_cast<int>(v);});
          }
        }
      } else {
        // Digits in normal mode
        // store evtnr, bcid,crc, size
        m_rodBCIDVec.at(drawerIndex) = digitsCollection->getRODBCID();
        m_sizeVec.at(drawerIndex) = digitsCollection->getFragSize();
        m_evtVec.at(drawerIndex) = m_evtNr;
 
        headerVec = digitsCollection->getFragChipHeaderWords();
        CRCmask = digitsCollection->getFragDMUMask(); //mask of FE+ROD DMU crc check (16bit+16bit) 0xffffffff == All ok
        fe_crc = CRCmask & 0xFFFF;
        rod_crc = CRCmask >> 16;
 
        int headsize = headerVec.size();
 
        for (int ih = 0; ih < headsize; ++ih) {
          m_bcidVec.at(drawerIndex)[ih] = (headerVec[ih] & 0xFFF);
          m_DMUheaderVec.at(drawerIndex)[ih] = headerVec[ih];   
          m_DMUformatErrVec.at(drawerIndex)[ih] = CheckDMUFormat(headerVec[ih]); 
          m_DMUparityErrVec.at(drawerIndex)[ih] = CheckDMUParity(headerVec[ih]); 
          m_DMUmemoryErrVec.at(drawerIndex)[ih] = (headerVec[ih] >> 25 & 0x1); 
          m_DMUSstrobeErrVec.at(drawerIndex)[ih] = (headerVec[ih] >> 24 & 0x1); 
          m_DMUDstrobeErrVec.at(drawerIndex)[ih] = (headerVec[ih] >> 23 & 0x1); 
          m_feCRCVec.at(drawerIndex)[ih] = (fe_crc >> ih & 0x1);
          m_rodCRCVec.at(drawerIndex)[ih] = (rod_crc >> ih & 0x1);
 
          if (dqStatus) {
            m_DMUHeadparityErrVec.at(drawerIndex)[ih] = dqStatus->checkHeaderParityErr(ros, drawer, ih, 0);
            m_DMUDataparityErrVec.at(drawerIndex)[ih] = dqStatus->checkSampleParityErr(ros, drawer, ih, 0);
          }
        }
 
        m_slinkCRCVec.at(drawerIndex)[0] = (digitsCollection->getFragCRC() >> 16) & 0xffff;
        m_dmuMaskVec.at(drawerIndex)[0] = (digitsCollection->getFragDMUMask() >> 16) & 0xffff;
        m_slinkCRCVec.at(drawerIndex)[1] = digitsCollection->getFragCRC() & 0xffff;
        m_dmuMaskVec.at(drawerIndex)[1] = digitsCollection->getFragDMUMask() & 0xffff;
 
        int nSamplesInDrawer = m_nSamplesInDrawer[drawerIndex];
 
        int dcnt = 0;
        // go through all TileDigits in collection
        for (const TileDigits* tile_digits : *digitsCollection) {
          emptyColl = false;
 
          HWIdentifier hwid = tile_digits->adc_HWID();
          // determine channel
          channel = m_tileHWID->channel(hwid);
          // cabling for testbeam
          if ((m_TBperiod < 2015 ||
               (m_TBperiod==2015 && fragType<3) ||
               ((m_TBperiod==2016 || m_TBperiod==2021) && ((fragId&0xFF)<4 && fragId != 0x201)) ||
               (m_TBperiod==2017 && ((fragId&0xFF)<4 && !(fragId == 0x201 || fragId == 0x203))) ||
               (m_TBperiod==2018 && ((fragId&0xFF)<4 && !(fragId == 0x201 || fragId == 0x402))) ||
               (m_TBperiod==2019 && ((fragId&0xFF)<5 && !(fragId == 0x201 || fragId == 0x203 || fragId >= 0x402))) ||
               (m_TBperiod==2022 && ((fragId&0xFF)<4 && !(fragId == 0x201 || (m_run >= 2210456 && fragId == 0x402)))))
              && fragType > 0 && m_pmtOrder)
            channel = digiChannel2PMT(fragType, channel);
 
          // gain
          m_gainVec.at(drawerIndex)[channel] = m_tileHWID->adc(hwid);
          ATH_MSG_DEBUG( "Storing TD for channel: " << channel
                         << " with gain " << m_tileHWID->adc(hwid) );
 
          // get digits
          sampleVec = tile_digits->samples();
          int siz = sampleVec.size();
          if (msgLvl(MSG::DEBUG)) {
            msg(MSG::DEBUG) << "Digits(" << siz << ")." << (dcnt++) << " {";
 
            for (int i = 0; i < siz; i++) {
              msg(MSG::DEBUG) << static_cast<int>(sampleVec[i]) << " ";
            }
 
            if (siz > nSamplesInDrawer) {
              msg(MSG::DEBUG) << "} ONLY " << nSamplesInDrawer << " digits saved to ntuple" << endmsg;
            } else {
              msg(MSG::DEBUG) << "}" << endmsg;
            }
          }
 
          if (siz > nSamplesInDrawer) siz = nSamplesInDrawer;
          std::transform(sampleVec.begin(), sampleVec.begin() + siz, &m_sampleVec.at(drawerIndex).get()[0] + nSamplesInDrawer * channel, [] (float v) {return static_cast<int>(v);});
        }
      }
    }
    // next container
  }
 
  if (emptyColl)
    return StatusCode::FAILURE;
  else
    return StatusCode::SUCCESS;
}

storeDigitsFlx()

StatusCode TileTBAANtuple::storeDigitsFlx ( const EventContext & ctx, const SG::ReadHandleKey< TileDigitsContainer > & containerKey )

private

Store ROD header info from collection (should be just one per ROD, but we don't know how many RODs we have, so store it for every collection)

Definition at line 1635 of file TileTBAANtuple.cxx.

                                                                                                                           {
 
  if (containerKey.empty()) { // empty name, nothing to do
    return StatusCode::FAILURE;
  }
 
  // Read Digits from TES
  const TileDigitsContainer* digitsCntFlx = SG::makeHandle (containerKey, ctx).get();
 
  bool emptyColl = true;
 
  // drawerIndex is 0 - m_nDrawersFlx-1, fragType is 1-4 B+/B-/EB+/EB-
  int drawerIndex, channel;
 
  std::vector<float> sampleVecLo;
  std::vector<float> sampleVecHi;
  // Go through all TileDigitsCollections
  for (const TileDigitsCollection* digitsCollection : *digitsCntFlx) {
    // determine type of frag
    int fragId = digitsCollection->identify();
    drawerMap_iterator itr = m_drawerFlxMap.find(fragId);
    if ( itr != m_drawerFlxMap.end() ) {
      drawerIndex = (*itr).second;
    } else {
      drawerIndex = -1;
    }
 
    if (drawerIndex < 0) {
      if ( !digitsCollection->empty() )
        ATH_MSG_DEBUG( "FELIX frag id 0x" << MSG::hex << fragId << MSG::dec <<" was not found among valid frag IDs when storing TRC!" );
 
    } else {
 
      ATH_MSG_DEBUG( "Event# " << m_evtNr
                     << " FELIX Frag id 0x" << MSG::hex << fragId << MSG::dec
                     << " index " << drawerIndex);
 
      ATH_MSG_DEBUG( "       Size=" << digitsCollection->getFragSize());
 
      ATH_MSG_DEBUG( "       Lvl1ID=" << digitsCollection->getLvl1Id()
                     << " EvBCID=" << digitsCollection->getRODBCID()
                     << " EvType=" << digitsCollection->getDetEvType() );
      ATH_MSG_DEBUG( "       Headers = "<< digitsCollection->getFragExtraWords() );
 
      if (m_completeNtuple) {
        int index = drawerIndex + m_nDrawers;
        m_l1ID.at(index) = digitsCollection->getLvl1Id();
        m_l1Type.at(index) = digitsCollection->getLvl1Type();
        m_evType.at(index) = digitsCollection->getDetEvType();
        m_evBCID.at(index) = digitsCollection->getRODBCID();
        // store FrBCID
        m_frBCID.at(index) = digitsCollection->getFragBCID();
      }
 
      // Digits in calib mode
      // check gain for first digits in collection
 
 
      std::vector<uint32_t> extraWords = digitsCollection->getFragExtraWords();
      if (extraWords.size() >= 10 * MAX_MINIDRAWER) {
 
        std::reference_wrapper<std::array<int,MAX_MINIDRAWER>>
          md[] = {m_mdL1idflxVec.at(drawerIndex), m_mdBcidflxVec.at(drawerIndex),
                  m_mdModuleflxVec.at(drawerIndex), m_mdRunTypeflxVec.at(drawerIndex),
                  m_mdRunflxVec.at(drawerIndex), m_mdPedLoflxVec.at(drawerIndex),
                  m_mdPedHiflxVec.at(drawerIndex), m_mdChargeflxVec.at(drawerIndex),
                  m_mdChargeTimeflxVec.at(drawerIndex), m_mdCapacitorflxVec.at(drawerIndex)};
 
        auto it = extraWords.begin();
        for (int i = 0; i < 10; ++i) {
          std::copy(it + i * MAX_MINIDRAWER, it + (i + 1) * MAX_MINIDRAWER, &md[i].get()[0]);
        }
 
      }
 
      int dcnt=0;
      // non empty collection
      if(!digitsCollection->empty()) {
        int drawerIndexHi = drawerIndex + m_nDrawersFlx;
        // store evtnr, bcid,crc, size
        if (m_bsInput) {
          m_rodBCIDflxVec.at(drawerIndex) = digitsCollection->getRODBCID();
          m_sizeflxVec.at(drawerIndex) = digitsCollection->getFragSize();
          m_evtflxVec.at(drawerIndex) = m_evtNr;
        }
 
        // go through all TileDigits in collection
        for (const TileDigits* tile_digits : *digitsCollection) {
          emptyColl = false;
          HWIdentifier hwid = tile_digits->adc_HWID();
          // determine gain
          int gain = m_tileHWID->adc(hwid);
          // add m_nDrawersFlx to index if hi gain
          int index = (gain == 1) ? drawerIndexHi : drawerIndex;
 
          // determine channel
          channel = m_tileHWID->channel(hwid);
          // cabling for testbeam (convert to pmt#-1)
 
          // gain determined for all digits in collection
          m_gainflxVec.at(index)[channel] = gain;
          ATH_MSG_DEBUG( "Storing TD for channel: " << channel
                         << " with gain " << m_tileHWID->adc(hwid)
                         << " index " << index );
 
          // get digits
          if (gain == 0) {
            sampleVecLo = tile_digits->samples();
          } else if (gain == 1) {
            sampleVecHi = tile_digits->samples();
          }
 
          int sizLo = sampleVecLo.size();
          int sizHi = sampleVecHi.size();
 
          int nSamplesInDrawer = m_nSamplesFlxInDrawer[index];
 
          if (msgLvl(MSG::DEBUG)) {
            if (sizLo > 0 ){
              msg(MSG::DEBUG) << "Low gain Digits(" << sizLo << ")." << (dcnt++) << " {";
              for (int i = 0; i < sizLo; i++) {
                msg(MSG::DEBUG) << static_cast<int>(sampleVecLo[i]) << " ";
              }
 
              if (sizLo > nSamplesInDrawer) {
                msg(MSG::DEBUG) << "} ONLY " << nSamplesInDrawer << " digits saved to ntuple" << endmsg;
              } else {
                msg(MSG::DEBUG) << "}" << endmsg;
              }
            }
 
            if (sizHi > 0 ){
              msg(MSG::DEBUG) << "High gain Digits(" << sizHi << ")." << (dcnt++) << " {";
              for (int i = 0; i < sizHi; i++) {
                msg(MSG::DEBUG) << static_cast<int>(sampleVecHi[i]) << " ";
              }
 
              if (sizHi > nSamplesInDrawer) {
                msg(MSG::DEBUG) << "} ONLY " << nSamplesInDrawer << " digits saved to ntuple" << endmsg;
              } else {
                msg(MSG::DEBUG) << "}" << endmsg;
              }
            }
 
 
          }
          if (sizLo > nSamplesInDrawer) sizLo = nSamplesInDrawer;
          if (sizHi > nSamplesInDrawer) sizHi = nSamplesInDrawer;
          std::transform(sampleVecLo.begin(), sampleVecLo.begin() + sizLo, &m_sampleflxVec.at(index).get()[0] + nSamplesInDrawer * channel, [] (float v) {return static_cast<int>(v);});
          std::transform(sampleVecHi.begin(), sampleVecHi.begin() + sizHi, &m_sampleflxVec.at(index).get()[0] + nSamplesInDrawer * channel, [] (float v) {return static_cast<int>(v);});
          sampleVecLo.clear();
          sampleVecHi.clear();
        }
      }
    }
    // next container
  }
 
  if (emptyColl)
    return StatusCode::FAILURE;
  else
    return StatusCode::SUCCESS;
}

storeHit()

void TileTBAANtuple::storeHit ( const TileHit * hit, int fragType, int fragId, std::array< float, MAX_CHAN > & ehitVec, std::array< float, MAX_CHAN > & thitVec, const TileSamplingFraction * samplingFraction )

private

Definition at line 1900 of file TileTBAANtuple.cxx.

                                                                            {
 
  // determine channel
  HWIdentifier hwid = cinp->pmt_HWID();
  int channel = m_tileHWID->channel(hwid);
 
  int size = cinp->size();
  if (msgLvl(MSG::VERBOSE)) {
    msg(MSG::VERBOSE) << "hit hwid="
                      << m_tileHWID->to_string(hwid, -1) << " ener=";
 
    for (int i = 0; i < size; ++i)
      msg(MSG::VERBOSE) << cinp->energy(i) << " ";
 
    msg(MSG::VERBOSE) << "time=";
    for (int i = 0; i < size; ++i)
      msg(MSG::VERBOSE) << cinp->time(i) << " ";
 
    msg(MSG::VERBOSE) << endmsg;
  }
 
  double ehit=0.0, thit=0.0;
  for(int i=0;i<size;++i) {
 
    double e = cinp->energy(i);
    double t = cinp->time(i);
 
    if (-75.<t && t<75.) {
      ehit += e;
      thit += e*t;
    }
  }
 
  if (ehit!=0) {
    thit /= ehit;
    // conversion factor from hit energy to final energy units 
    int drawerIdx = TileCalibUtils::getDrawerIdxFromFragId(fragId);
    ehit *= samplingFraction->getSamplingFraction(drawerIdx, channel);
    if (m_rchUnit != TileRawChannelUnit::MegaElectronVolts) {
      ehit /= m_tileToolEmscale->channelCalib(drawerIdx, channel, TileID::HIGHGAIN, 1.,
                                              m_rchUnit, TileRawChannelUnit::MegaElectronVolts);
    }
  } else {
    thit=0.0;
  }
 
  // cabling for testbeam
  if ((m_TBperiod < 2015 || 
       (m_TBperiod==2015 && fragType<3) ||
       ((m_TBperiod==2016 || m_TBperiod==2021) && ((fragId&0xFF)<4 && fragId != 0x201)) ||
       (m_TBperiod==2017 && ((fragId&0xFF)<4 && !(fragId == 0x201 || fragId == 0x203))) ||
       (m_TBperiod==2018 && ((fragId&0xFF)<4 && !(fragId == 0x201 || fragId == 0x402))) ||
       (m_TBperiod==2019 && ((fragId&0xFF)<5 && !(fragId == 0x201 || fragId == 0x203 || fragId >= 0x402))) ||
       (m_TBperiod==2022 && ((fragId&0xFF)<4 && !(fragId == 0x201 || (m_run >= 2210456 && fragId == 0x402)))))
      && fragType > 0 && m_pmtOrder)
    channel = digiChannel2PMT(fragType, channel);
 
  ehitVec[channel] = ehit;
  thitVec[channel] = thit;
 
  ATH_MSG_DEBUG( "HIT ene=" << ehit
                 << " time=" << thit
                 << " pmt-1=" << channel
                 << " index " << m_drawerMap.at(fragId) );
 
}

storeHitContainer()

StatusCode TileTBAANtuple::storeHitContainer ( const EventContext & ctx )

private

/ Fill Ntuple with MC truth info from simulation / Namely, hit energies corrected by photoelectron statistics and Birks' law

Definition at line 1853 of file TileTBAANtuple.cxx.

                                                                    {
 
  if (m_hitContainerKey.empty()) { // empty name, nothing to do
    return StatusCode::FAILURE;
  }
 
  // Read Hit Vector from TDS
  const TileHitContainer* hitCnt = SG::makeHandle (m_hitContainerKey, ctx).get();
 
  SG::ReadCondHandle<TileSamplingFraction> samplingFraction(m_samplingFractionKey, ctx);
  ATH_CHECK( samplingFraction.isValid() );
 
  bool emptyColl = true;
 
  // Go through all TileHitCollections
  for (const TileHitCollection* hitCollection : *hitCnt) {
 
    // determine type of frag
    int fragId = hitCollection->identify();
    drawerMap_iterator itr = m_drawerMap.find(fragId);
    int drawerIndex = ( itr != m_drawerMap.end() ) ? static_cast<int>((*itr).second) : -1;
 
    if (drawerIndex < 0) {
      if ( !hitCollection->empty() )
        ATH_MSG_WARNING( "frag id 0x" << MSG::hex << fragId << MSG::dec <<" was not found among valid frag IDs when storing HITS!" );
 
    } else {
      int fragType = m_drawerType[drawerIndex];
 
      ATH_MSG_DEBUG( "Event# " << m_evtNr
                     << " Frag id 0x" << MSG::hex << fragId << MSG::dec
                     << " index " << drawerIndex );
 
      if (emptyColl) emptyColl = hitCollection->empty();
      // go through all TileHit in collection
      for (const TileHit* cinp : *hitCollection) {
        storeHit(cinp,fragType,fragId,m_ehitCnt.at(drawerIndex),m_thitCnt.at(drawerIndex),*samplingFraction);
      }
    }
  }
 
  if (emptyColl)
    return StatusCode::FAILURE;
  else
    return StatusCode::SUCCESS;
}

storeHitVector()

StatusCode TileTBAANtuple::storeHitVector ( const EventContext & ctx )

private

Definition at line 1808 of file TileTBAANtuple.cxx.

                                                                 {
 
  if (m_hitVectorKey.empty()) { // empty name, nothing to do
    return StatusCode::FAILURE;
  }
 
  // Read Hit Vector from TDS
  const TileHitVector* hitVec = SG::makeHandle (m_hitVectorKey, ctx).get();
 
  ATH_MSG_DEBUG( "Event# " << m_evtNr << " reading Hit Vector");
 
  SG::ReadCondHandle<TileSamplingFraction> samplingFraction(m_samplingFractionKey, ctx);
  ATH_CHECK( samplingFraction.isValid() );
 
  // Go through all TileHit
  for (const TileHit& cinp : *hitVec) {
 
    // get hits
    HWIdentifier hwid = cinp.pmt_HWID();
 
    // determine type of frag
    int fragId = m_tileHWID->frag(hwid);
    drawerMap_iterator itr = m_drawerMap.find(fragId);
    int drawerIndex = ( itr != m_drawerMap.end() ) ? static_cast<int>((*itr).second) : -1;
 
    if (drawerIndex < 0) {
      ATH_MSG_WARNING( "frag id 0x" << MSG::hex << fragId << MSG::dec <<" was not found among valid frag IDs when storing HITS!" );
 
    } else {
      int fragType = m_drawerType[drawerIndex];
      storeHit(&cinp,fragType,fragId,m_ehitVec.at(drawerIndex),m_thitVec.at(drawerIndex), *samplingFraction);
    }
  }
 
  if (hitVec->empty())
    return StatusCode::FAILURE;
  else
    return StatusCode::SUCCESS;
}

storeLaser()

StatusCode TileTBAANtuple::storeLaser ( const EventContext & ctx )

private

Definition at line 455 of file TileTBAANtuple.cxx.

                                                             {
 
  if (m_laserObjectKey.empty()) { // empty name, nothing to do
    return StatusCode::FAILURE;
  }
 
  // Read Hit Vector from TDS
  const TileLaserObject* laserObj = SG::makeHandle (m_laserObjectKey, ctx).get();
 
  m_las_BCID   = laserObj->getBCID();
 
  m_las_Filt = laserObj->getFiltNumber();
  m_las_ReqAmp = laserObj->getDiodeCurrOrd();
  m_las_MeasAmp = laserObj->getDiodeCurrMeas();
 
 
  m_las_D1_ADC = laserObj->getDiodeADC(0);
  m_las_D2_ADC = laserObj->getDiodeADC(1);
  m_las_D3_ADC = laserObj->getDiodeADC(2);
  m_las_D4_ADC = laserObj->getDiodeADC(3);
 
  m_las_D1_Ped = laserObj->getDiodePedestal(0);
  m_las_D2_Ped = laserObj->getDiodePedestal(1);
  m_las_D3_Ped = laserObj->getDiodePedestal(2);
  m_las_D4_Ped = laserObj->getDiodePedestal(3);
 
  m_las_D1_Ped_RMS = laserObj->getDiodeSigmaPedestal(0);
  m_las_D2_Ped_RMS = laserObj->getDiodeSigmaPedestal(1);
  m_las_D3_Ped_RMS = laserObj->getDiodeSigmaPedestal(2);
  m_las_D4_Ped_RMS = laserObj->getDiodeSigmaPedestal(3);
 
  m_las_D1_Alpha = laserObj->getAlpha(0);
  m_las_D2_Alpha = laserObj->getAlpha(1);
  m_las_D3_Alpha = laserObj->getAlpha(2);
  m_las_D4_Alpha = laserObj->getAlpha(3);
 
  m_las_D1_Alpha_RMS = laserObj->getSigmaAlpha(0);
  m_las_D2_Alpha_RMS = laserObj->getSigmaAlpha(1);
  m_las_D3_Alpha_RMS = laserObj->getSigmaAlpha(2);
  m_las_D4_Alpha_RMS = laserObj->getSigmaAlpha(3);
 
  m_las_D1_AlphaPed = laserObj->getPedestalAlpha(0);
  m_las_D2_AlphaPed = laserObj->getPedestalAlpha(1);
  m_las_D3_AlphaPed = laserObj->getPedestalAlpha(2);
  m_las_D4_AlphaPed = laserObj->getPedestalAlpha(3);
 
  m_las_D1_AlphaPed_RMS = laserObj->getSigmaPedAlpha(0);
  m_las_D2_AlphaPed_RMS = laserObj->getSigmaPedAlpha(1);
  m_las_D3_AlphaPed_RMS = laserObj->getSigmaPedAlpha(2);
  m_las_D4_AlphaPed_RMS = laserObj->getSigmaPedAlpha(3);
 
  m_las_PMT1_ADC = laserObj->getPMADC(0);
  m_las_PMT2_ADC = laserObj->getPMADC(1);
 
  m_las_PMT1_TDC = laserObj->getTDC(0);
  m_las_PMT2_TDC = laserObj->getTDC(1);
 
  m_las_PMT1_Ped = laserObj->getPMPedestal(0);
  m_las_PMT2_Ped = laserObj->getPMPedestal(1);
 
  m_las_PMT1_Ped_RMS = laserObj->getPMSigmaPedestal(0);
  m_las_PMT2_Ped_RMS = laserObj->getPMSigmaPedestal(1);
 
  m_las_Temperature = laserObj->getPumpDiodeTemp();
 
  ATH_MSG_DEBUG( "storeLaser() completed" );
 
  return StatusCode::SUCCESS;
}

storeRawChannels()

StatusCode TileTBAANtuple::storeRawChannels ( const EventContext & ctx, const SG::ReadHandleKey< TileRawChannelContainer > & containerKey, bool calib_mode, std::vector< std::array< float, MAX_CHAN > > * eneVec, std::vector< std::array< float, MAX_CHAN > > * timeVec, std::vector< std::array< float, MAX_CHAN > > * chi2Vec, std::vector< std::array< float, MAX_CHAN > > * pedVec, bool saveDQstatus = false )

private

/ Fill ntuple with data from TRC.

/ Default TRC container contains flat filtered and / named container fitted/opt filtered

final calibration

Definition at line 1182 of file TileTBAANtuple.cxx.

{
 
 
  if (containerKey.empty()) {// empty name, nothing to do
    return StatusCode::FAILURE;
  }
 
  bool isFELIX = containerKey.key().find("Flx") != std::string::npos;
  int nDrawers = isFELIX ? m_nDrawersFlx : m_nDrawers;
  std::map<unsigned int, unsigned int, std::less<unsigned int>>& drawerMap = (isFELIX) ? m_drawerFlxMap : m_drawerMap;
 
  // get named container
  const TileRawChannelContainer* rcCnt = SG::makeHandle (containerKey, ctx).get();
 
  TileRawChannelUnit::UNIT rChUnit = rcCnt->get_unit();
  ATH_MSG_DEBUG( "RawChannel unit is " << rChUnit );
 
 
  if (rChUnit >= TileRawChannelUnit::OnlineADCcounts) { // this is container with DSP results
    m_dspUnit = rChUnit;
    m_dspFlags = rcCnt->get_bsflags() >> 16;
    ATH_MSG_DEBUG( "DSP flag is 0x" << MSG::hex << m_dspFlags << MSG::dec  << " DSP unit is " << m_dspUnit );
 
  } else if ((m_useDspUnits || m_finalUnit >= TileRawChannelUnit::OnlineADCcounts)
             && rChUnit != TileRawChannelUnit::ADCcounts) {
 
      ATH_MSG_ERROR( "RawChannel units are not ADC counts, can't apply DSP-like calibration" );
      return StatusCode::FAILURE;
  }
 
  if (m_calibrateEnergyThisEvent) {
    if (m_useDspUnits) { // calibrate a-la online
      m_rchUnit = m_dspUnit;
    } else { // convert to final units
      m_rchUnit = (TileRawChannelUnit::UNIT)m_finalUnit.value();
    }
  } else {
    m_rchUnit = rChUnit;
  }
 
  ATH_MSG_DEBUG( "Final RawChannel unit is " << m_rchUnit );
 
  // drawerIndex is 0 - m_nDrawers-1, fragType is 1-4 B+/B-/EB+/EB-
  int drawerIndex, fragType;
 
  // Go through all TileRawChannelCollections
  for (const TileRawChannelCollection* rawChannelCollection : *rcCnt) {
    // determine type
    int fragId = rawChannelCollection->identify();
    int drawerIdx = TileCalibUtils::getDrawerIdxFromFragId(fragId);
 
    drawerMap_iterator itr = drawerMap.find(fragId);
    if ( itr != drawerMap.end() ) {
      drawerIndex = (*itr).second;
    } else {
      drawerIndex= -1;
    }
 
    if (drawerIndex < 0) {
      if ( !rawChannelCollection->empty() )
        ATH_MSG_DEBUG( "frag id 0x" << MSG::hex << fragId << MSG::dec << " was not found among valid frag IDs when storing TRC!" );
    } else {
      fragType = isFELIX ? fragId >> 8 : m_drawerType[drawerIndex];
      ATH_MSG_DEBUG( "TRC (" << containerKey.key()
                    << ") Event# " << m_evtNr
                    << " Frag id 0x" << MSG::hex << fragId << MSG::dec
                    << " index "<< drawerIndex );
 
      // go through all TileRawChannels in collection
      for (const TileRawChannel* rch : *rawChannelCollection) {
        int index = drawerIndex;
        HWIdentifier hwid = rch->adc_HWID();
 
        // determine channel and gain
        int channel = m_tileHWID->channel(hwid);
        int gain = m_tileHWID->adc(hwid);
 
        if (calib_mode) {
          // gain, if hi add m_nDrawers to index
          if (gain == 1) index += nDrawers;
        }

        double energy = rch->amplitude();
        if (m_rchUnit != rChUnit) {
          if (m_rchUnit < TileRawChannelUnit::OnlineOffset)
            energy = m_tileToolEmscale->channelCalib(drawerIdx, channel, gain, energy, rChUnit, m_rchUnit);
          else
            energy = m_tileToolEmscale->channelCalibOnl(drawerIdx, channel, gain, energy, m_rchUnit);
        }
 
        // cabling for testbeam (convert to pmt#-1)
        if ((m_TBperiod < 2015 ||
             (m_TBperiod==2015 && fragType<3) ||
             ((m_TBperiod==2016 || m_TBperiod==2021) && ((fragId&0xFF)<4 && fragId != 0x201)) ||
             (m_TBperiod==2017 && ((fragId&0xFF)<4 && !(fragId == 0x201 || fragId == 0x203))) ||
             (m_TBperiod==2018 && ((fragId&0xFF)<4 && !(fragId == 0x201 || fragId == 0x402))) ||
             (m_TBperiod==2019 && ((fragId&0xFF)<5 && !(fragId == 0x201 || fragId == 0x203 || fragId >= 0x402))) ||
             (m_TBperiod==2022 && ((fragId&0xFF)<4 && !(fragId == 0x201 || (m_run >= 2210456 && fragId == 0x402)))))
            && fragType > 0 && m_pmtOrder)
          channel = digiChannel2PMT(fragType, channel);
 
        eneVec->at(index)[channel] = energy;
        timeVec->at(index)[channel] = rch->time();
        if (chi2Vec) {
          chi2Vec->at(index)[channel] = rch->quality();
        }
        if (pedVec) {
          pedVec->at(index)[channel] = rch->pedestal();
        }
 
        ATH_MSG_DEBUG( "TRC ene=" << energy
                       << " time=" << rch->time()
                       << " chi2=" << rch->quality()
                       << " ped=" << rch->pedestal()
                       << " pmt-1=" << channel
                       << " index " << index );
 
      }
 
      if (saveDQstatus && !isFELIX) {
 
        int index1 = drawerIndex, index2 = drawerIndex + 1;
        if (calib_mode) index2 += m_nDrawers;
 
        for (int index = index1; index < index2; index += m_nDrawers) {
 
          m_ROD_GlobalCRCVec.at(index) = rawChannelCollection->getFragGlobalCRC() & 1;
          m_ROD_DMUMaskVec.at(index)[0] = rawChannelCollection->getFragRODChipMask();
          m_ROD_DMUMaskVec.at(index)[1] = rawChannelCollection->getFragFEChipMask();
 
          for (unsigned int dmu = 0; dmu < MAX_DMU; ++dmu) {
 
            m_ROD_DMUBCIDVec.at(index)[dmu] = (rawChannelCollection->getFragBCID() >> dmu) & 1;
            m_ROD_DMUmemoryErrVec.at(index)[dmu] = (rawChannelCollection->getFragMemoryPar() >> dmu) & 1;
            m_ROD_DMUSstrobeErrVec.at(index)[dmu] = (rawChannelCollection->getFragSstrobe() >> dmu) & 1;
            m_ROD_DMUDstrobeErrVec.at(index)[dmu] = (rawChannelCollection->getFragDstrobe() >> dmu) & 1;
            m_ROD_DMUHeadformatErrVec.at(index)[dmu] = (rawChannelCollection->getFragHeaderBit() >> dmu) & 1;
            m_ROD_DMUDataformatErrVec.at(index)[dmu] = (rawChannelCollection->getFragSampleBit() >> dmu) & 1;
          }
        }
      }
    }
  }
 
  return StatusCode::SUCCESS;
}

sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

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

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

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

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

Definition at line 66 of file AthAlgorithm.cxx.

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

sysStart()

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

overridevirtualinherited

Handle START transition.

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

TRIGGER_addBranch()

void TileTBAANtuple::TRIGGER_addBranch ( void )

private

///////////////////////////////////////////////////////////////////////////

/ Variables Legenda / / - C : a character string terminated by the 0 character / - B : an 8 bit signed integer / - b : an 8 bit unsigned integer 2^8=256 / - S : a 16 bit signed integer (i.e. a "short") / - s : a 16 bit unsigned integer 2^16=65536 / - I : a 32 bit signed integer (i.e an "int") / - i : a 32 bit unsigned integer 2^32=4294967296 / - F : a 32 bit floating point (i.e. a "float") / - D : a 64 bit floating point (i.e. a "double") / - L : a 64 bit signed integer / - l : a 64 bit unsigned integer / - O : a boolean //////////////////////////////////////////////////////////////////////////// /Add TRIGGER variables to the Tree

////////////////////////////////////////////////////////////////////////////

Definition at line 2473 of file TileTBAANtuple.cxx.

{
 
  m_ntuplePtr->Branch("EvTime",&m_evTime,"EvTime/I");
  m_ntuplePtr->Branch("Run",&m_run,"Run/I");
  m_ntuplePtr->Branch("Evt",&m_evt,"Evt/I");
  m_ntuplePtr->Branch("Trig",&m_trigType,"Trig/S");
  m_ntuplePtr->Branch("DSPflags",&m_dspFlags,"DSPflags/I");
  m_ntuplePtr->Branch("DSPunits",&m_dspUnit,"DSPunits/S");
  m_ntuplePtr->Branch("OFLunits",&m_rchUnit,"OFLunits/S");
 
  if ( m_completeNtuple ) {
    int nDrawersAll = m_nDrawers.value() + m_nDrawersFlx.value();
    if (nDrawersAll > 0) {
      m_l1ID.resize(nDrawersAll + 1);
      m_l1Type.resize(nDrawersAll + 1);
      m_evType.resize(nDrawersAll + 1);
      m_evBCID.resize(nDrawersAll + 1);
      m_frBCID.resize(nDrawersAll);
    }
 
    // Info from ROD headers
    // m_nDrawers drawers separately (i.e. with duplications) + Beam ROD
    m_ntuplePtr->Branch("L1ID",&m_l1ID);
    m_ntuplePtr->Branch("L1Type",&m_l1Type);
    m_ntuplePtr->Branch("EvType",&m_evType);
    m_ntuplePtr->Branch("EvBCID",&m_evBCID);
    m_ntuplePtr->Branch("FrBCID",&m_frBCID);
  }
 
 
}

TRIGGER_clearBranch()

void TileTBAANtuple::TRIGGER_clearBranch ( void )

private

//////////////////////////////////////////////////////////////////////////// Clear Tree TRIGGER variables

////////////////////////////////////////////////////////////////////////////

Definition at line 2572 of file TileTBAANtuple.cxx.

                                             {
 
  m_evTime = 0;
  m_run = 0;
  m_evt = 0;
  m_trigType = 0;
 
  if (m_completeNtuple) {
    std::fill(m_l1ID.begin(), m_l1ID.end(), 0);
    std::fill(m_l1Type.begin(), m_l1Type.end(), 0);
    std::fill(m_evType.begin(), m_evType.end(), 0);
    std::fill(m_evBCID.begin(), m_evBCID.end(), 0);
    std::fill(m_frBCID.begin(), m_frBCID.end(), 0);
  }
}

updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

m_adder

std::vector<int*> TileTBAANtuple::m_adder

private

Definition at line 444 of file TileTBAANtuple.h.

m_adderFilterAlgTool

ToolHandle<TileRawChannelBuilderFlatFilter> TileTBAANtuple::m_adderFilterAlgTool

private

Initial value:

{this,
        "TileRawChannelBuilderFlatFilter", "TileRawChannelBuilderFlatFilter", "Tile raw channel builder tool"}

Definition at line 159 of file TileTBAANtuple.h.

                                                                    {this,
        "TileRawChannelBuilderFlatFilter", "TileRawChannelBuilderFlatFilter", "Tile raw channel builder tool"};

m_adderPayload

std::vector<int> TileTBAANtuple::m_adderPayload

private

Definition at line 445 of file TileTBAANtuple.h.

m_BarEne

std::array<float, 3> TileTBAANtuple::m_BarEne {}

private

Definition at line 581 of file TileTBAANtuple.h.

{};

m_bcidVec

std::vector<std::array<int, MAX_DMU> > TileTBAANtuple::m_bcidVec

private

Definition at line 515 of file TileTBAANtuple.h.

m_beamBC0X1

Gaudi::Property<float> TileTBAANtuple::m_beamBC0X1 {this, "BC0X1", 0.0, "Params for Beam TDC: Beam chamber: 0"}

private

Definition at line 198 of file TileTBAANtuple.h.

{this, "BC0X1", 0.0, "Params for Beam TDC: Beam chamber: 0"};

m_beamBC0X2

Gaudi::Property<float> TileTBAANtuple::m_beamBC0X2 {this, "BC0X2", 0.2, "Params for Beam TDC: Beam chamber: 0"}

private

Definition at line 199 of file TileTBAANtuple.h.

{this, "BC0X2", 0.2, "Params for Beam TDC: Beam chamber: 0"};

m_beamBC0Y1

Gaudi::Property<float> TileTBAANtuple::m_beamBC0Y1 {this, "BC0Y1", 0.0, "Params for Beam TDC: Beam chamber: 0"}

private

Definition at line 200 of file TileTBAANtuple.h.

{this, "BC0Y1", 0.0, "Params for Beam TDC: Beam chamber: 0"};

m_beamBC0Y2

Gaudi::Property<float> TileTBAANtuple::m_beamBC0Y2 {this, "BC0Y2", 0.2, "Params for Beam TDC: Beam chamber: 0"}

private

Definition at line 201 of file TileTBAANtuple.h.

{this, "BC0Y2", 0.2, "Params for Beam TDC: Beam chamber: 0"};

m_beamBC0Z

Gaudi::Property<float> TileTBAANtuple::m_beamBC0Z {this, "BC0Z", 17138.0, "Params for Beam TDC: Beam chamber: 0"}

private

Definition at line 202 of file TileTBAANtuple.h.

{this, "BC0Z", 17138.0, "Params for Beam TDC: Beam chamber: 0"};

m_beamBC1X1

Gaudi::Property<float> TileTBAANtuple::m_beamBC1X1 {this, "BC1X1", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}

private

Definition at line 204 of file TileTBAANtuple.h.

{this, "BC1X1", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"};

m_beamBC1X2

Gaudi::Property<float> TileTBAANtuple::m_beamBC1X2 {this, "BC1X2", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}

private

Definition at line 205 of file TileTBAANtuple.h.

{this, "BC1X2", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"};

m_beamBC1Y1

Gaudi::Property<float> TileTBAANtuple::m_beamBC1Y1 {this, "BC1Y1", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}

private

Definition at line 206 of file TileTBAANtuple.h.

{this, "BC1Y1", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"};

m_beamBC1Y2

Gaudi::Property<float> TileTBAANtuple::m_beamBC1Y2 {this, "BC1Y2", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}

private

Definition at line 207 of file TileTBAANtuple.h.

{this, "BC1Y2", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"};

m_beamBC1Z

Gaudi::Property<float> TileTBAANtuple::m_beamBC1Z {this, "BC1Z", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}

private

Definition at line 208 of file TileTBAANtuple.h.

{this, "BC1Z", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"};

m_beamBC1Z_0

Gaudi::Property<float> TileTBAANtuple::m_beamBC1Z_0 {this, "BC1Z_0", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}

private

Definition at line 209 of file TileTBAANtuple.h.

{this, "BC1Z_0", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"};

m_beamBC1Z_90

Gaudi::Property<float> TileTBAANtuple::m_beamBC1Z_90 {this, "BC1Z_90", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}

private

Definition at line 210 of file TileTBAANtuple.h.

{this, "BC1Z_90", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"};

m_beamBC1Z_min90

Gaudi::Property<float> TileTBAANtuple::m_beamBC1Z_min90 {this, "BC1Z_min90", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"}

private

Definition at line 211 of file TileTBAANtuple.h.

{this, "BC1Z_min90", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 1"};

m_beamBC2X1

Gaudi::Property<float> TileTBAANtuple::m_beamBC2X1 {this, "BC2X1", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}

private

Definition at line 213 of file TileTBAANtuple.h.

{this, "BC2X1", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"};

m_beamBC2X2

Gaudi::Property<float> TileTBAANtuple::m_beamBC2X2 {this, "BC2X2", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}

private

Definition at line 214 of file TileTBAANtuple.h.

{this, "BC2X2", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"};

m_beamBC2Y1

Gaudi::Property<float> TileTBAANtuple::m_beamBC2Y1 {this, "BC2Y1", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}

private

Definition at line 215 of file TileTBAANtuple.h.

{this, "BC2Y1", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"};

m_beamBC2Y2

Gaudi::Property<float> TileTBAANtuple::m_beamBC2Y2 {this, "BC2Y2", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}

private

Definition at line 216 of file TileTBAANtuple.h.

{this, "BC2Y2", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"};

m_beamBC2Z

Gaudi::Property<float> TileTBAANtuple::m_beamBC2Z {this, "BC2Z", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}

private

Definition at line 217 of file TileTBAANtuple.h.

{this, "BC2Z", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"};

m_beamBC2Z_0

Gaudi::Property<float> TileTBAANtuple::m_beamBC2Z_0 {this, "BC2Z_0", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}

private

Definition at line 218 of file TileTBAANtuple.h.

{this, "BC2Z_0", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"};

m_beamBC2Z_90

Gaudi::Property<float> TileTBAANtuple::m_beamBC2Z_90 {this, "BC2Z_90", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}

private

Definition at line 219 of file TileTBAANtuple.h.

{this, "BC2Z_90", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"};

m_beamBC2Z_min90

Gaudi::Property<float> TileTBAANtuple::m_beamBC2Z_min90 {this, "BC2Z_min90", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"}

private

Definition at line 220 of file TileTBAANtuple.h.

{this, "BC2Z_min90", NOT_SETUP - 1, "Params for Beam TDC: Beam chamber: 2"};

m_beamBN1X1

Gaudi::Property<float> TileTBAANtuple::m_beamBN1X1 {this, "BN1X1", 0.0, "Params for Beam TDC: Beam chamber: -1"}

private

Definition at line 193 of file TileTBAANtuple.h.

{this, "BN1X1", 0.0, "Params for Beam TDC: Beam chamber: -1"};

m_beamBN1X2

Gaudi::Property<float> TileTBAANtuple::m_beamBN1X2 {this, "BN1X2", 0.2, "Params for Beam TDC: Beam chamber: -1"}

private

Definition at line 194 of file TileTBAANtuple.h.

{this, "BN1X2", 0.2, "Params for Beam TDC: Beam chamber: -1"};

m_beamBN1Y1

Gaudi::Property<float> TileTBAANtuple::m_beamBN1Y1 {this, "BN1Y1", 0.0, "Params for Beam TDC: Beam chamber: -1"}

private

Definition at line 195 of file TileTBAANtuple.h.

{this, "BN1Y1", 0.0, "Params for Beam TDC: Beam chamber: -1"};

m_beamBN1Y2

Gaudi::Property<float> TileTBAANtuple::m_beamBN1Y2 {this, "BN1Y2", 0.2, "Params for Beam TDC: Beam chamber: -1"}

private

Definition at line 196 of file TileTBAANtuple.h.

{this, "BN1Y2", 0.2, "Params for Beam TDC: Beam chamber: -1"};

m_beamBN2X1

Gaudi::Property<float> TileTBAANtuple::m_beamBN2X1 {this, "BN2X1", 0.0, "Params for Beam TDC: Beam chamber: -2"}

private

Definition at line 188 of file TileTBAANtuple.h.

{this, "BN2X1", 0.0, "Params for Beam TDC: Beam chamber: -2"};

m_beamBN2X2

Gaudi::Property<float> TileTBAANtuple::m_beamBN2X2 {this, "BN2X2", 0.2, "Params for Beam TDC: Beam chamber: -2"}

private

Definition at line 189 of file TileTBAANtuple.h.

{this, "BN2X2", 0.2, "Params for Beam TDC: Beam chamber: -2"};

m_beamBN2Y1

Gaudi::Property<float> TileTBAANtuple::m_beamBN2Y1 {this, "BN2Y1", 0.0, "Params for Beam TDC: Beam chamber: -2"}

private

Definition at line 190 of file TileTBAANtuple.h.

{this, "BN2Y1", 0.0, "Params for Beam TDC: Beam chamber: -2"};

m_beamBN2Y2

Gaudi::Property<float> TileTBAANtuple::m_beamBN2Y2 {this, "BN2Y2", 0.2, "Params for Beam TDC: Beam chamber: -2"}

private

Definition at line 191 of file TileTBAANtuple.h.

{this, "BN2Y2", 0.2, "Params for Beam TDC: Beam chamber: -2"};

m_beamElemContainerKey

SG::ReadHandleKey<TileBeamElemContainer> TileTBAANtuple::m_beamElemContainerKey

private

Initial value:

{this,
        "TileBeamElemContainer", "TileBeamElemCnt", "Input Tile beam elements container"}

Definition at line 120 of file TileTBAANtuple.h.

                                                                 {this,
        "TileBeamElemContainer", "TileBeamElemCnt", "Input Tile beam elements container"};

m_beamFragList

Gaudi::Property<std::vector<std::string> > TileTBAANtuple::m_beamFragList {this, "beamFragList", {}, "List of beam frag IDs to store in the ntuple"}

private

Definition at line 186 of file TileTBAANtuple.h.

{this, "beamFragList", {}, "List of beam frag IDs to store in the ntuple"};

m_beamIdList

bool TileTBAANtuple::m_beamIdList[32]

private

Definition at line 607 of file TileTBAANtuple.h.

m_bsInput

Gaudi::Property<bool> TileTBAANtuple::m_bsInput {this, "BSInput", true, "Bytestream input"}

private

Definition at line 171 of file TileTBAANtuple.h.

{this, "BSInput", true, "Bytestream input"};

m_btdc

std::vector<std::vector<int> >* TileTBAANtuple::m_btdc {}

private

Definition at line 473 of file TileTBAANtuple.h.

{};

m_btdc1

std::array<int, 16> TileTBAANtuple::m_btdc1 {}

private

Definition at line 470 of file TileTBAANtuple.h.

{};

m_btdc2

std::array<int, 16> TileTBAANtuple::m_btdc2 {}

private

Definition at line 471 of file TileTBAANtuple.h.

{};

m_btdcNchMultiHit

int TileTBAANtuple::m_btdcNchMultiHit[2] {}

private

Definition at line 477 of file TileTBAANtuple.h.

{};

m_btdcNhit

int TileTBAANtuple::m_btdcNhit[16] {}

private

Definition at line 476 of file TileTBAANtuple.h.

{};

m_cabling

const TileCablingService* TileTBAANtuple::m_cabling {nullptr}

private

Definition at line 628 of file TileTBAANtuple.h.

{nullptr};

m_calibMode

Gaudi::Property<bool> TileTBAANtuple::m_calibMode {this, "CalibMode", false, "If data should be put in calib mode"}

private

Definition at line 168 of file TileTBAANtuple.h.

{this, "CalibMode", false, "If data should be put in calib mode"};

m_calibrateEnergy

Gaudi::Property<bool> TileTBAANtuple::m_calibrateEnergy {this, "CalibrateEnergy", true, "Calibrate energy"}

private

Definition at line 165 of file TileTBAANtuple.h.

{this, "CalibrateEnergy", true, "Calibrate energy"};

m_calibrateEnergyThisEvent

bool TileTBAANtuple::m_calibrateEnergyThisEvent {false}

private

Definition at line 620 of file TileTBAANtuple.h.

{false};

m_cellContainerKey

SG::ReadHandleKey<CaloCellContainer> TileTBAANtuple::m_cellContainerKey

private

Initial value:

{this,
        "CaloCellContainer", "AllCalo", "Input Calo cell container"}

Definition at line 153 of file TileTBAANtuple.h.

                                                         {this,
        "CaloCellContainer", "AllCalo", "Input Calo cell container"};

m_cher1

uint32_t TileTBAANtuple::m_cher1 {}

private

Definition at line 456 of file TileTBAANtuple.h.

{};

m_cher2

uint32_t TileTBAANtuple::m_cher2 {}

private

Definition at line 457 of file TileTBAANtuple.h.

{};

m_cher3

uint32_t TileTBAANtuple::m_cher3 {}

private

Definition at line 458 of file TileTBAANtuple.h.

{};

m_chi2cVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_chi2cVec

private

Definition at line 542 of file TileTBAANtuple.h.

m_chi2DspVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_chi2DspVec

private

Definition at line 549 of file TileTBAANtuple.h.

m_chi2FlatVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_chi2FlatVec

private

Definition at line 534 of file TileTBAANtuple.h.

m_chi2flxfitVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_chi2flxfitVec

private

Definition at line 553 of file TileTBAANtuple.h.

m_chi2flxoptVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_chi2flxoptVec

private

Definition at line 557 of file TileTBAANtuple.h.

m_chi2OptVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_chi2OptVec

private

Definition at line 546 of file TileTBAANtuple.h.

m_chi2Vec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_chi2Vec

private

Definition at line 538 of file TileTBAANtuple.h.

m_cispar

int TileTBAANtuple::m_cispar[16] {}

private

Definition at line 450 of file TileTBAANtuple.h.

{};

m_coincFlag1

int TileTBAANtuple::m_coincFlag1 {}

private

Definition at line 594 of file TileTBAANtuple.h.

{};

m_coincFlag2

int TileTBAANtuple::m_coincFlag2 {}

private

Definition at line 595 of file TileTBAANtuple.h.

{};

m_coincFlag3

int TileTBAANtuple::m_coincFlag3 {}

private

Definition at line 596 of file TileTBAANtuple.h.

{};

m_coincFlag4

int TileTBAANtuple::m_coincFlag4 {}

private

Definition at line 597 of file TileTBAANtuple.h.

{};

m_coincFlag5

int TileTBAANtuple::m_coincFlag5 {}

private

Definition at line 598 of file TileTBAANtuple.h.

{};

m_coincFlag6

int TileTBAANtuple::m_coincFlag6 {}

private

Definition at line 599 of file TileTBAANtuple.h.

{};

m_coincFlag7

int TileTBAANtuple::m_coincFlag7 {}

private

Definition at line 600 of file TileTBAANtuple.h.

{};

m_coincFlag8

int TileTBAANtuple::m_coincFlag8 {}

private

Definition at line 601 of file TileTBAANtuple.h.

{};

m_coincTrig1

std::array<unsigned int, 96> TileTBAANtuple::m_coincTrig1 {}

private

Definition at line 585 of file TileTBAANtuple.h.

{};

m_coincTrig2

std::array<unsigned int, 96> TileTBAANtuple::m_coincTrig2 {}

private

Definition at line 586 of file TileTBAANtuple.h.

{};

m_coincTrig3

std::array<unsigned int, 96> TileTBAANtuple::m_coincTrig3 {}

private

Definition at line 587 of file TileTBAANtuple.h.

{};

m_coincTrig4

std::array<unsigned int, 96> TileTBAANtuple::m_coincTrig4 {}

private

Definition at line 588 of file TileTBAANtuple.h.

{};

m_coincTrig5

std::array<unsigned int, 96> TileTBAANtuple::m_coincTrig5 {}

private

Definition at line 589 of file TileTBAANtuple.h.

{};

m_coincTrig6

std::array<unsigned int, 96> TileTBAANtuple::m_coincTrig6 {}

private

Definition at line 590 of file TileTBAANtuple.h.

{};

m_coincTrig7

std::array<unsigned int, 96> TileTBAANtuple::m_coincTrig7 {}

private

Definition at line 591 of file TileTBAANtuple.h.

{};

m_coincTrig8

std::array<unsigned int, 96> TileTBAANtuple::m_coincTrig8 {}

private

Definition at line 592 of file TileTBAANtuple.h.

{};

m_commonPU

int TileTBAANtuple::m_commonPU {}

private

Definition at line 441 of file TileTBAANtuple.h.

{};

m_completeNtuple

Gaudi::Property<bool> TileTBAANtuple::m_completeNtuple {this, "CompleteNtuple", true, "Complete the ntuple"}

private

Definition at line 170 of file TileTBAANtuple.h.

{this, "CompleteNtuple", true, "Complete the ntuple"};

m_detStore

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

privateinherited

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

m_digitsContainerFlxKey

SG::ReadHandleKey<TileDigitsContainer> TileTBAANtuple::m_digitsContainerFlxKey

private

Initial value:

{this,
        "TileDigitsContainerFlx", "", "Input Tile FELIX digits container"}

Definition at line 117 of file TileTBAANtuple.h.

                                                                {this,
        "TileDigitsContainerFlx", "", "Input Tile FELIX digits container"};

m_digitsContainerKey

SG::ReadHandleKey<TileDigitsContainer> TileTBAANtuple::m_digitsContainerKey

private

Initial value:

{this,
        "TileDigitsContainer", "TileDigitsCnt", "Input Tile digits container"}

Definition at line 114 of file TileTBAANtuple.h.

                                                             {this,
        "TileDigitsContainer", "TileDigitsCnt", "Input Tile digits container"};

m_DMUDataparityErrVec

std::vector<std::array<short, MAX_DMU> > TileTBAANtuple::m_DMUDataparityErrVec

private

Definition at line 523 of file TileTBAANtuple.h.

m_DMUDstrobeErrVec

std::vector<std::array<short, MAX_DMU> > TileTBAANtuple::m_DMUDstrobeErrVec

private

Definition at line 520 of file TileTBAANtuple.h.

m_DMUformatErrVec

std::vector<std::array<short, MAX_DMU> > TileTBAANtuple::m_DMUformatErrVec

private

Definition at line 517 of file TileTBAANtuple.h.

m_DMUheaderVec

std::vector<std::array<uint32_t, MAX_DMU> > TileTBAANtuple::m_DMUheaderVec

private

Definition at line 516 of file TileTBAANtuple.h.

m_DMUHeadparityErrVec

std::vector<std::array<short, MAX_DMU> > TileTBAANtuple::m_DMUHeadparityErrVec

private

Definition at line 522 of file TileTBAANtuple.h.

m_dmuMaskVec

std::vector<std::array<int, 2> > TileTBAANtuple::m_dmuMaskVec

private

Definition at line 524 of file TileTBAANtuple.h.

m_DMUmemoryErrVec

std::vector<std::array<short, MAX_DMU> > TileTBAANtuple::m_DMUmemoryErrVec

private

Definition at line 519 of file TileTBAANtuple.h.

m_DMUparityErrVec

std::vector<std::array<short, MAX_DMU> > TileTBAANtuple::m_DMUparityErrVec

private

Definition at line 518 of file TileTBAANtuple.h.

m_DMUSstrobeErrVec

std::vector<std::array<short, MAX_DMU> > TileTBAANtuple::m_DMUSstrobeErrVec

private

Definition at line 521 of file TileTBAANtuple.h.

m_dqStatusKey

SG::ReadHandleKey<TileDQstatus> TileTBAANtuple::m_dqStatusKey

private

Initial value:

{ this,
       "TileDQstatus", "TileDQstatus", "TileDQstatus key" }

Definition at line 162 of file TileTBAANtuple.h.

                                               { this,
       "TileDQstatus", "TileDQstatus", "TileDQstatus key" };

m_drawerFlxMap

std::map<unsigned int, unsigned int, std::less<unsigned int> > TileTBAANtuple::m_drawerFlxMap

private

Definition at line 604 of file TileTBAANtuple.h.

m_drawerList

Gaudi::Property<std::vector<std::string> > TileTBAANtuple::m_drawerList {this, "drawerList", {"-1"}, "List of frag IDs in correct order; Setup drawer list from data"}

private

Definition at line 184 of file TileTBAANtuple.h.

{this, "drawerList", {"-1"}, "List of frag IDs in correct order; Setup drawer list from data"};

m_drawerMap

std::map<unsigned int, unsigned int, std::less<unsigned int> > TileTBAANtuple::m_drawerMap

private

Definition at line 603 of file TileTBAANtuple.h.

m_drawerType

Gaudi::Property<std::vector<int> > TileTBAANtuple::m_drawerType {this, "drawerType", {}, "Type of every drawer 1-4: B+, B-, EB+, EB-; Take drawer type from Frag ID (doesn't work for 2003)"}

private

Definition at line 185 of file TileTBAANtuple.h.

{this, "drawerType", {}, "Type of every drawer 1-4: B+, B-, EB+, EB-; Take drawer type from Frag ID (doesn't work for 2003)"};

m_dspFlags

int TileTBAANtuple::m_dspFlags {}

private

Definition at line 356 of file TileTBAANtuple.h.

{};

m_dspRawChannelContainerKey

SG::ReadHandleKey<TileRawChannelContainer> TileTBAANtuple::m_dspRawChannelContainerKey

private

Initial value:

{this,
        "TileRawChannelContainerDsp", "", "Input Tile DSP raw channel container"}

Definition at line 132 of file TileTBAANtuple.h.

                                                                        {this,
        "TileRawChannelContainerDsp", "", "Input Tile DSP raw channel container"};

m_dspUnit

TileRawChannelUnit::UNIT TileTBAANtuple::m_dspUnit

private

Unit for TileRawChannels in DSP.

Definition at line 623 of file TileTBAANtuple.h.

m_ecal

std::array<float,8> TileTBAANtuple::m_ecal {}

private

Definition at line 372 of file TileTBAANtuple.h.

{};

m_eDspVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_eDspVec

private

Definition at line 547 of file TileTBAANtuple.h.

m_efitcVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_efitcVec

private

Definition at line 539 of file TileTBAANtuple.h.

m_efitVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_efitVec

private

Definition at line 535 of file TileTBAANtuple.h.

m_eflxfitVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_eflxfitVec

private

Definition at line 551 of file TileTBAANtuple.h.

m_eflxoptVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_eflxoptVec

private

Definition at line 555 of file TileTBAANtuple.h.

m_ehitCnt

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_ehitCnt

private

Definition at line 617 of file TileTBAANtuple.h.

m_ehitVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_ehitVec

private

Definition at line 615 of file TileTBAANtuple.h.

m_eneAdd

std::array<float, 16> TileTBAANtuple::m_eneAdd {}

private

Definition at line 446 of file TileTBAANtuple.h.

{};

m_eneVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_eneVec

private

Definition at line 531 of file TileTBAANtuple.h.

m_eOptVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_eOptVec

private

Definition at line 543 of file TileTBAANtuple.h.

m_eta

float TileTBAANtuple::m_eta {}

private

Definition at line 611 of file TileTBAANtuple.h.

{};

m_etaFileName

Gaudi::Property<std::string> TileTBAANtuple::m_etaFileName {this, "EtaFileName", "TileEtaCTB.txt", "File name with ETA, for CTB 2004 only"}

private

Definition at line 223 of file TileTBAANtuple.h.

{this, "EtaFileName", "TileEtaCTB.txt", "File name with ETA, for CTB 2004 only"};

m_evBCID

std::vector<int> TileTBAANtuple::m_evBCID

private

Definition at line 361 of file TileTBAANtuple.h.

m_eventsPerFile

Gaudi::Property<int> TileTBAANtuple::m_eventsPerFile {this, "EventsPerFile", 200000, "Number of events per file"}

private

Definition at line 178 of file TileTBAANtuple.h.

{this, "EventsPerFile", 200000, "Number of events per file"};

m_evt

int TileTBAANtuple::m_evt {}

private

Definition at line 354 of file TileTBAANtuple.h.

{};

m_evtflxVec

std::vector<int> TileTBAANtuple::m_evtflxVec

private

Definition at line 508 of file TileTBAANtuple.h.

m_evTime

int TileTBAANtuple::m_evTime {}

private

Definition at line 352 of file TileTBAANtuple.h.

{};

m_evtNr

int TileTBAANtuple::m_evtNr {-1}

private

Definition at line 349 of file TileTBAANtuple.h.

{-1};

m_evtStore

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

privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

m_evtVec

std::vector<int> TileTBAANtuple::m_evtVec

private

Definition at line 504 of file TileTBAANtuple.h.

m_evType

std::vector<int> TileTBAANtuple::m_evType

private

Definition at line 362 of file TileTBAANtuple.h.

m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects

privateinherited

Definition at line 114 of file AthAlgorithm.h.

m_ExtEne

std::array<float, 3> TileTBAANtuple::m_ExtEne {}

private

Definition at line 582 of file TileTBAANtuple.h.

{};

m_feCRCVec

std::vector<std::array<int, MAX_DMU> > TileTBAANtuple::m_feCRCVec

private

Definition at line 528 of file TileTBAANtuple.h.

m_finalUnit

Gaudi::Property<int> TileTBAANtuple::m_finalUnit {this, "OfflineUnits", TileRawChannelUnit::MegaElectronVolts, "Calibrate everything to this level"}

private

Definition at line 167 of file TileTBAANtuple.h.

{this, "OfflineUnits", TileRawChannelUnit::MegaElectronVolts, "Calibrate everything to this level"};

m_fitcRawChannelContainerKey

SG::ReadHandleKey<TileRawChannelContainer> TileTBAANtuple::m_fitcRawChannelContainerKey

private

Initial value:

{this,
        "TileRawChannelContainerFitCool", "", "Input Tile raw channel container reconstructed with Fit COOL method"}

Definition at line 135 of file TileTBAANtuple.h.

                                                                         {this,
        "TileRawChannelContainerFitCool", "", "Input Tile raw channel container reconstructed with Fit COOL method"};

m_fitRawChannelContainerKey

SG::ReadHandleKey<TileRawChannelContainer> TileTBAANtuple::m_fitRawChannelContainerKey

private

Initial value:

{this,
        "TileRawChannelContainerFit", "", "Input Tile raw channel container reconstructed with Fit method"}

Definition at line 126 of file TileTBAANtuple.h.

                                                                        {this,
        "TileRawChannelContainerFit", "", "Input Tile raw channel container reconstructed with Fit method"};

m_flatRawChannelContainerKey

SG::ReadHandleKey<TileRawChannelContainer> TileTBAANtuple::m_flatRawChannelContainerKey

private

Initial value:

{this,
        "TileRawChannelContainerFlat", "", "Input Tile raw channel container reconstructed with Flat method"}

Definition at line 123 of file TileTBAANtuple.h.

                                                                         {this,
        "TileRawChannelContainerFlat", "", "Input Tile raw channel container reconstructed with Flat method"};

m_flxFitRawChannelContainerKey

SG::ReadHandleKey<TileRawChannelContainer> TileTBAANtuple::m_flxFitRawChannelContainerKey

private

Initial value:

{this,
        "TileRawChannelContainerFitFlx", "", "Input Tile FELIX raw channel container reconstructed with Fit method"}

Definition at line 138 of file TileTBAANtuple.h.

                                                                           {this,
        "TileRawChannelContainerFitFlx", "", "Input Tile FELIX raw channel container reconstructed with Fit method"};

m_flxOptRawChannelContainerKey

SG::ReadHandleKey<TileRawChannelContainer> TileTBAANtuple::m_flxOptRawChannelContainerKey

private

Initial value:

{this,
        "TileRawChannelContainerOptFlx", "", "Input Tile FELIX raw channel container reconstructed with Opt2 method"}

Definition at line 141 of file TileTBAANtuple.h.

                                                                           {this,
        "TileRawChannelContainerOptFlx", "", "Input Tile FELIX raw channel container reconstructed with Opt2 method"};

m_frBCID

std::vector<int> TileTBAANtuple::m_frBCID

private

Definition at line 363 of file TileTBAANtuple.h.

m_gainflxVec

std::vector<std::array<int, MAX_CHAN> > TileTBAANtuple::m_gainflxVec

private

Definition at line 512 of file TileTBAANtuple.h.

m_gainVec

std::vector<std::array<int, MAX_CHAN> > TileTBAANtuple::m_gainVec

private

Definition at line 526 of file TileTBAANtuple.h.

m_GapEne

std::array<float, 3> TileTBAANtuple::m_GapEne {}

private

Definition at line 583 of file TileTBAANtuple.h.

{};

m_hitContainerKey

SG::ReadHandleKey<TileHitContainer> TileTBAANtuple::m_hitContainerKey

private

Initial value:

{this,
        "TileHitContainer", "TileHitCnt", "Input Tile hit container"}

Definition at line 147 of file TileTBAANtuple.h.

                                                       {this,
        "TileHitContainer", "TileHitCnt", "Input Tile hit container"};

m_hitVectorKey

SG::ReadHandleKey<TileHitVector> TileTBAANtuple::m_hitVectorKey

private

Initial value:

{this,
        "TileHitVector", "TileHitVec", "Input Tile hit vector"}

Definition at line 150 of file TileTBAANtuple.h.

                                                 {this,
        "TileHitVector", "TileHitVec", "Input Tile hit vector"};

m_l1ID

std::vector<int> TileTBAANtuple::m_l1ID

private

Definition at line 359 of file TileTBAANtuple.h.

m_l1Type

std::vector<int> TileTBAANtuple::m_l1Type

private

Definition at line 360 of file TileTBAANtuple.h.

m_LarEne

std::array<float, 4> TileTBAANtuple::m_LarEne {}

private

Definition at line 580 of file TileTBAANtuple.h.

{};

m_las0

float TileTBAANtuple::m_las0 {}

private

Definition at line 434 of file TileTBAANtuple.h.

{};

m_las1

float TileTBAANtuple::m_las1 {}

private

Definition at line 435 of file TileTBAANtuple.h.

{};

m_las2

float TileTBAANtuple::m_las2 {}

private

Definition at line 436 of file TileTBAANtuple.h.

{};

m_las3

float TileTBAANtuple::m_las3 {}

private

Definition at line 437 of file TileTBAANtuple.h.

{};

m_las_BCID

int TileTBAANtuple::m_las_BCID {}

private

Definition at line 378 of file TileTBAANtuple.h.

{};

m_las_D1_ADC

int TileTBAANtuple::m_las_D1_ADC {}

private

Definition at line 384 of file TileTBAANtuple.h.

{};

m_las_D1_Alpha

double TileTBAANtuple::m_las_D1_Alpha {}

private

Definition at line 399 of file TileTBAANtuple.h.

{};

m_las_D1_Alpha_RMS

double TileTBAANtuple::m_las_D1_Alpha_RMS {}

private

Definition at line 404 of file TileTBAANtuple.h.

{};

m_las_D1_AlphaPed

double TileTBAANtuple::m_las_D1_AlphaPed {}

private

Definition at line 409 of file TileTBAANtuple.h.

{};

m_las_D1_AlphaPed_RMS

double TileTBAANtuple::m_las_D1_AlphaPed_RMS {}

private

Definition at line 414 of file TileTBAANtuple.h.

{};

m_las_D1_Ped

double TileTBAANtuple::m_las_D1_Ped {}

private

Definition at line 389 of file TileTBAANtuple.h.

{};

m_las_D1_Ped_RMS

double TileTBAANtuple::m_las_D1_Ped_RMS {}

private

Definition at line 394 of file TileTBAANtuple.h.

{};

m_las_D2_ADC

int TileTBAANtuple::m_las_D2_ADC {}

private

Definition at line 385 of file TileTBAANtuple.h.

{};

m_las_D2_Alpha

double TileTBAANtuple::m_las_D2_Alpha {}

private

Definition at line 400 of file TileTBAANtuple.h.

{};

m_las_D2_Alpha_RMS

double TileTBAANtuple::m_las_D2_Alpha_RMS {}

private

Definition at line 405 of file TileTBAANtuple.h.

{};

m_las_D2_AlphaPed

double TileTBAANtuple::m_las_D2_AlphaPed {}

private

Definition at line 410 of file TileTBAANtuple.h.

{};

m_las_D2_AlphaPed_RMS

double TileTBAANtuple::m_las_D2_AlphaPed_RMS {}

private

Definition at line 415 of file TileTBAANtuple.h.

{};

m_las_D2_Ped

double TileTBAANtuple::m_las_D2_Ped {}

private

Definition at line 390 of file TileTBAANtuple.h.

{};

m_las_D2_Ped_RMS

double TileTBAANtuple::m_las_D2_Ped_RMS {}

private

Definition at line 395 of file TileTBAANtuple.h.

{};

m_las_D3_ADC

int TileTBAANtuple::m_las_D3_ADC {}

private

Definition at line 386 of file TileTBAANtuple.h.

{};

m_las_D3_Alpha

double TileTBAANtuple::m_las_D3_Alpha {}

private

Definition at line 401 of file TileTBAANtuple.h.

{};

m_las_D3_Alpha_RMS

double TileTBAANtuple::m_las_D3_Alpha_RMS {}

private

Definition at line 406 of file TileTBAANtuple.h.

{};

m_las_D3_AlphaPed

double TileTBAANtuple::m_las_D3_AlphaPed {}

private

Definition at line 411 of file TileTBAANtuple.h.

{};

m_las_D3_AlphaPed_RMS

double TileTBAANtuple::m_las_D3_AlphaPed_RMS {}

private

Definition at line 416 of file TileTBAANtuple.h.

{};

m_las_D3_Ped

double TileTBAANtuple::m_las_D3_Ped {}

private

Definition at line 391 of file TileTBAANtuple.h.

{};

m_las_D3_Ped_RMS

double TileTBAANtuple::m_las_D3_Ped_RMS {}

private

Definition at line 396 of file TileTBAANtuple.h.

{};

m_las_D4_ADC

int TileTBAANtuple::m_las_D4_ADC {}

private

Definition at line 387 of file TileTBAANtuple.h.

{};

m_las_D4_Alpha

double TileTBAANtuple::m_las_D4_Alpha {}

private

Definition at line 402 of file TileTBAANtuple.h.

{};

m_las_D4_Alpha_RMS

double TileTBAANtuple::m_las_D4_Alpha_RMS {}

private

Definition at line 407 of file TileTBAANtuple.h.

{};

m_las_D4_AlphaPed

double TileTBAANtuple::m_las_D4_AlphaPed {}

private

Definition at line 412 of file TileTBAANtuple.h.

{};

m_las_D4_AlphaPed_RMS

double TileTBAANtuple::m_las_D4_AlphaPed_RMS {}

private

Definition at line 417 of file TileTBAANtuple.h.

{};

m_las_D4_Ped

double TileTBAANtuple::m_las_D4_Ped {}

private

Definition at line 392 of file TileTBAANtuple.h.

{};

m_las_D4_Ped_RMS

double TileTBAANtuple::m_las_D4_Ped_RMS {}

private

Definition at line 397 of file TileTBAANtuple.h.

{};

m_las_Filt

int TileTBAANtuple::m_las_Filt {}

private

Definition at line 380 of file TileTBAANtuple.h.

{};

m_las_MeasAmp

double TileTBAANtuple::m_las_MeasAmp {}

private

Definition at line 382 of file TileTBAANtuple.h.

{};

m_las_PMT1_ADC

int TileTBAANtuple::m_las_PMT1_ADC {}

private

Definition at line 419 of file TileTBAANtuple.h.

{};

m_las_PMT1_Ped

double TileTBAANtuple::m_las_PMT1_Ped {}

private

Definition at line 425 of file TileTBAANtuple.h.

{};

m_las_PMT1_Ped_RMS

double TileTBAANtuple::m_las_PMT1_Ped_RMS {}

private

Definition at line 428 of file TileTBAANtuple.h.

{};

m_las_PMT1_TDC

int TileTBAANtuple::m_las_PMT1_TDC {}

private

Definition at line 422 of file TileTBAANtuple.h.

{};

m_las_PMT2_ADC

int TileTBAANtuple::m_las_PMT2_ADC {}

private

Definition at line 420 of file TileTBAANtuple.h.

{};

m_las_PMT2_Ped

double TileTBAANtuple::m_las_PMT2_Ped {}

private

Definition at line 426 of file TileTBAANtuple.h.

{};

m_las_PMT2_Ped_RMS

double TileTBAANtuple::m_las_PMT2_Ped_RMS {}

private

Definition at line 429 of file TileTBAANtuple.h.

{};

m_las_PMT2_TDC

int TileTBAANtuple::m_las_PMT2_TDC {}

private

Definition at line 423 of file TileTBAANtuple.h.

{};

m_las_ReqAmp

double TileTBAANtuple::m_las_ReqAmp {}

private

Definition at line 381 of file TileTBAANtuple.h.

{};

m_las_Temperature

double TileTBAANtuple::m_las_Temperature {}

private

Definition at line 431 of file TileTBAANtuple.h.

{};

m_laserObjectKey

SG::ReadHandleKey<TileLaserObject> TileTBAANtuple::m_laserObjectKey

private

Initial value:

{this,
        "TileLaserObj", "", "Input Tile laser object"}

Definition at line 144 of file TileTBAANtuple.h.

                                                     {this,
        "TileLaserObj", "", "Input Tile laser object"};

m_lasExtra

std::array<float, 4> TileTBAANtuple::m_lasExtra {}

private

Definition at line 438 of file TileTBAANtuple.h.

{};

m_lasFlag

int TileTBAANtuple::m_lasFlag {}

private

Definition at line 433 of file TileTBAANtuple.h.

{};

m_mdBcidflxVec

std::vector<std::array<int, MAX_MINIDRAWER> > TileTBAANtuple::m_mdBcidflxVec

private

Definition at line 570 of file TileTBAANtuple.h.

m_mdCapacitorflxVec

std::vector<std::array<int, MAX_MINIDRAWER> > TileTBAANtuple::m_mdCapacitorflxVec

private

Definition at line 578 of file TileTBAANtuple.h.

m_mdChargeflxVec

std::vector<std::array<int, MAX_MINIDRAWER> > TileTBAANtuple::m_mdChargeflxVec

private

Definition at line 576 of file TileTBAANtuple.h.

m_mdChargeTimeflxVec

std::vector<std::array<int, MAX_MINIDRAWER> > TileTBAANtuple::m_mdChargeTimeflxVec

private

Definition at line 577 of file TileTBAANtuple.h.

m_mdL1idflxVec

std::vector<std::array<int, MAX_MINIDRAWER> > TileTBAANtuple::m_mdL1idflxVec

private

Definition at line 569 of file TileTBAANtuple.h.

m_mdModuleflxVec

std::vector<std::array<int, MAX_MINIDRAWER> > TileTBAANtuple::m_mdModuleflxVec

private

Definition at line 571 of file TileTBAANtuple.h.

m_mdPedHiflxVec

std::vector<std::array<int, MAX_MINIDRAWER> > TileTBAANtuple::m_mdPedHiflxVec

private

Definition at line 574 of file TileTBAANtuple.h.

m_mdPedLoflxVec

std::vector<std::array<int, MAX_MINIDRAWER> > TileTBAANtuple::m_mdPedLoflxVec

private

Definition at line 573 of file TileTBAANtuple.h.

m_mdRunflxVec

std::vector<std::array<int, MAX_MINIDRAWER> > TileTBAANtuple::m_mdRunflxVec

private

Definition at line 575 of file TileTBAANtuple.h.

m_mdRunTypeflxVec

std::vector<std::array<int, MAX_MINIDRAWER> > TileTBAANtuple::m_mdRunTypeflxVec

private

Definition at line 572 of file TileTBAANtuple.h.

m_muBack

std::array<float,14> TileTBAANtuple::m_muBack {}

private

Definition at line 368 of file TileTBAANtuple.h.

{}; // MUON/MuBack

m_muBackHit

float TileTBAANtuple::m_muBackHit {}

private

Definition at line 366 of file TileTBAANtuple.h.

{};

m_muBackSum

float TileTBAANtuple::m_muBackSum {}

private

Definition at line 367 of file TileTBAANtuple.h.

{};

m_muCalib

std::array<float,2> TileTBAANtuple::m_muCalib {}

private

Definition at line 369 of file TileTBAANtuple.h.

{}; // MUON/MuCalib

m_muHalo

uint32_t TileTBAANtuple::m_muHalo {}

private

Definition at line 460 of file TileTBAANtuple.h.

{};

m_muTag

uint32_t TileTBAANtuple::m_muTag {}

private

Definition at line 459 of file TileTBAANtuple.h.

{};

m_muVeto

uint32_t TileTBAANtuple::m_muVeto {}

private

Definition at line 461 of file TileTBAANtuple.h.

{};

m_nDrawers

Gaudi::Property<unsigned int> TileTBAANtuple::m_nDrawers {this, "NDrawers", 6, "Number of drawers"}

private

Definition at line 175 of file TileTBAANtuple.h.

{this, "NDrawers", 6, "Number of drawers"};

m_nDrawersFlx

Gaudi::Property<unsigned int> TileTBAANtuple::m_nDrawersFlx {this, "NDrawersFelix", 0, "Number of drawers for FELIX"}

private

Definition at line 176 of file TileTBAANtuple.h.

{this, "NDrawersFelix", 0, "Number of drawers for FELIX"};

m_nSamples

Gaudi::Property<int> TileTBAANtuple::m_nSamples {this, "NSamples", NOT_SETUP, "Number of samples"}

private

Definition at line 173 of file TileTBAANtuple.h.

{this, "NSamples", NOT_SETUP, "Number of samples"};

m_nSamplesFlx

Gaudi::Property<int> TileTBAANtuple::m_nSamplesFlx {this, "NSamplesFelix", NOT_SETUP, "Number of samples for FELIX"}

private

Definition at line 174 of file TileTBAANtuple.h.

{this, "NSamplesFelix", NOT_SETUP, "Number of samples for FELIX"};

m_nSamplesFlxInDrawer

std::vector<int> TileTBAANtuple::m_nSamplesFlxInDrawer

private

Definition at line 633 of file TileTBAANtuple.h.

m_nSamplesFlxInDrawerMap

std::map<int, int> TileTBAANtuple::m_nSamplesFlxInDrawerMap

private

Definition at line 631 of file TileTBAANtuple.h.

m_nSamplesInDrawer

std::vector<int> TileTBAANtuple::m_nSamplesInDrawer

private

Definition at line 632 of file TileTBAANtuple.h.

m_nSamplesInDrawerMap

std::map<int, int> TileTBAANtuple::m_nSamplesInDrawerMap

private

Definition at line 630 of file TileTBAANtuple.h.

m_ntupleCreated

bool TileTBAANtuple::m_ntupleCreated {false}

private

Definition at line 346 of file TileTBAANtuple.h.

{false};

m_ntupleID

Gaudi::Property<std::string> TileTBAANtuple::m_ntupleID {this, "NTupleID", "h1000", "Name of the ntuple ID"}

private

Definition at line 181 of file TileTBAANtuple.h.

{this, "NTupleID", "h1000", "Name of the ntuple ID"};

m_ntuplePtr

TTree* TileTBAANtuple::m_ntuplePtr {nullptr}

private

Definition at line 345 of file TileTBAANtuple.h.

{nullptr};

m_optRawChannelContainerKey

SG::ReadHandleKey<TileRawChannelContainer> TileTBAANtuple::m_optRawChannelContainerKey

private

Initial value:

{this,
        "TileRawChannelContainerOpt", "TileRawChannelOpt2", "Input Tile raw channel container reconstructed with Opt2 method"}

Definition at line 129 of file TileTBAANtuple.h.

                                                                        {this,
        "TileRawChannelContainerOpt", "TileRawChannelOpt2", "Input Tile raw channel container reconstructed with Opt2 method"};

m_pedfitcVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_pedfitcVec

private

Definition at line 541 of file TileTBAANtuple.h.

m_pedfitVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_pedfitVec

private

Definition at line 537 of file TileTBAANtuple.h.

m_pedFlatVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_pedFlatVec

private

Definition at line 533 of file TileTBAANtuple.h.

m_pedflxfitVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_pedflxfitVec

private

Definition at line 554 of file TileTBAANtuple.h.

m_pedflxoptVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_pedflxoptVec

private

Definition at line 558 of file TileTBAANtuple.h.

m_pedOptVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_pedOptVec

private

Definition at line 545 of file TileTBAANtuple.h.

m_pmtOrder

Gaudi::Property<bool> TileTBAANtuple::m_pmtOrder {this, "PMTOrder", true, "Change channel ordering to pmt ordering in the ntuple"}

private

Definition at line 172 of file TileTBAANtuple.h.

{this, "PMTOrder", true, "Change channel ordering to pmt ordering in the ntuple"};

m_qdc

std::array<uint32_t,33> TileTBAANtuple::m_qdc {}

private

Definition at line 375 of file TileTBAANtuple.h.

{};

m_radius

Gaudi::Property<float> TileTBAANtuple::m_radius {this, "Radius", 2280.0, "Inner radius of calo, for CTB 2004 only"}

private

Definition at line 222 of file TileTBAANtuple.h.

{this, "Radius", 2280.0, "Inner radius of calo, for CTB 2004 only"};

m_rchUnit

TileRawChannelUnit::UNIT TileTBAANtuple::m_rchUnit

private

Unit for TileRawChannels (ADC, pCb, MeV).

Definition at line 622 of file TileTBAANtuple.h.

m_ROD_DMUBCIDVec

std::vector<std::array<short, MAX_DMU> > TileTBAANtuple::m_ROD_DMUBCIDVec

private

Definition at line 561 of file TileTBAANtuple.h.

m_ROD_DMUDataformatErrVec

std::vector<std::array<short, MAX_DMU> > TileTBAANtuple::m_ROD_DMUDataformatErrVec

private

Definition at line 566 of file TileTBAANtuple.h.

m_ROD_DMUDstrobeErrVec

std::vector<std::array<short, MAX_DMU> > TileTBAANtuple::m_ROD_DMUDstrobeErrVec

private

Definition at line 564 of file TileTBAANtuple.h.

m_ROD_DMUHeadformatErrVec

std::vector<std::array<short, MAX_DMU> > TileTBAANtuple::m_ROD_DMUHeadformatErrVec

private

Definition at line 565 of file TileTBAANtuple.h.

m_ROD_DMUMaskVec

std::vector<std::array<short, 2> > TileTBAANtuple::m_ROD_DMUMaskVec

private

Definition at line 567 of file TileTBAANtuple.h.

m_ROD_DMUmemoryErrVec

std::vector<std::array<short, MAX_DMU> > TileTBAANtuple::m_ROD_DMUmemoryErrVec

private

Definition at line 562 of file TileTBAANtuple.h.

m_ROD_DMUSstrobeErrVec

std::vector<std::array<short, MAX_DMU> > TileTBAANtuple::m_ROD_DMUSstrobeErrVec

private

Definition at line 563 of file TileTBAANtuple.h.

m_ROD_GlobalCRCVec

std::vector<short> TileTBAANtuple::m_ROD_GlobalCRCVec

private

Definition at line 560 of file TileTBAANtuple.h.

m_rodBCIDflxVec

std::vector<short> TileTBAANtuple::m_rodBCIDflxVec

private

Definition at line 509 of file TileTBAANtuple.h.

m_rodBCIDVec

std::vector<short> TileTBAANtuple::m_rodBCIDVec

private

Definition at line 505 of file TileTBAANtuple.h.

m_rodCRCVec

std::vector<std::array<int, MAX_DMU> > TileTBAANtuple::m_rodCRCVec

private

Definition at line 529 of file TileTBAANtuple.h.

m_rosName

Gaudi::Property<std::vector<std::string> > TileTBAANtuple::m_rosName {this, "rosName", {"B", "A", "C", "D", "E"}, "Name of arrays in ntuple for different ROSes"}

private

Definition at line 183 of file TileTBAANtuple.h.

{this, "rosName", {"B", "A", "C", "D", "E"}, "Name of arrays in ntuple for different ROSes"};

m_run

int TileTBAANtuple::m_run {}

private

Definition at line 353 of file TileTBAANtuple.h.

{};

m_runNumber

int TileTBAANtuple::m_runNumber {}

private

Definition at line 610 of file TileTBAANtuple.h.

{};

m_s1cou

uint32_t TileTBAANtuple::m_s1cou {}

private

Definition at line 453 of file TileTBAANtuple.h.

{};

m_s2cou

uint32_t TileTBAANtuple::m_s2cou {}

private

Definition at line 454 of file TileTBAANtuple.h.

{};

m_s2extra

int TileTBAANtuple::m_s2extra {}

private

Definition at line 463 of file TileTBAANtuple.h.

{};

m_s3cou

uint32_t TileTBAANtuple::m_s3cou {}

private

Definition at line 455 of file TileTBAANtuple.h.

{};

m_s3extra

int TileTBAANtuple::m_s3extra {}

private

Definition at line 464 of file TileTBAANtuple.h.

{};

m_sampleflxVec

std::vector<std::unique_ptr<int []> > TileTBAANtuple::m_sampleflxVec

private

Definition at line 513 of file TileTBAANtuple.h.

m_sampleVec

std::vector<std::unique_ptr<int []> > TileTBAANtuple::m_sampleVec

private

Definition at line 527 of file TileTBAANtuple.h.

m_samplingFractionKey

SG::ReadCondHandleKey<TileSamplingFraction> TileTBAANtuple::m_samplingFractionKey

private

Initial value:

{this,
        "TileSamplingFraction", "TileSamplingFraction", "Input Tile sampling fraction"}

Name of TileSamplingFraction in condition store.

Definition at line 111 of file TileTBAANtuple.h.

                                                                   {this,
        "TileSamplingFraction", "TileSamplingFraction", "Input Tile sampling fraction"};

m_saveFelixData

bool TileTBAANtuple::m_saveFelixData {false}

private

Definition at line 634 of file TileTBAANtuple.h.

{false};

m_sc1

int TileTBAANtuple::m_sc1 {}

private

Definition at line 466 of file TileTBAANtuple.h.

{};

m_sc2

int TileTBAANtuple::m_sc2 {}

private

Definition at line 467 of file TileTBAANtuple.h.

{};

m_scaler

std::array<int, 16> TileTBAANtuple::m_scaler {}

private

Definition at line 472 of file TileTBAANtuple.h.

{};

m_sizeflxVec

std::vector<short> TileTBAANtuple::m_sizeflxVec

private

Definition at line 510 of file TileTBAANtuple.h.

m_sizeVec

std::vector<short> TileTBAANtuple::m_sizeVec

private

Definition at line 506 of file TileTBAANtuple.h.

m_slinkCRCVec

std::vector<std::array<int, 2> > TileTBAANtuple::m_slinkCRCVec

private

Definition at line 525 of file TileTBAANtuple.h.

m_streamName

Gaudi::Property<std::string> TileTBAANtuple::m_streamName {this, "StreamName", "AANT", "Name of the output stream"}

private

Definition at line 180 of file TileTBAANtuple.h.

{this, "StreamName", "AANT", "Name of the output stream"};

m_TBperiod

Gaudi::Property<int> TileTBAANtuple::m_TBperiod {this, "TBperiod", 2016, "Tuned for 2016 testbeam by default"}

private

Definition at line 177 of file TileTBAANtuple.h.

{this, "TBperiod", 2016, "Tuned for 2016 testbeam by default"};

m_tDspVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_tDspVec

private

Definition at line 548 of file TileTBAANtuple.h.

m_tfitcVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_tfitcVec

private

Definition at line 540 of file TileTBAANtuple.h.

m_tfitVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_tfitVec

private

Definition at line 536 of file TileTBAANtuple.h.

m_tflxfitVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_tflxfitVec

private

Definition at line 552 of file TileTBAANtuple.h.

m_tflxoptVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_tflxoptVec

private

Definition at line 556 of file TileTBAANtuple.h.

m_theta

float TileTBAANtuple::m_theta {}

private

Definition at line 612 of file TileTBAANtuple.h.

{};

m_thistSvc

ServiceHandle<ITHistSvc> TileTBAANtuple::m_thistSvc

private

Definition at line 342 of file TileTBAANtuple.h.

m_thitCnt

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_thitCnt

private

Definition at line 618 of file TileTBAANtuple.h.

m_thitVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_thitVec

private

Definition at line 616 of file TileTBAANtuple.h.

m_tileHWID

const TileHWID* TileTBAANtuple::m_tileHWID {nullptr}

private

Definition at line 627 of file TileTBAANtuple.h.

{nullptr};

m_tileID

const TileID* TileTBAANtuple::m_tileID {nullptr}

private

Definition at line 626 of file TileTBAANtuple.h.

{nullptr};

m_tileToolEmscale

ToolHandle<TileCondToolEmscale> TileTBAANtuple::m_tileToolEmscale

private

Initial value:

{this,
        "TileCondToolEmscale", "TileCondToolEmscale", "Tile EMS conditions tool"}

Definition at line 156 of file TileTBAANtuple.h.

                                                     {this,
        "TileCondToolEmscale", "TileCondToolEmscale", "Tile EMS conditions tool"};

m_timeAdd

std::array<float, 16> TileTBAANtuple::m_timeAdd {}

private

Definition at line 447 of file TileTBAANtuple.h.

{};

m_timeVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_timeVec

private

Definition at line 532 of file TileTBAANtuple.h.

m_tjitter

int TileTBAANtuple::m_tjitter {}

private

Definition at line 474 of file TileTBAANtuple.h.

{};

m_tof

std::array<int, 16> TileTBAANtuple::m_tof {}

private

Definition at line 469 of file TileTBAANtuple.h.

{};

m_tOptVec

std::vector<std::array<float, MAX_CHAN> > TileTBAANtuple::m_tOptVec

private

Definition at line 544 of file TileTBAANtuple.h.

m_treeSize

Gaudi::Property<Long64_t> TileTBAANtuple::m_treeSize {this, "TreeSize", 16000000000LL, "Size of tree"}

private

Definition at line 179 of file TileTBAANtuple.h.

{this, "TreeSize", 16000000000LL, "Size of tree"};

m_trigType

int TileTBAANtuple::m_trigType {}

private

Definition at line 355 of file TileTBAANtuple.h.

{};

m_tscTOF

int TileTBAANtuple::m_tscTOF {}

private

Definition at line 475 of file TileTBAANtuple.h.

{};

m_unpackAdder

Gaudi::Property<bool> TileTBAANtuple::m_unpackAdder {this, "UnpackAdder", false, "Unpack adder"}

private

Definition at line 169 of file TileTBAANtuple.h.

{this, "UnpackAdder", false, "Unpack adder"};

m_useDspUnits

Gaudi::Property<bool> TileTBAANtuple::m_useDspUnits {this, "UseDspUnits", false, "Use DSP untis"}

private

Definition at line 166 of file TileTBAANtuple.h.

{this, "UseDspUnits", false, "Use DSP untis"};

m_varHandleArraysDeclared

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

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

m_xCha0

float TileTBAANtuple::m_xCha0 {}

private

Definition at line 483 of file TileTBAANtuple.h.

{};

m_xCha1

float TileTBAANtuple::m_xCha1 {}

private

Definition at line 486 of file TileTBAANtuple.h.

{};

m_xCha1_0

float TileTBAANtuple::m_xCha1_0 {}

private

Definition at line 490 of file TileTBAANtuple.h.

{};

m_xCha2

float TileTBAANtuple::m_xCha2 {}

private

Definition at line 488 of file TileTBAANtuple.h.

{};

m_xCha2_0

float TileTBAANtuple::m_xCha2_0 {}

private

Definition at line 492 of file TileTBAANtuple.h.

{};

m_xChN1

float TileTBAANtuple::m_xChN1 {}

private

Definition at line 481 of file TileTBAANtuple.h.

{};

m_xChN2

float TileTBAANtuple::m_xChN2 {}

private

Definition at line 479 of file TileTBAANtuple.h.

{};

m_xImp

float TileTBAANtuple::m_xImp {}

private

Definition at line 494 of file TileTBAANtuple.h.

{};

m_xImp_0

float TileTBAANtuple::m_xImp_0 {}

private

Definition at line 497 of file TileTBAANtuple.h.

{};

m_xImp_90

float TileTBAANtuple::m_xImp_90 {}

private

Definition at line 499 of file TileTBAANtuple.h.

{};

m_xImp_min90

float TileTBAANtuple::m_xImp_min90 {}

private

Definition at line 501 of file TileTBAANtuple.h.

{};

m_yCha0

float TileTBAANtuple::m_yCha0 {}

private

Definition at line 484 of file TileTBAANtuple.h.

{};

m_yCha1

float TileTBAANtuple::m_yCha1 {}

private

Definition at line 487 of file TileTBAANtuple.h.

{};

m_yCha1_0

float TileTBAANtuple::m_yCha1_0 {}

private

Definition at line 491 of file TileTBAANtuple.h.

{};

m_yCha2

float TileTBAANtuple::m_yCha2 {}

private

Definition at line 489 of file TileTBAANtuple.h.

{};

m_yCha2_0

float TileTBAANtuple::m_yCha2_0 {}

private

Definition at line 493 of file TileTBAANtuple.h.

{};

m_yChN1

float TileTBAANtuple::m_yChN1 {}

private

Definition at line 482 of file TileTBAANtuple.h.

{};

m_yChN2

float TileTBAANtuple::m_yChN2 {}

private

Definition at line 480 of file TileTBAANtuple.h.

{};

m_yImp

float TileTBAANtuple::m_yImp {}

private

Definition at line 495 of file TileTBAANtuple.h.

{};

m_yImp_0

float TileTBAANtuple::m_yImp_0 {}

private

Definition at line 498 of file TileTBAANtuple.h.

{};

m_yImp_90

float TileTBAANtuple::m_yImp_90 {}

private

Definition at line 500 of file TileTBAANtuple.h.

{};

m_yImp_min90

float TileTBAANtuple::m_yImp_min90 {}

private

Definition at line 502 of file TileTBAANtuple.h.

{};

---

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

• TileTBAANtuple.h
• TileTBAANtuple.cxx