diff --git a/benchmarks/rich/src/CherenkovPIDAnalysis.cc b/benchmarks/rich/src/CherenkovPIDAnalysis.cc index 795dbe814d9a7640b253fddec5f7570063a4ef47..3a7c55e20dd11901591e2034385d19aec1075925 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`