#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 final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const 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

bool	wantsStreamRuns () const

	~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(), and VarParsing.VarParsing::parseArguments().

◆ 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