sistrip::EnsembleCalibrationLA Class Reference

#include <EnsembleCalibrationLA.h>

Inheritance diagram for sistrip::EnsembleCalibrationLA:

Classes
struct	MethodCalibrations

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	endJob () override
void	endRun (const edm::Run &, const edm::EventSetup &) override
	EnsembleCalibrationLA (const edm::ParameterSet &)
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	write_calibrations () const
void	write_ensembles_plots (const Book &) const
void	write_ensembles_text (const Book &)
void	write_samples_plots (const Book &) const

Private Attributes
std::map< std::string, MethodCalibrations >	calibrations
const double	highBin
const std::string	inFileLocation
const std::vector< std::string >	inputFiles
const double	lowBin
const unsigned	maxEvents
const unsigned	nbins
const std::string	Prefix
const unsigned	samples
const TrackerTopology *	tTopo_
std::vector< int >	vMethods

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 9 of file EnsembleCalibrationLA.h.

Constructor & Destructor Documentation

EnsembleCalibrationLA()

sistrip::EnsembleCalibrationLA::EnsembleCalibrationLA ( const edm::ParameterSet &  conf )

explicit

Definition at line 13 of file EnsembleCalibrationLA.cc.

      : inputFiles(conf.getParameter<std::vector<std::string> >("InputFiles")),
        inFileLocation(conf.getParameter<std::string>("InFileLocation")),
        Prefix(conf.getUntrackedParameter<std::string>("Prefix", "")),
        maxEvents(conf.getUntrackedParameter<unsigned>("MaxEvents", 0)),
        samples(conf.getParameter<unsigned>("Samples")),
        nbins(conf.getParameter<unsigned>("NBins")),
        lowBin(conf.getParameter<double>("LowBin")),
        highBin(conf.getParameter<double>("HighBin")),
        vMethods(conf.getParameter<std::vector<int> >("Methods")) {}

Member Function Documentation

analyze()

void sistrip::EnsembleCalibrationLA::analyze ( const edm::Event &  , const edm::EventSetup &    )

inlineoverridevirtual

Implements edm::EDAnalyzer.

Definition at line 12 of file EnsembleCalibrationLA.h.

{}

endJob()

void sistrip::EnsembleCalibrationLA::endJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 24 of file EnsembleCalibrationLA.cc.

                                     {
    Book book("la_ensemble");
    TChain* const chain = new TChain("la_ensemble");
    for (auto const& file : inputFiles)
      chain->Add((file + inFileLocation).c_str());
 
    int methods = 0;
    for (unsigned int method : vMethods)
      methods |= method;
 
    LA_Filler_Fitter laff(methods, samples, nbins, lowBin, highBin, maxEvents, tTopo_);
    laff.fill(chain, book);
    laff.fit(book);
    laff.summarize_ensembles(book);
 
    write_ensembles_text(book);
    write_ensembles_plots(book);
    write_samples_plots(book);
    write_calibrations();
  }

References FrontierConditions_GlobalTag_cff::file, LA_Filler_Fitter::fill(), LA_Filler_Fitter::fit(), highBin, inFileLocation, inputFiles, lowBin, maxEvents, AlcaSiPixelAliHarvester0T_cff::method, nbins, samples, LA_Filler_Fitter::summarize_ensembles(), tTopo_, vMethods, write_calibrations(), write_ensembles_plots(), write_ensembles_text(), and write_samples_plots().

Referenced by o2olib.O2ORunMgr::executeJob().

endRun()

void sistrip::EnsembleCalibrationLA::endRun ( const edm::Run &  , const edm::EventSetup &  eSetup  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 45 of file EnsembleCalibrationLA.cc.

                                                                               {
    edm::ESHandle<TrackerTopology> tTopoHandle;
    eSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
    tTopo_ = tTopoHandle.product();
  }

References edm::EventSetup::get(), get, edm::ESHandle< T >::product(), and tTopo_.

write_calibrations()

void sistrip::EnsembleCalibrationLA::write_calibrations ( ) const

private

Definition at line 106 of file EnsembleCalibrationLA.cc.

                                                       {
    std::fstream file((Prefix + "calibrations.dat").c_str(), std::ios::out);
    for (auto const& cal : calibrations) {
      file << cal.first << std::endl << "\t slopes(";
      for (float i : cal.second.slopes)
        file << i << ",";
      file << ")" << std::endl << "\t offsets(";
      for (float i : cal.second.offsets)
        file << i << ",";
      file << ")" << std::endl << "\t pulls(";
      for (float i : cal.second.pulls)
        file << i << ",";
      file << ")" << std::endl;
    }
    file.close();
  }

References calibrations, FrontierConditions_GlobalTag_cff::file, mps_fire::i, MillePedeFileConverter_cfg::out, and Prefix.

Referenced by endJob().

write_ensembles_plots()

void sistrip::EnsembleCalibrationLA::write_ensembles_plots ( const Book &  book ) const

private

Definition at line 90 of file EnsembleCalibrationLA.cc.

                                                                          {
    TFile file((Prefix + "sampleFits.root").c_str(), "RECREATE");
    for (Book::const_iterator hist = book.begin(".*(profile|ratio|reconstruction|symm|symmchi2|_w\\d)");
         hist != book.end();
         ++hist)
      hist->second->Write();
    file.Close();
  }

References Book::begin(), Book::end(), FrontierConditions_GlobalTag_cff::file, compare::hist, and Prefix.

Referenced by endJob().

write_ensembles_text()

void sistrip::EnsembleCalibrationLA::write_ensembles_text ( const Book &  book )

private

Definition at line 51 of file EnsembleCalibrationLA.cc.

                                                                   {
    for (auto const& ensemble : LA_Filler_Fitter::ensemble_summary(book)) {
      std::fstream file((Prefix + ensemble.first + ".dat").c_str(), std::ios::out);
      for (auto const& summary : ensemble.second)
        file << summary << std::endl;
 
      const std::pair<std::pair<float, float>, std::pair<float, float> > line =
          LA_Filler_Fitter::offset_slope(ensemble.second);
      const float pull = LA_Filler_Fitter::pull(ensemble.second);
 
      unsigned index = 15;
      std::string label;
      {
        std::cout << ensemble.first << std::endl;
        boost::regex format(".*(T[IO]B)_layer(\\d)([as])_(.*)");
        if (boost::regex_match(ensemble.first, format)) {
          const bool TIB = "TIB" == boost::regex_replace(ensemble.first, format, "\\1");
          const bool stereo = "s" == boost::regex_replace(ensemble.first, format, "\\3");
          const unsigned layer = boost::lexical_cast<unsigned>(boost::regex_replace(ensemble.first, format, "\\2"));
          label = boost::regex_replace(ensemble.first, format, "\\4");
          index = LA_Filler_Fitter::layer_index(TIB, stereo, layer);
 
          calibrations[label].slopes[index] = line.second.first;
          calibrations[label].offsets[index] = line.first.first;
          calibrations[label].pulls[index] = pull;
        }
      }
 
      file << std::endl
           << std::endl
           << "# Best Fit Line: " << line.first.first << "(" << line.first.second << ")   +   x* " << line.second.first
           << "(" << line.second.second << ")" << std::endl
           << "# Pull (average sigma of (x_measure-x_truth)/e_measure): " << pull << std::endl
           << "LA_Calibration( METHOD_XXXXX , xxx, " << line.second.first << ", " << line.first.first << ", " << pull
           << ")," << std::endl;
      file.close();
    }
  }

References calibrations, gather_cfg::cout, LA_Filler_Fitter::ensemble_summary(), FrontierConditions_GlobalTag_cff::file, dqm-mbProfile::format, label, LA_Filler_Fitter::layer_index(), mps_splice::line, LA_Filler_Fitter::offset_slope(), MillePedeFileConverter_cfg::out, Prefix, LA_Filler_Fitter::pull(), AlCaHLTBitMon_QueryRunRegistry::string, edmLumisInFiles::summary, and sistrip::TIB.

Referenced by endJob().

write_samples_plots()

void sistrip::EnsembleCalibrationLA::write_samples_plots ( const Book &  book ) const

private

Definition at line 99 of file EnsembleCalibrationLA.cc.

                                                                        {
    TFile file((Prefix + "ensembleFits.root").c_str(), "RECREATE");
    for (Book::const_iterator hist = book.begin(".*(measure|merr|ensembleReco|pull)"); hist != book.end(); ++hist)
      hist->second->Write();
    file.Close();
  }

References Book::begin(), Book::end(), FrontierConditions_GlobalTag_cff::file, compare::hist, and Prefix.

Referenced by endJob().

Member Data Documentation

calibrations

std::map<std::string, MethodCalibrations> sistrip::EnsembleCalibrationLA::calibrations

private

Definition at line 35 of file EnsembleCalibrationLA.h.

Referenced by write_calibrations(), and write_ensembles_text().

highBin

const double sistrip::EnsembleCalibrationLA::highBin

private

Definition at line 25 of file EnsembleCalibrationLA.h.

Referenced by endJob().

inFileLocation

const std::string sistrip::EnsembleCalibrationLA::inFileLocation

private

Definition at line 23 of file EnsembleCalibrationLA.h.

Referenced by endJob().

inputFiles

const std::vector<std::string> sistrip::EnsembleCalibrationLA::inputFiles

private

Definition at line 22 of file EnsembleCalibrationLA.h.

Referenced by endJob(), and VarParsing.VarParsing::parseArguments().

lowBin

const double sistrip::EnsembleCalibrationLA::lowBin

private

Definition at line 25 of file EnsembleCalibrationLA.h.

Referenced by endJob().

maxEvents

const unsigned sistrip::EnsembleCalibrationLA::maxEvents

private

Definition at line 24 of file EnsembleCalibrationLA.h.

Referenced by endJob().

nbins

const unsigned sistrip::EnsembleCalibrationLA::nbins

private

Definition at line 24 of file EnsembleCalibrationLA.h.

Referenced by python.rootplot.utilities.Hist::__init_TGraph(), python.rootplot.utilities.Hist::__len__(), python.rootplot.root2matplotlib.Hist::_prepare_xaxis(), python.rootplot.root2matplotlib.Hist::_prepare_yaxis(), python.rootplot.utilities.Hist::av_xerr(), python.rootplot.utilities.Hist::av_yerr(), python.rootplot.utilities.Hist::delete_bin(), python.rootplot.utilities.Hist::divide(), endJob(), python.rootplot.utilities.Hist::TGraph(), and python.rootplot.utilities.Hist::TH1F().

Prefix

const std::string sistrip::EnsembleCalibrationLA::Prefix

private

Definition at line 23 of file EnsembleCalibrationLA.h.

Referenced by write_calibrations(), write_ensembles_plots(), write_ensembles_text(), and write_samples_plots().

samples

const unsigned sistrip::EnsembleCalibrationLA::samples

private

Definition at line 24 of file EnsembleCalibrationLA.h.

Referenced by endJob().

tTopo_

const TrackerTopology* sistrip::EnsembleCalibrationLA::tTopo_

private

Definition at line 37 of file EnsembleCalibrationLA.h.

Referenced by endJob(), and endRun().

vMethods

std::vector<int> sistrip::EnsembleCalibrationLA::vMethods

private

Definition at line 26 of file EnsembleCalibrationLA.h.

Referenced by endJob().