TileCellContainerCnv Class Reference

#include <TileCellContainerCnv.h>

Inheritance diagram for TileCellContainerCnv:

Collaboration diagram for TileCellContainerCnv:

Public Types
typedef T_AthenaPoolCustCnv< TRANS, PERS >	BaseType

Public Member Functions
	TileCellContainerCnv (ISvcLocator *svcloc)
virtual	~TileCellContainerCnv ()
virtual StatusCode	initialize () override
	initialization More...
virtual TileCellVec *	createPersistent (TileCellContainer *cont) override
virtual TileCellContainer *	createTransient () override
virtual StatusCode	finalize () override
virtual long	repSvcType () const override
virtual StatusCode	createObj (IOpaqueAddress *pAddr, DataObject *&pObj) override
	Create a transient object from a POOL persistent representation. More...
virtual StatusCode	createRep (DataObject *pObj, IOpaqueAddress *&pAddr) override
	Create a POOL persistent representation for a transient object. More...
virtual StatusCode	fillRepRefs (IOpaqueAddress *pAddr, DataObject *pObj) override
	Create a POOL persistent representation for a transient object. More...
void	addTrack (TIDA::Track *t)
const std::vector< TIDA::Track * > &	tracks () const
void	clear ()
void	selectTracks (const TrigInDetTrackCollection *trigtracks)
void	selectTracks (const Rec::TrackParticleContainer *trigtracks)
void	selectTracks (const Analysis::MuonContainer *muontracks)
bool	msgLvl (const MSG::Level lvl) const
	Test the output level. More...
MsgStream &	msg () const
	The standard message stream. More...
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream. More...
void	setLevel (MSG::Level lvl)
	Change the current logging level. More...

Static Public Member Functions
static const CLID &	classID ()
static long	storageType ()

Protected Types
typedef std::map< std::string, RootType >	ClassMap
typedef ClassMap::const_iterator	ClassMapIt
typedef std::mutex	CallMutex

Protected Member Functions
virtual PERS *	createPersistent (TRANS *obj)=0
virtual PERS *	createPersistentWithKey (TRANS *obj, const std::string &) override
	method to be implemented by the developer. More...
virtual TRANS *	createTransientWithKey (const std::string &) override
	method to be implemented by the developer. More...
template<class P >
P *	poolReadObject ()
	Read object of type P. More...
template<class P >
void	poolReadObject (TopLevelTPCnvBase &tlp_converter)
	Read object of type P (plus all extending objects) using the indicated top-level TP converter. More...
void	keepPoolObj (PERS *obj, const std::string &output)
	Remember the POOL object to be written out (will be deleted after commit) More...
virtual StatusCode	transToPers (TRANS *, PERS *&) override
	Obsolete methods replaced by createPersistent() and createTransient() obsolete. More...
virtual StatusCode	persToTrans (TRANS *&, PERS *) override
	obsolete More...
virtual StatusCode	DataObjectToPers (DataObject *pObj, IOpaqueAddress *&pAddr) override
	Convert an object into Persistent. More...
virtual StatusCode	DataObjectToPool (IOpaqueAddress *pAddr, DataObject *pObj) override
	Write an object into POOL. More...
virtual StatusCode	PoolToDataObject (DataObject *&pObj, const Token *token, const std::string &key) override
	Read an object from POOL. More...
virtual StatusCode	cleanUp (const std::string &output) override
	Callback from the CleanupSvc to delete persistent object in the local list. More...
template<class P >
Placement	setPlacementForP (P &p, const std::string &key, const std::string &output)
	Set POOL placement for type P. More...
template<class P >
StatusCode	objectToPool (P *pObj, Token *&token, const std::string &key, const std::string &output)
	Write an object into POOL returning its token. More...
template<class P >
StatusCode	poolToObject (const Token *&token, P *&pObj)
	Read an object from POOL. More...
virtual bool	compareClassGuid (const Guid &clid) const
	Compare POOL class GUID with the one of the object being read. More...
virtual void	setToken (const std::string &token)
	Set the token (in std::string representation) and classID for the object that will be read next. More...
virtual Placement	setPlacement (const std::string &key, const std::string &output)
	Set POOL placement. More...
virtual Placement	setPlacementWithType (const std::string &tname, const std::string &key, const std::string &output)
	Set POOL placement hint for a given type. More...
virtual const DataObject *	getDataObject () const
double	phiCorr (double phi)
void	ipCorr (double d0, double z0, double &d0c, double &z0c, double phi0, double eta, double pt)

Protected Attributes
std::map< std::string, std::vector< std::unique_ptr< PERS > > >	m_persObjLists
	Local cache for persistent objects created by this converter, grouped by processing slot These objects are deleted after a commit. More...
std::mutex	m_pListMutex
	protection mutex for m_persObjLists More...
Guid	m_classID
ServiceHandle< IAthenaPoolCnvSvc >	m_athenaPoolCnvSvc
RootType	m_classDesc
std::string	m_className
ClassMap	m_classDescs
std::string	m_containerPrefix
std::string	m_containerNameHint
std::string	m_branchNameHint
const DataObject *	m_dataObject
const Token *	m_i_poolToken
CallMutex	m_conv_mut
double	m_beamX
double	m_beamY
double	m_beamZ
std::vector< TIDA::Track * >	m_tracks

Private Member Functions
int	cell_index (int side, int phi, int eta) const
void	initIdToIndex ()
int	round32 (double x)
int	round16 (double x)
void	initMessaging () const
	Initialize our message level and MessageSvc. More...

Private Attributes
StoreGateSvc *	m_storeGate
const TileTBID *	m_tileTBID
const MbtsDetDescrManager *	m_mbtsMgr
int	m_version
Identifier	m_id [NCELLMBTS]
CaloDetDescrElement *	m_dde [NCELLMBTS]
int	m_gainIndex [17]
int	m_gain [8]
std::map< DataObject *, PERS * >	m_persMap
std::string	m_nm
	Message source name. More...
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels) More...
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer. More...
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level. More...
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging) More...

Static Private Attributes
static const int	NSIDE = 2
static const int	NPHI = 8
static const int	NETA = 2
static const int	NCELLMBTS = NSIDE*NPHI*NETA

Friends
class	CnvFactory< TileCellContainerCnv >

Detailed Description

Definition at line 18 of file TileCellContainerCnv.h.

Member Typedef Documentation

BaseType

template<class TRANS , class PERS >

typedef T_AthenaPoolCustCnv<TRANS, PERS> T_AthenaPoolCustomCnvWithKey< TRANS, PERS >::BaseType

inherited

Definition at line 41 of file T_AthenaPoolCustomCnv.h.

CallMutex

typedef std::mutex AthenaPoolConverter::CallMutex

protectedinherited

Definition at line 126 of file AthenaPoolConverter.h.

ClassMap

typedef std::map<std::string, RootType> AthenaPoolConverter::ClassMap

protectedinherited

Definition at line 114 of file AthenaPoolConverter.h.

ClassMapIt

typedef ClassMap::const_iterator AthenaPoolConverter::ClassMapIt

protectedinherited

Definition at line 115 of file AthenaPoolConverter.h.

Constructor & Destructor Documentation

TileCellContainerCnv()

TileCellContainerCnv::TileCellContainerCnv

(

ISvcLocator *

svcloc

)

Definition at line 18 of file TileCellContainerCnv.cxx.

  : TileCellContainerCnvBase::T_AthenaPoolCustomCnv(svcloc)
 // Must create DataVector that does NOT own elements
  , m_storeGate(0)
  , m_tileTBID(0)
  , m_mbtsMgr(0)
  , m_version(0)
  , m_id()
  , m_dde()
  , m_gainIndex()
  , m_gain()
{
}

~TileCellContainerCnv()

TileCellContainerCnv::~TileCellContainerCnv ( )

virtual

Definition at line 33 of file TileCellContainerCnv.cxx.

{
}

Member Function Documentation

addTrack()

void Converter::addTrack ( TIDA::Track *  t )

inlineinherited

Definition at line 45 of file Converter.h.

                              {
    m_tracks.push_back(t);
  }

cell_index()

int TileCellContainerCnv::cell_index ( int  side, int  phi, int  eta  ) const

inlineprivate

Definition at line 45 of file TileCellContainerCnv.h.

{ return (side*NPHI+phi)*NETA+eta; }

classID()

template<class TRANS , class PERS >

static const CLID& T_AthenaPoolCustCnv< TRANS, PERS >::classID ( )

staticinherited

Returns

class ID.

cleanUp()

template<class TRANS , class PERS >

virtual StatusCode T_AthenaPoolCustomCnvWithKey< TRANS, PERS >::cleanUp ( const std::string &  output )

overrideprotectedvirtualinherited

Callback from the CleanupSvc to delete persistent object in the local list.

Reimplemented from AthenaPoolConverter.

clear()

void Converter::clear ( )

inlineinherited

Definition at line 53 of file Converter.h.

{m_tracks.clear();} 

compareClassGuid()

template<class TRANS , class PERS >

virtual bool T_AthenaPoolCustCnv< TRANS, PERS >::compareClassGuid ( const Guid &  clid ) const

protectedvirtualinherited

Compare POOL class GUID with the one of the object being read.

Parameters

clid	[IN] GUID to compare

Returns

true if the type of the object about to be read matches the GUID

createObj()

StatusCode AthenaPoolConverter::createObj ( IOpaqueAddress *  pAddr, DataObject *&  pObj  )

overridevirtualinherited

Create a transient object from a POOL persistent representation.

Parameters

pAddr	[IN] IOpaqueAddress of POOL persistent representation.
pObj	[OUT] pointer to the transient object.

Reimplemented in T_AthenaPoolCoolMultChanCnv< COLL_T, ELEM_T, ELEM_P >, T_AthenaPoolCoolMultChanCnv< LArHVScaleCorrComplete, LArConditionsSubset< LArSingleFloatP >, LArSingleFloatSubset_p1 >, and T_AthenaPoolCoolMultChanCnv< TransientClass, LArConditionsSubset< LArSingleFloatP >, LArSingleFloatSubset_p1 >.

Definition at line 64 of file AthenaPoolConverter.cxx.

                                                                                  {
   TokenAddress* tokAddr = dynamic_cast<TokenAddress*>(pAddr);
   
   bool ownTokAddr = false;
   if (tokAddr == nullptr || tokAddr->getToken() == nullptr) {
      ownTokAddr = true;
      Token* token = new Token;
      token->fromString(*(pAddr->par()));
      GenericAddress* genAddr = dynamic_cast<GenericAddress*>(pAddr);
      if (not genAddr){
        ATH_MSG_ERROR("Dynamic cast failed in AthenaPoolConverter::createObj");
        //clean up
        delete token;
        return StatusCode::FAILURE;
      }
      tokAddr = new TokenAddress(*genAddr, token);
   }
   if( tokAddr->ipar()[0] > 0 and tokAddr->getToken()->auxString().empty() ) {
      char text[32];
      ::sprintf(text, "[CTXT=%08X]", static_cast<int>(*(pAddr->ipar())));
      tokAddr->getToken()->setAuxString(text);
   }
   ATH_MSG_VERBOSE("createObj: " << tokAddr->getToken()->toString() << ", CTX=" << tokAddr->ipar()[0]
                   << ", auxStr=" << tokAddr->getToken()->auxString() );
   std::lock_guard<CallMutex> lock(m_conv_mut);
   m_i_poolToken = tokAddr->getToken();
   try {
      std::string key = pAddr->par()[1];
      if (!PoolToDataObject(pObj, tokAddr->getToken(), key).isSuccess()) {
         ATH_MSG_ERROR("createObj PoolToDataObject() failed, Token = " << (tokAddr->getToken() ? tokAddr->getToken()->toString() : "NULL"));
         pObj = nullptr;
      }
   } catch (std::exception& e) {
      ATH_MSG_ERROR("createObj - caught exception: " << e.what());
      pObj = nullptr;
   }
   if (pObj == nullptr) {
      ATH_MSG_ERROR("createObj failed to get DataObject, Token = " << (tokAddr->getToken() ? tokAddr->getToken()->toString() : "NULL"));
   }
   if (ownTokAddr) {
      delete tokAddr; tokAddr = nullptr;
   }
   m_i_poolToken = nullptr;
   if (pObj == nullptr) {
      return(StatusCode::FAILURE);
   }
   return(StatusCode::SUCCESS);
}

createPersistent() [1/2]

TileCellVec * TileCellContainerCnv::createPersistent ( TileCellContainer *  cont )

overridevirtual

Definition at line 125 of file TileCellContainerCnv.cxx.

{
    // Convert every TileCell to 3 32-bit integers: ID,Ene, and (time,qual,qain)
 
    MsgStream log(msgSvc(),"TileCellContainerCnv" );
    MSG::Level logLevel = log.level();
    bool lDebug = (logLevel<=MSG::DEBUG);
    bool lVerbose = (logLevel<=MSG::VERBOSE);
 
    std::string name = m_storeGate->proxy(cont)->name();
    auto vecCell = std::make_unique<TileCellVec>();
    vecCell->reserve(NCELLMBTS);
 
    if (lDebug) log << MSG::DEBUG << "storing TileCells from " << name << " in POOL" << endmsg;
 
    vecCell->push_back(m_version);
    int nMBTSfound=0;
 
    std::vector<const TileCell *> allCells;
    
    switch (m_version) {
 
    case 1: // 3 words per cell, energy scale factor is 1000, time scale factor is 100
      for (const TileCell* cell : *cont) {
        if (lVerbose)
          log << MSG::VERBOSE 
              << "ene=" << cell->energy()
              << " time=" << cell->time() 
              << " qual=" << (int)cell->qual1()
              << " gain=" << (int)cell->gain()
              << endmsg;
        unsigned int id  = cell->ID().get_identifier32().get_compact();
        int          ene = round32(cell->energy() * 1000.);
        unsigned int tim = 0x8000 + round16(cell->time()*100.);
        unsigned int qua = std::max(0, std::min(0xFF, (int)cell->qual1()));
        unsigned int gai = std::max(0, std::min(0xFF,   0x80   + (int)(cell->gain())));
        unsigned int tqg = (tim<<16) | (qua<<8) | gai;
        vecCell->push_back(id);
        vecCell->push_back((unsigned int)ene);
        vecCell->push_back(tqg);
        if (lVerbose)
          log << MSG::VERBOSE << "packing cell in three words " 
              << MSG::hex << id << " " << ene << " " << tqg << MSG::dec << endmsg;
      }
      break;
 
    case 2: // 1 or 2 words for MBTS cells, 3 words for others, energy scale factor is 1000, time scale factor is 100
 
      // prepare vector with all cells first, expect at least 32 MBTS cells
      allCells.resize(NCELLMBTS);
      for (const TileCell* cell : *cont) {
        Identifier id = cell->ID();
        if (m_tileTBID->is_tiletb(id)) {
          int side = std::max(0,m_tileTBID->type(id));
          int phi  = m_tileTBID->module(id);
          int eta  = m_tileTBID->channel(id);
          int ind  = cell_index(side,phi,eta);
          if (eta<NETA && phi<NPHI && ind < NCELLMBTS) {
            allCells[ind] = cell;
            ++nMBTSfound;
          } else {
            allCells.push_back(cell);
          }
        } else {
          allCells.push_back(cell);
        }
      }
 
     if (nMBTSfound>0) {
 
      // save first 32 cells (MBTS) without identifiers, 2 words per cell, put zeros for empty cells
      // if MBTS energy is in pCb, then LSB corresponds to 1/12 ADC count of high gain
      for (int ind=0; ind<NCELLMBTS; ++ind) {
        int energy = 0;
        int time   = 0;
        int quality= 0;
        int gain   = m_gain[0]; // non-existing gain in CaloGain - to mark non-existing cells
        const TileCell* cell = allCells[ind];
        if (cell) {
          energy = round32(cell->energy() * 1000.);
          time   = round16(cell->time() * 100.);
          quality= cell->qual1();
          gain   = cell->gain();
 
          if (lVerbose)
            log << MSG::VERBOSE 
                << "ind="  << ind
                << " ene=" << cell->energy()
                << " time=" << cell->time() 
                << " qual=" << (int)cell->qual1()
               << " gain=" << (int)cell->gain()
                << endmsg;
        } 
        else {
          if (lVerbose)
            log << MSG::VERBOSE 
                << "ind="  << ind
                << " create MBTS cell with zero energy"
                << endmsg;
        }
        
        // put correct MBTS cells in one word
        if (time == 0 && // expect time to be equal to zero
            -0x10000 < energy && energy < 0xEFFFF && // expect energy within (-65,980) pCb 
            -17 < gain && gain < 0 ) { // expext only gains in TileCal range
 
          unsigned int ene = energy+0x10000;  // shift by 65 pCb (65*10^3 because of scaling)
          unsigned int qua = std::max(0, std::min(0xFF, quality)); // 8 bits for quality
          unsigned int gai = m_gainIndex[-gain];
          unsigned int gqe = (gai << 28) | (qua<<20) | ene; // upper most bit is always 1 here
          vecCell->push_back(gqe);
 
          if (lVerbose)
            log << MSG::VERBOSE << "packing cell " << ind << " in one word "
                << MSG::hex << gqe << MSG::dec << endmsg;
 
        } else { // cells with time, use 2 words for channel
                 // but make sure that upper most bit in energy word is zero
 
          unsigned int ene = std::max(0, std::min(0x7FFFFFFF, 0x40000000 + energy));
          unsigned int tim = std::max(0, std::min(0xFFFF, 0x8000 + time));
          unsigned int qua = std::max(0, std::min(0xFF, quality)); // 8 bits for quality
          unsigned int gai = std::max(0, std::min(0xFF,   0x80   + gain));
          unsigned int tqg = (tim<<16) | (qua<<8) | gai;
          vecCell->push_back(ene);
          vecCell->push_back(tqg);
 
          if (lVerbose)
            log << MSG::VERBOSE << "packing cell " << ind << " in two words "
                << MSG::hex << ene << " " << tqg << MSG::dec << endmsg;
        }
      }
 
     } else {
 
       (*vecCell)[0] = 1; // no MBTS found - use version 1 for packing
     }
 
      // keep all other cells (if any) with identifiers, 3 words per cell
      for (unsigned int ind=NCELLMBTS; ind<allCells.size(); ++ind) {
        
        const TileCell* cell = allCells[ind];
 
        if (lVerbose)
          log << MSG::VERBOSE 
              << "ind="  << ind
              << " ene=" << cell->energy()
              << " time=" << cell->time() 
              << " qual=" << (int)cell->qual1()
              << " gain=" << (int)cell->gain()
              << endmsg;
 
        unsigned int id  = cell->ID().get_identifier32().get_compact();
        int          ene = round32(cell->energy() * 1000.);
        unsigned int tim = 0x8000 + round16(cell->time()*100.);
        unsigned int qua = std::max(0, std::min(0xFF, (int)cell->qual1()));
        unsigned int gai = std::max(0, std::min(0xFF,   0x80   + (int)(cell->gain())));
        unsigned int tqg = (tim<<16) | (qua<<8) | gai;
        vecCell->push_back(id);
        vecCell->push_back((unsigned int)ene);
        vecCell->push_back(tqg);
 
        if (lVerbose)
          log << MSG::VERBOSE << "packing cell " << ind << " in three words "
              << MSG::hex << id << " " << ene << " " << tqg << MSG::dec << endmsg;
      }
      break;
 
    default:
 
      log << MSG::ERROR << "Unknown version of TileCellVec, ver="<<m_version << endmsg;
 
    }
 
    if (lDebug) log << MSG::DEBUG << "Storing data vector of size " << vecCell->size() << " with version " << vecCell->front() << endmsg;
 
    return vecCell.release();
}

createPersistent() [2/2]

template<class TRANS , class PERS >

virtual PERS* T_AthenaPoolCustomCnv< TRANS, PERS >::createPersistent ( TRANS *  obj )

protectedpure virtualinherited

Implemented in ZdcRawChannelCollectionCnv, TruthParticleContainerCnv, TruthEtIsolationsContainerCnv, TRTUncompressedHitCollectionCnv, TrigTauTracksInfoCnv, TrigTauClusterDetailsCnv, TrigTauClusterCnv, TrigTauCnv, TrigRoiDescriptorCollectionCnv, TrigRoiDescriptorCnv, TrigRNNOutputCnv, TrigOperationalInfoCollectionCnv, TrigOperationalInfoCnv, TrigMonEventCollectionCnv, TrigMonEventCnv, TrigMonConfigCollectionCnv, TrigMonConfigCnv, TrigEMClusterCnv, TrigCompositeCnv, T_AthenaPoolViewVectorCnv< DV >, T_AthenaPoolLegacyCnv< TRANS, PERS, LEGACY >::LegacyCnv, T_AthenaPoolLegacyCnv< TRANS, PERS, LEGACY >, TrackRecordCollectionCnv, TileTTL1ContainerCnv, TileRawChannelContainerCnv, TileMuonReceiverContainerCnv, TileMuContainerCnv, TileL2ContainerCnv, TileHitVectorCnv, TileHitContainerCnv, TileDigitsContainerCnv, TileCosmicMuonContainerCnv, TileBeamElemContainerCnv, LArRawEventContainerCnv< T, P, C >, MuonRdoContainerTPCnv< T, P, CONV >, SimulationHitCollectionCnv, SiHitCollectionCnv, RingerRingsCnv, RecoTimingObjCnv, MuonSpShowerContainerCnv, PhotonContainerCnv, MuonCaloEnergyContainerCnv, MissingEtTruthCnv, MissingEtCaloCnv, MissingETCnv, McEventCollectionCnv, LArTTL1ContainerCnv, LArHitFloatContainerCnv, LArHitContainerCnv, IParticleLinkContainerCnv, InDetLowBetaCandidateCnv, INav4MomLinkContainerCnv, INav4MomAssocsCnv, ElectronContainerCnv, egDetailContainerCnv, egammaContainerCnv, eflowObjectContainerCnv, CaloRingsContainerCnv, CaloCalibrationHitContainerCnv, BCM_RDO_ContainerCnv, AthExFatObjectCnv, AthExElephantinoCnv, AthExDecayCnv, MuonContainerCnv, TauJetContainerCnv, TauDetailsContainerCnv, ALFA_ODHitCollectionCnv, ALFA_ODDigitCollectionCnv, ALFA_LocRecODEvCollectionCnv, ALFA_LocRecEvCollectionCnv, ALFA_LocRecCorrODEvCollectionCnv, ALFA_LocRecCorrEvCollectionCnv, ALFA_HitCollectionCnv, ALFA_GloRecEvCollectionCnv, ALFA_DigitCollectionCnv, ALFA_CLinkEventCnv, AFP_TDSimHitCollectionCnv, AFP_TDLocRecoEvCollectionCnv, AFP_TDDigiCollectionCnv, AFP_SIDSimHitCollectionCnv, AFP_SIDLocRecoEvCollectionCnv, and AFP_SiDigiCollectionCnv.

createPersistentWithKey()

template<class TRANS , class PERS >

virtual PERS* T_AthenaPoolCustomCnv< TRANS, PERS >::createPersistentWithKey ( TRANS *  obj, const std::string &  key  )

overrideprotectedvirtualinherited

method to be implemented by the developer.

It should create the persistent representation of the object, using the default top-level TP converter.

Parameters

obj	[IN] transient object

Returns

the created persistent representation (by pointer)

Implements T_AthenaPoolCustomCnvWithKey< TRANS, PERS >.

createRep()

StatusCode AthenaPoolConverter::createRep ( DataObject *  pObj, IOpaqueAddress *&  pAddr  )

overridevirtualinherited

Create a POOL persistent representation for a transient object.

Parameters

pObj	[IN] pointer to the transient object.
pAddr	[OUT] IOpaqueAddress of POOL persistent representation.

Reimplemented in T_AthenaPoolCoolMultChanCnv< COLL_T, ELEM_T, ELEM_P >, T_AthenaPoolCoolMultChanCnv< LArHVScaleCorrComplete, LArConditionsSubset< LArSingleFloatP >, LArSingleFloatSubset_p1 >, and T_AthenaPoolCoolMultChanCnv< TransientClass, LArConditionsSubset< LArSingleFloatP >, LArSingleFloatSubset_p1 >.

Definition at line 113 of file AthenaPoolConverter.cxx.

                                                                                  {
   const SG::DataProxy* proxy = dynamic_cast<SG::DataProxy*>(pObj->registry());
   if (proxy == nullptr) {
      ATH_MSG_ERROR("AthenaPoolConverter CreateRep failed to cast DataProxy, key = " << pObj->name());
      return(StatusCode::FAILURE);
   }
   try {
      std::lock_guard<CallMutex> lock(m_conv_mut);
      if (!DataObjectToPers(pObj, pAddr).isSuccess()) {
         ATH_MSG_ERROR("CreateRep failed, key = " << pObj->name());
         return(StatusCode::FAILURE);
      }
   } catch (std::exception& e) {
      ATH_MSG_ERROR("createRep - caught exception: " << e.what());
      return(StatusCode::FAILURE);
   }
   const CLID clid = proxy->clID();
   if (pAddr == nullptr) {
      // Create a IOpaqueAddress for this object.
      pAddr = new TokenAddress(this->storageType(), clid, "", "", 0, 0);
   } else {
      GenericAddress* gAddr = dynamic_cast<GenericAddress*>(pAddr);
      if (gAddr != nullptr) {
         gAddr->setSvcType(this->storageType());
      }
   }
   return(StatusCode::SUCCESS);
}

createTransient()

TileCellContainer * TileCellContainerCnv::createTransient ( )

overridevirtual

Implements T_AthenaPoolCustomCnv< TRANS, PERS >.

Definition at line 304 of file TileCellContainerCnv.cxx.

{
    // Fill TileCellContainer from vector, creating cells from 3 integers 
 
    MsgStream log(msgSvc(), "TileCellContainerCnv");
    MSG::Level logLevel = log.level();
    bool lDebug = (logLevel<=MSG::DEBUG);
    bool lVerbose = (logLevel<=MSG::VERBOSE);
 
    std::unique_ptr<TileCellVec> vec(this->poolReadObject<TileCellVec>());
 
    if (lDebug) log << MSG::DEBUG << "Read TileCell Vec, size " << vec->size() << endmsg;
 
    // create the TileCellContainer
    auto cont = std::make_unique<TileCellContainer>();
 
    TileCellVec::const_iterator it   = vec->begin();
    TileCellVec::const_iterator last = vec->end();
 
    unsigned int version = *it++;
    int iCell = 0;
    
    switch (version) {
    case 1:
 
      for (; it != last; ) {
 
        Identifier id(Identifier32(*it++));
        int ene = (int)(*it++);
        unsigned int tqg = *it++;
      
        float ener = ene*1e-3;
        float time = ((int)(tqg>>16) - 0x8000 ) * 0.01;
        uint16_t qual = ((tqg>>8) & 0xFF);
        uint16_t qbit = TileCell::MASK_CMPC | TileCell::MASK_TIME;
        int   gain = (int)(tqg & 0xFF) - 0x80;
 
        if (lVerbose) {
          log << MSG::VERBOSE << "reading cell " << (iCell++) << " "
              << MSG::hex << id << MSG::dec << " " << ene << " " 
              << MSG::hex << tqg << MSG::dec << endmsg;
          log << MSG::VERBOSE << "ene=" << ener << " time=" << time 
              << " qual=" << qual << " gain=" << gain << endmsg;
        }
        
        TileCell * cell = new TileCell(NULL,id,ener,time,qual,qbit,(CaloGain::CaloGain)gain);
        cont->push_back(cell);
      }
      break;
 
    case 2:
      
      for (; it != last; ) {
 
        Identifier id;
        CaloDetDescrElement * dde = NULL;
        float ener = 0.0;
        float time = 0.0;
        uint16_t qual = 0;
        uint16_t qbit = TileCell::MASK_CMPC | TileCell::MASK_TIME;
        int gain = m_gain[0]; // non-existing gain in CaloGain - to mark non-existing cells
        
        if (lVerbose)
          log << MSG::VERBOSE << "reading cell " << iCell << " ";
 
        if (iCell < NCELLMBTS) { // first 32 cells are MBTS cells without identifier
 
          id = m_id[iCell]; // identifier is taken from array
          dde = m_dde[iCell]; // mbtsDDE is taken from array
 
          int ene = (int)(*it++); // first word is energy
 
          if (lVerbose)
            log << MSG::hex << id << " " << ene << " " << MSG::dec;
 
          if (ene < 0 ) { // upper most bit is set, it means that everything is packed in one word
 
            if (lVerbose)
              log << endmsg;
 
            time = 0.0;   // time was zero and it was not saved
            ener = ((ene & 0xFFFFF) - 0x10000) * 1e-3;
            qual = ((ene>>20) & 0xFF);
            gain = m_gain[((ene>>28) & 0x7)]; // gain is taken from array
            
          } else { // two words packing
          
            unsigned int tqg = *it++;
            if (lVerbose)
              log << MSG::hex << tqg << MSG::dec << endmsg;
 
            ener = (ene - 0x40000000) * 1e-3;
            time = ((int)(tqg>>16) - 0x8000 ) * 0.01;
            qual = ((tqg>>8) & 0xFF);
            gain = (int)(tqg & 0xFF) - 0x80;
          }
 
        } else { // three words packing for remaining cells
 
          id = Identifier(Identifier32(*it++));
          int ene = (int)(*it++);
          unsigned int tqg = *it++;
      
          if (lVerbose)
            log << MSG::hex << id << MSG::dec << " " << ene << " " 
                << MSG::hex << tqg << MSG::dec << endmsg;
 
          ener = ene*1e-3;
          time = ((int)(tqg>>16) - 0x8000 ) * 0.01;
          qual = ((tqg>>8) & 0xFF);
          gain = (int)(tqg & 0xFF) - 0x80;
        }
 
        if (lVerbose)
          log << MSG::VERBOSE << "ene=" << ener << " time=" << time 
              << " qual=" << qual << " gain=" << gain << endmsg;
 
        if (gain != m_gain[0]) { // don't create cells with non-existing gain
          TileCell * cell = new TileCell(dde,id,ener,time,qual,qbit,(CaloGain::CaloGain)gain);
          cont->push_back(cell);
        }
        else {
          if (lVerbose)
            log << MSG::VERBOSE << "Don't create MBTS cell with invalid gain" << endmsg;
        }
        ++iCell;
      }
      break;
      
    default:
 
      log << MSG::ERROR << "Unknown version of TileCellVec, ver="<<version << endmsg;
    }
    
    return cont.release();
}

createTransientWithKey()

template<class TRANS , class PERS >

virtual TRANS* T_AthenaPoolCustomCnv< TRANS, PERS >::createTransientWithKey ( const std::string &  key )

overrideprotectedvirtualinherited

method to be implemented by the developer.

It has to find out the type of the persistent object to be read (by comparing GUIDs), read it using poolReadObject<P>(), call TLP converter to create a transient representation and return it. if the version 1 of poolReadObject is used, the persistent object HAS TO BE DELETED manually.

Implements T_AthenaPoolCustomCnvWithKey< TRANS, PERS >.

DataObjectToPers()

template<class TRANS , class PERS >

virtual StatusCode T_AthenaPoolCustomCnvWithKey< TRANS, PERS >::DataObjectToPers ( DataObject *  pObj, IOpaqueAddress *&  pAddr  )

overrideprotectedvirtualinherited

Convert an object into Persistent.

Parameters

pObj	[IN] pointer to the transient object.
key	[IN] StoreGate key (string) - placement hint to generate POOL container name

Reimplemented from T_AthenaPoolCustCnv< TRANS, PERS >.

DataObjectToPool()

template<class TRANS , class PERS >

virtual StatusCode T_AthenaPoolCustomCnvWithKey< TRANS, PERS >::DataObjectToPool ( IOpaqueAddress *  pAddr, DataObject *  pObj  )

overrideprotectedvirtualinherited

Write an object into POOL.

Parameters

pObj	[IN] pointer to the transient object.
key	[IN] StoreGate key (string) - placement hint to generate POOL container name

Reimplemented from T_AthenaPoolCustCnv< TRANS, PERS >.

Reimplemented in T_AthenaPoolLegacyCnv< TRANS, PERS, LEGACY >, and DataHeaderCnv.

fillRepRefs()

StatusCode AthenaPoolConverter::fillRepRefs ( IOpaqueAddress *  pAddr, DataObject *  pObj  )

overridevirtualinherited

Create a POOL persistent representation for a transient object.

Parameters

pAddr	[IN] IOpaqueAddress of POOL persistent representation.
pObj	[IN] pointer to the transient object.

Reimplemented in T_AthenaPoolCoolMultChanCnv< COLL_T, ELEM_T, ELEM_P >, T_AthenaPoolCoolMultChanCnv< LArHVScaleCorrComplete, LArConditionsSubset< LArSingleFloatP >, LArSingleFloatSubset_p1 >, and T_AthenaPoolCoolMultChanCnv< TransientClass, LArConditionsSubset< LArSingleFloatP >, LArSingleFloatSubset_p1 >.

Definition at line 142 of file AthenaPoolConverter.cxx.

                                                                                   {
   std::lock_guard<CallMutex> lock(m_conv_mut);
   try {
      if (!DataObjectToPool(pAddr, pObj).isSuccess()) {
         ATH_MSG_ERROR("FillRepRefs failed, key = " << pObj->name());
         return(StatusCode::FAILURE);
      }
   } catch (std::exception& e) {
      ATH_MSG_ERROR("fillRepRefs - caught exception: " << e.what());
      return(StatusCode::FAILURE);
   }
   return(StatusCode::SUCCESS);
}

finalize()

StatusCode AthenaPoolConverter::finalize ( )

overridevirtualinherited

Definition at line 52 of file AthenaPoolConverter.cxx.

                                         {
   // Release AthenaPoolCnvSvc
   if (!m_athenaPoolCnvSvc.release().isSuccess()) {
      ATH_MSG_WARNING("Cannot release AthenaPoolCnvSvc.");
   }
   return(::Converter::finalize());
}

getDataObject()

const DataObject * AthenaPoolConverter::getDataObject ( ) const

protectedvirtualinherited

Returns

data object from the converter.

Definition at line 269 of file AthenaPoolConverter.cxx.

                                                           {
   return(m_dataObject);
}

initialize()

StatusCode TileCellContainerCnv::initialize ( )

overridevirtual

initialization

Reimplemented from T_AthenaPoolCustCnv< TRANS, PERS >.

Definition at line 37 of file TileCellContainerCnv.cxx.

{
    // Call base clase initialize
    if (!AthenaPoolConverter::initialize().isSuccess()) {
        return StatusCode::FAILURE;
    }
 
    // version 2 by default
    m_version = 2;
    
    // Get the messaging service, print where you are
    MsgStream log(msgSvc(), "TileCellContainerCnv");
    log << MSG::INFO << "TileCellContainerCnv::initialize(), packing format version " << m_version << endmsg;
 
    // get StoreGate service
    StatusCode sc=service("StoreGateSvc",m_storeGate);
    if (sc.isFailure()) {
      this->initIdToIndex();
      log << MSG::ERROR << "StoreGate service not found !" << endmsg;
      // log << MSG::FATAL << "StoreGate service not found !" << endmsg;
      // return StatusCode::FAILURE;
    }
 
    StoreGateSvc* detStore;
    sc=service("DetectorStore",detStore);
    if (sc.isFailure()) {
      this->initIdToIndex();
      log << MSG::ERROR << "DetectorStore service not found!" << endmsg;
      // log << MSG::FATAL << "DetectorStore service not found!" << endmsg;
      // return StatusCode::FAILURE;
    }
 
    sc = detStore->retrieve(m_tileTBID);
    if (sc.isFailure()) {
      this->initIdToIndex();
      log << MSG::ERROR << "No TileTBID helper" << endmsg;
      // log << MSG::FATAL << "No TileTBID helper" << endmsg;
      // return StatusCode::FAILURE;
    } else {
      for (int side=0; side<NSIDE; ++side) {
        for (int phi=0; phi<NPHI; ++phi) {
          for (int eta=0; eta<NETA; ++eta) {
            m_id[cell_index(side,phi,eta)] = m_tileTBID->channel_id((side>0)?1:-1,phi,eta);
          }
        }
      }
    }
    
    sc = detStore->retrieve(m_mbtsMgr);
    if (sc.isFailure()) {
      log << MSG::WARNING << "Unable to retrieve MbtsDetDescrManager from DetectorStore" << endmsg;
      memset(m_dde,0,sizeof(m_dde));
    } else {
      for (int side=0; side<NSIDE; ++side) {
        for (int phi=0; phi<NPHI; ++phi) {
          for (int eta=0; eta<NETA; ++eta) {
            m_dde[cell_index(side,phi,eta)] = m_mbtsMgr->get_element(m_id[cell_index(side,phi,eta)]);
          }
        }
      }
    }
    
 
    // set CaloGain <-> gain index mapping for all possible TileCal gains
    for (int i=0; i<17; ++i)  m_gainIndex[i] = 8;
    m_gain[0] = -2; // put non-existing gain here
    m_gainIndex[-CaloGain::TILELOWLOW]   = 8 + 1;  m_gain[1] = (int)CaloGain::TILELOWLOW;
    m_gainIndex[-CaloGain::TILELOWHIGH]  = 8 + 2;  m_gain[2] = (int)CaloGain::TILELOWHIGH;
    m_gainIndex[-CaloGain::TILEHIGHLOW]  = 8 + 3;  m_gain[3] = (int)CaloGain::TILEHIGHLOW;
    m_gainIndex[-CaloGain::TILEHIGHHIGH] = 8 + 4;  m_gain[4] = (int)CaloGain::TILEHIGHHIGH;
    m_gainIndex[-CaloGain::TILEONELOW]   = 8 + 5;  m_gain[5] = (int)CaloGain::TILEONELOW;
    m_gainIndex[-CaloGain::TILEONEHIGH]  = 8 + 6;  m_gain[6] = (int)CaloGain::TILEONEHIGH;
    m_gainIndex[-CaloGain::INVALIDGAIN]  = 8 + 7;  m_gain[7] = (int)CaloGain::INVALIDGAIN;
    
    return StatusCode::SUCCESS;
}

initIdToIndex()

void TileCellContainerCnv::initIdToIndex ( )

private

Definition at line 114 of file TileCellContainerCnv.cxx.

{
  for (int side=0; side<NSIDE; ++side) {
    for (int phi=0; phi<NPHI; ++phi) {
      for (int eta=0; eta<NETA; ++eta) {
        m_id[cell_index(side,phi,eta)] = (Identifier)cell_index(side,phi,eta);
      }
    }
  }
}

initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

{
  m_imsg = Athena::getMessageSvc();
  m_lvl = m_imsg ?
    static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
    MSG::INFO;
}

ipCorr()

void Converter::ipCorr ( double  d0, double  z0, double &  d0c, double &  z0c, double  phi0, double  eta, double  pt  )

inlineprotectedinherited

Definition at line 257 of file Converter.h.

                                                                                                  {
    
    double sn  = sin(phi0);
    double cs  = cos(phi0);
    double sd0 = (d0 != 0 ? d0/fabs(d0) : 1);
    double spt = (pt != 0 ? pt/fabs(pt) : 1);
    
    if (fabs(pt) >= 1*CLHEP::TeV) {
 
      d0c = d0 + m_beamX*sn - m_beamY*cs;
 
    } else {
      
      double rc = fabs(pt)*15.0/(9.0*1.042);
      
      double xc = (fabs(d0)-spt*sd0*rc)*cos(phi0+M_PI/2*sd0) - m_beamX;
      double yc = (fabs(d0)-spt*sd0*rc)*sin(phi0+M_PI/2*sd0) - m_beamY;
      
      double newphi;
      double xd01,yd01,xd02,yd02;
      
      if (xc == 0) {
    xd01 = 0; yd01 = rc+yc;
    xd02 = 0; yd02 = -rc+yc; 
      } else {
    xd01 = xc+yc/xc*yc+sqrt(pow((xc+yc/xc*yc),2)-xc*xc-yc*yc+rc*rc); yd01 = yc/xc*xd01;
    xd02 = xc+yc/xc*yc-sqrt(pow((xc+yc/xc*yc),2)-xc*xc-yc*yc+rc*rc); yd02 = yc/xc*xd02;
      }
      
      double r1 = sqrt(xd01*xd01+yd01*yd01);
      double r2 = sqrt(xd02*xd02+yd02*yd02);
      
      double phiV;
      
      if (r1 < r2)
    phiV = atan2(yd01,xd01); 
      else
    phiV = atan2(yd02,xd02); 
      
      double phi1 = phiCorr(phiV+M_PI/2);
      double phi2 = phiCorr(phiV-M_PI/2);
      
      if (fabs(phiCorr(phi1-phi0))<=fabs(phiCorr(phi2-phi0)))
    newphi = phi1;
      else
    newphi = phi2;
      
      d0c = fabs(sqrt(xc*xc+yc*yc)-rc)*sin(phiV-newphi);
      
      double theta=2*atan2(exp(-eta),1);
      double theta_save=theta;
      theta = theta - (1+spt)/2*M_PI;
      if (theta>0) theta = theta_save;
      
      double deltaz0= -spt*rc/tan(theta)*phiCorr(phiCorr(newphi)-phiCorr(phi0));
      z0c = z0 + deltaz0;
      
    }
  }

keepPoolObj()

template<class TRANS , class PERS >

void T_AthenaPoolCustomCnvWithKey< TRANS, PERS >::keepPoolObj ( PERS *  obj, const std::string &  output  )

protectedinherited

Remember the POOL object to be written out (will be deleted after commit)

Parameters

obj	[IN] persistent object

msg() [1/2]

MsgStream & AthMessaging::msg ( ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 164 of file AthMessaging.h.

{
  MsgStream* ms = m_msg_tls.get();
  if (!ms) {
    if (!m_initialized.test_and_set()) initMessaging();
    ms = new MsgStream(m_imsg,m_nm);
    m_msg_tls.reset( ms );
  }
 
  ms->setLevel (m_lvl);
  return *ms;
}

msg() [2/2]

MsgStream & AthMessaging::msg ( const MSG::Level  lvl ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 179 of file AthMessaging.h.

{ return msg() << lvl; }

msgLvl()

bool AthMessaging::msgLvl ( const MSG::Level  lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns

boolean Indicating if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

{
  if (!m_initialized.test_and_set()) initMessaging();
  if (m_lvl <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

objectToPool()

template<class TRANS , class PERS >

template<class P >

StatusCode T_AthenaPoolCustCnv< TRANS, PERS >::objectToPool ( P *  pObj, Token *&  token, const std::string &  key, const std::string &  output  )

protectedinherited

Write an object into POOL returning its token.

Parameters

pObj	[IN] pointer to the object to be written.
token	[OUT] POOL token of the persistent representation.
key	[IN] StoreGate key (string) - placement hint to generate POOL container name

persToTrans()

template<class TRANS , class PERS >

virtual StatusCode T_AthenaPoolCustomCnvWithKey< TRANS, PERS >::persToTrans ( TRANS *&  , PERS *    )

inlineoverrideprotectedvirtualinherited

obsolete

Implements T_AthenaPoolCustCnv< TRANS, PERS >.

Definition at line 96 of file T_AthenaPoolCustomCnv.h.

{ return(StatusCode::FAILURE); }

phiCorr()

double Converter::phiCorr ( double  phi )

inlineprotectedinherited

Definition at line 248 of file Converter.h.

                            {
    
    if (phi < -TMath::Pi()) phi += 2*TMath::Pi();
    if (phi >  TMath::Pi()) phi -= 2*TMath::Pi();
    
    return phi;
  }

poolReadObject() [1/2]

template<class TRANS , class PERS >

template<class P >

P* T_AthenaPoolCustomCnvWithKey< TRANS, PERS >::poolReadObject ( )

protectedinherited

Read object of type P.

This is an exception-throwing version of poolToObject() plus reading of all extending objects. Version 1 - (see createTransient() above)

Returns

object read from POOL (by pointer)

poolReadObject() [2/2]

template<class TRANS , class PERS >

template<class P >

void T_AthenaPoolCustomCnvWithKey< TRANS, PERS >::poolReadObject ( TopLevelTPCnvBase &  tlp_converter )

protectedinherited

Read object of type P (plus all extending objects) using the indicated top-level TP converter.

Version 2 - (see createTransient() above) NOTE: the TLP converter will delete the persistent object after createTransient()

Parameters

tlp_converter

[IN] top-level TP converter to be used when reading

PoolToDataObject()

template<class TRANS , class PERS >

virtual StatusCode T_AthenaPoolCustomCnvWithKey< TRANS, PERS >::PoolToDataObject ( DataObject *&  pObj, const Token *  token, const std::string &  key  )

overrideprotectedvirtualinherited

Read an object from POOL.

Parameters

pObj	[OUT] pointer to the transient object.
token	[IN] POOL token of the persistent representation.
key	[IN] SG key of the object being read.

Reimplemented from T_AthenaPoolCustCnv< TRANS, PERS >.

poolToObject()

template<class TRANS , class PERS >

template<class P >

StatusCode T_AthenaPoolCustCnv< TRANS, PERS >::poolToObject ( const Token *&  token, P *&  pObj  )

protectedinherited

Read an object from POOL.

Parameters

token	[IN] POOL token of the persistent representation.
pObj	[OUT] pointer to the object read.

repSvcType()

long AthenaPoolConverter::repSvcType ( ) const

overridevirtualinherited

Definition at line 60 of file AthenaPoolConverter.cxx.

                                           {
   return(POOL_StorageType);
}

round16()

int TileCellContainerCnv::round16 ( double  x )

inlineprivate

Definition at line 59 of file TileCellContainerCnv.h.

                                 { 
      if (x<-32767.) return -0x7FFF;
      else if (x>32767.) return 0x7FFF;
      else return (int)lround(x);
    }

round32()

int TileCellContainerCnv::round32 ( double  x )

inlineprivate

Definition at line 53 of file TileCellContainerCnv.h.

                                 { 
      if (x<-2147483647.) return -0x7FFFFFFF;
      else if (x>2147483647.) return 0x7FFFFFFF;
      else return (int)lround(x);
    }

selectTracks() [1/3]

void Converter::selectTracks ( const Analysis::MuonContainer *  muontracks )

inlineinherited

Definition at line 177 of file Converter.h.

                                                               {
    
    static int hpmap[20] = { 0, 1, 2,  7, 8, 9,  3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
    
    Analysis::MuonContainer::const_iterator muonitr = muontracks->begin();
    Analysis::MuonContainer::const_iterator muonend = muontracks->end();
    
    while(muonitr!=muonend) { 
      
      bool hasInDet = (*muonitr)->hasInDetTrackParticle();
      const Rec::TrackParticle *trackitr=NULL; 
      if(hasInDet) trackitr=(*muonitr)->inDetTrackParticle();
      muonitr++;
      if(!hasInDet) continue;
      
      const Trk::MeasuredPerigee* measPer = trackitr->measuredPerigee();
      
      CLHEP::HepVector perigeeParams = measPer->parameters();    
      double pT    = measPer->pT(); 
      double eta   = measPer->eta();
      double phi   = perigeeParams[Trk::phi0];
      double z0    = perigeeParams[Trk::z0];
      double d0    = perigeeParams[Trk::d0];
      //double theta = perigeeParams[Trk::theta];
      
      double deta = 1;
      double dphi = 1;
      double dz0  = 1;
      double dd0  = 1;
      double dpT  = 1;    
      
      // Check number of hits
      // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really 
      // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
      // stereo clusters making a spacepoint are two "hits"
      const Trk::TrackSummary *summary = trackitr->trackSummary();
      int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits); 
      int nPixelHits  = 2*summary->get(Trk::numberOfPixelHits);  
      int nSctHits    = summary->get(Trk::numberOfSCTHits); 
      int nStrawHits  = summary->get(Trk::numberOfTRTHits);
      int nTrHits     = summary->get(Trk::numberOfTRTHighThresholdHits);
      
      int nSiHits     = nPixelHits + nSctHits;
      
      const Trk::FitQuality *quality   = trackitr->fitQuality();
      double chi2 = quality->chiSquared();
      double dof  = quality->quality->numberDoF();
 
      unsigned bitmap = 0;
      
      for ( int ih=0 ; ih<20 ; ih++ ) {
    if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );    
      }
      
      // Shift d0 and z0 according to beam position
      ipCorr(d0, z0, d0, z0, phi, eta, pT);
      
      // Create and save Track  
      TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
                                 deta, dphi, dz0, dd0, dpT,
                                 nBlayerHits, nPixelHits, nSctHits, nSiHits,
                                 nStrawHits, nTrHits, bitmap, 0,
                                 -1) ;
      
      addTrack(t);
    }
    
  }

selectTracks() [2/3]

void Converter::selectTracks ( const Rec::TrackParticleContainer *  trigtracks )

inlineinherited

Definition at line 111 of file Converter.h.

                                                                   { 
    
    Rec::TrackParticleContainer::const_iterator  trackitr = trigtracks->begin();
    Rec::TrackParticleContainer::const_iterator  trackend = trigtracks->end();
    
    static int hpmap[20] = { 0, 1, 2,  7, 8, 9,  3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
 
    while ( trackitr!=trackend ) { 
        
      const Trk::MeasuredPerigee* measPer = (*trackitr)->measuredPerigee();
 
      CLHEP::HepVector perigeeParams = measPer->parameters();    
      double pT    = measPer->pT(); 
      double eta   = measPer->eta();
      double phi   = perigeeParams[Trk::phi0];
      double z0    = perigeeParams[Trk::z0];
      double d0    = perigeeParams[Trk::d0];
      //double theta = perigeeParams[Trk::theta];
    
      double deta = 1;
      double dphi = 1;
      double dz0  = 1;
      double dd0  = 1;
      double dpT  = 1;    
 
      // Check number of hits
      // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really 
      // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
      // stereo clusters making a spacepoint are two "hits"
      const Trk::TrackSummary *summary = (*trackitr)->trackSummary();
      int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits); 
      int nPixelHits  = 2*summary->get(Trk::numberOfPixelHits);  
      int nSctHits    = summary->get(Trk::numberOfSCTHits); 
      int nStrawHits  = summary->get(Trk::numberOfTRTHits);
      int nTrHits     = summary->get(Trk::numberOfTRTHighThresholdHits);
      int nSiHits     = nPixelHits + nSctHits;
      bool expectBL   = false; // Not stored for Rec::TrackParticle
 
      const Trk::FitQuality *quality   = (*trackitr)->fitQuality();
      double chi2 = quality->chiSquared();
      double dof  = quality->numberDoF();
 
      unsigned bitmap = 0;
      
      for ( int ih=0 ; ih<20 ; ih++ ) {
    if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );    
      }
 
      // Shift d0 and z0 according to beam position
      ipCorr(d0, z0, d0, z0, phi, eta, pT);
 
      // Create and save Track  
      TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof, 
                                 deta, dphi, dz0, dd0, dpT,
                                 nBlayerHits, nPixelHits, nSctHits, nSiHits,
                                 nStrawHits, nTrHits, bitmap, 0,
                                 -1,
                                 expectBL) ;
  
      addTrack(t);
      trackitr++;
    }
    
  }

selectTracks() [3/3]

void Converter::selectTracks ( const TrigInDetTrackCollection *  trigtracks )

inlineinherited

Definition at line 56 of file Converter.h.

                                                                {
    
    TrigInDetTrackCollection::const_iterator  trackitr = trigtracks->begin();
    TrigInDetTrackCollection::const_iterator  trackend = trigtracks->end();
 
    while ( trackitr!=trackend ) { 
 
    double eta    = (*trackitr)->param()->eta();
    double phi    = (*trackitr)->param()->phi0();
    double z0     = (*trackitr)->param()->z0(); 
    double pT     = (*trackitr)->param()->pT(); 
    double d0     = (*trackitr)->param()->a0();
    //double theta  = 2*atan2(exp(-(*trackitr)->param()->eta()),1);
 
    double deta    = (*trackitr)->param()->eeta();
    double dphi    = (*trackitr)->param()->ephi0();
    double dz0     = (*trackitr)->param()->ez0(); 
    double dpT     = (*trackitr)->param()->epT(); 
    double dd0     = (*trackitr)->param()->ea0();
 
    int   algoid  = (*trackitr)->algorithmId();           
 
    int nBlayerHits = ((*trackitr)->HitPattern() & 0x1);
    int nPixelHits  = 2*(*trackitr)->NPixelSpacePoints();  // NB: for comparison with offline 
    int nSctHits    = 2*(*trackitr)->NSCT_SpacePoints();   //     a spacepoint is 2 "hits"
    int nStrawHits  = (*trackitr)->NStrawHits();
    int nTrHits     = (*trackitr)->NTRHits();
               
    int nSiHits     = nPixelHits + nSctHits;
    bool expectBL   = false;                               //not available with TrigInDetTrack
 
    unsigned hitPattern = (*trackitr)->HitPattern();
    unsigned multiPattern = 0;
 
    double chi2    = (*trackitr)->chi2();
    double dof     = 0;
 
    // Shift d0 and z0 according to beam position
    ipCorr(d0, z0, d0, z0, phi, eta, pT);
    
    // Create and save Track  
    TIDA::Track* t = new TIDA::Track(eta,  phi,  z0,  d0,  pT, chi2, dof, 
                                   deta, dphi, dz0, dd0, dpT, 
                                   nBlayerHits, nPixelHits, nSctHits, nSiHits, 
                                   nStrawHits, nTrHits, hitPattern, multiPattern, 
                                   algoid,
                                   expectBL ) ; 
 
    addTrack(t);
    trackitr++; 
    }
  }

setLevel()

void AthMessaging::setLevel ( MSG::Level  lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

{
  m_lvl = lvl;
}

setPlacement()

virtual Placement T_AthenaPoolCnvBase< TRANS >::setPlacement ( const std::string &  key, const std::string &  output  )

protectedvirtualinherited

Set POOL placement.

Implements AthenaPoolConverter.

setPlacementForP()

template<class TRANS , class PERS >

template<class P >

Placement T_AthenaPoolCustCnv< TRANS, PERS >::setPlacementForP ( P &  p, const std::string &  key, const std::string &  output  )

protectedinherited

Set POOL placement for type P.

setPlacementWithType()

Placement AthenaPoolConverter::setPlacementWithType ( const std::string &  tname, const std::string &  key, const std::string &  output  )

protectedvirtualinherited

Set POOL placement hint for a given type.

Parameters

tname	[IN] type name.
key	[IN] SG key.

Definition at line 176 of file AthenaPoolConverter.cxx.

                                                                                                                         {
   Placement placement;
   // Override streaming parameters from StreamTool if requested.
   std::string::size_type pos1 = output.find('[');
   std::string outputConnectionSpec = output.substr(0, pos1);
   int tech = 0;
   m_athenaPoolCnvSvc->decodeOutputSpec(outputConnectionSpec, tech).ignore();
   // Set DB and Container names
   placement.setFileName(outputConnectionSpec);
 
   std::string containerPrefix = m_containerPrefix;
   std::string dhContainerPrefix = pool::ROOTRNTUPLE_StorageType.exactMatch(tech) ? APRDefaults::RNTupleNames::DataHeader : APRDefaults::TTreeNames::DataHeader;
   std::string containerName;
 
   // Get Technology from containerPrefix
   std::size_t colonPos = containerPrefix.find(':');
   if (colonPos != std::string::npos) {
      dhContainerPrefix = containerPrefix.substr(0, colonPos + 1) + dhContainerPrefix;
   }
 
   // Override streaming parameters from StreamTool if requested.
   std::string containerNameHint = m_containerNameHint;
   std::string branchNameHint = m_branchNameHint;
   std::string containerFriendPostfix;
   while (pos1 != std::string::npos) {
      const std::string::size_type pos2 = output.find('=', pos1);
      const std::string thisKey = output.substr(pos1 + 1, pos2 - pos1 - 1);
      const std::string::size_type pos3 = output.find(']', pos2);
      const std::string value = output.substr(pos2 + 1, pos3 - pos2 - 1);
      if (thisKey == "OutputCollection") {
         dhContainerPrefix = value;
      } else if (thisKey == "PoolContainerPrefix") {
         containerPrefix = value;
      } else if (thisKey == "TopLevelContainerName") {
         containerNameHint = value;
      } else if (thisKey == "SubLevelBranchName") {
         branchNameHint = value;
      } else if (thisKey == "PoolContainerFriendPostfix") {
         containerFriendPostfix = value;
      }
      pos1 = output.find('[', pos3);
   }
 
   // ---  Special types:   DataHeader & Form
   if( tname.compare(0, 10, "DataHeader") == 0 ) {
      if( tname.compare(10, 4, "Form") == 0 ) {
         containerName = dhContainerPrefix + "Form" + "(" + tname + ")";
      } else {
         if (key[key.size() - 1] == '/') {
            containerName = dhContainerPrefix + "(" + key + tname + ")";
         } else {
            containerName = dhContainerPrefix + "(" + tname + ")";
         }
      }
   }
   // AttributeList - writing attributes separately to EventTag container group
   else if (tname.compare(0, 13, "AttributeList") == 0) {
      // Find the right storage type and name for EventTag values
      if( pool::ROOTRNTUPLE_StorageType.exactMatch(tech) ) {
         containerName = std::string(APRDefaults::RNTupleNames::EventTag) + "(" + key + ")";
      } else {
         // no indexing needed (nothing points to Tags)
         // safe to set tech here - it will not be overwritten by decodeOutput
         tech = pool::ROOTTREE_StorageType.type();
         containerName = std::string(APRDefaults::TTreeNames::EventTag) + "(" + key + ")";
      }
   }
   // all other object types
   else {
      const std::string typeTok = "<type>", keyTok = "<key>";
      containerName = containerPrefix + containerFriendPostfix + containerNameHint;
      if (!branchNameHint.empty()) {
         containerName += "(" + branchNameHint + ")";
      }
      const std::size_t pos1 = containerName.find(typeTok);
      if (pos1 != std::string::npos) {
         containerName.replace(pos1, typeTok.size(), tname);
      }
      const std::size_t pos2 = containerName.find(keyTok);
      if (pos2 != std::string::npos) {
         if (key.empty()) {
            containerName.replace(pos2, keyTok.size(), tname);
         } else {
            containerName.replace(pos2, keyTok.size(), key);
         }
      }
   }
   m_athenaPoolCnvSvc->decodeOutputSpec(containerName, tech).ignore();
   placement.setContainerName(containerName);
   placement.setTechnology(tech);
   return(placement);
}

setToken()

template<class TRANS , class PERS >

virtual void T_AthenaPoolCustCnv< TRANS, PERS >::setToken ( const std::string &  token )

protectedvirtualinherited

Set the token (in std::string representation) and classID for the object that will be read next.

Required by compareClassGuid() and poolReadObject(). Only extending converter needs to call this method explicitely

Parameters

token

[IN] token for the object to be read

Reimplemented in T_AthenaPoolCoolMultChanCnv< COLL_T, ELEM_T, ELEM_P >, T_AthenaPoolCoolMultChanCnv< LArHVScaleCorrComplete, LArConditionsSubset< LArSingleFloatP >, LArSingleFloatSubset_p1 >, and T_AthenaPoolCoolMultChanCnv< TransientClass, LArConditionsSubset< LArSingleFloatP >, LArSingleFloatSubset_p1 >.

storageType()

long AthenaPoolConverter::storageType ( )

staticinherited

Returns

storage type.

Definition at line 156 of file AthenaPoolConverter.cxx.

                                      {
   return(POOL_StorageType);
}

tracks()

const std::vector<TIDA::Track*>& Converter::tracks ( ) const

inlineinherited

Definition at line 50 of file Converter.h.

{return m_tracks;}

transToPers()

template<class TRANS , class PERS >

virtual StatusCode T_AthenaPoolCustomCnvWithKey< TRANS, PERS >::transToPers ( TRANS *  , PERS *&    )

inlineoverrideprotectedvirtualinherited

Obsolete methods replaced by createPersistent() and createTransient() obsolete.

Implements T_AthenaPoolCustCnv< TRANS, PERS >.

Definition at line 94 of file T_AthenaPoolCustomCnv.h.

{ return(StatusCode::FAILURE); }

Friends And Related Function Documentation

CnvFactory< TileCellContainerCnv >

friend class CnvFactory< TileCellContainerCnv >

friend

Definition at line 16 of file TileCellContainerCnv.h.

Member Data Documentation

ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging)

Definition at line 141 of file AthMessaging.h.

m_athenaPoolCnvSvc

ServiceHandle<IAthenaPoolCnvSvc> AthenaPoolConverter::m_athenaPoolCnvSvc

protectedinherited

Definition at line 111 of file AthenaPoolConverter.h.

m_beamX

double Converter::m_beamX

protectedinherited

Definition at line 320 of file Converter.h.

m_beamY

double Converter::m_beamY

protectedinherited

Definition at line 320 of file Converter.h.

m_beamZ

double Converter::m_beamZ

protectedinherited

Definition at line 320 of file Converter.h.

m_branchNameHint

std::string AthenaPoolConverter::m_branchNameHint

protectedinherited

Definition at line 121 of file AthenaPoolConverter.h.

m_classDesc

RootType AthenaPoolConverter::m_classDesc

protectedinherited

Definition at line 112 of file AthenaPoolConverter.h.

m_classDescs

ClassMap AthenaPoolConverter::m_classDescs

protectedinherited

Definition at line 117 of file AthenaPoolConverter.h.

m_classID

template<class TRANS , class PERS >

Guid T_AthenaPoolCustCnv< TRANS, PERS >::m_classID

protectedinherited

Definition at line 96 of file T_AthenaPoolCustCnv.h.

m_className

std::string AthenaPoolConverter::m_className

protectedinherited

Definition at line 116 of file AthenaPoolConverter.h.

m_containerNameHint

std::string AthenaPoolConverter::m_containerNameHint

protectedinherited

Definition at line 120 of file AthenaPoolConverter.h.

m_containerPrefix

std::string AthenaPoolConverter::m_containerPrefix

protectedinherited

Definition at line 119 of file AthenaPoolConverter.h.

m_conv_mut

CallMutex AthenaPoolConverter::m_conv_mut

protectedinherited

Definition at line 127 of file AthenaPoolConverter.h.

m_dataObject

const DataObject* AthenaPoolConverter::m_dataObject

protectedinherited

Definition at line 123 of file AthenaPoolConverter.h.

m_dde

CaloDetDescrElement* TileCellContainerCnv::m_dde[NCELLMBTS]

private

Definition at line 49 of file TileCellContainerCnv.h.

m_gain

int TileCellContainerCnv::m_gain[8]

private

Definition at line 51 of file TileCellContainerCnv.h.

m_gainIndex

int TileCellContainerCnv::m_gainIndex[17]

private

Definition at line 50 of file TileCellContainerCnv.h.

m_i_poolToken

const Token* AthenaPoolConverter::m_i_poolToken

protectedinherited

Definition at line 124 of file AthenaPoolConverter.h.

m_id

Identifier TileCellContainerCnv::m_id[NCELLMBTS]

private

Definition at line 48 of file TileCellContainerCnv.h.

m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

m_mbtsMgr

const MbtsDetDescrManager* TileCellContainerCnv::m_mbtsMgr

private

Definition at line 36 of file TileCellContainerCnv.h.

m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

m_persMap

template<class TRANS , class PERS >

std::map<DataObject*, PERS*> T_AthenaPoolCustCnv< TRANS, PERS >::m_persMap

privateinherited

Definition at line 103 of file T_AthenaPoolCustCnv.h.

m_persObjLists

template<class TRANS , class PERS >

std::map<std::string, std::vector< std::unique_ptr<PERS> > > T_AthenaPoolCustomCnvWithKey< TRANS, PERS >::m_persObjLists

protectedinherited

Local cache for persistent objects created by this converter, grouped by processing slot These objects are deleted after a commit.

Definition at line 120 of file T_AthenaPoolCustomCnv.h.

m_pListMutex

template<class TRANS , class PERS >

std::mutex T_AthenaPoolCustomCnvWithKey< TRANS, PERS >::m_pListMutex

protectedinherited

protection mutex for m_persObjLists

Definition at line 123 of file T_AthenaPoolCustomCnv.h.

m_storeGate

StoreGateSvc* TileCellContainerCnv::m_storeGate

private

Definition at line 34 of file TileCellContainerCnv.h.

m_tileTBID

const TileTBID* TileCellContainerCnv::m_tileTBID

private

Definition at line 35 of file TileCellContainerCnv.h.

m_tracks

std::vector<TIDA::Track*> Converter::m_tracks

protectedinherited

Definition at line 321 of file Converter.h.

m_version

int TileCellContainerCnv::m_version

private

Definition at line 38 of file TileCellContainerCnv.h.

NCELLMBTS

const int TileCellContainerCnv::NCELLMBTS = NSIDE*NPHI*NETA

staticprivate

Definition at line 43 of file TileCellContainerCnv.h.

NETA

const int TileCellContainerCnv::NETA = 2

staticprivate

Definition at line 42 of file TileCellContainerCnv.h.

NPHI

const int TileCellContainerCnv::NPHI = 8

staticprivate

Definition at line 41 of file TileCellContainerCnv.h.

NSIDE

const int TileCellContainerCnv::NSIDE = 2

staticprivate

Definition at line 40 of file TileCellContainerCnv.h.

---

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

• TileCellContainerCnv.h
• TileCellContainerCnv.cxx