/data/refman/pasoursint/CMSSW_5_3_8_patch3/src/DataFormats/EgammaCandidates/interface/SiStripElectron.h

Go to the documentation of this file.

#ifndef EgammaCandidates_SiStripElectron_h
#define EgammaCandidates_SiStripElectron_h
// -*- C++ -*-
//
// Package:     EgammaCandidates
// Class  :     SiStripElectron
// 
//
// Original Author:  Jim Pivarski
//         Created:  Fri May 26 15:43:14 EDT 2006
// $Id: SiStripElectron.h,v 1.14 2008/11/03 17:33:44 nancy Exp $
//

// system include files

#include <vector>

// user include files

#include "DataFormats/RecoCandidate/interface/RecoCandidate.h"
#include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2DCollection.h"
#include "DataFormats/Common/interface/RefVector.h" 

// forward declarations

namespace reco {
  
  class SiStripElectron : public RecoCandidate {
  public:
    SiStripElectron() : RecoCandidate() { }
    SiStripElectron(const reco::SuperClusterRef& superCluster,
                    Charge q,
                    const std::vector<SiStripRecHit2D>& rphiRecHits,
                    const std::vector<SiStripRecHit2D>& stereoRecHits,
                    double superClusterPhiVsRSlope,
                    double phiVsRSlope,
                    double phiAtOrigin,
                    double chi2,
                    int ndof,
                    double pt,
                    double pz,
                    double zVsRSlope,
                    unsigned int numberOfStereoHits,
                    unsigned int numberOfBarrelRphiHits,
                    unsigned int numberOfEndcapZphiHits)
      : RecoCandidate(q, LorentzVector(pt*cos(phiAtOrigin), pt*sin(phiAtOrigin), pz, 
                                       sqrt( pt*pt + pz*pz + 0.000510*0.000510)), Point(0,0,0), -11 * q )
        , superCluster_(superCluster)
      , rphiRecHits_(rphiRecHits)
      , stereoRecHits_(stereoRecHits)
      , superClusterPhiVsRSlope_(superClusterPhiVsRSlope)
      , phiVsRSlope_(phiVsRSlope)
      , phiAtOrigin_(phiAtOrigin)
      , chi2_(chi2)
      , ndof_(ndof)
      , pt_(pt)
      , pz_(pz)
      , zVsRSlope_(zVsRSlope)
      , numberOfStereoHits_(numberOfStereoHits)
      , numberOfBarrelRphiHits_(numberOfBarrelRphiHits)
      , numberOfEndcapZphiHits_(numberOfEndcapZphiHits) { }
    
    SiStripElectron(const SiStripElectron& rhs)
      : RecoCandidate(rhs)
      , superCluster_(rhs.superCluster())
      , rphiRecHits_(rhs.rphiRecHits())
      , stereoRecHits_(rhs.stereoRecHits())
      , superClusterPhiVsRSlope_(rhs.superClusterPhiVsRSlope())
      , phiVsRSlope_(rhs.phiVsRSlope())
      , phiAtOrigin_(rhs.phiAtOrigin())
      , chi2_(rhs.chi2())
      , ndof_(rhs.ndof())
      , pt_(rhs.pt())
      , pz_(rhs.pz())
      , zVsRSlope_(rhs.zVsRSlope())
      , numberOfStereoHits_(rhs.numberOfStereoHits())
      , numberOfBarrelRphiHits_(rhs.numberOfBarrelRphiHits())
      , numberOfEndcapZphiHits_(rhs.numberOfEndcapZphiHits()) { }
    
    SiStripElectron( Charge q, const LorentzVector & p4, const Point & vtx = Point( 0, 0, 0 ) ) : 
      RecoCandidate( q, p4, vtx, -11 * q ) { }
    virtual ~SiStripElectron();
    virtual SiStripElectron * clone() const;
    virtual reco::SuperClusterRef superCluster() const;
    
    const std::vector<SiStripRecHit2D>& rphiRecHits() const { return rphiRecHits_; }
    const std::vector<SiStripRecHit2D>& stereoRecHits() const { return stereoRecHits_; }
    
    double superClusterPhiVsRSlope() const { return superClusterPhiVsRSlope_; }
    double phiVsRSlope() const { return phiVsRSlope_; }
    double phiAtOrigin() const { return phiAtOrigin_; }
    double chi2() const { return chi2_; }
    int ndof() const { return ndof_; }
    
    double pt() const { return pt_; }
    double pz() const { return pz_; }
    
    double zVsRSlope() const { return zVsRSlope_; }
    
    unsigned int numberOfStereoHits() const { return numberOfStereoHits_; }
    unsigned int numberOfBarrelRphiHits() const { return numberOfBarrelRphiHits_; }
    unsigned int numberOfEndcapZphiHits() const { return numberOfEndcapZphiHits_; }

    bool isElectron() const;
  private:
    virtual bool overlap( const Candidate & ) const;
    reco::SuperClusterRef superCluster_;
    std::vector<SiStripRecHit2D> rphiRecHits_;
    std::vector<SiStripRecHit2D> stereoRecHits_;
    
    double superClusterPhiVsRSlope_;
    double phiVsRSlope_;
    double phiAtOrigin_;
    double chi2_;
    int ndof_;
    
    double pt_;
    double pz_;
    
    double zVsRSlope_;
    
    unsigned int numberOfStereoHits_;
    unsigned int numberOfBarrelRphiHits_;
    unsigned int numberOfEndcapZphiHits_;
   };
}

#endif