include/N3AnaBase/cuts/energy_cut.h

// -*- mode: c++; -*-
// energy_cut.h

#ifndef __N3AnaBase__energy_cut_h
#define __N3AnaBase__energy_cut_h 1

#include <cstdlib>
#include <iostream>
#include <string>
#include <sstream>
#include <stdexcept>

#include <boost/numeric/interval.hpp>

#include <N3AnaBase/cuts/i_ana_event_cut.h>

namespace nemo3 {

  class energy_cut : public nemo3::i_ana_event_cut
  {
  public:
    enum mode_t
      {
        MODE_INVALID     = -1, 
        MODE_ENERGY_MIN  = 0,  
        MODE_ENERGY_MAX  = 1,  
        MODE_ENERGY_SUM  = 2,  
        MODE_ENERGY_EACH = 3,  

        MODE_UNASSOCIATED_ENERGY_MIN  = 4, 
        MODE_UNASSOCIATED_ENERGY_MAX  = 5, 
        MODE_UNASSOCIATED_ENERGY_SUM  = 6, 
        MODE_UNASSOCIATED_ENERGY_EACH = 7, 

        MODE_TRACK_ENERGY_MIN  = 8,  
        MODE_TRACK_ENERGY_MAX  = 9,  
        MODE_TRACK_ENERGY_SUM  = 10, 
        MODE_TRACK_ENERGY_EACH = 11, 

        MODE_GAMMA_TRACK_ENERGY_MIN  = 12, 
        MODE_GAMMA_TRACK_ENERGY_MAX  = 13, 
        MODE_GAMMA_TRACK_ENERGY_SUM  = 14, 
        MODE_GAMMA_TRACK_ENERGY_EACH = 15, 

        MODE_GAMMA_TRACKED_TRACK_ENERGY_MIN  = 16, 
        MODE_GAMMA_TRACKED_TRACK_ENERGY_MAX  = 17, 
        MODE_GAMMA_TRACKED_TRACK_ENERGY_SUM  = 18, 
        MODE_GAMMA_TRACKED_TRACK_ENERGY_EACH = 19, 

        MODE_ENERGY_CUT_LAST, 

      };

    static const double RANGE_MIN; 
    static const double RANGE_MAX; 

    static int get_mode_from_label (const std::string & label_); 

    static std::string get_label_from_mode (int mode_); 

    static bool is_mode_valid (int mode_); 


    typedef boost::numeric::interval<double> range_t; 

    static std::ostream & print_range (std::ostream & out_,
                                       const range_t & range_);

  private:
    bool    __debug;
    int     __mode;
    range_t __range;
    bool    __inverted; 
    double  __proba;    

  public:

    bool is_inverted () const; 

    bool is_mode_energy_min () const; 

    bool is_mode_energy_max () const; 

    bool is_mode_energy_sum () const; 

    bool is_mode_energy_each () const; 

    bool is_mode (std::string mode_) const; 

    energy_cut (bool debug_ = false); 

    virtual ~energy_cut (); 

    void set (int mode_,
              double min_,
              double max_,
              double proba_,
              bool inverted_ = false);

  public:

    bool is_debug () const; 

    void set_debug (bool); 

    void dump (std::ostream & out_ = std::clog) const; 

  public:

    bool accept ();

  protected:
    virtual void __prepare_cut (); 

  private:
    static creator_registration<energy_cut> __CR;

  public:

    static const std::string CUT_ID;

    virtual std::string cut_id () const;

    virtual cut_creator_t cut_creator () const;

    static nemo3::i_ana_event_cut * create (const nemocuts::ICutFactory::cut_properties_t & cut_properties_); 

  };

} // end of namespace nemo3

#endif /* __N3AnaBase_energy__cut_h */

// end of energy_cut.h