include/N3AnaBase/cuts/multiplicity_cut.h

// multiplicity_cut.h

/* A basic multiplicity cut on NEMO3 ANA event.
 * Cut is performed on some multiplicity among:
 * - number of scintillator hits
 * - number of Geiger hits
 * - number of prompt tracks
 * - number of delayed tracks
 *
 */

#ifndef __n3anabase__multiplicity_cut_h
#define __n3anabase__multiplicity_cut_h 1

#include <string>
#include <algorithm>

#include <boost/numeric/interval.hpp>
#include <N3AnaBase/measurements/cluster.h>

#include <N3AnaBase/cuts/i_ana_event_cut.h> 

namespace nemo3 {
  
  class multiplicity_cut : public i_ana_event_cut
    {
    public:
      enum mode_t
        {
          MODE_INVALID                    = -1,
          MODE_NB_SCIN                    = 0,
          MODE_NB_SCIN_FROM_INTERNAL_WALL = 1,
          MODE_NB_SCIN_FROM_EXTERNAL_WALL = 2,
          MODE_NB_SCIN_FROM_BOTTOM_CAP    = 3,
          MODE_NB_SCIN_FROM_TOP_CAP       = 4,
          MODE_NB_SCIN_ASSOCIATED         = 5,
          MODE_NB_SCIN_UNASSOCIATED       = 6,
          MODE_NB_CLUSTERS                = 7,
          MODE_NB_GG                      = 8,
          MODE_NB_GG_UNASSOCIATED         = 9,
          MODE_NB_DELAYED_HITS            = 10,
          MODE_NB_TRACKS                  = 11,
          MODE_NB_DELAYED_TRACKS          = 12
        };

      static const int32_t RANGE_MIN;
      static const int32_t RANGE_MAX;

      static const std::string CUT_ID;
      
      //typedef boost::numeric::interval<int> range_t;
      
    private:
      static creator_registration<multiplicity_cut> __CR;

    private:
      bool __debug;
      int  __mode;
      std::string __tracking_name;
      //range_t __range;
      int  __min_multiplicity;
      int  __max_multiplicity;
      bool __reversed;

    public:

      bool is_debug () const
        {
          return __debug;
        }

      bool is_reversed () const
        {
          return __reversed;
        }

      multiplicity_cut (bool debug_ = false);

      virtual ~multiplicity_cut ();

      void set (int mode_,
                int min_multiplicity_, 
                int max_multiplicity_,
                bool reverse_ = false);

      bool accept (); 

      static bool in_range (int value_, int min_, int max_);

      static bool not_in_range (int value_, int min_, int max_);


    protected:
      virtual void __prepare_cut ();
      
    public:

      // used in factory method:
      virtual std::string cut_id () const;

      virtual cut_creator_t cut_creator () const;

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

      static void test ();

    };
  
} // end of namespace nemo3

#endif // __n3anabase__multiplicity_cut_h

// end of multiplicity_cut.h