Incident angle in Benchmark

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`