CastorCluster.h

Go to the documentation of this file.

#ifndef CastorReco_CastorCluster_h
#define CastorReco_CastorCluster_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/CastorReco/interface/CastorTower.h"

namespace reco {

  class CastorCluster {
  public:
    CastorCluster()
        : energy_(0.),
          position_(ROOT::Math::XYZPoint(0., 0., 0.)),
          emEnergy_(0.),
          hadEnergy_(0.),
          fem_(0.),
          width_(0.),
          depth_(0.),
          fhot_(0.),
          sigmaz_(0.) {}

    CastorCluster(const double energy,
                  const ROOT::Math::XYZPoint& position,
                  const double emEnergy,
                  const double hadEnergy,
                  const double fem,
                  const double width,
                  const double depth,
                  const double fhot,
                  const double sigmaz,
                  const CastorTowerRefVector& usedTowers);

    virtual ~CastorCluster();

    double energy() const { return energy_; }

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

    double emEnergy() const { return emEnergy_; }

    double hadEnergy() const { return hadEnergy_; }

    double fem() const { return fem_; }

    double width() const { return width_; }

    double depth() const { return depth_; }

    double fhot() const { return fhot_; }

    double sigmaz() const { return sigmaz_; }

    CastorTowerRefVector getUsedTowers() const { return usedTowers_; }

    CastorTower_iterator towersBegin() const { return usedTowers_.begin(); }

    CastorTower_iterator towersEnd() const { return usedTowers_.end(); }

    size_t towersSize() const { return usedTowers_.size(); }

    void add(const CastorTowerRef& tower) { usedTowers_.push_back(tower); }

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

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

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

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

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

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

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

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

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

  private:
    double energy_;

    ROOT::Math::XYZPoint position_;

    double emEnergy_;

    double hadEnergy_;

    double fem_;

    double width_;

    double depth_;

    double fhot_;

    double sigmaz_;

    CastorTowerRefVector usedTowers_;
  };

  typedef std::vector<CastorCluster> CastorClusterCollection;

  // persistent reference to CastorCluster objects
  typedef edm::Ref<CastorClusterCollection> CastorClusterRef;

  typedef edm::RefVector<CastorClusterCollection> CastorClusterRefVector;

  typedef CastorClusterRefVector::iterator CastorCluster_iterator;
}  // namespace reco

#endif