MuonEtaRange.cc

Go to the documentation of this file.

#include "RecoMuon/Navigation/interface/MuonEtaRange.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include<iostream>


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

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

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

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

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

  if ( range.isInside(*this) ) {
    edm::LogInfo ("MuonEtaRange") << "MuonEtaRange: range is inside!" << std::endl;
    return *this;
  }
  if ( !range.isCompatible(*this) ) {
    edm::LogInfo ("MuonEtaRange") << "MuonEtaRange: 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 MuonEtaRange(max,min);
}