TauJetTraits.h

Go to the documentation of this file.

// emacs: this is -*- c++ -*-
 
#ifndef  IDTPM2_TAUJETTRAITS_H
#define  IDTPM2_TAUJETTRAITS_H
 
#include "xAODTau/TauJetContainer.h"
 
#include "TrackTraits.h"
 
#include <cmath>
#include <iostream>

 
template<>
struct TrackTraits<xAOD::TauJet> : public TrackTraitsBase<xAOD::TauJet> {
 
private:
  
  inline static const SG::ConstAccessor<xAOD::TauTrack_v1::TrackParticleLinks_t> m_trackAcc = {"trackLinks"};
  
public:
  
  static float    eta(const xAOD::TauJet& t) { return t.eta();    } 
  static float    phi(const xAOD::TauJet& t) { return t.phi();    } 

  static float     pt(const xAOD::TauJet& t) { return t.pt(); }
 
  static float  theta(const xAOD::TauJet& t) { return 2*std::atan(std::exp(-t.eta()));  }
  static float qOverP(const xAOD::TauJet& t) { return t.charge()*std::sin(theta(t))/std::fabs(t.pt()); }
  
  static float      E(const xAOD::TauJet& t) { return t.e(); }
 
  static float     ET(const xAOD::TauJet& t) { return t.e()*std::sin(theta(t)); }
 
  static float      m(const xAOD::TauJet& t) { return t.m(); }
 
  static float      p(const xAOD::TauJet& t) {
    double th = theta(t);
    if ( th==0 ) throw std::runtime_error( "divide by 0 - tau theta is 0" );
    return t.pt()/std::sin(th);
  }
 
  static float charge(const xAOD::TauJet& t) { return t.charge(); }

  template<typename U>
  static void addChildren( const xAOD::TauJet& t, std::vector<U>& v) {
    auto tt = t.tracks();
    for ( size_t i=0 ; i<tt.size() ; i++ ) {
      const auto& links = m_trackAcc(*tt[i]);
      if ( !links.empty() && links[0].isValid() ) v.emplace_back(tt[i]->track());
      //  else std::cerr << "duff TauTrack ElementLink for: " << tt[i] << std::endl;
    }
  }
  
};
 
 
 
 
#endif  // IDTPM2_TAUJETTRAITS_H