GammaZInterference.h

Go to the documentation of this file.

#ifndef PhysicsTools_Utilities_GammaZInterference_h
#define PhysicsTools_Utilities_GammaZInterference_h
#include "PhysicsTools/Utilities/interface/Parameter.h"
#include <boost/shared_ptr.hpp>
#include <cmath>

namespace funct {

  struct GammaZInterference {
    GammaZInterference(const Parameter& m, const Parameter& g): 
      mass(m.ptr()), width(g.ptr()) { }
    GammaZInterference(boost::shared_ptr<double> m, boost::shared_ptr<double> g): 
      mass(m), width(g) {}
    double operator()(double x) const { 
      double m2 = *mass * (*mass); 
      double g2 = *width * (*width);
      double g2OverM2 = g2/m2; 
      double s = x*x;
      double deltaS = s - m2;
      double interference = 0;
      if (fabs(deltaS/m2)<16) {
    double prop = deltaS*deltaS + s*s*g2OverM2;
    interference =  5*(*mass)*deltaS/prop;
      }
      return interference;
    }
    boost::shared_ptr<double> mass, width;
  };

}

#endif