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 final
 
bool wantsGlobalRuns () const final
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () const final
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~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< 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 50 of file MillePedeDQMModule.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum
Enumerator
SIZE_LG_STRUCTS 
SIZE_HG_STRUCTS 
SIZE_INDEX 

Definition at line 58 of file MillePedeDQMModule.h.

58 { 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  aliThrToken_(esConsumes<edm::Transition::BeginRun>()),
33  geomToken_(esConsumes<edm::Transition::BeginRun>()),
34  outputFolder_(config.getParameter<std::string>("outputFolder")),
35  mpReaderConfig_(config.getParameter<edm::ParameterSet>("MillePedeFileReader")),
36  isHG_(mpReaderConfig_.getParameter<bool>("isHG")) {
37  consumes<AlignmentToken, edm::InProcess>(config.getParameter<edm::InputTag>("alignmentTokenSrc"));
38 }
MillePedeDQMModule::mpReaderConfig_
const edm::ParameterSet mpReaderConfig_
Definition: MillePedeDQMModule.h:104
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
MillePedeDQMModule::gDetToken_
const edm::ESGetToken< GeometricDet, IdealGeometryRecord > gDetToken_
Definition: MillePedeDQMModule.h:98
config
Definition: config.py:1
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
MillePedeDQMModule::tTopoToken_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > tTopoToken_
Definition: MillePedeDQMModule.h:97
MillePedeDQMModule::aliThrToken_
const edm::ESGetToken< AlignPCLThresholdsHG, AlignPCLThresholdsHGRcd > aliThrToken_
Definition: MillePedeDQMModule.h:100
MillePedeDQMModule::ptpToken_
const edm::ESGetToken< PTrackerParameters, PTrackerParametersRcd > ptpToken_
Definition: MillePedeDQMModule.h:99
MillePedeDQMModule::outputFolder_
const std::string outputFolder_
Definition: MillePedeDQMModule.h:103
MillePedeDQMModule::geomToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
Definition: MillePedeDQMModule.h:101

◆ ~MillePedeDQMModule()

MillePedeDQMModule::~MillePedeDQMModule ( )
override

Definition at line 40 of file MillePedeDQMModule.cc.

40 {}

Member Function Documentation

◆ beginRun()

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

Definition at line 173 of file MillePedeDQMModule.cc.

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

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

◆ bookHistograms()

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

Definition at line 46 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().

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

◆ dqmEndJob()

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

Implements DQMEDHarvester.

Definition at line 127 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.

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

◆ 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 310 of file MillePedeDQMModule.cc.

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

Referenced by fillExpertHistos().

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

◆ 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 435 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_2024v10_cff::sigCut.

Referenced by fillExpertHistos_HG().

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

◆ fillExpertHistos()

void MillePedeDQMModule::fillExpertHistos ( )
private

Definition at line 246 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().

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

◆ fillExpertHistos_HG()

void MillePedeDQMModule::fillExpertHistos_HG ( )
private

Definition at line 353 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().

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

◆ fillStatusHisto()

void MillePedeDQMModule::fillStatusHisto ( MonitorElement *  statusHisto)
private

Definition at line 208 of file MillePedeDQMModule.cc.

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

Referenced by dqmEndJob().

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

◆ fillStatusHistoHG()

void MillePedeDQMModule::fillStatusHistoHG ( MonitorElement *  statusHisto)
private

Definition at line 226 of file MillePedeDQMModule.cc.

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

Referenced by dqmEndJob().

226  {
227  TH2F* histo_status = statusHisto->getTH2F();
228  auto& theResults = mpReader_->getResultsHG();
229  histo_status->GetXaxis()->SetBinLabel(1, "#DeltaX");
230  histo_status->GetXaxis()->SetBinLabel(2, "#DeltaY");
231  histo_status->GetXaxis()->SetBinLabel(3, "#DeltaZ");
232  histo_status->GetXaxis()->SetBinLabel(4, "#Delta#theta_{X}");
233  histo_status->GetXaxis()->SetBinLabel(5, "#Delta#theta_{Y}");
234  histo_status->GetXaxis()->SetBinLabel(6, "#Delta#theta_{Z}");
235 
236  int i = 0;
237  for (const auto& result : theResults) {
238  histo_status->GetYaxis()->SetBinLabel(i + 1, result.first.data());
239  for (std::size_t j = 0; j < result.second.size(); ++j) {
240  histo_status->SetBinContent(j + 1, i + 1, result.second[j]);
241  }
242  i++;
243  }
244 }
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:106

◆ getIndexFromString()

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

Definition at line 510 of file MillePedeDQMModule.cc.

References Exception.

Referenced by fillExpertHistos(), and fillExpertHistos_HG().

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

◆ setupChanged()

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

Definition at line 497 of file MillePedeDQMModule.cc.

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

Referenced by beginRun().

497  {
498  bool changed{false};
499 
500  if (watchIdealGeometryRcd_.check(setup))
501  changed = true;
502  if (watchTrackerTopologyRcd_.check(setup))
503  changed = true;
504  if (watchPTrackerParametersRcd_.check(setup))
505  changed = true;
506 
507  return changed;
508 }
MillePedeDQMModule::watchIdealGeometryRcd_
edm::ESWatcher< IdealGeometryRecord > watchIdealGeometryRcd_
Definition: MillePedeDQMModule.h:112
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:111
MillePedeDQMModule::watchPTrackerParametersRcd_
edm::ESWatcher< PTrackerParametersRcd > watchPTrackerParametersRcd_
Definition: MillePedeDQMModule.h:113

Member Data Documentation

◆ aliThrToken_

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

Definition at line 100 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ binariesAvalaible

MonitorElement* MillePedeDQMModule::binariesAvalaible
private

Definition at line 131 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and dqmEndJob().

◆ exitCode

MonitorElement* MillePedeDQMModule::exitCode
private

Definition at line 132 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and dqmEndJob().

◆ gDetToken_

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

Definition at line 98 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ geomToken_

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

Definition at line 101 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ h_xPos

MonitorElement* MillePedeDQMModule::h_xPos
private

Definition at line 116 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 123 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos_HG().

◆ h_xRot

MonitorElement* MillePedeDQMModule::h_xRot
private

Definition at line 117 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 124 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos_HG().

◆ h_yPos

MonitorElement* MillePedeDQMModule::h_yPos
private

Definition at line 118 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 125 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos_HG().

◆ h_yRot

MonitorElement* MillePedeDQMModule::h_yRot
private

Definition at line 119 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 126 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos_HG().

◆ h_zPos

MonitorElement* MillePedeDQMModule::h_zPos
private

Definition at line 120 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 127 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos_HG().

◆ h_zRot

MonitorElement* MillePedeDQMModule::h_zRot
private

Definition at line 121 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 128 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHistos_HG().

◆ isHG_

bool MillePedeDQMModule::isHG_
private

Definition at line 135 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and dqmEndJob().

◆ isVetoed

MonitorElement* MillePedeDQMModule::isVetoed
private

Definition at line 133 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and dqmEndJob().

◆ layerVec

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

Definition at line 109 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and fillExpertHisto_HG().

◆ mpReader_

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

Definition at line 106 of file MillePedeDQMModule.h.

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

◆ mpReaderConfig_

const edm::ParameterSet MillePedeDQMModule::mpReaderConfig_
private

Definition at line 104 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ outputFolder_

const std::string MillePedeDQMModule::outputFolder_
private

Definition at line 103 of file MillePedeDQMModule.h.

Referenced by bookHistograms().

◆ pixelTopologyMap_

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

Definition at line 107 of file MillePedeDQMModule.h.

Referenced by beginRun(), and bookHistograms().

◆ ptpToken_

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

Definition at line 99 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ statusResults

MonitorElement* MillePedeDQMModule::statusResults
private

Definition at line 130 of file MillePedeDQMModule.h.

Referenced by bookHistograms(), and dqmEndJob().

◆ tracker_

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

Definition at line 105 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ tTopoToken_

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

Definition at line 97 of file MillePedeDQMModule.h.

Referenced by beginRun().

◆ watchIdealGeometryRcd_

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

Definition at line 112 of file MillePedeDQMModule.h.

Referenced by setupChanged().

◆ watchPTrackerParametersRcd_

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

Definition at line 113 of file MillePedeDQMModule.h.

Referenced by setupChanged().

◆ watchTrackerTopologyRcd_

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

Definition at line 111 of file MillePedeDQMModule.h.

Referenced by setupChanged().

Generated for CMSSW Reference Manual by doxygen 1.8.14