|
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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
|
#include <algorithm>
#include "Calibration/Tools/interface/TrackDetMatchInfo.h"
double HTrackDetMatchInfo::ecalEnergyFromRecHits() {
double energy(0);
for (std::vector<EcalRecHit>::const_iterator hit = crossedEcalRecHits.begin(); hit != crossedEcalRecHits.end(); hit++)
energy += hit->energy();
return energy;
}
double HTrackDetMatchInfo::ecalConeEnergyFromRecHits() {
double energy(0);
for (std::vector<EcalRecHit>::const_iterator hit = coneEcalRecHits.begin(); hit != coneEcalRecHits.end(); hit++) {
energy += hit->energy();
// std::cout<< hit->detid().rawId()<<" "<<hit->energy()<<" "<<energy<<std::endl;
}
return energy;
}
double HTrackDetMatchInfo::ecalEnergyFromCaloTowers() {
double energy(0);
for (std::vector<CaloTower>::const_iterator hit = crossedTowers.begin(); hit != crossedTowers.end(); hit++) {
energy += hit->emEnergy();
}
return energy;
}
double HTrackDetMatchInfo::ecalConeEnergyFromCaloTowers() {
double energy(0);
for (std::vector<CaloTower>::const_iterator hit = coneTowers.begin(); hit != coneTowers.end(); hit++)
energy += hit->emEnergy();
return energy;
}
double HTrackDetMatchInfo::hcalEnergyFromRecHits() {
double energy(0);
for (std::vector<HBHERecHit>::const_iterator hit = crossedHcalRecHits.begin(); hit != crossedHcalRecHits.end(); hit++)
energy += hit->energy();
return energy;
}
double HTrackDetMatchInfo::hcalConeEnergyFromRecHits() {
double energy(0);
for (std::vector<HBHERecHit>::const_iterator hit = coneHcalRecHits.begin(); hit != coneHcalRecHits.end(); hit++) {
energy += hit->energy();
}
return energy;
}
double HTrackDetMatchInfo::hcalBoxEnergyFromRecHits() {
double energy(0);
for (std::vector<HBHERecHit>::const_iterator hit = boxHcalRecHits.begin(); hit != boxHcalRecHits.end(); hit++)
energy += hit->energy();
return energy;
}
double HTrackDetMatchInfo::hcalEnergyFromCaloTowers() {
double energy(0);
for (std::vector<CaloTower>::const_iterator tower = crossedTowers.begin(); tower != crossedTowers.end(); tower++)
energy += tower->hadEnergy();
return energy;
}
double HTrackDetMatchInfo::hcalConeEnergyFromCaloTowers() {
double energy(0);
for (std::vector<CaloTower>::const_iterator hit = coneTowers.begin(); hit != coneTowers.end(); hit++) {
energy += hit->hadEnergy();
}
return energy;
}
double HTrackDetMatchInfo::hcalBoxEnergyFromCaloTowers() {
double energy(0);
for (std::vector<CaloTower>::const_iterator hit = boxTowers.begin(); hit != boxTowers.end(); hit++)
energy += hit->hadEnergy();
return energy;
}
double HTrackDetMatchInfo::outerHcalEnergy() {
double energy(0);
for (std::vector<CaloTower>::const_iterator tower = crossedTowers.begin(); tower != crossedTowers.end(); tower++)
energy += tower->outerEnergy();
return energy;
}
|