Tools/src/TrackDetMatchInfo.cc

0001 #include <algorithm>
0002 #include "Calibration/Tools/interface/TrackDetMatchInfo.h"
0003
0004 double HTrackDetMatchInfo::ecalEnergyFromRecHits() {
0005   double energy(0);
0006   for (std::vector<EcalRecHit>::const_iterator hit = crossedEcalRecHits.begin(); hit != crossedEcalRecHits.end(); hit++)
0007     energy += hit->energy();
0008   return energy;
0009 }
0010
0011 double HTrackDetMatchInfo::ecalConeEnergyFromRecHits() {
0012   double energy(0);
0013   for (std::vector<EcalRecHit>::const_iterator hit = coneEcalRecHits.begin(); hit != coneEcalRecHits.end(); hit++) {
0014     energy += hit->energy();
0015     //     std::cout<< hit->detid().rawId()<<" "<<hit->energy()<<" "<<energy<<std::endl;
0016   }
0017   return energy;
0018 }
0019
0020 double HTrackDetMatchInfo::ecalEnergyFromCaloTowers() {
0021   double energy(0);
0022   for (std::vector<CaloTower>::const_iterator hit = crossedTowers.begin(); hit != crossedTowers.end(); hit++) {
0023     energy += hit->emEnergy();
0024   }
0025   return energy;
0026 }
0027
0028 double HTrackDetMatchInfo::ecalConeEnergyFromCaloTowers() {
0029   double energy(0);
0030   for (std::vector<CaloTower>::const_iterator hit = coneTowers.begin(); hit != coneTowers.end(); hit++)
0031     energy += hit->emEnergy();
0032   return energy;
0033 }
0034
0035 double HTrackDetMatchInfo::hcalEnergyFromRecHits() {
0036   double energy(0);
0037   for (std::vector<HBHERecHit>::const_iterator hit = crossedHcalRecHits.begin(); hit != crossedHcalRecHits.end(); hit++)
0038     energy += hit->energy();
0039   return energy;
0040 }
0041
0042 double HTrackDetMatchInfo::hcalConeEnergyFromRecHits() {
0043   double energy(0);
0044   for (std::vector<HBHERecHit>::const_iterator hit = coneHcalRecHits.begin(); hit != coneHcalRecHits.end(); hit++) {
0045     energy += hit->energy();
0046   }
0047   return energy;
0048 }
0049
0050 double HTrackDetMatchInfo::hcalBoxEnergyFromRecHits() {
0051   double energy(0);
0052   for (std::vector<HBHERecHit>::const_iterator hit = boxHcalRecHits.begin(); hit != boxHcalRecHits.end(); hit++)
0053     energy += hit->energy();
0054   return energy;
0055 }
0056
0057 double HTrackDetMatchInfo::hcalEnergyFromCaloTowers() {
0058   double energy(0);
0059   for (std::vector<CaloTower>::const_iterator tower = crossedTowers.begin(); tower != crossedTowers.end(); tower++)
0060     energy += tower->hadEnergy();
0061   return energy;
0062 }
0063
0064 double HTrackDetMatchInfo::hcalConeEnergyFromCaloTowers() {
0065   double energy(0);
0066   for (std::vector<CaloTower>::const_iterator hit = coneTowers.begin(); hit != coneTowers.end(); hit++) {
0067     energy += hit->hadEnergy();
0068   }
0069   return energy;
0070 }
0071
0072 double HTrackDetMatchInfo::hcalBoxEnergyFromCaloTowers() {
0073   double energy(0);
0074   for (std::vector<CaloTower>::const_iterator hit = boxTowers.begin(); hit != boxTowers.end(); hit++)
0075     energy += hit->hadEnergy();
0076   return energy;
0077 }
0078
0079 double HTrackDetMatchInfo::outerHcalEnergy() {
0080   double energy(0);
0081   for (std::vector<CaloTower>::const_iterator tower = crossedTowers.begin(); tower != crossedTowers.end(); tower++)
0082     energy += tower->outerEnergy();
0083   return energy;
0084 }