CaloJet.h

Go to the documentation of this file.

#ifndef JetReco_CaloJet_h
#define JetReco_CaloJet_h

#include "DataFormats/JetReco/interface/Jet.h"
#include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"


namespace reco {
class CaloJet : public Jet {
 public:
  struct Specific {
    Specific () :
      mMaxEInEmTowers (0),
     mMaxEInHadTowers (0),
     mHadEnergyInHO (0),
     mHadEnergyInHB (0),
     mHadEnergyInHF (0),
     mHadEnergyInHE (0),
     mEmEnergyInEB (0),
     mEmEnergyInEE (0),
     mEmEnergyInHF (0),
     mEnergyFractionHadronic(0),
     mEnergyFractionEm (0),
     mTowersArea (0)
    {}
    
    float mMaxEInEmTowers;
    float mMaxEInHadTowers;
    float mHadEnergyInHO;
    float mHadEnergyInHB;
    float mHadEnergyInHF;
    float mHadEnergyInHE;
    float mEmEnergyInEB;
    float mEmEnergyInEE;
    float mEmEnergyInHF;
    float mEnergyFractionHadronic;
    float mEnergyFractionEm;
    float mTowersArea;
  };
  
  CaloJet() {}
  
  CaloJet(const LorentzVector& fP4, const Point& fVertex, const Specific& fSpecific, 
      const Jet::Constituents& fConstituents);

  CaloJet(const LorentzVector& fP4, const Point& fVertex, const Specific& fSpecific);

  CaloJet(const LorentzVector& fP4, const Specific& fSpecific, 
      const Jet::Constituents& fConstituents);

  
  virtual ~CaloJet() {};
  
  float maxEInEmTowers() const {return m_specific.mMaxEInEmTowers;}
  float maxEInHadTowers() const {return m_specific.mMaxEInHadTowers;}
  float energyFractionHadronic () const {return m_specific.mEnergyFractionHadronic;}
  float emEnergyFraction() const {return m_specific.mEnergyFractionEm;}
  float hadEnergyInHB() const {return m_specific.mHadEnergyInHB;}
  float hadEnergyInHO() const {return m_specific.mHadEnergyInHO;}
  float hadEnergyInHE() const {return m_specific.mHadEnergyInHE;}
  float hadEnergyInHF() const {return m_specific.mHadEnergyInHF;}
  float emEnergyInEB() const {return m_specific.mEmEnergyInEB;}
  float emEnergyInEE() const {return m_specific.mEmEnergyInEE;}
  float emEnergyInHF() const {return m_specific.mEmEnergyInHF;}
  float towersArea() const {return m_specific.mTowersArea;}
  int n90() const {return nCarrying (0.9);}
  int n60() const {return nCarrying (0.6);}

    //  float physicsEtaQuick (float fZVertex) const;
    //float physicsEta (float fZVertex) const {return physicsEtaQuick (fZVertex);}
    //LorentzVector physicsP4 (float fZVertex) const;
  LorentzVector physicsP4 (const Particle::Point &vertex) const;
  LorentzVector detectorP4 () const;
  
    //float physicsEtaDetailed (float fZVertex) const;

    //float detectorEta () const {return eta();}


  virtual CaloTowerPtr getCaloConstituent (unsigned fIndex) const;
  virtual std::vector <CaloTowerPtr> getCaloConstituents () const;
  
  // block accessors
  
  const Specific& getSpecific () const {return m_specific;}

  virtual CaloJet* clone () const;

  virtual std::string print () const;

  std::vector<CaloTowerDetId> getTowerIndices() const;
 private:
  virtual bool overlap( const Candidate & ) const;
  
  //Variables specific to to the CaloJet class
  Specific m_specific;
};
}
// temporary fix before include_checcker runs globally
#include "DataFormats/JetReco/interface/CaloJetCollection.h" //INCLUDECHECKER:SKIP
#endif