TrackParticleTraits.h

Go to the documentation of this file.

// emacs: this is -*- c++ -*-
 
#ifndef  TRACKPARTICLETRAITS_H
#define  TRACKPARTICLETRAITS_H
 
#include "TrackTraits.h"
 
#include "xAODTracking/TrackParticleContainer.h"
 
 
template<>
struct TrackTraits<xAOD::TrackParticle> : public TrackTraitsBase<xAOD::TrackParticle> {
  
public:
  
  static uint8_t summary_value( const xAOD::TrackParticle& t, xAOD::SummaryType type ) {
    uint8_t n(0); t.summaryValue(n, type); return n;
  }
 
  static float  theta(const xAOD::TrackParticle& t) { return t.theta();  }
  static float qOverP(const xAOD::TrackParticle& t) { return t.qOverP(); }
  
  static float     z0(const xAOD::TrackParticle& t) { return t.z0()+t.vz(); }
  static float     d0(const xAOD::TrackParticle& t) { return t.d0(); }
 
  static float      E(const xAOD::TrackParticle& t) { return t.e(); }
  static float     ET(const xAOD::TrackParticle& t) { return t.e()*t.theta(); }
 
  static float      p(const xAOD::TrackParticle& t) { 
    if ( t.qOverP()==0 ) throw std::runtime_error( "divide by 0 - track qOverP() is 0" );
    return 1/t.qOverP();
  }
 
  
  static float     dphi(const xAOD::TrackParticle& t) { return std::sqrt(t.definingParametersCovMatrix()(Trk::phi0,Trk::phi0)); }
  static float   dtheta(const xAOD::TrackParticle& t) { return std::sqrt(t.definingParametersCovMatrix()(Trk::theta,Trk::theta)); }
  static float  dqOverP(const xAOD::TrackParticle& t) { return std::sqrt(t.definingParametersCovMatrix()(Trk::qOverP,Trk::qOverP)); }
  
  static float     deta(const xAOD::TrackParticle& t) { return 0.5*dtheta(t)/(std::cos(0.5*theta(t))*std::cos(0.5*theta(t))*std::tan(0.5*theta(t))); }  // ???? CHECK THIS <<--
 
  static float      dz0(const xAOD::TrackParticle& t) { return std::sqrt(t.definingParametersCovMatrix()(Trk::z0,Trk::z0)); }
  static float      dd0(const xAOD::TrackParticle& t) { return std::sqrt(t.definingParametersCovMatrix()(Trk::d0,Trk::d0)); }
  
  static float dpt(const xAOD::TrackParticle& t) {
    double th  = theta(t);
    double dth = dtheta(t);
    double sintheta = std::sin(th);
    double costheta = std::cos(th);
    double p_ = p(t);
    double dqovp = dqOverP(t);
    double dpt2 = (p_*p_*sintheta)*(p_*p_*sintheta)*dqovp*dqovp + (p_*costheta)*(p_*costheta)*dth*dth - 2*(p_*p_*sintheta)*(p_*costheta)*covthetaOvP(t);    
    if ( dpt2>0 ) return std::sqrt( dpt2 );
    return 0;
  }
  
  static float covthetaOvP(const xAOD::TrackParticle& t) { return t.definingParametersCovMatrix()(Trk::qOverP,Trk::theta); }; 
  
  static float chi2(const xAOD::TrackParticle& t) { return t.chiSquared(); }
  static float ndof(const xAOD::TrackParticle& t) { return t.numberDoF();  }  
  
  static uint8_t hasValidTime(const xAOD::TrackParticle& t) { return t.hasValidTime(); }
  static float           time(const xAOD::TrackParticle& t) { return t.time(); }

  //  static float      pTsig(const xAOD::TrackParticle& t) { return t.charge() ? std::copysign( t.pt(), t.charge() ) : 0.; }
  //  static float z0SinTheta(const xAOD::TrackParticle& t) { return z0( p ) * std::sin( theta( p ) ); }
  


 
  static int       nPixels( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfPixelHits );       }  
  static int nPixelsShared( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfPixelSharedHits ); }
  static int   nPixelHoles( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfPixelHoles );      }
 
  static int       nPixelsInner( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfInnermostPixelLayerHits ); }
  static int nPixelsInnerEndcap( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfInnermostPixelLayerEndcapHits ); }
 
  static int       nPixelsInnerShared( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfInnermostPixelLayerSharedHits ); }
  static int nPixelsInnerSharedEndcap( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfInnermostPixelLayerSharedEndcapHits ); }
  
  static int       nPixelsBlayer( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfNextToInnermostPixelLayerHits ); }
  static int nPixelsBlayerEndcap( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfNextToInnermostPixelLayerEndcapHits ); }

    
  static int       nSCT( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfSCTHits ); }
  static int nSCTShared( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfSCTSharedHits ); }
  static int  nSCTHoles( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfSCTHoles ); }
 
  static int      nSi( const xAOD::TrackParticle& t ) { return nPixels(t)+nSCT(t); }
  static int nSiHoles( const xAOD::TrackParticle& t ) { return nPixelHoles(t)+nSCTHoles(t); } 
  
  static int   nTRT( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfTRTHits ); }
  static int nTRTHi( const xAOD::TrackParticle& t ) { return summary_value( t, xAOD::numberOfTRTHighThresholdHitsTotal ); }
 
};
 
 
 
 
#endif  // TRACKPARTICLETRAITS_H