resolutionFunctionType18< T > Class Template Reference

#include <Functions.h>

Inheritance diagram for resolutionFunctionType18< T >:

Public Member Functions
	resolutionFunctionType18 ()
virtual void	setParameters (double *Start, double *Step, double *Mini, double *Maxi, int *ind, TString *parname, const T &parResol, const std::vector< int > &parResolOrder, const int muonType)
	This method is used to differentiate parameters among the different functions. More...
virtual double	sigmaCotgTh (const double &pt, const double &eta, const T &parval)
virtual double	sigmaPhi (const double &pt, const double &eta, const T &parval)
virtual double	sigmaPt (const double &pt, const double &eta, const T &parval)
Public Member Functions inherited from resolutionFunctionBase< T >
virtual double	covPt1Pt2 (const double &pt1, const double &eta1, const double &pt2, const double &eta2, const T &parval)
virtual int	parNum () const
	resolutionFunctionBase ()
virtual void	setParameters (double *Start, double *Step, double *Mini, double *Maxi, int *ind, TString *parname, const T &parResol, const std::vector< int > &parResolOrder, const std::vector< double > &parStep, const std::vector< double > &parMin, const std::vector< double > &parMax, const int muonType)
virtual double	sigmaPtError (const double &pt, const double &eta, const T &parval, const T &parError)
virtual	~resolutionFunctionBase ()=0

Additional Inherited Members
Protected Member Functions inherited from resolutionFunctionBase< T >
virtual void	setPar (double *Start, double *Step, double *Mini, double *Maxi, int *ind, TString *parname, const T &parResol, const std::vector< int > &parResolOrder, double *thisStep, double *thisMini, double *thisMaxi, TString *thisParName)
	This method sets the parameters. More...
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< ParameterSet > &parSet)
Protected Attributes inherited from resolutionFunctionBase< T >
int	parNum_

Detailed Description

template<class T>
class resolutionFunctionType18< T >

Definition at line 3898 of file Functions.h.

Constructor & Destructor Documentation

template<class T >

resolutionFunctionType18< T >::resolutionFunctionType18 ( )

inline

Definition at line 3900 of file Functions.h.

References resolutionFunctionBase< T >::parNum_.

{ this->parNum_ = 14; }

Member Function Documentation

template<class T >

virtual void resolutionFunctionType18< T >::setParameters ( double *  Start, double *  Step, double *  Mini, double *  Maxi, int *  ind, TString *  parname, const T &  parResol, const std::vector< int > &  parResolOrder, const int  muonType  )

inlinevirtual

This method is used to differentiate parameters among the different functions.

Reimplemented from resolutionFunctionBase< T >.

Definition at line 3929 of file Functions.h.

References resolutionFunctionBase< T >::setPar().

                                                                                                                                                                                                {
    double thisStep[] = { 0.01, 0.0001, 0.00001,
              0.001, 0.0001, 0.000001,
                          0.01, 0.001, 0.01, 0.001,
              0.01, 0.001, 0.01, 0.001};
    TString thisParName[] = { "offsetPt", "hyperbolicPt", "linearPt",
                              "offsetEtaCentral", "linaerEtaCentral", "parabEtaCentral",
                  "offsetEtaEndcapRight", "linearEtaRight", "rightParabCenter", "parabolicEtaRight",
                              "offsetEtaEndcapLeft", "linearEtaLeft", "leftParabCenter", "parabolicEtaLeft" };
    double thisMini[] = { -1.5, -0.001, 0.00005,
                          -0.05, -0.1, 0.0,
                          -0.6, -0.0009, 0., 0.0005,
              -0.6, -0.1, 1., 0.01
                        };

    if( muonType == 1 ) {
      double thisMaxi[] = { 1., 1., 1.,
                1., 1., 1.,
                            1., 1. ,1., 1.,
                1., 1., 1., 1.
                          };
      this->setPar( Start, Step, Mini, Maxi, ind, parname, parResol, parResolOrder, thisStep, thisMini, thisMaxi, thisParName );
    } else {
      double thisMaxi[] = { 1.1, 0.8, 0.005,
                0.05, 0.1, 0.08,
                            0.9, 0.5, 1.99, 0.15,
                0.9, 0.5, 1.99, 0.15
      };

      this->setPar( Start, Step, Mini, Maxi, ind, parname, parResol, parResolOrder, thisStep, thisMini, thisMaxi, thisParName );
    }
  }

template<class T >

virtual double resolutionFunctionType18< T >::sigmaCotgTh ( const double &  pt, const double &  eta, const T &  parval  )

inlinevirtual

Implements resolutionFunctionBase< T >.

Definition at line 3922 of file Functions.h.

                                                                                      {
    return( 0.004 );
  }

template<class T >

virtual double resolutionFunctionType18< T >::sigmaPhi ( const double &  pt, const double &  eta, const T &  parval  )

inlinevirtual

Implements resolutionFunctionBase< T >.

Definition at line 3926 of file Functions.h.

                                                                                   {
    return( 0.001 );
  }

template<class T >

virtual double resolutionFunctionType18< T >::sigmaPt ( const double &  pt, const double &  eta, const T &  parval  )

inlinevirtual

Implements resolutionFunctionBase< T >.

Definition at line 3902 of file Functions.h.

                                                                                  {
    double fabsEta = std::fabs(eta);
    double ptPart =  parval[0] + parval[1]*1./pt + pt*parval[2];

    if(fabsEta<=0.6)
      return( ptPart + parval[3]);
    else if((eta>0.6 && eta<=1.3) || (eta>=-1.3 && eta<-0.6)) {//eta in barrel + overlap
      double par = parval[3] - 0.6*parval[4] - 0.6*0.6*parval[5];
      return( ptPart + par + parval[4]*fabsEta + parval[5]*eta*eta );
    }
    else if (eta>1.3){//eta in right endcap
      double par = parval[3] - 0.6*parval[4] - 0.6*0.6*parval[5] + parval[4]*1.3 + parval[5]*1.3*1.3 - (parval[6] + parval[7]*std::fabs(1.3-parval[8]) + parval[9]*(1.3-parval[8])*(1.3-parval[8]));
      return( par +  ptPart + parval[6] + parval[7]*std::fabs((fabsEta-parval[8])) + parval[9]*(fabsEta-parval[8])*(fabsEta-parval[8]) );
    }
    else{//eta in left endcap
      double par = parval[3] - 0.6*parval[4] - 0.6*0.6*parval[5] + parval[4]*1.3 + parval[5]*1.3*1.3 - (parval[10] + parval[11]*std::fabs(1.3-parval[12]) + parval[13]*(1.3-parval[12])*(1.3-parval[12]));
      return( par + ptPart + parval[10] + parval[11]*std::fabs((fabsEta-parval[12])) + parval[13]*(fabsEta-parval[12])*(fabsEta-parval[12]) );
    }
  }