#include <MillePedeDQMModule.h>

Inheritance diagram for MillePedeDQMModule:

Public Member Functions
virtual void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override

	MillePedeDQMModule (const edm::ParameterSet &)

virtual	~MillePedeDQMModule ()

Public Member Functions inherited from DQMEDHarvester
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

virtual void	beginRun (edm::Run const &, edm::EventSetup const &) override

	DQMEDHarvester (void)

virtual void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	endJob () final

virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

void	endLuminosityBlockProduce (edm::LuminosityBlock &, edm::EventSetup const &) override final

virtual void	endRun (edm::Run const &, edm::EventSetup const &) override

virtual void	produce (edm::Event &, edm::EventSetup const &) override final

Public Member Functions inherited from edm::one::EDProducer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources, edm::EndLuminosityBlockProducer >
	EDProducer ()=default

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

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

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
virtual void	beginRun (const edm::Run &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &)

void	fillExpertHisto (MonitorElement *histo, const std::array< double, 6 > &cut, const std::array< double, 6 > &sigCut, const std::array< double, 6 > &maxMoveCut, const std::array< double, 6 > &maxErrorCut, const std::array< double, 6 > &obs, const std::array< double, 6 > &obsErr)

void	fillExpertHistos ()

int	getIndexFromString (const std::string &alignableId)

bool	setupChanged (const edm::EventSetup &)

Private Attributes
MonitorElement *	h_xPos

MonitorElement *	h_xRot

MonitorElement *	h_yPos

MonitorElement *	h_yRot

MonitorElement *	h_zPos

MonitorElement *	h_zRot

std::unique_ptr< MillePedeFileReader >	mpReader_

const edm::ParameterSet	mpReaderConfig_

std::unique_ptr< AlignableTracker >	tracker_

edm::ESWatcher< IdealGeometryRecord >	watchIdealGeometryRcd_

edm::ESWatcher< PTrackerParametersRcd >	watchPTrackerParametersRcd_

edm::ESWatcher< TrackerTopologyRcd >	watchTrackerTopologyRcd_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
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)

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 44 of file MillePedeDQMModule.h.

Constructor & Destructor Documentation

MillePedeDQMModule::MillePedeDQMModule ( const edm::ParameterSet & config )

Definition at line 36 of file MillePedeDQMModule.cc.

References ~MillePedeDQMModule().

                                                   :
   mpReaderConfig_(
     config.getParameter<edm::ParameterSet>("MillePedeFileReader")
   )
 {
 }

MillePedeDQMModule::~MillePedeDQMModule ( )

virtual

Definition at line 44 of file MillePedeDQMModule.cc.

References bookHistograms().

Referenced by MillePedeDQMModule().

45 {

46 }

Member Function Documentation

void MillePedeDQMModule::beginRun	(	const edm::Run &	,
		const edm::EventSetup &	setup
	)

overrideprivatevirtual

Definition at line 92 of file MillePedeDQMModule.cc.

References edm::ParameterSet::addUntrackedParameter(), TrackerGeomBuilderFromGeometricDet::build(), fillExpertHistos(), reco::get(), edm::EventSetup::get(), MillePedeAlignmentAlgorithm_cfi::pedeLabeler, edm::ESHandle< T >::product(), AlignPCLThresholds::setAlignPCLThresholds(), AlCaHLTBitMon_QueryRunRegistry::string, and sistrip::tracker_.

Referenced by dqmEndJob().

                                                     {
 
   if (!setupChanged(setup)) return;
 
   edm::ESHandle<TrackerTopology> tTopo;
   setup.get<TrackerTopologyRcd>().get(tTopo);
   edm::ESHandle<GeometricDet> geometricDet;
   setup.get<IdealGeometryRecord>().get(geometricDet);
   edm::ESHandle<PTrackerParameters> ptp;
   setup.get<PTrackerParametersRcd>().get(ptp);
 
   // take the thresholds from DB
   edm::ESHandle<AlignPCLThresholds> thresholdHandle;
   setup.get<AlignPCLThresholdsRcd>().get(thresholdHandle);
   auto thresholds_ = thresholdHandle.product();
 
   auto myThresholds = std::make_shared<AlignPCLThresholds>(); 
   myThresholds->setAlignPCLThresholds(thresholds_->getNrecords(),thresholds_->getThreshold_Map());
 
   TrackerGeomBuilderFromGeometricDet builder;
 
   const auto trackerGeometry = builder.build(&(*geometricDet), *ptp, &(*tTopo));
   tracker_ = std::make_unique<AlignableTracker>(trackerGeometry, &(*tTopo));
 
   const std::string labelerPlugin{"PedeLabeler"};
   edm::ParameterSet labelerConfig{};
   labelerConfig.addUntrackedParameter("plugin", labelerPlugin);
   labelerConfig.addUntrackedParameter("RunRangeSelection", edm::VParameterSet{});
 
   std::shared_ptr<PedeLabelerBase> pedeLabeler{
     PedeLabelerPluginFactory::get()
       ->create(labelerPlugin,
               PedeLabelerBase::TopLevelAlignables(tracker_.get(), nullptr, nullptr),
               labelerConfig)
   };
 
 
   mpReader_ = std::make_unique<MillePedeFileReader>(mpReaderConfig_, pedeLabeler, std::shared_ptr<const AlignPCLThresholds>(myThresholds)); 
 
 }

void MillePedeDQMModule::bookHistograms ( DQMStore::IBooker & booker )

private

Definition at line 53 of file MillePedeDQMModule.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::cd(), dqmEndJob(), and DQMStore::IBooker::setCurrentFolder().

Referenced by ~MillePedeDQMModule().

 {
   edm::LogInfo("MillePedeDQMModule") << "Booking histograms";
 
   booker.cd();
   booker.setCurrentFolder("AlCaReco/SiPixelAli/");
 
   h_xPos = booker.book1D("Xpos",   "Alignment fit #DeltaX;;#mum", 36, 0., 36.);
   h_xRot = booker.book1D("Xrot",   "Alignment fit #Delta#theta_{X};;#murad", 36, 0., 36.);
   h_yPos = booker.book1D("Ypos",   "Alignment fit #DeltaY;;#mum", 36, 0., 36.);
   h_yRot = booker.book1D("Yrot",   "Alignment fit #Delta#theta_{Y};;#murad", 36, 0., 36.);
   h_zPos = booker.book1D("Zpos",   "Alignment fit #DeltaZ;;#mum", 36, 0., 36.);
   h_zRot = booker.book1D("Zrot",   "Alignment fit #Delta#theta_{Z};;#murad", 36, 0., 36.);
 
   booker.cd();
 }

void MillePedeDQMModule::dqmEndJob	(	DQMStore::IBooker &	booker,
		DQMStore::IGetter &
	)

overridevirtual

Implements DQMEDHarvester.

Definition at line 72 of file MillePedeDQMModule.cc.

References beginRun(), bookHistograms(), and Exception.

Referenced by bookHistograms().

 {
   bookHistograms(booker);
   if (mpReader_) {
     mpReader_->read();
   } else {
     throw cms::Exception("LogicError")
       << "@SUB=MillePedeDQMModule::dqmEndJob\n"
       << "Try to read MillePede results before initializing MillePedeFileReader";
   }
   fillExpertHistos();
 }

void MillePedeDQMModule::fillExpertHisto	(	MonitorElement *	histo,
		const std::array< double, 6 > &	cut,
		const std::array< double, 6 > &	sigCut,
		const std::array< double, 6 > &	maxMoveCut,
		const std::array< double, 6 > &	maxErrorCut,
		const std::array< double, 6 > &	obs,
		const std::array< double, 6 > &	obsErr
	)

private

Definition at line 202 of file MillePedeDQMModule.cc.

References MonitorElement::getTH1F(), mps_fire::i, and setupChanged().

Referenced by fillExpertHistos().

 {
   TH1F* histo_0 = histo->getTH1F();
 
   double max_ = *std::max_element(maxMoveCut.begin(),maxMoveCut.end());
 
   histo_0->SetMinimum(-(max_));
   histo_0->SetMaximum(  max_);
 
   //  Schematics of the bin contents
   // 
   //  XX XX XX XX XX XX    OO OO OO OO    II II II II 
   // |--|--|--|--|--|--|--|--|--|--|--|--|--|--|--|--|--|
   // | 1| 2| 3| 4| 5| 6| 7| 8| 9|10|11|12|13|14|15|16|17|  ...
   //
   // |-----------------|  |-----------|  |-----------|  
   // |observed movement|  |thresholds1|  |thresholds2|     
 
   for (size_t i = 0; i < obs.size(); ++i) {
     
     // fist obs.size() bins for observed movements
     histo_0->SetBinContent(i+1, obs[i]);
     histo_0->SetBinError(i+1, obsErr[i]); 
     
     // then at bin 8,8+5,8+10,... for cutoffs
     // 5 bins is the space allocated for the 4 other thresholds + 1 empty separation bin
     histo_0->SetBinContent(8+i*5 , cut[i]);
     
     // then at bin 9,9+5,9+10,... for significances
     histo_0->SetBinContent(9+i*5 , sigCut[i]);
 
     // then at bin 10,10+5,10+10,... for maximum movements
     histo_0->SetBinContent(10+i*5, maxMoveCut[i]);
 
     // then at bin 11,11+5,11+10,... for maximum errors
     histo_0->SetBinContent(11+i*5, maxErrorCut[i]);
     
   }
 }

void MillePedeDQMModule::fillExpertHistos ( )

private

Definition at line 135 of file MillePedeDQMModule.cc.

References fillExpertHisto().

Referenced by beginRun().

 {
  
   std::array<double, 6> Xcut_,  sigXcut_,  maxMoveXcut_,  maxErrorXcut_; 
   std::array<double, 6> tXcut_, sigtXcut_, maxMovetXcut_, maxErrortXcut_;
                                                     
   std::array<double, 6> Ycut_,  sigYcut_,  maxMoveYcut_,  maxErrorYcut_; 
   std::array<double, 6> tYcut_, sigtYcut_, maxMovetYcut_, maxErrortYcut_;
                                                     
   std::array<double, 6> Zcut_,  sigZcut_,  maxMoveZcut_,  maxErrorZcut_; 
   std::array<double, 6> tZcut_, sigtZcut_, maxMovetZcut_, maxErrortZcut_;
 
   auto myMap = mpReader_->getThresholdMap();
 
   std::vector<std::string> alignablesList;
   for(auto it = myMap.begin(); it != myMap.end() ; ++it){
     alignablesList.push_back(it->first);
   }
 
   for (auto &alignable : alignablesList){
 
     int detIndex = getIndexFromString(alignable);
 
     Xcut_[detIndex]          = myMap[alignable].getXcut() ;
     sigXcut_[detIndex]       = myMap[alignable].getSigXcut() ;
     maxMoveXcut_[detIndex]   = myMap[alignable].getMaxMoveXcut() ;
     maxErrorXcut_[detIndex]  = myMap[alignable].getErrorXcut() ;
 
     Ycut_[detIndex]          = myMap[alignable].getYcut() ;
     sigYcut_[detIndex]       = myMap[alignable].getSigYcut() ;
     maxMoveYcut_[detIndex]   = myMap[alignable].getMaxMoveYcut() ; 
     maxErrorYcut_[detIndex]  = myMap[alignable].getErrorYcut() ;
  
     Zcut_[detIndex]          = myMap[alignable].getZcut() ;
     sigZcut_[detIndex]       = myMap[alignable].getSigZcut() ;
     maxMoveZcut_[detIndex]   = myMap[alignable].getMaxMoveZcut() ;
     maxErrorZcut_[detIndex]  = myMap[alignable].getErrorZcut() ;
 
     tXcut_[detIndex]         = myMap[alignable].getThetaXcut() ;
     sigtXcut_[detIndex]      = myMap[alignable].getSigThetaXcut() ; 
     maxMovetXcut_[detIndex]  = myMap[alignable].getMaxMoveThetaXcut() ;
     maxErrortXcut_[detIndex] = myMap[alignable].getErrorThetaXcut() ;
 
     tYcut_[detIndex]         = myMap[alignable].getThetaYcut() ;
     sigtYcut_[detIndex]      = myMap[alignable].getSigThetaYcut() ;
     maxMovetYcut_[detIndex]  = myMap[alignable].getMaxMoveThetaYcut() ;
     maxErrortYcut_[detIndex] = myMap[alignable].getErrorThetaYcut() ;
  
     tZcut_[detIndex]         = myMap[alignable].getThetaZcut() ;
     sigtZcut_[detIndex]      = myMap[alignable].getSigThetaYcut() ;
     maxMovetZcut_[detIndex]  = myMap[alignable].getMaxMoveThetaYcut() ;
     maxErrortZcut_[detIndex] = myMap[alignable].getErrorThetaYcut() ;
 
   }
 
   fillExpertHisto(h_xPos,  Xcut_, sigXcut_,  maxMoveXcut_,  maxErrorXcut_,  mpReader_->getXobs(),  mpReader_->getXobsErr());
   fillExpertHisto(h_xRot, tXcut_, sigtXcut_, maxMovetXcut_, maxErrortXcut_, mpReader_->getTXobs(), mpReader_->getTXobsErr());
 
   fillExpertHisto(h_yPos,  Ycut_, sigYcut_,  maxMoveYcut_,  maxErrorYcut_,  mpReader_->getYobs(),  mpReader_->getYobsErr());
   fillExpertHisto(h_yRot, tYcut_, sigtYcut_, maxMovetYcut_, maxErrortYcut_, mpReader_->getTYobs(), mpReader_->getTYobsErr());
 
   fillExpertHisto(h_zPos,  Zcut_, sigZcut_,  maxMoveZcut_,  maxErrorZcut_,  mpReader_->getZobs(),  mpReader_->getZobsErr());
   fillExpertHisto(h_zRot, tZcut_, sigtZcut_, maxMovetZcut_, maxErrortZcut_, mpReader_->getTZobs(), mpReader_->getTZobsErr());
 
 }

int MillePedeDQMModule::getIndexFromString ( const std::string & alignableId )

private

Definition at line 259 of file MillePedeDQMModule.cc.

References Exception.

Referenced by setupChanged().

 {
   
   if(alignableId == "TPBHalfBarrelXminus"){
     return 3;
   } else if(alignableId == "TPBHalfBarrelXplus"){
     return 2;
   } else if(alignableId == "TPEHalfCylinderXminusZminus") {
     return 1;
   } else if(alignableId == "TPEHalfCylinderXplusZminus") {
     return 0;
   } else if(alignableId == "TPEHalfCylinderXminusZplus") {
     return 5;
   } else if(alignableId == "TPEHalfCylinderXplusZplus") {
     return 4;
   } else{
     throw cms::Exception("LogicError") 
       << "@SUB=MillePedeDQMModule::getIndexFromString\n"
       << "Retrieving conversion for not supported Alignable partition" 
       << alignableId;
   }
 }

bool MillePedeDQMModule::setupChanged ( const edm::EventSetup & setup )

private

Definition at line 246 of file MillePedeDQMModule.cc.

References getIndexFromString().

Referenced by fillExpertHisto().

 {
   bool changed{false};
 
   if (watchIdealGeometryRcd_.check(setup)) changed = true;
   if (watchTrackerTopologyRcd_.check(setup)) changed = true;
   if (watchPTrackerParametersRcd_.check(setup)) changed = true;
 
   return changed;
 }

Member Data Documentation

MonitorElement* MillePedeDQMModule::h_xPos

private

Definition at line 86 of file MillePedeDQMModule.h.

MonitorElement* MillePedeDQMModule::h_xRot

private

Definition at line 87 of file MillePedeDQMModule.h.

MonitorElement* MillePedeDQMModule::h_yPos

private

Definition at line 88 of file MillePedeDQMModule.h.

MonitorElement* MillePedeDQMModule::h_yRot

private

Definition at line 89 of file MillePedeDQMModule.h.

MonitorElement* MillePedeDQMModule::h_zPos

private

Definition at line 90 of file MillePedeDQMModule.h.

MonitorElement* MillePedeDQMModule::h_zRot

private

Definition at line 91 of file MillePedeDQMModule.h.

std::unique_ptr<MillePedeFileReader> MillePedeDQMModule::mpReader_

private

Definition at line 79 of file MillePedeDQMModule.h.

const edm::ParameterSet MillePedeDQMModule::mpReaderConfig_

private

Definition at line 77 of file MillePedeDQMModule.h.

std::unique_ptr<AlignableTracker> MillePedeDQMModule::tracker_

private

Definition at line 78 of file MillePedeDQMModule.h.

edm::ESWatcher<IdealGeometryRecord> MillePedeDQMModule::watchIdealGeometryRcd_

private

Definition at line 82 of file MillePedeDQMModule.h.

edm::ESWatcher<PTrackerParametersRcd> MillePedeDQMModule::watchPTrackerParametersRcd_

private

Definition at line 83 of file MillePedeDQMModule.h.

edm::ESWatcher<TrackerTopologyRcd> MillePedeDQMModule::watchTrackerTopologyRcd_

private

Definition at line 81 of file MillePedeDQMModule.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation