|
|
|
#include <TrackAssociator.h>
Definition at line 47 of file TrackAssociator.h.
| HTrackAssociator::HTrackAssociator | ( | ) |
Definition at line 54 of file TrackAssociator.cc.
{
ivProp_ = 0;
defProp_ = 0;
debug_ = 0;
caloTowerMap_ = 0;
useDefaultPropagator_ = false;
}
| HTrackAssociator::~HTrackAssociator | ( | ) |
Definition at line 63 of file TrackAssociator.cc.
| void HTrackAssociator::addDataLabels | ( | const std::string | className, |
| const std::string | moduleLabel, | ||
| const std::string | productInstanceLabel = "" |
||
| ) |
specify names of EDProducts to use for different input data types
Definition at line 69 of file TrackAssociator.cc.
{
if (className == "EBRecHitCollection")
{
EBRecHitCollectionLabels.clear();
EBRecHitCollectionLabels.push_back(moduleLabel);
EBRecHitCollectionLabels.push_back(productInstanceLabel);
}
if (className == "EERecHitCollection")
{
EERecHitCollectionLabels.clear();
EERecHitCollectionLabels.push_back(moduleLabel);
EERecHitCollectionLabels.push_back(productInstanceLabel);
}
if (className == "HBHERecHitCollection")
{
HBHERecHitCollectionLabels.clear();
HBHERecHitCollectionLabels.push_back(moduleLabel);
HBHERecHitCollectionLabels.push_back(productInstanceLabel);
}
if (className == "CaloTowerCollection")
{
CaloTowerCollectionLabels.clear();
CaloTowerCollectionLabels.push_back(moduleLabel);
CaloTowerCollectionLabels.push_back(productInstanceLabel);
}
}
| HTrackDetMatchInfo HTrackAssociator::associate | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| const FreeTrajectoryState & | trackOrigin, | ||
| const HAssociatorParameters & | parameters | ||
| ) |
propagate a track across the whole detector and find associated objects. Association is done in two modes 1) an object is associated to a track only if crossed by track 2) an object is associated to a track if it is withing an eta-phi cone of some radius with respect to a track. (the cone origin is at (0,0,0))
Definition at line 141 of file TrackAssociator.cc.
References HTrackAssociator::HAssociatorParameters::dRCalo, HTrackAssociator::HAssociatorParameters::dREcal, HTrackAssociator::HAssociatorParameters::dRHcal, HTrackAssociator::HAssociatorParameters::idRCalo, HTrackAssociator::HAssociatorParameters::idREcal, HTrackAssociator::HAssociatorParameters::idRHcal, info, init, HTrackAssociator::HAssociatorParameters::useCalo, HTrackAssociator::HAssociatorParameters::useEcal, and HTrackAssociator::HAssociatorParameters::useHcal.
{
HTrackDetMatchInfo info;
using namespace edm;
HTimerStack timers;
init( iSetup );
FreeTrajectoryState currentPosition(trackOrigin);
if (parameters.useEcal) fillEcal( iEvent, iSetup, info, currentPosition,parameters.idREcal, parameters.dREcal);
if (parameters.useHcal) fillHcal( iEvent, iSetup, info, currentPosition,parameters.idRHcal,parameters.dRHcal);
if (parameters.useCalo) fillCaloTowers( iEvent, iSetup, info, currentPosition,parameters.idRCalo,parameters.dRCalo);
return info;
}
| std::vector< EcalRecHit > HTrackAssociator::associateEcal | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| const FreeTrajectoryState & | trackOrigin, | ||
| const double | dR = -1 |
||
| ) |
associate ECAL only and return RecHits negative dR means only crossed elements
Definition at line 162 of file TrackAssociator.cc.
References HTrackDetMatchInfo::coneEcalRecHits, HTrackDetMatchInfo::crossedEcalRecHits, HTrackAssociator::HAssociatorParameters::dREcal, info, Parameters::parameters, and HTrackAssociator::HAssociatorParameters::useHcal.
{
HAssociatorParameters parameters;
parameters.useHcal = false;
parameters.dREcal = dR;
HTrackDetMatchInfo info( associate(iEvent, iSetup, trackOrigin, parameters ));
if (dR>0)
return info.coneEcalRecHits;
else
return info.crossedEcalRecHits;
}
| std::vector< CaloTower > HTrackAssociator::associateHcal | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| const FreeTrajectoryState & | trackOrigin, | ||
| const double | dR = -1 |
||
| ) |
associate ECAL only and return RecHits negative dR means only crossed elements
Definition at line 194 of file TrackAssociator.cc.
References HTrackDetMatchInfo::coneTowers, HTrackDetMatchInfo::crossedTowers, HTrackAssociator::HAssociatorParameters::dRHcal, info, Parameters::parameters, and HTrackAssociator::HAssociatorParameters::useEcal.
{
HAssociatorParameters parameters;
parameters.useEcal = false;
parameters.dRHcal = dR;
HTrackDetMatchInfo info( associate(iEvent, iSetup, trackOrigin, parameters ));
if (dR>0)
return info.coneTowers;
else
return info.crossedTowers;
}
| void HTrackAssociator::fillCaloTowers | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| HTrackDetMatchInfo & | info, | ||
| const FreeTrajectoryState & | trajectoryPoint, | ||
| const int | idR, | ||
| const double | dR | ||
| ) | [private] |
Definition at line 312 of file TrackAssociator.cc.
References HTrackDetMatchInfo::boxTowers, HiRecoJets_cff::caloTowers, HTrackDetMatchInfo::coneTowers, gather_cfg::cout, HTrackDetMatchInfo::crossedTowers, edm::Event::getByLabel(), errorMatrix2Lands_multiChannel::id, HTrackDetMatchInfo::isGoodCalo, HTrackDetMatchInfo::isGoodEcal, HTrackDetMatchInfo::isGoodHcal, edm::HandleBase::isValid(), LogTrace, HTimerStack::pop_and_push(), HTimerStack::push(), HTrackDetMatchInfo::regionTowers, HTrackDetMatchInfo::trkGlobPosAtEcal, and HTrackDetMatchInfo::trkGlobPosAtHcal.
{
// ECAL hits are not used for the CaloTower identification
HTimerStack timers;
timers.push("HTrackAssociator::fillCaloTowers");
caloDetIdAssociator_.setGeometry(&*theCaloGeometry_);
// HCAL points (HB+HE)
timers.push("HTrackAssociator::fillCaloTowers::propagation");
std::vector<GlobalPoint> hcalPoints;
hcalPoints.push_back(GlobalPoint(135.,0,310.));
hcalPoints.push_back(GlobalPoint(150.,0,340.));
hcalPoints.push_back(GlobalPoint(170.,0,370.));
hcalPoints.push_back(GlobalPoint(190.,0,400.));
hcalPoints.push_back(GlobalPoint(240.,0,500.));
hcalPoints.push_back(GlobalPoint(280.,0,550.));
std::vector<GlobalPoint> hcalTrajectory = caloDetIdAssociator_.getTrajectory(trajectoryPoint, hcalPoints);
// if(hcalTrajectory.empty()) throw cms::Exception("FatalError") << "Failed to propagate the track to HCAL\n";
if(hcalTrajectory.empty()) {
LogTrace("HTrackAssociator::fillEcal") << "Failed to propagate a track to ECAL; moving on\n";
info.isGoodCalo = 0;
info.isGoodEcal = 0;
std::cout<<" HTrackAssociator::fillCaloTowers::Failed to propagate a track to ECAL "<<std::endl;
return;
}
info.isGoodCalo = 1;
info.isGoodEcal = 1;
info.trkGlobPosAtEcal = getPoint(hcalTrajectory[0]);
if(hcalTrajectory.size()<4) {
LogTrace("HTrackAssociator::fillEcal") << "Failed to propagate a track to HCAL; moving on\n";
info.isGoodHcal = 0;
}
info.isGoodHcal = 1;
info.trkGlobPosAtHcal = getPoint(hcalTrajectory[4]);
// find crossed CaloTowers
timers.pop_and_push("HTrackAssociator::fillCaloTowers::access::CaloTowers");
edm::Handle<CaloTowerCollection> caloTowers;
if (CaloTowerCollectionLabels.empty())
throw cms::Exception("FatalError") << "Module lable is not set for CaloTowers.\n";
else
iEvent.getByLabel (CaloTowerCollectionLabels[0], CaloTowerCollectionLabels[1], caloTowers);
if (!caloTowers.isValid()) throw cms::Exception("FatalError") << "Unable to find CaloTowers in event!\n";
timers.push("HTrackAssociator::fillCaloTowers::matching");
// first get DetIds in a predefined NxN region
// std::set<DetId> caloTowerIdsInBigRegion = caloDetIdAssociator_.getDetIdsCloseToAPoint(hcalTrajectory[0],idR+1);
std::set<DetId> caloTowerIdsInRegion = caloDetIdAssociator_.getDetIdsCloseToAPoint(hcalTrajectory[0],idR);
std::set<DetId> caloTowerIdsInACone;
std::set<DetId> crossedCaloTowerIds;
std::set<DetId> caloTowerIdsInBox;
caloTowerIdsInACone = caloDetIdAssociator_.getDetIdsInACone(caloTowerIdsInRegion, hcalTrajectory, dR);
// get DetId of the most energetic tower in that region
crossedCaloTowerIds = caloDetIdAssociator_.getMaxEDetId(caloTowerIdsInRegion, caloTowers);
// get DetIds of the towers surrounding the most energetic one
caloTowerIdsInBox = caloDetIdAssociator_.getDetIdsInACone(crossedCaloTowerIds, hcalTrajectory, -1.);
//
// Debug prints
//
// std::cout <<" Debug printout in CaloTowers "<<std::endl;
// std::cout <<" with position at outer layer:r,z,phi "<<trajectoryPoint.position().eta()<<
// " "<<trajectoryPoint.position().phi()<<
// " "<<trajectoryPoint.position().perp()<<
// " "<<trajectoryPoint.position().z()<<
// " "<<trajectoryPoint.charge()<<std::endl;
// std::cout <<" Trajectory point at ECAL surface:eta:phi:radius:z "<<(hcalTrajectory[0]).eta()<<
// " "<<(hcalTrajectory[0]).phi()<<
// " "<<(hcalTrajectory[0]).perp()<<
// " "<<(hcalTrajectory[0]).z()<<
// " momentum "<<trajectoryPoint.momentum().perp()<<std::endl;
//
// std::cout<<" Number of towers in the region "<<caloTowerIdsInRegion.size()<<" idR= "<<idR<<std::endl;
// add CaloTowers
timers.push("HTrackAssociator::fillCaloTowers::addCaloTowers");
for(std::set<DetId>::const_iterator itr=crossedCaloTowerIds.begin(); itr!=crossedCaloTowerIds.end();itr++)
{
DetId id(*itr);
CaloTowerCollection::const_iterator tower = (*caloTowers).find(id);
if(tower != (*caloTowers).end())
info.crossedTowers.push_back(*tower);
else
LogTrace("HTrackAssociator::fillEcal") << "CaloTower is not found for DetId: " << id.rawId() << "\n";
}
for(std::set<DetId>::const_iterator itr=caloTowerIdsInACone.begin(); itr!=caloTowerIdsInACone.end();itr++)
{
DetId id(*itr);
CaloTowerCollection::const_iterator tower = (*caloTowers).find(id);
if(tower != (*caloTowers).end()) {
info.coneTowers.push_back(*tower);
}
else
LogTrace("HTrackAssociator::fillEcal") << "CaloTower is not found for DetId: " << id.rawId() << "\n";
}
for(std::set<DetId>::const_iterator itr=caloTowerIdsInBox.begin(); itr!=caloTowerIdsInBox.end();itr++)
{
DetId id(*itr);
CaloTowerCollection::const_iterator tower = (*caloTowers).find(id);
if(tower != (*caloTowers).end()) {
info.boxTowers.push_back(*tower);
}
else
LogTrace("HTrackAssociator::fillEcal") << "CaloTower is not found for DetId: " << id.rawId() << "\n";
}
for(std::set<DetId>::const_iterator itr=caloTowerIdsInRegion.begin(); itr!=caloTowerIdsInRegion.end();itr++)
{
DetId id(*itr);
CaloTowerCollection::const_iterator tower = (*caloTowers).find(id);
if(tower != (*caloTowers).end()) {
info.regionTowers.push_back(*tower);
}
else
LogTrace("HTrackAssociator::fillEcal") << "CaloTower is not found for DetId: " << id.rawId() << "\n";
}
}
| void HTrackAssociator::fillEcal | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| HTrackDetMatchInfo & | info, | ||
| const FreeTrajectoryState & | trajectoryPoint, | ||
| const int | idR, | ||
| const double | dR | ||
| ) | [private] |
Definition at line 227 of file TrackAssociator.cc.
References HTrackDetMatchInfo::coneEcalRecHits, HTrackDetMatchInfo::crossedEcalRecHits, egHLT::errCodes::EBRecHits, egHLT::errCodes::EERecHits, edm::Event::getByLabel(), HTrackDetMatchInfo::isGoodEcal, LogTrace, HTimerStack::pop_and_push(), HTimerStack::push(), and HTrackDetMatchInfo::trkGlobPosAtEcal.
{
HTimerStack timers;
timers.push("HTrackAssociator::fillEcal");
ecalDetIdAssociator_.setGeometry(&*theCaloGeometry_);
timers.push("HTrackAssociator::fillEcal::propagation");
// ECAL points (EB+EE)
std::vector<GlobalPoint> ecalPoints;
ecalPoints.push_back(GlobalPoint(135.,0,310.));
ecalPoints.push_back(GlobalPoint(150.,0,340.));
ecalPoints.push_back(GlobalPoint(170.,0,370.));
std::vector<GlobalPoint> ecalTrajectory = ecalDetIdAssociator_.getTrajectory(trajectoryPoint, ecalPoints);
// if(ecalTrajectory.empty()) throw cms::Exception("FatalError") << "Failed to propagate a track to ECAL\n";
if(ecalTrajectory.empty()) {
LogTrace("HTrackAssociator::fillEcal") << "Failed to propagate a track to ECAL; moving on\n";
info.isGoodEcal = 0;
return;
}
info.isGoodEcal = 1;
info.trkGlobPosAtEcal = getPoint(ecalTrajectory[0]);
// Find ECAL crystals
timers.pop_and_push("HTrackAssociator::fillEcal::access::EcalBarrel");
edm::Handle<EBRecHitCollection> EBRecHits;
edm::Handle<EERecHitCollection> EERecHits;
// if (EBRecHitCollectionLabels.empty())
// throw cms::Exception("FatalError") << "Module lable is not set for EBRecHitCollection.\n";
// else
iEvent.getByLabel (EBRecHitCollectionLabels[0], EBRecHitCollectionLabels[1], EBRecHits);
// if (!EBRecHits.isValid()) throw cms::Exception("FatalError") << "Unable to find EBRecHitCollection in event!\n";
if (EERecHitCollectionLabels[1]!=EBRecHitCollectionLabels[1]) iEvent.getByLabel (EERecHitCollectionLabels[0], EERecHitCollectionLabels[1], EERecHits);
// if (!EERecHits.isValid()) throw cms::Exception("FatalError") << "Unable to find EERecHitCollection in event!\n";
timers.pop_and_push("HTrackAssociator::fillEcal::matching");
std::set<DetId> ecalIdsInRegion = ecalDetIdAssociator_.getDetIdsCloseToAPoint(ecalTrajectory[0],idR);
std::set<DetId> ecalIdsInACone = ecalDetIdAssociator_.getDetIdsInACone(ecalIdsInRegion, ecalTrajectory, dR);
std::set<DetId> crossedEcalIds = ecalDetIdAssociator_.getCrossedDetIds(ecalIdsInRegion, ecalTrajectory);
// add EcalRecHits
timers.pop_and_push("HTrackAssociator::fillEcal::addEcalRecHits");
for(std::set<DetId>::const_iterator itr=crossedEcalIds.begin(); itr!=crossedEcalIds.end();itr++) {
std::vector<EcalRecHit>::const_iterator hit = (*EBRecHits).find(*itr);
if(hit != (*EBRecHits).end())
info.crossedEcalRecHits.push_back(*hit);
else
LogTrace("HTrackAssociator::fillEcal") << "EcalRecHit is not found for DetId: " << itr->rawId() <<"\n";
}
for(std::set<DetId>::const_iterator itr=ecalIdsInACone.begin(); itr!=ecalIdsInACone.end();itr++) {
std::vector<EcalRecHit>::const_iterator hit = (*EBRecHits).find(*itr);
if(hit != (*EBRecHits).end()) {
info.coneEcalRecHits.push_back(*hit);
}
else
LogTrace("HTrackAssociator::fillEcal") << "EcalRecHit is not found for DetId: " << itr->rawId() <<"\n";
}
if (EERecHitCollectionLabels[1]==EBRecHitCollectionLabels[1])return;
for(std::set<DetId>::const_iterator itr=crossedEcalIds.begin(); itr!=crossedEcalIds.end();itr++) {
std::vector<EcalRecHit>::const_iterator hit = (*EERecHits).find(*itr);
if(hit != (*EERecHits).end())
info.crossedEcalRecHits.push_back(*hit);
else
LogTrace("HTrackAssociator::fillEcal") << "EcalRecHit is not found for DetId: " << itr->rawId() <<"\n";
}
for(std::set<DetId>::const_iterator itr=ecalIdsInACone.begin(); itr!=ecalIdsInACone.end();itr++) {
std::vector<EcalRecHit>::const_iterator hit = (*EERecHits).find(*itr);
if(hit != (*EERecHits).end()) {
info.coneEcalRecHits.push_back(*hit);
}
else
LogTrace("HTrackAssociator::fillEcal") << "EcalRecHit is not found for DetId: " << itr->rawId() <<"\n";
}
}
| void HTrackAssociator::fillHcal | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| HTrackDetMatchInfo & | info, | ||
| const FreeTrajectoryState & | trajectoryPoint, | ||
| const int | idR, | ||
| const double | dR | ||
| ) | [private] |
Definition at line 449 of file TrackAssociator.cc.
References HTrackDetMatchInfo::boxHcalRecHits, HTrackDetMatchInfo::coneHcalRecHits, HTrackDetMatchInfo::crossedHcalRecHits, edm::Event::getByLabel(), egHLT::errCodes::HBHERecHits, HTrackDetMatchInfo::isGoodHcal, edm::HandleBase::isValid(), LogTrace, HTimerStack::pop_and_push(), HTimerStack::push(), HTrackDetMatchInfo::regionHcalRecHits, and HTrackDetMatchInfo::trkGlobPosAtHcal.
{
HTimerStack timers;
timers.push("HTrackAssociator::fillHcal");
hcalDetIdAssociator_.setGeometry(&*theCaloGeometry_);
// HCAL points (HB+HE)
timers.push("HTrackAssociator::fillHcal::propagation");
std::vector<GlobalPoint> hcalPoints;
hcalPoints.push_back(GlobalPoint(190.,0,400.));
hcalPoints.push_back(GlobalPoint(240.,0,500.));
hcalPoints.push_back(GlobalPoint(280.,0,550.));
std::vector<GlobalPoint> hcalTrajectory = hcalDetIdAssociator_.getTrajectory(trajectoryPoint, hcalPoints);
// if(hcalTrajectory.empty()) throw cms::Exception("FatalError") << "Failed to propagate the track to HCAL\n";
if(hcalTrajectory.empty()) {
LogTrace("HTrackAssociator::fillHcal") << "Failed to propagate a track to HCAL; moving on\n";
info.isGoodHcal = 0;
// std::cout<<" HTrackAssociator::fillHcal::Failed to propagate a track to HCAL "<<std::endl;
return;
}
info.isGoodHcal = 1;
info.trkGlobPosAtHcal = getPoint(hcalTrajectory[0]);
timers.pop_and_push("HTrackAssociator::fillHcal::access::Hcal");
edm::Handle<HBHERecHitCollection> HBHERecHits;
// if (HBHERecHitCollectionLabels.empty())
// throw cms::Exception("FatalError") << "Module label is not set for HBHERecHitCollection.\n";
// else
iEvent.getByLabel (HBHERecHitCollectionLabels[0], HBHERecHitCollectionLabels[1], HBHERecHits);
if (!HBHERecHits.isValid()) throw cms::Exception("FatalError") << "Unable to find HBHERecHitCollection in event!\n";
timers.pop_and_push("HTrackAssociator::fillHcal::matching");
// first get DetIds in a predefined NxN region
// std::set<DetId> hcalIdsInBigRegion = hcalDetIdAssociator_.getDetIdsCloseToAPoint(hcalTrajectory[0],idR+1);
std::set<DetId> hcalIdsInRegion = hcalDetIdAssociator_.getDetIdsCloseToAPoint(hcalTrajectory[0],idR);
std::set<DetId> hcalIdsInACone;
std::set<DetId> crossedHcalIds;
std::set<DetId> hcalIdsInBox;
hcalIdsInACone = hcalDetIdAssociator_.getDetIdsInACone(hcalIdsInRegion, hcalTrajectory, dR);
// get DetId of the most energetic tower in that region
crossedHcalIds = hcalDetIdAssociator_.getMaxEDetId(hcalIdsInRegion, HBHERecHits);
// get DetIds of the towers surrounding the most energetic one
hcalIdsInBox = hcalDetIdAssociator_.getDetIdsInACone(crossedHcalIds, hcalTrajectory, -1.);
// add HcalRecHits
timers.pop_and_push("HTrackAssociator::fillHcal::addHcalRecHits");
for(std::set<DetId>::const_iterator itr=crossedHcalIds.begin(); itr!=crossedHcalIds.end();itr++) {
std::vector<HBHERecHit>::const_iterator hit = (*HBHERecHits).find(*itr);
if(hit != (*HBHERecHits).end())
info.crossedHcalRecHits.push_back(*hit);
else
LogTrace("HTrackAssociator::fillHcal") << "HcalRecHit is not found for DetId: " << itr->rawId() <<"\n";
}
for(std::set<DetId>::const_iterator itr=hcalIdsInACone.begin(); itr!=hcalIdsInACone.end();itr++) {
std::vector<HBHERecHit>::const_iterator hit = (*HBHERecHits).find(*itr);
if(hit != (*HBHERecHits).end())
info.coneHcalRecHits.push_back(*hit);
else
LogTrace("HTrackAssociator::fillHcal") << "HcalRecHit is not found for DetId: " << itr->rawId() <<"\n";
}
for(std::set<DetId>::const_iterator itr=hcalIdsInBox.begin(); itr!=hcalIdsInBox.end();itr++) {
std::vector<HBHERecHit>::const_iterator hit = (*HBHERecHits).find(*itr);
if(hit != (*HBHERecHits).end())
info.boxHcalRecHits.push_back(*hit);
else
LogTrace("HTrackAssociator::fillHcal") << "HcalRecHit is not found for DetId: " << itr->rawId() <<"\n";
}
for(std::set<DetId>::const_iterator itr=hcalIdsInRegion.begin(); itr!=hcalIdsInRegion.end();itr++) {
std::vector<HBHERecHit>::const_iterator hit = (*HBHERecHits).find(*itr);
if(hit != (*HBHERecHits).end())
info.regionHcalRecHits.push_back(*hit);
else
LogTrace("HTrackAssociator::fillHcal") << "HcalRecHit is not found for DetId: " << itr->rawId() <<"\n";
}
}
| void HTrackAssociator::fillHcalTowers | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| HTrackDetMatchInfo & | info, | ||
| const FreeTrajectoryState & | trajectoryPoint, | ||
| const int | idR, | ||
| const double | dR | ||
| ) | [private] |
Definition at line 537 of file TrackAssociator.cc.
References HiRecoJets_cff::caloTowers, HTrackDetMatchInfo::coneTowers, gather_cfg::cout, HTrackDetMatchInfo::crossedTowers, edm::Event::getByLabel(), errorMatrix2Lands_multiChannel::id, HTrackDetMatchInfo::isGoodCalo, edm::HandleBase::isValid(), LogTrace, HTimerStack::pop_and_push(), HTimerStack::push(), and HTrackDetMatchInfo::trkGlobPosAtHcal.
{
// ECAL hits are not used for the CaloTower identification
HTimerStack timers;
timers.push("HTrackAssociator::fillCaloTowers");
caloDetIdAssociator_.setGeometry(&*theCaloGeometry_);
// HCAL points (HB+HE)
timers.push("HTrackAssociator::fillCaloTowers::propagation");
std::vector<GlobalPoint> hcalPoints;
hcalPoints.push_back(GlobalPoint(190.,0,400.));
hcalPoints.push_back(GlobalPoint(240.,0,500.));
hcalPoints.push_back(GlobalPoint(280.,0,550.));
std::vector<GlobalPoint> hcalTrajectory = caloDetIdAssociator_.getTrajectory(trajectoryPoint, hcalPoints);
// if(hcalTrajectory.empty()) throw cms::Exception("FatalError") << "Failed to propagate the track to HCAL\n";
if(hcalTrajectory.empty()) {
LogTrace("HTrackAssociator::fillEcal") << "Failed to propagate a track to HCAL; moving on\n";
info.isGoodCalo = 0;
std::cout<<" HTrackAssociator::fillCaloTowers::Failed to propagate a track to HCAL "<<std::endl;
return;
}
info.isGoodCalo = 1;
info.trkGlobPosAtHcal = getPoint(hcalTrajectory[0]);
// find crossed CaloTowers
timers.pop_and_push("HTrackAssociator::fillCaloTowers::access::CaloTowers");
edm::Handle<CaloTowerCollection> caloTowers;
if (CaloTowerCollectionLabels.empty())
throw cms::Exception("FatalError") << "Module lable is not set for CaloTowers.\n";
else
iEvent.getByLabel (CaloTowerCollectionLabels[0], CaloTowerCollectionLabels[1], caloTowers);
if (!caloTowers.isValid()) throw cms::Exception("FatalError") << "Unable to find CaloTowers in event!\n";
timers.push("HTrackAssociator::fillCaloTowers::matching");
// first get DetIds in a predefined NxN region
std::set<DetId> caloTowerIdsInBigRegion = caloDetIdAssociator_.getDetIdsCloseToAPoint(hcalTrajectory[0],idR+1);
std::set<DetId> caloTowerIdsInRegion = caloDetIdAssociator_.getDetIdsCloseToAPoint(hcalTrajectory[0],idR);
std::set<DetId> caloTowerIdsInACone;
std::set<DetId> crossedCaloTowerIds;
std::set<DetId> caloTowerIdsInBox;
caloTowerIdsInACone = caloDetIdAssociator_.getDetIdsInACone(caloTowerIdsInBigRegion, hcalTrajectory, dR);
// get DetId of the most energetic tower in that region
crossedCaloTowerIds = caloDetIdAssociator_.getMaxEDetId(caloTowerIdsInRegion, caloTowers);
// get DetIds of the towers surrounding the most energetic one
caloTowerIdsInBox = caloDetIdAssociator_.getDetIdsInACone(crossedCaloTowerIds, hcalTrajectory, -1.);
// add CaloTowers
timers.push("HTrackAssociator::fillCaloTowers::addCaloTowers");
for(std::set<DetId>::const_iterator itr=crossedCaloTowerIds.begin(); itr!=crossedCaloTowerIds.end();itr++)
{
DetId id(*itr);
CaloTowerCollection::const_iterator tower = (*caloTowers).find(id);
if(tower != (*caloTowers).end())
info.crossedTowers.push_back(*tower);
else
LogTrace("HTrackAssociator::fillEcal") << "CaloTower is not found for DetId: " << id.rawId() << "\n";
}
for(std::set<DetId>::const_iterator itr=caloTowerIdsInACone.begin(); itr!=caloTowerIdsInACone.end();itr++)
{
DetId id(*itr);
CaloTowerCollection::const_iterator tower = (*caloTowers).find(id);
if(tower != (*caloTowers).end())
info.coneTowers.push_back(*tower);
else
LogTrace("HTrackAssociator::fillEcal") << "CaloTower is not found for DetId: " << id.rawId() << "\n";
}
}
| double HTrackAssociator::getEcalEnergy | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| const FreeTrajectoryState & | trackOrigin, | ||
| const double | dR = -1 |
||
| ) |
associate ECAL only and return energy negative dR means only crossed elements
Definition at line 178 of file TrackAssociator.cc.
References HTrackAssociator::HAssociatorParameters::dREcal, HTrackDetMatchInfo::ecalConeEnergyFromRecHits(), HTrackDetMatchInfo::ecalEnergyFromRecHits(), info, Parameters::parameters, and HTrackAssociator::HAssociatorParameters::useHcal.
{
HAssociatorParameters parameters;
parameters.useHcal = false;
parameters.dREcal = dR;
HTrackDetMatchInfo info = associate(iEvent, iSetup, trackOrigin, parameters );
if(dR>0)
return info.ecalConeEnergyFromRecHits();
else
return info.ecalEnergyFromRecHits();
}
| FreeTrajectoryState HTrackAssociator::getFreeTrajectoryState | ( | const edm::EventSetup & | iSetup, |
| const reco::Track & | track | ||
| ) |
get FreeTrajectoryState from different track representations
Definition at line 641 of file TrackAssociator.cc.
References ecalTB2006H4_GenSimDigiReco_cfg::bField, reco::TrackBase::charge(), alignCSCRings::e, edm::EventSetup::get(), reco::TrackBase::momentum(), point, and reco::TrackBase::vertex().
{
edm::ESHandle<MagneticField> bField;
iSetup.get<IdealMagneticFieldRecord>().get(bField);
GlobalVector vector( track.momentum().x(), track.momentum().y(), track.momentum().z() );
GlobalPoint point( track.vertex().x(), track.vertex().y(), track.vertex().z() );
GlobalTrajectoryParameters tPars(point, vector, track.charge(), &*bField);
// FIX THIS !!!
// need to convert from perigee to global or helix (curvilinear) frame
// for now just an arbitrary matrix.
AlgebraicSymMatrix66 covT= AlgebraicMatrixID(); covT *= 1e-6; // initialize to sigma=1e-3
CartesianTrajectoryError tCov(covT);
return FreeTrajectoryState(tPars, tCov);
}
| FreeTrajectoryState HTrackAssociator::getFreeTrajectoryState | ( | const edm::EventSetup & | iSetup, |
| const SimTrack & | track, | ||
| const SimVertex & | vertex | ||
| ) |
Definition at line 621 of file TrackAssociator.cc.
References ecalTB2006H4_GenSimDigiReco_cfg::bField, DeDxDiscriminatorTools::charge(), alignCSCRings::e, edm::EventSetup::get(), CoreSimTrack::momentum(), point, CoreSimVertex::position(), and CoreSimTrack::type().
{
edm::ESHandle<MagneticField> bField;
iSetup.get<IdealMagneticFieldRecord>().get(bField);
GlobalVector vector( track.momentum().x(), track.momentum().y(), track.momentum().z() );
// convert mm to cm
GlobalPoint point( vertex.position().x()*.1, vertex.position().y()*.1, vertex.position().z()*.1 );
int charge = track.type( )> 0 ? -1 : 1;
GlobalTrajectoryParameters tPars(point, vector, charge, &*bField);
AlgebraicSymMatrix66 covT= AlgebraicMatrixID(); covT *= 1e-6; // initialize to sigma=1e-3
CartesianTrajectoryError tCov(covT);
return FreeTrajectoryState(tPars, tCov);
}
| double HTrackAssociator::getHcalEnergy | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| const FreeTrajectoryState & | trackOrigin, | ||
| const double | dR = -1 |
||
| ) |
associate ECAL only and return energy negative dR means only crossed elements
Definition at line 211 of file TrackAssociator.cc.
References HTrackAssociator::HAssociatorParameters::dRHcal, HTrackDetMatchInfo::hcalConeEnergyFromRecHits(), HTrackDetMatchInfo::hcalEnergyFromRecHits(), info, Parameters::parameters, and HTrackAssociator::HAssociatorParameters::useEcal.
{
HAssociatorParameters parameters;
parameters.useEcal = false;
parameters.dRHcal = dR;
HTrackDetMatchInfo info( associate(iEvent, iSetup, trackOrigin, parameters ));
if (dR>0)
return info.hcalConeEnergyFromRecHits();
else
return info.hcalEnergyFromRecHits();
}
| math::XYZPoint HTrackAssociator::getPoint | ( | const GlobalPoint & | point | ) | [inline, private] |
Definition at line 173 of file TrackAssociator.h.
References PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{
return math::XYZPoint(point.x(),point.y(),point.z());
}
| math::XYZPoint HTrackAssociator::getPoint | ( | const LocalPoint & | point | ) | [inline, private] |
Definition at line 178 of file TrackAssociator.h.
References PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{
return math::XYZPoint(point.x(),point.y(),point.z());
}
| math::XYZVector HTrackAssociator::getVector | ( | const GlobalVector & | vec | ) | [inline, private] |
Definition at line 183 of file TrackAssociator.h.
References PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{
return math::XYZVector(vec.x(),vec.y(),vec.z());
}
| math::XYZVector HTrackAssociator::getVector | ( | const LocalVector & | vec | ) | [inline, private] |
Definition at line 188 of file TrackAssociator.h.
References PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{
return math::XYZVector(vec.x(),vec.y(),vec.z());
}
| void HTrackAssociator::init | ( | const edm::EventSetup & | iSetup | ) | [private] |
Definition at line 117 of file TrackAssociator.cc.
References anyDirection, SteppingHelixPropagator::applyRadX0Correction(), ecalTB2006H4_GenSimDigiReco_cfg::bField, edm::EventSetup::get(), SteppingHelixPropagator::setMaterialMode(), and SteppingHelixPropagator_cfi::SteppingHelixPropagator.
{
// access the calorimeter geometry
iSetup.get<CaloGeometryRecord>().get(theCaloGeometry_);
if (!theCaloGeometry_.isValid())
throw cms::Exception("FatalError") << "Unable to find IdealGeometryRecord in event!\n";
if (useDefaultPropagator_ && ! defProp_ ) {
// setup propagator
edm::ESHandle<MagneticField> bField;
iSetup.get<IdealMagneticFieldRecord>().get(bField);
SteppingHelixPropagator* prop = new SteppingHelixPropagator(&*bField,anyDirection);
prop->setMaterialMode(false);
prop->applyRadX0Correction(true);
// prop->setDebug(true); // tmp
defProp_ = prop;
setPropagator(defProp_);
}
}
| void HTrackAssociator::setPropagator | ( | Propagator * | ptr | ) |
use a user configured propagator
Definition at line 101 of file TrackAssociator.cc.
| void HTrackAssociator::useDefaultPropagator | ( | ) |
use the default propagator
Definition at line 110 of file TrackAssociator.cc.
{
useDefaultPropagator_ = true;
}
Definition at line 201 of file TrackAssociator.h.
std::vector<std::string> HTrackAssociator::CaloTowerCollectionLabels [private] |
CaloTowers.
Definition at line 212 of file TrackAssociator.h.
std::vector<std::vector<std::set<uint32_t> > >* HTrackAssociator::caloTowerMap_ [private] |
Definition at line 197 of file TrackAssociator.h.
int HTrackAssociator::debug_ [private] |
Definition at line 196 of file TrackAssociator.h.
Propagator* HTrackAssociator::defProp_ [private] |
Definition at line 194 of file TrackAssociator.h.
std::vector<std::string> HTrackAssociator::EBRecHitCollectionLabels [private] |
Labels of the detector EDProducts (empty by default) ECAL
Definition at line 207 of file TrackAssociator.h.
Definition at line 199 of file TrackAssociator.h.
std::vector<std::string> HTrackAssociator::EERecHitCollectionLabels [private] |
Definition at line 208 of file TrackAssociator.h.
std::vector<std::string> HTrackAssociator::HBHERecHitCollectionLabels [private] |
HCAL.
Definition at line 210 of file TrackAssociator.h.
Definition at line 200 of file TrackAssociator.h.
Propagator* HTrackAssociator::ivProp_ [private] |
Definition at line 193 of file TrackAssociator.h.
Definition at line 203 of file TrackAssociator.h.
bool HTrackAssociator::useDefaultPropagator_ [private] |
Definition at line 195 of file TrackAssociator.h.
1.7.3