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
const Definitions &	definitions () const
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

Private Member Functions
template<class Archive >
void	serialize (Archive &ar, const unsigned int version)
	is this a valid set? More...

Private Attributes
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 19 of file JetCorrectorParameters.h.

Constructor & Destructor Documentation

JetCorrectorParameters::JetCorrectorParameters ( )

inline

Definition at line 81 of file JetCorrectorParameters.h.

References AlCaHLTBitMon_QueryRunRegistry::string, and valid_.

{ valid_ = false;}

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

Definition at line 109 of file JetCorrectorParameters.cc.

References trackerTree::check(), definitions(), JetCorrectorParameters::Definitions::formula(), insertMaterial::getSection(), i, input, geometryCSVtoXML::line, mDefinitions, mRecords, JetCorrectorParameters::Definitions::nBinVar(), JetCorrectorParameters::Record::nParameters(), record(), vertexPlots::section, AlCaHLTBitMon_QueryRunRegistry::string, 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))
    {
      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()==""))
            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=="")
    handleError("JetCorrectorParameters","No definitions found!!!");
  if (mRecords.empty() && currentSection == "") mRecords.push_back(Record());
  if (mRecords.empty() && currentSection != "")
    {
      std::stringstream sserr;
      sserr<<"the requested section "<<fSection<<" doesn't exist!";
      handleError("JetCorrectorParameters",sserr.str());
    }
  std::sort(mRecords.begin(), mRecords.end());
  valid_ = true;
}

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

inline

Definition at line 83 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 252 of file JetCorrectorParameters.cc.

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

Referenced by 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 160 of file JetCorrectorParameters.cc.

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

Referenced by SimpleJetCorrector::correction(), JetResolution::parameterEtaEval(), size(), 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;
}

const Definitions& JetCorrectorParameters::definitions ( ) const

inline

Definition at line 88 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;   }

bool JetCorrectorParameters::isValid ( void  ) const

inline

Definition at line 96 of file JetCorrectorParameters.h.

References valid_.

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

{ return valid_; }

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

Definition at line 188 of file JetCorrectorParameters.cc.

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

Referenced by SimpleJetCorrector::correction(), and size().

{
  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 300 of file JetCorrectorParameters.cc.

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

Referenced by JetCorrectorDBReader::analyze(), and size().

{
  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 262 of file JetCorrectorParameters.cc.

References JetCorrectorParameters::Definitions::binVar(), gather_cfg::cout, definitions(), JetCorrectorParameters::Definitions::formula(), i, 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(), and size().

{
  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 87 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

is this a valid set?

unsigned JetCorrectorParameters::size ( void  ) const

inline

Definition at line 89 of file JetCorrectorParameters.h.

References binCenters(), binIndex(), mRecords, neighbourBin(), printFile(), printScreen(), and AlCaHLTBitMon_QueryRunRegistry::string.

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

{return mRecords.size();}

unsigned JetCorrectorParameters::size

(

unsigned

fVar

)

const

Definition at line 230 of file JetCorrectorParameters.cc.

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

Referenced by ntuple._Collection::__iter__(), and ntuple._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 104 of file JetCorrectorParameters.h.

template<typename CondSerializationT , typename Enabled >

friend struct cond::serialization::access

friend

Definition at line 104 of file JetCorrectorParameters.h.

Member Data Documentation

JetCorrectorParameters::Definitions JetCorrectorParameters::mDefinitions

private

Definition at line 100 of file JetCorrectorParameters.h.

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

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

private

Definition at line 101 of file JetCorrectorParameters.h.

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

bool JetCorrectorParameters::valid_

private

Definition at line 102 of file JetCorrectorParameters.h.

Referenced by isValid(), and JetCorrectorParameters().