sistrip::MeasureLA Class Reference

#include <MeasureLA.h>

Inheritance diagram for sistrip::MeasureLA:

Public Member Functions
	MeasureLA (const edm::ParameterSet &)
std::unique_ptr< SiStripLorentzAngle >	produce (const SiStripLorentzAngleRcd &)
Public Member Functions inherited from edm::ESProducer
	ESProducer ()
	ESProducer (const ESProducer &)=delete
ESProxyIndex const *	getTokenIndices (unsigned int iIndex) const
ESProducer const &	operator= (const ESProducer &)=delete
template<typename Record >
void	updateFromMayConsumes (unsigned int iIndex, const Record &iRecord)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &) final
	~ESProducer ()(false) override
Public Member Functions inherited from edm::ESProxyFactoryProducer
	ESProxyFactoryProducer ()
	ESProxyFactoryProducer (const ESProxyFactoryProducer &)=delete
const ESProxyFactoryProducer &	operator= (const ESProxyFactoryProducer &)=delete
	~ESProxyFactoryProducer () noexcept(false) override
Public Member Functions inherited from edm::eventsetup::DataProxyProvider
void	createKeyedProxies (EventSetupRecordKey const &key, unsigned int nConcurrentIOVs)
	DataProxyProvider ()
	DataProxyProvider (const DataProxyProvider &)=delete
const ComponentDescription &	description () const
void	fillRecordsNotAllowingConcurrentIOVs (std::set< EventSetupRecordKey > &recordsNotAllowingConcurrentIOVs) const
virtual void	initConcurrentIOVs (EventSetupRecordKey const &key, unsigned int nConcurrentIOVs)
bool	isUsingRecord (const EventSetupRecordKey &key) const
KeyedProxies &	keyedProxies (const EventSetupRecordKey &iRecordKey, unsigned int iovIndex=0)
const DataProxyProvider &	operator= (const DataProxyProvider &)=delete
void	setAppendToDataLabel (const edm::ParameterSet &)
void	setDescription (const ComponentDescription &iDescription)
std::set< EventSetupRecordKey >	usingRecords () const
virtual	~DataProxyProvider () noexcept(false)

Private Types
enum	GRANULARITY { LAYER = 0, MODULE = 1, MODULESUMMARY = 2 }

Private Member Functions
void	calibrate (const std::pair< unsigned, LA_Filler_Fitter::Method >, LA_Filler_Fitter::Result &) const
std::pair< unsigned, LA_Filler_Fitter::Method >	calibration_key (const std::string layer, const LA_Filler_Fitter::Method) const
std::pair< unsigned, LA_Filler_Fitter::Method >	calibration_key (const uint32_t detid, const LA_Filler_Fitter::Method) const
void	process_reports () const
void	store_calibrations ()
void	store_methods_and_granularity (const edm::VParameterSet &)
void	summarize_module_muH_byLayer (const LA_Filler_Fitter &)
void	write_report_plots (const std::string, const LA_Filler_Fitter::Method, const GRANULARITY) const
template<class T >
void	write_report_text (const std::string, const LA_Filler_Fitter::Method &, const std::map< T, LA_Filler_Fitter::Result > &) const
void	write_report_text_ms (const std::string, const LA_Filler_Fitter::Method) const

Static Private Member Functions
static std::string	granularity (int32_t g)

Private Attributes
Book	book
bool	byLayer
bool	byModule
const edm::VParameterSet	calibrations
std::map< std::pair< uint32_t, LA_Filler_Fitter::Method >, float >	error_scaling
const edm::FileInPath	fp_
const std::string	inFileLocation
const std::vector< std::string >	inputFiles
const float	localybin
const unsigned	maxEvents
const edm::VParameterSet	measurementPreferences
int32_t	methods
std::map< std::pair< uint32_t, LA_Filler_Fitter::Method >, float >	offset
const edm::VParameterSet	reports
std::map< std::pair< uint32_t, LA_Filler_Fitter::Method >, float >	slope
const unsigned	stripsperbin
TrackerTopology	tTopo_

Additional Inherited Members
Static Public Member Functions inherited from edm::eventsetup::DataProxyProvider
static void	prevalidate (ConfigurationDescriptions &)
Protected Types inherited from edm::ESProxyFactoryProducer
using	EventSetupRecordKey = eventsetup::EventSetupRecordKey
Protected Types inherited from edm::eventsetup::DataProxyProvider
using	KeyedProxiesVector = std::vector< std::pair< DataKey, std::shared_ptr< DataProxy >>>
Protected Member Functions inherited from edm::ESProducer
template<typename T >
auto	setWhatProduced (T *iThis, const es::Label &iLabel={})
template<typename T >
auto	setWhatProduced (T *iThis, const char *iLabel)
template<typename T >
auto	setWhatProduced (T *iThis, const std::string &iLabel)
template<typename T , typename TDecorator >
auto	setWhatProduced (T *iThis, const TDecorator &iDec, const es::Label &iLabel={})
template<typename T , typename TReturn , typename TRecord >
auto	setWhatProduced (T *iThis, TReturn(T::*iMethod)(const TRecord &), const es::Label &iLabel={})
template<typename T , typename TReturn , typename TRecord , typename TArg >
ESConsumesCollectorT< TRecord >	setWhatProduced (T *iThis, TReturn(T::*iMethod)(const TRecord &), const TArg &iDec, const es::Label &iLabel={})
Protected Member Functions inherited from edm::ESProxyFactoryProducer
template<class TFactory >
void	registerFactory (std::unique_ptr< TFactory > iFactory, const std::string &iLabel=std::string())
virtual void	registerFactoryWithKey (const EventSetupRecordKey &iRecord, std::unique_ptr< eventsetup::ProxyFactoryBase > iFactory, const std::string &iLabel=std::string())
KeyedProxiesVector	registerProxies (const EventSetupRecordKey &, unsigned int iovIndex) override
Protected Member Functions inherited from edm::eventsetup::DataProxyProvider
template<class T >
void	usingRecord ()
void	usingRecordWithKey (const EventSetupRecordKey &key)

Detailed Description

Definition at line 17 of file MeasureLA.h.

Member Enumeration Documentation

enum sistrip::MeasureLA::GRANULARITY

private

Enumerator
LAYER
MODULE
MODULESUMMARY

Definition at line 23 of file MeasureLA.h.

{ LAYER = 0, MODULE = 1, MODULESUMMARY = 2 };

Constructor & Destructor Documentation

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

explicit

Definition at line 20 of file MeasureLA.cc.

References book, byLayer, byModule, FrontierConditions_GlobalTag_cff::file, LA_Filler_Fitter::fill(), LA_Filler_Fitter::fit(), inFileLocation, inputFiles, localybin, maxEvents, measurementPreferences, methods, process_reports(), produce(), reports, edm::ESProducer::setWhatProduced(), store_calibrations(), store_methods_and_granularity(), stripsperbin, summarize_module_muH_byLayer(), and tTopo_.

      : inputFiles(conf.getParameter<std::vector<std::string> >("InputFiles")),
        inFileLocation(conf.getParameter<std::string>("InFileLocation")),
        fp_(conf.getParameter<edm::FileInPath>("SiStripDetInfo")),
        reports(conf.getParameter<edm::VParameterSet>("Reports")),
        measurementPreferences(conf.getParameter<edm::VParameterSet>("MeasurementPreferences")),
        calibrations(conf.getParameter<edm::VParameterSet>("Calibrations")),
        methods(0),
        byModule(false),
        byLayer(false),
        localybin(conf.getUntrackedParameter<double>("LocalYBin", 0.0)),
        stripsperbin(conf.getUntrackedParameter<unsigned>("StripsPerBin", 0)),
        maxEvents(conf.getUntrackedParameter<unsigned>("MaxEvents", 0)),
        tTopo_(StandaloneTrackerTopology::fromTrackerParametersXMLFile(
            conf.getParameter<edm::FileInPath>("TrackerParameters").fullPath())) {
    store_methods_and_granularity(reports);
    store_methods_and_granularity(measurementPreferences);
    store_calibrations();

    TChain* const chain = new TChain("la_data");
    for (auto const& file : inputFiles)
      chain->Add((file + inFileLocation).c_str());

    LA_Filler_Fitter laff(methods, byLayer, byModule, localybin, stripsperbin, maxEvents, &tTopo_);
    laff.fill(chain, book);
    laff.fit(book);
    summarize_module_muH_byLayer(laff);
    process_reports();

    setWhatProduced(this, &MeasureLA::produce);
  }

Member Function Documentation

void sistrip::MeasureLA::calibrate ( const std::pair< unsigned, LA_Filler_Fitter::Method >  key, LA_Filler_Fitter::Result &  result  ) const

inlineprivate

Definition at line 173 of file MeasureLA.cc.

References LA_Filler_Fitter::Result::calMeasured, error_scaling, LA_Filler_Fitter::Result::measured, offset, and slope.

Referenced by granularity(), summarize_module_muH_byLayer(), and write_report_text().

                                                                         {
    result.calMeasured = std::make_pair<float, float>(
        (result.measured.first - offset.find(key)->second) / slope.find(key)->second,
        result.measured.second * error_scaling.find(key)->second / slope.find(key)->second);
  }

std::pair< uint32_t, LA_Filler_Fitter::Method > sistrip::MeasureLA::calibration_key ( const std::string  layer, const LA_Filler_Fitter::Method  method  ) const

private

Definition at line 180 of file MeasureLA.cc.

References dqm-mbProfile::format, and LA_Filler_Fitter::layer_index().

Referenced by granularity(), summarize_module_muH_byLayer(), and write_report_text().

                                                                        {
    boost::regex format(".*(T[IO]B)_layer(\\d)([as]).*");
    const bool isTIB = "TIB" == boost::regex_replace(layer, format, "\\1");
    const bool stereo = "s" == boost::regex_replace(layer, format, "\\3");
    const unsigned layerNum = boost::lexical_cast<unsigned>(boost::regex_replace(layer, format, "\\2"));
    return std::make_pair(LA_Filler_Fitter::layer_index(isTIB, stereo, layerNum), method);
  }

std::pair< uint32_t, LA_Filler_Fitter::Method > sistrip::MeasureLA::calibration_key ( const uint32_t  detid, const LA_Filler_Fitter::Method  method  ) const

private

Definition at line 189 of file MeasureLA.cc.

References LA_Filler_Fitter::layer_index(), SiStripDetId::subDetector(), SiStripDetId::TIB, TrackerTopology::tibLayer(), TrackerTopology::tibStereo(), TrackerTopology::tobStereo(), and tTopo_.

                                                                       {
    const bool isTIB = SiStripDetId(detid).subDetector() == SiStripDetId::TIB;
    const bool stereo = isTIB ? tTopo_.tibStereo(detid) : tTopo_.tobStereo(detid);
    const unsigned layer = isTIB ? tTopo_.tibLayer(detid) : tTopo_.tobStereo(detid);

    return std::make_pair(LA_Filler_Fitter::layer_index(isTIB, stereo, layer), method);
  }

static std::string sistrip::MeasureLA::granularity ( int32_t  g )

inlinestaticprivate

Definition at line 24 of file MeasureLA.h.

References calibrate(), calibration_key(), LAYER, MODULE, MODULESUMMARY, process_reports(), store_calibrations(), store_methods_and_granularity(), AlCaHLTBitMon_QueryRunRegistry::string, summarize_module_muH_byLayer(), write_report_plots(), write_report_text(), and write_report_text_ms().

Referenced by summarize_module_muH_byLayer(), write_report_plots(), and write_report_text_ms().

                                            {
      switch (g) {
        case LAYER:
          return "_layer";
        case MODULE:
          return "_module";
        case MODULESUMMARY:
          return "_moduleSummary";
      }
      return "";
    };

void sistrip::MeasureLA::process_reports ( ) const

private

Definition at line 88 of file MeasureLA.cc.

References Book::begin(), book, Book::end(), LAYER, LA_Filler_Fitter::layer_results(), AlcaSiPixelAliHarvester0T_cff::method, MODULE, LA_Filler_Fitter::module_results(), MODULESUMMARY, Skims_PA_cff::name, AlCaHLTBitMon_ParallelJobs::p, reports, AlCaHLTBitMon_QueryRunRegistry::string, write_report_plots(), write_report_text(), and write_report_text_ms().

Referenced by granularity(), and MeasureLA().

                                        {
    for (auto const& p : reports) {
      const GRANULARITY gran = (GRANULARITY)p.getParameter<int32_t>("Granularity");
      const std::string name = p.getParameter<std::string>("ReportName");
      const LA_Filler_Fitter::Method method = (LA_Filler_Fitter::Method)p.getParameter<int32_t>("Method");

      write_report_plots(name, method, gran);
      switch (gran) {
        case LAYER:
          write_report_text(name, method, LA_Filler_Fitter::layer_results(book, method));
          break;
        case MODULE:
          write_report_text(name, method, LA_Filler_Fitter::module_results(book, method));
          break;
        case MODULESUMMARY:
          write_report_text_ms(name, method);
          break;
      }
    }

    {
      TFile widthsFile("widths.root", "RECREATE");
      for (Book::const_iterator it = book.begin(".*_width"); it != book.end(); it++)
        if (it->second)
          it->second->Write();
      widthsFile.Close();
    }
  }

std::unique_ptr< SiStripLorentzAngle > sistrip::MeasureLA::produce

(

const SiStripLorentzAngleRcd &

)

Definition at line 52 of file MeasureLA.cc.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), HTMLExport.HTMLExportStatic::export(), and MeasureLA().

                                                                                     {
    auto lorentzAngle = std::make_unique<SiStripLorentzAngle>();
    /*
  std::map<uint32_t,LA_Filler_Fitter::Result> 
    module_results = LA_Filler_Fitter::module_results(book, LA_Filler_Fitter::SQRTVAR);
  
  BOOST_FOREACH(const uint32_t& detid, SiStripDetInfoFileReader(fp_.fullPath()).getAllDetIds()) {
    float la = module_results[detid].measure / module_results[detid].field ;
    lorentzAngle->putLorentzAngle( detid, la );
  }
  */
    return lorentzAngle;
  }

void sistrip::MeasureLA::store_calibrations ( )

private

Definition at line 154 of file MeasureLA.cc.

References calibrations, error_scaling, mps_fire::i, crabWrapper::key, AlcaSiPixelAliHarvester0T_cff::method, offset, unpackBuffers-CaloStage1::offsets, AlCaHLTBitMon_ParallelJobs::p, and slope.

Referenced by granularity(), and MeasureLA().

                                     {
    for (auto const& p : calibrations) {
      LA_Filler_Fitter::Method method = (LA_Filler_Fitter::Method)p.getParameter<int32_t>("Method");
      std::vector<double> slopes(p.getParameter<std::vector<double> >("Slopes"));
      assert(slopes.size() == 14);
      std::vector<double> offsets(p.getParameter<std::vector<double> >("Offsets"));
      assert(offsets.size() == 14);
      std::vector<double> pulls(p.getParameter<std::vector<double> >("Pulls"));
      assert(pulls.size() == 14);

      for (unsigned i = 0; i < 14; i++) {
        const std::pair<unsigned, LA_Filler_Fitter::Method> key(i, method);
        offset[key] = offsets[i];
        slope[key] = slopes[i];
        error_scaling[key] = pulls[i];
      }
    }
  }

void sistrip::MeasureLA::store_methods_and_granularity ( const edm::VParameterSet &  vpset )

private

Definition at line 12 of file MeasureLA.cc.

References byLayer, byModule, methods, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by granularity(), and MeasureLA().

                                                                             {
    for (auto const& p : vpset) {
      methods |= p.getParameter<int32_t>("Method");
      byModule = byModule || p.getParameter<int32_t>("Granularity");
      byLayer = byLayer || !p.getParameter<int32_t>("Granularity");
    }
  }

void sistrip::MeasureLA::summarize_module_muH_byLayer ( const LA_Filler_Fitter &  laff )

private

Definition at line 66 of file MeasureLA.cc.

References Book::begin(), book, calibrate(), calibration_key(), Book::end(), Book::fill(), LA_Filler_Fitter::FIRST_METHOD, granularity(), label, LA_Filler_Fitter::LAST_METHOD, LA_Filler_Fitter::layerLabel(), visualization-live-secondInstance_cfg::m, AlcaSiPixelAliHarvester0T_cff::method, LA_Filler_Fitter::method(), LA_Filler_Fitter::module_results(), MODULESUMMARY, funct::pow(), mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by granularity(), and MeasureLA().

                                                                           {
    for (int m = LA_Filler_Fitter::FIRST_METHOD; m <= LA_Filler_Fitter::LAST_METHOD; m <<= 1) {
      const LA_Filler_Fitter::Method method = (LA_Filler_Fitter::Method)m;
      for (auto& result : LA_Filler_Fitter::module_results(book, method)) {
        calibrate(calibration_key(result.first, method), result.second);
        std::string label =
            laff.layerLabel(result.first) + granularity(MODULESUMMARY) + LA_Filler_Fitter::method(method);
        label = boost::regex_replace(label, boost::regex("layer"), "");

        const double mu_H = -result.second.calMeasured.first / result.second.field;
        const double sigma_mu_H = result.second.calMeasured.second / result.second.field;
        const double weight = pow(1. / sigma_mu_H, 2);

        book.fill(mu_H, label, 150, -0.05, 0.1, weight);
      }
      for (Book::iterator it = book.begin(".*" + granularity(MODULESUMMARY) + ".*"); it != book.end(); ++it) {
        if (it->second->GetEntries())
          it->second->Fit("gaus", "LLQ");
      }
    }
  }

void sistrip::MeasureLA::write_report_plots ( const std::string  name, const LA_Filler_Fitter::Method  method, const GRANULARITY  gran  ) const

private

Definition at line 117 of file MeasureLA.cc.

References Book::begin(), book, Book::end(), FrontierConditions_GlobalTag_cff::file, granularity(), compare::hist, crabWrapper::key, LA_Filler_Fitter::method(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by granularity(), and process_reports().

                                                                                                          {
    TFile file((name + ".root").c_str(), "RECREATE");
    const std::string key = ".*" + granularity(gran) + ".*(" + LA_Filler_Fitter::method(method) + "|" +
                            LA_Filler_Fitter::method(method, false) + ".*)";
    for (Book::const_iterator hist = book.begin(key); hist != book.end(); ++hist)
      if (hist->second)
        hist->second->Write();
    file.Close();
  }

template<class T >

void sistrip::MeasureLA::write_report_text ( const std::string  name, const LA_Filler_Fitter::Method &  _method, const std::map< T, LA_Filler_Fitter::Result > &  _results  ) const

private

Definition at line 128 of file MeasureLA.cc.

References calibrate(), calibration_key(), FrontierConditions_GlobalTag_cff::file, AlcaSiPixelAliHarvester0T_cff::method, MillePedeFileConverter_cfg::out, mps_fire::result, and bookConverter::results.

Referenced by granularity(), and process_reports().

                                                                                             {
    LA_Filler_Fitter::Method method = _method;
    std::map<T, LA_Filler_Fitter::Result> results = _results;
    std::fstream file((name + ".dat").c_str(), std::ios::out);
    for (auto& result : results) {
      calibrate(calibration_key(result.first, method), result.second);
      file << result.first << "\t" << result.second << std::endl;
    }
    file.close();
  }

void sistrip::MeasureLA::write_report_text_ms ( const std::string  name, const LA_Filler_Fitter::Method  method  ) const

private

Definition at line 141 of file MeasureLA.cc.

References Book::begin(), book, Book::end(), f, FrontierConditions_GlobalTag_cff::file, granularity(), crabWrapper::key, LA_Filler_Fitter::method(), MODULESUMMARY, MillePedeFileConverter_cfg::out, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by granularity(), and process_reports().

                                                                                          {
    std::fstream file((name + ".dat").c_str(), std::ios::out);
    const std::string key = ".*" + granularity(MODULESUMMARY) + LA_Filler_Fitter::method(method);
    for (Book::const_iterator it = book.begin(key); it != book.end(); ++it) {
      const TF1* const f = it->second->GetFunction("gaus");
      if (f) {
        file << it->first << "\t" << f->GetParameter(1) << "\t" << f->GetParError(1) << "\t" << f->GetParameter(2)
             << "\t" << f->GetParError(2) << std::endl;
      }
    }
    file.close();
  }

Member Data Documentation

Book sistrip::MeasureLA::book

private

Definition at line 63 of file MeasureLA.h.

Referenced by MeasureLA(), process_reports(), summarize_module_muH_byLayer(), write_report_plots(), and write_report_text_ms().

bool sistrip::MeasureLA::byLayer

private

Definition at line 60 of file MeasureLA.h.

Referenced by MeasureLA(), and store_methods_and_granularity().

bool sistrip::MeasureLA::byModule

private

Definition at line 60 of file MeasureLA.h.

Referenced by MeasureLA(), and store_methods_and_granularity().

const edm::VParameterSet sistrip::MeasureLA::calibrations

private

Definition at line 57 of file MeasureLA.h.

Referenced by store_calibrations().

std::map<std::pair<uint32_t, LA_Filler_Fitter::Method>, float> sistrip::MeasureLA::error_scaling

private

Definition at line 58 of file MeasureLA.h.

Referenced by calibrate(), and store_calibrations().

const edm::FileInPath sistrip::MeasureLA::fp_

private

Definition at line 56 of file MeasureLA.h.

const std::string sistrip::MeasureLA::inFileLocation

private

Definition at line 55 of file MeasureLA.h.

Referenced by MeasureLA().

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

private

Definition at line 54 of file MeasureLA.h.

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

const float sistrip::MeasureLA::localybin

private

Definition at line 61 of file MeasureLA.h.

Referenced by MeasureLA().

const unsigned sistrip::MeasureLA::maxEvents

private

Definition at line 62 of file MeasureLA.h.

Referenced by MeasureLA().

const edm::VParameterSet sistrip::MeasureLA::measurementPreferences

private

Definition at line 57 of file MeasureLA.h.

Referenced by MeasureLA().

int32_t sistrip::MeasureLA::methods

private

Definition at line 59 of file MeasureLA.h.

Referenced by MeasureLA(), and store_methods_and_granularity().

std::map<std::pair<uint32_t, LA_Filler_Fitter::Method>, float> sistrip::MeasureLA::offset

private

Definition at line 58 of file MeasureLA.h.

Referenced by calibrate(), store_calibrations(), and upgradeWorkflowComponents.UpgradeWorkflow::workflow_().

const edm::VParameterSet sistrip::MeasureLA::reports

private

Definition at line 57 of file MeasureLA.h.

Referenced by MeasureLA(), and process_reports().

std::map<std::pair<uint32_t, LA_Filler_Fitter::Method>, float> sistrip::MeasureLA::slope

private

Definition at line 58 of file MeasureLA.h.

Referenced by calibrate(), Vispa.Views.LineDecayView.DecayLine::containsPoint(), store_calibrations(), and Vispa.Views.LineDecayView.DecayLine::transform().

const unsigned sistrip::MeasureLA::stripsperbin

private

Definition at line 62 of file MeasureLA.h.

Referenced by MeasureLA().

TrackerTopology sistrip::MeasureLA::tTopo_

private

Definition at line 65 of file MeasureLA.h.

Referenced by calibration_key(), and MeasureLA().