#include <EnsembleCalibrationLA.h>

Inheritance diagram for sistrip::EnsembleCalibrationLA:

Classes
struct	MethodCalibrations

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginRun (const edm::Run &, const edm::EventSetup &) override

void	endJob () override

void	endRun (const edm::Run &, const edm::EventSetup &) override

	EnsembleCalibrationLA (const edm::ParameterSet &)

	~EnsembleCalibrationLA () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDAnalyzerBase () 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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices 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_

edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken_

std::vector< int >	vMethods

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, 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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 10 of file EnsembleCalibrationLA.h.

Constructor & Destructor Documentation

◆ EnsembleCalibrationLA()

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

explicit

Definition at line 11 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")),
         tTopoToken_(esConsumes<edm::Transition::EndRun>()) {}

◆ ~EnsembleCalibrationLA()

sistrip::EnsembleCalibrationLA::~EnsembleCalibrationLA ( )

override

Definition at line 23 of file EnsembleCalibrationLA.cc.

23 {};

Member Function Documentation

◆ analyze()

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

inlineoverridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 13 of file EnsembleCalibrationLA.h.

13 {}

◆ beginRun()

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

override

Definition at line 46 of file EnsembleCalibrationLA.cc.

46 {}

◆ endJob()

void sistrip::EnsembleCalibrationLA::endJob ( void )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 25 of file EnsembleCalibrationLA.cc.

References geometryDiff::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().

                                      {
     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();
   }

◆ endRun()

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

override

Definition at line 48 of file EnsembleCalibrationLA.cc.

References edm::EventSetup::getData(), tTopo_, and tTopoToken_.

                                                                                {
     tTopo_ = &eSetup.getData(tTopoToken_);
   }

◆ write_calibrations()

void sistrip::EnsembleCalibrationLA::write_calibrations ( ) const

private

Definition at line 107 of file EnsembleCalibrationLA.cc.

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

Referenced by endJob().

                                                        {
     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();
   }

◆ write_ensembles_plots()

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

private

Definition at line 91 of file EnsembleCalibrationLA.cc.

References Book::begin(), Book::end(), geometryDiff::file, compareTotals::hist, and Prefix.

Referenced by endJob().

                                                                           {
     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();
   }

◆ write_ensembles_text()

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

private

Definition at line 52 of file EnsembleCalibrationLA.cc.

References calibrations, gather_cfg::cout, LA_Filler_Fitter::ensemble_summary(), geometryDiff::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().

                                                                    {
     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;
         std::regex format(".*(T[IO]B)_layer(\\d)([as])_(.*)");
         if (std::regex_match(ensemble.first, format)) {
           const bool TIB = "TIB" == std::regex_replace(ensemble.first, format, "\\1");
           const bool stereo = "s" == std::regex_replace(ensemble.first, format, "\\3");
           const unsigned layer = std::stoul(std::regex_replace(ensemble.first, format, "\\2"));
           label = std::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();
     }
   }

◆ write_samples_plots()

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

private

Definition at line 100 of file EnsembleCalibrationLA.cc.

References Book::begin(), Book::end(), geometryDiff::file, compareTotals::hist, and Prefix.

Referenced by endJob().

                                                                         {
     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();
   }

Member Data Documentation

◆ calibrations

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

private

Definition at line 39 of file EnsembleCalibrationLA.h.

Referenced by write_calibrations(), and write_ensembles_text().

◆ highBin

const double sistrip::EnsembleCalibrationLA::highBin

private

Definition at line 29 of file EnsembleCalibrationLA.h.

Referenced by endJob().

◆ inFileLocation

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

private

Definition at line 27 of file EnsembleCalibrationLA.h.

Referenced by endJob().

◆ inputFiles

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

private

Definition at line 26 of file EnsembleCalibrationLA.h.

Referenced by endJob(), VarParsing.VarParsing::parseArguments(), and postprocessor.PostProcessor::run().

◆ lowBin

const double sistrip::EnsembleCalibrationLA::lowBin

private

Definition at line 29 of file EnsembleCalibrationLA.h.

Referenced by endJob().

◆ maxEvents

const unsigned sistrip::EnsembleCalibrationLA::maxEvents

private

Definition at line 28 of file EnsembleCalibrationLA.h.

Referenced by endJob().

◆ nbins

const unsigned sistrip::EnsembleCalibrationLA::nbins

private

Definition at line 28 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 27 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 28 of file EnsembleCalibrationLA.h.

Referenced by endJob().

◆ tTopo_

const TrackerTopology* sistrip::EnsembleCalibrationLA::tTopo_

private

Definition at line 42 of file EnsembleCalibrationLA.h.

Referenced by endJob(), and endRun().

◆ tTopoToken_

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> sistrip::EnsembleCalibrationLA::tTopoToken_

private

Definition at line 41 of file EnsembleCalibrationLA.h.

Referenced by endRun().

◆ vMethods

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

private

Definition at line 30 of file EnsembleCalibrationLA.h.

Referenced by endJob().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EnsembleCalibrationLA()

◆ ~EnsembleCalibrationLA()

Member Function Documentation

◆ analyze()

◆ beginRun()

◆ endJob()

◆ endRun()

◆ write_calibrations()

◆ write_ensembles_plots()

◆ write_ensembles_text()

◆ write_samples_plots()

Member Data Documentation

◆ calibrations

◆ highBin

◆ inFileLocation

◆ inputFiles

◆ lowBin

◆ maxEvents

◆ nbins

◆ Prefix

◆ samples

◆ tTopo_

◆ tTopoToken_

◆ vMethods