#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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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 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)

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)

Detailed Description

Definition at line 9 of file EnsembleCalibrationLA.h.

Constructor & Destructor Documentation

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

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

inlineoverride

Definition at line 12 of file EnsembleCalibrationLA.h.

References endJob(), endRun(), write_calibrations(), write_ensembles_plots(), write_ensembles_text(), and write_samples_plots().

12 {}

void sistrip::EnsembleCalibrationLA::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 24 of file EnsembleCalibrationLA.cc.

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 analyze(), and 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();
   }

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

override

Definition at line 45 of file EnsembleCalibrationLA.cc.

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

Referenced by analyze().

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

void sistrip::EnsembleCalibrationLA::write_calibrations ( ) const

private

Definition at line 106 of file EnsembleCalibrationLA.cc.

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

Referenced by analyze(), and 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();
   }

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

private

Definition at line 90 of file EnsembleCalibrationLA.cc.

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

Referenced by analyze(), and 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();
   }

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

private

Definition at line 51 of file EnsembleCalibrationLA.cc.

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 analyze(), and 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;
         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();
     }
   }

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

private

Definition at line 99 of file EnsembleCalibrationLA.cc.

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

Referenced by analyze(), and 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

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().

const double sistrip::EnsembleCalibrationLA::highBin

private

Definition at line 25 of file EnsembleCalibrationLA.h.

Referenced by endJob().

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

private

Definition at line 23 of file EnsembleCalibrationLA.h.

Referenced by endJob().

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

private

Definition at line 22 of file EnsembleCalibrationLA.h.

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

const double sistrip::EnsembleCalibrationLA::lowBin

private

Definition at line 25 of file EnsembleCalibrationLA.h.

Referenced by endJob().

const unsigned sistrip::EnsembleCalibrationLA::maxEvents

private

Definition at line 24 of file EnsembleCalibrationLA.h.

Referenced by endJob().

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().

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().

const unsigned sistrip::EnsembleCalibrationLA::samples

private

Definition at line 24 of file EnsembleCalibrationLA.h.

Referenced by endJob().

const TrackerTopology* sistrip::EnsembleCalibrationLA::tTopo_

private

Definition at line 37 of file EnsembleCalibrationLA.h.

Referenced by endJob(), and endRun().

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

private

Definition at line 26 of file EnsembleCalibrationLA.h.

Referenced by endJob().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation