CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
ElectronEPcombinator Class Reference

#include <ElectronEPcombinator.h>

Public Member Functions

void combine (SimpleElectron &electron)
 
 ElectronEPcombinator ()
 
void setCombinationMode (int mode)
 

Private Member Functions

void computeEPcombination ()
 

Private Attributes

double combinedMomentum_
 
double combinedMomentumError_
 
int elClass_
 
SimpleElectron electron_
 
int mode_
 
double scEnergy_
 
double scEnergyError_
 
double trackerMomentum_
 
double trackerMomentumError_
 

Detailed Description

Definition at line 8 of file ElectronEPcombinator.h.

Constructor & Destructor Documentation

ElectronEPcombinator::ElectronEPcombinator ( )
inline

Definition at line 11 of file ElectronEPcombinator.h.

11 {}

Member Function Documentation

void ElectronEPcombinator::combine ( SimpleElectron &  electron)

Definition at line 5 of file ElectronEPcombinator.cc.

References combinedMomentum_, combinedMomentumError_, computeEPcombination(), metsig::electron, electron_, SimpleElectron::setCombinedMomentum(), and SimpleElectron::setCombinedMomentumError().

6 {
7  electron_ = electron;
8  computeEPcombination();
9  electron.setCombinedMomentum(combinedMomentum_);
10  electron.setCombinedMomentumError(combinedMomentumError_);
11 }
SimpleElectron::setCombinedMomentum
void setCombinedMomentum(double combinedMomentum)
Definition: SimpleElectron.h:67
SimpleElectron::setCombinedMomentumError
void setCombinedMomentumError(double combinedMomentumError)
Definition: SimpleElectron.h:68
void ElectronEPcombinator::computeEPcombination ( )
private

Definition at line 14 of file ElectronEPcombinator.cc.

References combinedMomentum_, combinedMomentumError_, elClass_, electron_, SimpleElectron::getElClass(), SimpleElectron::getNewEnergy(), SimpleElectron::getNewEnergyError(), SimpleElectron::getRegEnergy(), SimpleElectron::getRegEnergyError(), SimpleElectron::getTrackerMomentum(), SimpleElectron::getTrackerMomentumError(), mode_, scEnergy_, scEnergyError_, mathSSE::sqrt(), trackerMomentum_, and trackerMomentumError_.

Referenced by combine().

15 {
16  if ( mode_==1 )
17  {
18  scEnergy_ = electron_.getNewEnergy();
19  scEnergyError_ = electron_.getNewEnergyError();
20  }
21  if ( mode_==2 )
22  {
23  scEnergy_ = electron_.getRegEnergy();
24  scEnergyError_ = electron_.getRegEnergyError();
25  }
26  trackerMomentum_ = electron_.getTrackerMomentum();
27  trackerMomentumError_ = electron_.getTrackerMomentumError();
28  elClass_ = electron_.getElClass();
29 
30  combinedMomentum_ = scEnergy_; // initial
31  combinedMomentumError_ = 999.;
32 
33  // first check for large errors
34 
35  if ( trackerMomentumError_/trackerMomentum_ > 0.5
36  && scEnergyError_/scEnergy_ <= 0.5 )
37  {
38  combinedMomentum_ = scEnergy_;
39  combinedMomentumError_ = scEnergyError_;
40  }
41  else if ( trackerMomentumError_/trackerMomentum_ <= 0.5
42  && scEnergyError_/scEnergy_ > 0.5 )
43  {
44  combinedMomentum_ = trackerMomentum_;
45  combinedMomentumError_ = trackerMomentumError_;
46  }
47  else if ( trackerMomentumError_/trackerMomentum_ > 0.5
48  && scEnergyError_/scEnergy_ > 0.5 )
49  {
50  if ( trackerMomentumError_/trackerMomentum_ < scEnergyError_/scEnergy_ )
51  {
52  combinedMomentum_ = trackerMomentum_;
53  combinedMomentumError_ = trackerMomentumError_;
54  }
55  else
56  {
57  combinedMomentum_ = scEnergy_;
58  combinedMomentumError_ = scEnergyError_;
59  }
60  }
61 
62  // then apply the combination algorithm
63  else
64  {
65  // calculate E/p and corresponding error
66  float eOverP = scEnergy_ / trackerMomentum_;
67  float errorEOverP = sqrt(
68  (scEnergyError_/trackerMomentum_)*(scEnergyError_/trackerMomentum_) +
69  (scEnergy_*trackerMomentumError_/trackerMomentum_/trackerMomentum_) *
70  (scEnergy_*trackerMomentumError_/trackerMomentum_/trackerMomentum_) );
71 
72  bool eleIsNotInCombination = false;
73  if ( ( eOverP > 1 + 2.5*errorEOverP )
74  || ( eOverP < 1 - 2.5*errorEOverP )
75  || ( eOverP < 0.8 ) || ( eOverP > 1.3 ) )
76  {
77  eleIsNotInCombination = true;
78  }
79 
80  if ( eleIsNotInCombination )
81  {
82  if ( eOverP > 1 )
83  {
84  combinedMomentum_ = scEnergy_;
85  combinedMomentumError_ = scEnergyError_;
86  }
87  else
88  {
89  if ( elClass_ == 0 ) // == reco::GsfElectron::GOLDEN)
90  {
91  combinedMomentum_ = scEnergy_;
92  combinedMomentumError_ = scEnergyError_;
93  }
94  if ( elClass_ == 1 ) //reco::GsfElectron::BIGBREM)
95  {
96  if ( scEnergy_ < 36 )
97  {
98  combinedMomentum_ = trackerMomentum_;
99  combinedMomentumError_ = trackerMomentumError_;
100  } else
101  {
102  combinedMomentum_ = scEnergy_;
103  combinedMomentumError_ = scEnergyError_;
104  }
105  }
106  if ( elClass_ == 2 ) // == reco::GsfElectron::BADTRACK)
107  {
108  combinedMomentum_ = scEnergy_;
109  combinedMomentumError_ = scEnergyError_;
110  }
111  if ( elClass_ == 3 ) //reco::GsfElectron::SHOWERING)
112  {
113  if ( scEnergy_ < 30 )
114  {
115  combinedMomentum_ = trackerMomentum_;
116  combinedMomentumError_ = trackerMomentumError_;
117  } else
118  {
119  combinedMomentum_ = scEnergy_;
120  combinedMomentumError_ = scEnergyError_;
121  }
122  }
123  if ( elClass_ == 4 ) //reco::GsfElectron::GAP)
124  {
125  if ( scEnergy_ < 60 )
126  {
127  combinedMomentum_ = trackerMomentum_;
128  combinedMomentumError_ = trackerMomentumError_;
129  } else
130  {
131  combinedMomentum_ = scEnergy_;
132  combinedMomentumError_ = scEnergyError_;
133  }
134  }
135  }
136  } else
137  {
138  // combination
139  combinedMomentum_ = (scEnergy_/scEnergyError_/scEnergyError_
140  + trackerMomentum_/trackerMomentumError_/trackerMomentumError_) /
141  (1/scEnergyError_/scEnergyError_ + 1/trackerMomentumError_/trackerMomentumError_);
142  float combinedMomentum_Variance = 1 / (1/scEnergyError_/scEnergyError_ + 1/trackerMomentumError_/trackerMomentumError_);
143  combinedMomentumError_ = sqrt(combinedMomentum_Variance);
144  }
145  }
146 }
SimpleElectron::getRegEnergyError
double getRegEnergyError()
Definition: SimpleElectron.h:54
SimpleElectron::getTrackerMomentum
double getTrackerMomentum()
Definition: SimpleElectron.h:55
SimpleElectron::getNewEnergyError
double getNewEnergyError()
Definition: SimpleElectron.h:46
SimpleElectron::getTrackerMomentumError
double getTrackerMomentumError()
Definition: SimpleElectron.h:56
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
SimpleElectron::getNewEnergy
double getNewEnergy()
Definition: SimpleElectron.h:45
SimpleElectron::getRegEnergy
double getRegEnergy()
Definition: SimpleElectron.h:53
void ElectronEPcombinator::setCombinationMode ( int  mode)
inline

Definition at line 13 of file ElectronEPcombinator.h.

References alignBH_cfg::mode, and mode_.

13 {mode_ = mode;}
alignBH_cfg.mode
list mode
Definition: alignBH_cfg.py:7

Member Data Documentation

double ElectronEPcombinator::combinedMomentum_
private

Definition at line 18 of file ElectronEPcombinator.h.

Referenced by combine(), and computeEPcombination().

double ElectronEPcombinator::combinedMomentumError_
private

Definition at line 19 of file ElectronEPcombinator.h.

Referenced by combine(), and computeEPcombination().

int ElectronEPcombinator::elClass_
private

Definition at line 24 of file ElectronEPcombinator.h.

Referenced by computeEPcombination().

SimpleElectron ElectronEPcombinator::electron_
private

Definition at line 16 of file ElectronEPcombinator.h.

Referenced by combine(), and computeEPcombination().

int ElectronEPcombinator::mode_
private

Definition at line 25 of file ElectronEPcombinator.h.

Referenced by computeEPcombination(), and setCombinationMode().

double ElectronEPcombinator::scEnergy_
private

Definition at line 20 of file ElectronEPcombinator.h.

Referenced by computeEPcombination().

double ElectronEPcombinator::scEnergyError_
private

Definition at line 21 of file ElectronEPcombinator.h.

Referenced by computeEPcombination().

double ElectronEPcombinator::trackerMomentum_
private

Definition at line 22 of file ElectronEPcombinator.h.

Referenced by computeEPcombination().

double ElectronEPcombinator::trackerMomentumError_
private

Definition at line 23 of file ElectronEPcombinator.h.

Referenced by computeEPcombination().

Generated for CMSSW Reference Manual by doxygen 1.8.5