FactorizedJetCorrectorCalculator Class Reference

#include <FactorizedJetCorrectorCalculator.h>

Classes
class	VariableValues

Public Types
enum	LevelTypes {   kL1, kL2, kL3, kL4,   kL5, kL6, kL7, kL1fj,   kL1JPT }
enum	VarTypes {   kJetPt, kJetEta, kJetPhi, kJetE,   kJetEMF, kRelLepPt, kPtRel, kNPV,   kJetA, kRho, kJPTrawE, kJPTrawEt,   kJPTrawPt, kJPTrawEta, kJPTrawOff }

Public Member Functions
	FactorizedJetCorrectorCalculator ()
	FactorizedJetCorrectorCalculator (const std::string &fLevels, const std::string &fTags, const std::string &fOptions="")
	FactorizedJetCorrectorCalculator (const std::vector< JetCorrectorParameters > &fParameters)
float	getCorrection (VariableValues &) const
std::vector< float >	getSubCorrections (VariableValues &) const
	~FactorizedJetCorrectorCalculator ()

Private Member Functions
void	checkConsistency (const std::vector< std::string > &fLevels, const std::vector< std::string > &fTags)
	FactorizedJetCorrectorCalculator (const FactorizedJetCorrectorCalculator &)
std::vector< float >	fillVector (const std::vector< VarTypes > &fVarTypes, const VariableValues &) const
float	getLepPt (const VariableValues &) const
float	getPtRel (const VariableValues &) const
float	getRelLepPt (const VariableValues &) const
void	initCorrectors (const std::string &fLevels, const std::string &fFiles, const std::string &fOptions)
std::vector< VarTypes >	mapping (const std::vector< std::string > &fNames) const
FactorizedJetCorrectorCalculator &	operator= (const FactorizedJetCorrectorCalculator &)
std::vector< std::string >	parseLevels (const std::string &ss) const
std::string	parseOption (const std::string &ss, const std::string &type) const
std::string	removeSpaces (const std::string &ss) const

Private Attributes
std::vector< std::vector < VarTypes > >	mBinTypes
std::vector < SimpleJetCorrector const * >	mCorrectors
std::vector< LevelTypes >	mLevels
std::vector< std::vector < VarTypes > >	mParTypes

Detailed Description

Definition at line 17 of file FactorizedJetCorrectorCalculator.h.

Member Enumeration Documentation

enum FactorizedJetCorrectorCalculator::LevelTypes

Enumerator
kL1
kL2
kL3
kL4
kL5
kL6
kL7
kL1fj
kL1JPT

Definition at line 79 of file FactorizedJetCorrectorCalculator.h.

{kL1,kL2,kL3,kL4,kL5,kL6,kL7,kL1fj,kL1JPT};

enum FactorizedJetCorrectorCalculator::VarTypes

Enumerator
kJetPt
kJetEta
kJetPhi
kJetE
kJetEMF
kRelLepPt
kPtRel
kNPV
kJetA
kRho
kJPTrawE
kJPTrawEt
kJPTrawPt
kJPTrawEta
kJPTrawOff

Definition at line 78 of file FactorizedJetCorrectorCalculator.h.

{kJetPt,kJetEta,kJetPhi,kJetE,kJetEMF,kRelLepPt,kPtRel,kNPV,kJetA,kRho,kJPTrawE,kJPTrawEt,kJPTrawPt,kJPTrawEta,kJPTrawOff};

Constructor & Destructor Documentation

FactorizedJetCorrectorCalculator::FactorizedJetCorrectorCalculator

(

)

Definition at line 20 of file FactorizedJetCorrectorCalculator.cc.

{
}

FactorizedJetCorrectorCalculator::FactorizedJetCorrectorCalculator	(	const std::string &	fLevels,
		const std::string &	fTags,
		const std::string &	fOptions = ""
	)

Definition at line 26 of file FactorizedJetCorrectorCalculator.cc.

References initCorrectors().

{
  initCorrectors(fLevels, fFiles, fOptions);
}

FactorizedJetCorrectorCalculator::FactorizedJetCorrectorCalculator

(

const std::vector< JetCorrectorParameters > &

fParameters

)

Definition at line 33 of file FactorizedJetCorrectorCalculator.cc.

References i, kL1, kL1fj, kL1JPT, kL2, kL3, kL4, kL5, kL6, kL7, mapping(), mBinTypes, mCorrectors, mLevels, mParTypes, Parameters::parameters, contentValuesCheck::ss, and AlCaHLTBitMon_QueryRunRegistry::string.

{
  for(unsigned i=0;i<fParameters.size();i++) {
    std::string ss = fParameters[i].definitions().level();
    if (ss == "L1Offset")
      mLevels.push_back(kL1);
    else if (ss == "L1JPTOffset")
      mLevels.push_back(kL1JPT);
    else if (ss == "L2Relative")
      mLevels.push_back(kL2);
    else if (ss == "L3Absolute")
      mLevels.push_back(kL3);
    else if (ss == "L4EMF")
      mLevels.push_back(kL4);
    else if (ss == "L5Flavor")
      mLevels.push_back(kL5);
    else if (ss == "L6SLB")
      mLevels.push_back(kL6);
    else if (ss == "L7Parton")
      mLevels.push_back(kL7);
    else if (ss == "L1FastJet")
      mLevels.push_back(kL1fj);
    mCorrectors.push_back(new SimpleJetCorrector(fParameters[i]));
    mBinTypes.push_back(mapping(mCorrectors[i]->parameters().definitions().binVar()));
    mParTypes.push_back(mapping(mCorrectors[i]->parameters().definitions().parVar()));
  }
}

FactorizedJetCorrectorCalculator::~FactorizedJetCorrectorCalculator

(

)

Definition at line 64 of file FactorizedJetCorrectorCalculator.cc.

References i, and mCorrectors.

{
  for(unsigned i=0;i<mCorrectors.size();i++)
    delete mCorrectors[i];
}

FactorizedJetCorrectorCalculator::FactorizedJetCorrectorCalculator ( const FactorizedJetCorrectorCalculator &  )

private

Member Function Documentation

void FactorizedJetCorrectorCalculator::checkConsistency ( const std::vector< std::string > &  fLevels, const std::vector< std::string > &  fTags  )

private

Definition at line 178 of file FactorizedJetCorrectorCalculator.cc.

References i.

Referenced by initCorrectors().

{
  //---- First check: the number of tags must be equal to the number of sub-corrections.
  if (fLevels.size() != fTags.size()) {
    std::stringstream sserr;
    sserr<<"number of correction levels: "<<fLevels.size()<<" doesn't match # of tags: "<<fTags.size();
    handleError("FactorizedJetCorrectorCalculator",sserr.str());
  }
  //---- Second check: each tag must contain the corresponding sub-correction level.
  for(unsigned int i=0;i<fTags.size();i++) {
    if ((int)fTags[i].find(fLevels[i])<0) {
      std::stringstream sserr;
      sserr<<"inconsistent tag: "<<fTags[i]<<" for "<<"the requested correction: "<<fLevels[i];
      handleError("FactorizedJetCorrectorCalculator",sserr.str());
    }
  }
}

std::vector< float > FactorizedJetCorrectorCalculator::fillVector ( const std::vector< VarTypes > &  fVarTypes, const VariableValues &  iValues  ) const

private

Definition at line 308 of file FactorizedJetCorrectorCalculator.cc.

References getPtRel(), getRelLepPt(), i, kJetA, kJetE, kJetEMF, kJetEta, kJetPhi, kJetPt, kJPTrawE, kJPTrawEt, kJPTrawEta, kJPTrawOff, kJPTrawPt, kNPV, kPtRel, kRelLepPt, kRho, FactorizedJetCorrectorCalculator::VariableValues::mIsAddLepToJetset, FactorizedJetCorrectorCalculator::VariableValues::mIsJetAset, FactorizedJetCorrectorCalculator::VariableValues::mIsJetEMFset, FactorizedJetCorrectorCalculator::VariableValues::mIsJetEset, FactorizedJetCorrectorCalculator::VariableValues::mIsJetEtaset, FactorizedJetCorrectorCalculator::VariableValues::mIsJetPhiset, FactorizedJetCorrectorCalculator::VariableValues::mIsJetPtset, FactorizedJetCorrectorCalculator::VariableValues::mIsJPTrawOFFset, FactorizedJetCorrectorCalculator::VariableValues::mIsJPTrawP4set, FactorizedJetCorrectorCalculator::VariableValues::mIsLepPxset, FactorizedJetCorrectorCalculator::VariableValues::mIsLepPyset, FactorizedJetCorrectorCalculator::VariableValues::mIsLepPzset, FactorizedJetCorrectorCalculator::VariableValues::mIsNPVset, FactorizedJetCorrectorCalculator::VariableValues::mIsRhoset, FactorizedJetCorrectorCalculator::VariableValues::mJetA, FactorizedJetCorrectorCalculator::VariableValues::mJetE, FactorizedJetCorrectorCalculator::VariableValues::mJetEMF, FactorizedJetCorrectorCalculator::VariableValues::mJetEta, FactorizedJetCorrectorCalculator::VariableValues::mJetPhi, FactorizedJetCorrectorCalculator::VariableValues::mJetPt, FactorizedJetCorrectorCalculator::VariableValues::mJPTrawE, FactorizedJetCorrectorCalculator::VariableValues::mJPTrawEt, FactorizedJetCorrectorCalculator::VariableValues::mJPTrawEta, FactorizedJetCorrectorCalculator::VariableValues::mJPTrawOff, FactorizedJetCorrectorCalculator::VariableValues::mJPTrawPt, FactorizedJetCorrectorCalculator::VariableValues::mNPV, FactorizedJetCorrectorCalculator::VariableValues::mRho, and query::result.

Referenced by getSubCorrections().

{
//  std::vector<VarTypes> fVarTypes = _fVarTypes;
  std::vector<float> result;
  for(unsigned i=0;i<fVarTypes.size();i++) {
    if (fVarTypes[i] == kJetEta) {
      if (!iValues.mIsJetEtaset)
        handleError("FactorizedJetCorrectorCalculator","jet eta is not set");
      result.push_back(iValues.mJetEta);
    }
    else if (fVarTypes[i] == kNPV) {
      if (!iValues.mIsNPVset)
        handleError("FactorizedJetCorrectorCalculator","number of primary vertices is not set");
      result.push_back(iValues.mNPV);
    }
    else if (fVarTypes[i] == kJetPt) {
      if (!iValues.mIsJetPtset)
        handleError("FactorizedJetCorrectorCalculator","jet pt is not set");
      result.push_back(iValues.mJetPt);
    }
    else if (fVarTypes[i] == kJetPhi) {
      if (!iValues.mIsJetPhiset)
        handleError("FactorizedJetCorrectorCalculator","jet phi is not set");
      result.push_back(iValues.mJetPhi);
    }
    else if (fVarTypes[i] == kJetE) {
      if (!iValues.mIsJetEset)
        handleError("FactorizedJetCorrectorCalculator","jet E is not set");
      result.push_back(iValues.mJetE);
    }
    else if (fVarTypes[i] == kJetEMF) {
      if (!iValues.mIsJetEMFset)
        handleError("FactorizedJetCorrectorCalculator","jet EMF is not set");
      result.push_back(iValues.mJetEMF);
    }
    else if (fVarTypes[i] == kJetA) {
      if (!iValues.mIsJetAset)
        handleError("FactorizedJetCorrectorCalculator","jet area is not set");
      result.push_back(iValues.mJetA);
    }
    else if (fVarTypes[i] == kRho) {
      if (!iValues.mIsRhoset)
        handleError("FactorizedJetCorrectorCalculator","fastjet density Rho is not set");
      result.push_back(iValues.mRho);
    }
    else if (fVarTypes[i] == kJPTrawE) {
      if (!iValues.mIsJPTrawP4set)
        handleError("FactorizedJetCorrectorCalculator","raw CaloJet P4 for JPT is not set");
      result.push_back(iValues.mJPTrawE);
    }
    else if (fVarTypes[i] == kJPTrawEt) {
      if (!iValues.mIsJPTrawP4set)
        handleError("FactorizedJetCorrectorCalculator","raw CaloJet P4 for JPT is not set");
      result.push_back(iValues.mJPTrawEt);
    }
    else if (fVarTypes[i] == kJPTrawPt) {
      if (!iValues.mIsJPTrawP4set)
        handleError("FactorizedJetCorrectorCalculator","raw CaloJet P4 for JPT is not set");
      result.push_back(iValues.mJPTrawPt);
    }
    else if (fVarTypes[i] == kJPTrawEta) {
      if (!iValues.mIsJPTrawP4set)
        handleError("FactorizedJetCorrectorCalculator","raw CaloJet P4 for JPT is not set");
      result.push_back(iValues.mJPTrawEta);
    }
    else if (fVarTypes[i] == kJPTrawOff) {
      if (!iValues.mIsJPTrawOFFset)
        handleError("FactorizedJetCorrectorCalculator","Offset correction for JPT is not set");
      result.push_back(iValues.mJPTrawOff);
    }
    else if (fVarTypes[i] == kRelLepPt) {
      if (!iValues.mIsJetPtset||!iValues.mIsAddLepToJetset||!iValues.mIsLepPxset||!iValues.mIsLepPyset)
        handleError("FactorizedJetCorrectorCalculator","can't calculate rel lepton pt");
      result.push_back(getRelLepPt(iValues));
    }
    else if (fVarTypes[i] == kPtRel) {
      if (!iValues.mIsJetPtset||!iValues.mIsJetEtaset||!iValues.mIsJetPhiset||!iValues.mIsJetEset||!iValues.mIsAddLepToJetset||!iValues.mIsLepPxset||!iValues.mIsLepPyset||!iValues.mIsLepPzset)
        handleError("FactorizedJetCorrectorCalculator","can't calculate ptrel");
      result.push_back(getPtRel(iValues));
    }
    else {
      std::stringstream sserr;
      sserr<<"unknown parameter "<<fVarTypes[i];
      handleError("FactorizedJetCorrectorCalculator",sserr.str());
    }
  }
  return result;
}

float FactorizedJetCorrectorCalculator::getCorrection

(

FactorizedJetCorrectorCalculator::VariableValues &

iValues

)

const

Definition at line 265 of file FactorizedJetCorrectorCalculator.cc.

References getSubCorrections().

Referenced by JetReCalibrator.JetReCalibrator::correct(), L1FastjetCorrector::correction(), L6SLBCorrector::correction(), and FactorizedJetCorrector::getCorrection().

{
  std::vector<float> && vv = getSubCorrections(iValues);
  return vv.back();
}

float FactorizedJetCorrectorCalculator::getLepPt ( const VariableValues &  iValues ) const

private

Definition at line 400 of file FactorizedJetCorrectorCalculator.cc.

References FactorizedJetCorrectorCalculator::VariableValues::mLepPx, FactorizedJetCorrectorCalculator::VariableValues::mLepPy, and mathSSE::sqrt().

Referenced by getRelLepPt(), and getSubCorrections().

{
  return std::sqrt(iValues.mLepPx*iValues.mLepPx + iValues.mLepPy*iValues.mLepPy);
}

float FactorizedJetCorrectorCalculator::getPtRel ( const VariableValues &  iValues ) const

private

Definition at line 415 of file FactorizedJetCorrectorCalculator.cc.

References metsig::jet, RazorAnalyzer::lep2, FactorizedJetCorrectorCalculator::VariableValues::mAddLepToJet, FactorizedJetCorrectorCalculator::VariableValues::mJetE, FactorizedJetCorrectorCalculator::VariableValues::mJetEta, FactorizedJetCorrectorCalculator::VariableValues::mJetPhi, FactorizedJetCorrectorCalculator::VariableValues::mJetPt, FactorizedJetCorrectorCalculator::VariableValues::mLepPx, FactorizedJetCorrectorCalculator::VariableValues::mLepPy, FactorizedJetCorrectorCalculator::VariableValues::mLepPz, and mathSSE::sqrt().

Referenced by fillVector().

{
  typedef ROOT::Math::LorentzVector<ROOT::Math::PtEtaPhiE4D<float> >
    PtEtaPhiELorentzVector;
  typedef ROOT::Math::DisplacementVector3D<ROOT::Math::Cartesian3D<float> >
    XYZVector;
  PtEtaPhiELorentzVector jet;
  XYZVector lep;
  jet.SetPt(iValues.mJetPt);
  jet.SetEta(iValues.mJetEta);
  jet.SetPhi(iValues.mJetPhi);
  jet.SetE(iValues.mJetE);
  lep.SetXYZ(iValues.mLepPx,iValues.mLepPy,iValues.mLepPz);
  float lj_x = (iValues.mAddLepToJet) ? lep.X()+jet.Px() : jet.Px();
  float lj_y = (iValues.mAddLepToJet) ? lep.Y()+jet.Py() : jet.Py();
  float lj_z = (iValues.mAddLepToJet) ? lep.Z()+jet.Pz() : jet.Pz();
  // absolute values squared
  float lj2  = lj_x*lj_x+lj_y*lj_y+lj_z*lj_z;
  if (lj2<=0) {
    std::stringstream sserr;
    sserr<<"lepton+jet momentum sq is not positive: "<<lj2;
    handleError("FactorizedJetCorrectorCalculator",sserr.str());
  }
  float lep2 = lep.X()*lep.X()+lep.Y()*lep.Y()+lep.Z()*lep.Z();
  // projection vec(mu) to lepjet axis
  float lepXlj = lep.X()*lj_x+lep.Y()*lj_y+lep.Z()*lj_z;
  // absolute value squared and normalized
  float pLrel2 = lepXlj*lepXlj/lj2;
  // lep2 = pTrel2 + pLrel2
  float pTrel2 = lep2-pLrel2;
  return (pTrel2 > 0) ? std::sqrt(pTrel2) : 0.0;
}

float FactorizedJetCorrectorCalculator::getRelLepPt ( const VariableValues &  iValues ) const

private

Definition at line 407 of file FactorizedJetCorrectorCalculator.cc.

References getLepPt(), FactorizedJetCorrectorCalculator::VariableValues::mAddLepToJet, and FactorizedJetCorrectorCalculator::VariableValues::mJetPt.

Referenced by fillVector().

{
  float lepPt = getLepPt(iValues);
  return (iValues.mAddLepToJet) ? lepPt/(iValues.mJetPt + lepPt) : lepPt/iValues.mJetPt;
}

std::vector< float > FactorizedJetCorrectorCalculator::getSubCorrections

(

FactorizedJetCorrectorCalculator::VariableValues &

iValues

)

const

Definition at line 273 of file FactorizedJetCorrectorCalculator.cc.

References fillVector(), getLepPt(), i, kL1, kL1fj, kL6, FactorizedJetCorrectorCalculator::VariableValues::mAddLepToJet, mBinTypes, mCorrectors, FactorizedJetCorrectorCalculator::VariableValues::mIsJPTrawOFFset, FactorizedJetCorrectorCalculator::VariableValues::mIsJPTrawP4set, FactorizedJetCorrectorCalculator::VariableValues::mJetE, FactorizedJetCorrectorCalculator::VariableValues::mJetPt, mLevels, mParTypes, FactorizedJetCorrectorCalculator::VariableValues::reset(), pileupReCalc_HLTpaths::scale, and FactorizedJetCorrectorCalculator::VariableValues::setJPTrawOff().

Referenced by getCorrection(), and FactorizedJetCorrector::getSubCorrections().

{
  float scale,factor;
  std::vector<float> factors;
  std::vector<float> vx,vy;
  factor = 1;
  for(unsigned int i=0;i<mLevels.size();i++) {
    vx = fillVector(mBinTypes[i],iValues);
    vy = fillVector(mParTypes[i],iValues);
    //if (mLevels[i]==kL2 || mLevels[i]==kL6)
      //mCorrectors[i]->setInterpolation(true);
    scale = mCorrectors[i]->correction(vx,vy);
    //----- For JPT jets, the offset is stored in order to be used later by the the L1JPTOffset
    if ((mLevels[i]==kL1 || mLevels[i]==kL1fj) && iValues.mIsJPTrawP4set && !iValues.mIsJPTrawOFFset) {
      iValues.setJPTrawOff(scale);
    }
    else if (mLevels[i]==kL6 && iValues.mAddLepToJet) {
      scale  *= 1.0 + getLepPt(iValues) / iValues.mJetPt;
      iValues.mJetE  *= scale;
      iValues.mJetPt *= scale;
      factor *= scale;
    }
    else {
      iValues.mJetE  *= scale;
      iValues.mJetPt *= scale;
      factor *= scale;
    }
    factors.push_back(factor);
  }
  iValues.reset();
  return factors;
}

void FactorizedJetCorrectorCalculator::initCorrectors ( const std::string &  fLevels, const std::string &  fFiles, const std::string &  fOptions  )

private

Definition at line 72 of file FactorizedJetCorrectorCalculator.cc.

References checkConsistency(), i, kL1, kL1fj, kL1JPT, kL2, kL3, kL4, kL5, kL6, kL7, mapping(), mBinTypes, mCorrectors, mLevels, mParTypes, Parameters::parameters, parseLevels(), parseOption(), removeSpaces(), AlCaHLTBitMon_QueryRunRegistry::string, and tmp.

Referenced by FactorizedJetCorrectorCalculator().

{
  //---- Read the CorrectionLevels string and parse the requested sub-correction levels.
  std::vector<std::string> tmp = parseLevels(removeSpaces(fLevels));
  for(unsigned i=0;i<tmp.size();i++) {
    if (tmp[i] == "L1Offset")
      mLevels.push_back(kL1);
    else if (tmp[i] == "L1JPTOffset")
      mLevels.push_back(kL1JPT);
    else if (tmp[i] == "L2Relative")
      mLevels.push_back(kL2);
    else if (tmp[i] == "L3Absolute")
      mLevels.push_back(kL3);
    else if (tmp[i] == "L4EMF")
      mLevels.push_back(kL4);
    else if (tmp[i] == "L5Flavor")
      mLevels.push_back(kL5);
    else if (tmp[i] == "L6SLB")
      mLevels.push_back(kL6);
    else if (tmp[i] == "L7Parton")
      mLevels.push_back(kL7);
    else if (tmp[i] == "L1FastJet")
      mLevels.push_back(kL1fj);
    else {
      std::stringstream sserr;
      sserr<<"unknown correction level "<<tmp[i];
      handleError("FactorizedJetCorrectorCalculator",sserr.str());
    }
  }
  //---- Read the parameter filenames string and parse the requested sub-correction tags.
  std::vector<std::string> Files = parseLevels(removeSpaces(fFiles));
  //---- Read the Options string and define the FlavorOption and PartonOption.
  std::string FlavorOption = parseOption(removeSpaces(fOptions),"L5Flavor");
  std::string PartonOption = parseOption(removeSpaces(fOptions),"L7Parton");
  //---- Check the consistency between tags and requested sub-corrections.
  checkConsistency(tmp,Files);
  //---- Create instances of the requested sub-correctors.
  for(unsigned i=0;i<mLevels.size();i++) {
    if (mLevels[i]==kL1||mLevels[i]==kL1JPT||mLevels[i]==kL2||mLevels[i]==kL3||mLevels[i]==kL4||mLevels[i]==kL6||mLevels[i]==kL1fj)
      mCorrectors.push_back(new SimpleJetCorrector(Files[i]));
    else if (mLevels[i]==kL5 && FlavorOption.length()==0)
      handleError("FactorizedJetCorrectorCalculator","must specify flavor option when requesting L5Flavor correction!");
    else if (mLevels[i]==kL5 && FlavorOption.length()>0)
      mCorrectors.push_back(new SimpleJetCorrector(Files[i],FlavorOption));
    else if (mLevels[i]==kL7 && PartonOption.length()==0)
      handleError("FactorizedJetCorrectorCalculator","must specify parton option when requesting L7Parton correction!");
    else if (mLevels[i]==kL7 && PartonOption.length()>0)
      mCorrectors.push_back(new SimpleJetCorrector(Files[i],PartonOption));
    else {
      std::stringstream sserr;
      sserr<<"unknown correction level "<<tmp[i];
      handleError("FactorizedJetCorrectorCalculator",sserr.str());
    }
    mBinTypes.push_back(mapping(mCorrectors[i]->parameters().definitions().binVar()));
    mParTypes.push_back(mapping(mCorrectors[i]->parameters().definitions().parVar()));
  }
}

std::vector< FactorizedJetCorrectorCalculator::VarTypes > FactorizedJetCorrectorCalculator::mapping ( const std::vector< std::string > &  fNames ) const

private

Definition at line 132 of file FactorizedJetCorrectorCalculator.cc.

References i, kJetA, kJetE, kJetEMF, kJetEta, kJetPhi, kJetPt, kJPTrawE, kJPTrawEt, kJPTrawEta, kJPTrawOff, kJPTrawPt, kNPV, kPtRel, kRelLepPt, kRho, query::result, contentValuesCheck::ss, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by FactorizedJetCorrectorCalculator(), and initCorrectors().

{
  std::vector<VarTypes> result;
  for(unsigned i=0;i<fNames.size();i++) {
    std::string ss = fNames[i];
    if (ss=="JetPt")
      result.push_back(kJetPt);
    else if (ss=="JetEta")
      result.push_back(kJetEta);
    else if (ss=="JetPhi")
      result.push_back(kJetPhi);
    else if (ss=="JetE")
      result.push_back(kJetE);
    else if (ss=="JetEMF")
      result.push_back(kJetEMF);
    else if (ss=="RelLepPt")
      result.push_back(kRelLepPt);
    else if (ss=="PtRel")
      result.push_back(kPtRel);
    else if (ss=="NPV")
      result.push_back(kNPV);
    else if (ss=="JetA")
      result.push_back(kJetA);
    else if (ss=="Rho")
      result.push_back(kRho);
    else if (ss=="JPTrawE")
      result.push_back(kJPTrawE);
    else if (ss=="JPTrawEt")
      result.push_back(kJPTrawEt);
    else if (ss=="JPTrawPt")
      result.push_back(kJPTrawPt);
    else if (ss=="JPTrawEta")
      result.push_back(kJPTrawEta);
    else if (ss=="JPTrawOff")
      result.push_back(kJPTrawOff);
    else {
      std::stringstream sserr;
      sserr<<"unknown parameter name: "<<ss;
      handleError("FactorizedJetCorrectorCalculator",sserr.str());
    }
  }
  return result;
}

FactorizedJetCorrectorCalculator& FactorizedJetCorrectorCalculator::operator= ( const FactorizedJetCorrectorCalculator &  )

private

std::vector< std::string > FactorizedJetCorrectorCalculator::parseLevels ( const std::string &  ss ) const

private

Definition at line 198 of file FactorizedJetCorrectorCalculator.cc.

References i, j, query::result, AlCaHLTBitMon_QueryRunRegistry::string, and tmp.

Referenced by initCorrectors().

{
  std::vector<std::string> result;
  unsigned int pos(0),j,newPos;
  int i;
  std::string tmp;
  //---- The ss string must be of the form: "LX:LY:...:LZ"
  while (pos<ss.length()) {
    tmp = "";
    i = ss.find(":" , pos);
    if (i<0 && pos==0) {
      result.push_back(ss);
      pos = ss.length();
    }
    else if (i<0 && pos>0) {
      for(j=pos;j<ss.length();j++)
        tmp+=ss[j];
      result.push_back(tmp);
      pos = ss.length();
    }
    else {
      newPos = i;
      for(j=pos;j<newPos;j++)
        tmp+=ss[j];
      result.push_back(tmp);
      pos = newPos+1;
    }
  }
  return result;
}

std::string FactorizedJetCorrectorCalculator::parseOption ( const std::string &  ss, const std::string &  type  ) const

private

Definition at line 231 of file FactorizedJetCorrectorCalculator.cc.

References query::result, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by initCorrectors().

{
  std::string result;
  int pos1(-1),pos2(-1);
  //---- The ss string must be of the form: "type1:option1&type2:option2&..."
  pos1 = ss.find(type+":");
  if (pos1<0)
    result = "";
  else {
    pos2 = ss.find("&",pos1+type.length()+1);
    if (pos2<0)
      result = ss.substr(pos1+type.length()+1,ss.length()-pos1-type.length()-1);
    else
      result = ss.substr(pos1+type.length()+1,pos2-pos1-type.length()-1);
  }
  return result;
}

std::string FactorizedJetCorrectorCalculator::removeSpaces ( const std::string &  ss ) const

private

Definition at line 251 of file FactorizedJetCorrectorCalculator.cc.

References i, query::result, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by initCorrectors().

{
  std::string result("");
  std::string aChar;
  for(unsigned int i=0;i<ss.length();i++) {
    aChar = ss.substr(i,1);
    if (aChar != " ")
      result+=aChar;
  }
  return result;
}

Member Data Documentation

std::vector<std::vector<VarTypes> > FactorizedJetCorrectorCalculator::mBinTypes

private

Definition at line 104 of file FactorizedJetCorrectorCalculator.h.

Referenced by FactorizedJetCorrectorCalculator(), getSubCorrections(), and initCorrectors().

std::vector<SimpleJetCorrector const*> FactorizedJetCorrectorCalculator::mCorrectors

private

Definition at line 105 of file FactorizedJetCorrectorCalculator.h.

Referenced by FactorizedJetCorrectorCalculator(), getSubCorrections(), initCorrectors(), and ~FactorizedJetCorrectorCalculator().

std::vector<LevelTypes> FactorizedJetCorrectorCalculator::mLevels

private

Definition at line 103 of file FactorizedJetCorrectorCalculator.h.

Referenced by FactorizedJetCorrectorCalculator(), getSubCorrections(), and initCorrectors().

std::vector<std::vector<VarTypes> > FactorizedJetCorrectorCalculator::mParTypes

private

Definition at line 104 of file FactorizedJetCorrectorCalculator.h.

Referenced by FactorizedJetCorrectorCalculator(), getSubCorrections(), and initCorrectors().