Public Member Functions
	ParametrizationTauJet (int ptype, std::vector< double > x, double u)

double	value (double, double) const

Private Attributes
std::map< int, double >	theEtabound

std::map< int, std::vector < double > >	theParam

int	type

Detailed Description

Definition at line 28 of file TauJetCorrector.h.

Constructor & Destructor Documentation

TauJetCorrector::ParametrizationTauJet::ParametrizationTauJet	(	int	ptype,
		std::vector< double >	x,
		double	u
	)

inline

Definition at line 30 of file TauJetCorrector.h.

References theEtabound, theParam, type, and vdt::x.

                                                                     {
             type=ptype;
             theParam[type] = x;
             theEtabound[type] = u;
             //cout<<"ParametrizationTauJet "<<type<<" "<<u<<endl;
     };

Member Function Documentation

double TauJetCorrector::ParametrizationTauJet::value	(	double	et,
		double	eta
	)		const

Definition at line 15 of file TauJetCorrector.cc.

References edm::second(), mathSSE::sqrt(), and vdt::x.

                                                                             {
 
   double x=et;
   double etnew(et);
   double etabound = (*theEtabound.find(type)).second;
   std::vector<double> taus = (*theParam.find(type)).second;
   
        if ( fabs(eta) > etabound) {
           //cout << " ===> ETA outside of range - CORRECTION NOT DONE ***" << endl;
           //cout << " eta = " << eta <<" pt = " << et << " etabound "<<etabound<<endl;
           return et;
         }
   //cout << "correction parameters: " << taus[0] << " " << taus[1] << " " << taus[2] << endl;
   switch(type){
   case 1:
     {
       etnew = 2*(x-taus[0])/(taus[1]+sqrt(taus[1]*taus[1]-4*taus[0]*taus[2]+4*x*taus[2]));
       break;
     }
   case 2:
     {
       etnew = 2*(x-taus[0])/(taus[1]+sqrt(taus[1]*taus[1]-4*taus[0]*taus[2]+4*x*taus[2]));
       break;
     }
   case 3:
     {
       etnew = 2*(x-taus[0])/(taus[1]+sqrt(taus[1]*taus[1]-4*taus[0]*taus[2]+4*x*taus[2]));
       break;
     }
 
   default:
     edm::LogError("TauJetCorrector: Error: unknown parametrization type ") << type << " in TauJetCorrector. No correction applied" << endl;
     //cerr<<"TauJetCorrector: Error: unknown parametrization type '"<<type<<"' in TauJetCorrector. No correction applied"<<endl;
     break;
   }
   return etnew;
 }