heppy::scaleFunction50< T > Class Template Reference

#include <MuScleFitCorrector_Functions.h>

Inheritance diagram for heppy::scaleFunction50< T >:

Public Member Functions
virtual double	scale (const double &pt, const double &eta, const double &phi, const int chg, const T &parScale) const
	scaleFunction50 ()
Public Member Functions inherited from heppy::scaleFunctionBase< T >
virtual int	parNum () const
virtual	~scaleFunctionBase ()=0

Additional Inherited Members
Protected Member Functions inherited from heppy::scaleFunctionBase< T >
virtual void	setPar (double *Start, double *Step, double *Mini, double *Maxi, int *ind, TString *parname, const T &parResol, const std::vector< int > &parResolOrder, const std::vector< ParSet > &parSet)
Protected Attributes inherited from heppy::scaleFunctionBase< T >
int	parNum_

Detailed Description

template<class T>
class heppy::scaleFunction50< T >

Definition at line 82 of file MuScleFitCorrector_Functions.h.

Constructor & Destructor Documentation

template<class T >

heppy::scaleFunction50< T >::scaleFunction50 ( )

inline

Definition at line 84 of file MuScleFitCorrector_Functions.h.

References heppy::scaleFunctionBase< T >::parNum_.

{ this->parNum_ = 27; }

Member Function Documentation

template<class T >

virtual double heppy::scaleFunction50< T >::scale ( const double &  pt, const double &  eta, const double &  phi, const int  chg, const T &  parScale  ) const

inlinevirtual

Implements heppy::scaleFunctionBase< T >.

Definition at line 85 of file MuScleFitCorrector_Functions.h.

References PVValHelper::eta, PVValHelper::phase2, and funct::sin().

                                                                                                        {
      double ampl(0), phase(0), twist(0), ampl2(0), freq2(0), phase2(0);

      // very bwd bin
      if (eta < parScale[4]) {
        ampl = parScale[1];
        phase = parScale[2];
        ampl2 = parScale[21];
        freq2 = parScale[22];
        phase2 = parScale[23];
        twist =
            parScale[3] * (eta - parScale[4]) + parScale[7] * (parScale[4] - parScale[8]) + parScale[11] * parScale[8];
        // bwd bin
      } else if (parScale[4] <= eta && eta < parScale[8]) {
        ampl = parScale[5];
        phase = parScale[6];
        twist = parScale[7] * (eta - parScale[8]) + parScale[11] * parScale[8];
        // barrel bin
      } else if (parScale[8] <= eta && eta < parScale[12]) {
        ampl = parScale[9];
        phase = parScale[10];
        twist = parScale[11] * eta;
        // fwd bin
      } else if (parScale[12] <= eta && eta < parScale[16]) {
        ampl = parScale[13];
        phase = parScale[14];
        twist = parScale[15] * (eta - parScale[12]) + parScale[11] * parScale[12];
        // very fwd bin
      } else if (parScale[16] < eta) {
        ampl = parScale[17];
        phase = parScale[18];
        ampl2 = parScale[24];
        freq2 = parScale[25];
        phase2 = parScale[26];
        twist = parScale[19] * (eta - parScale[16]) + parScale[15] * (parScale[16] - parScale[12]) +
                parScale[11] * parScale[12];
      }

      // apply the correction
      double curv = (1. + parScale[0]) * ((double)chg / pt - twist - ampl * sin(phi + phase) -
                                          ampl2 * sin(freq2 * phi + phase2) - 0.5 * parScale[20]);
      return 1. / ((double)chg * curv);
    }