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

#endif