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/HcalRecHitFwd.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);

    virtual ~CastorTower();

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

}

#endif