JetCorrectorParameters Class Reference

#include <JetCorrectorParameters.h>

Classes
class	Definitions
class	Record

Public Member Functions
std::vector< float >	binCenters (unsigned fVar) const
int	binIndex (const std::vector< float > &fX) const
int	binIndexN (const std::vector< float > &fX) const
const Definitions &	definitions () const
void	init ()
bool	isValid () const
	JetCorrectorParameters ()
	JetCorrectorParameters (const std::string &fFile, const std::string &fSection="")
	JetCorrectorParameters (const JetCorrectorParameters::Definitions &fDefinitions, const std::vector< JetCorrectorParameters::Record > &fRecords)
int	neighbourBin (unsigned fIndex, unsigned fVar, bool fNext) const
void	printFile (const std::string &fFileName) const
void	printScreen () const
const Record &	record (unsigned fBin) const
unsigned	size () const
unsigned	size (unsigned fVar) const

Static Public Attributes
static const int	MAX_SIZE_DIMENSIONALITY = 3

Private Member Functions
template<class Archive >
void	serialize (Archive &ar, const unsigned int version)

Private Attributes
std::shared_ptr < JetCorrectorParametersHelper >	helper
	is this a valid set? More...
JetCorrectorParameters::Definitions	mDefinitions
std::vector < JetCorrectorParameters::Record >	mRecords
bool	valid_

Friends
class	boost::serialization::access
template<typename CondSerializationT , typename Enabled >
struct	cond::serialization::access

Detailed Description

Definition at line 25 of file JetCorrectorParameters.h.

Constructor & Destructor Documentation

JetCorrectorParameters::JetCorrectorParameters ( )

inline

Definition at line 113 of file JetCorrectorParameters.h.

References valid_.

{ valid_ = false; }

JetCorrectorParameters::JetCorrectorParameters	(	const std::string &	fFile,
		const std::string &	fSection = ""
	)

Definition at line 113 of file JetCorrectorParameters.cc.

References trackerTree::check(), definitions(), JetCorrectorParameters::Definitions::formula(), mps_fire::i, init(), input, geometryCSVtoXML::line, MAX_SIZE_DIMENSIONALITY, mDefinitions, mRecords, JetCorrectorParameters::Definitions::nBinVar(), JetCorrectorParameters::Record::nParameters(), record(), vertexPlots::section, AlCaHLTBitMon_QueryRunRegistry::string, createJobs::tmp, valid_, JetCorrectorParameters::Record::xMax(), and JetCorrectorParameters::Record::xMin().

                                                                                                {
  std::ifstream input(fFile.c_str());
  std::string currentSection = "";
  std::string line;
  std::string currentDefinitions = "";
  while (std::getline(input, line)) {
    if (line.empty())
      continue;
    std::string section = getSection(line);
    std::string tmp = getDefinitions(line);
    if (!section.empty() && tmp.empty()) {
      currentSection = section;
      continue;
    }
    if (currentSection == fSection) {
      if (!tmp.empty()) {
        currentDefinitions = tmp;
        continue;
      }
      Definitions definitions(currentDefinitions);
      if (!(definitions.nBinVar() == 0 && definitions.formula().empty()))
        mDefinitions = definitions;
      Record record(line, mDefinitions.nBinVar());
      bool check(true);
      for (unsigned i = 0; i < mDefinitions.nBinVar(); ++i)
        if (record.xMin(i) == 0 && record.xMax(i) == 0)
          check = false;
      if (record.nParameters() == 0)
        check = false;
      if (check)
        mRecords.push_back(record);
    }
  }
  if (currentDefinitions.empty())
    handleError("JetCorrectorParameters", "No definitions found!!!");
  if (mRecords.empty() && currentSection.empty())
    mRecords.push_back(Record());
  if (mRecords.empty() && !currentSection.empty()) {
    std::stringstream sserr;
    sserr << "the requested section " << fSection << " doesn't exist!";
    handleError("JetCorrectorParameters", sserr.str());
  }
  std::sort(mRecords.begin(), mRecords.end());
  valid_ = true;

  if (mDefinitions.nBinVar() <= MAX_SIZE_DIMENSIONALITY) {
    init();
  } else {
    std::stringstream sserr;
    sserr << "since the binned dimensionality is greater than " << MAX_SIZE_DIMENSIONALITY
          << " the SimpleJetCorrector will default to using the legacy binIndex function!";
    handleError("JetCorrectorParameters", sserr.str());
  }
}

JetCorrectorParameters::JetCorrectorParameters ( const JetCorrectorParameters::Definitions &  fDefinitions, const std::vector< JetCorrectorParameters::Record > &  fRecords  )

inline

Definition at line 115 of file JetCorrectorParameters.h.

References valid_.

      : mDefinitions(fDefinitions), mRecords(fRecords) {
    valid_ = true;
  }

Member Function Documentation

std::vector< float > JetCorrectorParameters::binCenters

(

unsigned

fVar

)

const

Definition at line 269 of file JetCorrectorParameters.cc.

References mps_fire::i, record(), mps_fire::result, and size().

                                                                       {
  std::vector<float> result;
  for (unsigned i = 0; i < size(); ++i)
    result.push_back(record(i).xMiddle(fVar));
  return result;
}

int JetCorrectorParameters::binIndex

(

const std::vector< float > &

fX

)

const

Definition at line 178 of file JetCorrectorParameters.cc.

References mps_fire::i, dqmiolumiharvest::j, mDefinitions, N, JetCorrectorParameters::Definitions::nBinVar(), record(), mps_fire::result, size(), createJobs::tmp, JetCorrectorParameters::Record::xMax(), and JetCorrectorParameters::Record::xMin().

Referenced by SimpleJetCorrector::correction(), JetResolution::parameterEtaEval(), and SimpleJetCorrectionUncertainty::uncertainty().

                                                                     {
  int result = -1;
  unsigned N = mDefinitions.nBinVar();
  if (N != fX.size()) {
    std::stringstream sserr;
    sserr << "# bin variables " << N << " doesn't correspont to requested #: " << fX.size();
    handleError("JetCorrectorParameters", sserr.str());
  }
  unsigned tmp;
  for (unsigned i = 0; i < size(); ++i) {
    tmp = 0;
    for (unsigned j = 0; j < N; j++)
      if (fX[j] >= record(i).xMin(j) && fX[j] < record(i).xMax(j))
        tmp += 1;
    if (tmp == N) {
      result = i;
      break;
    }
  }
  return result;
}

int JetCorrectorParameters::binIndexN

(

const std::vector< float > &

fX

)

const

Definition at line 202 of file JetCorrectorParameters.cc.

References helper, and mRecords.

Referenced by SimpleJetCorrector::correction().

                                                                      {
  if (helper->size() > 0) {
    return helper->binIndexN(fX, mRecords);
  } else {
    return -1;
  }
}

const Definitions& JetCorrectorParameters::definitions ( ) const

inline

Definition at line 122 of file JetCorrectorParameters.h.

References mDefinitions.

Referenced by SimpleJetCorrector::correction(), SimpleJetCorrector::correctionBin(), SimpleJetCorrector::findInvertVar(), JetCorrectionUncertainty::getUncertainty(), JetResolution::initialize(), JetCorrectorParameters(), L1FastjetCorrector::L1FastjetCorrector(), L1JPTOffsetCorrector::L1JPTOffsetCorrector(), L1OffsetCorrector::L1OffsetCorrector(), LXXXCorrector::LXXXCorrector(), JetResolution::parameterEtaEval(), printFile(), printScreen(), and SimpleJetCorrector::SimpleJetCorrector().

{ return mDefinitions; }

void JetCorrectorParameters::init

(

void

)

Definition at line 170 of file JetCorrectorParameters.cc.

References helper, mDefinitions, and mRecords.

Referenced by JetCorrectorParameters().

                                  {
  std::sort(mRecords.begin(), mRecords.end());
  helper = std::make_shared<JetCorrectorParametersHelper>();
  helper->init(mDefinitions, mRecords);
}

bool JetCorrectorParameters::isValid ( void  ) const

inline

Definition at line 131 of file JetCorrectorParameters.h.

References valid_.

Referenced by ntupleDataFormat._Object::_checkIsValid(), and core.AutoHandle.AutoHandle::ReallyLoad().

{ return valid_; }

int JetCorrectorParameters::neighbourBin	(	unsigned	fIndex,
		unsigned	fVar,
		bool	fNext
	)		const

Definition at line 212 of file JetCorrectorParameters.cc.

References mps_fire::i, dqmiolumiharvest::j, mDefinitions, N, JetCorrectorParameters::Definitions::nBinVar(), record(), mps_fire::result, size(), and createJobs::tmp.

Referenced by SimpleJetCorrector::correction().

                                                                                         {
  int result = -1;
  unsigned N = mDefinitions.nBinVar();
  if (fVar >= N) {
    std::stringstream sserr;
    sserr << "# of bin variables " << N << " doesn't correspond to requested #: " << fVar;
    handleError("JetCorrectorParameters", sserr.str());
  }
  unsigned tmp;
  for (unsigned i = 0; i < size(); ++i) {
    tmp = 0;
    for (unsigned j = 0; j < fVar; j++)
      if (fabs(record(i).xMin(j) - record(fIndex).xMin(j)) < 0.0001)
        tmp += 1;
    for (unsigned j = fVar + 1; j < N; j++)
      if (fabs(record(i).xMin(j) - record(fIndex).xMin(j)) < 0.0001)
        tmp += 1;
    if (tmp < N - 1)
      continue;
    if (tmp == N - 1) {
      if (fNext)
        if (fabs(record(i).xMin(fVar) - record(fIndex).xMax(fVar)) < 0.0001)
          tmp += 1;
      if (!fNext)
        if (fabs(record(i).xMax(fVar) - record(fIndex).xMin(fVar)) < 0.0001)
          tmp += 1;
    }
    if (tmp == N) {
      result = i;
      break;
    }
  }
  return result;
}

void JetCorrectorParameters::printFile

(

const std::string &

fFileName

)

const

Definition at line 316 of file JetCorrectorParameters.cc.

References JetCorrectorParameters::Definitions::binVar(), definitions(), EcalElecEmulTccOnlyExample_cfg::formula, JetCorrectorParameters::Definitions::formula(), mps_fire::i, dqmiolumiharvest::j, JetCorrectorParameters::Definitions::level(), JetCorrectorParameters::Definitions::nBinVar(), JetCorrectorParameters::Record::nParameters(), JetCorrectorParameters::Definitions::nParVar(), JetCorrectorParameters::Record::parameter(), JetCorrectorParameters::Definitions::parVar(), record(), size(), JetCorrectorParameters::Record::xMax(), and JetCorrectorParameters::Record::xMin().

Referenced by JetCorrectorDBReader::analyze().

                                                                       {
  std::ofstream txtFile;
  txtFile.open(fFileName.c_str());
  txtFile.setf(std::ios::right);
  txtFile << "{" << definitions().nBinVar() << std::setw(15);
  for (unsigned i = 0; i < definitions().nBinVar(); i++)
    txtFile << definitions().binVar(i) << std::setw(15);
  txtFile << definitions().nParVar() << std::setw(15);
  for (unsigned i = 0; i < definitions().nParVar(); i++)
    txtFile << definitions().parVar(i) << std::setw(15);
  txtFile << std::setw(definitions().formula().size() + 15) << definitions().formula() << std::setw(15);
  if (definitions().isResponse())
    txtFile << "Response" << std::setw(15);
  else
    txtFile << "Correction" << std::setw(15);
  txtFile << definitions().level() << "}"
          << "\n";
  for (unsigned i = 0; i < size(); i++) {
    for (unsigned j = 0; j < definitions().nBinVar(); j++)
      txtFile << record(i).xMin(j) << std::setw(15) << record(i).xMax(j) << std::setw(15);
    txtFile << record(i).nParameters() << std::setw(15);
    for (unsigned j = 0; j < record(i).nParameters(); j++)
      txtFile << record(i).parameter(j) << std::setw(15);
    txtFile << "\n";
  }
  txtFile.close();
}

void JetCorrectorParameters::printScreen

(

)

const

Definition at line 278 of file JetCorrectorParameters.cc.

References JetCorrectorParameters::Definitions::binVar(), gather_cfg::cout, definitions(), JetCorrectorParameters::Definitions::formula(), mps_fire::i, dqmiolumiharvest::j, JetCorrectorParameters::Definitions::level(), JetCorrectorParameters::Definitions::nBinVar(), JetCorrectorParameters::Record::nParameters(), JetCorrectorParameters::Definitions::nParVar(), JetCorrectorParameters::Record::parameter(), JetCorrectorParameters::Definitions::parVar(), record(), size(), JetCorrectorParameters::Record::xMax(), and JetCorrectorParameters::Record::xMin().

Referenced by JetCorrectorDBReader::analyze().

                                               {
  std::cout << "--------------------------------------------" << std::endl;
  std::cout << "////////  PARAMETERS: //////////////////////" << std::endl;
  std::cout << "--------------------------------------------" << std::endl;
  std::cout << "Number of binning variables:   " << definitions().nBinVar() << std::endl;
  std::cout << "Names of binning variables:    ";
  for (unsigned i = 0; i < definitions().nBinVar(); i++)
    std::cout << definitions().binVar(i) << " ";
  std::cout << std::endl;
  std::cout << "--------------------------------------------" << std::endl;
  std::cout << "Number of parameter variables: " << definitions().nParVar() << std::endl;
  std::cout << "Names of parameter variables:  ";
  for (unsigned i = 0; i < definitions().nParVar(); i++)
    std::cout << definitions().parVar(i) << " ";
  std::cout << std::endl;
  std::cout << "--------------------------------------------" << std::endl;
  std::cout << "Parametrization Formula:       " << definitions().formula() << std::endl;
  if (definitions().isResponse())
    std::cout << "Type (Response or Correction): "
              << "Response" << std::endl;
  else
    std::cout << "Type (Response or Correction): "
              << "Correction" << std::endl;
  std::cout << "Correction Level:              " << definitions().level() << std::endl;
  std::cout << "--------------------------------------------" << std::endl;
  std::cout << "------- Bin contents -----------------------" << std::endl;
  for (unsigned i = 0; i < size(); i++) {
    for (unsigned j = 0; j < definitions().nBinVar(); j++)
      std::cout << record(i).xMin(j) << " " << record(i).xMax(j) << " ";
    std::cout << record(i).nParameters() << " ";
    for (unsigned j = 0; j < record(i).nParameters(); j++)
      std::cout << record(i).parameter(j) << " ";
    std::cout << std::endl;
  }
}

const Record& JetCorrectorParameters::record ( unsigned  fBin ) const

inline

Definition at line 121 of file JetCorrectorParameters.h.

References mRecords.

Referenced by binCenters(), binIndex(), SimpleJetCorrector::correction(), SimpleJetCorrector::correctionBin(), JetCorrectorParameters(), neighbourBin(), JetResolution::parameterEtaEval(), printFile(), printScreen(), size(), and SimpleJetCorrectionUncertainty::uncertaintyBin().

{ return mRecords[fBin]; }

template<class Archive >

void JetCorrectorParameters::serialize ( Archive &  ar, const unsigned int  version  )

private

unsigned JetCorrectorParameters::size ( void  ) const

inline

Definition at line 123 of file JetCorrectorParameters.h.

References mRecords.

Referenced by ntupleDataFormat._Collection::__iter__(), ntupleDataFormat._Collection::__len__(), binCenters(), binIndex(), SimpleJetCorrector::correctionBin(), neighbourBin(), JetResolution::parameterEtaEval(), printFile(), printScreen(), size(), and SimpleJetCorrectionUncertainty::uncertaintyBin().

{ return mRecords.size(); }

unsigned JetCorrectorParameters::size

(

unsigned

fVar

)

const

Definition at line 249 of file JetCorrectorParameters.cc.

References mps_fire::i, mDefinitions, JetCorrectorParameters::Definitions::nBinVar(), record(), mps_fire::result, size(), JetCorrectorParameters::Record::xMax(), and JetCorrectorParameters::Record::xMin().

Referenced by ntupleDataFormat._Collection::__iter__(), and ntupleDataFormat._Collection::__len__().

                                                         {
  if (fVar >= mDefinitions.nBinVar()) {
    std::stringstream sserr;
    sserr << "requested bin variable index " << fVar << " is greater than number of variables "
          << mDefinitions.nBinVar();
    handleError("JetCorrectorParameters", sserr.str());
  }
  unsigned result = 0;
  float tmpMin(-9999), tmpMax(-9999);
  for (unsigned i = 0; i < size(); ++i)
    if (record(i).xMin(fVar) > tmpMin && record(i).xMax(fVar) > tmpMax) {
      result++;
      tmpMin = record(i).xMin(fVar);
      tmpMax = record(i).xMax(fVar);
    }
  return result;
}

Friends And Related Function Documentation

friend class boost::serialization::access

friend

Definition at line 144 of file JetCorrectorParameters.h.

template<typename CondSerializationT , typename Enabled >

friend struct cond::serialization::access

friend

Definition at line 144 of file JetCorrectorParameters.h.

Member Data Documentation

std::shared_ptr<JetCorrectorParametersHelper> JetCorrectorParameters::helper

private

is this a valid set?

Definition at line 142 of file JetCorrectorParameters.h.

Referenced by binIndexN(), and init().

const int JetCorrectorParameters::MAX_SIZE_DIMENSIONALITY = 3

static

Definition at line 134 of file JetCorrectorParameters.h.

Referenced by JetCorrectorParametersHelper::binIndexN(), JetCorrectorParametersHelper::init(), JetCorrectorParametersHelper::initTransientMaps(), and JetCorrectorParameters().

JetCorrectorParameters::Definitions JetCorrectorParameters::mDefinitions

private

Definition at line 138 of file JetCorrectorParameters.h.

Referenced by binIndex(), definitions(), init(), JetCorrectorParameters(), neighbourBin(), and size().

std::vector<JetCorrectorParameters::Record> JetCorrectorParameters::mRecords

private

Definition at line 139 of file JetCorrectorParameters.h.

Referenced by binIndexN(), init(), JetCorrectorParameters(), record(), and size().

bool JetCorrectorParameters::valid_

private

Definition at line 140 of file JetCorrectorParameters.h.

Referenced by isValid(), and JetCorrectorParameters().