CheckRegion.hh

Go to the documentation of this file.

#ifndef CHECKREGION_h 
#define CHECKREGION_h 
#include <stdexcept>
#include "ProcessingCut.hh"
namespace ProcessingCuts{
  
  class CheckRegion : public ProcessingCut{
  public:
    CheckRegion();
    ~CheckRegion() {}
    
    static std::string GetCutName(){ return "CheckRegion";}
    
    bool Process(EventDataPtr event);
    int AddDependenciesToModule(BaseModule* mod);
    enum variable_t {MIN, MAX, INTEGRAL, NPE, AMPLITUDE};
  private:
    bool ProcessChannel(ChannelData* chdata);
    int channel;                 
    int region;                  
    double minimum;                 
    double maximum;                 
    variable_t variable;         
    bool default_pass;           
  };

  inline
  std::ostream& operator<<(std::ostream& out, 
                           const CheckRegion::variable_t& var){
    if(var == CheckRegion::MIN)
      return out<<"min";
    else if(var == CheckRegion::MAX)
      return out<<"max";
    else if(var == CheckRegion::INTEGRAL)
      return out<<"integral";
    else if(var == CheckRegion::NPE)
      return out<<"npe";
    else if(var == CheckRegion::AMPLITUDE)
      return out<<"amplitude";
    return out;
  }
  inline
  std::istream& operator>>(std::istream& in, CheckRegion::variable_t& var){
    std::string vname;
    in>>vname;
    if(vname == "min") var = CheckRegion::MIN;
    else if (vname == "max") var = CheckRegion::MAX;
    else if (vname == "integral") var = CheckRegion::INTEGRAL;
    else if (vname == "npe") var = CheckRegion::NPE;
    else if (vname == "amplitude") var = CheckRegion::AMPLITUDE;
    else{
      std::cerr<<"Unknown variable type "<<vname<<std::endl;
      throw std::invalid_argument("CheckRegion::variable_t has no entry"+vname);
    }
    return in;
  }
}

#endif