diff --git a/benchmarks/zdc_sigma/analysis/gen_sigma_decay.cxx b/benchmarks/zdc_sigma/analysis/gen_sigma_decay.cxx index 4e808434441f6115b9ebfbdf3c83672c3769f9af..ee708209b3f5caa229c212595258f3f6338d6be3 100644 --- a/benchmarks/zdc_sigma/analysis/gen_sigma_decay.cxx +++ b/benchmarks/zdc_sigma/analysis/gen_sigma_decay.cxx @@ -291,11 +291,11 @@ void gen_sigma_decay(int n_events = 100000, UInt_t seed = 0, char* out_fname = " TVector3 extrap_gamma2=lambda_decay_position+gamma2_lab.Vect()*((zdc_z-pbeam_dir.Dot(lambda_decay_position))/(pbeam_dir.Dot(gamma2_lab.Vect()))); if (extrap_neutron.Angle(pbeam_dir)<0.004 && extrap_gamma1.Angle(pbeam_dir)<0.004 && extrap_gamma2.Angle(pbeam_dir)<0.004 && extrap_gamma.Angle(pbeam_dir)<0.004 && lambda_decay_position.Dot(pbeam_dir)<zdc_z){ hepmc_output.write_event(evt); + if (accepted_events % 1000 == 0) { + std::cout << "Event: " << accepted_events << std::endl; + } accepted_events ++; } - if (accepted_events % 1000 == 0) { - std::cout << "Event: " << accepted_events << std::endl; - } evt.clear(); } hepmc_output.close();