---

Jet.cxx

Go to the documentation of this file.

// Jet implementation
//

#include "AtlfastEvent/Jet.h"
#include "AtlfastEvent/ParticleCodes.h"
//#include "AtlfastUtils/HepMC_helper/IMCselector.h"

#include "GaudiKernel/ObjectList.h"
#include "GaudiKernel/ObjectVector.h"
#include "GaudiKernel/ObjectFactory.h"

static ContainedObjectFactory < Atlfast::Jet > s_AtlfastJetFactory;
const IFactory& AtlfastJetFactory = s_AtlfastJetFactory;

namespace Atlfast {
  Jet::Jet( const Cluster& cluster):
    ContainedObject(), 
    AssociationManager(), 
    
    m_momentum(cluster.momentum()),
    m_flavor(ParticleCodes::UQUARK),
    m_tagged(true){
    
    this->associate(&cluster);
  }

  Jet::Jet( const HepLorentzVector& vec, Cluster& cluster ): 
    ContainedObject(), 
    AssociationManager(), 
    
    m_momentum(vec),
    m_flavor(ParticleCodes::UQUARK),
    m_tagged(true){
    // Here we fill the related cells vector with pointers to the cells
    // from which it was constructed.
    this->associate( &cluster );
  }

  void Jet::setMomentum(const HepLorentzVector& vec) {m_momentum=vec;}

  HepLorentzVector Jet::momentum() const {return m_momentum; } 
  double  Jet::eta() const { return m_momentum.pseudoRapidity() ;}
  double  Jet::phi() const { return m_momentum.phi() ;  } 
  double  Jet::pT()  const { return m_momentum.perp() ;}
  double  Jet::eT()  const { 
    return m_momentum.e()*m_momentum.perp()/m_momentum.rho() ;} 
  double  Jet::mT()  const { return m_momentum.mt(); } 
  double  Jet::px()  const { return m_momentum.px(); } 
  double  Jet::py()  const { return m_momentum.py(); } 
  double  Jet::pz()  const { return m_momentum.pz(); } 
      
  double  Jet::e()   const { return m_momentum.e(); } 
      
  //methods to set and interogate flavor tag  
  void Jet::setBTag() {
    m_flavor=ParticleCodes::BQUARK;
    m_tagged=true;
  }
  void Jet::setCTag() {
    m_flavor=ParticleCodes::CQUARK;
    m_tagged=true;
  }

 void Jet::setLightTag() {
    m_flavor=ParticleCodes::UQUARK;
    m_tagged=true;
  }

  void Jet::setTauTag(int tauTag){
    m_flavor=tauTag;
    m_tagged=true;
  }
  void Jet::setUntagged() {
    m_flavor=0;
    m_tagged=false;
  }

  bool Jet::isBTag()   const {
    return m_flavor==ParticleCodes::BQUARK && m_tagged;
  }

  bool Jet::isCTag()   const{
    return m_flavor==ParticleCodes::CQUARK && m_tagged;
  }

  bool Jet::isLightTag()   const{
    return m_flavor==ParticleCodes::UQUARK && m_tagged;
  }

  bool Jet::isTauTag() const {
    return m_flavor==ParticleCodes::TAU    && m_tagged;
  }
  bool Jet::isTagged() const {
    return m_tagged;
  }

  int Jet::pdg_id()    const {return m_flavor;}

} // end of namespace bracket

---