Public Member Functions
	ParametrizationJetParton (int thePartonMixture, const vector< double > &x, const vector< double > &y, const vector< double > &z)

double	value (double arg1, double arg2) const

Private Attributes
std::vector< double >	pg

std::vector< double >	pq

std::vector< double >	pqcd

int	type

Detailed Description

Definition at line 25 of file JetPartonCorrector.cc.

Constructor & Destructor Documentation

◆ ParametrizationJetParton()

JetPartonNamespace::ParametrizationJetParton::ParametrizationJetParton	(	int	thePartonMixture,
		const vector< double > &	x,
		const vector< double > &	y,
		const vector< double > &	z
	)

inline

Definition at line 27 of file JetPartonCorrector.cc.

                                                       {
       type = thePartonMixture;
       pq = x;
       pg = y;
       pqcd = z;
     }

Member Function Documentation

◆ value()

double JetPartonNamespace::ParametrizationJetParton::value	(	double	arg1,
		double	arg2
	)		const

Definition at line 46 of file JetPartonCorrector.cc.

                                                                    {
     double enew(e);
     double x = e;
     double y = eta;
  
     if (abs(x - pq[2]) < pq[1] / pq[0] || abs(x - pg[2]) < pg[1] / pg[0] || abs(x - pqcd[2]) < pqcd[1] / pqcd[0]) {
       return enew;
     }
  
     double kjetq = pq[0] - pq[1] / (x - pq[2]);
     double kjetg = pg[0] - pg[1] / (x - pg[2]);
     double kjetqcd = pqcd[0] - pqcd[1] / (x - pqcd[2]);
  
     switch (type) {
       case 1: {
         if (abs(kjetq) > 0.0001)
           enew = e / kjetq;
         break;
       }
       case 2: {
         if (abs(kjetg) > 0.0001)
           enew = e / kjetg;
         break;
       }
       case 3: {
         if (abs(kjetqcd) > 0.0001)
           enew = e / kjetqcd;
         break;
       }
       case 4: {
         cout << "[Jets] JetPartonCorrector: Warning! Calibration to b-quark - does not implemented yet. Light quark "
                 "calibration is used instead "
              << endl;
         if (abs(kjetq) > 0.0001)
           enew = e / kjetq;
         break;
       }
       case 100: {
         UserPartonMixture upm;
         double f = upm.mixt(x, y);
         double kjet = (f * kjetq + kjetg) / (f + 1);
         if (abs(kjet) > 0.0001)
           enew = e / kjet;
         break;
       }
  
       default:
         cerr << "[Jets] JetPartonCorrector: Error! unknown parametrization type = " << type
              << " No correction applied ..." << endl;
         break;
     }
     return enew;
   }