CMS 3D CMS Logo
List of all members | Public Types | Public Member Functions | Private Member Functions | Private Attributes
MillePedeDQMModule Class Reference

#include <MillePedeDQMModule.h>

Inheritance diagram for MillePedeDQMModule:
DQMEDHarvester edm::one::EDProducer< edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::EndProcessBlockProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources, edm::Accumulator > edm::one::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Types

enum  { SIZE_LG_STRUCTS = 6, SIZE_HG_STRUCTS = 820, SIZE_INDEX = 8 }
 
- Public Types inherited from DQMEDHarvester
typedef dqm::harvesting::DQMStore DQMStore
 
typedef dqm::harvesting::MonitorElement MonitorElement
 
- Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase ModuleType
 
- Public Types inherited from edm::ProducerBase
template<typename T >
using BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >
 
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > >
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 

Public Member Functions

void dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override
 
 MillePedeDQMModule (const edm::ParameterSet &)
 
 ~MillePedeDQMModule () override
 
- Public Member Functions inherited from DQMEDHarvester
void accumulate (edm::Event const &ev, edm::EventSetup const &es) final
 
void beginJob () override
 
void beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &) final
 
void beginRun (edm::Run const &run, edm::EventSetup const &) override
 
virtual void dqmAnalyze (DQMStore::IBooker &, DQMStore::IGetter &, edm::Event const &, edm::EventSetup const &)
 
 DQMEDHarvester (edm::ParameterSet const &iConfig)
 
 DQMEDHarvester ()
 
virtual void dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)
 
virtual void dqmEndRun (DQMStore::IBooker &, DQMStore::IGetter &, edm::Run const &, edm::EventSetup const &)
 
void endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
 
void endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &es) final
 
void endProcessBlockProduce (edm::ProcessBlock &) final
 
void endRun (edm::Run const &, edm::EventSetup const &) override
 
void endRunProduce (edm::Run &run, edm::EventSetup const &es) final
 
 ~DQMEDHarvester () override=default
 
- Public Member Functions inherited from edm::one::EDProducer< edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::EndProcessBlockProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources, edm::Accumulator >
 EDProducer ()=default
 
 EDProducer (const EDProducer &)=delete
 
SerialTaskQueue * globalLuminosityBlocksQueue () final
 
SerialTaskQueue * globalRunsQueue () final
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
const EDProducer & operator= (const EDProducer &)=delete
 
bool wantsGlobalLuminosityBlocks () const noexcept final
 
bool wantsGlobalRuns () const noexcept final
 
bool wantsInputProcessBlocks () const noexcept final
 
bool wantsProcessBlocks () const noexcept final
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const noexcept
 
bool wantsStreamRuns () const noexcept
 
 ~EDProducerBase () override
 
- 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 ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
std::vector< bool > const & recordProvenanceList () const
 
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, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
TypeLabelList const & typeLabelList () const
 used by the fwk to register the list of products of this module More...
 
 ~ProducerBase () noexcept(false) 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
 
ESResolverIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESResolverIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

void beginRun (const edm::Run &, const edm::EventSetup &) override
 
void bookHistograms (DQMStore::IBooker &)
 
void fillExpertHisto (MonitorElement *histo, const std::array< double, SIZE_INDEX > &cut, const std::array< double, SIZE_INDEX > &sigCut, const std::array< double, SIZE_INDEX > &maxMoveCut, const std::array< double, SIZE_INDEX > &maxErrorCut, const std::array< double, SIZE_LG_STRUCTS > &obs, const std::array< double, SIZE_LG_STRUCTS > &obsErr)
 
void fillExpertHisto_HG (std::map< std::string, MonitorElement *> &histo_map, const std::array< double, SIZE_INDEX > &cut, const std::array< double, SIZE_INDEX > &sigCut, const std::array< double, SIZE_INDEX > &maxMoveCut, const std::array< double, SIZE_INDEX > &maxErrorCut, const std::array< double, SIZE_HG_STRUCTS > &obs, const std::array< double, SIZE_HG_STRUCTS > &obsErr)
 
void fillExpertHistos ()
 
void fillExpertHistos_HG ()
 
void fillStatusHisto (MonitorElement *statusHisto)
 
void fillStatusHistoHG (MonitorElement *statusHisto)
 
int getIndexFromString (const std::string &alignableId)
 
bool setupChanged (const edm::EventSetup &)
 

Private Attributes

const edm::ESGetToken< AlignPCLThresholdsHG, AlignPCLThresholdsHGRcd > aliThrToken_
 
MonitorElement * binariesAvalaible
 
MonitorElement * exitCode
 
const edm::ESGetToken< GeometricDet, IdealGeometryRecord > gDetToken_
 
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
 
MonitorElement * h_xPos
 
std::map< std::string, MonitorElement * > h_xPos_HG
 
MonitorElement * h_xRot
 
std::map< std::string, MonitorElement * > h_xRot_HG
 
MonitorElement * h_yPos
 
std::map< std::string, MonitorElement * > h_yPos_HG
 
MonitorElement * h_yRot
 
std::map< std::string, MonitorElement * > h_yRot_HG
 
MonitorElement * h_zPos
 
std::map< std::string, MonitorElement * > h_zPos_HG
 
MonitorElement * h_zRot
 
std::map< std::string, MonitorElement * > h_zRot_HG
 
bool isHG_
 
MonitorElement * isVetoed
 
std::vector< std::pair< std::string, int > > layerVec
 
std::unique_ptr< MillePedeFileReader > mpReader_
 
const edm::ParameterSet mpReaderConfig_
 
const std::string outputFolder_
 
std::shared_ptr< PixelTopologyMap > pixelTopologyMap_
 
const edm::ESGetToken< PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd > ptitpToken_
 
const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd > ptpToken_
 
MonitorElement * statusResults
 
std::unique_ptr< AlignableTracker > tracker_
 
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > tTopoToken_
 
edm::ESWatcher< IdealGeometryRecord > watchIdealGeometryRcd_
 
edm::ESWatcher< PTrackerParametersRcd > watchPTrackerParametersRcd_
 
edm::ESWatcher< TrackerTopologyRcd > watchTrackerTopologyRcd_
 

Additional Inherited Members

- 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::ProducerBase
template<Transition Tr = Transition::Event>
auto produces (std::string instanceName) noexcept
 declare what type of product will make and with which optional label More...
 
template<Transition B>
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
template<BranchType B>
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<class ProductType >
BranchAliasSetterT< ProductType > produces ()
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces ()
 
template<class ProductType >
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces ()
 
template<Transition Tr = Transition::Event>
auto produces () noexcept
 
ProducesCollector producesCollector ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B > consumes (edm::InputTag tag) noexcept
 
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
constexpr auto esConsumes ()
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag)
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
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)
 
void resetItemsToGetFrom (BranchType iType)
 
- Protected Attributes inherited from DQMEDHarvester
DQMStore * dqmstore_
 
edm::GetterOfProducts< DQMToken > jobmegetter_
 
edm::EDPutTokenT< DQMToken > jobToken_
 
edm::GetterOfProducts< DQMToken > lumimegetter_
 
edm::EDPutTokenT< DQMToken > lumiToken_
 
edm::GetterOfProducts< DQMToken > runmegetter_
 
edm::EDPutTokenT< DQMToken > runToken_
 

Detailed Description

Definition at line 51 of file MillePedeDQMModule.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum
Enumerator
SIZE_LG_STRUCTS 
SIZE_HG_STRUCTS 
SIZE_INDEX 

Definition at line 59 of file MillePedeDQMModule.h.

59 { SIZE_LG_STRUCTS = 6, SIZE_HG_STRUCTS = 820, SIZE_INDEX = 8 };

Constructor & Destructor Documentation

◆ MillePedeDQMModule()

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

Definition at line 28 of file MillePedeDQMModule.cc.

29  : tTopoToken_(esConsumes<edm::Transition::BeginRun>()),
30  gDetToken_(esConsumes<edm::Transition::BeginRun>()),
31  ptpToken_(esConsumes<edm::Transition::BeginRun>()),
32  ptitpToken_(esConsumes<edm::Transition::BeginRun>()),
33  aliThrToken_(esConsumes<edm::Transition::BeginRun>()),
34  geomToken_(esConsumes<edm::Transition::BeginRun>()),
35  outputFolder_(config.getParameter<std::string>("outputFolder")),
36  mpReaderConfig_(config.getParameter<edm::ParameterSet>("MillePedeFileReader")),
37  isHG_(mpReaderConfig_.getParameter<bool>("isHG")) {
38  consumes<AlignmentToken, edm::InProcess>(config.getParameter<edm::InputTag>("alignmentTokenSrc"));
39 }
MillePedeDQMModule::mpReaderConfig_
const edm::ParameterSet mpReaderConfig_
Definition: MillePedeDQMModule.h:106
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
MillePedeDQMModule::gDetToken_
const edm::ESGetToken< GeometricDet, IdealGeometryRecord > gDetToken_
Definition: MillePedeDQMModule.h:99
config
Definition: config.py:1
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
MillePedeDQMModule::ptitpToken_
const edm::ESGetToken< PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd > ptitpToken_
Definition: MillePedeDQMModule.h:101
MillePedeDQMModule::tTopoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > tTopoToken_
Definition: MillePedeDQMModule.h:98
MillePedeDQMModule::aliThrToken_
const edm::ESGetToken< AlignPCLThresholdsHG, AlignPCLThresholdsHGRcd > aliThrToken_
Definition: MillePedeDQMModule.h:102
MillePedeDQMModule::ptpToken_
const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd > ptpToken_
Definition: MillePedeDQMModule.h:100
MillePedeDQMModule::outputFolder_
const std::string outputFolder_
Definition: MillePedeDQMModule.h:105
MillePedeDQMModule::geomToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
Definition: MillePedeDQMModule.h:103

◆ ~MillePedeDQMModule()

MillePedeDQMModule::~MillePedeDQMModule ( )
override

Definition at line 41 of file MillePedeDQMModule.cc.

41 {}

Member Function Documentation

◆ beginRun()

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

Definition at line 174 of file MillePedeDQMModule.cc.

References edm::ParameterSet::addUntrackedParameter(), aliThrToken_, TrackerGeomBuilderFromGeometricDet::build(), gDetToken_, relativeConstraints::geom, geomToken_, get, mpReader_, mpReaderConfig_, MillePedeAlignmentAlgorithm_cfi::pedeLabeler, pixelTopologyMap_, ptitpToken_, ptpToken_, singleTopDQM_cfi::setup, setupChanged(), AlCaHLTBitMon_QueryRunRegistry::string, tracker_, and tTopoToken_.

174  {
175  if (!setupChanged(setup))
176  return;
177 
178  const TrackerTopology* const tTopo = &setup.getData(tTopoToken_);
179  const GeometricDet* geometricDet = &setup.getData(gDetToken_);
180  const PTrackerParameters* ptp = &setup.getData(ptpToken_);
181  const PTrackerAdditionalParametersPerDet* ptitp = &setup.getData(ptitpToken_);
182  const TrackerGeometry* geom = &setup.getData(geomToken_);
183 
184  pixelTopologyMap_ = std::make_shared<PixelTopologyMap>(geom, tTopo);
185 
186  // take the thresholds from DB
187  const auto& thresholds_ = &setup.getData(aliThrToken_);
188 
189  auto myThresholds = std::make_shared<AlignPCLThresholdsHG>();
190  myThresholds->setAlignPCLThresholds(thresholds_->getNrecords(), thresholds_->getThreshold_Map());
191  myThresholds->setFloatMap(thresholds_->getFloatMap());
192 
193  TrackerGeomBuilderFromGeometricDet builder;
194 
195  const auto trackerGeometry = builder.build(geometricDet, ptitp, *ptp, tTopo);
196  tracker_ = std::make_unique<AlignableTracker>(trackerGeometry, tTopo);
197 
198  const std::string labelerPlugin{"PedeLabeler"};
199  edm::ParameterSet labelerConfig{};
200  labelerConfig.addUntrackedParameter("plugin", labelerPlugin);
201  labelerConfig.addUntrackedParameter("RunRangeSelection", edm::VParameterSet{});
202 
203  std::shared_ptr<PedeLabelerBase> pedeLabeler{PedeLabelerPluginFactory::get()->create(
204  labelerPlugin, PedeLabelerBase::TopLevelAlignables(tracker_.get(), nullptr, nullptr), labelerConfig)};
205 
206  mpReader_ = std::make_unique<MillePedeFileReader>(
207  mpReaderConfig_, pedeLabeler, std::shared_ptr<const AlignPCLThresholdsHG>(myThresholds), pixelTopologyMap_);
208 }
MillePedeDQMModule::mpReaderConfig_
const edm::ParameterSet mpReaderConfig_
Definition: MillePedeDQMModule.h:106
edm::VParameterSet
std::vector< ParameterSet > VParameterSet
Definition: ParameterSet.h:35
TrackerGeomBuilderFromGeometricDet::build
TrackerGeometry * build(const GeometricDet *gd, const PTrackerAdditionalParametersPerDet *ptitp, const PTrackerParameters &ptp, const TrackerTopology *tTopo)
Definition: TrackerGeomBuilderFromGeometricDet.cc:46
MillePedeDQMModule::setupChanged
bool setupChanged(const edm::EventSetup &)
Definition: MillePedeDQMModule.cc:499
MillePedeDQMModule::gDetToken_
const edm::ESGetToken< GeometricDet, IdealGeometryRecord > gDetToken_
Definition: MillePedeDQMModule.h:99
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
MillePedeDQMModule::ptitpToken_
const edm::ESGetToken< PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd > ptitpToken_
Definition: MillePedeDQMModule.h:101
singleTopDQM_cfi.setup
setup
SETUP
Definition: singleTopDQM_cfi.py:37
MillePedeDQMModule::tTopoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > tTopoToken_
Definition: MillePedeDQMModule.h:98
MillePedeDQMModule::pixelTopologyMap_
std::shared_ptr< PixelTopologyMap > pixelTopologyMap_
Definition: MillePedeDQMModule.h:109
edm::ParameterSet::addUntrackedParameter
void addUntrackedParameter(std::string const &name, T const &value)
Definition: ParameterSet.h:193
MillePedeDQMModule::aliThrToken_
const edm::ESGetToken< AlignPCLThresholdsHG, AlignPCLThresholdsHGRcd > aliThrToken_
Definition: MillePedeDQMModule.h:102
MillePedeAlignmentAlgorithm_cfi.pedeLabeler
pedeLabeler
Definition: MillePedeAlignmentAlgorithm_cfi.py:38
MillePedeDQMModule::ptpToken_
const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd > ptpToken_
Definition: MillePedeDQMModule.h:100
MillePedeDQMModule::tracker_
std::unique_ptr< AlignableTracker > tracker_
Definition: MillePedeDQMModule.h:107
MillePedeDQMModule::mpReader_
std::unique_ptr< MillePedeFileReader > mpReader_
Definition: MillePedeDQMModule.h:108
get
#define get
MillePedeDQMModule::geomToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
Definition: MillePedeDQMModule.h:103

◆ bookHistograms()

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

Definition at line 47 of file MillePedeDQMModule.cc.

References binariesAvalaible, dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookInt(), dqm::implementation::IBooker::bookString(), dqm::implementation::NavigatorBase::cd(), Exception, exitCode, h_xPos, h_xPos_HG, h_xRot, h_xRot_HG, h_yPos, h_yPos_HG, h_yRot, h_yRot_HG, h_zPos, h_zPos_HG, h_zRot, h_zRot_HG, isHG_, isVetoed, nano_mu_digi_cff::layer, layerVec, outputFolder_, pixelTopologyMap_, dqm::implementation::NavigatorBase::setCurrentFolder(), and statusResults.

Referenced by dqmEndJob().

47  {
48  edm::LogInfo("MillePedeDQMModule") << "Booking histograms";
49 
50  booker.cd();
51  if (!isHG_) {
52  if (outputFolder_.find("HG") != std::string::npos) {
53  throw cms::Exception("LogicError")
54  << "MillePedeDQMModule is configured as Low Granularity but the outputfolder is for High Granularity";
55  }
56 
57  booker.setCurrentFolder(outputFolder_);
58  h_xPos = booker.book1D("Xpos", "Alignment fit #DeltaX;;#mum", 36, 0., 36.);
59  h_xRot = booker.book1D("Xrot", "Alignment fit #Delta#theta_{X};;#murad", 36, 0., 36.);
60  h_yPos = booker.book1D("Ypos", "Alignment fit #DeltaY;;#mum", 36, 0., 36.);
61  h_yRot = booker.book1D("Yrot", "Alignment fit #Delta#theta_{Y};;#murad", 36, 0., 36.);
62  h_zPos = booker.book1D("Zpos", "Alignment fit #DeltaZ;;#mum", 36, 0., 36.);
63  h_zRot = booker.book1D("Zrot", "Alignment fit #Delta#theta_{Z};;#murad", 36, 0., 36.);
64  statusResults = booker.book2D("statusResults", "Status of SiPixelAli PCL workflow;;", 6, 0., 6., 1, 0., 1.);
65  } else {
66  if (outputFolder_.find("HG") == std::string::npos) {
67  throw cms::Exception("LogicError")
68  << "MillePedeDQMModule is configured as High Granularity but the outputfolder is for Low Granularity";
69  }
70 
71  booker.setCurrentFolder(outputFolder_);
72 
73  layerVec = {{"Layer1", pixelTopologyMap_->getPXBLadders(1)},
74  {"Layer2", pixelTopologyMap_->getPXBLadders(2)},
75  {"Layer3", pixelTopologyMap_->getPXBLadders(3)},
76  {"Layer4", pixelTopologyMap_->getPXBLadders(4)},
77  {"Disk-3", pixelTopologyMap_->getPXFBlades(-3) * 2},
78  {"Disk-2", pixelTopologyMap_->getPXFBlades(-2) * 2},
79  {"Disk-1", pixelTopologyMap_->getPXFBlades(-1) * 2},
80  {"Disk1", pixelTopologyMap_->getPXFBlades(1) * 2},
81  {"Disk2", pixelTopologyMap_->getPXFBlades(2) * 2},
82  {"Disk3", pixelTopologyMap_->getPXFBlades(3) * 2}};
83 
84  for (const auto& layer : layerVec) {
85  h_xPos_HG[layer.first] = booker.book1D("Xpos_HG_" + layer.first,
86  "Alignment fit #DeltaX for " + layer.first + ";;#mum",
87  layer.second + 5,
88  0.,
89  layer.second + 5);
90  h_xRot_HG[layer.first] = booker.book1D("Xrot_HG_" + layer.first,
91  "Alignment fit #Delta#theta_{X} for " + layer.first + ";;#murad",
92  layer.second + 5,
93  0.,
94  layer.second + 5);
95  h_yPos_HG[layer.first] = booker.book1D("Ypos_HG_" + layer.first,
96  "Alignment fit #DeltaY for " + layer.first + ";;#mum",
97  layer.second + 5,
98  0.,
99  layer.second + 5);
100  h_yRot_HG[layer.first] = booker.book1D("Yrot_HG_" + layer.first,
101  "Alignment fit #Delta#theta_{Y} for " + layer.first + ";;#murad",
102  layer.second + 5,
103  0.,
104  layer.second + 5);
105  h_zPos_HG[layer.first] = booker.book1D("Zpos_HG_" + layer.first,
106  "Alignment fit #DeltaZ for " + layer.first + ";;#mum",
107  layer.second + 5,
108  0.,
109  layer.second + 5);
110  h_zRot_HG[layer.first] = booker.book1D("Zrot_HG_" + layer.first,
111  "Alignment fit #Delta#theta_{Z} for " + layer.first + ";;#murad",
112  layer.second + 5,
113  0.,
114  layer.second + 5);
115  }
116 
117  statusResults =
118  booker.book2D("statusResults", "Fraction threshold check for SiPixelAliHG PCL;;", 6, 0., 6., 10, 0., 10.);
119  }
120 
121  binariesAvalaible = booker.bookInt("BinariesFound");
122  exitCode = booker.bookString("PedeExitCode", "");
123  isVetoed = booker.bookString("IsVetoed", "");
124 
125  booker.cd();
126 }
MillePedeDQMModule::h_xRot_HG
std::map< std::string, MonitorElement * > h_xRot_HG
Definition: MillePedeDQMModule.h:126
MillePedeDQMModule::layerVec
std::vector< std::pair< std::string, int > > layerVec
Definition: MillePedeDQMModule.h:111
MillePedeDQMModule::h_yPos_HG
std::map< std::string, MonitorElement * > h_yPos_HG
Definition: MillePedeDQMModule.h:127
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:36
MillePedeDQMModule::h_zPos
MonitorElement * h_zPos
Definition: MillePedeDQMModule.h:122
MillePedeDQMModule::h_zRot_HG
std::map< std::string, MonitorElement * > h_zRot_HG
Definition: MillePedeDQMModule.h:130
MillePedeDQMModule::binariesAvalaible
MonitorElement * binariesAvalaible
Definition: MillePedeDQMModule.h:133
dqm::implementation::IBooker::bookString
MonitorElement * bookString(TString const &name, TString const &value, FUNC onbooking=NOOP())
Definition: DQMStore.h:87
MillePedeDQMModule::isVetoed
MonitorElement * isVetoed
Definition: MillePedeDQMModule.h:135
MillePedeDQMModule::h_yRot_HG
std::map< std::string, MonitorElement * > h_yRot_HG
Definition: MillePedeDQMModule.h:128
MillePedeDQMModule::h_xPos_HG
std::map< std::string, MonitorElement * > h_xPos_HG
Definition: MillePedeDQMModule.h:125
MillePedeDQMModule::h_xPos
MonitorElement * h_xPos
Definition: MillePedeDQMModule.h:118
MillePedeDQMModule::h_yPos
MonitorElement * h_yPos
Definition: MillePedeDQMModule.h:120
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131
nano_mu_digi_cff.layer
layer
Definition: nano_mu_digi_cff.py:28
MillePedeDQMModule::pixelTopologyMap_
std::shared_ptr< PixelTopologyMap > pixelTopologyMap_
Definition: MillePedeDQMModule.h:109
MillePedeDQMModule::statusResults
MonitorElement * statusResults
Definition: MillePedeDQMModule.h:132
dqm::implementation::IBooker::book2D
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:221
dqm::implementation::IBooker::bookInt
MonitorElement * bookInt(TString const &name, FUNC onbooking=NOOP())
Definition: DQMStore.h:73
MillePedeDQMModule::h_xRot
MonitorElement * h_xRot
Definition: MillePedeDQMModule.h:119
dqm::implementation::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
MillePedeDQMModule::h_zRot
MonitorElement * h_zRot
Definition: MillePedeDQMModule.h:123
MillePedeDQMModule::exitCode
MonitorElement * exitCode
Definition: MillePedeDQMModule.h:134
MillePedeDQMModule::h_zPos_HG
std::map< std::string, MonitorElement * > h_zPos_HG
Definition: MillePedeDQMModule.h:129
MillePedeDQMModule::outputFolder_
const std::string outputFolder_
Definition: MillePedeDQMModule.h:105
MillePedeDQMModule::h_yRot
MonitorElement * h_yRot
Definition: MillePedeDQMModule.h:121

◆ dqmEndJob()

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

Implements DQMEDHarvester.

Definition at line 128 of file MillePedeDQMModule.cc.

References binariesAvalaible, bookHistograms(), Exception, exitCode, dqm::impl::MonitorElement::Fill(), fillExpertHistos(), fillExpertHistos_HG(), fillStatusHisto(), fillStatusHistoHG(), isHG_, isVetoed, mpReader_, statusResults, and AlCaHLTBitMon_QueryRunRegistry::string.

128  {
129  bookHistograms(booker);
130  if (mpReader_) {
131  mpReader_->read();
132  } else {
133  throw cms::Exception("LogicError") << "@SUB=MillePedeDQMModule::dqmEndJob\n"
134  << "Try to read MillePede results before initializing MillePedeFileReader";
135  }
136  if (!isHG_) {
137  fillExpertHistos();
138  fillStatusHisto(statusResults);
139  } else {
140  fillExpertHistos_HG();
141  fillStatusHistoHG(statusResults);
142  }
143  binariesAvalaible->Fill(mpReader_->binariesAmount());
144  auto theResults = mpReader_->getResults();
145  std::string exitCodeStr = theResults.getExitMessage();
146 
147  std::string vetoStr{};
148  if (mpReader_->storeAlignments()) {
149  vetoStr = "DB Updated!"; /* easy peasy, fait accompli an alignment is there */
150  } else {
151  if (theResults.isHighGranularity()) { /* HG case */
152  if (theResults.getDBVetoed() && theResults.getDBUpdated()) {
153  vetoStr = "DB Update Vetoed"; /* this can happen in the HG PCL case */
154  } else {
155  vetoStr = "N/A";
156  }
157  } else { /* LG case */
158  if (theResults.exceedsCutoffs()) {
159  vetoStr = "DB Update Vetoed"; /* this can happen in the LG PCL case */
160  } else {
161  vetoStr = "N/A";
162  } // if the alignment exceeds the cutoffs
163  } // LG case
164  } // if the alignment was not stored
165 
166  exitCode->Fill(exitCodeStr);
167  isVetoed->Fill(vetoStr);
168 }
MillePedeDQMModule::fillStatusHisto
void fillStatusHisto(MonitorElement *statusHisto)
Definition: MillePedeDQMModule.cc:210
MillePedeDQMModule::binariesAvalaible
MonitorElement * binariesAvalaible
Definition: MillePedeDQMModule.h:133
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:292
MillePedeDQMModule::isVetoed
MonitorElement * isVetoed
Definition: MillePedeDQMModule.h:135
MillePedeDQMModule::bookHistograms
void bookHistograms(DQMStore::IBooker &)
Definition: MillePedeDQMModule.cc:47
MillePedeDQMModule::statusResults
MonitorElement * statusResults
Definition: MillePedeDQMModule.h:132
MillePedeDQMModule::mpReader_
std::unique_ptr< MillePedeFileReader > mpReader_
Definition: MillePedeDQMModule.h:108
MillePedeDQMModule::exitCode
MonitorElement * exitCode
Definition: MillePedeDQMModule.h:134
MillePedeDQMModule::fillStatusHistoHG
void fillStatusHistoHG(MonitorElement *statusHisto)
Definition: MillePedeDQMModule.cc:228

◆ fillExpertHisto()

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

Definition at line 312 of file MillePedeDQMModule.cc.

References TkAlMuonSelectors_cfi::cut, timingPdfMaker::histo, mps_fire::i, and HLT_2024v13_cff::sigCut.

Referenced by fillExpertHistos().

318  {
319  TH1F* histo_0 = histo->getTH1F();
320 
321  double max_ = *std::max_element(maxMoveCut.begin(), maxMoveCut.end());
322 
323  histo_0->SetMinimum(-(max_));
324  histo_0->SetMaximum(max_);
325 
326  // Schematics of the bin contents
327  //
328  // XX XX XX XX XX XX OO OO OO OO II II II II
329  // |--|--|--|--|--|--|--|--|--|--|--|--|--|--|--|--|--|
330  // | 1| 2| 3| 4| 5| 6| 7| 8| 9|10|11|12|13|14|15|16|17| ...
331  //
332  // |-----------------| |-----------| |-----------|
333  // |observed movement| |thresholds1| |thresholds2|
334 
335  for (size_t i = 0; i < obs.size(); ++i) {
336  // fist obs.size() bins for observed movements
337  histo_0->SetBinContent(i + 1, obs[i]);
338  histo_0->SetBinError(i + 1, obsErr[i]);
339 
340  // then at bin 8,8+5,8+10,... for cutoffs
341  // 5 bins is the space allocated for the 4 other thresholds + 1 empty separation bin
342  histo_0->SetBinContent(8 + i * 5, cut[i]);
343 
344  // then at bin 9,9+5,9+10,... for significances
345  histo_0->SetBinContent(9 + i * 5, sigCut[i]);
346 
347  // then at bin 10,10+5,10+10,... for maximum movements
348  histo_0->SetBinContent(10 + i * 5, maxMoveCut[i]);
349 
350  // then at bin 11,11+5,11+10,... for maximum errors
351  histo_0->SetBinContent(11 + i * 5, maxErrorCut[i]);
352  }
353 }
mps_fire.i
i
Definition: mps_fire.py:429
TkAlMuonSelectors_cfi.cut
cut
Definition: TkAlMuonSelectors_cfi.py:25
HLT_2024v13_cff.sigCut
sigCut
Definition: HLT_2024v13_cff.py:38898

◆ fillExpertHisto_HG()

void MillePedeDQMModule::fillExpertHisto_HG ( std::map< std::string, MonitorElement *> &  histo_map,
const std::array< double, SIZE_INDEX > &  cut,
const std::array< double, SIZE_INDEX > &  sigCut,
const std::array< double, SIZE_INDEX > &  maxMoveCut,
const std::array< double, SIZE_INDEX > &  maxErrorCut,
const std::array< double, SIZE_HG_STRUCTS > &  obs,
const std::array< double, SIZE_HG_STRUCTS > &  obsErr 
)
private

Definition at line 437 of file MillePedeDQMModule.cc.

References funct::abs(), newFWLiteAna::bin, mps_splice::currentStart, TkAlMuonSelectors_cfi::cut, mps_fire::i, nano_mu_digi_cff::layer, layerVec, SiStripPI::max, and HLT_2024v13_cff::sigCut.

Referenced by fillExpertHistos_HG().

443  {
444  int currentStart = 0;
445  int bin = 0;
446  double max_ = 0;
447 
448  for (const auto& layer : layerVec) {
449  TH1F* histo_0 = histo_map[layer.first]->getTH1F();
450 
451  max_ = -1;
452  for (int i = currentStart; i < (currentStart + layer.second); ++i) {
453  // first obs.size() bins for observed movements
454  bin = i - currentStart + 1;
455 
456  // fill observed values
457  histo_0->SetBinContent(bin, obs[i]);
458  histo_0->SetBinError(bin, obsErr[i]);
459 
460  if (std::abs(obs[i]) > max_) {
461  max_ = std::abs(obs[i]);
462  }
463  }
464 
465  // five extra bins at the end, one empty, one with threshold, one with sigCut, one with maxMoveCut, one with MaxErrorCut
466  histo_0->SetBinContent(bin + 1, 0);
467  histo_0->SetBinError(bin + 1, 0);
468 
469  int detIndex;
470  if (layer.first.find("Disk") != std::string::npos) {
471  // 7 is the detId for panels, see getIndexFromString
472  detIndex = 7;
473  histo_0->GetXaxis()->SetTitle("Panel");
474  } else {
475  // 6 is the detId for ladders, see getIndexFromString
476  detIndex = 6;
477  histo_0->GetXaxis()->SetTitle("Ladder");
478  }
479 
480  histo_0->SetBinContent(bin + 2, cut[detIndex]);
481  histo_0->SetBinError(bin + 2, 0);
482  histo_0->SetBinContent(bin + 3, sigCut[detIndex]);
483  histo_0->SetBinError(bin + 3, 0);
484  histo_0->SetBinContent(bin + 4, maxMoveCut[detIndex]);
485  histo_0->SetBinError(bin + 4, 0);
486  histo_0->SetBinContent(bin + 5, maxErrorCut[detIndex]);
487  histo_0->SetBinError(bin + 5, 0);
488 
489  // always scale so the cutoff is visible
490  max_ = std::max(cut[detIndex] * 1.2, max_);
491 
492  histo_0->SetMinimum(-(max_) * 1.2);
493  histo_0->SetMaximum(max_ * 1.2);
494 
495  currentStart += layer.second;
496  }
497 }
mps_fire.i
i
Definition: mps_fire.py:429
MillePedeDQMModule::layerVec
std::vector< std::pair< std::string, int > > layerVec
Definition: MillePedeDQMModule.h:111
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
nano_mu_digi_cff.layer
layer
Definition: nano_mu_digi_cff.py:28
mps_splice.currentStart
int currentStart
Definition: mps_splice.py:66
TkAlMuonSelectors_cfi.cut
cut
Definition: TkAlMuonSelectors_cfi.py:25
HLT_2024v13_cff.sigCut
sigCut
Definition: HLT_2024v13_cff.py:38898

◆ fillExpertHistos()

void MillePedeDQMModule::fillExpertHistos ( )
private

Definition at line 248 of file MillePedeDQMModule.cc.

References fillExpertHisto(), getIndexFromString(), h_xPos, h_xRot, h_yPos, h_yRot, h_zPos, h_zRot, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, and mpReader_.

Referenced by dqmEndJob().

248  {
249  std::array<double, SIZE_INDEX> Xcut_, sigXcut_, maxMoveXcut_, maxErrorXcut_;
250  std::array<double, SIZE_INDEX> tXcut_, sigtXcut_, maxMovetXcut_, maxErrortXcut_;
251 
252  std::array<double, SIZE_INDEX> Ycut_, sigYcut_, maxMoveYcut_, maxErrorYcut_;
253  std::array<double, SIZE_INDEX> tYcut_, sigtYcut_, maxMovetYcut_, maxErrortYcut_;
254 
255  std::array<double, SIZE_INDEX> Zcut_, sigZcut_, maxMoveZcut_, maxErrorZcut_;
256  std::array<double, SIZE_INDEX> tZcut_, sigtZcut_, maxMovetZcut_, maxErrortZcut_;
257 
258  auto myMap = mpReader_->getThresholdMap();
259 
260  std::vector<std::string> alignablesList;
261  for (auto it = myMap.begin(); it != myMap.end(); ++it) {
262  alignablesList.push_back(it->first);
263  }
264 
265  for (auto& alignable : alignablesList) {
266  int detIndex = getIndexFromString(alignable);
267 
268  Xcut_[detIndex] = myMap[alignable].getXcut();
269  sigXcut_[detIndex] = myMap[alignable].getSigXcut();
270  maxMoveXcut_[detIndex] = myMap[alignable].getMaxMoveXcut();
271  maxErrorXcut_[detIndex] = myMap[alignable].getErrorXcut();
272 
273  Ycut_[detIndex] = myMap[alignable].getYcut();
274  sigYcut_[detIndex] = myMap[alignable].getSigYcut();
275  maxMoveYcut_[detIndex] = myMap[alignable].getMaxMoveYcut();
276  maxErrorYcut_[detIndex] = myMap[alignable].getErrorYcut();
277 
278  Zcut_[detIndex] = myMap[alignable].getZcut();
279  sigZcut_[detIndex] = myMap[alignable].getSigZcut();
280  maxMoveZcut_[detIndex] = myMap[alignable].getMaxMoveZcut();
281  maxErrorZcut_[detIndex] = myMap[alignable].getErrorZcut();
282 
283  tXcut_[detIndex] = myMap[alignable].getThetaXcut();
284  sigtXcut_[detIndex] = myMap[alignable].getSigThetaXcut();
285  maxMovetXcut_[detIndex] = myMap[alignable].getMaxMoveThetaXcut();
286  maxErrortXcut_[detIndex] = myMap[alignable].getErrorThetaXcut();
287 
288  tYcut_[detIndex] = myMap[alignable].getThetaYcut();
289  sigtYcut_[detIndex] = myMap[alignable].getSigThetaYcut();
290  maxMovetYcut_[detIndex] = myMap[alignable].getMaxMoveThetaYcut();
291  maxErrortYcut_[detIndex] = myMap[alignable].getErrorThetaYcut();
292 
293  tZcut_[detIndex] = myMap[alignable].getThetaZcut();
294  sigtZcut_[detIndex] = myMap[alignable].getSigThetaZcut();
295  maxMovetZcut_[detIndex] = myMap[alignable].getMaxMoveThetaZcut();
296  maxErrortZcut_[detIndex] = myMap[alignable].getErrorThetaZcut();
297  }
298 
299  fillExpertHisto(h_xPos, Xcut_, sigXcut_, maxMoveXcut_, maxErrorXcut_, mpReader_->getXobs(), mpReader_->getXobsErr());
300  fillExpertHisto(
301  h_xRot, tXcut_, sigtXcut_, maxMovetXcut_, maxErrortXcut_, mpReader_->getTXobs(), mpReader_->getTXobsErr());
302 
303  fillExpertHisto(h_yPos, Ycut_, sigYcut_, maxMoveYcut_, maxErrorYcut_, mpReader_->getYobs(), mpReader_->getYobsErr());
304  fillExpertHisto(
305  h_yRot, tYcut_, sigtYcut_, maxMovetYcut_, maxErrortYcut_, mpReader_->getTYobs(), mpReader_->getTYobsErr());
306 
307  fillExpertHisto(h_zPos, Zcut_, sigZcut_, maxMoveZcut_, maxErrorZcut_, mpReader_->getZobs(), mpReader_->getZobsErr());
308  fillExpertHisto(
309  h_zRot, tZcut_, sigtZcut_, maxMovetZcut_, maxErrortZcut_, mpReader_->getTZobs(), mpReader_->getTZobsErr());
310 }
MillePedeDQMModule::fillExpertHisto
void fillExpertHisto(MonitorElement *histo, const std::array< double, SIZE_INDEX > &cut, const std::array< double, SIZE_INDEX > &sigCut, const std::array< double, SIZE_INDEX > &maxMoveCut, const std::array< double, SIZE_INDEX > &maxErrorCut, const std::array< double, SIZE_LG_STRUCTS > &obs, const std::array< double, SIZE_LG_STRUCTS > &obsErr)
Definition: MillePedeDQMModule.cc:312
MillePedeDQMModule::h_zPos
MonitorElement * h_zPos
Definition: MillePedeDQMModule.h:122
MillePedeDQMModule::h_xPos
MonitorElement * h_xPos
Definition: MillePedeDQMModule.h:118
MillePedeDQMModule::h_yPos
MonitorElement * h_yPos
Definition: MillePedeDQMModule.h:120
MillePedeDQMModule::h_xRot
MonitorElement * h_xRot
Definition: MillePedeDQMModule.h:119
ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it
auto & it
Definition: splitVertices.h:48
MillePedeDQMModule::mpReader_
std::unique_ptr< MillePedeFileReader > mpReader_
Definition: MillePedeDQMModule.h:108
MillePedeDQMModule::getIndexFromString
int getIndexFromString(const std::string &alignableId)
Definition: MillePedeDQMModule.cc:512
MillePedeDQMModule::h_zRot
MonitorElement * h_zRot
Definition: MillePedeDQMModule.h:123
MillePedeDQMModule::h_yRot
MonitorElement * h_yRot
Definition: MillePedeDQMModule.h:121

◆ fillExpertHistos_HG()

void MillePedeDQMModule::fillExpertHistos_HG ( )
private

Definition at line 355 of file MillePedeDQMModule.cc.

References fillExpertHisto_HG(), getIndexFromString(), h_xPos_HG, h_xRot_HG, h_yPos_HG, h_yRot_HG, h_zPos_HG, h_zRot_HG, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, and mpReader_.

Referenced by dqmEndJob().

355  {
356  std::array<double, SIZE_INDEX> Xcut_, sigXcut_, maxMoveXcut_, maxErrorXcut_;
357  std::array<double, SIZE_INDEX> tXcut_, sigtXcut_, maxMovetXcut_, maxErrortXcut_;
358 
359  std::array<double, SIZE_INDEX> Ycut_, sigYcut_, maxMoveYcut_, maxErrorYcut_;
360  std::array<double, SIZE_INDEX> tYcut_, sigtYcut_, maxMovetYcut_, maxErrortYcut_;
361 
362  std::array<double, SIZE_INDEX> Zcut_, sigZcut_, maxMoveZcut_, maxErrorZcut_;
363  std::array<double, SIZE_INDEX> tZcut_, sigtZcut_, maxMovetZcut_, maxErrortZcut_;
364 
365  auto myMap = mpReader_->getThresholdMap();
366 
367  std::vector<std::string> alignablesList;
368  for (auto it = myMap.begin(); it != myMap.end(); ++it) {
369  alignablesList.push_back(it->first);
370  }
371 
372  for (auto& alignable : alignablesList) {
373  int detIndex = getIndexFromString(alignable);
374 
375  Xcut_[detIndex] = myMap[alignable].getXcut();
376  sigXcut_[detIndex] = myMap[alignable].getSigXcut();
377  maxMoveXcut_[detIndex] = myMap[alignable].getMaxMoveXcut();
378  maxErrorXcut_[detIndex] = myMap[alignable].getErrorXcut();
379 
380  Ycut_[detIndex] = myMap[alignable].getYcut();
381  sigYcut_[detIndex] = myMap[alignable].getSigYcut();
382  maxMoveYcut_[detIndex] = myMap[alignable].getMaxMoveYcut();
383  maxErrorYcut_[detIndex] = myMap[alignable].getErrorYcut();
384 
385  Zcut_[detIndex] = myMap[alignable].getZcut();
386  sigZcut_[detIndex] = myMap[alignable].getSigZcut();
387  maxMoveZcut_[detIndex] = myMap[alignable].getMaxMoveZcut();
388  maxErrorZcut_[detIndex] = myMap[alignable].getErrorZcut();
389 
390  tXcut_[detIndex] = myMap[alignable].getThetaXcut();
391  sigtXcut_[detIndex] = myMap[alignable].getSigThetaXcut();
392  maxMovetXcut_[detIndex] = myMap[alignable].getMaxMoveThetaXcut();
393  maxErrortXcut_[detIndex] = myMap[alignable].getErrorThetaXcut();
394 
395  tYcut_[detIndex] = myMap[alignable].getThetaYcut();
396  sigtYcut_[detIndex] = myMap[alignable].getSigThetaYcut();
397  maxMovetYcut_[detIndex] = myMap[alignable].getMaxMoveThetaYcut();
398  maxErrortYcut_[detIndex] = myMap[alignable].getErrorThetaYcut();
399 
400  tZcut_[detIndex] = myMap[alignable].getThetaZcut();
401  sigtZcut_[detIndex] = myMap[alignable].getSigThetaZcut();
402  maxMovetZcut_[detIndex] = myMap[alignable].getMaxMoveThetaZcut();
403  maxErrortZcut_[detIndex] = myMap[alignable].getErrorThetaZcut();
404  }
405 
406  fillExpertHisto_HG(
407  h_xPos_HG, Xcut_, sigXcut_, maxMoveXcut_, maxErrorXcut_, mpReader_->getXobs_HG(), mpReader_->getXobsErr_HG());
408  fillExpertHisto_HG(h_xRot_HG,
409  tXcut_,
410  sigtXcut_,
411  maxMovetXcut_,
412  maxErrortXcut_,
413  mpReader_->getTXobs_HG(),
414  mpReader_->getTXobsErr_HG());
415 
416  fillExpertHisto_HG(
417  h_yPos_HG, Ycut_, sigYcut_, maxMoveYcut_, maxErrorYcut_, mpReader_->getYobs_HG(), mpReader_->getYobsErr_HG());
418  fillExpertHisto_HG(h_yRot_HG,
419  tYcut_,
420  sigtYcut_,
421  maxMovetYcut_,
422  maxErrortYcut_,
423  mpReader_->getTYobs_HG(),
424  mpReader_->getTYobsErr_HG());
425 
426  fillExpertHisto_HG(
427  h_zPos_HG, Zcut_, sigZcut_, maxMoveZcut_, maxErrorZcut_, mpReader_->getZobs_HG(), mpReader_->getZobsErr_HG());
428  fillExpertHisto_HG(h_zRot_HG,
429  tZcut_,
430  sigtZcut_,
431  maxMovetZcut_,
432  maxErrortZcut_,
433  mpReader_->getTZobs_HG(),
434  mpReader_->getTZobsErr_HG());
435 }
MillePedeDQMModule::h_xRot_HG
std::map< std::string, MonitorElement * > h_xRot_HG
Definition: MillePedeDQMModule.h:126
MillePedeDQMModule::h_yPos_HG
std::map< std::string, MonitorElement * > h_yPos_HG
Definition: MillePedeDQMModule.h:127
MillePedeDQMModule::h_zRot_HG
std::map< std::string, MonitorElement * > h_zRot_HG
Definition: MillePedeDQMModule.h:130
MillePedeDQMModule::h_yRot_HG
std::map< std::string, MonitorElement * > h_yRot_HG
Definition: MillePedeDQMModule.h:128
MillePedeDQMModule::h_xPos_HG
std::map< std::string, MonitorElement * > h_xPos_HG
Definition: MillePedeDQMModule.h:125
ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it
auto & it
Definition: splitVertices.h:48
MillePedeDQMModule::mpReader_
std::unique_ptr< MillePedeFileReader > mpReader_
Definition: MillePedeDQMModule.h:108
MillePedeDQMModule::getIndexFromString
int getIndexFromString(const std::string &alignableId)
Definition: MillePedeDQMModule.cc:512
MillePedeDQMModule::h_zPos_HG
std::map< std::string, MonitorElement * > h_zPos_HG
Definition: MillePedeDQMModule.h:129
MillePedeDQMModule::fillExpertHisto_HG
void fillExpertHisto_HG(std::map< std::string, MonitorElement *> &histo_map, const std::array< double, SIZE_INDEX > &cut, const std::array< double, SIZE_INDEX > &sigCut, const std::array< double, SIZE_INDEX > &maxMoveCut, const std::array< double, SIZE_INDEX > &maxErrorCut, const std::array< double, SIZE_HG_STRUCTS > &obs, const std::array< double, SIZE_HG_STRUCTS > &obsErr)
Definition: MillePedeDQMModule.cc:437

◆ fillStatusHisto()

void MillePedeDQMModule::fillStatusHisto ( MonitorElement *  statusHisto)
private

Definition at line 210 of file MillePedeDQMModule.cc.

References dqm::legacy::MonitorElement::getTH2F(), and mpReader_.

Referenced by dqmEndJob().

210  {
211  TH2F* histo_status = statusHisto->getTH2F();
212  auto theResults = mpReader_->getResults();
213  theResults.print();
214  histo_status->SetBinContent(1, 1, theResults.getDBUpdated());
215  histo_status->GetXaxis()->SetBinLabel(1, "DB updated");
216  histo_status->SetBinContent(2, 1, theResults.exceedsCutoffs());
217  histo_status->GetXaxis()->SetBinLabel(2, "significant movement");
218  histo_status->SetBinContent(3, 1, theResults.getDBVetoed());
219  histo_status->GetXaxis()->SetBinLabel(3, "DB update vetoed");
220  histo_status->SetBinContent(4, 1, !theResults.exceedsThresholds());
221  histo_status->GetXaxis()->SetBinLabel(4, "within max movement");
222  histo_status->SetBinContent(5, 1, !theResults.exceedsMaxError());
223  histo_status->GetXaxis()->SetBinLabel(5, "within max error");
224  histo_status->SetBinContent(6, 1, !theResults.belowSignificance());
225  histo_status->GetXaxis()->SetBinLabel(6, "above significance");
226 }
dqm::legacy::MonitorElement::getTH2F
virtual TH2F * getTH2F() const
Definition: MonitorElement.h:499
MillePedeDQMModule::mpReader_
std::unique_ptr< MillePedeFileReader > mpReader_
Definition: MillePedeDQMModule.h:108

◆ fillStatusHistoHG()

void MillePedeDQMModule::fillStatusHistoHG ( MonitorElement *  statusHisto)
private

Definition at line 228 of file MillePedeDQMModule.cc.

References dqm::legacy::MonitorElement::getTH2F(), mps_fire::i, dqmiolumiharvest::j, mpReader_, and mps_fire::result.

Referenced by dqmEndJob().

228  {
229  TH2F* histo_status = statusHisto->getTH2F();
230  auto& theResults = mpReader_->getResultsHG();
231  histo_status->GetXaxis()->SetBinLabel(1, "#DeltaX");
232  histo_status->GetXaxis()->SetBinLabel(2, "#DeltaY");
233  histo_status->GetXaxis()->SetBinLabel(3, "#DeltaZ");
234  histo_status->GetXaxis()->SetBinLabel(4, "#Delta#theta_{X}");
235  histo_status->GetXaxis()->SetBinLabel(5, "#Delta#theta_{Y}");
236  histo_status->GetXaxis()->SetBinLabel(6, "#Delta#theta_{Z}");
237 
238  int i = 0;
239  for (const auto& result : theResults) {
240  histo_status->GetYaxis()->SetBinLabel(i + 1, result.first.data());
241  for (std::size_t j = 0; j < result.second.size(); ++j) {
242  histo_status->SetBinContent(j + 1, i + 1, result.second[j]);
243  }
244  i++;
245  }
246 }
mps_fire.i
i
Definition: mps_fire.py:429
mps_fire.result
result
Definition: mps_fire.py:311
dqm::legacy::MonitorElement::getTH2F
virtual TH2F * getTH2F() const
Definition: MonitorElement.h:499
MillePedeDQMModule::mpReader_
std::unique_ptr< MillePedeFileReader > mpReader_
Definition: MillePedeDQMModule.h:108

◆ getIndexFromString()

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

Definition at line 512 of file MillePedeDQMModule.cc.

References Exception.

Referenced by fillExpertHistos(), and fillExpertHistos_HG().

512  {
513  if (alignableId == "TPBHalfBarrelXminus") {
514  return 3;
515  } else if (alignableId == "TPBHalfBarrelXplus") {
516  return 2;
517  } else if (alignableId == "TPEHalfCylinderXminusZminus") {
518  return 1;
519  } else if (alignableId == "TPEHalfCylinderXplusZminus") {
520  return 0;
521  } else if (alignableId == "TPEHalfCylinderXminusZplus") {
522  return 5;
523  } else if (alignableId == "TPEHalfCylinderXplusZplus") {
524  return 4;
525  } else if (alignableId.rfind("TPBLadder", 0) == 0) {
526  return 6;
527  } else if (alignableId.rfind("TPEPanel", 0) == 0) {
528  return 7;
529  } else {
530  throw cms::Exception("LogicError") << "@SUB=MillePedeDQMModule::getIndexFromString\n"
531  << "Retrieving conversion for not supported Alignable partition" << alignableId;
532  }
533 }

◆ setupChanged()

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

Definition at line 499 of file MillePedeDQMModule.cc.

References edm::ESWatcher< T >::check(), singleTopDQM_cfi::setup, watchIdealGeometryRcd_, watchPTrackerParametersRcd_, and watchTrackerTopologyRcd_.

Referenced by beginRun().

499  {
500  bool changed{false};
501 
502  if (watchIdealGeometryRcd_.check(setup))
503  changed = true;
504  if (watchTrackerTopologyRcd_.check(setup))
505  changed = true;
506  if (watchPTrackerParametersRcd_.check(setup))
507  changed = true;
508 
509  return changed;
510 }
MillePedeDQMModule::watchIdealGeometryRcd_
edm::ESWatcher< IdealGeometryRecord > watchIdealGeometryRcd_
Definition: MillePedeDQMModule.h:114
singleTopDQM_cfi.setup
setup
SETUP
Definition: singleTopDQM_cfi.py:37
edm::ESWatcher::check
bool check(const edm::EventSetup &iSetup)
Definition: ESWatcher.h:57
MillePedeDQMModule::watchTrackerTopologyRcd_
edm::ESWatcher< TrackerTopologyRcd > watchTrackerTopologyRcd_
Definition: MillePedeDQMModule.h:113
MillePedeDQMModule::watchPTrackerParametersRcd_
edm::ESWatcher< PTrackerParametersRcd > watchPTrackerParametersRcd_
Definition: MillePedeDQMModule.h:115

Member Data Documentation

◆ aliThrToken_

const edm::ESGetToken<AlignPCLThresholdsHG, AlignPCLThresholdsHGRcd> MillePedeDQMModule::aliThrToken_
private

Definition at line 102 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ binariesAvalaible

MonitorElement* MillePedeDQMModule::binariesAvalaible
private

Definition at line 133 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and dqmEndJob().

◆ exitCode

MonitorElement* MillePedeDQMModule::exitCode
private

Definition at line 134 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and dqmEndJob().

◆ gDetToken_

const edm::ESGetToken<GeometricDet, IdealGeometryRecord> MillePedeDQMModule::gDetToken_
private

Definition at line 99 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ geomToken_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> MillePedeDQMModule::geomToken_
private

Definition at line 103 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ h_xPos

MonitorElement* MillePedeDQMModule::h_xPos
private

Definition at line 118 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos().

◆ h_xPos_HG

std::map<std::string, MonitorElement*> MillePedeDQMModule::h_xPos_HG
private

Definition at line 125 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos_HG().

◆ h_xRot

MonitorElement* MillePedeDQMModule::h_xRot
private

Definition at line 119 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos().

◆ h_xRot_HG

std::map<std::string, MonitorElement*> MillePedeDQMModule::h_xRot_HG
private

Definition at line 126 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos_HG().

◆ h_yPos

MonitorElement* MillePedeDQMModule::h_yPos
private

Definition at line 120 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos().

◆ h_yPos_HG

std::map<std::string, MonitorElement*> MillePedeDQMModule::h_yPos_HG
private

Definition at line 127 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos_HG().

◆ h_yRot

MonitorElement* MillePedeDQMModule::h_yRot
private

Definition at line 121 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos().

◆ h_yRot_HG

std::map<std::string, MonitorElement*> MillePedeDQMModule::h_yRot_HG
private

Definition at line 128 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos_HG().

◆ h_zPos

MonitorElement* MillePedeDQMModule::h_zPos
private

Definition at line 122 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos().

◆ h_zPos_HG

std::map<std::string, MonitorElement*> MillePedeDQMModule::h_zPos_HG
private

Definition at line 129 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos_HG().

◆ h_zRot

MonitorElement* MillePedeDQMModule::h_zRot
private

Definition at line 123 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos().

◆ h_zRot_HG

std::map<std::string, MonitorElement*> MillePedeDQMModule::h_zRot_HG
private

Definition at line 130 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos_HG().

◆ isHG_

bool MillePedeDQMModule::isHG_
private

Definition at line 137 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and dqmEndJob().

◆ isVetoed

MonitorElement* MillePedeDQMModule::isVetoed
private

Definition at line 135 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and dqmEndJob().

◆ layerVec

std::vector<std::pair<std::string, int> > MillePedeDQMModule::layerVec
private

Definition at line 111 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHisto_HG().

◆ mpReader_

std::unique_ptr<MillePedeFileReader> MillePedeDQMModule::mpReader_
private

Definition at line 108 of file MillePedeDQMModule.h.

Referenced by beginRun(), dqmEndJob(), fillExpertHistos(), fillExpertHistos_HG(), fillStatusHisto(), and fillStatusHistoHG().

◆ mpReaderConfig_

const edm::ParameterSet MillePedeDQMModule::mpReaderConfig_
private

Definition at line 106 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ outputFolder_

const std::string MillePedeDQMModule::outputFolder_
private

Definition at line 105 of file MillePedeDQMModule.h.

Referenced by bookHistograms().

◆ pixelTopologyMap_

std::shared_ptr<PixelTopologyMap> MillePedeDQMModule::pixelTopologyMap_
private

Definition at line 109 of file MillePedeDQMModule.h.

Referenced by beginRun(), and bookHistograms().

◆ ptitpToken_

const edm::ESGetToken<PTrackerAdditionalParametersPerDet, PTrackerAdditionalParametersPerDetRcd> MillePedeDQMModule::ptitpToken_
private

Definition at line 101 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ ptpToken_

const edm::ESGetToken<PTrackerParameters, PTrackerParametersRcd> MillePedeDQMModule::ptpToken_
private

Definition at line 100 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ statusResults

MonitorElement* MillePedeDQMModule::statusResults
private

Definition at line 132 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and dqmEndJob().

◆ tracker_

std::unique_ptr<AlignableTracker> MillePedeDQMModule::tracker_
private

Definition at line 107 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ tTopoToken_

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> MillePedeDQMModule::tTopoToken_
private

Definition at line 98 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ watchIdealGeometryRcd_

edm::ESWatcher<IdealGeometryRecord> MillePedeDQMModule::watchIdealGeometryRcd_
private

Definition at line 114 of file MillePedeDQMModule.h.

Referenced by setupChanged().

◆ watchPTrackerParametersRcd_

edm::ESWatcher<PTrackerParametersRcd> MillePedeDQMModule::watchPTrackerParametersRcd_
private

Definition at line 115 of file MillePedeDQMModule.h.

Referenced by setupChanged().

◆ watchTrackerTopologyRcd_

edm::ESWatcher<TrackerTopologyRcd> MillePedeDQMModule::watchTrackerTopologyRcd_
private

Definition at line 113 of file MillePedeDQMModule.h.

Referenced by setupChanged().

Generated for CMSSW Reference Manual by doxygen 1.8.14