LocalError.h

Go to the documentation of this file.

#ifndef Geom_LocalError_H
#define Geom_LocalError_H

#include "DataFormats/GeometrySurface/interface/TrivialError.h"
#include <cmath>
#include <iosfwd>

class LocalError {
public:
  LocalError() : thexx(0), thexy(0), theyy(0) {}
  LocalError(InvalidError) : thexx(-9999.e10f), thexy(0), theyy(-9999.e10f) {}

  LocalError( float xx, float xy, float yy) :
    thexx(xx), thexy(xy), theyy(yy) {}

  bool invalid() const { return thexx<-1.e10f;}
  bool valid() const { return !invalid();}


  float xx() const { return thexx;}
  float xy() const { return thexy;}
  float yy() const { return theyy;}

  LocalError scale(float s) const { 
    float s2 = s*s;
    return LocalError(s2*xx(), s2*xy(), s2*yy());
  }

  LocalError rotate(float x, float y) const { 
    return rotateCosSin( x, y, 1.f/(x*x+y*y) );
  }

  LocalError rotate(float phi) const { 
    return rotateCosSin( cos(phi), sin(phi));
  }

  LocalError rotateCosSin( float c, float s, float mag2i=1.f) const {
    return LocalError( mag2i*( (c*c)*xx() + (s*s)*yy() - 2.f*(c*s)*xy()),
               mag2i*( (c*s)*(xx() - yy()) +  (c*c-s*s)*xy()) ,
               mag2i*( (s*s)*xx() + (c*c)*yy() + 2.f*(c*s)*xy())
                     );
  }

private:

  float thexx;
  float thexy;
  float theyy;

};  

std::ostream & operator<<( std::ostream& s, const LocalError& err) ;

#endif // LocalError_H