GsfElectronCore.cc

Go to the documentation of this file.

#include "DataFormats/EgammaCandidates/interface/GsfElectronCore.h"
#include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
#include "DataFormats/EgammaReco/interface/ElectronSeedFwd.h"
#include "DataFormats/EgammaReco/interface/ElectronSeed.h"

#include <cmath>

using namespace reco ;

GsfElectronCore::GsfElectronCore
 ()
 : ctfGsfOverlap_(0.), isEcalDrivenSeed_(false), isTrackerDrivenSeed_(false)
 {}

 GsfElectronCore::GsfElectronCore
  ( const GsfTrackRef & gsfTrack )
  : gsfTrack_(gsfTrack), ctfGsfOverlap_(0.), isEcalDrivenSeed_(false), isTrackerDrivenSeed_(false)
  {
   edm::RefToBase<TrajectorySeed> seed = gsfTrack_->extra()->seedRef() ;
   if (seed.isNull())
    { edm::LogError("GsfElectronCore")<<"The GsfTrack has no seed ?!" ; }
   else
    {
     ElectronSeedRef elseed = seed.castTo<ElectronSeedRef>() ;
     if (elseed.isNull())
      { edm::LogError("GsfElectronCore")<<"The GsfTrack seed is not an ElectronSeed ?!" ; }
     else
      {
       if (elseed->isEcalDriven()) isEcalDrivenSeed_ = true ;
       if (elseed->isTrackerDriven()) isTrackerDrivenSeed_ = true ;
      }
    }
  }

GsfElectronCore * GsfElectronCore::clone() const
 { return new GsfElectronCore(*this) ; }