CBNTAA_TBInfo Class Reference

#include <CBNTAA_TBInfo.h>

Inheritance diagram for CBNTAA_TBInfo:

Collaboration diagram for CBNTAA_TBInfo:

Public Member Functions
	CBNTAA_TBInfo (const std::string &name, ISvcLocator *pSvcLocator)
virtual	~CBNTAA_TBInfo ()
virtual StatusCode	CBNT_initialize () override
virtual StatusCode	CBNT_execute () override
virtual StatusCode	CBNT_finalize () override
virtual StatusCode	initialize () override
virtual StatusCode	execute () override
virtual StatusCode	finalize () override
virtual StatusCode	pre_execute ()
virtual StatusCode	sysInitialize () override
	Override sysInitialize. More...
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies. More...
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
virtual StatusCode	sysStart () override
	Handle START transition. More...
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
MsgStream &	msg (const MSG::Level lvl) const
bool	msgLvl (const MSG::Level lvl) const

Protected Types
enum	{ NOT_VALID = -999 }

Protected Member Functions
template<class T >
void	addBranch (const std::string &branchname, T &obj, const std::string &leaflist)
template<class T >
void	addBranch (const std::string &branchname, T *&obj)
template<class T >
void	addBranch (const std::string &branchname, T *&obj, int bufferSize, int splitLevel)
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...
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. More...

Protected Attributes
std::string	m_ntpath
std::string	m_ntTitle
TTree *	m_nt
MsgStream *	m_log
const LArEM_ID *	m_emId
const LArHEC_ID *	m_hecId
const LArFCAL_ID *	m_fcalId
const LArOnlineID *	m_onlineId

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
float	GetEnergy (float CurrB3, float CurrB4)
float	GetBdl3 (float current)
float	GetBdl4 (float current)
float	GetErrColl (float coll3_down, float coll3_up, float coll9_down, float coll9_up)
float	GetErrAbsE (float energy)
float	GetErrCurr (float currB3, float currB4)
float	SyncLoss (float energy, float currB3, float currB4)
float	SyncLossErr (float loss)
void	calculateAll (int runNumber, int eventNumber)
virtual StatusCode	CBNT_clear () override
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
ToolHandle< ICaloCoordinateTool >	m_coord {this, "CaloCoordinateTool", "TBCaloCoordinate"}
ToolHandle< TBCondRunParTool >	m_runpar {this, "CondRunParTool", "TBCondRunParTool"}
bool	m_DBRead
int	m_beamener
int	m_beamtype
bool	m_computeBeamEnergy
bool	m_dumpBeamLine
float	m_B8_Bdl
bool	m_is_VLE
int	m_prevrunNum
float	m_prevEnergy
float	m_preverrAbsEnergy
float	m_preverrCollimators
float	m_preverrCurrents
float	m_prevsycLoss
float	m_preverrSycLoss
float	m_prevB3
float	m_prevB4
float	m_prevB8
double	m_runEta_nt
long	m_runParticleType_nt
long	m_runEnergy_nt
std::vector< float > *	m_quad_file
std::vector< float > *	m_quad_equip
std::vector< float > *	m_bend_file
std::vector< float > *	m_bend_equip
std::vector< float > *	m_trim_file
std::vector< float > *	m_trim_equip
std::vector< float > *	m_coll_file
std::vector< float > *	m_coll_equip
float	m_energy
float	m_errAbsEnergy
float	m_errCollimators
float	m_errCurrents
float	m_sycLoss
float	m_errSycLoss
bool	m_initialized
DataObjIDColl	m_extendedExtraObjects
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 17 of file CBNTAA_TBInfo.h.

Member Typedef Documentation

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

protectedinherited

Enumerator
NOT_VALID

Definition at line 56 of file CBNT_TBRecBase.h.

{NOT_VALID = -999};

Constructor & Destructor Documentation

CBNTAA_TBInfo()

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

Definition at line 19 of file CBNTAA_TBInfo.cxx.

                                                                            :
CBNT_TBRecBase(name, pSvcLocator) 
{
  declareProperty("DBRead",     m_DBRead=false);
  declareProperty("BeamEnergy", m_beamener);
  declareProperty("BeamType",   m_beamtype);
  declareProperty("ComputeBeamEnergy", m_computeBeamEnergy=false);
  declareProperty("DumpBeamLine",      m_dumpBeamLine=true);
  m_quad_file = NULL;
  m_quad_equip = NULL;
  m_bend_file = NULL;
  m_bend_equip = NULL;
  m_trim_file = NULL;
  m_trim_equip = NULL;
  m_coll_file = NULL;
  m_coll_equip = NULL;  
  m_B8_Bdl = 0;
  m_energy = 0;
  m_errAbsEnergy = 0;
  m_errCollimators = 0;
  m_errCurrents = 0;
  m_errSycLoss = 0;
  m_is_VLE = 0;
  m_prevB3 = 0;
  m_prevB4 = 0;
  m_prevB8 = 0;
  m_prevEnergy = 0;
  m_preverrAbsEnergy = 0;
  m_preverrCollimators = 0;
  m_preverrCurrents = 0;
  m_preverrSycLoss = 0;
  m_prevrunNum = 0;
  m_prevsycLoss = 0;
  m_runEnergy_nt = 0;
  m_runEta_nt = 0;
  m_runParticleType_nt = 0;
  m_sycLoss = 0;
}

~CBNTAA_TBInfo()

CBNTAA_TBInfo::~CBNTAA_TBInfo ( )

virtual

Definition at line 58 of file CBNTAA_TBInfo.cxx.

{
   if(m_quad_file) delete m_quad_file;
   if(m_quad_equip) delete m_quad_equip;
   if(m_bend_file) delete m_bend_file;
   if(m_bend_equip) delete m_bend_equip;
   if(m_trim_file) delete m_trim_file;
   if(m_trim_equip) delete m_trim_equip;
   if(m_coll_file) delete m_coll_file;
   if(m_coll_equip) delete m_coll_equip;   
}

Member Function Documentation

addBranch() [1/3]

template<class T >

void CBNT_TBRecBase::addBranch ( const std::string &  branchname, T &  obj, const std::string &  leaflist  )

inlineprotectedinherited

Definition at line 44 of file CBNT_TBRecBase.h.

                                                                                                   {
     m_nt->Branch(branchname.c_str(), &obj, leaflist.c_str());
 }

addBranch() [2/3]

template<class T >

void CBNT_TBRecBase::addBranch ( const std::string &  branchname, T *&  obj  )

inlineprotectedinherited

Definition at line 47 of file CBNT_TBRecBase.h.

                                                                         {
     obj = new T();
     m_nt->Branch(branchname.c_str(), &obj);    
  }

addBranch() [3/3]

template<class T >

void CBNT_TBRecBase::addBranch ( const std::string &  branchname, T *&  obj, int  bufferSize, int  splitLevel  )

inlineprotectedinherited

Definition at line 51 of file CBNT_TBRecBase.h.

                                                                                                          {
     obj = new T();
     m_nt->Branch(branchname.c_str(), &obj, bufferSize, splitLevel);    
  }

calculateAll()

void CBNTAA_TBInfo::calculateAll ( int  runNumber, int  eventNumber  )

private

Definition at line 307 of file CBNTAA_TBInfo.cxx.

                                                              {
 
  int run = runNumber; 
  int event = eventNumber;
 
  ATH_MSG_VERBOSE( "run   " << run  );
 
  //check if the run is a VLE run then set the flag
  if (((run < 2101225) && (run >= 2101022)) || ((run < 2102165) && (run >= 2102003)))
    m_is_VLE=1;
 
  //starting from event=2 we dont have to recalculate the variables unless the run number changes
  if (event > 1)
    {      
      if (run == m_prevrunNum)
    {      
      m_energy = m_prevEnergy;
      m_errAbsEnergy = m_preverrAbsEnergy;
      m_errCollimators = m_preverrCollimators;
      m_errCurrents = m_preverrCurrents;
      m_sycLoss = m_prevsycLoss;
      m_errSycLoss = m_preverrSycLoss;
 
      return;
    }
      else
    m_prevrunNum=run;
    }
 
  m_runEnergy_nt= m_runpar->beam_energy();
  m_runParticleType_nt = m_runpar->beam_type();
  
  float bend_equip[9];
  float coll_equip[24];
  
  int nc_quad = 22;
  int nc_bend = 9;
  int nc_trim = 10;
  int nc_coll = 24;
 
  const char* folder;
  if (m_dumpBeamLine == true)
    {
//      const GenericDbTable* dcstbl_quad_file[22] = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
//      const GenericDbTable* dcstbl_quad_equip[22] = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
//      const GenericDbTable* dcstbl_bend_file[9] = {0,0,0,0,0,0,0,0,0};
//      const GenericDbTable* dcstbl_bend_equip[9] = {0,0,0,0,0,0,0,0,0};
//      const GenericDbTable* dcstbl_trim_file[10] = {0,0,0,0,0,0,0,0,0,0}; 
//      const GenericDbTable* dcstbl_trim_equip[10] = {0,0,0,0,0,0,0,0,0,0}; 
//      const GenericDbTable* dcstbl_coll_file[24] = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0}; 
//      const GenericDbTable* dcstbl_coll_equip[24] = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0}; 
      
//      int ncol_quad_file[22], nrow_quad_file[22];
//      int ncol_bend_file[9], nrow_bend_file[9];
//      int ncol_trim_file[10], nrow_trim_file[10];
//      int ncol_coll_file[24], nrow_coll_file[24];
//      int ncol_quad_equip[22], nrow_quad_equip[22]; 
//      int ncol_bend_equip[9], nrow_bend_equip[9];
//      int ncol_trim_equip[10], nrow_trim_equip[10];
//      int ncol_coll_equip[24], nrow_coll_equip[24];
      
//      std::vector<std::string> names_quad_file[22], rows_quad_file[22];
//      std::vector<std::string> names_bend_file[9], rows_bend_file[9];
//      std::vector<std::string> names_trim_file[10], rows_trim_file[10];
//      std::vector<std::string> names_coll_file[24], rows_coll_file[24];
//      std::vector<std::string> names_quad_equip[22], rows_quad_equip[22];
//      std::vector<std::string> names_bend_equip[9], rows_bend_equip[9];
//      std::vector<std::string> names_trim_equip[10], rows_trim_equip[10];
//      std::vector<std::string> names_coll_equip[24], rows_coll_equip[24];
      
      float quad_file[22];
      float bend_file[9];
      float trim_file[10];
      float coll_file[24];
      float quad_equip[22];
      float trim_equip[10];
      
      m_quad_file->reserve(nc_quad);
      m_quad_equip->reserve(nc_quad);
      for (int i=0; i<nc_quad; i++)
    {
      quad_file[i] = -1;
      quad_equip[i] = -1;
          unsigned int val1;
          unsigned int val2;
 
#if 0     
      if (run < 1000454)
        {
          folder = "/TILE/DCS/SYSTEM1/BEAM";
        }
      else
#endif
        {
          folder = "/TILE/DCS/SYSTEM1/BEAM";
        }
          val1 = i+6;
          val2 = i+48;
      
      if (StatusCode::SUCCESS!=m_runpar->getVal(folder,val1,quad_file[i])) {
            ATH_MSG_ERROR( "Cannot find val" << val1  );
          }
      if (StatusCode::SUCCESS!=m_runpar->getVal(folder,val2,quad_equip[i])){
        ATH_MSG_ERROR( "Cannot find val" << val2  );
          }
      
      (*m_quad_file)[i] = quad_file[i];
      (*m_quad_equip)[i] = quad_equip[i];
    }
    
 
      m_bend_file->reserve(nc_bend);
      m_bend_equip->reserve(nc_bend);
      for (int i=0; i<nc_bend; i++)
    {
      bend_file[i] = -1;
      bend_equip[i] = -1;
          unsigned int val1;
          unsigned int val2;
      
      if (run < 1000454)
        {
          folder = "/TILE/DCS/SYSTEM1/BEAM";
        }
      else
        {
          folder = "/TILE/DCS/TILE_LV_62/BEAM";
        }
          val1 = i+28;
          val2 = i+70;
      
      if (StatusCode::SUCCESS!=m_runpar->getVal(folder,val1,bend_file[i])) {
        ATH_MSG_ERROR( "Cannot find val" << val1  );
          }
      if (StatusCode::SUCCESS!=m_runpar->getVal(folder,val2,bend_equip[i])) {
        ATH_MSG_ERROR( "Cannot find val" << val2  );
          }
      
      (*m_bend_file)[i] = bend_file[i];
      (*m_bend_equip)[i] = bend_equip[i];
 
      //added on 12.03.2006 by H.Hakobyan for caching purposes
      if (!m_is_VLE)
        {
          //if the currents B3 and B4 are the same as in previous event then quit..
          if((m_prevB3 == bend_equip[2]) && (m_prevB4 == bend_equip[3]))
        return;
        }
      else
        {         
          //if the current B8 is the same as in previous event then quit..
          if(m_prevB8 == bend_equip[7])
        return; 
        }
        
    }
      
      m_trim_file->reserve(nc_trim);
      m_trim_equip->reserve(nc_trim);
      for (int i=0; i<nc_trim; i++)
    {
      trim_file[i] = -1;
      trim_equip[i] = -1;
          unsigned int val1;
          unsigned int val2;
      
      if (run < 1000454)
        {
          folder = "/TILE/DCS/SYSTEM1/BEAM";
        }
      else
        {
          folder = "/TILE/DCS/TILE_LV_62/BEAM";
        }
          val1 = i+37;
          val2 = i+79;
      
      if (StatusCode::SUCCESS!=m_runpar->getVal(folder,val1,trim_file[i])) {
        ATH_MSG_ERROR( "Cannot find val" << val1  );
          }
      if (StatusCode::SUCCESS!=m_runpar->getVal(folder,val2,trim_equip[i])) {
        ATH_MSG_ERROR( "Cannot find val" << val2  );
          }
      
      (*m_trim_file)[i] = trim_file[i];
      (*m_trim_equip)[i] = trim_equip[i];
    }
      
      
      m_coll_file->reserve(nc_coll);
      m_coll_equip->reserve(nc_coll);
      for (int i=0; i<nc_coll; i++)
    {
      coll_file[i] = -1;
      coll_equip[i] = -1;
          unsigned int val1;
          unsigned int val2;
      
      
      if (run < 1000454)
        {
          folder = "/TILE/DCS/SYSTEM1/BEAM";
        }
      else
        {
          folder = "/TILE/DCS/TILE_LV_62/BEAM";
        }
          val1 = i+90;
          val2 = i+115;
      
      if (StatusCode::SUCCESS!=m_runpar->getVal(folder,val1,coll_file[i])) {
        ATH_MSG_ERROR( "Cannot find val" << val1  );
          }
      if (StatusCode::SUCCESS!=m_runpar->getVal(folder,val2,coll_equip[i])) {
        ATH_MSG_ERROR( "Cannot find val" << val2  );
          }
      
      (*m_coll_file)[i] = coll_file[i];
      (*m_coll_equip)[i] = coll_equip[i];
    }
      
      
      for (int i=0; i<nc_quad; i++)
    {
      ATH_MSG_DEBUG( "Q" << i+1 << " from file   " << (*m_quad_file)[i]  );
    }
      for (int i=0; i<nc_bend; i++)
    {
      ATH_MSG_DEBUG( "B" << i+1 << " from file   " << (*m_bend_file)[i]  );
    }
      for (int i=0; i<nc_trim; i++)
    {
      ATH_MSG_DEBUG( "Trim" << i+1 << " from file   " << (*m_trim_file)[i]  );
    }
      for (int i=0; i<nc_coll; i++)
    {
      int j = i/2;
      if (i == 2*j)
        {
          ATH_MSG_DEBUG( "Coll" << j+1 << " from file   (up or right position)     " << (*m_coll_file)[i]  );
          ATH_MSG_DEBUG( "Coll" << j+1 << " from file   (down or left position)    " << (*m_coll_file)[i+1]  );
        }
    }
      
      
      for (int i=0; i<nc_quad; i++)
    {
      ATH_MSG_DEBUG( "Q" << i+1 << " from equipement   " << (*m_quad_equip)[i]  );
    }
      for (int i=0; i<nc_bend; i++)
    {
      ATH_MSG_DEBUG( "B" << i+1 << " from equipement   " << (*m_bend_equip)[i]  );
    }
      for (int i=0; i<nc_trim; i++)
    {
      ATH_MSG_DEBUG( "Trim" << i+1 << " from equipement   " << (*m_trim_equip)[i]  );
    }
      for (int i=0; i<nc_coll; i++)
    {
      int j = i/2;
      if (i == 2*j)
        {
          ATH_MSG_DEBUG( "Coll" << j+1 << " from equipement  (up or right position)   " << (*m_coll_equip)[i]  );
          ATH_MSG_DEBUG( "Coll" << j+1 << " from file   (down or left position)       " << (*m_coll_equip)[i+1]  );
        }
    }//for
    }//end of if (m_dumpBeamLine = true)
  else // if(m_dumpBeamLine == false)
    {
//      const GenericDbTable* dcstbl_bend_equip[9] = {0,0,0,0,0,0,0,0,0};
//      const GenericDbTable* dcstbl_coll_equip[24] = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0}; 
//      int ncol_bend_equip[9], nrow_bend_equip[9];
//      int ncol_coll_equip[9], nrow_coll_equip[9];
//      std::vector<std::string> names_bend_equip[9], rows_bend_equip[9];
//      std::vector<std::string> names_coll_equip[9], rows_coll_equip[9];
      
      for (int i=0; i<nc_bend; i++)
    {
      bend_equip[i] = -1;
      unsigned int val2;
 
      if (run < 1000454)
        {
          folder = "/TILE/DCS/SYSTEM1/BEAM";
        }
      else
        {
          folder = "/TILE/DCS/TILE_LV_62/BEAM";
        }
          val2 = i+70;
      if (StatusCode::SUCCESS!=m_runpar->getVal(folder,val2,bend_equip[i])) {
        ATH_MSG_ERROR( "Cannot find val" << val2  );
          }
    }
      
      
      for (int i=0; i<nc_coll; i++)
    {
      coll_equip[i] = -1;
      unsigned int val2;
      
      if (run < 1000454)
        {
          folder = "/TILE/DCS/SYSTEM1/BEAM";
        }
      else
        {
          folder = "/TILE/DCS/TILE_LV_62/BEAM";
        }
          val2 = i+115;
      if (StatusCode::SUCCESS!=m_runpar->getVal(folder,val2,coll_equip[i])) {
        ATH_MSG_ERROR( "Cannot find val" << val2  );
          }
    }
      
    }//end of if(m_dumpBeamLine = false)
  
  //check if the run is a VLE run then set the flag
  //if (((run < 2101225) && (run >= 2101022)) || ((run < 2102165) && (run >= 2102003)))
  //  m_is_VLE=1;
  
  // --- modifying the code introducing two cases - for VLE and HLE, with different beam energy calculations ---
  if(m_is_VLE)//run is a VLE run
    {   
      
      //ATH_MSG_INFO( "THIS is a VLE run..."  );
      //ATH_MSG_INFO( "B8 current=" << bend_equip[7]  );
      
      //a polynomial function giving the values of Bdl for any given B8 current value...MIGRAD method for 4th order polynome ...:)
      //    m_B8_Bdl = 0.00525593*bend_equip[7] + 0.00000410442*bend_equip[7]*bend_equip[7] - 
      //      0.00000000906592*bend_equip[7]*bend_equip[7]*bend_equip[7]+ 
      //      0.00000000000345019*bend_equip[7]*bend_equip[7]*bend_equip[7]*bend_equip[7];
      
      //a polynomial function giving the values of Bdl for any given B8 current value...MINOS method ...for 5th order polynome:)
      m_B8_Bdl = 0.00530982*bend_equip[7] + 0.00000328502*bend_equip[7]*bend_equip[7] - 
    0.00000000532397*bend_equip[7]*bend_equip[7]*bend_equip[7]- 
    0.00000000000276483*bend_equip[7]*bend_equip[7]*bend_equip[7]*bend_equip[7]+
    0.000000000000003368*bend_equip[7]*bend_equip[7]*bend_equip[7]*bend_equip[7]*bend_equip[7];
      
      //beam energy calculation basedon a famous formula..
      m_energy = 2.49826 * m_B8_Bdl;
      m_errAbsEnergy = GetErrAbsE(m_energy);
      
      ATH_MSG_DEBUG( "B8dl=" << m_B8_Bdl  );
      ATH_MSG_DEBUG( "m_energy=" << m_energy  );
      
      //C12 opening calculations... C12 is the onlycollimator which is responsible for VLE runs
      float coll12_opening = (fabs(coll_equip[11]) + fabs(coll_equip[23]))/2.0;
      m_errCollimators = sqrt(coll12_opening*coll12_opening)*(1./27);
      
      //calculation of errors of magnet currents..
      m_errCurrents = GetErrCurr(bend_equip[6], bend_equip[7]);
      m_sycLoss = SyncLoss(m_energy, bend_equip[6], bend_equip[7]);
      m_errSycLoss = SyncLossErr(m_sycLoss);
 
      //for caching purposes
      m_prevB8 = bend_equip[7];
    }
  else//run is a HLE run
    {
      m_energy = GetEnergy(bend_equip[2], bend_equip[3]);
      m_errAbsEnergy = GetErrAbsE(m_energy);
      m_errCollimators = GetErrColl(coll_equip[4], coll_equip[5], coll_equip[16], coll_equip[17]);
      m_errCurrents = GetErrCurr(bend_equip[2], bend_equip[3]);
      m_sycLoss = SyncLoss(m_energy, bend_equip[2], bend_equip[3]);
      m_errSycLoss = SyncLossErr(m_sycLoss);
 
      //for caching purposes
      m_prevB3 = bend_equip[2];
      m_prevB4 = bend_equip[3];
    }
 
  //assigning the previous vars for caching..
  m_prevEnergy = m_energy;
  m_preverrAbsEnergy = m_errAbsEnergy;
  m_preverrCollimators = m_errCollimators;
  m_preverrCurrents = m_errCurrents;
  m_prevsycLoss = m_sycLoss;
  m_preverrSycLoss = m_errSycLoss;
 
 
}

CBNT_clear()

StatusCode CBNTAA_TBInfo::CBNT_clear ( )

overrideprivatevirtual

Reimplemented from CBNT_TBRecBase.

Definition at line 292 of file CBNTAA_TBInfo.cxx.

{
   if(m_quad_file) m_quad_file->clear();
   if(m_quad_equip) m_quad_equip->clear();
   if(m_bend_file) m_bend_file->clear();
   if(m_bend_equip) m_bend_equip->clear();
   if(m_trim_file) m_trim_file->clear();
   if(m_trim_equip) m_trim_equip->clear();
   if(m_coll_file) m_coll_file->clear();
   if(m_coll_equip) m_coll_equip->clear();   
 
 
   return StatusCode::SUCCESS;
}

CBNT_execute()

StatusCode CBNTAA_TBInfo::CBNT_execute ( )

overridevirtual

Reimplemented from CBNT_TBRecBase.

Definition at line 124 of file CBNTAA_TBInfo.cxx.

                                      {
 
  const EventContext& ctx = Gaudi::Hive::currentContext();
 
  m_runEnergy_nt=m_beamener;
  m_runParticleType_nt=m_beamtype;
    
  if ( m_DBRead ) {
 
    int run = ctx.eventID().run_number(); 
    uint64_t event = ctx.eventID().event_number();
 
    ATH_MSG_VERBOSE( "run   " << run  );
 
    if((event > 0) && (m_prevrunNum==-1))
      m_prevrunNum=run;//initializing the prev run
 
    // This tool handles the conversion between local and ctb coordinates
    ATH_CHECK( m_coord.retrieve() );
    m_runEta_nt = m_coord->beam_local_eta();
    
    ATH_CHECK( m_runpar.retrieve() );
    
    if ( m_computeBeamEnergy ) {
 
    //below the calculation of all energies and their errors...
    calculateAll(run,event);
 
    ATH_MSG_DEBUG( "m_energy  " << m_energy  );
    ATH_MSG_DEBUG( "m_errAbsEnergy  " << m_errAbsEnergy  );
    ATH_MSG_DEBUG( "m_errCollimators  " << m_errCollimators  );
    ATH_MSG_DEBUG( "m_errCurrents  " << m_errCurrents  );
    ATH_MSG_DEBUG( "m_sycLoss  " << m_sycLoss  );
    ATH_MSG_DEBUG( "m_errSycLoss  " << m_errSycLoss  );
 
    }
 
  }//end of if (m_DBRead)
  
  return StatusCode::SUCCESS;
}

CBNT_finalize()

StatusCode CBNTAA_TBInfo::CBNT_finalize ( )

overridevirtual

Reimplemented from CBNT_TBRecBase.

Definition at line 285 of file CBNTAA_TBInfo.cxx.

                                        {
 
  return StatusCode::SUCCESS;
}

CBNT_initialize()

StatusCode CBNTAA_TBInfo::CBNT_initialize ( )

overridevirtual

Reimplemented from CBNT_TBRecBase.

Definition at line 70 of file CBNTAA_TBInfo.cxx.

                                         {
 
  StatusCode sc;
 
  addBranch("EVENT_RunEta", m_runEta_nt,"RunEta/D");
  addBranch("EVENT_RunPartType", m_runParticleType_nt,"RunPartType/l"); 
  addBranch("EVENT_RunEnergy", m_runEnergy_nt,"RunEnergy/l");
  addBranch("EVENT_Energy", m_energy,"Energy/f");
 
  if ( m_computeBeamEnergy ) {
  
     addBranch("EVENT_ErrAbsEnergy", m_errAbsEnergy,"ErrAbsEnergy/f");
     addBranch("EVENT_ErrCollimators", m_errCollimators,"ErrCollimators/f");
     addBranch("EVENT_ErrCurrents", m_errCurrents,"ErrCurrents/f");
     addBranch("EVENT_SyncLoss", m_sycLoss,"SyncLoss/f");
     addBranch("EVENT_ErrSyncLoss", m_errSycLoss,"ErrSyncLoss/f");
 
     if ( m_dumpBeamLine ) {
  
         addBranch("EVENT_Quad_file", m_quad_file); 
         addBranch("EVENT_Quad_equip", m_quad_equip);
     
         addBranch("EVENT_Bend_file", m_bend_file); 
         addBranch("EVENT_Bend_equip", m_bend_equip);
     
         addBranch("EVENT_Trim_file", m_trim_file); 
         addBranch("EVENT_Trim_equip", m_trim_equip);
      
         addBranch("EVENT_Coll_file", m_coll_file); 
         addBranch("EVENT_Coll_equip", m_coll_equip);
     }
 
  }
 
  //VLE flag set to false
  m_is_VLE=0;
  //Bdl value of B8 magnet is also set to 0 - the case for HLE runs..
  m_B8_Bdl=0.0;
  //set a previous vars to -1
  m_prevrunNum = -1;
  m_prevEnergy = -1.0;
  m_preverrAbsEnergy = -1.0;
  m_preverrCollimators = -1.0;
  m_preverrCurrents = -1.0;
  m_prevsycLoss = -1.0;
  m_preverrSycLoss = -1.0;
  m_prevB3 = -1.0;
  m_prevB4 = -1.0;
  m_prevB8 = -1.0;
 
  return StatusCode::SUCCESS;
}

declareGaudiProperty() [1/4]

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

inlineprivateinherited

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

Definition at line 170 of file AthCommonDataStore.h.

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

declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  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());
 
  }

declareGaudiProperty() [3/4]

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

inlineprivateinherited

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

Definition at line 184 of file AthCommonDataStore.h.

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

declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T > &  t, const SG::NotHandleType &    )

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(t);
  }

declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleBase &  hndl, const std::string &  doc, const SG::VarHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

  {
    this->declare(hndl.vhKey());
    hndl.vhKey().setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleKey &  hndl, const std::string &  doc, const SG::VarHandleKeyType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

  {
    this->declare(hndl);
    hndl.setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleKeyArray &  hndArr, const std::string &  doc, const SG::VarHandleKeyArrayType &    )

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

  {
 
    // std::ostringstream ost;
    // ost << Algorithm::name() << " VHKA declareProp: " << name 
    //     << " size: " << hndArr.keys().size() 
    //     << " mode: " << hndArr.mode() 
    //     << "  vhka size: " << m_vhka.size()
    //     << "\n";
    // debug() << ost.str() << endmsg;
 
    hndArr.setOwner(this);
    m_vhka.push_back(&hndArr);
 
    Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
    if (p != 0) {
      p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
    } else {
      ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                    << name);
    }
 
    return p;
 
  }

declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc, const SG::NotHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(name, property, doc);
  }

declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc = "none"  )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

  {
    typedef typename SG::HandleClassifier<T>::type htype;
    return declareProperty (name, property, doc, htype());
  }

declareProperty() [6/6]

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

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

detStore()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore() [1/2]

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

evtStore() [2/2]

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

inlineinherited

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

Definition at line 90 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode CBNT_TBRecBase::execute ( )

overridevirtualinherited

Definition at line 154 of file CBNT_TBRecBase.cxx.

                                   {
 
  StatusCode sc;
 
   // run pre-execution
  try {
            sc = this->pre_execute();
  }
        
  catch( const std::exception& Exception ) {
   *m_log << MSG::ERROR << " Standard exception " 
               << Exception.what() 
               << " caught from sub-algorithm::CBNT_pre_execute (). Disable !" << endmsg ;
    sc = this->setProperty(BooleanProperty( "Enable",false ) );
    return StatusCode::FAILURE;
   
  }
  catch (...) {
    *m_log << MSG::ERROR << " Unknown exception " 
         << " caught from sub-algorithm::CBNT_pre_execute (). Disable !" << endmsg ;
     sc = this->setProperty(BooleanProperty( "Enable",false ) );
     return StatusCode::FAILURE;
  }
 
  if (sc.isFailure()) {
    *m_log << MSG::ERROR << "CBNT_pre_execute() failed. Disable !" << endmsg;
     sc = this->setProperty(BooleanProperty( "Enable",false ) );
     return sc;
  }
           
  // now subalgorithm execution
  try {
     sc = this->CBNT_execute();
  }
 
  catch( const std::exception& Exception ) {
    *m_log << MSG::ERROR << " Standard exception " 
         << Exception.what() 
         << " caught from sub-algorithm::CBNT_execute () :" << endmsg ;
     return StatusCode::FAILURE;
  }
  catch (...) {
    *m_log << MSG::ERROR << " Unknown exception " 
         << " caught from sub-algorithm::CBNT_execute () :" << endmsg ;
     return StatusCode::FAILURE;
  }
 
  return sc;
}

extraDeps_update_handler()

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

protectedinherited

Add StoreName to extra input/output deps as needed.

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

extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 50 of file AthAlgorithm.cxx.

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

finalize()

StatusCode CBNT_TBRecBase::finalize ( )

overridevirtualinherited

Definition at line 205 of file CBNT_TBRecBase.cxx.

                                    {
 
  StatusCode sc;
 
  //now subalgorithm finalisation
  try {
    sc = this->CBNT_finalize();
  }
  catch( const std::exception& Exception ) {
    *m_log << MSG::ERROR << " Standard exception " 
         << Exception.what() 
         << " caught from sub-algorithm::CBNT_finalize () :" << endmsg ;
  }
  catch (...) {
     *m_log << MSG::ERROR << " Unknown exception " 
          << " caught from sub-algorithm::CBNT_finalize () :" << endmsg ;
  }
 
  return sc;
}

GetBdl3()

float CBNTAA_TBInfo::GetBdl3 ( float  current )

private

Definition at line 180 of file CBNTAA_TBInfo.cxx.

{
  float Bend3_low[2] = {0.0079347, 0.021813};
  float Bend3_high[4] = {0.437329, 0.0212849, 0.00000355516, -0.0000000147542};
 
  float Bdl3 = 0;
  
  current = fabs(current);
 
  if (current <= 824.12) Bdl3 = Bend3_low[0] + Bend3_low[1]*current;
  else Bdl3 = Bend3_high[0] + Bend3_high[1]*current + Bend3_high[2]*(current-824.12)*(current-824.12) + Bend3_high[3]*(current-824.12)*(current-824.12)*(current-824.12);
  
  return Bdl3;
}

GetBdl4()

float CBNTAA_TBInfo::GetBdl4 ( float  current )

private

Definition at line 196 of file CBNTAA_TBInfo.cxx.

{
  float Bend4_low[2] = {0.00786406, 0.021814};
  float Bend4_high[4] = {0.434258, 0.0212899, 0.00000356113, -0.0000000147379};
 
  float Bdl4 = 0;
 
  current = fabs(current);
 
  if (current <= 824.12) Bdl4 = Bend4_low[0] + Bend4_low[1]*current;
  else  Bdl4 = Bend4_high[0] + Bend4_high[1]*current + Bend4_high[2]*(current-824.12)*(current-824.12) + Bend4_high[3]*(current-824.12)*(current-824.12)*(current-824.12);
  
  return Bdl4;
}

GetEnergy()

float CBNTAA_TBInfo::GetEnergy ( float  CurrB3, float  CurrB4  )

private

Definition at line 169 of file CBNTAA_TBInfo.cxx.

{ 
  float  Bdltot = GetBdl3(currB3) + GetBdl4(currB4);
  
  float energy = (0.3/41.)*Bdltot * 1000;
  
  return energy;
}

GetErrAbsE()

float CBNTAA_TBInfo::GetErrAbsE ( float  energy )

private

Definition at line 222 of file CBNTAA_TBInfo.cxx.

{ 
  float fac = 25./energy;
  float err_abs = sqrt(fac*fac + 0.5*0.5);
 
  return err_abs;
}

GetErrColl()

float CBNTAA_TBInfo::GetErrColl ( float  coll3_down, float  coll3_up, float  coll9_down, float  coll9_up  )

private

Definition at line 212 of file CBNTAA_TBInfo.cxx.

{ 
  float coll3_opening = (fabs(coll3_down) + fabs(coll3_up))/2.0;
  float coll9_opening = (fabs(coll9_down) + fabs(coll9_up))/2.0;
 
  float errcoll = sqrt(coll3_opening*coll3_opening + coll9_opening*coll9_opening)*(1./27);
 
  return errcoll;
}

GetErrCurr()

float CBNTAA_TBInfo::GetErrCurr ( float  currB3, float  currB4  )

private

Definition at line 231 of file CBNTAA_TBInfo.cxx.

{
  float Bdl3 = GetBdl3(currB3);
  float Bdl4 = GetBdl4(currB4);
 
  float Bdl_tot = Bdl3 + Bdl4;
 
  float errBdl3 = (GetBdl3(currB3 + 0.1) - GetBdl3(currB3 - 0.1))/2;
  float errBdl4 = (GetBdl4(currB4 + 0.1) - GetBdl4(currB4 - 0.1))/2;
  const float fac3 = errBdl3/Bdl3;
  const float fac4 = errBdl4/Bdl4;
  float errBdl = sqrt(fac3*fac3 + fac4*fac4)*Bdl_tot;
 
  float resField = sqrt(0.0010*0.0010 + 0.0010*0.0010);
 
  errBdl = sqrt(errBdl*errBdl + resField*resField + Bdl_tot*0.0003 * Bdl_tot*0.0003)/Bdl_tot*100;
 
  return errBdl;
}

initialize()

StatusCode CBNT_TBRecBase::initialize ( )

overridevirtualinherited

Definition at line 22 of file CBNT_TBRecBase.cxx.

                                      {
  m_log=new MsgStream(msgSvc(), name());
  
  *m_log << MSG::DEBUG << "Initializing CBNT_TBRecBase base class" << endmsg;
  
  const CaloCell_ID* idHelper = nullptr;
  ATH_CHECK( detStore()->retrieve (idHelper, "CaloCell_ID") );
  m_emId=idHelper->em_idHelper();
  m_fcalId=idHelper->fcal_idHelper();
  m_hecId=idHelper->hec_idHelper();
 
 
  if (!m_emId) {
    (*m_log) << MSG::ERROR << "Could not access lar EM ID helper" << endmsg;
    return StatusCode::FAILURE;
  }
  if (!m_fcalId) {
    (*m_log) << MSG::ERROR << "Could not access lar FCAL ID helper" << endmsg;
    return StatusCode::FAILURE;
  }
  if (!m_hecId) {
    (*m_log) << MSG::ERROR << "Could not access lar HEC ID helper" << endmsg;
    return StatusCode::FAILURE;
  }
 
  StatusCode sc = detStore()->retrieve(m_onlineId, "LArOnlineID");
  if (sc.isFailure()) {
    (*m_log) << MSG::ERROR << "Could not get LArOnlineID helper !" << endmsg;
    return StatusCode::FAILURE;
  }
    else {
      (*m_log) << MSG::DEBUG << " Found the LArOnlineID helper. " << endmsg;
    }
 
 
  m_initialized=true;
  return StatusCode::SUCCESS;
}

inputHandles()

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

overridevirtualinherited

Return this algorithm's input handles.

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

msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msg() [2/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( const MSG::Level  lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                  {
    return this->msgStream(lvl);
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

outputHandles()

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

pre_execute()

StatusCode CBNT_TBRecBase::pre_execute ( )

virtualinherited

Definition at line 61 of file CBNT_TBRecBase.cxx.

                                       {
 
  StatusCode sc;
 
  if(!m_initialized) {
 
    if (m_ntpath.size()==0 || m_ntTitle.size()==0) {
      *m_log << MSG::ERROR << "Need to set variable 'm_ntpath' and 'm_ntTitle' in constructor of deriving class!" << endmsg;
      return StatusCode::FAILURE;
    }
 
    size_t i=m_ntpath.rfind('/');
    if (i==std::string::npos) {
      *m_log << MSG::ERROR << "Expected at least on '/' in path " << m_ntpath << endmsg;
      return StatusCode::FAILURE;
    }
    std::string basepath(m_ntpath.begin(),m_ntpath.begin()+i);
    //std::cout << "Basepath" << basepath << std::endl;
 
    // retrive pointer to THistSvc
    ServiceHandle<ITHistSvc> tHistSvc("THistSvc", name());
    ATH_CHECK( tHistSvc.retrieve() );
 
    // get TTree
    sc = tHistSvc->getTree(m_ntpath,m_nt);
    if (sc.isFailure()) {
       *m_log << MSG::ERROR << "Unable to retrieve TTree : " << m_ntpath << endmsg;
       return sc;
    }
    /*
    NTupleFilePtr file1(ntupleSvc(),basepath);
    if (!file1){
      (*m_log)  << MSG::ERROR << "Could not get NTupleFilePtr with path " << basepath << " failed" << endmsg;
      return StatusCode::FAILURE;
    }
    NTuplePtr nt(ntupleSvc(),m_ntpath);
    if (!nt) {
      nt=ntupleSvc()->book(m_ntpath,CLID_ColumnWiseTuple,m_ntTitle);
    }
    if (!nt){
      (*m_log)  << MSG::ERROR << "Booking of NTuple at "<< m_ntpath << " and name " << m_ntTitle << " failed" << endmsg;
      return StatusCode::FAILURE; 
    }
    m_nt=nt;
    */ 
    //std::cout << "Ntuple ptr:" << m_nt << std::endl;
 
    // subalgorithm initialisation
    try {
      sc = this->CBNT_initialize();
    }
    catch( const std::exception& Exception ) {
      *m_log << MSG::ERROR << " Standard exception " 
                  << Exception.what() 
                  << " caught from sub-algorithm::CBNT_initialize () :" << endmsg ;
    }
    catch (...) {
              *m_log << MSG::ERROR << " Unknown exception " 
         << " caught from sub-algorithm::CBNT_initialize () :" << endmsg ;
   }
 
   if (sc.isFailure())
     {
              *m_log << MSG::ERROR << "CBNT_initialize() failed" << endmsg;
               return sc;
     }
 
    m_initialized=true;
  }
  
  // clear data members
  try {
            sc = this->CBNT_clear();
  }
  catch( const std::exception& Exception ) {
   *m_log << MSG::ERROR << " Standard exception " 
        << Exception.what() 
        << " caught from sub-algorithm::CBNT_clear () :" << endmsg ;
  }
  catch (...) {
   *m_log << MSG::ERROR << " Unknown exception " 
               << " caught from sub-algorithm::CBNT_clear () :" << endmsg ;
  }
 
  if (sc.isFailure()) {
      *m_log << MSG::ERROR << "CBNT_clear() failed" << endmsg;
       return sc;
  }
 
  return StatusCode::SUCCESS ;
}

renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

SyncLoss()

float CBNTAA_TBInfo::SyncLoss ( float  energy, float  currB3, float  currB4  )

private

Definition at line 253 of file CBNTAA_TBInfo.cxx.

{
  float Lmag = 5.127;
 
  float Bdl_B3 = GetBdl3(currB3);
  float Bdl_B4 = GetBdl4(currB4);
 
  float B2dl_B3 = (Bdl_B3 * Bdl_B3)*(1./(Lmag*9));
  float B2dl_B4 = (Bdl_B4 * Bdl_B4)*(1./(Lmag*9));
  float B2dl_tot = B2dl_B3 + B2dl_B4;
  
  float B5suB34 = 0.19;
  float B6suB34 = 0.22;
  float B56lB34 = 0.25;
 
  float B2dl_B5 = (B2dl_B3 + B2dl_B4)*B5suB34*B5suB34*B56lB34;
  float B2dl_B6 = (B2dl_B3 + B2dl_B4)*B6suB34*B6suB34*B56lB34;
 
  float loss = 0.00000127*energy*energy*(3*B2dl_tot/2 + B2dl_B5 + B2dl_B6);
  return loss;
}

SyncLossErr()

float CBNTAA_TBInfo::SyncLossErr ( float  loss )

private

Definition at line 277 of file CBNTAA_TBInfo.cxx.

{
  float errloss = loss * 0.02;
  return errloss;
}

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< Algorithm > >.

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

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< Algorithm > >::sysStart ( )

overridevirtualinherited

Handle START transition.

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

updateVHKA()

void AthCommonDataStore< AthCommonMsg< 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_B8_Bdl

float CBNTAA_TBInfo::m_B8_Bdl

private

Definition at line 45 of file CBNTAA_TBInfo.h.

m_beamener

int CBNTAA_TBInfo::m_beamener

private

Definition at line 39 of file CBNTAA_TBInfo.h.

m_beamtype

int CBNTAA_TBInfo::m_beamtype

private

Definition at line 40 of file CBNTAA_TBInfo.h.

m_bend_equip

std::vector<float>* CBNTAA_TBInfo::m_bend_equip

private

Definition at line 72 of file CBNTAA_TBInfo.h.

m_bend_file

std::vector<float>* CBNTAA_TBInfo::m_bend_file

private

Definition at line 71 of file CBNTAA_TBInfo.h.

m_coll_equip

std::vector<float>* CBNTAA_TBInfo::m_coll_equip

private

Definition at line 78 of file CBNTAA_TBInfo.h.

m_coll_file

std::vector<float>* CBNTAA_TBInfo::m_coll_file

private

Definition at line 77 of file CBNTAA_TBInfo.h.

m_computeBeamEnergy

bool CBNTAA_TBInfo::m_computeBeamEnergy

private

Definition at line 41 of file CBNTAA_TBInfo.h.

m_coord

ToolHandle<ICaloCoordinateTool> CBNTAA_TBInfo::m_coord {this, "CaloCoordinateTool", "TBCaloCoordinate"}

private

Definition at line 36 of file CBNTAA_TBInfo.h.

m_DBRead

bool CBNTAA_TBInfo::m_DBRead

private

Definition at line 38 of file CBNTAA_TBInfo.h.

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_dumpBeamLine

bool CBNTAA_TBInfo::m_dumpBeamLine

private

Definition at line 42 of file CBNTAA_TBInfo.h.

m_emId

const LArEM_ID* CBNT_TBRecBase::m_emId

protectedinherited

Definition at line 65 of file CBNT_TBRecBase.h.

m_energy

float CBNTAA_TBInfo::m_energy

private

Definition at line 80 of file CBNTAA_TBInfo.h.

m_errAbsEnergy

float CBNTAA_TBInfo::m_errAbsEnergy

private

Definition at line 81 of file CBNTAA_TBInfo.h.

m_errCollimators

float CBNTAA_TBInfo::m_errCollimators

private

Definition at line 82 of file CBNTAA_TBInfo.h.

m_errCurrents

float CBNTAA_TBInfo::m_errCurrents

private

Definition at line 83 of file CBNTAA_TBInfo.h.

m_errSycLoss

float CBNTAA_TBInfo::m_errSycLoss

private

Definition at line 85 of file CBNTAA_TBInfo.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects

privateinherited

Definition at line 79 of file AthAlgorithm.h.

m_fcalId

const LArFCAL_ID* CBNT_TBRecBase::m_fcalId

protectedinherited

Definition at line 67 of file CBNT_TBRecBase.h.

m_hecId

const LArHEC_ID* CBNT_TBRecBase::m_hecId

protectedinherited

Definition at line 66 of file CBNT_TBRecBase.h.

m_initialized

bool CBNT_TBRecBase::m_initialized

privateinherited

Definition at line 40 of file CBNT_TBRecBase.h.

m_is_VLE

bool CBNTAA_TBInfo::m_is_VLE

private

Definition at line 46 of file CBNTAA_TBInfo.h.

m_log

MsgStream* CBNT_TBRecBase::m_log

protectedinherited

Definition at line 63 of file CBNT_TBRecBase.h.

m_nt

TTree* CBNT_TBRecBase::m_nt

protectedinherited

Definition at line 61 of file CBNT_TBRecBase.h.

m_ntpath

std::string CBNT_TBRecBase::m_ntpath

protectedinherited

Definition at line 58 of file CBNT_TBRecBase.h.

m_ntTitle

std::string CBNT_TBRecBase::m_ntTitle

protectedinherited

Definition at line 58 of file CBNT_TBRecBase.h.

m_onlineId

const LArOnlineID* CBNT_TBRecBase::m_onlineId

protectedinherited

Definition at line 68 of file CBNT_TBRecBase.h.

m_prevB3

float CBNTAA_TBInfo::m_prevB3

private

Definition at line 56 of file CBNTAA_TBInfo.h.

m_prevB4

float CBNTAA_TBInfo::m_prevB4

private

Definition at line 57 of file CBNTAA_TBInfo.h.

m_prevB8

float CBNTAA_TBInfo::m_prevB8

private

Definition at line 58 of file CBNTAA_TBInfo.h.

m_prevEnergy

float CBNTAA_TBInfo::m_prevEnergy

private

Definition at line 50 of file CBNTAA_TBInfo.h.

m_preverrAbsEnergy

float CBNTAA_TBInfo::m_preverrAbsEnergy

private

Definition at line 51 of file CBNTAA_TBInfo.h.

m_preverrCollimators

float CBNTAA_TBInfo::m_preverrCollimators

private

Definition at line 52 of file CBNTAA_TBInfo.h.

m_preverrCurrents

float CBNTAA_TBInfo::m_preverrCurrents

private

Definition at line 53 of file CBNTAA_TBInfo.h.

m_preverrSycLoss

float CBNTAA_TBInfo::m_preverrSycLoss

private

Definition at line 55 of file CBNTAA_TBInfo.h.

m_prevrunNum

int CBNTAA_TBInfo::m_prevrunNum

private

Definition at line 49 of file CBNTAA_TBInfo.h.

m_prevsycLoss

float CBNTAA_TBInfo::m_prevsycLoss

private

Definition at line 54 of file CBNTAA_TBInfo.h.

m_quad_equip

std::vector<float>* CBNTAA_TBInfo::m_quad_equip

private

Definition at line 69 of file CBNTAA_TBInfo.h.

m_quad_file

std::vector<float>* CBNTAA_TBInfo::m_quad_file

private

Definition at line 68 of file CBNTAA_TBInfo.h.

m_runEnergy_nt

long CBNTAA_TBInfo::m_runEnergy_nt

private

Definition at line 66 of file CBNTAA_TBInfo.h.

m_runEta_nt

double CBNTAA_TBInfo::m_runEta_nt

private

Definition at line 64 of file CBNTAA_TBInfo.h.

m_runpar

ToolHandle<TBCondRunParTool> CBNTAA_TBInfo::m_runpar {this, "CondRunParTool", "TBCondRunParTool"}

private

Definition at line 37 of file CBNTAA_TBInfo.h.

m_runParticleType_nt

long CBNTAA_TBInfo::m_runParticleType_nt

private

Definition at line 65 of file CBNTAA_TBInfo.h.

m_sycLoss

float CBNTAA_TBInfo::m_sycLoss

private

Definition at line 84 of file CBNTAA_TBInfo.h.

m_trim_equip

std::vector<float>* CBNTAA_TBInfo::m_trim_equip

private

Definition at line 75 of file CBNTAA_TBInfo.h.

m_trim_file

std::vector<float>* CBNTAA_TBInfo::m_trim_file

private

Definition at line 74 of file CBNTAA_TBInfo.h.

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

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

• CBNTAA_TBInfo.h
• CBNTAA_TBInfo.cxx