TrigSerializeConverter< DATA > Class Template Reference

#include <TrigSerializeConverter.h>

Inheritance diagram for TrigSerializeConverter< DATA >:

Collaboration diagram for TrigSerializeConverter< DATA >:

Public Member Functions
StatusCode	createObj (IOpaqueAddress *iAddr, DataObject *&pO)
StatusCode	createRep (DataObject *pObj, IOpaqueAddress *&pAddr)
StatusCode	initialize ()
virtual long int	repSvcType () const
	TrigSerializeConverter (ISvcLocator *svcloc)
	~TrigSerializeConverter ()
void	addTrack (TIDA::Track *t)
const std::vector< TIDA::Track * > &	tracks () const
void	clear ()
void	selectTracks (const TrigInDetTrackCollection *trigtracks)
void	selectTracks (const Rec::TrackParticleContainer *trigtracks)
void	selectTracks (const Analysis::MuonContainer *muontracks)

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

Static Public Attributes
static const bool	typeIsxAOD = isxAOD<DATA>::value

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< ITrigSerializeConvHelper >	m_convHelper
ServiceHandle< StoreGateSvc >	m_sgHandle
MsgStream *	m_log
StoreGateSvc *	m_sgsvc

Detailed Description

template<typename DATA>
class TrigSerializeConverter< DATA >

Definition at line 139 of file TrigSerializeConverter.h.

Constructor & Destructor Documentation

TrigSerializeConverter()

template<typename DATA>

TrigSerializeConverter< DATA >::TrigSerializeConverter ( ISvcLocator * svcloc )

inline

Definition at line 280 of file TrigSerializeConverter.h.

      : Converter( storageType(), classID(), svcloc ),
        m_convHelper( "TrigSerializeConvHelper/TrigSerializeConvHelper" ),
        m_sgHandle( "StoreGateSvc", "TrigSerializeConverter" ),
        m_log( 0 ), m_sgsvc( 0 ) {
 
      m_log = new MsgStream( msgSvc(), "TrigSerializeConverter" );
   }

~TrigSerializeConverter()

template<typename DATA>

TrigSerializeConverter< DATA >::~TrigSerializeConverter ( )

inline

Definition at line 289 of file TrigSerializeConverter.h.

                            {
      delete m_log;
   }

Member Function Documentation

addTrack()

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

inlineinherited

Definition at line 45 of file Converter.h.

                              {
    m_tracks.push_back(t);
  }

classID()

template<typename DATA>

const CLID & TrigSerializeConverter< DATA >::classID ( )

inlinestatic

Definition at line 275 of file TrigSerializeConverter.h.

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

clear()

void Converter::clear ( )

inlineinherited

Definition at line 53 of file Converter.h.

{m_tracks.clear();} 

createObj()

template<typename DATA>

StatusCode TrigSerializeConverter< DATA >::createObj ( IOpaqueAddress * iAddr, DataObject *& pO )

inline

Definition at line 145 of file TrigSerializeConverter.h.

                                                                   {
 
      const std::string clname = ClassID_traits< DATA >::typeName();
 
      std::string normalized = SG::normalizedTypeinfoName( typeid(DATA) ); // normalized string representation (matching the names used in the root dictionary)
 
 
      void *ptr( 0 );
      if( m_log->level() <= MSG::DEBUG ) {
         *m_log << MSG::DEBUG << "In createObj for : " << clname << " normalized to " << normalized << endmsg;
      }
 
      StatusCode sc = m_convHelper->createObj( normalized, iAddr, ptr, typeIsxAOD);
      if( sc.isFailure() ){
         *m_log << MSG::WARNING << "m_convHelper->createObj failed for "
                << clname << endmsg;
      }
 
      DATA *nObj = ( DATA* ) ptr;
      nObj = TrigSerialization::finishRead (nObj);
      if (!nObj){
        *m_log << MSG::ERROR << "TrigSerializeConverter::createObj: nObj is null."<<  endmsg;
        return StatusCode::FAILURE;
      }
      pO = SG::asStorable( nObj );
      if( m_log->level() <= MSG::DEBUG ) {
         *m_log << MSG::DEBUG << "IOpaq: " << iAddr 
                << " created nObj: " << nObj << endmsg;
      }
 
      
      
 
      if( m_sgsvc){
         TrigStreamAddress *addr = dynamic_cast< TrigStreamAddress* >( iAddr );
         if( addr ) {
            if( m_sgsvc->contains< DATA >( addr->sgkey() ) ) {
               *m_log << MSG::WARNING << "TrigSerializeConverter::createObj object "
                      << clname << " / " << addr->sgkey()
                      << " is already in the store; not overwriting"
                      <<  endmsg;
            } else {
               sc = m_sgsvc->record( nObj, addr->sgkey() );
            }
            if( sc.isFailure() ) {
               *m_log << MSG::ERROR << "SG::record failed for " << addr->sgkey()
                      <<  endmsg;
            } else if( m_log->level() <= MSG::DEBUG ) {
               *m_log << MSG::DEBUG << "SG::record key: "
                      << addr->sgkey() << " class: " << clname << endmsg;
            }
         } else {
            *m_log << MSG::WARNING << "createObj cast failed" << endmsg;
         }
      } else if( m_log->level() <= MSG::DEBUG ) {
         *m_log << MSG::DEBUG << "did not put an object into SG" << endmsg;
      }
 
      return sc;
   }

createRep()

template<typename DATA>

StatusCode TrigSerializeConverter< DATA >::createRep ( DataObject * pObj, IOpaqueAddress *& pAddr )

inline

Definition at line 207 of file TrigSerializeConverter.h.

                                                                    {
 
      // const std::string clname = ClassID_traits<DATA>::typeName(); // old style, does not contain the _vX of the actual class type
 
      // const std::type_info& typeId = ClassID_traits<DATA>::typeId();  // damn, this line just does not compile
 
      std::string classname = SG::normalizedTypeinfoName( typeid(DATA) ); // normalized string representation (matching the names used in the root dictionary)
 
      DATA *d( 0 );
      SG::fromStorable( pObj, d );
      if( m_log->level() <= MSG::DEBUG ) {
         *m_log << MSG::DEBUG << "My createRep for " << classID() << " "
                << classname << endmsg;
      }
 
      std::unique_ptr<DATA> holder;
      void *serptr = TrigSerialization::prepareForWrite (d, holder);
      //
      std::vector< uint32_t > ser;
 
      StatusCode sc = m_convHelper->createRep( classname, serptr, ser, typeIsxAOD);
      if( m_log->level() <= MSG::DEBUG ) {
         *m_log << MSG::DEBUG << "convHelper " << ser.size() << endmsg;
      }
 
      TrigStreamAddress *addr = new TrigStreamAddress( classID(), classname, "" );
      addr->add( ser );
 
      pAddr = addr;
 
      return sc;
   }

initialize()

template<typename DATA>

StatusCode TrigSerializeConverter< DATA >::initialize ( )

inline

Definition at line 240 of file TrigSerializeConverter.h.

                          {
 
      if( m_log ) {
         m_log->setLevel( msgSvc()->outputLevel( "TrigSerializeConverter" ) );
      } else {
         return StatusCode::FAILURE;
      }
 
      StatusCode sc = m_convHelper.retrieve();
      if( ! sc.isSuccess() ) {
         if( m_log->level() <= MSG::DEBUG ) {
            *m_log << MSG::DEBUG << m_convHelper << "not retrieved" << endmsg;
         }
      } else {
         if( m_log->level() <= MSG::DEBUG ) {
            *m_log << MSG::DEBUG << m_convHelper << "m_serializer retrieved"
                   << endmsg;
         }
         if( ! m_convHelper->initialize().isSuccess() ) {
            if( m_log->level() <= MSG::DEBUG ) {
               *m_log << MSG::DEBUG << m_convHelper << " not initialized"
                      << endmsg;
            }
         }
      }
 
      StatusCode scsg = m_sgHandle.retrieve();
      if (scsg.isFailure())
         *m_log << MSG::ERROR << "cannot access SG" << endmsg;
      else {
         m_sgsvc = m_sgHandle.operator->();
      }
      return sc;
   }

ipCorr()

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

inlineprotectedinherited

Definition at line 257 of file Converter.h.

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

phiCorr()

double Converter::phiCorr ( double phi )

inlineprotectedinherited

Definition at line 248 of file Converter.h.

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

repSvcType()

template<typename DATA>

long int TrigSerializeConverter< DATA >::repSvcType ( ) const

virtual

Definition at line 306 of file TrigSerializeConverter.h.

                                                        {
   return storageType();
}

selectTracks() [1/3]

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

inlineinherited

Definition at line 177 of file Converter.h.

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

selectTracks() [2/3]

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

inlineinherited

Definition at line 111 of file Converter.h.

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

selectTracks() [3/3]

void Converter::selectTracks ( const TrigInDetTrackCollection * trigtracks )

inlineinherited

Definition at line 56 of file Converter.h.

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

storageType()

template<typename DATA>

long int TrigSerializeConverter< DATA >::storageType ( )

static

Definition at line 302 of file TrigSerializeConverter.h.

                                                   {
   return TrigSerializeCnvSvc::storageType();
}

tracks()

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

inlineinherited

Definition at line 50 of file Converter.h.

{return m_tracks;}

Member Data Documentation

m_beamX

double Converter::m_beamX

protectedinherited

Definition at line 320 of file Converter.h.

m_beamY

double Converter::m_beamY

protectedinherited

Definition at line 320 of file Converter.h.

m_beamZ

double Converter::m_beamZ

protectedinherited

Definition at line 320 of file Converter.h.

m_convHelper

template<typename DATA>

ToolHandle< ITrigSerializeConvHelper > TrigSerializeConverter< DATA >::m_convHelper

private

Definition at line 294 of file TrigSerializeConverter.h.

m_log

template<typename DATA>

MsgStream* TrigSerializeConverter< DATA >::m_log

private

Definition at line 296 of file TrigSerializeConverter.h.

m_sgHandle

template<typename DATA>

ServiceHandle< StoreGateSvc > TrigSerializeConverter< DATA >::m_sgHandle

private

Definition at line 295 of file TrigSerializeConverter.h.

m_sgsvc

template<typename DATA>

StoreGateSvc* TrigSerializeConverter< DATA >::m_sgsvc

private

Definition at line 297 of file TrigSerializeConverter.h.

m_tracks

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

protectedinherited

Definition at line 321 of file Converter.h.

typeIsxAOD

template<typename DATA>

const bool TrigSerializeConverter< DATA >::typeIsxAOD = isxAOD<DATA>::value

static

Definition at line 142 of file TrigSerializeConverter.h.

---

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

• TrigSerializeConverter.h