---

Jet.h

Go to the documentation of this file.

// ================================================
// Jet class description
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS STANDARD FORMAT
//
// Namespace Atlfast::
//
// class: Jet
//
// Authors: P.Clarke, H.Phillips, E.Richter-Was, P.Sherwood, R.Steward
//
// Description: 
//
//  Class to represent Jets constructed from Clusters
//
//..................................................


#ifndef __Atlfast_Jet__
#define __Atlfast_Jet__


// Kinematic classes.
#include "CLHEP/Vector/LorentzVector.h"
#include "AtlfastCode/IKinematic.h"

// Gaudi includes
#include "GaudiKernel/ContainedObject.h"
#include "GaudiKernel/ObjectVector.h"
#include "GaudiKernel/SmartRefVector.h"  

// Atlfast includes
#include "AtlfastCode/Cluster.h"
#include "AtlfastCode/HepMC_helper.h"

//Class ID for Jet class
static const CLID CLID_ATLFAST_Jet=2304 ;


namespace Atlfast {

  //........................................
  class Jet : 
    public IKinematic, 
    public ContainedObject, 
    public AssociationManager
    {

    private:

      HepLorentzVector m_momentum ;
      int m_flavor;

      /*#  Cluster lnkCluster; */
    public:
      
      // Constructors
      
      Jet() : ContainedObject(), AssociationManager() {}

      Jet( const Cluster& cluster):
        ContainedObject(), 
        AssociationManager(), 

        m_momentum(cluster.momentum()),
        m_flavor(HepMC_helper::UQUARK){

            this->associate(&cluster);
      }

      Jet( const HepLorentzVector& vec, Cluster& cluster ): 
        ContainedObject(), 
        AssociationManager(), 

        m_momentum(vec),
        m_flavor(HepMC_helper::UQUARK)

        {
          // Here we fill the related cells vector with pointers to the cells
          // from which it was constructed.
          this->associate( &cluster );
        }
      void setMomentum(const HepLorentzVector& vec) {m_momentum=vec;}
      // Other Athena mandatory methods
      static const CLID& classID() { return CLID_ATLFAST_Jet; }
      virtual const CLID& clID() const  {  return CLID_ATLFAST_Jet; }

      // To implement IKinematic
      virtual HepLorentzVector momentum() const {return m_momentum; } 
      virtual double  eta() const { return m_momentum.pseudoRapidity() ;}
      virtual double  phi() const { return m_momentum.phi() ;  } 
      virtual double  pT()  const { return m_momentum.perp() ;}
      //  FIXME: not sure if this is correct
      virtual double  eT()  const { 
        return m_momentum.e()*m_momentum.perp()/m_momentum.rho() ;} 
      virtual double  mT()  const { return m_momentum.mt(); } 
      virtual double  px()  const { return m_momentum.px(); } 
      virtual double  py()  const { return m_momentum.py(); } 
      virtual double  pz()  const { return m_momentum.pz(); } 
      
      double  e()   const { return m_momentum.e(); } 

      //methods to set and interogate flavor tag  
      void setBTag()  {m_flavor=HepMC_helper::BQUARK;}
      void setCTag()  {m_flavor=HepMC_helper::CQUARK;}
      void setTauTag(int tauTag){m_flavor=tauTag;}
      bool isBTag()   {return m_flavor==HepMC_helper::BQUARK;}
      bool isCTag()   {return m_flavor==HepMC_helper::CQUARK;}
      bool isTauTag() {return m_flavor==HepMC_helper::TAU;}
      int pdg_id() {return m_flavor;}

      // Input streamer
      virtual StreamBuffer& serialize(StreamBuffer& s);

      // Output streamer
      virtual StreamBuffer& serialize(StreamBuffer& s) const;
    };

  // ..............................
  // This typedef is used by JetMaker to store Jets
  //
  // It is put here rather than JetMaker.h so that any algorithm wishing 
  // to use Jets from the TES need only #include Jet.h
  
  typedef ObjectVector<Jet> t_jetCollection ;
  
  
} // end of namespace 


#endif

---