#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

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)

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 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 14 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 14 of file EnsembleCalibrationLA.h.

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

14 {}

void sistrip::EnsembleCalibrationLA::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 26 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, AlCaHLTBitMon_QueryRunRegistry::string, 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"); 
   BOOST_FOREACH(std::string file, inputFiles) chain->Add((file+inFileLocation).c_str());
 
   int methods = 0;  BOOST_FOREACH(unsigned 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 110 of file EnsembleCalibrationLA.cc.

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

Referenced by analyze(), endJob(), and write_samples_plots().

                            {
   std::fstream file((Prefix+"calibrations.dat").c_str(),std::ios::out);
   std::pair<std::string,MethodCalibrations> cal;
   BOOST_FOREACH(cal,calibrations) {
     file << cal.first << std::endl
      << "\t slopes(";    BOOST_FOREACH(float i, cal.second.slopes) file << i<< ","; file << ")" << std::endl
      << "\t offsets(";   BOOST_FOREACH(float i, cal.second.offsets) file << i<< ","; file << ")" << std::endl
      << "\t pulls(";     BOOST_FOREACH(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 94 of file EnsembleCalibrationLA.cc.

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

Referenced by analyze(), endJob(), and write_ensembles_text().

                                               {
   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 52 of file EnsembleCalibrationLA.cc.

References calibrations, gather_cfg::cout, LA_Filler_Fitter::ensemble_summary(), FrontierConditions_GlobalTag_cff::file, cmsPerfStripChart::format, diffTwoXMLs::label, LA_Filler_Fitter::layer_index(), geometryCSVtoXML::line, LA_Filler_Fitter::offset_slope(), MillePedeFileConverter_cfg::out, Prefix, LA_Filler_Fitter::pull(), AlCaHLTBitMon_QueryRunRegistry::string, heppy_report::summary, sistrip::TIB, and write_ensembles_plots().

Referenced by analyze(), and endJob().

 {
   std::pair<std::string, std::vector<LA_Filler_Fitter::EnsembleSummary> > ensemble;
   BOOST_FOREACH(ensemble, LA_Filler_Fitter::ensemble_summary(book)) {
     std::fstream file((Prefix+ensemble.first+".dat").c_str(),std::ios::out);
     BOOST_FOREACH(LA_Filler_Fitter::EnsembleSummary 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 102 of file EnsembleCalibrationLA.cc.

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

Referenced by analyze(), endJob(), and write_ensembles_plots().

                                             {
   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 40 of file EnsembleCalibrationLA.h.

Referenced by write_calibrations(), and write_ensembles_text().

const double sistrip::EnsembleCalibrationLA::highBin

private

Definition at line 28 of file EnsembleCalibrationLA.h.

Referenced by endJob().

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

private

Definition at line 26 of file EnsembleCalibrationLA.h.

Referenced by endJob().

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

private

Definition at line 25 of file EnsembleCalibrationLA.h.

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

const double sistrip::EnsembleCalibrationLA::lowBin

private

Definition at line 28 of file EnsembleCalibrationLA.h.

Referenced by endJob().

const unsigned sistrip::EnsembleCalibrationLA::maxEvents

private

Definition at line 27 of file EnsembleCalibrationLA.h.

Referenced by endJob().

const unsigned sistrip::EnsembleCalibrationLA::nbins

private

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

Referenced by endJob().

const TrackerTopology* sistrip::EnsembleCalibrationLA::tTopo_

private

Definition at line 42 of file EnsembleCalibrationLA.h.

Referenced by endJob(), and endRun().

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

private

Definition at line 29 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