From bdb9e715d8fde16cc7f5a3fbf503eb1edefb0555 Mon Sep 17 00:00:00 2001
From: Fangyi Guo <guofangyi@ihep.ac.cn>
Date: Wed, 4 Dec 2024 16:01:35 +0800
Subject: [PATCH] Debug: assign correct Hcal hit energy in digitization
---
Digitization/DigiCalo/src/HcalDigiAlg.cpp | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/Digitization/DigiCalo/src/HcalDigiAlg.cpp b/Digitization/DigiCalo/src/HcalDigiAlg.cpp
index 4650ebe5..35332c9e 100644
--- a/Digitization/DigiCalo/src/HcalDigiAlg.cpp
+++ b/Digitization/DigiCalo/src/HcalDigiAlg.cpp
@@ -266,7 +266,8 @@ StatusCode HcalDigiAlg::execute()
//printf(" Step #%d: En %.2f, rotate angle %.2f, tile pos after rotation (%.2f, %.2f, %.2f), ", iCont, conb.getEnergy(), rotPhi, rot_tilepos.x(), rot_tilepos.y(), rot_tilepos.z());
//printf("rel pos after rotation (%.2f, %.2f, %.2f) \n", (rot_steppos-rot_tilepos).x(), (rot_steppos-rot_tilepos).y(), (rot_steppos-rot_tilepos).z());
//printf(" Project to bin (%d, %d), LY %.3f \n", ibinx, ibiny, GSTileResMap->GetBinContent( ibinx, ibiny ));
-
+
+ Npe_att += conb.getEnergy() / _MIPCali * GSTileResMap->GetBinContent( ibinx, ibiny ) * fLY_tempScale;
m_step_LY.push_back(GSTileResMap->GetBinContent( ibinx, ibiny ));
}
Ehit = Npe_att / _MIPLY / fLY_tempScale * _MIPCali;
--
GitLab