CastorTower.h

Go to the documentation of this file.

#ifndef CastorReco_CastorTower_h
#define CastorReco_CastorTower_h

#include <vector>
#include <memory>
#include "DataFormats/Math/interface/Point3D.h"

#include "DataFormats/Common/interface/RefProd.h"
#include "DataFormats/Common/interface/Ref.h"
#include "DataFormats/Common/interface/RefVector.h"
#include "DataFormats/Common/interface/SortedCollection.h"

#include "DataFormats/HcalRecHit/interface/CastorRecHit.h"
//#include "DataFormats/HcalRecHit/interface/HcalRecHitDefs.h"
#include "DataFormats/Candidate/interface/LeafCandidate.h"

namespace reco {

  class CastorTower : public LeafCandidate {
  public:
    typedef edm::SortedCollection<CastorRecHit> CastorRecHitCollection;
    typedef edm::Ref<CastorRecHitCollection> CastorRecHitRef;
    typedef edm::RefVector<CastorRecHitCollection> CastorRecHitRefs;

    // default constructor. Set energy and position to zero

    CastorTower()
        : energy_(0.),
          position_(ROOT::Math::XYZPoint(0., 0., 0.)),
          emEnergy_(0.),
          hadEnergy_(0.),
          fem_(0.),
          depth_(0.),
          fhot_(0.) {}

    // constructor from values
    CastorTower(const double energy,
                const ROOT::Math::XYZPoint& position,
                const double emEnergy,
                const double hadEnergy,
                const double fem,
                const double depth,
                const double fhot,
                const CastorRecHitRefs& usedRecHits);

    ~CastorTower() override;

    ROOT::Math::XYZPoint position() const { return position_; }

    double emEnergy() const { return emEnergy_; }

    double hadEnergy() const { return hadEnergy_; }

    double fem() const { return fem_; }

    double depth() const { return depth_; }

    double fhot() const { return fhot_; }

    CastorRecHitRefs getUsedRecHits() const { return usedRecHits_; }

    CastorRecHitRefs::iterator rechitsBegin() const { return usedRecHits_.begin(); }

    CastorRecHitRefs::iterator rechitsEnd() const { return usedRecHits_.end(); }

    size_t rechitsSize() const { return usedRecHits_.size(); }

    void add(const CastorRecHitRef& rechit) { usedRecHits_.push_back(rechit); }

    bool operator>=(const CastorTower& rhs) const { return (energy_ >= rhs.energy_); }

    bool operator>(const CastorTower& rhs) const { return (energy_ > rhs.energy_); }

    bool operator<=(const CastorTower& rhs) const { return (energy_ <= rhs.energy_); }

    bool operator<(const CastorTower& rhs) const { return (energy_ < rhs.energy_); }

    double xPos() const { return position_.x(); }

    double yPos() const { return position_.y(); }

    double rho() const { return position_.rho(); }

  private:
    double energy_;

    ROOT::Math::XYZPoint position_;

    double emEnergy_;

    double hadEnergy_;

    double fem_;

    double depth_;

    double fhot_;

    CastorRecHitRefs usedRecHits_;
  };

  typedef std::vector<CastorTower> CastorTowerCollection;

  // persistent reference to CastorTower objects
  typedef edm::Ref<CastorTowerCollection> CastorTowerRef;

  typedef edm::RefVector<CastorTowerCollection> CastorTowerRefVector;

  typedef CastorTowerRefVector::iterator CastorTower_iterator;

}  // namespace reco

#endif