include/N3AnaBase/cuts/track_scin_association_cut.h

// track_scin_association_cut.h

#ifndef __N3AnaBase__track_scin_association_cut_h
#define __N3AnaBase__track_scin_association_cut_h 1

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

#include <boost/numeric/interval.hpp>

#include <N3AnaBase/cuts/i_ana_event_cut.h>

using namespace std;

namespace nemo3 {

  class track_scin_association_cut : public nemo3::i_ana_event_cut
    {
    public:
      enum mode_t
        {
          MODE_EACH      = 0,
          MODE_MULTIPLE  = 1,
          MODE_MORE_THAN = 2,
          MODE_LESS_THAN = 3,
          MODE_INVALID   = -1
        };

      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_);


    private:
      bool    __debug;
      int     __mode;
      int     __number_of_tracks;
      string  __tracking_name;

    public:
      
      bool is_mode_each () const;

      bool is_mode_multiple () const;
   
      bool is_mode_more_than () const;

      bool is_mode_less_than () const;

      track_scin_association_cut (bool debug_ = false); // default ctor

      virtual ~track_scin_association_cut (); // dtor

      void set (int mode_, int nb_track_, const string & tracking_name_ = "");

 
   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 ();


    // used by the cuts' factory stuff:
    private:
      static creator_registration<track_scin_association_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 */