MTDEtaRange.cc

Go to the documentation of this file.

#include "RecoMTD/Navigation/interface/MTDEtaRange.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include<iostream>


MTDEtaRange::MTDEtaRange() : 
   theMin(0), theMax(0) {}
   
MTDEtaRange::MTDEtaRange(float max, float min) {
 if ( max < min ) {
   edm::LogWarning ("MTDEtaRange") << "Warning MTDEtaRange:: max < min!! correcting" <<std::endl;
   float tmp(min);
   min = max;
   max = tmp;
 }
 theMax = max;
 theMin = min;
}

MTDEtaRange::MTDEtaRange(const MTDEtaRange& range) :
    theMin(range.theMin), theMax(range.theMax) {}
MTDEtaRange& MTDEtaRange::operator=(const MTDEtaRange& range) {

  if ( this != &range ) {
    theMin = range.theMin;
    theMax = range.theMax;
  }
  return *this;
}

bool MTDEtaRange::isInside(float eta, float error) const {

  if ( (eta+error) > max() || (eta-error) < min() ) return false;
  return true;
}
bool MTDEtaRange::isInside(const MTDEtaRange& range) const {
  if ( min() > range.min() && max() < range.max() ) return true;
  return false;
}
bool MTDEtaRange::isCompatible(const MTDEtaRange& range) const {
  if ( range.min() > max() || range.max() < min() ) return false; 
  return true;
}
MTDEtaRange MTDEtaRange::add(const MTDEtaRange& range) const {
  float max = ( theMax > range.theMax ) ? theMax : range.theMax;
  float min = ( theMin < range.theMin ) ? theMin : range.theMin;
  return MTDEtaRange(max,min);
}
MTDEtaRange MTDEtaRange::subtract(const MTDEtaRange& range) const {

  if ( range.isInside(*this) ) {
    edm::LogInfo ("MTDEtaRange") << "MTDEtaRange: range is inside!" << std::endl;
    return *this;
  }
  if ( !range.isCompatible(*this) ) {
    edm::LogInfo ("MTDEtaRange") << "MTDEtaRange: no overlap between ranges" << std::endl;
    return *this;
  }

  float max = isInside(range.theMin) ? range.theMin : theMax;
  float min = isInside(range.theMax) ? range.theMax : theMin;
  return MTDEtaRange(max,min);
}