IMASelector.h

Go to the documentation of this file.

#ifndef UtilAlgos_IMASelector_h
#define UtilAlgos_IMASelector_h
#include "CommonTools/UtilAlgos/interface/ParameterAdapter.h"
#include <iostream>


struct IMASelector {
  IMASelector( double ESCOPinMin, double ESeedOPoutMin, double PinMPoutOPinMin,
               double ESCOPinMax, double ESeedOPoutMax, double PinMPoutOPinMax,
           double EMPoutMin, double EMPoutMax) : 
    ESCOPinMin_ (ESCOPinMin),
    ESeedOPoutMin_ (ESeedOPoutMin),
    PinMPoutOPinMin_ (PinMPoutOPinMin),
    ESCOPinMax_ (ESCOPinMax),
    ESeedOPoutMax_ (ESeedOPoutMax),
    PinMPoutOPinMax_ (PinMPoutOPinMax),
    EMPoutMin_ (EMPoutMin),
    EMPoutMax_ (EMPoutMax) {}

  template<typename T>
  bool operator()( const T & t ) const { 
    double pin = t.trackMomentumAtVtx ().R () ;
    double poMpiOpi = (pin - t.trackMomentumOut ().R ()) / pin ;
    double ESC = t.energy () ;     
    double pOut = t.trackMomentumOut().R();
    double EseedOPout = t.eSeedClusterOverPout () ;
    double EoPin = t.eSuperClusterOverP () ;
    double EoPout = (ESC)/pOut;
    bool result;
        result=EoPout > EMPoutMin_ ;
        result = EoPout<EMPoutMax_;
    return (poMpiOpi > PinMPoutOPinMin_ && poMpiOpi < PinMPoutOPinMax_ &&
            EseedOPout > ESeedOPoutMin_ && EseedOPout < ESeedOPoutMax_ &&
            EoPin > ESCOPinMin_ && EoPin < ESCOPinMax_ &&
        EoPout > EMPoutMin_ && EoPout < EMPoutMax_) ;
  }
private:
  double ESCOPinMin_, ESeedOPoutMin_, PinMPoutOPinMin_ ;
  double ESCOPinMax_, ESeedOPoutMax_, PinMPoutOPinMax_ ;
  double EMPoutMin_, EMPoutMax_;
};


namespace reco { 
  namespace modules { 
    template<> 
    struct ParameterAdapter<IMASelector> { 
      static IMASelector make(const edm::ParameterSet & cfg) { 
        return IMASelector(cfg.getParameter<double>("ESCOPinMin"), 
                           cfg.getParameter<double>("ESeedOPoutMin"), 
                           cfg.getParameter<double>("PinMPoutOPinMin"), 
                           cfg.getParameter<double>("ESCOPinMax"), 
                           cfg.getParameter<double>("ESeedOPoutMax"), 
                           cfg.getParameter<double>("PinMPoutOPinMax"), 
               cfg.getParameter<double>("EMPoutMin"),
               cfg.getParameter<double>("EMPoutMax")
               );
      } 
    }; 
  } 
} 

#endif