From ed97fd06af1767a813b5c94a257a6d34db43219d Mon Sep 17 00:00:00 2001 From: matteoceschia Date: Mon, 21 Sep 2020 15:15:39 +0100 Subject: [PATCH 1/2] Added a (conservative) one-line fix to take into account NaN radii --- SNFitterModule.cpp | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) diff --git a/SNFitterModule.cpp b/SNFitterModule.cpp index a166655..8737acd 100644 --- a/SNFitterModule.cpp +++ b/SNFitterModule.cpp @@ -122,13 +122,15 @@ FALAISE_REGISTER_MODULE(SNFitterModule) for (const auto &hit : gg_hits_col) { // work with geiger hits as members of a given cluster - GeigerRing ring{hit->get_r(), hit->get_x(), hit->get_y(), + if ( !std::isnan(hit->get_r()) ){ //take into account NaN radii: discard them + GeigerRing ring{hit->get_r(), hit->get_x(), hit->get_y(), hit->get_z(), hit->get_sigma_r(), hit->get_sigma_z()}; - MetaInfo mi{hit->get_id(), hit->get_side(), hit->get_row(), hit->get_layer()}; + MetaInfo mi{hit->get_id(), hit->get_side(), hit->get_row(), hit->get_layer()}; - rings.emplace_back(TrackerHit{ncl, ring, mi}); + rings.emplace_back(TrackerHit{ncl, ring, mi}); + } } - + std::cout<<"***** NaN vec length: "< hres = snf.fithelix(); From 2ee9101d4eeaa8b20b9d525afe2e58c794827e2f Mon Sep 17 00:00:00 2001 From: matteoceschia Date: Mon, 21 Sep 2020 15:17:56 +0100 Subject: [PATCH 2/2] Added a (conservative) one-line fix to take into account NaN radii --- SNFitterModule.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/SNFitterModule.cpp b/SNFitterModule.cpp index 8737acd..1362e5b 100644 --- a/SNFitterModule.cpp +++ b/SNFitterModule.cpp @@ -130,7 +130,7 @@ FALAISE_REGISTER_MODULE(SNFitterModule) rings.emplace_back(TrackerHit{ncl, ring, mi}); } } - std::cout<<"***** NaN vec length: "< hres = snf.fithelix();