TileL2ContByteStreamCnv Class Reference

This AthConstConverter class provides conversion between ByteStream and TileL2Container. More...

#include <TileL2ContByteStreamCnv.h>

Inheritance diagram for TileL2ContByteStreamCnv:

Collaboration diagram for TileL2ContByteStreamCnv:

Public Types
typedef TileL2ContByteStreamTool	BYTESTREAMTOOL

Public Member Functions
	TileL2ContByteStreamCnv (ISvcLocator *svcloc)
virtual StatusCode	initialize () override
virtual StatusCode	createObjConst (IOpaqueAddress *pAddr, DataObject *&pObj) const override
	Create the transient representation of an object.
virtual StatusCode	createRepConst (DataObject *pObj, IOpaqueAddress *&pAddr) const override
	Convert the transient object to the requested representation.
virtual StatusCode	finalize () override
virtual long	repSvcType () const override
	Storage type and class ID.
virtual StatusCode	createObj (IOpaqueAddress *pAddress, DataObject *&refpObject) override final
	Create the transient representation of an object.
virtual StatusCode	createRep (DataObject *pObject, IOpaqueAddress *&refpAddress) override final
	Convert the transient object to the requested representation.
const std::string &	name () const
void	addTrack (TIDA::Track *t)
const std::vector< TIDA::Track * > &	tracks () const
void	clear ()
void	selectTracks (const Rec::TrackParticleContainer *trigtracks)
void	selectTracks (const Analysis::MuonContainer *muontracks)
bool	msgLvl (const MSG::Level lvl) const
	Test the output level.
MsgStream &	msg () const
	The standard message stream.
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream.
void	setLevel (MSG::Level lvl)
	Change the current logging level.

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

Protected Member Functions
double	phiCorr (double phi)
void	ipCorr (double d0, double z0, double &d0c, double &z0c, double phi0, double eta, double pt)

Protected Attributes
double	m_beamX
double	m_beamY
double	m_beamZ
std::vector< TIDA::Track * >	m_tracks

Private Member Functions
void	initMessaging () const
	Initialize our message level and MessageSvc.

Private Attributes
ToolHandle< BYTESTREAMTOOL >	m_tool
ServiceHandle< IByteStreamCnvSvc >	m_byteStreamCnvSvc
ServiceHandle< StoreGateSvc >	m_storeGate
	Pointer to StoreGateSvc.
ServiceHandle< IROBDataProviderSvc >	m_robSvc
	Pointer to IROBDataProviderSvc.
ToolHandle< TileROD_Decoder >	m_decoder
	Pointer to TileROD_Decoder.
const TileHid2RESrcID *	m_hid2re
	Pointer to TileHid2RESrcID.
Athena::RecyclableDataQueue< TileRecyclableL2Container > m_queue	ATLAS_THREAD_SAFE
	Queue of data objects to recycle.
std::string	m_name
std::string	m_nm
	Message source name.
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels).
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer.
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level.

Detailed Description

This AthConstConverter class provides conversion between ByteStream and TileL2Container.

Author

Aranzazu Ruiz

This class provides methods to convert the bytestream data into TileL2 objects (muon and transverse energy info) and vice versa. The TileL2Container contains a TileL2 object per superdrawer.

Definition at line 66 of file TileL2ContByteStreamCnv.h.

Member Typedef Documentation

BYTESTREAMTOOL

typedef TileL2ContByteStreamTool TileL2ContByteStreamCnv::BYTESTREAMTOOL

Definition at line 72 of file TileL2ContByteStreamCnv.h.

Constructor & Destructor Documentation

TileL2ContByteStreamCnv()

TileL2ContByteStreamCnv::TileL2ContByteStreamCnv

(

ISvcLocator *

svcloc

)

Definition at line 60 of file TileL2ContByteStreamCnv.cxx.

  : AthConstConverter(storageType(), classID(), svcloc, "TileL2ContByteStreamCnv")
  , m_tool("TileL2ContByteStreamTool")
  , m_byteStreamCnvSvc("ByteStreamCnvSvc", name())
  , m_storeGate("StoreGateSvc", name())
  , m_robSvc("ROBDataProviderSvc", name())
  , m_decoder("TileROD_Decoder")
  , m_hid2re(0)
{
}

Member Function Documentation

addTrack()

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

inlineinherited

Definition at line 43 of file Converter.h.

                              {
    m_tracks.push_back(t);
  }

classID()

CLID TileL2ContByteStreamCnv::classID ( )

static

Definition at line 71 of file TileL2ContByteStreamCnv.cxx.

{ return ClassID_traits<TileL2Container>::ID(); }

clear()

void Converter::clear ( )

inlineinherited

Definition at line 51 of file Converter.h.

{m_tracks.clear();} 

createObj()

StatusCode AthConstConverter::createObj ( IOpaqueAddress * pAddress, DataObject *& refpObject )

finaloverridevirtualinherited

Create the transient representation of an object.

Non-const version; just calls the const version.

Definition at line 44 of file AthConstConverter.cxx.

{
  return createObjConst (pAddress, refpObject);
}

createObjConst()

StatusCode TileL2ContByteStreamCnv::createObjConst ( IOpaqueAddress * pAddress, DataObject *& refpObject ) const

overridevirtual

Create the transient representation of an object.

The default implementation is a no-op.

Reimplemented from AthConstConverter.

Definition at line 98 of file TileL2ContByteStreamCnv.cxx.

{
  ATH_MSG_DEBUG( " Executing createObj method" );
 
  ByteStreamAddress* pRE_Addr{};
  ATH_CHECK( (pRE_Addr = dynamic_cast<ByteStreamAddress*>(pAddr)) != nullptr );
 
  const EventContext& ctx = pRE_Addr->getEventContext();
  const RawEvent* re = m_robSvc->getEvent(ctx);
  if (!re) {
    ATH_MSG_ERROR( "Could not get raw event from ByteStreamInputSvc" );
    return StatusCode::FAILURE;
  }
 
  TileL2Container* cont = m_queue.get (*m_decoder);
 
  if (!m_decoder->convert(re, cont).isSuccess()) {
    ATH_MSG_WARNING( "Conversion tool returned an error. TileL2Container might be empty." );
  }
 
  pObj = SG::asStorable( cont ) ; 
 
  return StatusCode::SUCCESS;  
}

createRep()

StatusCode AthConstConverter::createRep ( DataObject * pObject, IOpaqueAddress *& refpAddress )

finaloverridevirtualinherited

Convert the transient object to the requested representation.

Non-const version; just calls the const version.

Definition at line 56 of file AthConstConverter.cxx.

{
  return createRepConst (pObject, refpAddress);
}

createRepConst()

StatusCode TileL2ContByteStreamCnv::createRepConst ( DataObject * pObject, IOpaqueAddress *& refpAddress ) const

overridevirtual

Convert the transient object to the requested representation.

The default implementation is a no-op.

Reimplemented from AthConstConverter.

Definition at line 123 of file TileL2ContByteStreamCnv.cxx.

{
  // convert TileL2s in the container into ByteStream
 
  ATH_MSG_DEBUG( " Executing createRep method" );
 
  // get Full Event Assembler
  FullEventAssembler<TileHid2RESrcID>* fea = 0;
  std::string key("Tile");
  
  ATH_CHECK( m_byteStreamCnvSvc->getFullEventAssembler(fea, key) );
 
  // create TileL2Container
  TileL2Container* l2cont(0);
  SG::fromStorable(pObj, l2cont);
  if(!l2cont){
    ATH_MSG_ERROR( " Can not cast to TileL2Container " );
    return StatusCode::FAILURE;
  }
 
  std::string name = pObj->registry()->name();
  ByteStreamAddress* addr = new ByteStreamAddress(classID(), name, "");
 
  pAddr = addr;
 
  // call TileL2ContByteStreamTool
  ATH_CHECK( m_tool->convert(l2cont, fea) );
 
  return StatusCode::SUCCESS;
}

finalize()

StatusCode TileL2ContByteStreamCnv::finalize ( )

overridevirtual

Definition at line 155 of file TileL2ContByteStreamCnv.cxx.

{
  return Converter::finalize(); 
}

initialize()

StatusCode TileL2ContByteStreamCnv::initialize ( )

overridevirtual

Definition at line 75 of file TileL2ContByteStreamCnv.cxx.

                                               {
 
  ATH_CHECK(Converter::initialize());
 
  ATH_MSG_DEBUG(" initialize ");
 
  // Get ByteStreamCnvSvc
  ATH_CHECK( m_byteStreamCnvSvc.retrieve() );
 
  // retrieve Tool
  ATH_CHECK( m_decoder.retrieve() );
  m_hid2re = m_decoder->getHid2re();
 
  ATH_CHECK( m_tool.retrieve() );
 
  ATH_CHECK( m_robSvc.retrieve() );
 
  ATH_CHECK( m_storeGate.retrieve() );
 
  return StatusCode::SUCCESS ;
}

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();
  // If user did not set an explicit level, set a default
  if (m_lvl == MSG::NIL) {
    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 200 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;
      
    }
  }

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 167 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 182 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 user did not set explicit message level we have to initialize
  // the messaging and retrieve the default via the MessageSvc.
  if (m_lvl==MSG::NIL && !m_initialized.test_and_set()) initMessaging();
 
  if (m_lvl <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

name()

const std::string & AthConstConverter::name ( ) const

inlineinherited

Definition at line 67 of file AthConstConverter.h.

{ return m_name; }

phiCorr()

double Converter::phiCorr ( double phi )

inlineprotectedinherited

Definition at line 191 of file Converter.h.

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

repSvcType()

virtual long TileL2ContByteStreamCnv::repSvcType ( ) const

inlineoverridevirtual

Storage type and class ID.

Definition at line 80 of file TileL2ContByteStreamCnv.h.

{ return i_repSvcType(); }

selectTracks() [1/2]

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

inlineinherited

Definition at line 120 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/2]

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

inlineinherited

Definition at line 54 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++;
    }
    
  }

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;
}

storageType()

long TileL2ContByteStreamCnv::storageType ( )

static

Definition at line 73 of file TileL2ContByteStreamCnv.cxx.

{ return ByteStreamAddress::storageType(); }

tracks()

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

inlineinherited

Definition at line 48 of file Converter.h.

{return m_tracks;}

Member Data Documentation

ATLAS_THREAD_SAFE

Athena::RecyclableDataQueue<TileRecyclableL2Container> m_queue TileL2ContByteStreamCnv::ATLAS_THREAD_SAFE

mutableprivate

Queue of data objects to recycle.

Definition at line 103 of file TileL2ContByteStreamCnv.h.

m_beamX

double Converter::m_beamX

protectedinherited

Definition at line 263 of file Converter.h.

m_beamY

double Converter::m_beamY

protectedinherited

Definition at line 263 of file Converter.h.

m_beamZ

double Converter::m_beamZ

protectedinherited

Definition at line 263 of file Converter.h.

m_byteStreamCnvSvc

ServiceHandle<IByteStreamCnvSvc> TileL2ContByteStreamCnv::m_byteStreamCnvSvc

private

Definition at line 88 of file TileL2ContByteStreamCnv.h.

m_decoder

ToolHandle<TileROD_Decoder> TileL2ContByteStreamCnv::m_decoder

private

Pointer to TileROD_Decoder.

Definition at line 97 of file TileL2ContByteStreamCnv.h.

m_hid2re

const TileHid2RESrcID* TileL2ContByteStreamCnv::m_hid2re

private

Pointer to TileHid2RESrcID.

Definition at line 100 of file TileL2ContByteStreamCnv.h.

m_imsg

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

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

{ nullptr };

m_lvl

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

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

{ MSG::NIL };

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_name

std::string AthConstConverter::m_name

privateinherited

Definition at line 71 of file AthConstConverter.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

m_robSvc

ServiceHandle<IROBDataProviderSvc> TileL2ContByteStreamCnv::m_robSvc

private

Pointer to IROBDataProviderSvc.

Definition at line 94 of file TileL2ContByteStreamCnv.h.

m_storeGate

ServiceHandle<StoreGateSvc> TileL2ContByteStreamCnv::m_storeGate

private

Pointer to StoreGateSvc.

Definition at line 91 of file TileL2ContByteStreamCnv.h.

m_tool

ToolHandle<BYTESTREAMTOOL> TileL2ContByteStreamCnv::m_tool

private

Definition at line 86 of file TileL2ContByteStreamCnv.h.

m_tracks

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

protectedinherited

Definition at line 264 of file Converter.h.

---

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

• TileL2ContByteStreamCnv.h
• TileL2ContByteStreamCnv.cxx