CastorEgamma.h

Go to the documentation of this file.

#ifndef CastorReco_CastorEgamma_h
#define CastorReco_CastorEgamma_h

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

namespace reco {

  class CastorEgamma : public CastorCluster {
  public:

    CastorEgamma() : energycal_(0.) { }

    CastorEgamma(const double energycal, const CastorClusterRef& usedCluster);

    virtual ~CastorEgamma();

    double energy() const { return (*usedCluster_).energy(); }

    double energycal() const { return energycal_; }

    ROOT::Math::XYZPoint position() const { return (*usedCluster_).position(); }
    
    CastorClusterRef getUsedCluster() const { return usedCluster_; }

    bool operator >=(const CastorEgamma& rhs) const { return (energycal_>=rhs.energycal_); }

    bool operator > (const CastorEgamma& rhs) const { return (energycal_> rhs.energycal_); }

    bool operator <=(const CastorEgamma& rhs) const { return (energycal_<=rhs.energycal_); }

    bool operator < (const CastorEgamma& rhs) const { return (energycal_< rhs.energycal_); }

    double emEnergy() const { return (*usedCluster_).emEnergy(); }

    double hadEnergy() const { return (*usedCluster_).hadEnergy(); }

    double fem() const { return (*usedCluster_).fem(); }

    double width() const { return (*usedCluster_).width(); }

    double depth() const { return (*usedCluster_).depth(); }

    double fhot() const { return (*usedCluster_).fhot(); }

    double sigmaz() const { return (*usedCluster_).sigmaz(); }

    double eta() const { return (*usedCluster_).eta(); }

    double phi() const { return (*usedCluster_).phi(); }

    double x() const { return (*usedCluster_).x(); }

    double y() const { return (*usedCluster_).y(); }

    double rho() const { return (*usedCluster_).rho(); }

  private:

    double energycal_;

    CastorClusterRef usedCluster_;
  };
  
  // define CastorEgammaCollection
  typedef std::vector<CastorEgamma> CastorEgammaCollection;

}

#endif