#include <TRT_HitCollectionCnv_p3.h>

Inheritance diagram for TRT_HitCollectionCnv_p3:

Collaboration diagram for TRT_HitCollectionCnv_p3:

Public Types
typedef TRANS	Trans_t
typedef PERS	Pers_t
typedef PERS	PersBase_t
typedef TRANS	TransBase_t
typedef ITPConverterFor< TRANS >	PolyCnvBase_t
typedef Gaudi::PluginService::Factory< ITPCnvBase *()>	Factory

Public Member Functions
	TRT_HitCollectionCnv_p3 ()
virtual TRTUncompressedHitCollection *	createTransient (const TRT_HitCollection_p3 *persObj, MsgStream &log)
virtual void	persToTrans (const TRT_HitCollection_p3 persCont, TRTUncompressedHitCollection transCont, MsgStream &log)
virtual void	transToPers (const TRTUncompressedHitCollection transCont, TRT_HitCollection_p3 persCont, MsgStream &log)
virtual TPObjRef	virt_toPersistent (const TRANS *trans, MsgStream &log)
	Internal interface method that is used to invoke the real conversion method (toPersistent_impl) in the derived converter.
virtual TPObjRef	virt_toPersistentWithKey (const TRANS *trans, const std::string &key, MsgStream &log)
	Internal interface method that is used to invoke the real conversion method (toPersistent_impl) in the derived converter.
virtual void	pstoreToTrans (unsigned index, TRANS *trans, MsgStream &log)
	Convert persistent representation stored in the storage vector of the top-level object to transient.
virtual TRANS *	createTransient (const PERS *persObj, MsgStream &log)
	Create transient representation of a persistent object.
virtual TRANS *	createTransientWithKey (const PERS *persObj, const std::string &key, MsgStream &log)
	Create transient representation of a persistent object, with SG key.
virtual TRANS *	virt_createTransFromPStore (unsigned index, MsgStream &log)
	Internal interface method that is used to invoke the real conversion method (createTransient)
virtual TRANS *	virt_createTransFromPStoreWithKey (unsigned index, const std::string &key, MsgStream &log)
	Internal interface method that is used to invoke the real conversion method (createTransient)
virtual void	persToTrans (const PERS persObj, TRANS transObj, MsgStream &log)=0
	Convert persistent representation to transient one.
virtual void	transToPers (const TRANS transObj, PERS persObj, MsgStream &log)=0
	Convert transient representation to persistent one.
virtual void	persToTransWithKey (const PERS persObj, TRANS transObj, const std::string &, MsgStream &log)
	Convert persistent representation to transient one.
virtual void	transToPersWithKey (const TRANS transObj, PERS persObj, const std::string &, MsgStream &log)
	Convert transient representation to persistent one.
virtual void	persToTransUntyped (const void pers, void trans, MsgStream &log)
	Convert persistent object representation to transient.
virtual void	transToPersUntyped (const void trans, void pers, MsgStream &log)
	Convert transient object representation to persistent.
virtual void	persToTransWithKeyUntyped (const void pers, void trans, const std::string &key, MsgStream &log)
	Convert persistent object representation to transient.
virtual void	transToPersWithKeyUntyped (const void trans, void pers, const std::string &key, MsgStream &log)
	Convert transient object representation to persistent.
virtual PERS *	createPersistent (const TRANS *transObj, MsgStream &log)
	Create persistent representation of a transient object.
virtual PERS *	createPersistentWithKey (const TRANS *transObj, const std::string &key, MsgStream &log)
	Create persistent representation of a transient object, with SG key.
TPObjRef	toPersistentWithKey_impl (const TRANS *trans, const std::string &key, MsgStream &log)
	Convert transient object to persistent representation.
virtual const std::type_info &	transientTInfo () const
	return C++ type id of the transient class this converter is for
virtual const std::type_info &	persistentTInfo () const
	return C++ type id of the persistent class this converter is for
void	setPStorage (std::vector< PERS > *storage)
	Tell this converter which storage vector it should use to store or retrieve persistent representations.
void	setRecursive (bool flag=true)
	Tell the converter if it should work in recursive mode slower but it can safely handle recursion.
void	ignoreRecursion (bool flag=false)
	Tell the converter to ignore recursion (do not throw errors) even when recurion is detected.
virtual void	reservePStorage (size_t size)
	Reserve 'size' elements for persistent storage.
template<class CNV>
CNV *	converterForType (CNV *cnv, const std::type_info &t_info, MsgStream &log) const
	Find converter for a given C++ type ID, that is or ihnerits from CNV type.
template<class CNV>
CNV *	converterForRef (CNV *cnv, const TPObjRef &ref, MsgStream &log) const
	Find converter for a TP type ID (passed in a TP Ref), that is or ihnerits from CNV type.
template<class CNV>
TPObjRef	baseToPersistent (CNV *cnv, const typename CNV::Trans_t transObj, MsgStream &log) const
	Persistify bass class of a given object and store the persistent represenation in the storage vector of the top-level persistent object.
template<class CNV>
TPObjRef	toPersistent (CNV *cnv, const typename CNV::TransBase_t transObj, MsgStream &log) const
	Persistify an object and store the persistent represenation in the storage vector of the top-level persistent object.
template<class CNV, class TRANS_T>
void	fillTransFromPStore (CNV *cnv, const TPObjRef &ref, TRANS_T trans, MsgStream &log) const
	Convert persistent object, stored in the the top-level persistent object and referenced by the TP Ref, to transient representation.
template<class CNV>
CNV::Trans_t *	createTransFromPStore (CNV **cnv, const TPObjRef &ref, MsgStream &log) const
	Create transient representation of a persistent object, stored in the the top-level persistent object and referenced by the TP Ref.
virtual void	initPrivateConverters (TopLevelTPCnvBase *)
virtual TopLevelTPCnvBase *	topConverter ()
	return the top-level converter for this elemental TP converter
virtual const TopLevelTPCnvBase *	topConverter () const
	return the top-level converter for this elemental TP converter
const std::type_info &	transBaseTInfo () const
	return C++ type id of the common base transient type for all converters for a group of polymorphic types
virtual const TPObjRef::typeID_t &	typeID () const
	Return TP typeID for persistent objects produced by this converter.
unsigned	typeIDvalue () const
	inlined non-virtual version to get the typeID value fast
virtual void	setRuntimeTopConverter (TopLevelTPCnvBase *topConverter)
	Set runtime top-level converter - usually it is the owning TL converter, but in case of extended objects it will be the TL converter of the extended object.
virtual void	setTopConverter (TopLevelTPCnvBase *topConverter, const TPObjRef::typeID_t &TPtypeID)
	Set which top-level converter owns this elemental converter, and what TPtypeID was assigned to the persistent objects it produces.
void	setReadingFlag ()
void	clearReadingFlag ()
bool	wasUsedForReading ()
virtual void	converterNotFound (const std::type_info &converterType, ITPConverter *c, const std::string &typeName, MsgStream &log) const
	method called when the right TP converter was not found during writing
virtual void	converterNotFound (unsigned typeID, ITPConverter *c, const std::string &typeName, MsgStream &log) const
	method called when the right TP converter was not found during reading

Protected Attributes
std::vector< PERS > *	m_pStorage
	the address of the storage vector for persistent representations
int	m_curRecLevel
	count recursive invocations, to detect recursion
bool	m_recursive
	if true, work in recursion-safe way (slower)
bool	m_ignoreRecursion
	if true, do not throw errors in case of recursion.
TPObjRef::typeID_t	m_pStorageTID
	TP Ref typeID for the persistent objects this converter is creating.
unsigned	m_pStorageTIDvalue
	m_pStorageTID converted to integer value
TopLevelTPCnvBase *	m_topConverter
	top level converter that owns this elemental TP converter it also holds the storage object
TopLevelTPCnvBase *	m_topConverterRuntime
	top level converter "owning" this TP converter at runtime (different from m_topConverter in case the top-level converter and object have extensions)
bool	m_wasUsedForReading
	flag set when using this converter for reading triggers search for a new converter before writing, to prevent possible use of old version

Detailed Description

Definition at line 16 of file TRT_HitCollectionCnv_p3.h.

Member Typedef Documentation

◆ Factory

typedef Gaudi::PluginService::Factory<ITPCnvBase*()> ITPCnvBase::Factory

inherited

Definition at line 26 of file ITPCnvBase.h.

◆ Pers_t

typedef PERS TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::Pers_t

inherited

Definition at line 335 of file TPConverter.h.

◆ PersBase_t

typedef PERS TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::PersBase_t

inherited

Definition at line 336 of file TPConverter.h.

◆ PolyCnvBase_t

template<class TRANS>

typedef ITPConverterFor< TRANS > ITPConverterFor< TRANS >::PolyCnvBase_t

inherited

Definition at line 41 of file TPConverter.h.

◆ Trans_t

typedef TRANS TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::Trans_t

inherited

Definition at line 334 of file TPConverter.h.

◆ TransBase_t

template<class TRANS>

typedef TRANS ITPConverterFor< TRANS >::TransBase_t

inherited

Definition at line 39 of file TPConverter.h.

Constructor & Destructor Documentation

◆ TRT_HitCollectionCnv_p3()

TRT_HitCollectionCnv_p3::TRT_HitCollectionCnv_p3 ( )

inline

Definition at line 20 of file TRT_HitCollectionCnv_p3.h.

20{};

Member Function Documentation

◆ baseToPersistent()

template<class TRANS>

template<class CNV>

TPObjRef ITPConverterFor< TRANS >::baseToPersistent	(	CNV **	cnv,
		const typename CNV::Trans_t *	transObj,
		MsgStream &	log ) const

inlineinherited

Persistify bass class of a given object and store the persistent represenation in the storage vector of the top-level persistent object.

The converter is located using the transient type from the CNV parameter, not from the object itself (because we need the base type, not the actual type)

Parameters

cnv	[IN/OUT] type of this parameter decides which converter will be used. Once the converter is found, this pointer will be set so the search is done only once
transObj	[IN] transient object
log	[IN] output message stream

Returns: TPObjRef TP reference to the persistent representation stored in the storage vector of the top-level persistent object

Definition at line 97 of file TPConverter.h.

                                                                                           {
     if( !*cnv || (*cnv)->wasUsedForReading() ) {
        // don't trust the converter if it was used for reading, find again
    *cnv = converterForType( *cnv, typeid(typename CNV::Trans_t), log );
    if( !*cnv )  return TPObjRef();
        (*cnv)->clearReadingFlag();
     }
//     return (**cnv).toPersistent_impl(transObj, log); 
     return (**cnv).virt_toPersistent(transObj, log); 
  }

◆ clearReadingFlag()

template<class TRANS>

void ITPConverterFor< TRANS >::clearReadingFlag ( )

inlineinherited

Definition at line 235 of file TPConverter.h.

235{ m_wasUsedForReading = false; }

ITPConverterFor::m_wasUsedForReading

bool m_wasUsedForReading

flag set when using this converter for reading triggers search for a new converter before writing,...

Definition TPConverter.h:306

◆ converterForRef()

template<class TRANS>

template<class CNV>

CNV * ITPConverterFor< TRANS >::converterForRef	(	CNV *	cnv,
		const TPObjRef &	ref,
		MsgStream &	log ) const

inlineinherited

Find converter for a TP type ID (passed in a TP Ref), that is or ihnerits from CNV type.

Parameters

cnv	[IN] parameter specifying the converter type
ref	[IN] TP Ref to an object for which a converter is sought
log	[IN] output message stream

Returns: CNV* pointer to the converter, if found

Definition at line 74 of file TPConverter.h.

                                                                         {
     ITPConverter *c = m_topConverterRuntime->converterForRef( ref );
     cnv = dynamic_cast<CNV*>(c);
     if( !cnv )
        this->converterNotFound( ref.typeID(), c, typeid(CNV).name(), log );
     return cnv;
  }

◆ converterForType()

template<class TRANS>

template<class CNV>

CNV * ITPConverterFor< TRANS >::converterForType	(	CNV *	cnv,
		const std::type_info &	t_info,
		MsgStream &	log ) const

inlineinherited

Find converter for a given C++ type ID, that is or ihnerits from CNV type.

Parameters

cnv	[IN] parameter specifying the converter type
t_info	[IN] C++ type id for which a converter is sought
log	[IN] output message stream

Returns: CNV* pointer to the converter, if found

Definition at line 58 of file TPConverter.h.

                                                                                 {
     ITPConverter *c = m_topConverterRuntime->converterForType( t_info );
     cnv = dynamic_cast< CNV* >( c );
     if( !cnv )
    this->converterNotFound( typeid(CNV), c, t_info.name(), log );
     return cnv;
  }

◆ converterNotFound() [1/2]

void ITPConverter::converterNotFound	(	const std::type_info &	converterType,
		ITPConverter *	c,
		const std::string &	typeName,
		MsgStream &	log ) const

virtualinherited

method called when the right TP converter was not found during writing

useful as a debugging hook, prints a detailed error message

Parameters

converterType	[IN] converterType that was requested
c	[IN] converter that was actually found (0 if not)
typeName	[IN] the C++ type name of the type for which converter was searched for
log	[IN] output message stream

Definition at line 22 of file ITPConverter.cxx.

{
   log << MSG::ERROR << ">>>>>> in parent TP converter " << typeid(*this).name()
       << ": could not find matching TP converter for type " << typeName << endmsg;
   if( c ) {
      log << MSG::ERROR << " - found incompatible converter of type "
      << typeid(*c).name() << endmsg;
   }
   log << MSG::ERROR << "  Converter handle type was " << converterType.name() << endmsg;
   errorHandler();
}

◆ converterNotFound() [2/2]

void ITPConverter::converterNotFound	(	unsigned	typeID,
		ITPConverter *	c,
		const std::string &	typeName,
		MsgStream &	log ) const

virtualinherited

method called when the right TP converter was not found during reading

useful as a debugging hook, prints a detailed error message

Parameters

typeID	[IN] converter ID that was requested
c	[IN] converter that was actually found (0 if not)
typeName	[IN] the C++ type name of the type for which converter was searched for
log	[IN] output message stream

Definition at line 40 of file ITPConverter.cxx.

{
   log << MSG::ERROR << ">>>>>> in parent TP converter " <<  typeid(*this).name()
       << "   requested TP converter for TP type ID " << typeID << " not found " << endmsg;
   if( c ) {
      log << MSG::ERROR << " - found converter " <<  typeid(*c).name()
          << " for " << c->transientTInfo().name()
          << " with an incompatible base type " << c->transBaseTInfo().name()
          << endmsg;
   }
   log << MSG::ERROR << "  Converter handle type was " << reqCnvTypeName <<  endmsg;
   errorHandler();
}

◆ createPersistent()

virtual PERS * TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::createPersistent	(	const TRANS *	transObj,
		MsgStream &	log )

virtualinherited

Create persistent representation of a transient object.

Simply creates a new persistent object and calls transToPers()

Parameters

transObj	[IN] transient object
log	[IN] output message stream

Returns: the created persistent representation

◆ createPersistentWithKey()

virtual PERS * TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::createPersistentWithKey	(	const TRANS *	transObj,
		const std::string &	key,
		MsgStream &	log )

virtualinherited

Create persistent representation of a transient object, with SG key.

Simply creates a new persistent object and calls transToPersWithKey()

Parameters

transObj	[IN] transient object
key	[IN] SG key of object being written
log	[IN] output message stream

Returns: the created persistent representation

◆ createTransFromPStore()

template<class TRANS>

template<class CNV>

CNV::Trans_t * ITPConverterFor< TRANS >::createTransFromPStore	(	CNV **	cnv,
		const TPObjRef &	ref,
		MsgStream &	log ) const

inlineinherited

Create transient representation of a persistent object, stored in the the top-level persistent object and referenced by the TP Ref.

If a TP converter is not specified, it will be found based on the Ref type.

Parameters

cnv	[IN][OUT] pointer to the converter, usually 0 at the start. Once the right converter is found, this pointer will be set so the search is done only once
ref	[IN] TP Ref to the persistent object to be converted
log	[IN] output message stream

Returns: pointer to the created transient represention

Definition at line 172 of file TPConverter.h.

                                                                               {
     if( ref.isNull() )  return 0;
     CNV *temp_cnv_p = 0;
     if( !cnv ) cnv = &temp_cnv_p;
     // see if we already have a converter and if it is the right one
     if( !*cnv || (*cnv)->typeID().value() != ref.typeID() ) {
        // we don't - find the right converter for ref.typeID()
    *cnv = converterForRef( *cnv, ref, log );
    if( !*cnv ) return 0;
        (*cnv)->setReadingFlag();
     }
     return (**cnv).virt_createTransFromPStore( ref.index(), log );
  }

◆ createTransient() [1/2]

virtual TRANS * TPPolyCnvBase< TRANS, TRANS, PERS >::createTransient	(	const PERS *	persObj,
		MsgStream &	log )

virtualinherited

Create transient representation of a persistent object.

Simply creates a new transient object and calls persToTrans()

Parameters

persObj	[IN] persistent object
log	[IN] output message stream

Returns: the created transient object

◆ createTransient() [2/2]

TRTUncompressedHitCollection * TRT_HitCollectionCnv_p3::createTransient	(	const TRT_HitCollection_p3 *	persObj,
		MsgStream &	log )

virtual

Definition at line 281 of file TRT_HitCollectionCnv_p3.cxx.

                                                                                                                          {
  std::unique_ptr<TRTUncompressedHitCollection> trans(std::make_unique<TRTUncompressedHitCollection>("DefaultCollectionName",persObj->m_nHits.size()));
  persToTrans(persObj, trans.get(), log);
  return(trans.release());
} //createTransient

◆ createTransientWithKey()

virtual TRANS * TPPolyCnvBase< TRANS, TRANS, PERS >::createTransientWithKey	(	const PERS *	persObj,
		const std::string &	key,
		MsgStream &	log )

virtualinherited

Create transient representation of a persistent object, with SG key.

Simply creates a new transient object and calls persToTransWithKey()

Parameters

persObj	[IN] persistent object
key	[IN] SG key of object being read
log	[IN] output message stream

Returns: the created transient object

◆ fillTransFromPStore()

template<class TRANS>

template<class CNV, class TRANS_T>

void ITPConverterFor< TRANS >::fillTransFromPStore	(	CNV **	cnv,
		const TPObjRef &	ref,
		TRANS_T *	trans,
		MsgStream &	log ) const

inlineinherited

Convert persistent object, stored in the the top-level persistent object and referenced by the TP Ref, to transient representation.

An empty transient object to be filled in is provided. If converter is not given, it will be found based on the Ref type.

Parameters

cnv	[IN][OUT] pointer to the converter, usually 0 at the start. Once the right converter is found, this pointer will be set so the search is done only once
ref	[IN] TP Ref to the persistent object to be converted
trans	[IN] pointer to the empty transient object
log	[IN] output message stream

Definition at line 145 of file TPConverter.h.

                                                                                                    {
    if( ref.isNull() )  return;
    CNV *temp_cnv_p = 0;
    if( !cnv ) cnv = &temp_cnv_p;
    // see if we already have a converter and if it is the right one
    if( !*cnv || (*cnv)->typeID().value() != ref.typeID() ) {
       // we don't - find the right converter for ref.typeID()
       *cnv = converterForRef( *cnv, ref, log );
       if( !*cnv ) return;
       (*cnv)->setReadingFlag();
     }
    (**cnv).pstoreToTrans( ref.index(), trans, log );
  }

◆ ignoreRecursion()

void TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::ignoreRecursion ( bool flag = false )

inlineinherited

Tell the converter to ignore recursion (do not throw errors) even when recurion is detected.

UNSAFE! use only if you are sure you preallocated enough persistent storage

Definition at line 568 of file TPConverter.h.

                                           {
     m_ignoreRecursion = flag;
  }

◆ initPrivateConverters()

template<class TRANS>

virtual void ITPConverterFor< TRANS >::initPrivateConverters ( TopLevelTPCnvBase * )

inlinevirtualinherited

Reimplemented in JetCnv_p4, MeasuredAtaSurfaceCnv_p1< SURFACE_CNV, ATA_SURFACE >, MeasuredAtaSurfaceCnv_p1< AtaCylinderCnv_p1, Trk::AtaCylinder >, MeasuredAtaSurfaceCnv_p1< AtaDiscCnv_p1, Trk::AtaDisc >, MeasuredAtaSurfaceCnv_p1< AtaPlaneCnv_p1, Trk::AtaPlane >, MeasuredAtaSurfaceCnv_p1< AtaStraightLineCnv_p1, Trk::AtaStraightLine >, MeasuredAtaSurfaceCnv_p1< CurvilinearParametersCnv_p1, Trk::CurvilinearParameters >, MeasuredNeutralAtaSurfaceCnv_p1< SURFACE_CNV, ATA_SURFACE >, MeasuredNeutralAtaSurfaceCnv_p1< NeutralAtaCylinderCnv_p1, Trk::NeutralAtaCylinder >, MeasuredNeutralAtaSurfaceCnv_p1< NeutralAtaDiscCnv_p1, Trk::NeutralAtaDisc >, MeasuredNeutralAtaSurfaceCnv_p1< NeutralAtaPlaneCnv_p1, Trk::NeutralAtaPlane >, MeasuredNeutralAtaSurfaceCnv_p1< NeutralAtaStraightLineCnv_p1, Trk::NeutralAtaStraightLine >, SegmentCnv_p1, TauPi0CandidateCnv_p1, TauPi0ClusterCnv_p1, TauPi0DetailsCnv_p1, TauPi0DetailsCnv_p2, TrackCnv_p12, TrackCnv_p1, TrackCnv_p2, TrackCnv_p3, TrackCnv_p4, TrackParticleBaseCnv_p1, TrackParticleCnv_p2, TrackParticleCnv_p3, TrigInDetTrackCollectionCnv_p1, TrigVertexCnv_p1, TrigVertexCnv_p2, V0CandidateCnv_p1, VxCandidateCnv_p1, and VxCandidateCnv_p2.

Definition at line 187 of file TPConverter.h.

187{}

◆ persistentTInfo()

virtual const std::type_info & TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::persistentTInfo ( ) const

inlinevirtualinherited

return C++ type id of the persistent class this converter is for

Returns: std::type_info&

Implements ITPCnvBase.

Definition at line 482 of file TPConverter.h.

482{ return typeid(PERS); }

◆ persToTrans() [1/2]

virtual void TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::persToTrans	(	const PERS *	persObj,
		TRANS *	transObj,
		MsgStream &	log )

pure virtualinherited

Convert persistent representation to transient one.

Copies data members from persistent object to an existing transient one. Needs to be implemented by the developer on the actual converter.

Parameters

persObj	[IN] persistent object
transObj	[IN] transient object
log	[IN] output message stream

◆ persToTrans() [2/2]

void TRT_HitCollectionCnv_p3::persToTrans	(	const TRT_HitCollection_p3 *	persCont,
		TRTUncompressedHitCollection *	transCont,
		MsgStream &	log )

virtual

Definition at line 289 of file TRT_HitCollectionCnv_p3.cxx.

{
 
  //    if (log.level() <= MSG::DEBUG) log << MSG::DEBUG << "In TRT_HitCollectionCnv_p3::persToTrans()" << endmsg;
 
  // some values are read less than once per hit, these need counters.
  unsigned int meanTimeofCount=0, startRCount=0, endRCount=0, hitEneCount=0;
  unsigned int idxBC=0, idxId=0, endHit=0, endBC=0, endId=0;
 
  // Assume that all Hits should be linked to the hard-scatter GenEvent
  IProxyDict *sg = SG::CurrentEventStore::store();
  const int event_number = HepMcParticleLink::getEventNumberAtPosition (0, sg);
 
  //
  // loop over strings - index [i]
  //
 
  for ( unsigned int i = 0; i < persCont->m_nHits.size(); i++ ) {
 
    if ( persCont->m_nHits[i] ) { // at least one hit in the string
 
      const unsigned int startHit = endHit;
      endHit += persCont->m_nHits[i];
 
      //
      // string strawId
      //
      const unsigned int i1 = persCont->m_strawId1b[i]; //  8 bits
      const unsigned int i2 = persCont->m_strawId2b[i]; // 16 bits
      const unsigned int strawId = i2*256+i1;           // => 24 bits (0 to 16,777,215)
 
      //
      // string startPhi
      //
      const unsigned int istartPhi = persCont->m_startPhi[i];         // 8 bits
      const double startPhi = -M_PI + (istartPhi+0.5)*2.0*M_PI/256.0; // rad (min = -pi, max = +pi)
 
      //
      // string startZ
      //
      const unsigned int istartZ = persCont->m_startZ[i] >> 1; // 4 bits
      double startZ = -365.0*CLHEP::mm + (istartZ+0.5)*730.0*CLHEP::mm/16.0; // (min = -365 mm, max = +365 mm)
 
      //
      // start Rflag
      //
      const unsigned int istartRflag = persCont->m_startZ[i] & 1; // 1 bit
 
      //
      // string startR
      //
      double startR;
      if ( istartRflag == 1 ) {
        startR = 2.0*CLHEP::mm; // 1 bit
      }
      else {
        const unsigned int istartR = persCont->m_startR[startRCount++]; // 8 bits
        startR = (istartR+0.5)*2.0*CLHEP::mm/256.0;                            // (range 0 - 2 mm)
        if ( startR < 0.0155*CLHEP::mm ) startR = 0.0155*CLHEP::mm;                   // The wire radius
      }
 
      //
      // string startX, startY (derived from R,Phi)
      //
      double startX = startR*cos(startPhi);
      double startY = startR*sin(startPhi);
      /*
      // Validation output
      std::cout.precision(15);
      std::cout << "AJBPtoTstrawId "    << strawId              << std::endl;
      std::cout << "AJBPtoTstartR "     << startR               << std::endl;
      std::cout << "AJBPtoTstartPhi "   << startPhi             << std::endl;
      std::cout << "AJBPtoTstartX "     << startX               << std::endl;
      std::cout << "AJBPtoTstartY "     << startY               << std::endl;
      std::cout << "AJBPtoTstartZ "     << startZ               << std::endl;
      std::cout << "AJBPtoTnHits "      << persCont->m_nHits[i] << std::endl;
      */
      //
      // loop over end hits in the string - index [j]
      //
 
      for ( unsigned int j = startHit; j < endHit; j++ ) {
 
        if ( j >= endBC + persCont->m_nBC[idxBC] ) endBC += persCont->m_nBC[idxBC++];
        if ( j >= endId + persCont->m_nId[idxId] ) endId += persCont->m_nId[idxId++];
 
        //
        // hit meanTime
        //
        const unsigned int imeanTime = persCont->m_meanTime[j] >> 2; // 10 bits
        double meanTime = (imeanTime+0.5)*75.0*CLHEP::ns/1024.0;            // (min = 0.0 ns, max = 75.0 ns)
        if ( imeanTime == 1023 ) meanTime = (double)persCont->m_meanTimeof[meanTimeofCount++]; // ns, 32-bit float overflow
 
        //
        // dZ sign
        //
        const unsigned int idZsign = (persCont->m_meanTime[j] >> 1 ) & 1; // 1 bit
 
        //
        // endR flag
        //
        const unsigned int iendRflag =  persCont->m_meanTime[j] & 1;         // 1 bit
 
        //
        // hit energy deposited in keV (only relevant for photons) 32-bit float
        //
        const double hitEne = ( persCont->m_id[idxId] == 22 ||
                                (int)(abs(persCont->m_id[idxId])/100000) == 41 ||
                                (int)(abs(persCont->m_id[idxId])/10000000) == 1
                                ) ? (double)persCont->m_hitEne[hitEneCount++] : 0.0;
 
        //
        // hit endPhi (can be modified later during "steplength preservation")
        //
        const unsigned int iendPhi =  persCont->m_endPhi[j];  // 8 bits
        double endPhi = -M_PI + (iendPhi+0.5)*2.0*M_PI/256.0; // rad (min = -pi, max = +pi)
 
        //
        // string endR (can be modified later during "steplength preservation")
        //
        double endR;
        if ( iendRflag==1 ) {
          endR = 2.0*CLHEP::mm;  // 1 bit
        }
        else {
          const unsigned int iendR = persCont->m_endR[endRCount++];
          endR = (iendR+0.5)*2.0*CLHEP::mm/256.0;          // 8 bits
          if ( endR < 0.0155*CLHEP::mm ) endR = 0.0155*CLHEP::mm; // the wire radius
        }
 
        //
        // hit endX, endY (derived from R,Phi)
        //
        double endX = endR*cos(endPhi); // can be modified later during "steplength preservation"
        double endY = endR*sin(endPhi); // can be modified later during "steplength preservation"
 
        // Save the (o)riginal endX, endY values for the next hit start because
        // they might get shrunk to fit the g4 steplength of the current hit.
        double endXo = endX;
        double endYo = endY;
 
        //
        // g4 step length of the hit, m_steplength, and
        // kinetic energy of the hit, m_kinEne, are both 15-bit short floats.
        // Note: a rare condition causes a 16-bit short float (mantissa=512).
        //
        const int kmantissa = persCont->m_kinEne[j]     >> 6;   // 9 bits (expected)
        const int smantissa = persCont->m_steplength[j] >> 6;
        const int kexponent = persCont->m_kinEne[j]     & 0x3F; // 6 bits
        const int sexponent = persCont->m_steplength[j] & 0x3F;
        const double kinEne = (kmantissa+512.5)/1024 * pow(2.0,kexponent) / 1.0e9; // MeV
        double g4steplength = (smantissa+512.5)/1024 * pow(2.0,sexponent) / 1.0e9; // mm
        if ( idZsign==0 ) g4steplength = -g4steplength;
 
        //
        // Preserving the steplength of the hit by setting endZ or shrinking dX,dY.
        //
        double dX = endX-startX;
        double dY = endY-startY;
        double dZ;
        double dXY2 = dX*dX+dY*dY;
        double dL2 = g4steplength*g4steplength;
        if ( dL2 > dXY2 ) {         // define dZ such that steplength = g4steplength
          dZ = sqrt(dL2-dXY2);
          if (g4steplength<0.0) dZ=-dZ;
        }
        else {    // dL2 < dXY2   // shrink dX,dY such that dXY = g4steplength
          dX = dX * sqrt(dL2/dXY2); // this includes the cases where dL2=0!
          dY = dY * sqrt(dL2/dXY2);
          dZ = 0.0*CLHEP::mm;
          endX = startX + dX;
          endY = startY + dY;
          //endR = sqrt( endX*endX + endY*endY ); // for validation information
          //endPhi = atan2(endY,endX);            // for validation information
        }
        double endZ = startZ + dZ;
        //dX = endX-startX; // for validation information
        //dY = endY-startY; // for validation information
        /*
        // Validation output
        std::cout.precision(15);
        std::cout << "AJBPtoTendR "     << endR     << std::endl;
        std::cout << "AJBPtoTendPhi "   << endPhi   << std::endl;
        std::cout << "AJBPtoTendX "     << endX     << std::endl;
        std::cout << "AJBPtoTendY "     << endY     << std::endl;
        std::cout << "AJBPtoTendZ "     << endZ     << std::endl;
        std::cout << "AJBPtoTmeanTime " << meanTime << std::endl;
        std::cout << "AJBPtoTkinEne "   << kinEne   << std::endl;
        std::cout << "AJBPtoThitEne "   << hitEne   << std::endl;
        std::cout << "AJBPtoTsteplength " << sqrt(dX*dX+dY*dY+dZ*dZ) << std::endl;
        */
        //
        // Notes:
        // - All units are CLHEP, except hitEne which is in keV.
        // - For charged particles kinEne is *zero*!
        //
 
        HepMcParticleLink partLink(persCont->m_barcode[idxBC], event_number, HepMcParticleLink::IS_EVENTNUM, HepMcParticleLink::IS_BARCODE, sg);
        if ( HepMC::BarcodeBased::is_truth_suppressed_pileup(static_cast<int>(persCont->m_barcode[idxBC])) ) {
          partLink.setTruthSuppressionType(EBC_PU_SUPPRESSED);
        }
        transCont->Emplace( strawId, partLink, persCont->m_id[idxId],
                            kinEne, hitEne, startX, startY, startZ,
                            endX, endY, endZ, meanTime );
        //
        // End of this hit becomes the start of the next;
        // use the original (uncorrected) values for X,Y
        // but the derived value for Z.
        //
        startX = endXo; startY = endYo; startZ = endZ;
 
      }
    } // nhits>0
  } // straw loop
} // persToTrans

◆ persToTransUntyped()

virtual void TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::persToTransUntyped	(	const void *	pers,
		void *	trans,
		MsgStream &	log )

inlinevirtualinherited

Convert persistent object representation to transient.

Parameters

pers	[IN] void* pointer to the persistent object
trans	[OUT] void* pointer to the empty transient object
log	[IN] output message stream

Implements ITPCnvBase.

Definition at line 400 of file TPConverter.h.

  {
    persToTrans (reinterpret_cast<const PERS*> (pers),
                 reinterpret_cast<TRANS*> (trans),
                 log);
  }

◆ persToTransWithKey()

virtual void TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::persToTransWithKey	(	const PERS *	persObj,
		TRANS *	transObj,
		const std::string &	,
		MsgStream &	log )

inlinevirtualinherited

Convert persistent representation to transient one.

Copies data members from persistent object to an existing transient one. Needs to be implemented by the developer on the actual converter.

Parameters

persObj	[IN] persistent object
transObj	[IN] transient object
log	[IN] output message stream

Reimplemented in AthExParticlesCnv_p1, CaloCellContainerCnv_p1, CaloCellLinkContainerCnv_p1, CaloCellLinkContainerCnv_p2, CaloClusterCellLinkContainerCnv_p1, TPConverterWithKeyBase< TRANS, PERS >, and xAODTauJetAuxContainerCnv_v2.

Definition at line 376 of file TPConverter.h.

  {
    return persToTrans (persObj, transObj, log);
  }

◆ persToTransWithKeyUntyped()

virtual void TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::persToTransWithKeyUntyped	(	const void *	pers,
		void *	trans,
		const std::string &	key,
		MsgStream &	log )

inlinevirtualinherited

Convert persistent object representation to transient.

Parameters

pers	[IN] void* pointer to the persistent object
trans	[OUT] void* pointer to the empty transient object
key	[IN] SG key of object being read.
log	[IN] output message stream

Reimplemented from ITPCnvBase.

Definition at line 420 of file TPConverter.h.

  {
    persToTransWithKey (reinterpret_cast<const PERS*> (pers),
                        reinterpret_cast<TRANS*> (trans),
                        key,
                        log);
  }

◆ pstoreToTrans()

template<class TRANS, class PERS>

virtual void TPConverterBase< TRANS, PERS >::pstoreToTrans	(	unsigned	index,
		TRANS *	trans,
		MsgStream &	log )

inlinevirtualinherited

Convert persistent representation stored in the storage vector of the top-level object to transient.

Internal.

Parameters

index	[IN] index of the persistent representation in the storage vector
trans	[IN] empty transient object
log	[IN] output message stream

Reimplemented from TPAbstractPolyCnvBase< TRANS, TRANS, PERS >.

Definition at line 760 of file TPConverter.h.

                                                                              {
     assert (index < this->m_pStorage->size());
     this->persToTrans( &(*this->m_pStorage)[index], trans, log );
  }

◆ reservePStorage()

virtual void TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::reservePStorage ( size_t size )

inlinevirtualinherited

Reserve 'size' elements for persistent storage.

Implements ITPConverter.

Definition at line 573 of file TPConverter.h.

                                               {
     m_pStorage->reserve( size );
  }

◆ setPStorage()

void TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::setPStorage ( std::vector< PERS > * storage )

inlineinherited

Tell this converter which storage vector it should use to store or retrieve persistent representations.

Parameters

storage [IN] the address of the storage vector

Definition at line 551 of file TPConverter.h.

                                                  {
     m_pStorage = storage;
     m_curRecLevel = 0;
  }

◆ setReadingFlag()

template<class TRANS>

void ITPConverterFor< TRANS >::setReadingFlag ( )

inlineinherited

Definition at line 234 of file TPConverter.h.

234{ m_wasUsedForReading = true; }

◆ setRecursive()

void TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::setRecursive ( bool flag = true )

inlineinherited

Tell the converter if it should work in recursive mode slower but it can safely handle recursion.

Definition at line 559 of file TPConverter.h.

                                       {
     m_recursive = flag;
  }

◆ setRuntimeTopConverter()

template<class TRANS>

virtual void ITPConverterFor< TRANS >::setRuntimeTopConverter ( TopLevelTPCnvBase * topConverter )

inlinevirtualinherited

Set runtime top-level converter - usually it is the owning TL converter, but in case of extended objects it will be the TL converter of the extended object.

Parameters

topConverter [IN] runtime top-level converter for this converter

Implements ITPConverter.

Definition at line 215 of file TPConverter.h.

                                                                          {
     m_topConverterRuntime = topConverter;
     initPrivateConverters( topConverter );
  }

◆ setTopConverter()

template<class TRANS>

virtual void ITPConverterFor< TRANS >::setTopConverter	(	TopLevelTPCnvBase *	topConverter,
		const TPObjRef::typeID_t &	TPtypeID )

inlinevirtualinherited

Set which top-level converter owns this elemental converter, and what TPtypeID was assigned to the persistent objects it produces.

Parameters

topConverter	[IN] the top-level converter owning this converter
TPtypeID	[IN] TP type id for persistent objects (used in TP refs)

Implements ITPConverter.

Definition at line 221 of file TPConverter.h.

  {
     m_topConverter = m_topConverterRuntime = topConverter;
     m_pStorageTID = TPtypeID;
     m_pStorageTIDvalue = TPtypeID.value();
     initPrivateConverters( topConverter );
  }

◆ topConverter() [1/2]

template<class TRANS>

virtual TopLevelTPCnvBase * ITPConverterFor< TRANS >::topConverter ( )

inlinevirtualinherited

return the top-level converter for this elemental TP converter

Returns: TopLevelTPCnvBas

Reimplemented from ITPConverter.

Definition at line 191 of file TPConverter.h.

                                               {
     return m_topConverter;
  }

◆ topConverter() [2/2]

template<class TRANS>

virtual const TopLevelTPCnvBase * ITPConverterFor< TRANS >::topConverter ( ) const

inlinevirtualinherited

return the top-level converter for this elemental TP converter

Returns: TopLevelTPCnvBas

Reimplemented from ITPConverter.

Definition at line 196 of file TPConverter.h.

                                                         {
     return m_topConverter;
  }

◆ toPersistent()

template<class TRANS>

template<class CNV>

TPObjRef ITPConverterFor< TRANS >::toPersistent	(	CNV **	cnv,
		const typename CNV::TransBase_t *	transObj,
		MsgStream &	log ) const

inlineinherited

Persistify an object and store the persistent represenation in the storage vector of the top-level persistent object.

The correct converter is located using the actual object type.

Parameters

cnv	[IN/OUT] pointer to the converter, usually 0 at the start. Once the right converter is found, this pointer will be set so the search is done only once
transObj	[IN] transient object
log	[IN] output message stream

Returns: TPObjRef TP reference to the persistent representation stored in the storage vector of the top-level persistent object

Definition at line 119 of file TPConverter.h.

                                                                                                   {
     if( !transObj ) return TPObjRef();
     CNV *temp_cnv_p = 0;
     if( !cnv ) cnv = &temp_cnv_p;
     if( !*cnv || (*cnv)->wasUsedForReading() ) {
        // don't trust the converter if it was used for reading, find again
    *cnv = converterForType( *cnv, typeid(*transObj), log );
    if( !*cnv )  return TPObjRef();
        (*cnv)->clearReadingFlag();
     }
     return (**cnv).virt_toPersistent(transObj, log);
  }

◆ toPersistentWithKey_impl()

TPObjRef TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::toPersistentWithKey_impl	(	const TRANS *	trans,
		const std::string &	key,
		MsgStream &	log )

inherited

Convert transient object to persistent representation.

Stores the result in the storage vector of the top-level object and returns a TP Ref to it.

Parameters

trans	[IN] transient object
key	[IN] SG key of object being converted
log	[IN] output message stream

Returns: TP reference to the persistent representation

◆ transBaseTInfo()

template<class TRANS>

const std::type_info & ITPConverterFor< TRANS >::transBaseTInfo ( ) const

inlinevirtualinherited

return C++ type id of the common base transient type for all converters for a group of polymorphic types

Returns: std::type_info& this method is not overwritten in the subclasses like transientTInfo()

Implements ITPConverter.

Definition at line 205 of file TPConverter.h.

205{ return typeid(TRANS); }

◆ transientTInfo()

virtual const std::type_info & TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::transientTInfo ( ) const

inlinevirtualinherited

return C++ type id of the transient class this converter is for

Returns: std::type_info&

Reimplemented from ITPConverterFor< TRANS >.

Definition at line 479 of file TPConverter.h.

479{ return typeid(TRANS); }

◆ transToPers() [1/2]

virtual void TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::transToPers	(	const TRANS *	transObj,
		PERS *	persObj,
		MsgStream &	log )

pure virtualinherited

Convert transient representation to persistent one.

Copies data members from transient object to an existing persistent one. Needs to be implemented by the developer on the actual converter.

Parameters

transObj	[IN] transient object
persObj	[IN] persistent object
log	[IN] output message stream

◆ transToPers() [2/2]

void TRT_HitCollectionCnv_p3::transToPers	(	const TRTUncompressedHitCollection *	transCont,
		TRT_HitCollection_p3 *	persCont,
		MsgStream &	log )

virtual

Definition at line 24 of file TRT_HitCollectionCnv_p3.cxx.

{
 
  /*
    Spring 2009
    Andrew Beddall - lossy TRT G4hit compression [p3]
 
    In p1, p2 versions, GEANT hits are persistified on disk as floats.
    In this p3 version, floats are compressed to "integers"/"short-floats" before persistifying.
    The saving is about 75%; see http://cern.ch/beddall/TRThitCompression/
 
    Spring 2008
    Rob Duxfield - lossless TRT G4hit compression [p2]
 
    Finds hits belonging to a "string" (in which the end point of one hit is
    the same as the start point of the next) and persistifies the end point
    of each hit plus the start point of the first hit in each string.
  */
 
  //  The original units from the hit simulation are indicated in comments;
  //  they are all in CLHEP units except for hitEne which is in keV.
  //  I sometimes make use of CLHEP scales *CLHEP::mm and *CLHEP::ns (both=1) for clarity (I hope!).
  //  See also https://twiki.cern.ch/twiki/bin/view/Atlas/TrtSoftware#Production_of_Hits
 
  static const double dRcut = 1.0e-7*CLHEP::mm;
  static const double dTcut = 1.0*CLHEP::ns; // redundant?
 
  //    if (log.level() <= MSG::DEBUG) log << MSG::DEBUG << "In TRT_HitCollectionCnv_p3::transToPers()" << endmsg;
 
  int lastBarcode = -1;
  int lastId = -1;
  double lastT = 0.0*CLHEP::ns;
  unsigned int idx = 0;
  unsigned int endBC = 0;
  unsigned int endId = 0;
  unsigned int endHit = 0;
  HepGeom::Point3D<double> lastEnd(0.0, 0.0, 0.0); // mm
 
  for (TRTUncompressedHitCollection::const_iterator it = transCont->begin(); it != transCont->end(); ++it) {
 
    TRTUncompressedHitCollection::const_iterator trtHit = it;
 
    if ( trtHit->particleLink().barcode() != lastBarcode  ||  idx - endBC > 65500) {   // max unsigned short =  65535;
      // store barcode once for set of consecutive hits with same barcode
      lastBarcode = trtHit->particleLink().barcode();
      using barcodeType = decltype(persCont->m_barcode)::value_type;
      persCont->m_barcode.push_back(static_cast<barcodeType>(lastBarcode));
      if ( idx > 0 ) {
        persCont->m_nBC.push_back(idx - endBC);
        endBC = idx;
      }
    }
 
    if ( (int)trtHit->GetParticleEncoding() != lastId || idx - endId >  65500) { // max unsigned short =  65535;
      // store id once for set of consecutive hits with same id
      lastId = trtHit->GetParticleEncoding();
      persCont->m_id.push_back(lastId);
      if ( idx > 0 ) {
        persCont->m_nId.push_back(idx - endId);
        endId = idx;
      }
    }
 
    const HepGeom::Point3D<double> hitStart(trtHit->GetPreStepX(), trtHit->GetPreStepY(), trtHit->GetPreStepZ()); // mm
 
    const double meanTime = trtHit->GetGlobalTime(); // ns  // Time of flight from the I.P. to the center of the hit.
    const double dTLast = fabs(meanTime - lastT);      // |d(meantime)| between the previous hit and the current one.
    const double dRLast = lastEnd.distance(hitStart);  // Distance between end of previous hit and start of current one;
    // this is zero if the hit is a continuation of the same particle in the same straw.
 
    // Begin a new string if the current and previous hits are disconnected;
    // it looks like dTcut is redundant (but not sure about this).
    if ( dRLast >= dRcut || dTLast >= dTcut ) {
 
      // if ( dRLast < dRcut) std::cout << "AJBdTLastTriggeredNewString " << dRLast << " " << dTLast << std::endl;
      // new hit string //
 
      //
      // Persistify string *strawId* using 24 bits.
      // Assumes 0 <= strawId <= 16,777,215 (strawId appears to be < 4,000,000)
      //
      const unsigned int strawId = trtHit->GetHitID();
      persCont->m_strawId1b.push_back(  (unsigned char)(strawId % 256) ); //  8 bits
      persCont->m_strawId2b.push_back( (unsigned short)(strawId / 256) ); // 16 bits
      if ( strawId>16777215 )
        log << MSG::WARNING << "TRT_HitCollectionCnv: strawId > 2^24-1 cannot be persistified correctly! " << endmsg;
 
      //
      // Persistify string start radius using 1 bit (istartRflag) or 8 bits (startR)
      // Note that the smallest value of R is the wire radius (0.0155 mm)
      //
      // R will be flagged as 2 mm if it is within 0.1 um of the straw wall => max error = 0.1 um,
      // otherwise compress with 8 bits => max error = 3.9 um (0.078 ns), RMS error = 1.1 um (0.022 ns)
      //
      const double startR = sqrt( hitStart.x()*hitStart.x() + hitStart.y()*hitStart.y() ); // mm
      unsigned short istartRflag;
      if ( startR > 1.9999*CLHEP::mm ) {
        istartRflag=1; // persistify as a 1-bit flag
      }
      else {
        istartRflag=0; // compress to 8 bits with a span of 2 mm
        persCont->m_startR.push_back( (unsigned char)(startR/(2.0*CLHEP::mm)*256.0) );
      }
 
      //
      // Persistify string *startPhi* using 8 bits (min=-pi, max=+pi)
      // Max. error = 12 mrad (< 24 um, 0.48 ns); RMS error = 7 mrad (< 14 um, 0.28 ns)
      //
      const double startPhi = atan2( hitStart.y(), hitStart.x() ); // returns range -pi to +pi rad
      persCont->m_startPhi.push_back( (unsigned char)( (startPhi+M_PI)/(2.0*M_PI)*256.0 ) );
 
      //
      // Persistify *startZ* using a 4 bits (min = -365 mm, max= +365 mm)
      // Max. error = 25 mm (25e-3/(0.75c) = 0.111 ns * 2 reflect = 0.222 ns)
      // RMS  error = 14 mm (14e-3/(0.75c) = 0.062 ns * 2 reflect = 0.124 ns)
      // Also the 1-bit *istartRflag* is packed into this variable.
      //
      // Note:
      // In the digi code we need to allow for something like 22.5 mm outside straw.
      // Also because we have short straws,
      // short straws are about < +-180 mm, long straws are about < +-350 mm
      // The following compressions can give a large "out of straw" value;
      // *don't* use these: (2.0), 32.0, 128.0, 256.0.
 
      unsigned char istartZ = (unsigned char)( (hitStart.z()+365.0*CLHEP::mm)/(730.0*CLHEP::mm)*16.0 );
      istartZ = (istartZ << 1) | istartRflag;
      persCont->m_startZ.push_back( istartZ );
 
      if ( idx > 0 ) {
        persCont->m_nHits.push_back( idx - endHit );
        endHit = idx;
      }
      /*
      // Validation output
      std::cout.precision(15);
      std::cout << "AJBTtoPstrawId "    << strawId      << std::endl;
      std::cout << "AJBTtoPstartR "     << startR       << std::endl;
      std::cout << "AJBTtoPstartPhi "   << startPhi     << std::endl;
      std::cout << "AJBTtoPstartX "     << hitStart.x() << std::endl;
      std::cout << "AJBTtoPstartY "     << hitStart.y() << std::endl;
      std::cout << "AJBTtoPstartZ "     << hitStart.z() << std::endl;
      */
    } // end of "begin new hit string"
    // Now for the end hits //
 
    const HepGeom::Point3D<double> hitEnd(trtHit->GetPostStepX(), trtHit->GetPostStepY(), trtHit->GetPostStepZ()); // mm
    const HepGeom::Point3D<double> hitLength = (hitEnd - hitStart);
 
    //
    // Here both *kinEne* (kinetic energy of the particle causing the hit) and
    // *steplength* (g4hit length) are persistified using a 15-bit "short float"
    // (9 bit unsigned mantissa, 6 bit unsigned exponent).
    // This stores values in the range 0.51*2^0 = 0.51 to 1.00*2^63 = 9.2e18.
    // I enforce the limits 1.0 and 9.0e18; see below.
    // Max relative error = 0.0010, RMS = 0.0004
    //
    // Notes:
    //
    //  - G4 gives kinEne in MeV; I sometimes see values ~ 1e-7 MeV (100 meV) [float round-off?]
    //    So I multiply by 1e9 and store in units of meV => range 1.0 meV to 9.0e18 meV (9000 TeV!)
    //  - About 1 in 10000 hits have steplength ~ 1e-7 mm [float round-off?]
    //    so again I multiply by 1e9 and store in units of pm => range 1.0 pm to 9.0e18 pm (9000 km)
    //  - The mantissa has maximum 9 bits, the exponent has maximum 6 bits,
    //    Note: a rare condition causes an 10-bit mantissa (mantissa=512).
    //
    double kinEne = trtHit->GetKineticEnergy()    * 1.0e9;  // nano Mev = meV.
    double steplength = hitLength.distance() * 1.0e9;  // nano mm  = pm.
    if ( kinEne     < 1.0 )        kinEne=1.0;         // Keep  the value
    if ( steplength < 1.0 )    steplength=1.0;         // well within the
    if ( kinEne     > 9.0e18 )     kinEne=9.0e18;      // range  of   the
    if ( steplength > 9.0e18 ) steplength=9.0e18;      // short    float.
    const unsigned int kexponent = (unsigned int)ceil(log10(kinEne)/0.30102999566398);
    const unsigned int sexponent = (unsigned int)ceil(log10(steplength)/0.30102999566398);
    const unsigned int kmantissa = (unsigned int)(kinEne/pow(2.0,kexponent)*1024) - 512;
    const unsigned int smantissa = (unsigned int)(steplength/pow(2.0,sexponent)*1024) - 512;
    persCont->m_kinEne.push_back(     (kmantissa << 6) | kexponent );
    persCont->m_steplength.push_back( (smantissa << 6) | sexponent );
 
    //
    // Persistify hit end radius using 1 bit (iendRflag) or 8 bits (endR).
    // Note that the smallest value of R is the wire radius (0.0155 mm)
    //
    // R will be flagged as 2 mm if it is within 0.1 um of the straw wall => max error = 0.1 um,
    // otherwise compress with 8 bits. The errors are as for startR, but can increased greatly
    // after steplength preservation in PtoT.
    //
    const double endR = sqrt( hitEnd.x()*hitEnd.x() + hitEnd.y()*hitEnd.y() ); // mm
    unsigned short iendRflag;
    if ( endR > 1.9999*CLHEP::mm ) {
      iendRflag=1; // persistify as a 1-bit flag
    }
    else {
      iendRflag=0; // compress to 8 bits with a span of 2 mm
      persCont->m_endR.push_back( (unsigned char)(endR/(2.0*CLHEP::mm)*256.0) );
    }
 
    //
    // Persistify string *endPhi* using 8 bits (min=-pi, max=+pi)
    // The errors are as for startPhi, but are very different after steplength
    // preservation in PtoT.
    //
    const double endPhi = atan2( hitEnd.y(), hitEnd.x() ); // returns range -pi to +pi rad
    persCont->m_endPhi.push_back( (unsigned char)( (endPhi+M_PI)/(2.0*M_PI)*256.0 ) );
 
    //
    // Persistify hit *meanTime* using 10 bits (min=0.,span=75 ns)
    // with float overflow for meanTime >= 75ns (the tail of the distribution).
    // Max. error = 0.037 ns; RMS error = 0.021 ns.
    // Also the 1-bit *iendRflag* and 1-bit *idZsign* are packed into this variable.
    //
    unsigned short idZsign = (hitLength.z()>0.0) ? 1 : 0; // flag the sign of dZ
    unsigned short imeanTime = ( meanTime < 75.0*CLHEP::ns ) ? (unsigned short)(meanTime/(75.0*CLHEP::ns)*1024.0) : 1023;
    if ( imeanTime == 1023 ) persCont->m_meanTimeof.push_back( (float)meanTime ); // "overflow flag"
    imeanTime = (imeanTime << 2) | (idZsign << 1) | iendRflag;
    persCont->m_meanTime.push_back( imeanTime );
 
    //
    // Persistify hit *hitEne* (the energy deposited by the hit in keV) using a float but only for photons
    // (relatively very few of these). Digitisation does not use hitEne for charged particles.
    //
    if ( lastId == 22 ||
         (int)(abs(lastId)/100000) == 41 ||
         (int)(abs(lastId)/10000000) == 1
         ) persCont->m_hitEne.push_back( (float)(trtHit->GetEnergyDeposit()) );  // keV
 
    lastEnd = hitEnd;
    lastT = meanTime;
    ++idx;
    /*
    // Validation output
    std::cout.precision(15);
    std::cout << "AJBTtoPendR "     << endR                  << std::endl;
    std::cout << "AJBTtoPendPhi "   << endPhi                << std::endl;
    std::cout << "AJBTtoPendX "     << hitEnd.x()            << std::endl;
    std::cout << "AJBTtoPendY "     << hitEnd.y()            << std::endl;
    std::cout << "AJBTtoPendZ "     << hitEnd.z()            << std::endl;
    std::cout << "AJBTtoPmeanTime " << meanTime              << std::endl;
    std::cout << "AJBTtoPkinEne "   << trtHit->kineticEnergy << std::endl;
    std::cout << "AJBTtoPhitEne "   << trtHit->energyDeposit << std::endl;
    std::cout << "AJBTtoPsteplength " << hitLength.distance() << std::endl;
    */
  }
 
  persCont->m_nBC.push_back(idx - endBC);
  persCont->m_nId.push_back(idx - endId);
  persCont->m_nHits.push_back( idx - endHit );
 
} // transToPers

◆ transToPersUntyped()

virtual void TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::transToPersUntyped	(	const void *	trans,
		void *	pers,
		MsgStream &	log )

inlinevirtualinherited

Convert transient object representation to persistent.

Parameters

trans	[IN] void* pointer to the transient object
pers	[OUT] void* pointer to the empty persistent object
log	[IN] output message stream

Implements ITPCnvBase.

Definition at line 410 of file TPConverter.h.

  {
    transToPers (reinterpret_cast<const TRANS*> (trans),
                 reinterpret_cast<PERS*> (pers),
                 log);
  }

◆ transToPersWithKey()

virtual void TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::transToPersWithKey	(	const TRANS *	transObj,
		PERS *	persObj,
		const std::string &	,
		MsgStream &	log )

inlinevirtualinherited

Convert transient representation to persistent one.

Copies data members from transient object to an existing persistent one. Needs to be implemented by the developer on the actual converter.

Parameters

transObj	[IN] transient object
persObj	[IN] persistent object
key	[IN] SG key of object being written.
log	[IN] output message stream

Reimplemented in AthExParticlesCnv_p1, CaloCellContainerCnv_p1, CaloCellLinkContainerCnv_p1, CaloCellLinkContainerCnv_p2, CaloClusterCellLinkContainerCnv_p1, TPConverterWithKeyBase< TRANS, PERS >, and xAODTauJetAuxContainerCnv_v2.

Definition at line 392 of file TPConverter.h.

  {
    return transToPers (transObj, persObj, log);
  }

◆ transToPersWithKeyUntyped()

virtual void TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::transToPersWithKeyUntyped	(	const void *	trans,
		void *	pers,
		const std::string &	key,
		MsgStream &	log )

inlinevirtualinherited

Convert transient object representation to persistent.

Parameters

trans	[IN] void* pointer to the transient object
pers	[OUT] void* pointer to the empty persistent object
key	[IN] SG key of object being written.
log	[IN] output message stream

Reimplemented from ITPCnvBase.

Definition at line 432 of file TPConverter.h.

  {
    transToPersWithKey (reinterpret_cast<const TRANS*> (trans),
                        reinterpret_cast<PERS*> (pers),
                        key,
                        log);
  }

◆ typeID()

template<class TRANS>

virtual const TPObjRef::typeID_t & ITPConverterFor< TRANS >::typeID ( ) const

inlinevirtualinherited

Return TP typeID for persistent objects produced by this converter.

Returns: TPObjRef::typeID_t&

Implements ITPConverter.

Definition at line 208 of file TPConverter.h.

208{ return m_pStorageTID; }

◆ typeIDvalue()

template<class TRANS>

unsigned ITPConverterFor< TRANS >::typeIDvalue ( ) const

inlineinherited

inlined non-virtual version to get the typeID value fast

Definition at line 211 of file TPConverter.h.

211{ return m_pStorageTIDvalue; }

◆ virt_createTransFromPStore()

virtual TRANS * TPPolyCnvBase< TRANS, TRANS, PERS >::virt_createTransFromPStore	(	unsigned	index,
		MsgStream &	log )

inlinevirtualinherited

Internal interface method that is used to invoke the real conversion method (createTransient)

Parameters

index	[IN] index of the persistent object in the storage vector
log	[IN] output message stream

Returns: Created transient object (by pointer)

Reimplemented from TPAbstractPolyCnvBase< TRANS, TRANS, PERS >.

Definition at line 706 of file TPConverter.h.

                                                                               {
     assert (index < this->m_pStorage->size());
     return createTransient( &(*this->m_pStorage)[index], log );
  }

◆ virt_createTransFromPStoreWithKey()

virtual TRANS * TPPolyCnvBase< TRANS, TRANS, PERS >::virt_createTransFromPStoreWithKey	(	unsigned	index,
		const std::string &	key,
		MsgStream &	log )

inlinevirtualinherited

Internal interface method that is used to invoke the real conversion method (createTransient)

Parameters

index	[IN] index of the persistent object in the storage vector
key	[IN] SG key of the object being converted
log	[IN] output message stream

Returns: Created transient object (by pointer)

Reimplemented from TPAbstractPolyCnvBase< TRANS, TRANS, PERS >.

Definition at line 718 of file TPConverter.h.

  {
     assert (index < this->m_pStorage->size());
     return createTransientWithKey( &(*this->m_pStorage)[index], key, log );
  }

◆ virt_toPersistent()

template<class TRANS, class PERS>

virtual TPObjRef TPConverterBase< TRANS, PERS >::virt_toPersistent	(	const TRANS *	trans,
		MsgStream &	log )

inlinevirtualinherited

Internal interface method that is used to invoke the real conversion method (toPersistent_impl) in the derived converter.

Parameters

trans	[IN] transient object
log	[IN] output message stream

Returns: TPObjRef TP reference to the persistent representation stored in the storage vector of the top-level persistent object Here toPersistent_impl is invoked with the dynamic cast of the transient type pointer to it's actual type

Reimplemented from TPAbstractPolyCnvBase< TRANS, TRANS, PERS >.

Definition at line 747 of file TPConverter.h.

                                                                           {
    return this->toPersistentWithKey_impl( trans, "", log);
  }

◆ virt_toPersistentWithKey()

template<class TRANS, class PERS>

virtual TPObjRef TPConverterBase< TRANS, PERS >::virt_toPersistentWithKey	(	const TRANS *	trans,
		const std::string &	key,
		MsgStream &	log )

inlinevirtualinherited

Internal interface method that is used to invoke the real conversion method (toPersistent_impl) in the derived converter.

Parameters

trans	[IN] transient object
key	[IN] SG key of the object being converted.
log	[IN] output message stream

Returns: TPObjRef TP reference to the persistent representation stored in the storage vector of the top-level persistent object Here toPersistentWithKey_impl is invoked with the dynamic cast of the transient type pointer to it's actual type

Reimplemented from TPAbstractPolyCnvBase< TRANS, TRANS, PERS >.

Definition at line 752 of file TPConverter.h.

  {
     return this->toPersistentWithKey_impl( trans, key, log);
  }

◆ wasUsedForReading()

template<class TRANS>

bool ITPConverterFor< TRANS >::wasUsedForReading ( )

inlineinherited

Definition at line 236 of file TPConverter.h.

236{ return m_wasUsedForReading; }

Member Data Documentation

◆ m_curRecLevel

int TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::m_curRecLevel

protectedinherited

count recursive invocations, to detect recursion

Definition at line 582 of file TPConverter.h.

◆ m_ignoreRecursion

bool TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::m_ignoreRecursion

protectedinherited

if true, do not throw errors in case of recursion.

Definition at line 588 of file TPConverter.h.

◆ m_pStorage

std::vector< PERS >* TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::m_pStorage

protectedinherited

the address of the storage vector for persistent representations

Definition at line 579 of file TPConverter.h.

◆ m_pStorageTID

template<class TRANS>

TPObjRef::typeID_t ITPConverterFor< TRANS >::m_pStorageTID

protectedinherited

TP Ref typeID for the persistent objects this converter is creating.

Definition at line 292 of file TPConverter.h.

◆ m_pStorageTIDvalue

template<class TRANS>

unsigned ITPConverterFor< TRANS >::m_pStorageTIDvalue

protectedinherited

m_pStorageTID converted to integer value

Definition at line 295 of file TPConverter.h.

◆ m_recursive

bool TPAbstractPolyCnvBase< TRANS, TRANS, PERS >::m_recursive

protectedinherited

if true, work in recursion-safe way (slower)

Definition at line 585 of file TPConverter.h.

◆ m_topConverter

template<class TRANS>

TopLevelTPCnvBase* ITPConverterFor< TRANS >::m_topConverter

protectedinherited

top level converter that owns this elemental TP converter it also holds the storage object

Definition at line 299 of file TPConverter.h.

◆ m_topConverterRuntime

template<class TRANS>

TopLevelTPCnvBase* ITPConverterFor< TRANS >::m_topConverterRuntime

protectedinherited

top level converter "owning" this TP converter at runtime (different from m_topConverter in case the top-level converter and object have extensions)

Definition at line 302 of file TPConverter.h.

◆ m_wasUsedForReading

template<class TRANS>

bool ITPConverterFor< TRANS >::m_wasUsedForReading

protectedinherited

flag set when using this converter for reading triggers search for a new converter before writing, to prevent possible use of old version

Definition at line 306 of file TPConverter.h.

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

Public Types

Public Member Functions

Protected Attributes

Detailed Description

Member Typedef Documentation

◆ Factory

◆ Pers_t

◆ PersBase_t

◆ PolyCnvBase_t

◆ Trans_t

◆ TransBase_t

Constructor & Destructor Documentation

◆ TRT_HitCollectionCnv_p3()

Member Function Documentation

◆ baseToPersistent()

◆ clearReadingFlag()

◆ converterForRef()

◆ converterForType()

◆ converterNotFound() [1/2]

◆ converterNotFound() [2/2]

◆ createPersistent()

◆ createPersistentWithKey()

◆ createTransFromPStore()

◆ createTransient() [1/2]

◆ createTransient() [2/2]

◆ createTransientWithKey()

◆ fillTransFromPStore()

◆ ignoreRecursion()

◆ initPrivateConverters()

◆ persistentTInfo()

◆ persToTrans() [1/2]

◆ persToTrans() [2/2]

◆ persToTransUntyped()

◆ persToTransWithKey()

◆ persToTransWithKeyUntyped()

◆ pstoreToTrans()

◆ reservePStorage()

◆ setPStorage()

◆ setReadingFlag()

◆ setRecursive()

◆ setRuntimeTopConverter()

◆ setTopConverter()

◆ topConverter() [1/2]

◆ topConverter() [2/2]

◆ toPersistent()

◆ toPersistentWithKey_impl()

◆ transBaseTInfo()

◆ transientTInfo()

◆ transToPers() [1/2]

◆ transToPers() [2/2]

◆ transToPersUntyped()

◆ transToPersWithKey()

◆ transToPersWithKeyUntyped()

◆ typeID()

◆ typeIDvalue()

◆ virt_createTransFromPStore()

◆ virt_createTransFromPStoreWithKey()

◆ virt_toPersistent()

◆ virt_toPersistentWithKey()

◆ wasUsedForReading()

Member Data Documentation

◆ m_curRecLevel

◆ m_ignoreRecursion

◆ m_pStorage

◆ m_pStorageTID

◆ m_pStorageTIDvalue

◆ m_recursive

◆ m_topConverter

◆ m_topConverterRuntime

◆ m_wasUsedForReading