HLT::HLTResultByteStreamCnv Class Reference

Gaudi bytestream Converter for the HLTResult class. More...

#include <HLTResultByteStreamCnv.h>

Inheritance diagram for HLT::HLTResultByteStreamCnv:

Collaboration diagram for HLT::HLTResultByteStreamCnv:

Public Member Functions
	HLTResultByteStreamCnv (ISvcLocator *svcloc)
	std Gaudi converter constructor
virtual StatusCode	initialize ()
	Gaudi initialize.
virtual StatusCode	finalize ()
	Gaudi finalize.
virtual StatusCode	createObj (IOpaqueAddress *pAddr, DataObject *&pObj)
	create RDO (HLTResult) from bytestream
virtual StatusCode	createRep (DataObject *pObj, IOpaqueAddress *&pAddr)
	create bytestream from RDOs (HLTResult)
long	repSvcType () const
	return repSvcType
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 ()
	storageType
static CLID	classID ()
	CLID.

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
MsgStream *	m_log
	std message stream
ToolHandle< HLTResultByteStreamTool >	m_tool
	tool for conversion work
ServiceHandle< IByteStreamEventAccess >	m_ByteStreamEventAccess
	Gaudi helper.
ServiceHandle< IROBDataProviderSvc >	m_robDataProvider
	Gaudi ROB helper.

Detailed Description

Gaudi bytestream Converter for the HLTResult class.

An AlgTool to provide conversion of Run-1&2 HLTResult from/to ByteStream.

This class essentially provides two methods:

1. createObj(IOpaqueAddress*, DataObject*&)
2. createRep(DataObject*, IOpaqueAddress*& )

The first one creates a DataObject from persistent information, captured in IOpaqueAddress, and the second is supposed to write the DataObject to a persistent representation. Here, DataObject will be casted to HLTResult

see here for a general ByteStream converter description: http://atlas-sw.cern.ch/cgi-bin/viewcvs-atlas.cgi/offline/Event/ByteStreamCnvSvc/doc/ByteStream.txt

beginning of this file:
ByteStream conversion Service (ByteStreamCnvSvc) provides a means to convert the raw data in the ByteStream format to and from the transient representation, Raw Data Objects (RDO). In Gaudi's terminology, ByteStream can be considered as a persistency for reading and storing raw data. In practice, the source of the ByteStream can be a file, or the online data source (or its emulation) through network, for example. Together with the IdentifiableContainer, used to contain the Raw Data Objects, ByteStreamCnvSvc provides a mechanism for creating the RDOs on demand.

This tool is used by the Gaudi converter class: HLTResultByteStreamCnv.

Definition at line 69 of file HLTResultByteStreamCnv.h.

Constructor & Destructor Documentation

HLTResultByteStreamCnv()

HLT::HLTResultByteStreamCnv::HLTResultByteStreamCnv

(

ISvcLocator *

svcloc

)

std Gaudi converter constructor

Constructor.

Definition at line 34 of file HLTResultByteStreamCnv.cxx.

                                                                     :
  Converter(storageType(), classID(),svcloc),
  m_log(0),
  m_tool("HLT::HLTResultByteStreamTool"),
  m_ByteStreamEventAccess("ByteStreamCnvSvc", "HLTResultByteStreamCnv"),
  m_robDataProvider("ROBDataProviderSvc", "HLTResultByteStreamCnv")
{ }

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

static

CLID.

Definition at line 45 of file HLTResultByteStreamCnv.cxx.

{
  return ClassID_traits<HLT::HLTResult>::ID() ;
}

clear()

void Converter::clear ( )

inlineinherited

Definition at line 51 of file Converter.h.

{m_tracks.clear();} 

createObj()

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

virtual

create RDO (HLTResult) from bytestream

createObj should create the RDO from bytestream.

Definition at line 148 of file HLTResultByteStreamCnv.cxx.

{
  HLTResult* result = new HLTResult();
 
  ByteStreamAddress *pBS_Addr = dynamic_cast<ByteStreamAddress*>(pAddr);
  if(pBS_Addr==nullptr) {
    (*m_log) << MSG::ERROR << " Can not cast to ByteStreamAddress " << endmsg ;
    delete result;
    return StatusCode::FAILURE;
  }
 
  const std::string nm = *(pBS_Addr->par()) ;
  std::vector<ROBData> vRobData;
 
  StatusCode sc = m_tool->convert(*m_robDataProvider, result, nm);
  if (sc != StatusCode::SUCCESS) {
    (*m_log) << MSG::ERROR << "Failed to create object " << nm << endmsg;
    return sc;
  }
 
  pObj = SG::asStorable(result);
  return sc;
}

createRep()

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

virtual

create bytestream from RDOs (HLTResult)

createRep should create the bytestream from RDOs.

Definition at line 175 of file HLTResultByteStreamCnv.cxx.

{
 
   HLTResult* result = 0;
   SG::fromStorable(pObj, result);
 
   if (!result) {
     (*m_log) << MSG::ERROR << " Cannot cast to HLTResult " << endmsg ;
     return StatusCode::FAILURE;
   }
 
   std::string nm = pObj->registry()->name();
 
   RawEventWrite* re = m_ByteStreamEventAccess->getRawEvent();
  // Convert to ByteStream
   StatusCode sc = m_tool->convert(result, re, nm);
 
   if ( pAddr != nullptr ) pAddr->release();
   ByteStreamAddress* addr = new ByteStreamAddress(classID(), nm, "");
   pAddr = addr;
   pAddr->addRef();
 
   return sc;
 
}

finalize()

StatusCode HLT::HLTResultByteStreamCnv::finalize ( )

virtual

Gaudi finalize.

Definition at line 55 of file HLTResultByteStreamCnv.cxx.

{
  delete m_log;
  return StatusCode::SUCCESS;
}

initialize()

StatusCode HLT::HLTResultByteStreamCnv::initialize ( )

virtual

Gaudi initialize.

Init method gets all necessary services etc.

Definition at line 65 of file HLTResultByteStreamCnv.cxx.

{
   StatusCode sc = Converter::initialize();
   if(StatusCode::SUCCESS!=sc)
   {
    return sc;
   }
 
   m_log = new MsgStream(msgSvc(), "HLTResultByteStreamCnv");
 
   (*m_log) << MSG::DEBUG << "HLTResultByteStreamCnv in initialize() " <<endmsg;
 
   //Get ByteStreamCnvSvc
   if ( m_ByteStreamEventAccess.retrieve().isFailure() ) {
     (*m_log) << MSG::FATAL << "failed to retrieve ByteStreamEventAccess service: "
          << m_ByteStreamEventAccess << endmsg;
     return  StatusCode::FAILURE;
   } else {
     (*m_log) << MSG::DEBUG << "successfully retrieved ByteStreamEventAccess service: "
          << m_ByteStreamEventAccess << endmsg;
   }
//    if(StatusCode::SUCCESS != service("ByteStreamCnvSvc",m_ByteStreamEventAccess)){
//      (*m_log) << MSG::ERROR << " Can't get ByteStreamEventAccess interface " << endmsg;
//      return StatusCode::FAILURE;
//    }
 
 
 
   // Retrieve Tool
   //   IToolSvc* toolSvc;
   //   if(StatusCode::SUCCESS != service("ToolSvc",toolSvc)){
   //     (*m_log) << MSG::ERROR << " Can't get ToolSvc " << endmsg;
   //     return StatusCode::FAILURE;
   //   }
 
//    std::string toolType = "HLT::HLTResultByteStreamTool" ;
//    if(StatusCode::SUCCESS !=toolSvc->retrieveTool(toolType,m_tool))
//    {
//      (*m_log) << MSG::ERROR << " Can't get ByteStreamTool " << endmsg;
//      return StatusCode::FAILURE;
//    }
   // get the converter tool:
   if ( m_tool.retrieve().isFailure() ) {
     (*m_log) << MSG::FATAL << "failed to retrieve converter tool : "
          << m_tool << endmsg;
     return  StatusCode::FAILURE;
   } else {
     (*m_log) << MSG::DEBUG << "successfully retrieved converter tool: "
          << m_tool << endmsg;
   }
 
 
   // Get ROBDataProvider
   if (m_robDataProvider.retrieve().isFailure()) {
    (*m_log) << MSG::FATAL << "failed to retrieve ROBDataProviderSvc service: "
         << m_robDataProvider << endmsg;
    return StatusCode::FAILURE;
   }
   (*m_log) << MSG::DEBUG << "successfully retrieved ROBDataProviderSvc service: "
        << m_robDataProvider << endmsg;
 
//    IService* robSvc ;
//    sc = serviceLocator()->getService("ROBDataProviderSvc", robSvc);
//    if(sc != StatusCode::SUCCESS ) {
//      (*m_log)<<MSG::ERROR << " Cant get ROBDataProviderSvc " <<endmsg;
 
//      //DBG:
//      //     return sc ;
//      return StatusCode::SUCCESS;
//    }
 
//    m_robDataProvider = dynamic_cast<IROBDataProviderSvc*> (robSvc);
//    if(m_robDataProvider == 0 ) {
//      (*m_log)<<MSG::ERROR<< " Cant cast to ROBDataProviderSvc " <<endmsg;
//      // return StatusCode::FAILURE ;
//    }
 
   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()

long HLT::HLTResultByteStreamCnv::repSvcType ( ) const

inline

return repSvcType

Definition at line 85 of file HLTResultByteStreamCnv.h.

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 HLT::HLTResultByteStreamCnv::storageType ( )

static

storageType

Definition at line 50 of file HLTResultByteStreamCnv.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

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 > HLT::HLTResultByteStreamCnv::m_ByteStreamEventAccess

private

Gaudi helper.

Definition at line 93 of file HLTResultByteStreamCnv.h.

m_log

MsgStream* HLT::HLTResultByteStreamCnv::m_log

private

std message stream

Definition at line 89 of file HLTResultByteStreamCnv.h.

m_robDataProvider

ServiceHandle< IROBDataProviderSvc > HLT::HLTResultByteStreamCnv::m_robDataProvider

private

Gaudi ROB helper.

Definition at line 95 of file HLTResultByteStreamCnv.h.

m_tool

ToolHandle< HLTResultByteStreamTool > HLT::HLTResultByteStreamCnv::m_tool

private

tool for conversion work

Definition at line 91 of file HLTResultByteStreamCnv.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:

• HLTResultByteStreamCnv.h
• HLTResultByteStreamCnv.cxx