|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
#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 "FWCore/MessageLogger/interface/MessageLogger.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); }
|