From 2a44633dc77d6dd7115c432c7c2f94bcebef49bc Mon Sep 17 00:00:00 2001
From: Chandradoy Chatterjee <chandradoy.chatterjee@ts.infn.it>
Date: Fri, 28 Jul 2023 15:15:04 +0200
Subject: [PATCH] Incident angle in Benchmark
---
benchmarks/rich/src/CherenkovPIDAnalysis.cc | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/benchmarks/rich/src/CherenkovPIDAnalysis.cc b/benchmarks/rich/src/CherenkovPIDAnalysis.cc
index 795dbe81..3a7c55e2 100644
--- a/benchmarks/rich/src/CherenkovPIDAnalysis.cc
+++ b/benchmarks/rich/src/CherenkovPIDAnalysis.cc
@@ -148,7 +148,7 @@ namespace benchmarks {
auto thrown_pdg = part->GetPDG();
// loop over `CherenkovParticleID` objects
- for(const auto& cherenkov_pid : cherenkov_pids) {
+ for( const auto& cherenkov_pid : cherenkov_pids) {
// skip if NPE==0
if(cherenkov_pid.getNpe() == 0) {
@@ -179,7 +179,9 @@ namespace benchmarks {
theta_rec += theta;
theta_rec /= cherenkov_pid.getNpe();
auto theta_rec_mrad = theta_rec * 1e3; // [rad] -> [mrad]
-
+ for(const auto& iAngle : cherenkov_pid.getIncidentAngle())
+ printf("--->%lf\n",iAngle);
+ //auto nn = cherenkov_pid.getIncidentAngle(cherenkov_pid.incidentAngle_size());
// calculate expected Cherenkov angle `theta_exp` and residual `theta_resid`,
// using refractive index from MC truth
auto mc_rindex = cherenkov_pid.getRefractiveIndex(); // average refractive index for photons used in this `cherenkov_pid`
--
GitLab