Skip to content

Commit

Permalink
Cell position cast to float
Browse files Browse the repository at this point in the history
  • Loading branch information
kjvbrt committed Apr 5, 2024
1 parent 0d670f1 commit 361e1da
Show file tree
Hide file tree
Showing 2 changed files with 10 additions and 7 deletions.
3 changes: 2 additions & 1 deletion RecCalorimeter/src/components/LayeredCaloTowerTool.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -81,7 +81,8 @@ void LayeredCaloTowerTool::towersNumber(int& nEta, int& nPhi) {
nPhi = m_nPhiTower;
}

uint LayeredCaloTowerTool::buildTowers(std::vector<std::vector<float>>& aTowers, bool fillTowerCells) {
uint LayeredCaloTowerTool::buildTowers(std::vector<std::vector<float>>& aTowers,
[[maybe_unused]] bool fillTowerCells) {
// Get the input collection with cells from simulation + digitisation (after
// calibration and with noise)
const edm4hep::CalorimeterHitCollection* cells = m_cells.get();
Expand Down
14 changes: 8 additions & 6 deletions RecFCCeeCalorimeter/src/components/CaloTopoClusterFCCee.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -189,9 +189,10 @@ StatusCode CaloTopoClusterFCCee::execute() {
else
warning() << "No cell positions tool found for system id " << systemId << ". " << endmsg;

newCell.setPosition(edm4hep::Vector3f{posCell.X() / dd4hep::mm,
posCell.Y() / dd4hep::mm,
posCell.Z() / dd4hep::mm});
newCell.setPosition(edm4hep::Vector3f{
static_cast<float>(posCell.X() / dd4hep::mm),
static_cast<float>(posCell.Y() / dd4hep::mm),
static_cast<float>(posCell.Z() / dd4hep::mm)});
posX += posCell.X() * newCell.getEnergy();
posY += posCell.Y() * newCell.getEnergy();
posZ += posCell.Z() * newCell.getEnergy();
Expand All @@ -209,9 +210,10 @@ StatusCode CaloTopoClusterFCCee::execute() {
info() << "Problem in erasing cell ID from map." << endmsg;
}
cluster.setEnergy(energy);
cluster.setPosition(edm4hep::Vector3f{(posX / energy) / dd4hep::mm,
(posY / energy) / dd4hep::mm,
(posZ / energy) / dd4hep::mm});
cluster.setPosition(edm4hep::Vector3f{
static_cast<float>((posX / energy) / dd4hep::mm),
static_cast<float>((posY / energy) / dd4hep::mm),
static_cast<float>((posZ / energy) / dd4hep::mm)});
// store deltaR of cluster in time for the moment..
sumPhi = sumPhi / energy;
sumTheta = sumTheta / energy;
Expand Down

0 comments on commit 361e1da

Please sign in to comment.