#include <PixelLumiDQM.h>

Inheritance diagram for PixelLumiDQM:

Classes
class	PixelClusterCount

Public Member Functions
	PixelLumiDQM (const edm::ParameterSet &)

	~PixelLumiDQM () override

Public Member Functions inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) override

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

	DQMOneEDAnalyzer ()

void	endRun (edm::Run const &, edm::EventSetup const &) final

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... >
	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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Static Public Attributes
static constexpr double	CM2_TO_NANOBARN = 1.0 / 1.e-33

static constexpr double	FREQ_ORBIT = 11245.5

static const unsigned int	lastBunchCrossing = 3564

static constexpr double	rXSEC_PIXEL_CLUSTER = 9.4e-24

static constexpr double	rXSEC_PIXEL_CLUSTER_UNC = 0.119e-24

static constexpr double	SECONDS_PER_LS = double(0x40000) / double(FREQ_ORBIT)

static constexpr double	XSEC_PIXEL_CLUSTER = 10.08e-24

static constexpr double	XSEC_PIXEL_CLUSTER_UNC = 0.17e-24

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

unsigned int	calculateBunchMask (MonitorElement *, std::vector< bool > &)

unsigned int	calculateBunchMask (std::vector< float > &, unsigned int, std::vector< bool > &)

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

Private Attributes
std::vector< bool >	bunchTriggerMask

UInt_t	fBXNo

Double_t	fChi2

std::vector< uint32_t >	fDeadModules

UInt_t	fEvtNo

unsigned int	fFillNumber

MonitorElement *	fHistBunchCrossings

MonitorElement *	fHistBunchCrossingsLastLumi

MonitorElement *	fHistClusByLS

MonitorElement *	fHistClusterCountByBxCumulative

MonitorElement *	fHistClusterCountByBxLastLumi

std::map< std::string, MonitorElement * >	fHistContainerThisRun

MonitorElement *	fHistnBClusVsLS [3]

MonitorElement *	fHistnFMClusVsLS [2]

MonitorElement *	fHistnFPClusVsLS [2]

MonitorElement *	fHistRecordedByBxCumulative

MonitorElement *	fHistRecordedByBxLastLumi

MonitorElement *	fHistTotalRecordedLumiByLS

unsigned int	filledAndUnmaskedBunches

bool	fIncludePixelClusterInfo

bool	fIncludePixelQualCheckHistos

bool	fIncludeVertexInfo

MonitorElement *	fIntActiveCrossingsFromDB

std::string	fLogFileName_

UInt_t	fLSNo

double	fMinClusterCharge

int	fMinPixelsPerCluster

Double_t	fNdof

std::map< int, PixelClusterCount >	fNumPixelClusters

UInt_t	fNumTrkPerVtx

UInt_t	fNumVtx

UInt_t	fNumVtxGood

edm::EDGetTokenT< edmNew::DetSetVector< SiPixelCluster > >	fPixelClusterLabel

int	fResetIntervalInLumiSections

UInt_t	fRunNo

UInt_t	fTimestamp

Double_t	fVtxX

Double_t	fVtxY

Double_t	fVtxZ

std::ofstream	logFile_

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tkGeomToken_

bool	useInnerBarrelLayer

Static Private Attributes
static constexpr size_t	kNumDisks = 12

static constexpr size_t	kNumLayers = 5

static constexpr size_t	kOffsetDisks = 4

static constexpr size_t	kOffsetLayers = 0

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::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

Protected Member Functions inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)

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 DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

PixelLumi/PixelLumiDQM/plugins/PixelLumiDQM.h

Description: DQM Module producing Pixel Cluster Count Luminosity

Implementation notes: 1) Filling scheme is put in by 'hand' for now. Can obtain it from the cluster count but need higher rate in trigger; Best would be to have filling scheme in the DB. 2) Afterglow correction is put in by hand. 3) Currently barrel layer 0 is excluded, but a version using all pixel layers and disks is also in place. 4) A stable beam flag should be obtained from the DB to turn on actual checks. (Pixel Lumi does not make sense outside stable beams.) 5) Need a top module to correlate the info from the three trigger categories. NB: at present the module only uses the ZeroBias trigger providing about 85 Hz.

Definition at line 41 of file PixelLumiDQM.h.

Constructor & Destructor Documentation

◆ PixelLumiDQM()

PixelLumiDQM::PixelLumiDQM ( const edm::ParameterSet & iConfig )

explicit

Definition at line 47 of file PixelLumiDQM.cc.

References fDeadModules, fIncludePixelClusterInfo, fIncludePixelQualCheckHistos, fMinClusterCharge, fMinPixelsPerCluster, edm::ParameterSet::getUntrackedParameter(), ProducerED_cfi::InputTag, and ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it.

     : fPixelClusterLabel(consumes<edmNew::DetSetVector<SiPixelCluster>>(
           iConfig.getUntrackedParameter<edm::InputTag>("pixelClusterLabel", edm::InputTag("siPixelClusters")))),
       tkGeomToken_(esConsumes()),
       fIncludePixelClusterInfo(iConfig.getUntrackedParameter<bool>("includePixelClusterInfo", true)),
       fIncludePixelQualCheckHistos(iConfig.getUntrackedParameter<bool>("includePixelQualCheckHistos", true)),
       fResetIntervalInLumiSections(iConfig.getUntrackedParameter<int>("resetEveryNLumiSections", 1)),
       fDeadModules(iConfig.getUntrackedParameter<std::vector<uint32_t>>("deadModules", std::vector<uint32_t>())),
       fMinPixelsPerCluster(iConfig.getUntrackedParameter<int>("minNumPixelsPerCluster", 0)),
       fMinClusterCharge(iConfig.getUntrackedParameter<double>("minChargePerCluster", 0)),
       bunchTriggerMask(lastBunchCrossing + 1, false),
       filledAndUnmaskedBunches(0),
       useInnerBarrelLayer(iConfig.getUntrackedParameter<bool>("useInnerBarrelLayer", false)),
       fLogFileName_(iConfig.getUntrackedParameter<std::string>("logFileName", "/tmp/pixel_lumi.txt")) {
   edm::LogInfo("Configuration") << "PixelLumiDQM looking for pixel clusters in '"
                                 << iConfig.getUntrackedParameter<edm::InputTag>("pixelClusterLabel",
                                                                                 edm::InputTag("siPixelClusters"))
                                 << "'";
   edm::LogInfo("Configuration") << "PixelLumiDQM storing pixel cluster info? " << fIncludePixelClusterInfo;
   edm::LogInfo("Configuration") << "PixelLumiDQM storing pixel cluster quality check histograms? "
                                 << fIncludePixelQualCheckHistos;
 
   if (fDeadModules.empty()) {
     edm::LogInfo("Configuration") << "No pixel modules specified to be ignored";
   } else {
     edm::LogInfo("Configuration") << fDeadModules.size() << " pixel modules specified to be ignored:";
     for (std::vector<uint32_t>::const_iterator it = fDeadModules.begin(); it != fDeadModules.end(); ++it) {
       edm::LogInfo("Configuration") << "  " << *it;
     }
   }
   edm::LogInfo("Configuration") << "Ignoring pixel clusters with less than " << fMinPixelsPerCluster << " pixels";
   edm::LogInfo("Configuration") << "Ignoring pixel clusters with charge less than " << fMinClusterCharge;
 }

◆ ~PixelLumiDQM()

PixelLumiDQM::~PixelLumiDQM ( )

override

Definition at line 81 of file PixelLumiDQM.cc.

81 {}

Member Function Documentation

◆ analyze()

void PixelLumiDQM::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 89 of file PixelLumiDQM.cc.

References cms::cuda::assert(), edmNew::DetSetVector< T >::begin(), ALCARECOTkAlJpsiMuMu_cff::charge, hcalRecHitTable_cff::detId, TrackerGeometry::dets(), edmNew::DetSetVector< T >::end(), fBXNo, fDeadModules, fEvtNo, fHistBunchCrossings, fHistBunchCrossingsLastLumi, fHistContainerThisRun, dqm::impl::MonitorElement::Fill(), fIncludePixelClusterInfo, fIncludePixelQualCheckHistos, spr::find(), fLSNo, fMinClusterCharge, fMinPixelsPerCluster, fNumPixelClusters, fPixelClusterLabel, fRunNo, fTimestamp, edm::EventSetup::getData(), mps_fire::i, TrackerGeometry::idToDet(), iEvent, GeomDetEnumerators::isTrackerPixel(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, kNumDisks, kNumLayers, kOffsetDisks, kOffsetLayers, nano_mu_digi_cff::layer, PixelBarrelNameUpgrade::layerName(), Topology::localPosition(), PixelEndcapNameUpgrade::mI, PixelEndcapNameUpgrade::mO, PV3DBase< T, PVType, FrameType >::phi(), PixelEndcapNameUpgrade::pI, PixelSubdetector::PixelBarrel, DQMOffline_LumiMonitoring_cff::pixelClusters, PixelSubdetector::PixelEndcap, PixelEndcapNameUpgrade::pO, findQualityFiles::size, PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, GeomDet::surface(), tkGeomToken_, to_string(), Surface::toGlobal(), x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                               {
   // Collect all bookkeeping information.
   fRunNo = iEvent.id().run();
   fEvtNo = iEvent.id().event();
   fLSNo = iEvent.getLuminosityBlock().luminosityBlock();
   fBXNo = iEvent.bunchCrossing();
   fTimestamp = iEvent.time().unixTime();
   fHistBunchCrossings->Fill(float(fBXNo));
   fHistBunchCrossingsLastLumi->Fill(float(fBXNo));
   // This serves as event counter to compute luminosity from cluster counts.
   std::map<int, PixelClusterCount>::iterator it = fNumPixelClusters.find(fBXNo);
   if (it == fNumPixelClusters.end())
     fNumPixelClusters[fBXNo] = PixelClusterCount();
   // Find tracker geometry.
   const TrackerGeometry *trackerGeo = &iSetup.getData(tkGeomToken_);
   if (fIncludePixelClusterInfo) {
     // Find pixel clusters.
     edm::Handle<edmNew::DetSetVector<SiPixelCluster>> pixelClusters;
     iEvent.getByToken(fPixelClusterLabel, pixelClusters);
 
     // Loop over entire tracker geometry.
     for (TrackerGeometry::DetContainer::const_iterator i = trackerGeo->dets().begin(); i != trackerGeo->dets().end();
          ++i) {
       // See if this is a pixel unit(?).
 
       if (GeomDetEnumerators::isTrackerPixel((*i)->subDetector())) {
         DetId detId = (*i)->geographicalId();
         // Find all clusters on this detector module.
         edmNew::DetSetVector<SiPixelCluster>::const_iterator iSearch = pixelClusters->find(detId);
         if (iSearch != pixelClusters->end()) {
           // Count the number of clusters with at least a minimum
           // number of pixels per cluster and at least a minimum charge.
           size_t numClusters = 0;
           for (edmNew::DetSet<SiPixelCluster>::const_iterator itClus = iSearch->begin(); itClus != iSearch->end();
                ++itClus) {
             if ((itClus->size() >= fMinPixelsPerCluster) && (itClus->charge() >= fMinClusterCharge)) {
               ++numClusters;
             }
           }
           // DEBUG DEBUG DEBUG
           assert(numClusters <= iSearch->size());
           // DEBUG DEBUG DEBUG end
 
           // Add up the cluster count based on the position of this detector
           // element.
           if (detId.subdetId() == PixelSubdetector::PixelBarrel) {
             PixelBarrelNameUpgrade detName = PixelBarrelNameUpgrade(detId);
             int layer = detName.layerName() - kOffsetLayers;
             fNumPixelClusters[fBXNo].numB.at(layer) += numClusters;
             fNumPixelClusters[fBXNo].dnumB.at(layer) += sqrt(numClusters);
           } else {
             // DEBUG DEBUG DEBUG
             assert(detId.subdetId() == PixelSubdetector::PixelEndcap);
             // DEBUG DEBUG DEBUG end
 
             PixelEndcapNameUpgrade detName = PixelEndcapNameUpgrade(detId);
             PixelEndcapNameUpgrade::HalfCylinder halfCylinder = detName.halfCylinder();
             int disk = detName.diskName() - kOffsetDisks;
             switch (halfCylinder) {
               case PixelEndcapNameUpgrade::mO:
               case PixelEndcapNameUpgrade::mI:
                 fNumPixelClusters[fBXNo].numFM.at(disk) += numClusters;
                 fNumPixelClusters[fBXNo].dnumFM.at(disk) += sqrt(numClusters);
                 break;
               case PixelEndcapNameUpgrade::pO:
               case PixelEndcapNameUpgrade::pI:
                 fNumPixelClusters[fBXNo].numFP.at(disk) += numClusters;
                 fNumPixelClusters[fBXNo].dnumFP.at(disk) += sqrt(numClusters);
                 break;
               default:
                 assert(false);
                 break;
             }
           }
         }
       }
     }
   }
   // ----------
 
   // Fill some pixel cluster quality check histograms if requested.
   if (fIncludePixelQualCheckHistos) {
     // Find pixel clusters.
     edm::Handle<edmNew::DetSetVector<SiPixelCluster>> pixelClusters;
     iEvent.getByToken(fPixelClusterLabel, pixelClusters);
 
     bool filterDeadModules = (!fDeadModules.empty());
     std::vector<uint32_t>::const_iterator deadModulesBegin = fDeadModules.begin();
     std::vector<uint32_t>::const_iterator deadModulesEnd = fDeadModules.end();
 
     // Loop over entire tracker geometry.
     for (TrackerGeometry::DetContainer::const_iterator i = trackerGeo->dets().begin(); i != trackerGeo->dets().end();
          ++i) {
       // See if this is a pixel module.
       if (GeomDetEnumerators::isTrackerPixel((*i)->subDetector())) {
         DetId detId = (*i)->geographicalId();
 
         // Skip this module if it's on the list of modules to be ignored.
         if (filterDeadModules && find(deadModulesBegin, deadModulesEnd, detId()) != deadModulesEnd) {
           continue;
         }
 
         // Find all clusters in this module.
         edmNew::DetSetVector<SiPixelCluster>::const_iterator iSearch = pixelClusters->find(detId);
 
         // Loop over all clusters in this module.
         if (iSearch != pixelClusters->end()) {
           for (edmNew::DetSet<SiPixelCluster>::const_iterator clus = iSearch->begin(); clus != iSearch->end(); ++clus) {
             if ((clus->size() >= fMinPixelsPerCluster) && (clus->charge() >= fMinClusterCharge)) {
               PixelGeomDetUnit const *theGeomDet = dynamic_cast<PixelGeomDetUnit const *>(trackerGeo->idToDet(detId));
               PixelTopology const *topol = &(theGeomDet->specificTopology());
               double x = clus->x();
               double y = clus->y();
               LocalPoint clustLP = topol->localPosition(MeasurementPoint(x, y));
               GlobalPoint clustGP = theGeomDet->surface().toGlobal(clustLP);
               double charge = clus->charge() / 1.e3;
               int size = clus->size();
 
               if (detId.subdetId() == PixelSubdetector::PixelBarrel) {
                 PixelBarrelNameUpgrade detName = PixelBarrelNameUpgrade(detId);
                 unsigned int layer = detName.layerName() - kOffsetLayers;
                 if (layer < kNumLayers) {
                   std::string histName;
                   histName = "clusPosBarrel" + std::to_string(layer);
                   fHistContainerThisRun[histName]->Fill(clustGP.z(), clustGP.phi());
                   histName = "clusChargeBarrel" + std::to_string(layer);
                   fHistContainerThisRun[histName]->Fill(iEvent.bunchCrossing(), charge);
                   histName = "clusSizeBarrel" + std::to_string(layer);
                   fHistContainerThisRun[histName]->Fill(iEvent.bunchCrossing(), size);
                 } else {
                   edm::LogWarning("pixelLumi") << "higher layer number, " << layer << ", than layers";
                 }
               } else {
                 // DEBUG DEBUG DEBUG
                 assert(detId.subdetId() == PixelSubdetector::PixelEndcap);
                 // DEBUG DEBUG DEBUG end
 
                 PixelEndcapNameUpgrade detName = PixelEndcapNameUpgrade(detId);
                 unsigned int disk = detName.diskName() - kOffsetDisks;
                 if (disk < kNumDisks) {
                   std::string histName;
                   histName = "clusPosEndCap" + std::to_string(disk);
                   fHistContainerThisRun[histName]->Fill(clustGP.x(), clustGP.y());
                   histName = "clusChargeEndCap" + std::to_string(disk);
                   fHistContainerThisRun[histName]->Fill(iEvent.bunchCrossing(), charge);
                   histName = "clusSizeEndCap" + std::to_string(disk);
                   fHistContainerThisRun[histName]->Fill(iEvent.bunchCrossing(), size);
                 } else {
                   edm::LogWarning("pixelLumi")
                       << "higher disk number, " << disk << ", than disks," << kNumDisks << std::endl;
                 }
               }
             }
           }
         }
       }
     }
   }
 }

◆ beginLuminosityBlock()

void PixelLumiDQM::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumiBlock,
		edm::EventSetup const &
	)

overrideprivate

Definition at line 374 of file PixelLumiDQM.cc.

References fHistBunchCrossingsLastLumi, fHistClusterCountByBxLastLumi, fHistRecordedByBxLastLumi, fResetIntervalInLumiSections, eostools::ls(), edm::LuminosityBlockAuxiliary::luminosityBlock(), edm::LuminosityBlock::luminosityBlockAuxiliary(), and dqm::impl::MonitorElement::Reset().

                                                                                                   {
   // Only reset and fill every fResetIntervalInLumiSections (default is 1 LS)
   // Return unless the PREVIOUS LS was at the right modulo value
   // (e.g. is resetinterval = 5 the rest will only be executed at LS=6
   // NB: reset is done here so the histograms by LS are sent before resetting.
   // NB: not being used for now since default is 1 LS. There is a bug here.
 
   unsigned int ls = lumiBlock.luminosityBlockAuxiliary().luminosityBlock();
 
   if ((ls - 1) % fResetIntervalInLumiSections == 0) {
     fHistBunchCrossingsLastLumi->Reset();
     fHistClusterCountByBxLastLumi->Reset();
     fHistRecordedByBxLastLumi->Reset();
   }
 }

◆ bookHistograms()

void PixelLumiDQM::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	run,
		edm::EventSetup const &
	)

overrideprivatevirtual

Implements DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 249 of file PixelLumiDQM.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), fHistBunchCrossings, fHistBunchCrossingsLastLumi, fHistClusByLS, fHistClusterCountByBxCumulative, fHistClusterCountByBxLastLumi, fHistContainerThisRun, fHistnBClusVsLS, fHistnFMClusVsLS, fHistnFPClusVsLS, fHistRecordedByBxCumulative, fHistRecordedByBxLastLumi, fHistTotalRecordedLumiByLS, fIncludePixelQualCheckHistos, printsummarytable::folder, mps_fire::i, submitPVResolutionJobs::key, kNumDisks, kNumLayers, lastBunchCrossing, Skims_PA_cff::name, Geom::pi(), writedatasetfile::run, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and runGCPTkAlMap::title.

                                                            {
   edm::LogInfo("Status") << "Starting processing of run #" << run.id().run();
 
   // Top folder containing high-level information about pixel and HF lumi.
   std::string folder = "PixelLumi/";
   ibooker.setCurrentFolder(folder);
 
   fHistTotalRecordedLumiByLS = ibooker.book1D("totalPixelLumiByLS", "Pixel Lumi in nb vs LS", 8000, 0.5, 8000.5);
   fHistRecordedByBxCumulative = ibooker.book1D("PXLumiByBXsum",
                                                "Pixel Lumi in nb by BX Cumulative vs LS",
                                                lastBunchCrossing,
                                                0.5,
                                                float(lastBunchCrossing) + 0.5);
 
   std::string subfolder = folder + "lastLS/";
   ibooker.setCurrentFolder(subfolder);
   fHistRecordedByBxLastLumi = ibooker.book1D(
       "PXByBXLastLumi", "Pixel By BX Last Lumi", lastBunchCrossing + 1, -0.5, float(lastBunchCrossing) + 0.5);
 
   subfolder = folder + "ClusterCountingDetails/";
   ibooker.setCurrentFolder(subfolder);
 
   fHistnBClusVsLS[0] = ibooker.book1D("nBClusVsLS_0", "Fraction of Clusters vs LS Barrel layer 0", 8000, 0.5, 8000.5);
   fHistnBClusVsLS[1] = ibooker.book1D("nBClusVsLS_1", "Fraction of Clusters vs LS Barrel layer 1", 8000, 0.5, 8000.5);
   fHistnBClusVsLS[2] = ibooker.book1D("nBClusVsLS_2", "Fraction of Clusters vs LS Barrel layer 2", 8000, 0.5, 8000.5);
   fHistnFPClusVsLS[0] = ibooker.book1D("nFPClusVsLS_0", "Fraction of Clusters vs LS Barrel layer 0", 8000, 0.5, 8000.5);
   fHistnFPClusVsLS[1] = ibooker.book1D("nFPClusVsLS_1", "Fraction of Clusters vs LS Barrel layer 1", 8000, 0.5, 8000.5);
   fHistnFMClusVsLS[0] = ibooker.book1D("nFMClusVsLS_0", "Fraction of Clusters vs LS Barrel layer 0", 8000, 0.5, 8000.5);
   fHistnFMClusVsLS[1] = ibooker.book1D("nFMClusVsLS_1", "Fraction of Clusters vs LS Barrel layer 1", 8000, 0.5, 8000.5);
   fHistBunchCrossings = ibooker.book1D(
       "BunchCrossings", "Cumulative Bunch Crossings", lastBunchCrossing, 0.5, float(lastBunchCrossing) + 0.5);
   fHistBunchCrossingsLastLumi = ibooker.book1D(
       "BunchCrossingsLL", "Bunch Crossings Last Lumi", lastBunchCrossing, 0.5, float(lastBunchCrossing) + 0.5);
   fHistClusterCountByBxLastLumi = ibooker.book1D(
       "ClusterCountByBxLL", "Cluster Count by BX Last Lumi", lastBunchCrossing, 0.5, float(lastBunchCrossing) + 0.5);
   fHistClusterCountByBxCumulative = ibooker.book1D(
       "ClusterCountByBxSum", "Cluster Count by BX Cumulative", lastBunchCrossing, 0.5, float(lastBunchCrossing) + 0.5);
   fHistClusByLS = ibooker.book1D("totalClusByLS", "Number of Clusters all dets vs LS", 8000, 0.5, 8000.5);
 
   // Add some pixel cluster quality check histograms (in a subfolder).
   subfolder = folder + "qualityChecks/";
   ibooker.setCurrentFolder(subfolder);
 
   if (fIncludePixelQualCheckHistos) {
     // Create histograms for this run if not already present in our list.
     edm::LogInfo("Status") << "Creating histograms for run #" << run.id().run();
 
     // Pixel cluster positions in the barrel - (z, phi).
     for (size_t i = 0; i <= kNumLayers; ++i) {
       std::stringstream key;
       key << "clusPosBarrel" << i;
       std::stringstream name;
       name << key.str() << "_" << run.run();
       std::stringstream title;
       title << "Pixel cluster position - barrel layer " << i;
       fHistContainerThisRun[key.str()] =
           ibooker.book2D(name.str().c_str(), title.str().c_str(), 100, -30., 30., 64, -Geom::pi(), Geom::pi());
     }
 
     // Pixel cluster positions in the endcaps (x, y).
     for (size_t i = 0; i <= kNumDisks; ++i) {
       std::stringstream key;
       key << "clusPosEndCap" << i;
       std::stringstream name;
       name << key.str() << "_" << run.run();
       std::stringstream title;
       title << "Pixel cluster position - endcap disk " << i;
       fHistContainerThisRun[key.str()] =
           ibooker.book2D(name.str().c_str(), title.str().c_str(), 100, -20., 20., 100, -20., 20.);
     }
 
     // Pixel cluster charge in the barrel, per bx.
     for (size_t i = 0; i <= kNumLayers; ++i) {
       std::stringstream key;
       key << "clusChargeBarrel" << i;
       std::stringstream name;
       name << key.str() << "_" << run.run();
       std::stringstream title;
       title << "Pixel cluster charge - barrel layer " << i;
       fHistContainerThisRun[key.str()] =
           ibooker.book2D(name.str().c_str(), title.str().c_str(), 3564, .5, 3564.5, 100, 0., 100.);
     }
 
     // Pixel cluster charge in the endcaps, per bx.
     for (size_t i = 0; i <= kNumDisks; ++i) {
       std::stringstream key;
       key << "clusChargeEndCap" << i;
       std::stringstream name;
       name << key.str() << "_" << run.run();
       std::stringstream title;
       title << "Pixel cluster charge - endcap disk " << i;
       fHistContainerThisRun[key.str()] =
           ibooker.book2D(name.str().c_str(), title.str().c_str(), 3564, .5, 3564.5, 100, 0., 100.);
     }
 
     // Pixel cluster size in the barrel, per bx.
     for (size_t i = 0; i <= kNumLayers; ++i) {
       std::stringstream key;
       key << "clusSizeBarrel" << i;
       std::stringstream name;
       name << key.str() << "_" << run.run();
       std::stringstream title;
       title << "Pixel cluster size - barrel layer " << i;
       fHistContainerThisRun[key.str()] =
           ibooker.book2D(name.str().c_str(), title.str().c_str(), 3564, .5, 3564.5, 100, 0., 100.);
     }
 
     // Pixel cluster size in the endcaps, per bx.
     for (size_t i = 0; i <= kNumDisks; ++i) {
       std::stringstream key;
       key << "clusSizeEndCap" << i;
       std::stringstream name;
       name << key.str() << "_" << run.run();
       std::stringstream title;
       title << "Pixel cluster size - endcap disk " << i;
       fHistContainerThisRun[key.str()] =
           ibooker.book2D(name.str().c_str(), title.str().c_str(), 3564, .5, 3564.5, 100, 0., 100.);
     }
   }
 }

◆ calculateBunchMask() [1/2]

unsigned int PixelLumiDQM::calculateBunchMask	(	MonitorElement *	e,
		std::vector< bool > &	mask
	)

private

Definition at line 549 of file PixelLumiDQM.cc.

References MillePedeFileConverter_cfg::e, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::mask, and LaserClient_cfi::nbins.

Referenced by endLuminosityBlock().

                                                                                       {
   unsigned int nbins = e->getNbinsX();
   std::vector<float> ar(nbins + 1, 0.);
   for (unsigned int i = 1; i <= nbins; i++) {
     ar[i] = e->getBinContent(i);
   }
   return calculateBunchMask(ar, nbins, mask);
 }

◆ calculateBunchMask() [2/2]

unsigned int PixelLumiDQM::calculateBunchMask	(	std::vector< float > &	e,
		unsigned int	nbins,
		std::vector< bool > &	mask
	)

private

Definition at line 557 of file PixelLumiDQM.cc.

References funct::abs(), newFWLiteAna::bin, MillePedeFileConverter_cfg::e, l1tPhase2CaloJetEmulator_cfi::fit, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::mask, SiStripPI::mean, and LaserClient_cfi::nbins.

                                                                                                             {
   // Take the cumulative cluster count histogram and find max and average of
   // non-empty bins.
   unsigned int active_count = 0;
   double maxc = 0.0;
   double ave = 0.0;  // Average of non-empty bins
   unsigned int non_empty_bins = 0;
 
   for (unsigned int i = 1; i <= nbins; i++) {
     double bin = e[i];
     if (bin != 0.0) {
       if (maxc < bin)
         maxc = bin;
       ave += bin;
       non_empty_bins++;
     }
   }
 
   ave /= non_empty_bins;
   edm::LogWarning("pixelLumi") << "Bunch mask finder - non empty bins " << non_empty_bins
                                << " average of non empty bins " << ave << " max content of one bin " << maxc;
   double mean = 0.;
   double sigma = 0.;
   if (non_empty_bins < 50) {
     mean = maxc;
     sigma = (maxc) / 20;
   } else {
     TH1F dist("dist", "dist", 500, 0., maxc + (maxc / 500.) * 20.);
     for (unsigned int i = 1; i <= nbins; i++) {
       double bin = e[i];
       dist.Fill(bin);
     }
     TF1 fit("plgaus", "gaus");
     dist.Fit(&fit, "", "", fmax(0., ave - (maxc - ave) / 5.), maxc);
     mean = fit.GetParameter("Mean");
     sigma = fit.GetParameter("Sigma");
   }
   edm::LogWarning("pixelLumi") << "Bunch mask will use mean" << mean << " sigma " << sigma;
   // Active BX defined as those which have nclus within fixed standard
   // deviations of peak.
   for (unsigned int i = 1; i <= nbins; i++) {
     double bin = e[i];
     if (bin > 0. && std::abs(bin - mean) < 5. * (sigma)) {
       mask[i] = true;
       active_count++;
     }
   }
   edm::LogWarning("pixelLumi") << "Bunch mask finds " << active_count << " active bunch crossings ";
   //   edm::LogWarning("pixelLumi") << "this is the full bx mask " ;
   //   for(unsigned int i = 1; i<= nbins; i++)
   //     edm::LogWarning("pixelLumi") << ((mask[i]) ? 1:0);
   return active_count;
 }

◆ endLuminosityBlock()

void PixelLumiDQM::endLuminosityBlock	(	edm::LuminosityBlock const &	lumiBlock,
		edm::EventSetup const &	es
	)

overrideprivate

Definition at line 392 of file PixelLumiDQM.cc.

References bunchTriggerMask, calculateBunchMask(), CM2_TO_NANOBARN, submitPVResolutionJobs::err, fHistBunchCrossingsLastLumi, fHistClusByLS, fHistClusterCountByBxCumulative, fHistClusterCountByBxLastLumi, fHistnBClusVsLS, fHistnFMClusVsLS, fHistnFPClusVsLS, fHistRecordedByBxCumulative, fHistRecordedByBxLastLumi, fHistTotalRecordedLumiByLS, filledAndUnmaskedBunches, fLogFileName_, fNumPixelClusters, FREQ_ORBIT, fResetIntervalInLumiSections, fRunNo, dqm::impl::MonitorElement::getBinContent(), dqm::impl::MonitorElement::getBinError(), mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, lastBunchCrossing, logFile_, eostools::ls(), edm::LuminosityBlockAuxiliary::luminosityBlock(), edm::LuminosityBlock::luminosityBlockAuxiliary(), rXSEC_PIXEL_CLUSTER, rXSEC_PIXEL_CLUSTER_UNC, SECONDS_PER_LS, dqm::impl::MonitorElement::setBinContent(), dqm::impl::MonitorElement::setBinError(), mathSSE::sqrt(), dqmMemoryStats::total, pileupReCalc_HLTpaths::trunc, useInnerBarrelLayer, XSEC_PIXEL_CLUSTER, and XSEC_PIXEL_CLUSTER_UNC.

                                                                                                   {
   unsigned int ls = lumiBlock.luminosityBlockAuxiliary().luminosityBlock();
 
   // Only fill every fResetIntervalInLumiSections (default is 1 LS)
   if (ls % fResetIntervalInLumiSections != 0)
     return;
 
   printf("Lumi Block = %d\n", ls);
 
   if ((ls - 1) % fResetIntervalInLumiSections == 0) {
   }
 
   unsigned int nBClus[3] = {0, 0, 0};
   unsigned int nFPClus[2] = {0, 0};
   unsigned int nFMClus[2] = {0, 0};
 
   double total_recorded = 0.;
   double total_recorded_unc_square = 0.;
 
   // Obtain bunch-by-bunch cluster counts and compute totals for lumi
   // calculation.
   double totalcounts = 0.0;
   double totalevents = 0.0;
   double lumi_factor_per_bx = 0.0;
   if (useInnerBarrelLayer)
     lumi_factor_per_bx = FREQ_ORBIT * SECONDS_PER_LS * fResetIntervalInLumiSections / XSEC_PIXEL_CLUSTER;
   else
     lumi_factor_per_bx = FREQ_ORBIT * SECONDS_PER_LS * fResetIntervalInLumiSections / rXSEC_PIXEL_CLUSTER;
 
   for (std::map<int, PixelClusterCount>::iterator it = fNumPixelClusters.begin(); it != fNumPixelClusters.end(); it++) {
     // Sum all clusters for this BX.
     unsigned int total = (*it).second.numB.at(1) + (*it).second.numB.at(2) + (*it).second.numFP.at(0) +
                          (*it).second.numFP.at(1) + (*it).second.numFM.at(0) + (*it).second.numFM.at(1);
     if (useInnerBarrelLayer)
       total += (*it).second.numB.at(0);
     totalcounts += total;
     double etotal = useInnerBarrelLayer
                         ? (*it).second.dnumB.at(0)
                         : (*it).second.dnumB.at(1) + (*it).second.dnumB.at(2) + (*it).second.dnumFP.at(0) +
                               (*it).second.dnumFP.at(1) + (*it).second.dnumFM.at(0) + (*it).second.dnumFM.at(1);
     etotal = sqrt(etotal);
 
     fHistClusterCountByBxLastLumi->setBinContent((*it).first, total);
     fHistClusterCountByBxLastLumi->setBinError((*it).first, etotal);
     fHistClusterCountByBxCumulative->setBinContent((*it).first,
                                                    fHistClusterCountByBxCumulative->getBinContent((*it).first) + total);
 
     unsigned int events_per_bx = fHistBunchCrossingsLastLumi->getBinContent((*it).first);
     totalevents += events_per_bx;
     double average_cluster_count = events_per_bx != 0 ? double(total) / events_per_bx : 0.;
     double average_cluster_count_unc = events_per_bx != 0 ? etotal / events_per_bx : 0.;
     double pixel_bx_lumi_per_ls = lumi_factor_per_bx * average_cluster_count / CM2_TO_NANOBARN;
     double pixel_bx_lumi_per_ls_unc = 0.0;
     if (useInnerBarrelLayer)
       pixel_bx_lumi_per_ls_unc = sqrt(lumi_factor_per_bx * lumi_factor_per_bx *
                                       (average_cluster_count_unc * average_cluster_count_unc +
                                        (average_cluster_count * XSEC_PIXEL_CLUSTER_UNC / XSEC_PIXEL_CLUSTER) *
                                            (average_cluster_count * XSEC_PIXEL_CLUSTER / XSEC_PIXEL_CLUSTER))) /
                                  CM2_TO_NANOBARN;
     else
       pixel_bx_lumi_per_ls_unc = sqrt(lumi_factor_per_bx * lumi_factor_per_bx *
                                       (average_cluster_count_unc * average_cluster_count_unc +
                                        (average_cluster_count * rXSEC_PIXEL_CLUSTER_UNC / rXSEC_PIXEL_CLUSTER) *
                                            (average_cluster_count * rXSEC_PIXEL_CLUSTER / rXSEC_PIXEL_CLUSTER))) /
                                  CM2_TO_NANOBARN;
 
     fHistRecordedByBxLastLumi->setBinContent((*it).first, pixel_bx_lumi_per_ls);
     fHistRecordedByBxLastLumi->setBinError((*it).first, pixel_bx_lumi_per_ls_unc);
 
     fHistRecordedByBxCumulative->setBinContent(
         (*it).first, fHistRecordedByBxCumulative->getBinContent((*it).first) + pixel_bx_lumi_per_ls);
 
     /*
       if(fHistRecordedByBxLastLumi->getBinContent((*it).first)!=0.)
       fHistRecordedByBxLastLumi->getBinContent((*it).first));
       if(fHistRecordedByBxCumulative->getBinContent((*it).first)!=0.)
       fHistRecordedByBxCumulative->getBinContent((*it).first));
     */
 
     nBClus[0] += (*it).second.numB.at(0);
     nBClus[1] += (*it).second.numB.at(1);
     nBClus[2] += (*it).second.numB.at(2);
     nFPClus[0] += (*it).second.numFP.at(0);
     nFPClus[1] += (*it).second.numFP.at(1);
     nFMClus[0] += (*it).second.numFM.at(0);
     nFMClus[1] += (*it).second.numFM.at(1);
 
     // Reset counters
     (*it).second.Reset();
 
     // edm::LogWarning("pixelLumi") << "bx="<< (*it).first << " clusters=" <<
     // (*it).second.numB.at(0));
   }
 
   if ((filledAndUnmaskedBunches = calculateBunchMask(fHistClusterCountByBxCumulative, bunchTriggerMask)) != 0) {
     for (unsigned int i = 0; i <= lastBunchCrossing; i++) {
       if (bunchTriggerMask[i]) {
         double err = fHistRecordedByBxLastLumi->getBinError(i);
         total_recorded += fHistRecordedByBxLastLumi->getBinContent(i);
         total_recorded_unc_square += err * err;
       }
     }
 
     // Replace the total obtained by summing over BXs with the average per BX
     // from the total cluster count and rescale
 
     if (totalevents > 10) {
       total_recorded = lumi_factor_per_bx * totalcounts / totalevents / CM2_TO_NANOBARN;
     } else
       total_recorded = 0.0;
 
     edm::LogWarning("pixelLumi") << " Total recorded " << total_recorded;
     fHistTotalRecordedLumiByLS->setBinContent(ls, total_recorded);
     fHistTotalRecordedLumiByLS->setBinError(ls, sqrt(total_recorded_unc_square));
   }
   // fill cluster counts by detector regions for sanity checks
   unsigned int all_detectors_counts = 0;
   for (unsigned int i = 0; i < 3; i++) {
     all_detectors_counts += nBClus[i];
   }
   for (unsigned int i = 0; i < 2; i++) {
     all_detectors_counts += nFPClus[i];
   }
   for (unsigned int i = 0; i < 2; i++) {
     all_detectors_counts += nFMClus[i];
   }
 
   fHistClusByLS->setBinContent(ls, all_detectors_counts);
 
   for (unsigned int i = 0; i < 3; i++) {
     fHistnBClusVsLS[i]->setBinContent(ls, float(nBClus[i]) / float(all_detectors_counts));
   }
   for (unsigned int i = 0; i < 2; i++) {
     fHistnFPClusVsLS[i]->setBinContent(ls, float(nFPClus[i]) / float(all_detectors_counts));
   }
   for (unsigned int i = 0; i < 2; i++) {
     fHistnFMClusVsLS[i]->setBinContent(ls, float(nFMClus[i]) / float(all_detectors_counts));
   }
 
   logFile_.open(fLogFileName_.c_str(), std::ios_base::trunc);
 
   timeval tv;
   gettimeofday(&tv, nullptr);
   tm *ts = gmtime(&tv.tv_sec);
   char datestring[256];
   strftime(datestring, sizeof(datestring), "%Y.%m.%d %T GMT %s", ts);
   logFile_ << "RunNumber " << fRunNo << std::endl;
   logFile_ << "EndTimeOfFit " << datestring << std::endl;
   logFile_ << "LumiRange " << ls << "-" << ls << std::endl;
   logFile_ << "Fill " << -99 << std::endl;
   logFile_ << "ActiveBunchCrossings " << filledAndUnmaskedBunches << std::endl;
   logFile_ << "PixelLumi " << fHistTotalRecordedLumiByLS->getBinContent(ls) * 0.98 << std::endl;
   logFile_ << "HFLumi " << -99 << std::endl;
   logFile_ << "Ratio " << -99 << std::endl;
   logFile_.close();
 }

◆ fillDescriptions()

void PixelLumiDQM::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 83 of file PixelLumiDQM.cc.

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

                                                                               {
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

Member Data Documentation

◆ bunchTriggerMask

std::vector<bool> PixelLumiDQM::bunchTriggerMask

private

Definition at line 167 of file PixelLumiDQM.h.

Referenced by endLuminosityBlock().

◆ CM2_TO_NANOBARN

constexpr double PixelLumiDQM::CM2_TO_NANOBARN = 1.0 / 1.e-33

static

Definition at line 55 of file PixelLumiDQM.h.

Referenced by endLuminosityBlock().

◆ fBXNo

UInt_t PixelLumiDQM::fBXNo

private

Definition at line 121 of file PixelLumiDQM.h.

Referenced by analyze().

◆ fChi2

Double_t PixelLumiDQM::fChi2

private

Definition at line 129 of file PixelLumiDQM.h.

◆ fDeadModules

std::vector<uint32_t> PixelLumiDQM::fDeadModules

private

Definition at line 142 of file PixelLumiDQM.h.

Referenced by analyze(), and PixelLumiDQM().

◆ fEvtNo

UInt_t PixelLumiDQM::fEvtNo

private

Definition at line 119 of file PixelLumiDQM.h.

Referenced by analyze().

◆ fFillNumber

unsigned int PixelLumiDQM::fFillNumber

private

Definition at line 170 of file PixelLumiDQM.h.

◆ fHistBunchCrossings

MonitorElement* PixelLumiDQM::fHistBunchCrossings

private

Definition at line 158 of file PixelLumiDQM.h.

Referenced by analyze(), and bookHistograms().

◆ fHistBunchCrossingsLastLumi

MonitorElement* PixelLumiDQM::fHistBunchCrossingsLastLumi

private

Definition at line 159 of file PixelLumiDQM.h.

Referenced by analyze(), beginLuminosityBlock(), bookHistograms(), and endLuminosityBlock().

◆ fHistClusByLS

MonitorElement* PixelLumiDQM::fHistClusByLS

private

Definition at line 162 of file PixelLumiDQM.h.

Referenced by bookHistograms(), and endLuminosityBlock().

◆ fHistClusterCountByBxCumulative

MonitorElement* PixelLumiDQM::fHistClusterCountByBxCumulative

private

Definition at line 161 of file PixelLumiDQM.h.

Referenced by bookHistograms(), and endLuminosityBlock().

◆ fHistClusterCountByBxLastLumi

MonitorElement* PixelLumiDQM::fHistClusterCountByBxLastLumi

private

Definition at line 160 of file PixelLumiDQM.h.

Referenced by beginLuminosityBlock(), bookHistograms(), and endLuminosityBlock().

◆ fHistContainerThisRun

std::map<std::string, MonitorElement *> PixelLumiDQM::fHistContainerThisRun

private

Definition at line 138 of file PixelLumiDQM.h.

Referenced by analyze(), and bookHistograms().

◆ fHistnBClusVsLS

MonitorElement* PixelLumiDQM::fHistnBClusVsLS[3]

private

Definition at line 155 of file PixelLumiDQM.h.

Referenced by bookHistograms(), and endLuminosityBlock().

◆ fHistnFMClusVsLS

MonitorElement* PixelLumiDQM::fHistnFMClusVsLS[2]

private

Definition at line 157 of file PixelLumiDQM.h.

Referenced by bookHistograms(), and endLuminosityBlock().

◆ fHistnFPClusVsLS

MonitorElement* PixelLumiDQM::fHistnFPClusVsLS[2]

private

Definition at line 156 of file PixelLumiDQM.h.

Referenced by bookHistograms(), and endLuminosityBlock().

◆ fHistRecordedByBxCumulative

MonitorElement* PixelLumiDQM::fHistRecordedByBxCumulative

private

Definition at line 165 of file PixelLumiDQM.h.

Referenced by bookHistograms(), and endLuminosityBlock().

◆ fHistRecordedByBxLastLumi

MonitorElement* PixelLumiDQM::fHistRecordedByBxLastLumi

private

Definition at line 164 of file PixelLumiDQM.h.

Referenced by beginLuminosityBlock(), bookHistograms(), and endLuminosityBlock().

◆ fHistTotalRecordedLumiByLS

MonitorElement* PixelLumiDQM::fHistTotalRecordedLumiByLS

private

Definition at line 163 of file PixelLumiDQM.h.

Referenced by bookHistograms(), and endLuminosityBlock().

◆ filledAndUnmaskedBunches

unsigned int PixelLumiDQM::filledAndUnmaskedBunches

private

Definition at line 168 of file PixelLumiDQM.h.

Referenced by endLuminosityBlock().

◆ fIncludePixelClusterInfo

bool PixelLumiDQM::fIncludePixelClusterInfo

private

Definition at line 134 of file PixelLumiDQM.h.

Referenced by analyze(), and PixelLumiDQM().

◆ fIncludePixelQualCheckHistos

bool PixelLumiDQM::fIncludePixelQualCheckHistos

private

Definition at line 135 of file PixelLumiDQM.h.

Referenced by analyze(), bookHistograms(), and PixelLumiDQM().

◆ fIncludeVertexInfo

bool PixelLumiDQM::fIncludeVertexInfo

private

Definition at line 133 of file PixelLumiDQM.h.

◆ fIntActiveCrossingsFromDB

MonitorElement* PixelLumiDQM::fIntActiveCrossingsFromDB

private

Definition at line 154 of file PixelLumiDQM.h.

◆ fLogFileName_

std::string PixelLumiDQM::fLogFileName_

private

Definition at line 171 of file PixelLumiDQM.h.

Referenced by endLuminosityBlock().

◆ fLSNo

UInt_t PixelLumiDQM::fLSNo

private

Definition at line 120 of file PixelLumiDQM.h.

Referenced by analyze().

◆ fMinClusterCharge

double PixelLumiDQM::fMinClusterCharge

private

Definition at line 150 of file PixelLumiDQM.h.

Referenced by analyze(), and PixelLumiDQM().

◆ fMinPixelsPerCluster

int PixelLumiDQM::fMinPixelsPerCluster

private

Definition at line 147 of file PixelLumiDQM.h.

Referenced by analyze(), and PixelLumiDQM().

◆ fNdof

Double_t PixelLumiDQM::fNdof

private

Definition at line 130 of file PixelLumiDQM.h.

◆ fNumPixelClusters

std::map<int, PixelClusterCount> PixelLumiDQM::fNumPixelClusters

private

Definition at line 131 of file PixelLumiDQM.h.

Referenced by analyze(), and endLuminosityBlock().

◆ fNumTrkPerVtx

UInt_t PixelLumiDQM::fNumTrkPerVtx

private

Definition at line 125 of file PixelLumiDQM.h.

◆ fNumVtx

UInt_t PixelLumiDQM::fNumVtx

private

Definition at line 123 of file PixelLumiDQM.h.

◆ fNumVtxGood

UInt_t PixelLumiDQM::fNumVtxGood

private

Definition at line 124 of file PixelLumiDQM.h.

◆ fPixelClusterLabel

edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster> > PixelLumiDQM::fPixelClusterLabel

private

Definition at line 115 of file PixelLumiDQM.h.

Referenced by analyze().

◆ FREQ_ORBIT

constexpr double PixelLumiDQM::FREQ_ORBIT = 11245.5

static

Definition at line 45 of file PixelLumiDQM.h.

Referenced by endLuminosityBlock().

◆ fResetIntervalInLumiSections

int PixelLumiDQM::fResetIntervalInLumiSections

private

Definition at line 136 of file PixelLumiDQM.h.

Referenced by beginLuminosityBlock(), and endLuminosityBlock().

◆ fRunNo

UInt_t PixelLumiDQM::fRunNo

private

Definition at line 118 of file PixelLumiDQM.h.

Referenced by analyze(), and endLuminosityBlock().

◆ fTimestamp

UInt_t PixelLumiDQM::fTimestamp

private

Definition at line 122 of file PixelLumiDQM.h.

Referenced by analyze().

◆ fVtxX

Double_t PixelLumiDQM::fVtxX

private

Definition at line 126 of file PixelLumiDQM.h.

◆ fVtxY

Double_t PixelLumiDQM::fVtxY

private

Definition at line 127 of file PixelLumiDQM.h.

◆ fVtxZ

Double_t PixelLumiDQM::fVtxZ

private

Definition at line 128 of file PixelLumiDQM.h.

◆ kNumDisks

constexpr size_t PixelLumiDQM::kNumDisks = 12

staticprivate

Definition at line 76 of file PixelLumiDQM.h.

Referenced by analyze(), and bookHistograms().

◆ kNumLayers

constexpr size_t PixelLumiDQM::kNumLayers = 5

staticprivate

Definition at line 75 of file PixelLumiDQM.h.

Referenced by analyze(), and bookHistograms().

◆ kOffsetDisks

constexpr size_t PixelLumiDQM::kOffsetDisks = 4

staticprivate

Definition at line 78 of file PixelLumiDQM.h.

Referenced by analyze().

◆ kOffsetLayers

constexpr size_t PixelLumiDQM::kOffsetLayers = 0

staticprivate

Definition at line 77 of file PixelLumiDQM.h.

Referenced by analyze().

◆ lastBunchCrossing

const unsigned int PixelLumiDQM::lastBunchCrossing = 3564

static

Definition at line 56 of file PixelLumiDQM.h.

Referenced by bookHistograms(), and endLuminosityBlock().

◆ logFile_

std::ofstream PixelLumiDQM::logFile_

private

Definition at line 173 of file PixelLumiDQM.h.

Referenced by endLuminosityBlock().

◆ rXSEC_PIXEL_CLUSTER

constexpr double PixelLumiDQM::rXSEC_PIXEL_CLUSTER = 9.4e-24

static

Definition at line 53 of file PixelLumiDQM.h.

Referenced by endLuminosityBlock().

◆ rXSEC_PIXEL_CLUSTER_UNC

constexpr double PixelLumiDQM::rXSEC_PIXEL_CLUSTER_UNC = 0.119e-24

static

Definition at line 54 of file PixelLumiDQM.h.

Referenced by endLuminosityBlock().

◆ SECONDS_PER_LS

constexpr double PixelLumiDQM::SECONDS_PER_LS = double(0x40000) / double(FREQ_ORBIT)

static

Definition at line 46 of file PixelLumiDQM.h.

Referenced by endLuminosityBlock().

◆ tkGeomToken_

edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> PixelLumiDQM::tkGeomToken_

private

Definition at line 116 of file PixelLumiDQM.h.

Referenced by analyze().

◆ useInnerBarrelLayer

bool PixelLumiDQM::useInnerBarrelLayer

private

Definition at line 169 of file PixelLumiDQM.h.

Referenced by endLuminosityBlock().

◆ XSEC_PIXEL_CLUSTER

constexpr double PixelLumiDQM::XSEC_PIXEL_CLUSTER = 10.08e-24

static

Definition at line 49 of file PixelLumiDQM.h.

Referenced by endLuminosityBlock().

◆ XSEC_PIXEL_CLUSTER_UNC

constexpr double PixelLumiDQM::XSEC_PIXEL_CLUSTER_UNC = 0.17e-24

static

Definition at line 50 of file PixelLumiDQM.h.

Referenced by endLuminosityBlock().

Classes

Public Member Functions

Static Public Member Functions

Static Public Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PixelLumiDQM()

◆ ~PixelLumiDQM()

Member Function Documentation

◆ analyze()

◆ beginLuminosityBlock()

◆ bookHistograms()

◆ calculateBunchMask() [1/2]

◆ calculateBunchMask() [2/2]

◆ endLuminosityBlock()

◆ fillDescriptions()

Member Data Documentation

◆ bunchTriggerMask

◆ CM2_TO_NANOBARN

◆ fBXNo

◆ fChi2

◆ fDeadModules

◆ fEvtNo

◆ fFillNumber

◆ fHistBunchCrossings

◆ fHistBunchCrossingsLastLumi

◆ fHistClusByLS

◆ fHistClusterCountByBxCumulative

◆ fHistClusterCountByBxLastLumi

◆ fHistContainerThisRun

◆ fHistnBClusVsLS

◆ fHistnFMClusVsLS

◆ fHistnFPClusVsLS

◆ fHistRecordedByBxCumulative

◆ fHistRecordedByBxLastLumi

◆ fHistTotalRecordedLumiByLS

◆ filledAndUnmaskedBunches

◆ fIncludePixelClusterInfo

◆ fIncludePixelQualCheckHistos

◆ fIncludeVertexInfo

◆ fIntActiveCrossingsFromDB

◆ fLogFileName_

◆ fLSNo

◆ fMinClusterCharge

◆ fMinPixelsPerCluster

◆ fNdof

◆ fNumPixelClusters

◆ fNumTrkPerVtx

◆ fNumVtx

◆ fNumVtxGood

◆ fPixelClusterLabel

◆ FREQ_ORBIT

◆ fResetIntervalInLumiSections

◆ fRunNo

◆ fTimestamp

◆ fVtxX

◆ fVtxY

◆ fVtxZ

◆ kNumDisks

◆ kNumLayers

◆ kOffsetDisks

◆ kOffsetLayers

◆ lastBunchCrossing

◆ logFile_

◆ rXSEC_PIXEL_CLUSTER

◆ rXSEC_PIXEL_CLUSTER_UNC

◆ SECONDS_PER_LS

◆ tkGeomToken_

◆ useInnerBarrelLayer

◆ XSEC_PIXEL_CLUSTER

◆ XSEC_PIXEL_CLUSTER_UNC