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:
    typedef CaloTowerPtr ConstituentTypePtr;
    typedef CaloTowerFwdPtr ConstituentTypeFwdPtr;

    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);

    ~CaloJet() override{};

    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; }

    CaloJet* clone() const override;

    std::string print() const override;

    std::vector<CaloTowerDetId> getTowerIndices() const;

  private:
    bool overlap(const Candidate&) const override;

    //Variables specific to to the CaloJet class
    Specific m_specific;
  };
}  // namespace reco
// temporary fix before include_checcker runs globally
#include "DataFormats/JetReco/interface/CaloJetCollection.h"  //INCLUDECHECKER:SKIP
#endif