ClusterParams.h

Go to the documentation of this file.

#ifndef UTIL_CLUSTERPARAMS_H
#define UTIL_CLUSTERPARAMS_H

#include "lardata/Utilities/PxUtils.h"
#include "larreco/RecoAlg/ClusterRecoUtil/Polygon2D.h"

namespace cluster {

  class cluster_params {
  public:
    cluster_params()
    {
      //Initialize all values.
      Clear();
    }

    Polygon2D PolyObject; 

    util::PxPoint start_point; 
    util::PxPoint end_point;   

    double sum_charge;       
    double mean_charge;      
    double rms_charge;       
    double sum_ADC;          
    double mean_ADC;         
    double rms_ADC;          
    double mean_x;           
    double mean_y;           
    double rms_x;            
    double rms_y;            
    double charge_wgt_x;     
    double charge_wgt_y;     
    double cluster_angle_2d; 
    double angle_2d;         
    double opening_angle;    
    double opening_angle_charge_wgt; 
    double closing_angle;            
    double closing_angle_charge_wgt; 
    double start_charge;             
    double end_charge;               
    double eigenvalue_principal;     
    double eigenvalue_secondary;     
    double verticalness;             

    double length;
    double width;

    double hit_density_1D;
    double hit_density_2D;
    double multi_hit_wires;
    double N_Wires;
    double modified_hit_density;
    double N_Hits;
    double N_Hits_HC;

    double modmeancharge;
    double RMS_charge;

    int direction;

    double showerness;   
    double trackness;    
    double offaxis_hits; 

    void Clear()
    {
      start_point.Clear();
      end_point.Clear();
      sum_charge = -999.999;
      mean_charge = -999.999;
      rms_charge = -999.999;
      sum_ADC = -999.999;
      mean_ADC = -999.999;
      rms_ADC = -999.999;
      mean_x = -999.999;
      mean_y = -999.999;
      rms_x = -999.999;
      rms_y = -999.999;
      charge_wgt_x = -999.999;
      charge_wgt_y = -999.999;
      cluster_angle_2d = -999.999;
      angle_2d = -999.999;
      eigenvalue_principal = -999.999;
      eigenvalue_secondary = -999.999;
      multi_hit_wires = -999.999;
      N_Wires = -999.999;
      verticalness = -999.999;
      opening_angle = -999.999;
      opening_angle_charge_wgt = -999.999;
      closing_angle = -999.999;
      closing_angle_charge_wgt = -999.999;
      start_charge = -999.999;
      end_charge = -999.999;
      offaxis_hits = -999.999;
      modmeancharge = -999.999;
      RMS_charge = -999.999;
      modified_hit_density = -999.999;
      N_Hits = -999.999;
      N_Hits_HC = -999.999;
      length = -999.999;
      width = -999.999;
      hit_density_1D = -999.999;
      hit_density_2D = -999.999;
      direction = 0;
      showerness = -999.999;
      trackness = -999.999;
    }

    template <typename Stream>
    void Report(Stream& os) const
    {

      os << "ClusterParams report: \n"
         << "   start_point  .............. : (" << start_point.w << "," << start_point.t << ")\n"
         << "   end_point  ................ : (" << end_point.w << "," << end_point.t << ")\n"
         << "   sum_charge  ............... : " << sum_charge << "\n"
         << "   mean_charge   ............. : " << mean_charge << "\n"
         << "   mean_x  ................... : " << mean_x << "\n"
         << "   mean_y  ................... : " << mean_y << "\n"
         << "   rms_x  .................... : " << rms_x << "\n"
         << "   rms_y  .................... : " << rms_y << "\n"
         << "   charge_wgt_x  ............. : " << charge_wgt_x << "\n"
         << "   charge_wgt_y  ............. : " << charge_wgt_y << "\n"
         << "   cluster_angle_2d  ......... : " << cluster_angle_2d << "\n"
         << "   angle_2d  ................. : " << angle_2d << "\n"
         << "   eigenvalue_principal  ..... : " << eigenvalue_principal << "\n"
         << "   eigenvalue_secondary  ..... : " << eigenvalue_secondary << "\n"
         << "   multi_hit_wires  .......... : " << multi_hit_wires << "\n"
         << "   N_Wires  .................. : " << N_Wires << "\n"
         << "   N_Hits  ................... : " << N_Hits << "\n"
         << "   N_Hits_HC  ................ : " << N_Hits_HC << "\n"
         << "   verticalness  ............. : " << verticalness << "\n"
         << "   opening_angle  ............ : " << opening_angle << "\n"
         << "   opening_angle_charge_wgt  . : " << opening_angle_charge_wgt << "\n"
         << "   closing_angle  ............ : " << closing_angle << "\n"
         << "   closing_angle_charge_wgt  . : " << closing_angle_charge_wgt << "\n"
         << "   start_charge  ............. : " << start_charge << "\n"
         << "   end_charge  ............... : " << end_charge << "\n"
         << "   offaxis_hits  ............. : " << offaxis_hits << "\n"
         << "   modified_hit_density  ..... : " << modified_hit_density << "\n"
         << "   modified mean charge ...... : " << modmeancharge << "\n"
         << "   charge_RMS ................ : " << RMS_charge << "\n"
         << "   length  ................... : " << length << "\n"
         << "   width  .................... : " << width << "\n"
         << "   hit_density_1D  ........... : " << hit_density_1D << "\n"
         << "   hit_density_2D  ........... : " << hit_density_2D << "\n"
         << "   showerness  ............... : " << showerness << "\n"
         << "   trackness  ................ : " << trackness << "\n"
         << "   direction  ................ : " << direction << "\n";
    }

  }; // class cluster_params
} // namespace cluster

#endif