Inheritance diagram for LumiMonitor:

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

	~LumiMonitor () override=default

Public Member Functions inherited from edm::global::EDAnalyzer< edm::RunCache< Histograms >, Args... >
	EDAnalyzer ()=default

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

	~EDAnalyzerBase () 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

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

virtual	~EDConsumerBase () noexcept(false)

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

static void	fillHistoLSPSetDescription (edm::ParameterSetDescription &pset)

static void	fillHistoPSetDescription (edm::ParameterSetDescription &pset)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	bookHistograms (DQMStore::ConcurrentBooker &, edm::Run const &, edm::EventSetup const &, Histograms &) const override

void	dqmAnalyze (edm::Event const &, edm::EventSetup const &, Histograms const &) const override

Static Private Member Functions
static MEbinning	getHistoLSPSet (const edm::ParameterSet &pset)

static MEbinning	getHistoPSet (const edm::ParameterSet &pset)

Private Attributes
bool	doPixelLumi_

std::string	folderName_

MEbinning	ls_binning_

MEbinning	lumi_binning_

float	lumi_factor_per_bx_

edm::EDGetTokenT< LumiScalersCollection >	lumiScalersToken_

edm::EDGetTokenT< LumiSummary >	lumiSummaryToken_

int	minNumberOfPixelsPerCluster_

float	minPixelClusterCharge_

MEbinning	pixelCluster_binning_

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

MEbinning	pixellumi_binning_

MEbinning	pu_binning_

bool	useBPixLayer1_

Additional Inherited Members
Public Types inherited from edm::global::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ESProduct , Transition Tr = Transition::Event>
auto	esConsumes (eventsetup::EventSetupRecordKey const &, ESInputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 45 of file LumiMonitor.cc.

Constructor & Destructor Documentation

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

Definition at line 88 of file LumiMonitor.cc.

References doPixelLumi_, GetLumi::FREQ_ORBIT, getHistoPSet(), edm::ParameterSet::getParameter(), lumi_factor_per_bx_, pixellumi_binning_, GetLumi::rXSEC_PIXEL_CLUSTER, GetLumi::SECONDS_PER_LS, useBPixLayer1_, and GetLumi::XSEC_PIXEL_CLUSTER.

                                                       :
   folderName_             ( config.getParameter<std::string>("FolderName") )
   , lumiScalersToken_     ( consumes<LumiScalersCollection>(config.getParameter<edm::InputTag>("scalers") ) )
   , lumi_binning_         ( getHistoPSet  (config.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("lumiPSet"))         )
   , pu_binning_           ( getHistoPSet  (config.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("puPSet"))           )
   , ls_binning_           ( getHistoLSPSet(config.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("lsPSet"))           )
   , doPixelLumi_          ( config.getParameter<bool>("doPixelLumi") )
   , pixelClustersToken_           ( doPixelLumi_ ? consumes<edmNew::DetSetVector<SiPixelCluster>>(config.getParameter<edm::InputTag>("pixelClusters") ) : edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster>> () )
   , useBPixLayer1_                ( doPixelLumi_ ? config.getParameter<bool>   ("useBPixLayer1") : false            )
   , minNumberOfPixelsPerCluster_  ( doPixelLumi_ ? config.getParameter<int>    ("minNumberOfPixelsPerCluster") : -1 )
   , minPixelClusterCharge_        ( doPixelLumi_ ? config.getParameter<double> ("minPixelClusterCharge") : -1.      )
   , pixelCluster_binning_         ( doPixelLumi_ ? getHistoPSet(config.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("pixelClusterPSet")) : MEbinning {} )
   , pixellumi_binning_            ( doPixelLumi_ ? getHistoPSet(config.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("pixellumiPSet")) : MEbinning {} )
 {
   if(useBPixLayer1_) {
     lumi_factor_per_bx_ = GetLumi::FREQ_ORBIT * GetLumi::SECONDS_PER_LS / GetLumi::XSEC_PIXEL_CLUSTER  ;
   } else {
     lumi_factor_per_bx_ = GetLumi::FREQ_ORBIT * GetLumi::SECONDS_PER_LS / GetLumi::rXSEC_PIXEL_CLUSTER  ;
   }
 }

LumiMonitor::~LumiMonitor ( )

overridedefault

Member Function Documentation

void LumiMonitor::bookHistograms	(	DQMStore::ConcurrentBooker &	booker,
		edm::Run const &	run,
		edm::EventSetup const &	setup,
		Histograms &	histograms
	)		const

overrideprivatevirtual

Implements DQMGlobalEDAnalyzer< Histograms >.

Definition at line 127 of file LumiMonitor.cc.

References DQMStore::ConcurrentBooker::book1D(), DQMStore::ConcurrentBooker::bookProfile(), doPixelLumi_, folderName_, ls_binning_, lumi_binning_, pixelCluster_binning_, pixellumi_binning_, pu_binning_, and DQMStore::IBooker::setCurrentFolder().

 {
   booker.setCurrentFolder(folderName_);
 
   if (doPixelLumi_) {
     histograms.numberOfPixelClustersVsLS = booker.book1D(
         "numberOfPixelClustersVsLS",
         "number of pixel clusters vs LS",
         ls_binning_.nbins, ls_binning_.xmin, ls_binning_.xmax);
     histograms.numberOfPixelClustersVsLS.setAxisTitle("LS", 1);
     histograms.numberOfPixelClustersVsLS.setAxisTitle("number of pixel clusters", 2);
 
     histograms.numberOfPixelClustersVsLumi = booker.bookProfile(
         "numberOfPixelClustersVsLumi",
         "number of pixel clusters vs scal lumi",
         lumi_binning_.nbins, lumi_binning_.xmin, lumi_binning_.xmax,
         pixelCluster_binning_.xmin, pixelCluster_binning_.xmax);
     histograms.numberOfPixelClustersVsLumi.setAxisTitle("scal inst lumi E30 [Hz cm^{-2}]", 1);
     histograms.numberOfPixelClustersVsLumi.setAxisTitle("number of pixel clusters", 2);
 
     histograms.pixelLumiVsLS = booker.bookProfile(
         "pixelLumiVsLS",
         "pixel-lumi vs LS",
         ls_binning_.nbins, ls_binning_.xmin, ls_binning_.xmax,
         pixellumi_binning_.xmin, pixellumi_binning_.xmax);
     histograms.pixelLumiVsLS.setAxisTitle("LS", 1);
     histograms.pixelLumiVsLS.setAxisTitle("pixel-based inst lumi E30 [Hz cm^{-2}]", 2);
 
     histograms.pixelLumiVsLumi = booker.bookProfile(
         "pixelLumiVsLumi",
         "pixel-lumi vs scal lumi",
         lumi_binning_.nbins, lumi_binning_.xmin, lumi_binning_.xmax,
         pixellumi_binning_.xmin, lumi_binning_.xmax);
     histograms.pixelLumiVsLumi.setAxisTitle("scal inst lumi E30 [Hz cm^{-2}]", 1);
     histograms.pixelLumiVsLumi.setAxisTitle("pixel-based inst lumi E30 [Hz cm^{-2}]", 2);
   }
 
   histograms.lumiVsLS = booker.bookProfile(
       "lumiVsLS",
       "scal lumi vs LS",
       ls_binning_.nbins, ls_binning_.xmin, ls_binning_.xmax,
       lumi_binning_.xmin, lumi_binning_.xmax);
   histograms.lumiVsLS.setAxisTitle("LS", 1);
   histograms.lumiVsLS.setAxisTitle("scal inst lumi E30 [Hz cm^{-2}]", 2);
 
   histograms.puVsLS = booker.bookProfile(
       "puVsLS",
       "scal PU vs LS",
       ls_binning_.nbins, ls_binning_.xmin, ls_binning_.xmax,
       pu_binning_.xmin, pu_binning_.xmax);
   histograms.puVsLS.setAxisTitle("LS", 1);
   histograms.puVsLS.setAxisTitle("scal PU", 2);
 
 }

void LumiMonitor::dqmAnalyze	(	edm::Event const &	event,
		edm::EventSetup const &	setup,
		Histograms const &	histograms
	)		const

overrideprivate

Definition at line 182 of file LumiMonitor.cc.

References GetLumi::CM2_TO_NANOBARN, doPixelLumi_, edm::HandleBase::isValid(), eostools::ls(), lumi_factor_per_bx_, lumiScalersToken_, minNumberOfPixelsPerCluster_, minPixelClusterCharge_, PixelSubdetector::PixelBarrel, HiRegitMuonDetachedTripletStep_cff::pixelClusters, pixelClustersToken_, GeneralSetup::setup(), and DetId::subdetId().

 {
   int ls = event.id().luminosityBlock();
 
   float scal_lumi = -1.;
   float scal_pu   = -1.;
   edm::Handle<LumiScalersCollection> lumiScalers;
   event.getByToken(lumiScalersToken_, lumiScalers);
   if (lumiScalers.isValid() and not lumiScalers->empty()) {
     auto scalit = lumiScalers->begin();
     scal_lumi = scalit->instantLumi();
     scal_pu   = scalit->pileup();
   } else {
     scal_lumi = -1.;
     scal_pu   = -1.;
   }
   histograms.lumiVsLS.fill(ls, scal_lumi);
   histograms.puVsLS.fill(ls, scal_pu);
 
   if (doPixelLumi_) {
     size_t pixel_clusters = 0;
     float  pixel_lumi = -1.;
     edm::Handle<edmNew::DetSetVector<SiPixelCluster>> pixelClusters;
     event.getByToken(pixelClustersToken_, pixelClusters);
     if (pixelClusters.isValid()) {
       auto const& tTopo = edm::get<TrackerTopology, TrackerTopologyRcd>(setup);
 
       // Count the number of clusters with at least a minimum
       // number of pixels per cluster and at least a minimum charge.
       size_t tot = 0;
       for (auto pixCluDet = pixelClusters->begin(); pixCluDet != pixelClusters->end(); ++pixCluDet) {
         
         DetId detid = pixCluDet->detId();
         size_t subdetid = detid.subdetId();
         if (subdetid == (int) PixelSubdetector::PixelBarrel) {
           if (tTopo.layer(detid)==1) {
             continue;
           }
         }
         
         for (auto pixClu = pixCluDet->begin(); pixClu != pixCluDet->end(); ++pixClu) {
           ++tot;
           if ((pixClu->size()   >= minNumberOfPixelsPerCluster_) and
               (pixClu->charge() >= minPixelClusterCharge_      )) {
             ++pixel_clusters;
           }
         }
       }
       pixel_lumi = lumi_factor_per_bx_ * pixel_clusters / GetLumi::CM2_TO_NANOBARN ; // ?!?!
     } else {
       pixel_lumi = -1.;
     }
 
     histograms.numberOfPixelClustersVsLS.fill(ls, pixel_clusters);
     histograms.numberOfPixelClustersVsLumi.fill(scal_lumi, pixel_clusters);
     histograms.pixelLumiVsLS.fill(ls, pixel_lumi);
     histograms.pixelLumiVsLumi.fill(scal_lumi, pixel_lumi);
   }
 
 }

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

static

Definition at line 255 of file LumiMonitor.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), DEFINE_FWK_MODULE, funct::false, fillHistoLSPSetDescription(), fillHistoPSetDescription(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>( "pixelClusters",               edm::InputTag("hltSiPixelClusters") );
   desc.add<edm::InputTag>( "scalers",                     edm::InputTag("hltScalersRawToDigi"));
   desc.add<std::string>  ( "FolderName",                  "HLT/LumiMonitoring"                );
   desc.add<bool>         ( "doPixelLumi",                 false                               );
   desc.add<bool>         ( "useBPixLayer1",               false                               );
   desc.add<int>          ( "minNumberOfPixelsPerCluster", 2                                   ); // from DQM/PixelLumi/python/PixelLumiDQM_cfi.py
   desc.add<double>       ( "minPixelClusterCharge",       15000.                              );
 
   edm::ParameterSetDescription histoPSet;
   edm::ParameterSetDescription pixelClusterPSet;
   LumiMonitor::fillHistoPSetDescription(pixelClusterPSet);
   histoPSet.add("pixelClusterPSet", pixelClusterPSet);
 
   edm::ParameterSetDescription lumiPSet;
   fillHistoPSetDescription(lumiPSet);
   histoPSet.add<edm::ParameterSetDescription>("lumiPSet", lumiPSet);
 
   edm::ParameterSetDescription puPSet;
   fillHistoPSetDescription(puPSet);
   histoPSet.add<edm::ParameterSetDescription>("puPSet", puPSet);
 
   edm::ParameterSetDescription pixellumiPSet;
   fillHistoPSetDescription(pixellumiPSet);
   histoPSet.add<edm::ParameterSetDescription>("pixellumiPSet", pixellumiPSet);
 
   edm::ParameterSetDescription lsPSet;
   fillHistoLSPSetDescription(lsPSet);
   histoPSet.add<edm::ParameterSetDescription>("lsPSet", lsPSet);
 
   desc.add<edm::ParameterSetDescription>("histoPSet", histoPSet);
 
   descriptions.add("lumiMonitor", desc);
 }

void LumiMonitor::fillHistoLSPSetDescription ( edm::ParameterSetDescription & pset )

static

Definition at line 250 of file LumiMonitor.cc.

References edm::ParameterSetDescription::add().

Referenced by fillDescriptions().

 {
   pset.add<int>   ( "nbins", 2500);
 }

void LumiMonitor::fillHistoPSetDescription ( edm::ParameterSetDescription & pset )

static

Definition at line 243 of file LumiMonitor.cc.

References edm::ParameterSetDescription::add().

Referenced by fillDescriptions().

 {
   pset.add<int>   ( "nbins");
   pset.add<double>( "xmin" );
   pset.add<double>( "xmax" );
 }

MEbinning LumiMonitor::getHistoLSPSet ( const edm::ParameterSet & pset )

staticprivate

Definition at line 118 of file LumiMonitor.cc.

References edm::ParameterSet::getParameter().

 {
   return MEbinning{
     pset.getParameter<int32_t>("nbins"),
     -0.5,
     pset.getParameter<int32_t>("nbins") - 0.5
   };
 }

MEbinning LumiMonitor::getHistoPSet ( const edm::ParameterSet & pset )

staticprivate

Definition at line 109 of file LumiMonitor.cc.

References edm::ParameterSet::getParameter().

Referenced by LumiMonitor().

 {
   return MEbinning{
     pset.getParameter<int32_t>("nbins"),
     pset.getParameter<double>("xmin"),
     pset.getParameter<double>("xmax"),
   };
 }