Converter Class Reference

#include <Converter.h>

Inherited by TBDataCnv< TBTDC >, TBDataCnv< TBTDCRawCont >, TBDataCnv< TBADCRawCont >, TBDataCnv< TBLArDigitContainer >, TBDataCnv< TBLArCalibDigitContainer >, TBDataCnv< TBBPCRawCont >, TBDataCnv< TBScintillatorRawCont >, TBDataCnv< TBTriggerPatternUnit >, TBDataCnv< TBTailCatcherRaw >, TBDataCnv< TBEventInfo >, ToyConverter< Foo >, ToyConverter< Bar >, AthConstConverter, Athena::LeafCnv, Athena::RootCnv, Athena::xAODCnv, AthenaAttributeListCnv, AthenaPoolConverter, CTPByteStreamCnv, CondAttrListCollCnv, CondAttrListVecCnv, DetDescrConverter, EventInfoByteStreamAuxCnv, EventInfoByteStreamxAODCnv, EventInfoMuonCalibStreamCnv, HLT::HLTResultByteStreamCnv, HLT::HLTResultMTByteStreamCnv, L1TriggerResultByteStreamCnv, LUCID_DigitByteStreamCnv, LVL1BS::CmxCpHitsByteStreamAuxCnv, LVL1BS::CmxCpHitsByteStreamxAODCnv, LVL1BS::CmxCpTobByteStreamAuxCnv, LVL1BS::CmxCpTobByteStreamxAODCnv, LVL1BS::CmxEtSumsByteStreamAuxCnv, LVL1BS::CmxEtSumsByteStreamxAODCnv, LVL1BS::CmxJetHitsByteStreamAuxCnv, LVL1BS::CmxJetHitsByteStreamxAODCnv, LVL1BS::CmxJetTobByteStreamAuxCnv, LVL1BS::CmxJetTobByteStreamxAODCnv, LVL1BS::CmxRoIByteStreamAuxCnv, LVL1BS::CmxRoIByteStreamxAODCnv, LVL1BS::CpmTobRoiByteStreamAuxCnv, LVL1BS::CpmTobRoiByteStreamxAODCnv, LVL1BS::CpmTowerByteStreamAuxCnv, LVL1BS::CpmTowerByteStreamxAODCnv, LVL1BS::JemEtSumsByteStreamAuxCnv, LVL1BS::JemEtSumsByteStreamxAODCnv, LVL1BS::JemTobRoiByteStreamAuxCnv, LVL1BS::JemTobRoiByteStreamxAODCnv, LVL1BS::JetElementByteStreamAuxCnv, LVL1BS::JetElementByteStreamxAODCnv, LVL1BS::L1CaloByteStreamAuxCnv< ContainerT, AuxContainerT >, LVL1BS::L1TopoByteStreamAuxCnv, LVL1BS::L1TopoByteStreamxAODCnv, LVL1BS::PpmByteStreamAuxCnv, LVL1BS::PpmByteStreamxAODCnv, LVL1BS::RodHeaderByteStreamAuxCnv, LVL1BS::RodHeaderByteStreamxAODCnv, McEventCnv, MuCTPIByteStreamCnv, RecCTPByteStreamCnv, RecMuCTPIByteStreamCnv, RoIBResultByteStreamCnv< ROBF >, TBDataCnv< DATATYPE >, ToyConverter< DATA >, TrigSerializeConverter< DATA >, and xAODPrivate::Loader.

Collaboration diagram for Converter:

Public Member Functions
	Converter ()
	Converter (double x, double y, double z)
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)

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

Detailed Description

Definition at line 27 of file Converter.h.

Constructor & Destructor Documentation

Converter() [1/2]

Converter::Converter ( )

inline

Definition at line 31 of file Converter.h.

              {
    m_beamX = 0;
    m_beamY = 0;
    m_beamZ = 0;
  }

Converter() [2/2]

Converter::Converter ( double x, double y, double z )

inline

Definition at line 37 of file Converter.h.

                                          {
    m_beamX = x;
    m_beamY = y;
    m_beamZ = z;
  }

Member Function Documentation

addTrack()

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

inline

Definition at line 45 of file Converter.h.

                              {
    m_tracks.push_back(t);
  }

clear()

void Converter::clear ( )

inline

Definition at line 53 of file Converter.h.

{m_tracks.clear();} 

ipCorr()

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

inlineprotected

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 )

inlineprotected

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

selectTracks() [1/3]

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

inline

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 )

inline

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 )

inline

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

tracks()

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

inline

Definition at line 50 of file Converter.h.

{return m_tracks;}

Member Data Documentation

m_beamX

double Converter::m_beamX

protected

Definition at line 320 of file Converter.h.

m_beamY

double Converter::m_beamY

protected

Definition at line 320 of file Converter.h.

m_beamZ

double Converter::m_beamZ

protected

Definition at line 320 of file Converter.h.

m_tracks

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

protected

Definition at line 321 of file Converter.h.

---

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

• Converter.h