#include <TGC_IDDetDescrCnv.h>

Inheritance diagram for TGC_IDDetDescrCnv:

Collaboration diagram for TGC_IDDetDescrCnv:

Public Member Functions
	TGC_IDDetDescrCnv (ISvcLocator *svcloc)
virtual StatusCode	initialize () override
virtual StatusCode	createObj (IOpaqueAddress pAddr, DataObject &pObj) override
virtual long int	repSvcType () const override
virtual StatusCode	fillObjRefs (IOpaqueAddress pAddr, DataObject pObj)
virtual StatusCode	createRep (DataObject pObj, IOpaqueAddress &pAddr)
virtual StatusCode	fillRepRefs (IOpaqueAddress pAddr, DataObject pObj)
const ServiceHandle< StoreGateSvc > &	detStore () const
	Handle to DetectorStore.
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
ServiceHandle< StoreGateSvc >	m_detStore
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
	MessageSvc pointer.
std::atomic< MSG::Level >	m_lvl
	Current logging level.
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE
	Messaging initialized (initMessaging).

Detailed Description

Converter for TgcIdHelper.

Definition at line 15 of file TGC_IDDetDescrCnv.h.

Constructor & Destructor Documentation

◆ TGC_IDDetDescrCnv()

TGC_IDDetDescrCnv::TGC_IDDetDescrCnv ( ISvcLocator * svcloc )

inline

Definition at line 17 of file TGC_IDDetDescrCnv.h.

17 :

18 T_Muon_IDDetDescrCnv(svcloc, "TGC_IDDetDescrCnv") {}

T_Muon_IDDetDescrCnv< TgcIdHelper >::T_Muon_IDDetDescrCnv

T_Muon_IDDetDescrCnv(ISvcLocator *svcloc, const char *name)

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

inlinestaticinherited

Definition at line 34 of file T_Muon_IDDetDescrCnv.h.

34{ return ClassID_traits<IDHELPER>::ID(); }

ClassID_traits::ID

static constexpr CLID ID()

Definition Control/AthenaKernel/AthenaKernel/ClassID_traits.h:44

◆ clear()

void Converter::clear ( )

inlineinherited

Definition at line 51 of file Converter.h.

51{m_tracks.clear();}

◆ createObj()

virtual StatusCode T_Muon_IDDetDescrCnv< TgcIdHelper >::createObj	(	IOpaqueAddress *	pAddr,
		DataObject *&	pObj )

overridevirtualinherited

Implements DetDescrConverter.

◆ createRep()

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

virtualinherited

Definition at line 29 of file DetDescrConverter.cxx.

                                                            {
    return StatusCode::SUCCESS;
}

◆ detStore()

const ServiceHandle< StoreGateSvc > & DetDescrConverter::detStore ( ) const

inlineinherited

Handle to DetectorStore.

Definition at line 33 of file DetDescrConverter.h.

33{ return m_detStore; }

T_Muon_IDDetDescrCnv< TgcIdHelper >::m_detStore

ServiceHandle< StoreGateSvc > m_detStore

Definition DetDescrConverter.h:42

◆ fillObjRefs()

StatusCode DetDescrConverter::fillObjRefs	(	IOpaqueAddress *	pAddr,
		DataObject *	pObj )

virtualinherited

Definition at line 28 of file DetDescrConverter.cxx.

                                                         {
    return StatusCode::SUCCESS;
}

◆ fillRepRefs()

StatusCode DetDescrConverter::fillRepRefs	(	IOpaqueAddress *	pAddr,
		DataObject *	pObj )

virtualinherited

Definition at line 30 of file DetDescrConverter.cxx.

                                                         {
    return StatusCode::SUCCESS;
}

◆ initialize()

virtual StatusCode T_Muon_IDDetDescrCnv< TgcIdHelper >::initialize ( )

overridevirtualinherited

◆ initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 121 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 92 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 99 of file AthMessaging.h.

183{ return msg() << lvl; }

T_Muon_IDDetDescrCnv< TgcIdHelper >::msg

MsgStream & msg() const

Definition AthMessaging.h:167

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

◆ 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 int T_Muon_IDDetDescrCnv< TgcIdHelper >::repSvcType ( ) const

inlineoverridevirtualinherited

Definition at line 30 of file T_Muon_IDDetDescrCnv.h.

30{ return storageType(); }

T_Muon_IDDetDescrCnv< TgcIdHelper >::storageType

static long storageType()

Definition T_Muon_IDDetDescrCnv.h:33

◆ 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 108 of file AthMessaging.cxx.

{
  m_lvl = lvl;
}

◆ storageType()

long T_Muon_IDDetDescrCnv< TgcIdHelper >::storageType ( )

inlinestaticinherited

Definition at line 33 of file T_Muon_IDDetDescrCnv.h.

33{ return DetDescr_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

◆ ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE

mutableprivateinherited

Messaging initialized (initMessaging).

Definition at line 141 of file AthMessaging.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_detStore

ServiceHandle<StoreGateSvc> DetDescrConverter::m_detStore

privateinherited

Definition at line 42 of file DetDescrConverter.h.

◆ m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

135{ nullptr };

◆ m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

138{ 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_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.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 file:

TGC_IDDetDescrCnv.h

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TGC_IDDetDescrCnv()

Member Function Documentation

◆ addTrack()

◆ classID()

◆ clear()

◆ createObj()

◆ createRep()

◆ detStore()

◆ fillObjRefs()

◆ fillRepRefs()

◆ initialize()

◆ initMessaging()

◆ ipCorr()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ phiCorr()

◆ repSvcType()

◆ selectTracks() [1/2]

◆ selectTracks() [2/2]

◆ setLevel()

◆ storageType()

◆ tracks()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ m_beamX

◆ m_beamY

◆ m_beamZ

◆ m_detStore

◆ m_imsg

◆ m_lvl

◆ m_msg_tls

◆ m_nm

◆ m_tracks