CFAlgoQRatio.cxx

Go to the documentation of this file.

#include "larreco/RecoAlg/CMTool/CMTAlgMatch/CFAlgoQRatio.h"

namespace cmtool {

  //-------------------------------------------------------
  CFAlgoQRatio::CFAlgoQRatio() : CFloatAlgoBase()
  //-------------------------------------------------------
  {
    _qratio_cut = 0.1; // Preliminary cuts
  }

  //-----------------------------
  void
  CFAlgoQRatio::Reset()
  //-----------------------------
  {}

  //-----------------------------------------------------------------------------------------
  float
  CFAlgoQRatio::Float(util::GeometryUtilities const&,
                      const std::vector<const cluster::ClusterParamsAlg*>& clusters)
  //-----------------------------------------------------------------------------------------
  {

    // Code-block by Kazu starts
    // This ensures the algorithm works only if # clusters is > 2 (and not =2)
    // You may take out this block if you want to allow matching using clusters from only 2 planes.
    if (clusters.size() == 2) return -1;
    // Code-block by Kazu ends

    double q_max = -1;
    float ratio = 1;

    // Search for max charge & return if any cluster has negative charge
    for (auto const& c : clusters) {

      //Ignore all clusters with less than 40 hits

      if (c->GetParams().sum_charge < 0) {
        if (_verbose)
          std::cout << "\033[00m Found a cluster with negative charge!\033[00m ... aborting "
                    << __FUNCTION__ << std::endl;
        return -1;
      }

      if (q_max < c->GetParams().sum_charge) q_max = c->GetParams().sum_charge;
    }

    // Compute ratio
    for (auto const& c : clusters) {
      ratio *= c->GetParams().sum_charge / q_max;

      if (_verbose) {
        std::cout << "Sum Charge: " << c->GetParams().sum_charge;
        std::cout << "Q max : " << q_max << std::endl;
        std::cout << "Ratio is: " << ratio << std::endl;
      }
    }

    if (_verbose) {
      if (ratio > _qratio_cut)
        std::cout << " ... pass!" << std::endl;
      else
        std::cout << " ... below cut value: " << _qratio_cut << std::endl;
    }

    return (ratio > _qratio_cut ? ratio : -1);
  }

  //------------------------------
  void
  CFAlgoQRatio::Report()
  //------------------------------
  {}

}