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 ()
ESProxyIndex const *	getTokenIndices (unsigned int iIndex) const
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &) final
	~ESProducer ()(false) override
Public Member Functions inherited from edm::ESProxyFactoryProducer
	ESProxyFactoryProducer ()
void	newInterval (const eventsetup::EventSetupRecordKey &iRecordType, const ValidityInterval &iInterval) override
	overrides DataProxyProvider method More...
	~ESProxyFactoryProducer () noexcept(false) override
Public Member Functions inherited from edm::eventsetup::DataProxyProvider
	DataProxyProvider ()
const ComponentDescription &	description () const
bool	isUsingRecord (const EventSetupRecordKey &) const
const KeyedProxies &	keyedProxies (const EventSetupRecordKey &iRecordKey) const
void	resetProxies (const EventSetupRecordKey &iRecordType)
void	resetProxiesIfTransient (const EventSetupRecordKey &iRecordType)
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
Public Types inherited from edm::eventsetup::DataProxyProvider
typedef std::vector< std::pair< DataKey, edm::propagate_const< std::shared_ptr< DataProxy > > > >	KeyedProxies
typedef std::vector< EventSetupRecordKey >	Keys
typedef std::map< EventSetupRecordKey, KeyedProxies >	RecordProxies
Static Public Member Functions inherited from edm::eventsetup::DataProxyProvider
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::ESProducer
	ESProducer (const ESProducer &)=delete
ESProducer const &	operator= (const ESProducer &)=delete
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 eventsetup::EventSetupRecordKey &iRecord, std::unique_ptr< eventsetup::ProxyFactoryBase > iFactory, const std::string &iLabel=std::string())
void	registerProxies (const eventsetup::EventSetupRecordKey &iRecord, KeyedProxies &aProxyList) override
	override DataProxyProvider method More...
Protected Member Functions inherited from edm::eventsetup::DataProxyProvider
void	eraseAll (const EventSetupRecordKey &iRecordKey)
	deletes all the Proxies in aStream More...
void	invalidateProxies (const EventSetupRecordKey &iRecordKey)
template<class T >
void	usingRecord ()
void	usingRecordWithKey (const EventSetupRecordKey &)

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 26 of file MeasureLA.h.

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

Constructor & Destructor Documentation

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

explicit

Definition at line 22 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 170 of file MeasureLA.cc.

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

Referenced by granularity(), store_calibrations(), 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 176 of file MeasureLA.cc.

References cmsPerfStripChart::format, and LA_Filler_Fitter::layer_index().

Referenced by calibrate(), 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 185 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 27 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 93 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, dataset::name, AlCaHLTBitMon_ParallelJobs::p, reports, AlCaHLTBitMon_QueryRunRegistry::string, write_report_plots(), write_report_text(), and write_report_text_ms().

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

                        {
  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 53 of file MeasureLA.cc.

References summarize_module_muH_byLayer().

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 152 of file MeasureLA.cc.

References calibrate(), calibrations, error_scaling, mps_fire::i, crabWrapper::key, AlcaSiPixelAliHarvester0T_cff::method, offset, AlCaHLTBitMon_ParallelJobs::p, and slope.

Referenced by granularity(), MeasureLA(), and write_report_text_ms().

                     {
  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 14 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 71 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(), funct::m, AlcaSiPixelAliHarvester0T_cff::method, LA_Filler_Fitter::method(), LA_Filler_Fitter::module_results(), MODULESUMMARY, funct::pow(), process_reports(), mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by granularity(), MeasureLA(), and produce().

                                                           {
  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 113 of file MeasureLA.cc.

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

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 123 of file MeasureLA.cc.

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

Referenced by granularity(), process_reports(), and write_report_plots().

                                                                                                                                 {
  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 135 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, store_calibrations(), and AlCaHLTBitMon_QueryRunRegistry::string.

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

                                                                          {
  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 56 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 53 of file MeasureLA.h.

Referenced by MeasureLA(), and store_methods_and_granularity().

bool sistrip::MeasureLA::byModule

private

Definition at line 53 of file MeasureLA.h.

Referenced by MeasureLA(), and store_methods_and_granularity().

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

private

Definition at line 50 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 51 of file MeasureLA.h.

Referenced by calibrate(), and store_calibrations().

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

private

Definition at line 49 of file MeasureLA.h.

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

private

Definition at line 48 of file MeasureLA.h.

Referenced by MeasureLA().

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

private

Definition at line 47 of file MeasureLA.h.

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

const float sistrip::MeasureLA::localybin

private

Definition at line 54 of file MeasureLA.h.

Referenced by MeasureLA().

const unsigned sistrip::MeasureLA::maxEvents

private

Definition at line 55 of file MeasureLA.h.

Referenced by MeasureLA().

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

private

Definition at line 50 of file MeasureLA.h.

Referenced by MeasureLA().

int32_t sistrip::MeasureLA::methods

private

Definition at line 52 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 51 of file MeasureLA.h.

Referenced by calibrate(), and store_calibrations().

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

private

Definition at line 50 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 51 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 55 of file MeasureLA.h.

Referenced by MeasureLA().

TrackerTopology sistrip::MeasureLA::tTopo_

private

Definition at line 58 of file MeasureLA.h.

Referenced by calibration_key(), and MeasureLA().