ByteStream converter for the MuCTPI_RDO object. More...

#include <MuCTPIByteStreamCnv.h>

Inheritance diagram for MuCTPIByteStreamCnv:

Collaboration diagram for MuCTPIByteStreamCnv:

Public Member Functions
	MuCTPIByteStreamCnv (ISvcLocator *svcloc)
	Standard constructor.
virtual StatusCode	initialize () override
	Function connecting to all the needed services/tools.
virtual StatusCode	createObj (IOpaqueAddress pAddr, DataObject &pObj) override
	Function creating the MuCTPI_RDO object from a MuCTPI ROB fragment.
virtual StatusCode	createRep (DataObject pObj, IOpaqueAddress &pAddr) override
	Function creating the MuCTPI ROB fragment from a MuCTPI_RDO object.
virtual long	repSvcType () const override
	Function needed by the framework.
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)

Static Public Member Functions
static long	storageType ()
	Function needed by the framework.
static CLID	classID ()
	Function needed by the framework.

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 Attributes
ToolHandle< MuCTPIByteStreamTool >	m_tool
	Tool doing the actual conversion.
MuCTPISrcIdMap	m_srcIdMap
	Object storing the various IDs of the MuCTPI fragment.
ServiceHandle< IROBDataProviderSvc >	m_robDataProvider
	Service used when reading the BS data.
ServiceHandle< IByteStreamEventAccess >	m_ByteStreamEventAccess
	Service used when writing the BS data.

Detailed Description

ByteStream converter for the MuCTPI_RDO object.

   This is a quite standard BS converter for the output sent from the
   "Muon to Central Trigger Processor Interface" to the DAQ. It uses
   an external tool (MuCTPIByteStreamTool) to do the actual work of
   the conversion, this converter is "only" supposed to communicate
   with the framework.

See also: MuCTPIByteStreamTool

Author: Tadashi Maeno; Attila Krasznahorkay

Definition at line 47 of file MuCTPIByteStreamCnv.h.

Constructor & Destructor Documentation

◆ MuCTPIByteStreamCnv()

MuCTPIByteStreamCnv::MuCTPIByteStreamCnv ( ISvcLocator * svcloc )

Standard constructor.

The constructor sets up all the ToolHandle and ServiceHandle objects and initialises the base class in the correct way.

Definition at line 33 of file MuCTPIByteStreamCnv.cxx.

  : Converter( storageType(), classID(), svcloc ),
    m_tool( "MuCTPIByteStreamTool" ),
#ifdef CTP_MUCTPI_HAVE_SAME_ROS
    m_ctp_tool( "CTPByteStreamTool" ),
    m_ctp_rec_tool( "RecCTPByteStreamTool" ),
    m_storeGate( "StoreGateSvc", "MuCTPIByteStreamCnv" ),
#endif
    m_robDataProvider( "ROBDataProviderSvc", "MuCTPIByteStreamCnv" ),
    m_ByteStreamEventAccess( "ByteStreamCnvSvc", "MuCTPIByteStreamCnv" ) {
 
}

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 MuCTPIByteStreamCnv::classID ( )

static

Function needed by the framework.

Function telling the framework the Class ID of the object that this converter is for (MuCTPI_RDO).

Definition at line 50 of file MuCTPIByteStreamCnv.cxx.

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

◆ clear()

void Converter::clear ( )

inlineinherited

Definition at line 51 of file Converter.h.

51{m_tracks.clear();}

◆ createObj()

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

overridevirtual

Function creating the MuCTPI_RDO object from a MuCTPI ROB fragment.

This function creates the MuCTPI_RDO object from the BS data.

It requests the ROB fragment with the ROB Id of the MuCTPI and gives this fragment to the MuCTPIByteStreamTool for conversion.

Definition at line 116 of file MuCTPIByteStreamCnv.cxx.

                                                                                    {
 
  MsgStream log( msgSvc(), "MuCTPIByteStreamCnv" );
  log << MSG::DEBUG << "createObj() called" << endmsg;
 
  ByteStreamAddress* pBS_Addr{};
  ATH_CHECK( (pBS_Addr = dynamic_cast<ByteStreamAddress*>(pAddr)) != nullptr );
 
  const EventContext& ctx = pBS_Addr->getEventContext();
 
  log << MSG::DEBUG << " Creating Objects " << *( pBS_Addr->par() ) << endmsg;
 
  //
  // Get SourceID:
  //
  const uint32_t robId = m_srcIdMap.getRobID( m_srcIdMap.getRodID() );
 
  log << MSG::DEBUG << " expected ROB sub-detector ID: " << std::hex 
      << robId << std::dec << endmsg;  
 
  std::vector< uint32_t > vID;
  vID.push_back( robId );
 
  //
  // Get ROB fragment:
  //
  IROBDataProviderSvc::VROBFRAG robFrags;
  m_robDataProvider->getROBData(ctx, vID, robFrags );
 
  //
  // Size check:
  //
  if( robFrags.size() == 0 ) {
    
    // This is a hack - for some time CTP and MUCTPI were plugged in
    // the same ROB of the same ROS PC ...
    // 0x76000b : MUCTPI
    // 0x76000a : CTP
    // Note that the MUCTPI ROBID must come first, cos further down we
    // just use the first element of robFrags to create the RDO
    vID[0] = ( robId | 0x00000b );
#ifdef CTP_MUCTPI_HAVE_SAME_ROS
    vID.push_back( robId | 0x00000a );
#endif
    m_robDataProvider->getROBData(ctx, vID, robFrags );
 
    if( robFrags.size() == 0 ) {
      // yet another temporary hack!!! - for some time the ROB ID was 0!
      vID[0] = 0x00000b;
#ifdef CTP_MUCTPI_HAVE_SAME_ROS
      vID[1] = 0x00000a;
#endif
      m_robDataProvider->getROBData(ctx,  vID, robFrags );
      
      if( robFrags.size() == 0 ) {
        log << MSG::WARNING << "No MuCTPI ROB fragments found!" << endmsg;
 
        MuCTPI_RDO* result = new MuCTPI_RDO;
        pObj = SG::asStorable( result );
        return StatusCode::SUCCESS;
      }
    }
  }
 
  IROBDataProviderSvc::VROBFRAG::const_iterator it = robFrags.begin();
  MuCTPI_RDO* result = 0;
  StatusCode sc = m_tool->convert( ROBData( *it ).getROBFragment(), result );
  if( sc.isFailure() ) {
    log << MSG::ERROR << " Failed to create Objects: " << *( pBS_Addr->par() ) << endmsg;
    if (result!=0){
      delete result;
      result=0;
    }
    return sc;
  }
  pObj = SG::asStorable( result ) ;
 
#ifdef CTP_MUCTPI_HAVE_SAME_ROS
  // If CTP and MuCTPI send to the same ROS, the MuCTPI ROS, we get
  // the CTP fragment in the MuCTPI ROS fragment. In order to be able
  // to produce a CTP_RDO, we call the CTP ByteStream AlgTool from
  // here whenever we think we have a CTP fragment.
 
  if( robFrags.size() > 1 ) {
 
    ++it;
    if( it != robFrags.end() ) {
    
      CTP_RDO* ctp_result = 0;
      CTP_RIO* ctp_rio_result = 0;
 
      ATH_CHECK(  m_ctp_tool->convert( ROBData( *it ).getROBFragment(), ctp_result ) );
 
      SG::asStorable( ctp_result ) ;
      ATH_CHECK(  m_storeGate->record( ctp_result, "CTP_RDO" ) );
 
      ATH_CHECK(  m_ctp_rec_tool->convert( ROBData( *it ).getROBFragment(), ctp_rio_result, log ) );
 
      SG::asStorable( ctp_rio_result ) ;
      ATH_CHECK(  m_storeGate->record( ctp_rio_result, "CTP_RIO" ) );
    }
 
  }
#endif
 
  return sc;
}

◆ createRep()

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

overridevirtual

Function creating the MuCTPI ROB fragment from a MuCTPI_RDO object.

This function receives a MuCTPI_RDO object as input, and adds the ROB fragment of the MuCTPI to the current raw event using the IByteStreamEventAccess interface.

Definition at line 229 of file MuCTPIByteStreamCnv.cxx.

                                                                                    {
 
  MsgStream log( msgSvc(), "MuCTPIByteStreamCnv" );
  log << MSG::DEBUG << "createRep() called" << endmsg;
 
  RawEventWrite* re = m_ByteStreamEventAccess->getRawEvent();
 
  MuCTPI_RDO* result;
  if( ! SG::fromStorable( pObj, result ) ) {
    log << MSG::ERROR << " Cannot cast to MuCTPI_RDO" << endmsg;
    return StatusCode::FAILURE;
  }
 
  ByteStreamAddress* addr = new ByteStreamAddress( classID(), pObj->registry()->name(), "" );
 
  pAddr = addr;
 
  // Convert to ByteStream
  return m_tool->convert( result, re );
}

◆ initialize()

StatusCode MuCTPIByteStreamCnv::initialize ( )

overridevirtual

Function connecting to all the needed services/tools.

Init method gets all necessary services etc.

Definition at line 63 of file MuCTPIByteStreamCnv.cxx.

                                           {
 
  //
  // Initialise the base class:
  //
  ATH_CHECK(  Converter::initialize() );
 
  MsgStream log( msgSvc(), "MuCTPIByteStreamCnv" );
  log << MSG::DEBUG << "MuCTPIByteStreamCnv in initialize() " << endmsg;
 
  //
  // Get ByteStreamCnvSvc:
  //
  ATH_CHECK(  m_ByteStreamEventAccess.retrieve() );
 
  //
  // Get MuCTPIByteStreamTool:
  //
  ATH_CHECK(  m_tool.retrieve() );
  log << MSG::DEBUG << "Connected to MuCTPIByteStreamTool" << endmsg;
 
#ifdef CTP_MUCTPI_HAVE_SAME_ROS
  //
  // Get CTPByteStreamTool:
  //
  ATH_CHECK( m_ctp_tool.retrieve() );
  log << MSG::DEBUG << "Connected to CTPByteStreamTool" << endmsg;
 
  //
  // Get RecCTPByteStreamTool:
  //
  ATH_CHECK(  m_ctp_rec_tool.retrieve() );
  log << MSG::DEBUG << "Connected to RecCTPByteStreamTool" << endmsg;
#endif
 
  //
  // Get ROBDataProvider:
  //
  if( m_robDataProvider.retrieve().isFailure() ) {
    log << MSG::WARNING << "Can't get ROBDataProviderSvc" << endmsg;
    // return is disabled for Write BS which does not requre ROBDataProviderSvc
  } else {
    log << MSG::DEBUG << "Connected to ROBDataProviderSvc" << endmsg;
  }
 
  return StatusCode::SUCCESS;
}

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

◆ 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 MuCTPIByteStreamCnv::repSvcType ( ) const

inlineoverridevirtual

Function needed by the framework.

Definition at line 61 of file MuCTPIByteStreamCnv.h.

61{ 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++;
    }
    
  }

◆ storageType()

long MuCTPIByteStreamCnv::storageType ( )

static

Function needed by the framework.

Definition at line 56 of file MuCTPIByteStreamCnv.cxx.

                                      {
  return ByteStreamAddress::storageType();
}

◆ tracks()

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

inlineinherited

Definition at line 48 of file Converter.h.

48{return m_tracks;}

Member Data Documentation

◆ 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_ByteStreamEventAccess

ServiceHandle< IByteStreamEventAccess > MuCTPIByteStreamCnv::m_ByteStreamEventAccess

private

Service used when writing the BS data.

Definition at line 82 of file MuCTPIByteStreamCnv.h.

◆ m_robDataProvider

ServiceHandle< IROBDataProviderSvc > MuCTPIByteStreamCnv::m_robDataProvider

private

Service used when reading the BS data.

Definition at line 80 of file MuCTPIByteStreamCnv.h.

◆ m_srcIdMap

MuCTPISrcIdMap MuCTPIByteStreamCnv::m_srcIdMap

private

Object storing the various IDs of the MuCTPI fragment.

Definition at line 71 of file MuCTPIByteStreamCnv.h.

◆ m_tool

ToolHandle< MuCTPIByteStreamTool > MuCTPIByteStreamCnv::m_tool

private

Tool doing the actual conversion.

Definition at line 69 of file MuCTPIByteStreamCnv.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:

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MuCTPIByteStreamCnv()

Member Function Documentation

◆ addTrack()

◆ classID()

◆ clear()

◆ createObj()

◆ createRep()

◆ initialize()

◆ ipCorr()

◆ phiCorr()

◆ repSvcType()

◆ selectTracks() [1/2]

◆ selectTracks() [2/2]

◆ storageType()

◆ tracks()

Member Data Documentation

◆ m_beamX

◆ m_beamY

◆ m_beamZ

◆ m_ByteStreamEventAccess

◆ m_robDataProvider

◆ m_srcIdMap

◆ m_tool

◆ m_tracks