BeamMonitor Class Reference

#include <BeamMonitor.h>

Inheritance diagram for BeamMonitor:

Public Member Functions
	BeamMonitor (const edm::ParameterSet &)
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
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
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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > 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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override
void	beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &context) override
void	bookHistograms (DQMStore::IBooker &i, const edm::Run &r, const edm::EventSetup &c) override
void	dqmEndRun (const edm::Run &r, const edm::EventSetup &c) override
void	endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c) override
Protected Member Functions inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
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)
template<Transition B>
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 (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Private Member Functions
void	FitAndFill (const edm::LuminosityBlock &lumiSeg, int &, int &, int &)
void	formatFitTime (char *, const std::time_t &)
void	RestartFitting ()
void	scrollTH1 (TH1 *, std::time_t)
bool	testScroll (std::time_t &, std::time_t &)

Private Attributes
int	beginLumiOfBSFit_
int	beginLumiOfPVFit_
edm::EDGetTokenT< reco::BeamSpot >	bsSrc_
int	countEvt_
int	countGapLumi_
int	countLumi_
MonitorElement *	cutFlowTable
int	DBloggerReturn_
bool	debug_
double	deltaSigCut_
const int	dxBin_
const double	dxMax_
const double	dxMin_
const int	dzBin_
const double	dzMax_
const double	dzMin_
int	endLumiOfBSFit_
int	endLumiOfPVFit_
int	firstAverageFit_
int	fitNLumi_
int	fitPVNLumi_
MonitorElement *	fitResults
int	frun
edm::TimeValue_t	ftimestamp
MonitorElement *	h_d0_phi0
MonitorElement *	h_nTrk_lumi
MonitorElement *	h_nVtx
MonitorElement *	h_nVtx_lumi
MonitorElement *	h_nVtx_lumi_all
MonitorElement *	h_nVtx_st
MonitorElement *	h_PVx [2]
MonitorElement *	h_PVxz
MonitorElement *	h_PVy [2]
MonitorElement *	h_PVyz
MonitorElement *	h_PVz [2]
MonitorElement *	h_sigmaX0
MonitorElement *	h_sigmaY0
MonitorElement *	h_sigmaZ0
MonitorElement *	h_trk_z0
MonitorElement *	h_trkPt
MonitorElement *	h_trkVz
MonitorElement *	h_vx_dz
MonitorElement *	h_vx_vy
MonitorElement *	h_vy_dz
MonitorElement *	h_x0
MonitorElement *	h_y0
MonitorElement *	h_z0
edm::EDGetTokenT< edm::TriggerResults >	hltSrc_
std::vector< MonitorElement * >	hs
int	intervalInSec_
std::vector< std::string >	jetTrigger_
int	lastlumi_
int	lastNZbin
std::map< int, int >	mapBeginBSLS
std::map< int, std::time_t >	mapBeginBSTime
std::map< int, int >	mapBeginPVLS
std::map< int, std::time_t >	mapBeginPVTime
std::map< int, std::size_t >	mapLSBSTrkSize
std::map< int, TH1F >	mapLSCF
std::map< int, size_t >	mapLSPVStoreSize
std::map< int, std::vector< int > >	mapNPV
std::map< int, std::vector< float > >	mapPVx
std::map< int, std::vector< float > >	mapPVy
std::map< int, std::vector< float > >	mapPVz
double	maxZ_
unsigned int	min_Ntrks_
unsigned int	minNrVertices_
double	minVtxNdf_
double	minVtxWgt_
std::string	monitorName_
int	nextlumi_
int	nFitElements_
int	nFits_
unsigned int	nthBSTrk_
edm::Service< cond::service::OnlineDBOutputService >	onlineDbService_
bool	onlineMode_
const int	phiBin_
const double	phiMax_
const double	phiMin_
reco::BeamSpot	preBS
bool	processed_
MonitorElement *	pvResults
edm::EDGetTokenT< reco::VertexCollection >	pvSrc_
std::string	recordName_
reco::BeamSpot	refBS
std::time_t	refBStime [2]
std::time_t	refPVtime [2]
std::time_t	refTime
MonitorElement *	reportSummary
Float_t	reportSummary_
MonitorElement *	reportSummaryContents [3]
MonitorElement *	reportSummaryMap
int	resetFitNLumi_
bool	resetHistos_
int	resetPVNLumi_
bool	StartAverage_
std::time_t	startTime
Float_t	summaryContent_ [3]
Float_t	summarySum_
std::unique_ptr< BeamFitter >	theBeamFitter
std::time_t	tmpTime
edm::EDGetTokenT< reco::TrackCollection >	tracksLabel_
const int	vxBin_
const double	vxMax_
const double	vxMin_

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

Detailed Description

• Author

  Geng-yuan Jeng/UC Riverside Francisco Yumiceva/FNAL

Definition at line 34 of file BeamMonitor.h.

Constructor & Destructor Documentation

BeamMonitor()

BeamMonitor::BeamMonitor

(

const edm::ParameterSet &

ps

)

Definition at line 94 of file BeamMonitor.cc.

    : dxBin_(ps.getParameter<int>("dxBin")),
      dxMin_(ps.getParameter<double>("dxMin")),
      dxMax_(ps.getParameter<double>("dxMax")),
 
      vxBin_(ps.getParameter<int>("vxBin")),
      vxMin_(ps.getParameter<double>("vxMin")),
      vxMax_(ps.getParameter<double>("vxMax")),
 
      phiBin_(ps.getParameter<int>("phiBin")),
      phiMin_(ps.getParameter<double>("phiMin")),
      phiMax_(ps.getParameter<double>("phiMax")),
 
      dzBin_(ps.getParameter<int>("dzBin")),
      dzMin_(ps.getParameter<double>("dzMin")),
      dzMax_(ps.getParameter<double>("dzMax")),
 
      countEvt_(0),
      countLumi_(0),
      nthBSTrk_(0),
      nFitElements_(3),
      resetHistos_(false),
      StartAverage_(false),
      firstAverageFit_(0),
      countGapLumi_(0) {
  monitorName_ = ps.getUntrackedParameter<string>("monitorName", "YourSubsystemName");
  recordName_ = ps.getUntrackedParameter<string>("recordName");
  bsSrc_ = consumes<reco::BeamSpot>(ps.getUntrackedParameter<InputTag>("beamSpot"));
  tracksLabel_ = consumes<reco::TrackCollection>(
      ps.getParameter<ParameterSet>("BeamFitter").getUntrackedParameter<InputTag>("TrackCollection"));
  pvSrc_ = consumes<reco::VertexCollection>(ps.getUntrackedParameter<InputTag>("primaryVertex"));
  hltSrc_ = consumes<TriggerResults>(ps.getParameter<InputTag>("hltResults"));
  intervalInSec_ = ps.getUntrackedParameter<int>("timeInterval", 920);  //40 LS X 23"
  fitNLumi_ = ps.getUntrackedParameter<int>("fitEveryNLumi", -1);
  resetFitNLumi_ = ps.getUntrackedParameter<int>("resetEveryNLumi", -1);
  fitPVNLumi_ = ps.getUntrackedParameter<int>("fitPVEveryNLumi", -1);
  resetPVNLumi_ = ps.getUntrackedParameter<int>("resetPVEveryNLumi", -1);
  deltaSigCut_ = ps.getUntrackedParameter<double>("deltaSignificanceCut", 15);
  debug_ = ps.getUntrackedParameter<bool>("Debug");
  onlineMode_ = ps.getUntrackedParameter<bool>("OnlineMode");
  jetTrigger_ = ps.getUntrackedParameter<std::vector<std::string> >("jetTrigger");
  min_Ntrks_ = ps.getParameter<ParameterSet>("BeamFitter").getUntrackedParameter<int>("MinimumInputTracks");
  maxZ_ = ps.getParameter<ParameterSet>("BeamFitter").getUntrackedParameter<double>("MaximumZ");
  minNrVertices_ = ps.getParameter<ParameterSet>("PVFitter").getUntrackedParameter<unsigned int>("minNrVerticesForFit");
  minVtxNdf_ = ps.getParameter<ParameterSet>("PVFitter").getUntrackedParameter<double>("minVertexNdf");
  minVtxWgt_ = ps.getParameter<ParameterSet>("PVFitter").getUntrackedParameter<double>("minVertexMeanWeight");
 
  if (!monitorName_.empty())
    monitorName_ = monitorName_ + "/";
 
  theBeamFitter = std::make_unique<BeamFitter>(ps, consumesCollector());
  theBeamFitter->resetTrkVector();
  theBeamFitter->resetLSRange();
  theBeamFitter->resetRefTime();
  theBeamFitter->resetPVFitter();
 
  if (fitNLumi_ <= 0)
    fitNLumi_ = 1;
  nFits_ = beginLumiOfBSFit_ = endLumiOfBSFit_ = beginLumiOfPVFit_ = endLumiOfPVFit_ = 0;
  refBStime[0] = refBStime[1] = refPVtime[0] = refPVtime[1] = 0;
  maxZ_ = std::fabs(maxZ_);
  lastlumi_ = 0;
  nextlumi_ = 0;
  processed_ = false;
}

References beginLumiOfBSFit_, beginLumiOfPVFit_, bsSrc_, edm::EDConsumerBase::consumesCollector(), debug_, deltaSigCut_, endLumiOfBSFit_, endLumiOfPVFit_, fitNLumi_, fitPVNLumi_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hltSrc_, intervalInSec_, jetTrigger_, lastlumi_, maxZ_, min_Ntrks_, minNrVertices_, minVtxNdf_, minVtxWgt_, monitorName_, nextlumi_, nFits_, onlineMode_, processed_, pvSrc_, recordName_, refBStime, refPVtime, resetFitNLumi_, resetPVNLumi_, theBeamFitter, and tracksLabel_.

Member Function Documentation

analyze()

void BeamMonitor::analyze ( const edm::Event &  e, const edm::EventSetup &  c  )

overrideprotectedvirtual

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

Definition at line 646 of file BeamMonitor.cc.

                                                                       {
  const int nthlumi = iEvent.luminosityBlock();
  if (onlineMode_ && (nthlumi < nextlumi_)) {
    edm::LogInfo("BeamMonitor") << "analyze::  Spilt event from previous lumi section!" << std::endl;
    return;
  }
  if (onlineMode_ && (nthlumi > nextlumi_)) {
    edm::LogInfo("BeamMonitor") << "analyze::  Spilt event from next lumi section!!!" << std::endl;
    return;
  }
 
  countEvt_++;
  theBeamFitter->readEvent(
      iEvent);  //Remember when track fitter read the event in the same place the PVFitter read the events !!!!!!!!!
 
  Handle<reco::BeamSpot> recoBeamSpotHandle;
  iEvent.getByToken(bsSrc_, recoBeamSpotHandle);
  refBS = *recoBeamSpotHandle;
 
  //------Cut Flow Table filled every event!--------------------------------------
  {
    // Make a copy of the cut flow table from the beam fitter.
    auto tmphisto = static_cast<TH1F*>(theBeamFitter->getCutFlow());
    cutFlowTable->getTH1()->SetBins(
        tmphisto->GetNbinsX(), tmphisto->GetXaxis()->GetXmin(), tmphisto->GetXaxis()->GetXmax());
    // Update the bin labels
    if (countEvt_ == 1)  // SetLabel just once
      for (int n = 0; n < tmphisto->GetNbinsX(); n++)
        cutFlowTable->setBinLabel(n + 1, tmphisto->GetXaxis()->GetBinLabel(n + 1), 1);
    cutFlowTable->Reset();
    cutFlowTable->getTH1()->Add(tmphisto);
  }
 
  //----Reco tracks -------------------------------------
  Handle<reco::TrackCollection> TrackCollection;
  iEvent.getByToken(tracksLabel_, TrackCollection);
  const reco::TrackCollection* tracks = TrackCollection.product();
  for (reco::TrackCollection::const_iterator track = tracks->begin(); track != tracks->end(); ++track) {
    h_trkPt->Fill(track->pt());  //no need to change  here for average bs
    h_trkVz->Fill(track->vz());
  }
 
  //-------HLT Trigger --------------------------------
  edm::Handle<TriggerResults> triggerResults;
  bool JetTrigPass = false;
  if (iEvent.getByToken(hltSrc_, triggerResults)) {
    const edm::TriggerNames& trigNames = iEvent.triggerNames(*triggerResults);
    for (unsigned int i = 0; i < triggerResults->size(); i++) {
      const std::string& trigName = trigNames.triggerName(i);
 
      if (JetTrigPass)
        continue;
 
      for (size_t t = 0; t < jetTrigger_.size(); ++t) {
        if (JetTrigPass)
          continue;
 
        string string_search(jetTrigger_[t]);
        size_t found = trigName.find(string_search);
 
        if (found != string::npos) {
          int thisTrigger_ = trigNames.triggerIndex(trigName);
          if (triggerResults->accept(thisTrigger_))
            JetTrigPass = true;
        }  //if trigger found
      }    //for(t=0;..)
    }      //for(i=0; ..)
  }        //if trigger colleciton exist)
 
  //------ Primary Vertices-------
  edm::Handle<reco::VertexCollection> PVCollection;
 
  if (iEvent.getByToken(pvSrc_, PVCollection)) {
    int nPVcount = 0;
    int nPVcount_ST = 0;  //For Single Trigger(hence ST)
 
    for (reco::VertexCollection::const_iterator pv = PVCollection->begin(); pv != PVCollection->end(); ++pv) {
      //--- vertex selection
      if (pv->isFake() || pv->tracksSize() == 0)
        continue;
      nPVcount++;  // count non fake pv:
 
      if (JetTrigPass)
        nPVcount_ST++;  //non-fake pv with a specific trigger
 
      if (pv->ndof() < minVtxNdf_ || (pv->ndof() + 3.) / pv->tracksSize() < 2 * minVtxWgt_)
        continue;
 
      //Fill this map to store xyx for pv so that later we can remove the first one for run aver
      mapPVx[countLumi_].push_back(pv->x());
      mapPVy[countLumi_].push_back(pv->y());
      mapPVz[countLumi_].push_back(pv->z());
 
      if (!StartAverage_) {  //for first N LS
        h_PVx[0]->Fill(pv->x());
        h_PVy[0]->Fill(pv->y());
        h_PVz[0]->Fill(pv->z());
        h_PVxz->Fill(pv->z(), pv->x());
        h_PVyz->Fill(pv->z(), pv->y());
      }  //for first N LiS
      else {
        h_PVxz->Fill(pv->z(), pv->x());
        h_PVyz->Fill(pv->z(), pv->y());
      }
 
    }  //loop over pvs
 
    h_nVtx->Fill(nPVcount * 1.);  //no need to change it for average BS
 
    mapNPV[countLumi_].push_back((nPVcount_ST));
 
    if (!StartAverage_) {
      h_nVtx_st->Fill(nPVcount_ST * 1.);
    }
 
  }  //if pv collection is availaable
 
  if (StartAverage_) {
    map<int, std::vector<float> >::iterator itpvx = mapPVx.begin();
    map<int, std::vector<float> >::iterator itpvy = mapPVy.begin();
    map<int, std::vector<float> >::iterator itpvz = mapPVz.begin();
 
    map<int, std::vector<int> >::iterator itbspvinfo = mapNPV.begin();
 
    if ((int)mapPVx.size() > resetFitNLumi_) {  //sometimes the events is not there but LS is there!
      mapPVx.erase(itpvx);
      mapPVy.erase(itpvy);
      mapPVz.erase(itpvz);
      mapNPV.erase(itbspvinfo);
    }  //loop over Last N lumi collected
 
  }  //StartAverage==true
 
  processed_ = true;
}

References bsSrc_, countEvt_, countLumi_, cutFlowTable, dqm::impl::MonitorElement::Fill(), newFWLiteAna::found, dqm::impl::MonitorElement::getTH1(), h_nVtx, h_nVtx_st, h_PVx, h_PVxz, h_PVy, h_PVyz, h_PVz, h_trkPt, h_trkVz, hltSrc_, mps_fire::i, iEvent, jetTrigger_, mapNPV, mapPVx, mapPVy, mapPVz, minVtxNdf_, minVtxWgt_, dqmiodumpmetadata::n, nextlumi_, onlineMode_, processed_, MetAnalyzer::pv(), L1TEGammaOffline_cfi::PVCollection, pvSrc_, refBS, dqm::impl::MonitorElement::Reset(), resetFitNLumi_, dqm::impl::MonitorElement::setBinLabel(), StartAverage_, AlCaHLTBitMon_QueryRunRegistry::string, submitPVValidationJobs::t, theBeamFitter, HLT_FULL_cff::track, tracks, tracksLabel_, triggerResults, EgHLTOffTrigSelection_cfi::trigName, and trigNames.

beginLuminosityBlock()

void BeamMonitor::beginLuminosityBlock ( const edm::LuminosityBlock &  lumiSeg, const edm::EventSetup &  context  )

overrideprotected

Definition at line 517 of file BeamMonitor.cc.

                                                                                                {
  // start DB logger
  DBloggerReturn_ = 0;
  if (onlineDbService_.isAvailable()) {
    onlineDbService_->logger().start();
    onlineDbService_->logger().logInfo() << "BeamMonitor::beginLuminosityBlock - LS: " << lumiSeg.luminosityBlock();
  }
 
  int nthlumi = lumiSeg.luminosityBlock();
  const edm::TimeValue_t fbegintimestamp = lumiSeg.beginTime().value();
  const std::time_t ftmptime = fbegintimestamp >> 32;
 
  if (countLumi_ == 0 && (!processed_)) {
    beginLumiOfBSFit_ = beginLumiOfPVFit_ = nthlumi;
    refBStime[0] = refPVtime[0] = ftmptime;
    mapBeginBSLS[countLumi_] = nthlumi;
    mapBeginPVLS[countLumi_] = nthlumi;
    mapBeginBSTime[countLumi_] = ftmptime;
    mapBeginPVTime[countLumi_] = ftmptime;
  }  //for the first record
 
  if (nthlumi > nextlumi_) {
    if (processed_) {
      countLumi_++;
      //store here them will need when we remove the first one of Last N LS
      mapBeginBSLS[countLumi_] = nthlumi;
      mapBeginPVLS[countLumi_] = nthlumi;
      mapBeginBSTime[countLumi_] = ftmptime;
      mapBeginPVTime[countLumi_] = ftmptime;
    }  //processed passed but not the first lumi
    if ((!processed_) && countLumi_ != 0) {
      mapBeginBSLS[countLumi_] = nthlumi;
      mapBeginPVLS[countLumi_] = nthlumi;
      mapBeginBSTime[countLumi_] = ftmptime;
      mapBeginPVTime[countLumi_] = ftmptime;
    }  //processed fails for last lumi
  }    //nthLumi > nextlumi
 
  if (StartAverage_) {
    //Just Make sure it get rest
    refBStime[0] = 0;
    refPVtime[0] = 0;
    beginLumiOfPVFit_ = 0;
    beginLumiOfBSFit_ = 0;
 
    if (debug_)
      edm::LogInfo("BeamMonitor") << " beginLuminosityBlock:  Size of mapBeginBSLS before =  " << mapBeginBSLS.size()
                                  << endl;
    if (nthlumi >
        nextlumi_) {  //this make sure that it does not take into account this lumi for fitting and only look forward for new lumi
      //as countLumi also remains the same so map value  get overwritten once return to normal running.
      //even if few LS are misssing and DQM module do not sees them then it catchs up again
      map<int, int>::iterator itbs = mapBeginBSLS.begin();
      map<int, int>::iterator itpv = mapBeginPVLS.begin();
      map<int, std::time_t>::iterator itbstime = mapBeginBSTime.begin();
      map<int, std::time_t>::iterator itpvtime = mapBeginPVTime.begin();
 
      if (processed_) {  // otherwise if false then LS range of fit get messed up because we don't remove trk/pvs but we remove LS begin value . This prevent it as it happened if LS is there but no event are processed for some reason
        mapBeginBSLS.erase(itbs);
        mapBeginPVLS.erase(itpv);
        mapBeginBSTime.erase(itbstime);
        mapBeginPVTime.erase(itpvtime);
      }
      /*//not sure if want this or not ??
            map<int, int>::iterator itgapb=mapBeginBSLS.begin();
            map<int, int>::iterator itgape=mapBeginBSLS.end(); itgape--;
            countGapLumi_ = ( (itgape->second) - (itgapb->second) );
            //if we see Gap more than then 2*resetNFitLumi !!!!!!!
            //for example if 10-15 is fitted and if 16-25 are missing then we next fit will be for range 11-26 but BS can change in between
            // so better start  as fresh  and reset everything like starting in the begining!
            if(countGapLumi_ >= 2*resetFitNLumi_){RestartFitting(); mapBeginBSLS[countLumi_]   = nthlumi;}
            */
    }
 
    if (debug_)
      edm::LogInfo("BeamMonitor") << " beginLuminosityBlock::  Size of mapBeginBSLS After = " << mapBeginBSLS.size()
                                  << endl;
 
    map<int, int>::iterator bbs = mapBeginBSLS.begin();
    map<int, int>::iterator bpv = mapBeginPVLS.begin();
    map<int, std::time_t>::iterator bbst = mapBeginBSTime.begin();
    map<int, std::time_t>::iterator bpvt = mapBeginPVTime.begin();
 
    if (beginLumiOfPVFit_ == 0)
      beginLumiOfPVFit_ = bpv->second;  //new begin time after removing the LS
    if (beginLumiOfBSFit_ == 0)
      beginLumiOfBSFit_ = bbs->second;
    if (refBStime[0] == 0)
      refBStime[0] = bbst->second;
    if (refPVtime[0] == 0)
      refPVtime[0] = bpvt->second;
 
  }  //same logic for average fit as above commented line
 
  map<int, std::time_t>::iterator nbbst = mapBeginBSTime.begin();
  map<int, std::time_t>::iterator nbpvt = mapBeginPVTime.begin();
 
  if (onlineMode_ && (nthlumi < nextlumi_))
    return;
 
  if (onlineMode_) {
    if (nthlumi > nextlumi_) {
      if (countLumi_ != 0 && processed_)
        FitAndFill(lumiSeg, lastlumi_, nextlumi_, nthlumi);
      nextlumi_ = nthlumi;
      edm::LogInfo("BeamMonitor") << "beginLuminosityBlock:: Next Lumi to Fit: " << nextlumi_ << endl;
      if ((StartAverage_) && refBStime[0] == 0)
        refBStime[0] = nbbst->second;
      if ((StartAverage_) && refPVtime[0] == 0)
        refPVtime[0] = nbpvt->second;
    }
  } else {
    if (processed_)
      FitAndFill(lumiSeg, lastlumi_, nextlumi_, nthlumi);
    nextlumi_ = nthlumi;
    edm::LogInfo("BeamMonitor") << " beginLuminosityBlock:: Next Lumi to Fit: " << nextlumi_ << endl;
    if ((StartAverage_) && refBStime[0] == 0)
      refBStime[0] = nbbst->second;
    if ((StartAverage_) && refPVtime[0] == 0)
      refPVtime[0] = nbpvt->second;
  }
 
  //countLumi_++;
  if (processed_)
    processed_ = false;
  edm::LogInfo("BeamMonitor") << " beginLuminosityBlock::  Begin of Lumi: " << nthlumi << endl;
}

References beginLumiOfBSFit_, beginLumiOfPVFit_, edm::LuminosityBlockBase::beginTime(), countLumi_, DBloggerReturn_, debug_, FitAndFill(), edm::Service< T >::isAvailable(), lastlumi_, cond::service::PoolDBOutputService::logger(), cond::persistency::Logger::logInfo(), edm::LuminosityBlockBase::luminosityBlock(), mapBeginBSLS, mapBeginBSTime, mapBeginPVLS, mapBeginPVTime, nextlumi_, onlineDbService_, onlineMode_, processed_, refBStime, refPVtime, cond::persistency::Logger::start(), StartAverage_, and edm::Timestamp::value().

bookHistograms()

void BeamMonitor::bookHistograms ( DQMStore::IBooker &  i, const edm::Run &  r, const edm::EventSetup &  c  )

overrideprotectedvirtual

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

Definition at line 206 of file BeamMonitor.cc.

                                                                                                          {
  frun = iRun.run();
  ftimestamp = iRun.beginTime().value();
  tmpTime = ftimestamp >> 32;
  startTime = refTime = tmpTime;
  char eventTime[64];
  formatFitTime(eventTime, tmpTime);
  edm::LogInfo("BeamMonitor") << "TimeOffset = " << eventTime << std::endl;
  TDatime da(eventTime);
  if (debug_) {
    edm::LogInfo("BeamMonitor") << "TimeOffset = ";
    da.Print();
  }
  auto daTime = da.Convert(kTRUE);
 
  // book some histograms here
 
  // create and cd into new folder
  iBooker.setCurrentFolder(monitorName_ + "Fit");
 
  h_nTrk_lumi = iBooker.book1D("nTrk_lumi", "Num. of selected tracks vs lumi (Fit)", 20, 0.5, 20.5);
  h_nTrk_lumi->setAxisTitle("Lumisection", 1);
  h_nTrk_lumi->setAxisTitle("Num of Tracks for Fit", 2);
 
  //store vtx vs lumi for monitoring why fits fail
  h_nVtx_lumi = iBooker.book1D("nVtx_lumi", "Num. of selected Vtx vs lumi (Fit)", 20, 0.5, 20.5);
  h_nVtx_lumi->setAxisTitle("Lumisection", 1);
  h_nVtx_lumi->setAxisTitle("Num of Vtx for Fit", 2);
 
  h_nVtx_lumi_all = iBooker.book1D("nVtx_lumi_all", "Num. of selected Vtx vs lumi (Fit) all", 20, 0.5, 20.5);
  h_nVtx_lumi_all->getTH1()->SetCanExtend(TH1::kAllAxes);
  h_nVtx_lumi_all->setAxisTitle("Lumisection", 1);
  h_nVtx_lumi_all->setAxisTitle("Num of Vtx for Fit", 2);
 
  h_d0_phi0 = iBooker.bookProfile(
      "d0_phi0", "d_{0} vs. #phi_{0} (Selected Tracks)", phiBin_, phiMin_, phiMax_, dxBin_, dxMin_, dxMax_, "");
  h_d0_phi0->setAxisTitle("#phi_{0} (rad)", 1);
  h_d0_phi0->setAxisTitle("d_{0} (cm)", 2);
 
  h_vx_vy = iBooker.book2D(
      "trk_vx_vy", "Vertex (PCA) position of selected tracks", vxBin_, vxMin_, vxMax_, vxBin_, vxMin_, vxMax_);
  h_vx_vy->setOption("COLZ");
  //   h_vx_vy->getTH1()->SetCanExtend(TH1::kAllAxes);
  h_vx_vy->setAxisTitle("x coordinate of input track at PCA (cm)", 1);
  h_vx_vy->setAxisTitle("y coordinate of input track at PCA (cm)", 2);
 
  {
    TDatime* da = new TDatime();
    gStyle->SetTimeOffset(da->Convert(kTRUE));
  }
 
  const int nvar_ = 6;
  string coord[nvar_] = {"x", "y", "z", "sigmaX", "sigmaY", "sigmaZ"};
  string label[nvar_] = {
      "x_{0} (cm)", "y_{0} (cm)", "z_{0} (cm)", "#sigma_{X_{0}} (cm)", "#sigma_{Y_{0}} (cm)", "#sigma_{Z_{0}} (cm)"};
 
  hs.reserve(kNumHists);
  for (int i = 0; i < 4; i++) {
    iBooker.setCurrentFolder(monitorName_ + "Fit");
    for (int ic = 0; ic < nvar_; ++ic) {
      TString histName(coord[ic]);
      TString histTitle(coord[ic]);
      string ytitle(label[ic]);
      string xtitle("");
      string options("E1");
      bool createHisto = true;
      switch (i) {
        case 1:  // BS vs time
          histName += "0_time";
          xtitle = "Time [UTC]";
          if (ic < 3)
            histTitle += " coordinate of beam spot vs time (Fit)";
          else
            histTitle = histTitle.Insert(5, " ") + " of beam spot vs time (Fit)";
          break;
        case 2:  // PV vs lumi
          if (ic < 3) {
            iBooker.setCurrentFolder(monitorName_ + "PrimaryVertex");
            histName.Insert(0, "PV");
            histName += "_lumi";
            histTitle.Insert(0, "Avg. ");
            histTitle += " position of primary vtx vs lumi";
            xtitle = "Lumisection";
            ytitle.insert(0, "PV");
            ytitle += " #pm #sigma_{PV";
            ytitle += coord[ic];
            ytitle += "} (cm)";
          } else
            createHisto = false;
          break;
        case 3:  // PV vs time
          if (ic < 3) {
            iBooker.setCurrentFolder(monitorName_ + "PrimaryVertex");
            histName.Insert(0, "PV");
            histName += "_time";
            histTitle.Insert(0, "Avg. ");
            histTitle += " position of primary vtx vs time";
            xtitle = "Time [UTC]";
            ytitle.insert(0, "PV");
            ytitle += " #pm #sigma_{PV";
            ytitle += coord[ic];
            ytitle += "} (cm)";
          } else
            createHisto = false;
          break;
        default:  // BS vs lumi
          histName += "0_lumi";
          xtitle = "Lumisection";
          if (ic < 3)
            histTitle += " coordinate of beam spot vs lumi (Fit)";
          else
            histTitle = histTitle.Insert(5, " ") + " of beam spot vs lumi (Fit)";
          break;
      }
      if (createHisto) {
        edm::LogInfo("BeamMonitor") << "hitsName = " << histName << "; histTitle = " << histTitle << std::endl;
        auto tmpHs = iBooker.book1D(histName, histTitle, 40, 0.5, 40.5);
        hs.push_back(tmpHs);
        tmpHs->setAxisTitle(xtitle, 1);
        tmpHs->setAxisTitle(ytitle, 2);
        tmpHs->getTH1()->SetOption("E1");
        if (histName.Contains("time")) {
          //int nbins = (intervalInSec_/23 > 0 ? intervalInSec_/23 : 40);
          tmpHs->getTH1()->SetBins(intervalInSec_, 0.5, intervalInSec_ + 0.5);
          tmpHs->setAxisTimeDisplay(1);
          tmpHs->setAxisTimeFormat("%H:%M:%S", 1);
          tmpHs->getTH1()->GetXaxis()->SetTimeOffset(daTime);
        }
        histName += "_all";
        histTitle += " all";
        tmpHs = iBooker.book1D(histName, histTitle, 40, 0.5, 40.5);
        hs.push_back(tmpHs);
        tmpHs->getTH1()->SetCanExtend(TH1::kAllAxes);
        tmpHs->setAxisTitle(xtitle, 1);
        tmpHs->setAxisTitle(ytitle, 2);
        tmpHs->getTH1()->SetOption("E1");
        if (histName.Contains("time")) {
          //int nbins = (intervalInSec_/23 > 0 ? intervalInSec_/23 : 40);
          tmpHs->getTH1()->SetBins(intervalInSec_, 0.5, intervalInSec_ + 0.5);
          tmpHs->setAxisTimeDisplay(1);
          tmpHs->setAxisTimeFormat("%H:%M:%S", 1);
          tmpHs->getTH1()->GetXaxis()->SetTimeOffset(daTime);
        }
      }
    }
  }
  assert(hs.size() == kNumHists);
  assert(0 == strcmp(hs[k_sigmaY0_time]->getTH1()->GetName(), "sigmaY0_time"));
  assert(0 == strcmp(hs[k_PVz_lumi_all]->getTH1()->GetName(), "PVz_lumi_all"));
 
  iBooker.setCurrentFolder(monitorName_ + "Fit");
 
  h_trk_z0 = iBooker.book1D("trk_z0", "z_{0} of selected tracks", dzBin_, dzMin_, dzMax_);
  h_trk_z0->setAxisTitle("z_{0} of selected tracks (cm)", 1);
 
  h_vx_dz = iBooker.bookProfile(
      "vx_dz", "v_{x} vs. dz of selected tracks", dzBin_, dzMin_, dzMax_, dxBin_, dxMin_, dxMax_, "");
  h_vx_dz->setAxisTitle("dz (cm)", 1);
  h_vx_dz->setAxisTitle("x coordinate of input track at PCA (cm)", 2);
 
  h_vy_dz = iBooker.bookProfile(
      "vy_dz", "v_{y} vs. dz of selected tracks", dzBin_, dzMin_, dzMax_, dxBin_, dxMin_, dxMax_, "");
  h_vy_dz->setAxisTitle("dz (cm)", 1);
  h_vy_dz->setAxisTitle("y coordinate of input track at PCA (cm)", 2);
 
  h_x0 = iBooker.book1D("BeamMonitorFeedBack_x0", "x coordinate of beam spot (Fit)", 100, -0.01, 0.01);
  h_x0->setAxisTitle("x_{0} (cm)", 1);
  h_x0->getTH1()->SetCanExtend(TH1::kAllAxes);
 
  h_y0 = iBooker.book1D("BeamMonitorFeedBack_y0", "y coordinate of beam spot (Fit)", 100, -0.01, 0.01);
  h_y0->setAxisTitle("y_{0} (cm)", 1);
  h_y0->getTH1()->SetCanExtend(TH1::kAllAxes);
 
  h_z0 = iBooker.book1D("BeamMonitorFeedBack_z0", "z coordinate of beam spot (Fit)", dzBin_, dzMin_, dzMax_);
  h_z0->setAxisTitle("z_{0} (cm)", 1);
  h_z0->getTH1()->SetCanExtend(TH1::kAllAxes);
 
  h_sigmaX0 = iBooker.book1D("BeamMonitorFeedBack_sigmaX0", "sigma x0 of beam spot (Fit)", 100, 0, 0.05);
  h_sigmaX0->setAxisTitle("#sigma_{X_{0}} (cm)", 1);
  h_sigmaX0->getTH1()->SetCanExtend(TH1::kAllAxes);
 
  h_sigmaY0 = iBooker.book1D("BeamMonitorFeedBack_sigmaY0", "sigma y0 of beam spot (Fit)", 100, 0, 0.05);
  h_sigmaY0->setAxisTitle("#sigma_{Y_{0}} (cm)", 1);
  h_sigmaY0->getTH1()->SetCanExtend(TH1::kAllAxes);
 
  h_sigmaZ0 = iBooker.book1D("BeamMonitorFeedBack_sigmaZ0", "sigma z0 of beam spot (Fit)", 100, 0, 10);
  h_sigmaZ0->setAxisTitle("#sigma_{Z_{0}} (cm)", 1);
  h_sigmaZ0->getTH1()->SetCanExtend(TH1::kAllAxes);
 
  // Histograms of all reco tracks (without cuts):
  h_trkPt = iBooker.book1D("trkPt", "p_{T} of all reco'd tracks (no selection)", 200, 0., 50.);
  h_trkPt->setAxisTitle("p_{T} (GeV/c)", 1);
 
  h_trkVz = iBooker.book1D("trkVz", "Z coordinate of PCA of all reco'd tracks (no selection)", dzBin_, dzMin_, dzMax_);
  h_trkVz->setAxisTitle("V_{Z} (cm)", 1);
 
  cutFlowTable = iBooker.book1D("cutFlowTable", "Cut flow table of track selection", 9, 0, 9);
 
  // Results of previous good fit:
  fitResults = iBooker.book2D("fitResults", "Results of previous good beam fit", 2, 0, 2, 8, 0, 8);
  fitResults->setAxisTitle("Fitted Beam Spot (cm)", 1);
  fitResults->setBinLabel(8, "x_{0}", 2);
  fitResults->setBinLabel(7, "y_{0}", 2);
  fitResults->setBinLabel(6, "z_{0}", 2);
  fitResults->setBinLabel(5, "#sigma_{Z}", 2);
  fitResults->setBinLabel(4, "#frac{dx}{dz} (rad)", 2);
  fitResults->setBinLabel(3, "#frac{dy}{dz} (rad)", 2);
  fitResults->setBinLabel(2, "#sigma_{X}", 2);
  fitResults->setBinLabel(1, "#sigma_{Y}", 2);
  fitResults->setBinLabel(1, "Mean", 1);
  fitResults->setBinLabel(2, "Stat. Error", 1);
  fitResults->getTH1()->SetOption("text");
 
  // Histos of PrimaryVertices:
  iBooker.setCurrentFolder(monitorName_ + "PrimaryVertex");
 
  h_nVtx = iBooker.book1D("vtxNbr", "Reconstructed Vertices(non-fake) in all Event", 60, -0.5, 59.5);
  h_nVtx->setAxisTitle("Num. of reco. vertices", 1);
 
  //For one Trigger only
  h_nVtx_st = iBooker.book1D("vtxNbr_SelectedTriggers", "Reconstructed Vertices(non-fake) in Events", 60, -0.5, 59.5);
  //h_nVtx_st->setAxisTitle("Num. of reco. vertices for Un-Prescaled Jet Trigger",1);
 
  // Monitor only the PV with highest sum pt of assoc. trks:
  h_PVx[0] = iBooker.book1D("PVX", "x coordinate of Primary Vtx", 50, -0.01, 0.01);
  h_PVx[0]->setAxisTitle("PVx (cm)", 1);
  h_PVx[0]->getTH1()->SetCanExtend(TH1::kAllAxes);
 
  h_PVy[0] = iBooker.book1D("PVY", "y coordinate of Primary Vtx", 50, -0.01, 0.01);
  h_PVy[0]->setAxisTitle("PVy (cm)", 1);
  h_PVy[0]->getTH1()->SetCanExtend(TH1::kAllAxes);
 
  h_PVz[0] = iBooker.book1D("PVZ", "z coordinate of Primary Vtx", dzBin_, dzMin_, dzMax_);
  h_PVz[0]->setAxisTitle("PVz (cm)", 1);
 
  h_PVx[1] = iBooker.book1D("PVXFit", "x coordinate of Primary Vtx (Last Fit)", 50, -0.01, 0.01);
  h_PVx[1]->setAxisTitle("PVx (cm)", 1);
  h_PVx[1]->getTH1()->SetCanExtend(TH1::kAllAxes);
 
  h_PVy[1] = iBooker.book1D("PVYFit", "y coordinate of Primary Vtx (Last Fit)", 50, -0.01, 0.01);
  h_PVy[1]->setAxisTitle("PVy (cm)", 1);
  h_PVy[1]->getTH1()->SetCanExtend(TH1::kAllAxes);
 
  h_PVz[1] = iBooker.book1D("PVZFit", "z coordinate of Primary Vtx (Last Fit)", dzBin_, dzMin_, dzMax_);
  h_PVz[1]->setAxisTitle("PVz (cm)", 1);
 
  h_PVxz = iBooker.bookProfile("PVxz", "PVx vs. PVz", dzBin_ / 2, dzMin_, dzMax_, dxBin_ / 2, dxMin_, dxMax_, "");
  h_PVxz->setAxisTitle("PVz (cm)", 1);
  h_PVxz->setAxisTitle("PVx (cm)", 2);
 
  h_PVyz = iBooker.bookProfile("PVyz", "PVy vs. PVz", dzBin_ / 2, dzMin_, dzMax_, dxBin_ / 2, dxMin_, dxMax_, "");
  h_PVyz->setAxisTitle("PVz (cm)", 1);
  h_PVyz->setAxisTitle("PVy (cm)", 2);
 
  // Results of previous good fit:
  pvResults = iBooker.book2D("pvResults", "Results of fitting Primary Vertices", 2, 0, 2, 6, 0, 6);
  pvResults->setAxisTitle("Fitted Primary Vertex (cm)", 1);
  pvResults->setBinLabel(6, "PVx", 2);
  pvResults->setBinLabel(5, "PVy", 2);
  pvResults->setBinLabel(4, "PVz", 2);
  pvResults->setBinLabel(3, "#sigma_{X}", 2);
  pvResults->setBinLabel(2, "#sigma_{Y}", 2);
  pvResults->setBinLabel(1, "#sigma_{Z}", 2);
  pvResults->setBinLabel(1, "Mean", 1);
  pvResults->setBinLabel(2, "Stat. Error", 1);
  pvResults->getTH1()->SetOption("text");
 
  // Summary plots:
  iBooker.setCurrentFolder(monitorName_ + "EventInfo");
 
  reportSummary = iBooker.bookFloat("reportSummary");
  if (reportSummary)
    reportSummary->Fill(std::numeric_limits<double>::quiet_NaN());
 
  char histo[20];
  iBooker.setCurrentFolder(monitorName_ + "EventInfo/reportSummaryContents");
  for (int n = 0; n < nFitElements_; n++) {
    switch (n) {
      case 0:
        sprintf(histo, "x0_status");
        break;
      case 1:
        sprintf(histo, "y0_status");
        break;
      case 2:
        sprintf(histo, "z0_status");
        break;
    }
    reportSummaryContents[n] = iBooker.bookFloat(histo);
  }
 
  for (int i = 0; i < nFitElements_; i++) {
    summaryContent_[i] = 0.;
    reportSummaryContents[i]->Fill(std::numeric_limits<double>::quiet_NaN());
  }
 
  iBooker.setCurrentFolder(monitorName_ + "EventInfo");
 
  reportSummaryMap = iBooker.book2D("reportSummaryMap", "Beam Spot Summary Map", 1, 0, 1, 3, 0, 3);
  reportSummaryMap->setAxisTitle("", 1);
  reportSummaryMap->setAxisTitle("Fitted Beam Spot", 2);
  reportSummaryMap->setBinLabel(1, " ", 1);
  reportSummaryMap->setBinLabel(1, "x_{0}", 2);
  reportSummaryMap->setBinLabel(2, "y_{0}", 2);
  reportSummaryMap->setBinLabel(3, "z_{0}", 2);
  for (int i = 0; i < nFitElements_; i++) {
    reportSummaryMap->setBinContent(1, i + 1, -1.);
  }
}

References cms::cuda::assert(), edm::RunBase::beginTime(), dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookFloat(), dqm::implementation::IBooker::bookProfile(), cutFlowTable, Calorimetry_cff::da, debug_, dxBin_, dxMax_, dxMin_, dzBin_, dzMax_, dzMin_, dqm::impl::MonitorElement::Fill(), fitResults, formatFitTime(), frun, ftimestamp, dqm::impl::MonitorElement::getTH1(), h_d0_phi0, h_nTrk_lumi, h_nVtx, h_nVtx_lumi, h_nVtx_lumi_all, h_nVtx_st, h_PVx, h_PVxz, h_PVy, h_PVyz, h_PVz, h_sigmaX0, h_sigmaY0, h_sigmaZ0, h_trk_z0, h_trkPt, h_trkVz, h_vx_dz, h_vx_vy, h_vy_dz, h_x0, h_y0, h_z0, timingPdfMaker::histo, B2GTnPMonitor_cfi::histTitle, hs, mps_fire::i, intervalInSec_, label, monitorName_, dqmiodumpmetadata::n, nFitElements_, AlcaSiPixelAliHarvester0T_cff::options, phiBin_, phiMax_, phiMin_, pvResults, refTime, reportSummary, reportSummaryContents, reportSummaryMap, edm::RunBase::run(), dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinContent(), dqm::impl::MonitorElement::setBinLabel(), dqm::implementation::NavigatorBase::setCurrentFolder(), dqm::impl::MonitorElement::setOption(), startTime, summaryContent_, tmpTime, edm::Timestamp::value(), vxBin_, vxMax_, vxMin_, hgcalPlots::xtitle, and hgcalPlots::ytitle.

dqmEndRun()

void BeamMonitor::dqmEndRun ( const edm::Run &  r, const edm::EventSetup &  c  )

overrideprotectedvirtual

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

Definition at line 1532 of file BeamMonitor.cc.

                                                                   {
  if (debug_)
    edm::LogInfo("BeamMonitor") << "dqmEndRun:: Clearing all the Maps " << endl;
  //Clear all the Maps here
  mapPVx.clear();
  mapPVy.clear();
  mapPVz.clear();
  mapNPV.clear();
  mapBeginBSLS.clear();
  mapBeginPVLS.clear();
  mapBeginBSTime.clear();
  mapBeginPVTime.clear();
  mapLSBSTrkSize.clear();
  mapLSPVStoreSize.clear();
  mapLSCF.clear();
}

References debug_, mapBeginBSLS, mapBeginBSTime, mapBeginPVLS, mapBeginPVTime, mapLSBSTrkSize, mapLSCF, mapLSPVStoreSize, mapNPV, mapPVx, mapPVy, and mapPVz.

endLuminosityBlock()

void BeamMonitor::endLuminosityBlock ( const edm::LuminosityBlock &  lumiSeg, const edm::EventSetup &  c  )

overrideprotected

Definition at line 783 of file BeamMonitor.cc.

                                                                                             {
  int nthlumi = lumiSeg.id().luminosityBlock();
  edm::LogInfo("BeamMonitor") << "endLuminosityBlock:: Lumi of the last event before endLuminosityBlock: " << nthlumi
                              << endl;
 
  if (onlineMode_ && nthlumi < nextlumi_)
    return;
  const edm::TimeValue_t fendtimestamp = lumiSeg.endTime().value();
  const std::time_t fendtime = fendtimestamp >> 32;
  tmpTime = refBStime[1] = refPVtime[1] = fendtime;
 
  // end DB logger
  if (onlineDbService_.isAvailable()) {
    onlineDbService_->logger().logInfo() << "BeamMonitor::endLuminosityBlock";
    onlineDbService_->logger().end(DBloggerReturn_);
  }
}

References DBloggerReturn_, cond::persistency::Logger::end(), edm::LuminosityBlockBase::endTime(), edm::LuminosityBlockBase::id(), edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::logger(), cond::persistency::Logger::logInfo(), edm::LuminosityBlockID::luminosityBlock(), nextlumi_, onlineDbService_, onlineMode_, refBStime, refPVtime, tmpTime, and edm::Timestamp::value().

FitAndFill()

void BeamMonitor::FitAndFill ( const edm::LuminosityBlock &  lumiSeg, int &  lastlumi, int &  nextlumi, int &  nthlumi  )

private

Definition at line 802 of file BeamMonitor.cc.

                                                                                                       {
  if (onlineMode_ && (nthlumi <= nextlumi))
    return;
 
  //set the correct run number when no event in the LS for fake output
  if ((processed_) && theBeamFitter->getRunNumber() != frun)
    theBeamFitter->setRun(frun);
 
  int currentlumi = nextlumi;
  edm::LogInfo("BeamMonitor") << "FitAndFill::  Lumi of the current fit: " << currentlumi << endl;
  lastlumi = currentlumi;
  endLumiOfBSFit_ = currentlumi;
  endLumiOfPVFit_ = currentlumi;
 
  //---------Fix for Runninv average-------------
  mapLSPVStoreSize[countLumi_] = theBeamFitter->getPVvectorSize();
 
  if (StartAverage_) {
    std::map<int, std::size_t>::iterator rmLSPVi = mapLSPVStoreSize.begin();
    size_t SizeToRemovePV = rmLSPVi->second;
    for (std::map<int, std::size_t>::iterator rmLSPVe = mapLSPVStoreSize.end(); ++rmLSPVi != rmLSPVe;)
      rmLSPVi->second -= SizeToRemovePV;
 
    theBeamFitter->resizePVvector(SizeToRemovePV);
 
    map<int, std::size_t>::iterator tmpItpv = mapLSPVStoreSize.begin();
    mapLSPVStoreSize.erase(tmpItpv);
  }
  if (debug_)
    edm::LogInfo("BeamMonitor") << "FitAndFill:: Size of thePVvector After removing the PVs = "
                                << theBeamFitter->getPVvectorSize() << endl;
 
  //lets filt the PV for GUI here: It was in analyzer in preivous versiton but moved here due to absence of event in some lumis, works OK
  bool resetHistoFlag_ = false;
  if ((int)mapPVx.size() >= resetFitNLumi_ && (StartAverage_)) {
    h_PVx[0]->Reset();
    h_PVy[0]->Reset();
    h_PVz[0]->Reset();
    h_nVtx_st->Reset();
    resetHistoFlag_ = true;
  }
 
  int MaxPVs = 0;
  int countEvtLastNLS_ = 0;
  int countTotPV_ = 0;
 
  std::map<int, std::vector<int> >::iterator mnpv = mapNPV.begin();
  std::map<int, std::vector<float> >::iterator mpv2 = mapPVy.begin();
  std::map<int, std::vector<float> >::iterator mpv3 = mapPVz.begin();
 
  for (std::map<int, std::vector<float> >::iterator mpv1 = mapPVx.begin(); mpv1 != mapPVx.end();
       ++mpv1, ++mpv2, ++mpv3, ++mnpv) {
    std::vector<float>::iterator mpvs2 = (mpv2->second).begin();
    std::vector<float>::iterator mpvs3 = (mpv3->second).begin();
    for (std::vector<float>::iterator mpvs1 = (mpv1->second).begin(); mpvs1 != (mpv1->second).end();
         ++mpvs1, ++mpvs2, ++mpvs3) {
      if (resetHistoFlag_) {
        h_PVx[0]->Fill(*mpvs1);  //these histogram are reset after StartAverage_ flag is ON
        h_PVy[0]->Fill(*mpvs2);
        h_PVz[0]->Fill(*mpvs3);
      }
    }  //loop over second
 
    //Do the same here for nPV distr.
    for (std::vector<int>::iterator mnpvs = (mnpv->second).begin(); mnpvs != (mnpv->second).end(); ++mnpvs) {
      if ((*mnpvs > 0) && (resetHistoFlag_))
        h_nVtx_st->Fill((*mnpvs) * (1.0));
      countEvtLastNLS_++;
      countTotPV_ += (*mnpvs);
      if ((*mnpvs) > MaxPVs)
        MaxPVs = (*mnpvs);
    }  //loop over second of mapNPV
 
  }  //loop over last N lumis
 
  char tmpTitlePV[100];
  sprintf(tmpTitlePV, "%s %i %s %i", "Num. of reco. vertices for LS: ", beginLumiOfPVFit_, " to ", endLumiOfPVFit_);
  h_nVtx_st->setAxisTitle(tmpTitlePV, 1);
 
  std::vector<float> DipPVInfo_;
  DipPVInfo_.clear();
 
  if (countTotPV_ != 0) {
    DipPVInfo_.push_back((float)countEvtLastNLS_);
    DipPVInfo_.push_back(h_nVtx_st->getMean());
    DipPVInfo_.push_back(h_nVtx_st->getMeanError());
    DipPVInfo_.push_back(h_nVtx_st->getRMS());
    DipPVInfo_.push_back(h_nVtx_st->getRMSError());
    DipPVInfo_.push_back((float)MaxPVs);
    DipPVInfo_.push_back((float)countTotPV_);
    MaxPVs = 0;
  } else {
    for (size_t i = 0; i < 7; i++) {
      if (i > 0) {
        DipPVInfo_.push_back(0.);
      } else {
        DipPVInfo_.push_back((float)countEvtLastNLS_);
      }
    }
  }
  theBeamFitter->SetPVInfo(DipPVInfo_);
  countEvtLastNLS_ = 0;
 
  if (onlineMode_) {  // filling LS gap
    // FIXME: need to add protection for the case if the gap is at the resetting LS!
    const int countLS_bs = hs[k_x0_lumi]->getTH1()->GetEntries();
    const int countLS_pv = hs[k_PVx_lumi]->getTH1()->GetEntries();
    edm::LogInfo("BeamMonitor") << "FitAndFill:: countLS_bs = " << countLS_bs << " ; countLS_pv = " << countLS_pv
                                << std::endl;
    int LSgap_bs = currentlumi / fitNLumi_ - countLS_bs;
    int LSgap_pv = currentlumi / fitPVNLumi_ - countLS_pv;
    if (currentlumi % fitNLumi_ == 0)
      LSgap_bs--;
    if (currentlumi % fitPVNLumi_ == 0)
      LSgap_pv--;
    edm::LogInfo("BeamMonitor") << "FitAndFill::  LSgap_bs = " << LSgap_bs << " ; LSgap_pv = " << LSgap_pv << std::endl;
    // filling previous fits if LS gap ever exists
    for (int ig = 0; ig < LSgap_bs; ig++) {
      hs[k_x0_lumi]->ShiftFillLast(0., 0., fitNLumi_);  //x0 , x0err, fitNLumi_;  see DQMCore....
      hs[k_y0_lumi]->ShiftFillLast(0., 0., fitNLumi_);
      hs[k_z0_lumi]->ShiftFillLast(0., 0., fitNLumi_);
      hs[k_sigmaX0_lumi]->ShiftFillLast(0., 0., fitNLumi_);
      hs[k_sigmaY0_lumi]->ShiftFillLast(0., 0., fitNLumi_);
      hs[k_sigmaZ0_lumi]->ShiftFillLast(0., 0., fitNLumi_);
      h_nVtx_lumi->ShiftFillLast(0., 0., fitNLumi_);
    }
    for (int ig = 0; ig < LSgap_pv; ig++) {
      hs[k_PVx_lumi]->ShiftFillLast(0., 0., fitPVNLumi_);
      hs[k_PVy_lumi]->ShiftFillLast(0., 0., fitPVNLumi_);
      hs[k_PVz_lumi]->ShiftFillLast(0., 0., fitPVNLumi_);
    }
    const int previousLS = h_nTrk_lumi->getTH1()->GetEntries();
    for (int i = 1; i < (currentlumi - previousLS);
         i++)  //if (current-previoius)= 1 then never go inside the for loop!!!!!!!!!!!
      h_nTrk_lumi->ShiftFillLast(nthBSTrk_);
  }
 
  edm::LogInfo("BeamMonitor") << "FitAndFill:: Time lapsed since last scroll = " << tmpTime - refTime << std::endl;
 
  if (testScroll(tmpTime, refTime)) {
    scrollTH1(hs[k_x0_time]->getTH1(), refTime);
    scrollTH1(hs[k_y0_time]->getTH1(), refTime);
    scrollTH1(hs[k_z0_time]->getTH1(), refTime);
    scrollTH1(hs[k_sigmaX0_time]->getTH1(), refTime);
    scrollTH1(hs[k_sigmaY0_time]->getTH1(), refTime);
    scrollTH1(hs[k_sigmaZ0_time]->getTH1(), refTime);
    scrollTH1(hs[k_PVx_time]->getTH1(), refTime);
    scrollTH1(hs[k_PVy_time]->getTH1(), refTime);
    scrollTH1(hs[k_PVz_time]->getTH1(), refTime);
  }
 
  bool doPVFit = false;
 
  if (fitPVNLumi_ > 0) {
    if (onlineMode_) {
      if (currentlumi % fitPVNLumi_ == 0)
        doPVFit = true;
    } else if (countLumi_ % fitPVNLumi_ == 0)
      doPVFit = true;
  } else
    doPVFit = true;
 
  if (doPVFit) {
    edm::LogInfo("BeamMonitor") << "FitAndFill:: Do PV Fitting for LS = " << beginLumiOfPVFit_ << " to "
                                << endLumiOfPVFit_ << std::endl;
    // Primary Vertex Fit:
    if (h_PVx[0]->getTH1()->GetEntries() > minNrVertices_) {
      pvResults->Reset();
      char tmpTitle[50];
      sprintf(
          tmpTitle, "%s %i %s %i", "Fitted Primary Vertex (cm) of LS: ", beginLumiOfPVFit_, " to ", endLumiOfPVFit_);
      pvResults->setAxisTitle(tmpTitle, 1);
 
      std::unique_ptr<TF1> fgaus{new TF1("fgaus", "gaus")};
      double mean, width, meanErr, widthErr;
      fgaus->SetLineColor(4);
      h_PVx[0]->getTH1()->Fit(fgaus.get(), "QLM0");
      mean = fgaus->GetParameter(1);
      width = fgaus->GetParameter(2);
      meanErr = fgaus->GetParError(1);
      widthErr = fgaus->GetParError(2);
 
      hs[k_PVx_lumi]->ShiftFillLast(mean, width, fitPVNLumi_);
      hs[k_PVx_lumi_all]->setBinContent(currentlumi, mean);
      hs[k_PVx_lumi_all]->setBinError(currentlumi, width);
      int nthBin = tmpTime - refTime;
      if (nthBin < 0)
        edm::LogInfo("BeamMonitor") << "FitAndFill::  Event time outside current range of time histograms!"
                                    << std::endl;
      if (nthBin > 0) {
        hs[k_PVx_time]->setBinContent(nthBin, mean);
        hs[k_PVx_time]->setBinError(nthBin, width);
      }
      int jthBin = tmpTime - startTime;
      if (jthBin > 0) {
        hs[k_PVx_time_all]->setBinContent(jthBin, mean);
        hs[k_PVx_time_all]->setBinError(jthBin, width);
      }
      pvResults->setBinContent(1, 6, mean);
      pvResults->setBinContent(1, 3, width);
      pvResults->setBinContent(2, 6, meanErr);
      pvResults->setBinContent(2, 3, widthErr);
 
      {
        // snap shot of the fit
        auto tmphisto = h_PVx[0]->getTH1F();
        h_PVx[1]->getTH1()->SetBins(
            tmphisto->GetNbinsX(), tmphisto->GetXaxis()->GetXmin(), tmphisto->GetXaxis()->GetXmax());
        h_PVx[1]->Reset();
        h_PVx[1]->getTH1()->Add(tmphisto);
        h_PVx[1]->getTH1()->Fit(fgaus.get(), "QLM");
      }
 
      h_PVy[0]->getTH1()->Fit(fgaus.get(), "QLM0");
      mean = fgaus->GetParameter(1);
      width = fgaus->GetParameter(2);
      meanErr = fgaus->GetParError(1);
      widthErr = fgaus->GetParError(2);
      hs[k_PVy_lumi]->ShiftFillLast(mean, width, fitPVNLumi_);
      hs[k_PVy_lumi_all]->setBinContent(currentlumi, mean);
      hs[k_PVy_lumi_all]->setBinError(currentlumi, width);
      if (nthBin > 0) {
        hs[k_PVy_time]->setBinContent(nthBin, mean);
        hs[k_PVy_time]->setBinError(nthBin, width);
      }
      if (jthBin > 0) {
        hs[k_PVy_time_all]->setBinContent(jthBin, mean);
        hs[k_PVy_time_all]->setBinError(jthBin, width);
      }
      pvResults->setBinContent(1, 5, mean);
      pvResults->setBinContent(1, 2, width);
      pvResults->setBinContent(2, 5, meanErr);
      pvResults->setBinContent(2, 2, widthErr);
      // snap shot of the fit
      {
        auto tmphisto = h_PVy[0]->getTH1F();
        h_PVy[1]->getTH1()->SetBins(
            tmphisto->GetNbinsX(), tmphisto->GetXaxis()->GetXmin(), tmphisto->GetXaxis()->GetXmax());
        h_PVy[1]->Reset();
        h_PVy[1]->getTH1()->Add(tmphisto);
        h_PVy[1]->getTH1()->Fit(fgaus.get(), "QLM");
      }
 
      h_PVz[0]->getTH1()->Fit(fgaus.get(), "QLM0");
      mean = fgaus->GetParameter(1);
      width = fgaus->GetParameter(2);
      meanErr = fgaus->GetParError(1);
      widthErr = fgaus->GetParError(2);
      hs[k_PVz_lumi]->ShiftFillLast(mean, width, fitPVNLumi_);
      hs[k_PVz_lumi_all]->setBinContent(currentlumi, mean);
      hs[k_PVz_lumi_all]->setBinError(currentlumi, width);
      if (nthBin > 0) {
        hs[k_PVz_time]->setBinContent(nthBin, mean);
        hs[k_PVz_time]->setBinError(nthBin, width);
      }
      if (jthBin > 0) {
        hs[k_PVz_time_all]->setBinContent(jthBin, mean);
        hs[k_PVz_time_all]->setBinError(jthBin, width);
      }
      pvResults->setBinContent(1, 4, mean);
      pvResults->setBinContent(1, 1, width);
      pvResults->setBinContent(2, 4, meanErr);
      pvResults->setBinContent(2, 1, widthErr);
      {
        // snap shot of the fit
        auto tmphisto = h_PVz[0]->getTH1F();
        h_PVz[1]->getTH1()->SetBins(
            tmphisto->GetNbinsX(), tmphisto->GetXaxis()->GetXmin(), tmphisto->GetXaxis()->GetXmax());
        h_PVz[1]->Reset();
        h_PVz[1]->getTH1()->Add(tmphisto);
        h_PVz[1]->getTH1()->Fit(fgaus.get(), "QLM");
      }
    }  //check if found min Vertices
  }    //do PVfit
 
  if ((resetPVNLumi_ > 0 && countLumi_ == resetPVNLumi_) || StartAverage_) {
    beginLumiOfPVFit_ = 0;
    refPVtime[0] = 0;
  }
 
  //---------Readjustment of theBSvector, RefTime, beginLSofFit---------
  vector<BSTrkParameters> theBSvector1 = theBeamFitter->getBSvector();
  mapLSBSTrkSize[countLumi_] = (theBSvector1.size());
  size_t PreviousRecords = 0;  //needed to fill nth record of tracks in GUI
 
  if (StartAverage_) {
    size_t SizeToRemove = 0;
    std::map<int, std::size_t>::iterator rmls = mapLSBSTrkSize.begin();
    SizeToRemove = rmls->second;
    if (debug_)
      edm::LogInfo("BeamMonitor") << "  The size to remove is =  " << SizeToRemove << endl;
    int changedAfterThis = 0;
    for (std::map<int, std::size_t>::iterator rmLS = mapLSBSTrkSize.begin(); rmLS != mapLSBSTrkSize.end();
         ++rmLS, ++changedAfterThis) {
      if (changedAfterThis > 0) {
        (rmLS->second) = (rmLS->second) - SizeToRemove;
        if ((mapLSBSTrkSize.size() - (size_t)changedAfterThis) == 2)
          PreviousRecords = (rmLS->second);
      }
    }
 
    theBeamFitter->resizeBSvector(SizeToRemove);
 
    map<int, std::size_t>::iterator tmpIt = mapLSBSTrkSize.begin();
    mapLSBSTrkSize.erase(tmpIt);
 
    std::pair<int, int> checkfitLS = theBeamFitter->getFitLSRange();
    std::pair<time_t, time_t> checkfitTime = theBeamFitter->getRefTime();
    theBeamFitter->setFitLSRange(beginLumiOfBSFit_, checkfitLS.second);
    theBeamFitter->setRefTime(refBStime[0], checkfitTime.second);
  }
 
  //Fill the track for this fit
  vector<BSTrkParameters> theBSvector = theBeamFitter->getBSvector();
  h_nTrk_lumi->ShiftFillLast(theBSvector.size());
 
  if (debug_)
    edm::LogInfo("BeamMonitor") << "FitAndFill::   Size of  theBSViector.size()  After =" << theBSvector.size() << endl;
 
  bool countFitting = false;
  if (theBSvector.size() >= PreviousRecords && theBSvector.size() >= min_Ntrks_) {
    countFitting = true;
  }
 
  //---Fix for Cut Flow Table for Running average in a same way//the previous code  has problem for resetting!!!
  mapLSCF[countLumi_] = *theBeamFitter->getCutFlow();
  if (StartAverage_ && !mapLSCF.empty()) {
    const TH1F& cutFlowToSubtract = mapLSCF.begin()->second;
    // Subtract the last cut flow from all of the others.
    std::map<int, TH1F>::iterator cf = mapLSCF.begin();
    // Start on second entry
    for (; cf != mapLSCF.end(); ++cf) {
      cf->second.Add(&cutFlowToSubtract, -1);
    }
    theBeamFitter->subtractFromCutFlow(&cutFlowToSubtract);
    // Remove the obsolete lumi section
    mapLSCF.erase(mapLSCF.begin());
  }
 
  if (resetHistos_) {
    h_d0_phi0->Reset();
    h_vx_vy->Reset();
    h_vx_dz->Reset();
    h_vy_dz->Reset();
    h_trk_z0->Reset();
    resetHistos_ = false;
  }
 
  if (StartAverage_)
    nthBSTrk_ = PreviousRecords;  //after average proccess is ON//for 2-6 LS fit PreviousRecords is size from 2-5 LS
 
  edm::LogInfo("BeamMonitor") << " The Previous Recored for this fit is  =" << nthBSTrk_ << endl;
 
  unsigned int itrk = 0;
  for (vector<BSTrkParameters>::const_iterator BSTrk = theBSvector.begin(); BSTrk != theBSvector.end();
       ++BSTrk, ++itrk) {
    if (itrk >= nthBSTrk_) {  //fill for this record only !!
      h_d0_phi0->Fill(BSTrk->phi0(), BSTrk->d0());
      double vx = BSTrk->vx();
      double vy = BSTrk->vy();
      double z0 = BSTrk->z0();
      h_vx_vy->Fill(vx, vy);
      h_vx_dz->Fill(z0, vx);
      h_vy_dz->Fill(z0, vy);
      h_trk_z0->Fill(z0);
    }
  }
 
  nthBSTrk_ = theBSvector.size();  // keep track of num of tracks filled so far
 
  edm::LogInfo("BeamMonitor") << " The Current Recored for this fit is  =" << nthBSTrk_ << endl;
 
  if (countFitting)
    edm::LogInfo("BeamMonitor") << "FitAndFill::  Num of tracks collected = " << nthBSTrk_ << endl;
 
  if (fitNLumi_ > 0) {
    if (onlineMode_) {
      if (currentlumi % fitNLumi_ != 0) {
        //  for (std::map<TString,MonitorElement*>::iterator itAll = hs.begin();
        //       itAll != hs.end(); ++itAll) {
        //    if ((*itAll).first.Contains("all")) {
        //      (*itAll).second->setBinContent(currentlumi,0.);
        //      (*itAll).second->setBinError(currentlumi,0.);
        //    }
        //  }
        return;
      }
    } else if (countLumi_ % fitNLumi_ != 0)
      return;
  }
 
  edm::LogInfo("BeamMonitor") << "FitAndFill:: [DebugTime] refBStime[0] = " << refBStime[0]
                              << "; address =  " << &refBStime[0] << std::endl;
  edm::LogInfo("BeamMonitor") << "FitAndFill:: [DebugTime] refBStime[1] = " << refBStime[1]
                              << "; address =  " << &refBStime[1] << std::endl;
 
  //Fill for all LS even if fit fails
  h_nVtx_lumi->ShiftFillLast((theBeamFitter->getPVvectorSize()), 0., fitNLumi_);
  h_nVtx_lumi_all->setBinContent(currentlumi, (theBeamFitter->getPVvectorSize()));
 
  if (countFitting) {
    nFits_++;
    std::pair<int, int> fitLS = theBeamFitter->getFitLSRange();
    edm::LogInfo("BeamMonitor") << "FitAndFill::  [BeamFitter] Do BeamSpot Fit for LS = " << fitLS.first << " to "
                                << fitLS.second << std::endl;
    edm::LogInfo("BeamMonitor") << "FitAndFill::  [BeamMonitor] Do BeamSpot Fit for LS = " << beginLumiOfBSFit_
                                << " to " << endLumiOfBSFit_ << std::endl;
 
    //Now Run the PV and Track Fitter over the collected tracks and pvs
    if (theBeamFitter->runPVandTrkFitter()) {
      reco::BeamSpot bs = theBeamFitter->getBeamSpot();
      if (bs.type() > 0)  // with good beamwidth fit
        preBS = bs;       // cache good fit results
 
      edm::LogInfo("BeamMonitor") << "\n RESULTS OF DEFAULT FIT:" << endl;
      edm::LogInfo("BeamMonitor") << bs << endl;
      edm::LogInfo("BeamMonitor") << "[BeamFitter] fitting done \n" << endl;
 
      hs[k_x0_lumi]->ShiftFillLast(bs.x0(), bs.x0Error(), fitNLumi_);
      hs[k_y0_lumi]->ShiftFillLast(bs.y0(), bs.y0Error(), fitNLumi_);
      hs[k_z0_lumi]->ShiftFillLast(bs.z0(), bs.z0Error(), fitNLumi_);
      hs[k_sigmaX0_lumi]->ShiftFillLast(bs.BeamWidthX(), bs.BeamWidthXError(), fitNLumi_);
      hs[k_sigmaY0_lumi]->ShiftFillLast(bs.BeamWidthY(), bs.BeamWidthYError(), fitNLumi_);
      hs[k_sigmaZ0_lumi]->ShiftFillLast(bs.sigmaZ(), bs.sigmaZ0Error(), fitNLumi_);
      hs[k_x0_lumi_all]->setBinContent(currentlumi, bs.x0());
      hs[k_x0_lumi_all]->setBinError(currentlumi, bs.x0Error());
      hs[k_y0_lumi_all]->setBinContent(currentlumi, bs.y0());
      hs[k_y0_lumi_all]->setBinError(currentlumi, bs.y0Error());
      hs[k_z0_lumi_all]->setBinContent(currentlumi, bs.z0());
      hs[k_z0_lumi_all]->setBinError(currentlumi, bs.z0Error());
      hs[k_sigmaX0_lumi_all]->setBinContent(currentlumi, bs.BeamWidthX());
      hs[k_sigmaX0_lumi_all]->setBinError(currentlumi, bs.BeamWidthXError());
      hs[k_sigmaY0_lumi_all]->setBinContent(currentlumi, bs.BeamWidthY());
      hs[k_sigmaY0_lumi_all]->setBinError(currentlumi, bs.BeamWidthYError());
      hs[k_sigmaZ0_lumi_all]->setBinContent(currentlumi, bs.sigmaZ());
      hs[k_sigmaZ0_lumi_all]->setBinError(currentlumi, bs.sigmaZ0Error());
 
      int nthBin = tmpTime - refTime;
      if (nthBin > 0) {
        hs[k_x0_time]->setBinContent(nthBin, bs.x0());
        hs[k_y0_time]->setBinContent(nthBin, bs.y0());
        hs[k_z0_time]->setBinContent(nthBin, bs.z0());
        hs[k_sigmaX0_time]->setBinContent(nthBin, bs.BeamWidthX());
        hs[k_sigmaY0_time]->setBinContent(nthBin, bs.BeamWidthY());
        hs[k_sigmaZ0_time]->setBinContent(nthBin, bs.sigmaZ());
        hs[k_x0_time]->setBinError(nthBin, bs.x0Error());
        hs[k_y0_time]->setBinError(nthBin, bs.y0Error());
        hs[k_z0_time]->setBinError(nthBin, bs.z0Error());
        hs[k_sigmaX0_time]->setBinError(nthBin, bs.BeamWidthXError());
        hs[k_sigmaY0_time]->setBinError(nthBin, bs.BeamWidthYError());
        hs[k_sigmaZ0_time]->setBinError(nthBin, bs.sigmaZ0Error());
      }
 
      int jthBin = tmpTime - startTime;
      if (jthBin > 0) {
        hs[k_x0_time_all]->setBinContent(jthBin, bs.x0());
        hs[k_y0_time_all]->setBinContent(jthBin, bs.y0());
        hs[k_z0_time_all]->setBinContent(jthBin, bs.z0());
        hs[k_sigmaX0_time_all]->setBinContent(jthBin, bs.BeamWidthX());
        hs[k_sigmaY0_time_all]->setBinContent(jthBin, bs.BeamWidthY());
        hs[k_sigmaZ0_time_all]->setBinContent(jthBin, bs.sigmaZ());
        hs[k_x0_time_all]->setBinError(jthBin, bs.x0Error());
        hs[k_y0_time_all]->setBinError(jthBin, bs.y0Error());
        hs[k_z0_time_all]->setBinError(jthBin, bs.z0Error());
        hs[k_sigmaX0_time_all]->setBinError(jthBin, bs.BeamWidthXError());
        hs[k_sigmaY0_time_all]->setBinError(jthBin, bs.BeamWidthYError());
        hs[k_sigmaZ0_time_all]->setBinError(jthBin, bs.sigmaZ0Error());
      }
 
      h_x0->Fill(bs.x0());
      h_y0->Fill(bs.y0());
      h_z0->Fill(bs.z0());
      if (bs.type() > 0) {  // with good beamwidth fit
        h_sigmaX0->Fill(bs.BeamWidthX());
        h_sigmaY0->Fill(bs.BeamWidthY());
      }
      h_sigmaZ0->Fill(bs.sigmaZ());
 
      if (nthBSTrk_ >= 2 * min_Ntrks_) {
        double amp = std::sqrt(bs.x0() * bs.x0() + bs.y0() * bs.y0());
        double alpha = std::atan2(bs.y0(), bs.x0());
        std::unique_ptr<TF1> f1{new TF1("f1", "[0]*sin(x-[1])", -3.14, 3.14)};
        f1->SetParameters(amp, alpha);
        f1->SetParLimits(0, amp - 0.1, amp + 0.1);
        f1->SetParLimits(1, alpha - 0.577, alpha + 0.577);
        f1->SetLineColor(4);
        h_d0_phi0->getTProfile()->Fit(f1.get(), "QR");
 
        double mean = bs.z0();
        double width = bs.sigmaZ();
        std::unique_ptr<TF1> fgaus{new TF1("fgaus", "gaus")};
        fgaus->SetParameters(mean, width);
        fgaus->SetLineColor(4);
        h_trk_z0->getTH1()->Fit(fgaus.get(), "QLRM", "", mean - 3 * width, mean + 3 * width);
      }
 
      fitResults->Reset();
      std::pair<int, int> LSRange = theBeamFitter->getFitLSRange();
      char tmpTitle[50];
      sprintf(tmpTitle, "%s %i %s %i", "Fitted Beam Spot (cm) of LS: ", LSRange.first, " to ", LSRange.second);
      fitResults->setAxisTitle(tmpTitle, 1);
      fitResults->setBinContent(1, 8, bs.x0());
      fitResults->setBinContent(1, 7, bs.y0());
      fitResults->setBinContent(1, 6, bs.z0());
      fitResults->setBinContent(1, 5, bs.sigmaZ());
      fitResults->setBinContent(1, 4, bs.dxdz());
      fitResults->setBinContent(1, 3, bs.dydz());
      if (bs.type() > 0) {  // with good beamwidth fit
        fitResults->setBinContent(1, 2, bs.BeamWidthX());
        fitResults->setBinContent(1, 1, bs.BeamWidthY());
      } else {  // fill cached widths
        fitResults->setBinContent(1, 2, preBS.BeamWidthX());
        fitResults->setBinContent(1, 1, preBS.BeamWidthY());
      }
 
      fitResults->setBinContent(2, 8, bs.x0Error());
      fitResults->setBinContent(2, 7, bs.y0Error());
      fitResults->setBinContent(2, 6, bs.z0Error());
      fitResults->setBinContent(2, 5, bs.sigmaZ0Error());
      fitResults->setBinContent(2, 4, bs.dxdzError());
      fitResults->setBinContent(2, 3, bs.dydzError());
      if (bs.type() > 0) {  // with good beamwidth fit
        fitResults->setBinContent(2, 2, bs.BeamWidthXError());
        fitResults->setBinContent(2, 1, bs.BeamWidthYError());
      } else {  // fill cached width errors
        fitResults->setBinContent(2, 2, preBS.BeamWidthXError());
        fitResults->setBinContent(2, 1, preBS.BeamWidthYError());
      }
 
      // count good fit
      //     if (std::fabs(refBS.x0()-bs.x0())/bs.x0Error() < deltaSigCut_) { // disabled temporarily
      summaryContent_[0] += 1.;
      //     }
      //     if (std::fabs(refBS.y0()-bs.y0())/bs.y0Error() < deltaSigCut_) { // disabled temporarily
      summaryContent_[1] += 1.;
      //     }
      //     if (std::fabs(refBS.z0()-bs.z0())/bs.z0Error() < deltaSigCut_) { // disabled temporarily
      summaryContent_[2] += 1.;
      //     }
 
      // Create the BeamSpotOnlineObjects object
      BeamSpotOnlineObjects* BSOnline = new BeamSpotOnlineObjects();
      BSOnline->SetLastAnalyzedLumi(fitLS.second);
      BSOnline->SetLastAnalyzedRun(theBeamFitter->getRunNumber());
      BSOnline->SetLastAnalyzedFill(0);  // To be updated with correct LHC Fill number
      BSOnline->SetPosition(bs.x0(), bs.y0(), bs.z0());
      BSOnline->SetSigmaZ(bs.sigmaZ());
      BSOnline->SetBeamWidthX(bs.BeamWidthX());
      BSOnline->SetBeamWidthY(bs.BeamWidthY());
      BSOnline->SetBeamWidthXError(bs.BeamWidthXError());
      BSOnline->SetBeamWidthYError(bs.BeamWidthYError());
      BSOnline->Setdxdz(bs.dxdz());
      BSOnline->Setdydz(bs.dydz());
      BSOnline->SetType(bs.type());
      BSOnline->SetEmittanceX(bs.emittanceX());
      BSOnline->SetEmittanceY(bs.emittanceY());
      BSOnline->SetBetaStar(bs.betaStar());
      for (int i = 0; i < 7; ++i) {
        for (int j = 0; j < 7; ++j) {
          BSOnline->SetCovariance(i, j, bs.covariance(i, j));
        }
      }
      BSOnline->SetNumTracks(theBeamFitter->getNTracks());
      BSOnline->SetNumPVs(theBeamFitter->getNPVs());
      auto creationTime =
          std::chrono::duration_cast<std::chrono::microseconds>(std::chrono::system_clock::now().time_since_epoch())
              .count();
      BSOnline->SetCreationTime(creationTime);
 
      edm::LogInfo("BeamMonitor") << "FitAndFill::[PayloadCreation] BeamSpotOnline object created: \n" << std::endl;
      edm::LogInfo("BeamMonitor") << *BSOnline << std::endl;
 
      // Create the payload for BeamSpotOnlineObjects object
      if (onlineDbService_.isAvailable()) {
        edm::LogInfo("BeamMonitor") << "FitAndFill::[PayloadCreation] onlineDbService available \n" << std::endl;
        onlineDbService_->logger().logInfo() << "BeamMonitor::FitAndFill - Lumi of the current fit: " << currentlumi;
        onlineDbService_->logger().logInfo()
            << "BeamMonitor::FitAndFill - Do PV Fitting for LS = " << beginLumiOfPVFit_ << " to " << endLumiOfPVFit_;
        onlineDbService_->logger().logInfo()
            << "BeamMonitor::FitAndFill - [BeamFitter] Do BeamSpot Fit for LS = " << fitLS.first << " to "
            << fitLS.second;
        onlineDbService_->logger().logInfo()
            << "BeamMonitor::FitAndFill - [BeamMonitor] Do BeamSpot Fit for LS = " << beginLumiOfBSFit_ << " to "
            << endLumiOfBSFit_;
        onlineDbService_->logger().logInfo() << "BeamMonitor - RESULTS OF DEFAULT FIT:";
        onlineDbService_->logger().logInfo() << "\n" << bs;
        onlineDbService_->logger().logInfo()
            << "BeamMonitor::FitAndFill - [PayloadCreation] BeamSpotOnline object created:";
        onlineDbService_->logger().logInfo() << "\n" << *BSOnline;
        onlineDbService_->logger().logInfo() << "BeamMonitor::FitAndFill - [PayloadCreation] onlineDbService available";
        onlineDbService_->logger().logInfo()
            << "BeamMonitor::FitAndFill - [PayloadCreation] SetCreationTime: " << creationTime
            << " [epoch in microseconds]";
        try {
          onlineDbService_->writeForNextLumisection(BSOnline, recordName_);
          onlineDbService_->logger().logInfo()
              << "BeamMonitor::FitAndFill - [PayloadCreation] writeForNextLumisection executed correctly";
        } catch (const std::exception& e) {
          onlineDbService_->logger().logError() << "BeamMonitor - Error writing record: " << recordName_
                                                << " for Run: " << frun << " - Lumi: " << fitLS.second;
          onlineDbService_->logger().logError() << "Error is: " << e.what();
          onlineDbService_->logger().logError() << "RESULTS OF DEFAULT FIT WAS:";
          onlineDbService_->logger().logError() << "\n" << bs;
          DBloggerReturn_ = 2;
        }
      }
      edm::LogInfo("BeamMonitor") << "FitAndFill::[PayloadCreation] BeamSpotOnline payload created \n" << std::endl;
 
    }       //if (theBeamFitter->runPVandTrkFitter())
    else {  // beam fit fails
      reco::BeamSpot bs = theBeamFitter->getBeamSpot();
      edm::LogInfo("BeamMonitor") << "FitAndFill::   [BeamMonitor] Beam fit fails!!! \n" << endl;
      edm::LogInfo("BeamMonitor") << "FitAndFill::   [BeamMonitor] Output beam spot for DIP \n" << endl;
      edm::LogInfo("BeamMonitor") << bs << endl;
 
      if (onlineDbService_.isAvailable()) {
        onlineDbService_->logger().logInfo() << "BeamMonitor::FitAndFill - Beam fit fails!!!";
        onlineDbService_->logger().logInfo() << "BeamMonitor::FitAndFill - Output beam spot for DIP";
        onlineDbService_->logger().logInfo() << "\n" << bs;
        DBloggerReturn_ = 2;
      }
 
      hs[k_sigmaX0_lumi]->ShiftFillLast(bs.BeamWidthX(), bs.BeamWidthXError(), fitNLumi_);
      hs[k_sigmaY0_lumi]->ShiftFillLast(bs.BeamWidthY(), bs.BeamWidthYError(), fitNLumi_);
      hs[k_sigmaZ0_lumi]->ShiftFillLast(bs.sigmaZ(), bs.sigmaZ0Error(), fitNLumi_);
      hs[k_x0_lumi]->ShiftFillLast(bs.x0(), bs.x0Error(), fitNLumi_);
      hs[k_y0_lumi]->ShiftFillLast(bs.y0(), bs.y0Error(), fitNLumi_);
      hs[k_z0_lumi]->ShiftFillLast(bs.z0(), bs.z0Error(), fitNLumi_);
    }  // end of beam fit fails
 
  }       //-------- end of countFitting------------------------------------------
  else {  // no fit
    // Overwrite Fit LS and fit time when no event processed or no track selected
    theBeamFitter->setFitLSRange(beginLumiOfBSFit_, endLumiOfBSFit_);
    theBeamFitter->setRefTime(refBStime[0], refBStime[1]);
    if (theBeamFitter->runPVandTrkFitter()) {
    }  // Dump fake beam spot for DIP
    reco::BeamSpot bs = theBeamFitter->getBeamSpot();
    edm::LogInfo("BeamMonitor") << "FitAndFill::  [BeamMonitor] No fitting \n" << endl;
    edm::LogInfo("BeamMonitor") << "FitAndFill::  [BeamMonitor] Output fake beam spot for DIP \n" << endl;
    edm::LogInfo("BeamMonitor") << bs << endl;
 
    if (onlineDbService_.isAvailable()) {
      onlineDbService_->logger().logInfo() << "BeamMonitor::FitAndFill - No fitting";
      onlineDbService_->logger().logInfo() << "BeamMonitor::FitAndFill - Output fake beam spot for DIP";
      onlineDbService_->logger().logInfo() << "\n" << bs;
      DBloggerReturn_ = 2;
    }
 
    hs[k_sigmaX0_lumi]->ShiftFillLast(bs.BeamWidthX(), bs.BeamWidthXError(), fitNLumi_);
    hs[k_sigmaY0_lumi]->ShiftFillLast(bs.BeamWidthY(), bs.BeamWidthYError(), fitNLumi_);
    hs[k_sigmaZ0_lumi]->ShiftFillLast(bs.sigmaZ(), bs.sigmaZ0Error(), fitNLumi_);
    hs[k_x0_lumi]->ShiftFillLast(bs.x0(), bs.x0Error(), fitNLumi_);
    hs[k_y0_lumi]->ShiftFillLast(bs.y0(), bs.y0Error(), fitNLumi_);
    hs[k_z0_lumi]->ShiftFillLast(bs.z0(), bs.z0Error(), fitNLumi_);
  }
 
  // Fill summary report
  if (countFitting) {
    for (int n = 0; n < nFitElements_; n++) {
      reportSummaryContents[n]->Fill(summaryContent_[n] / (float)nFits_);
    }
 
    summarySum_ = 0;
    for (int ii = 0; ii < nFitElements_; ii++) {
      summarySum_ += summaryContent_[ii];
    }
    reportSummary_ = summarySum_ / (nFitElements_ * nFits_);
    if (reportSummary)
      reportSummary->Fill(reportSummary_);
 
    for (int bi = 0; bi < nFitElements_; bi++) {
      reportSummaryMap->setBinContent(1, bi + 1, summaryContent_[bi] / (float)nFits_);
    }
  }
 
  if ((resetFitNLumi_ > 0 &&
       ((onlineMode_ &&
         countLumi_ == resetFitNLumi_) ||  //OR it should be currentLumi_ (if in sequence then does not mattar)
        (!onlineMode_ && countLumi_ == resetFitNLumi_))) ||
      (StartAverage_)) {
    edm::LogInfo("BeamMonitor") << "FitAndFill:: The flag is ON for running average Beam Spot fit" << endl;
    StartAverage_ = true;
    firstAverageFit_++;
    resetHistos_ = true;
    nthBSTrk_ = 0;
    beginLumiOfBSFit_ = 0;
    refBStime[0] = 0;
  }
}

References alpha, reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthXError(), reco::BeamSpot::BeamWidthY(), reco::BeamSpot::BeamWidthYError(), beginLumiOfBSFit_, beginLumiOfPVFit_, cms::cuda::bs, countLumi_, BeamSpotPI::creationTime, DBloggerReturn_, debug_, MillePedeFileConverter_cfg::e, mps_fire::end, endLumiOfBSFit_, endLumiOfPVFit_, cppFunctionSkipper::exception, DeadROC_duringRun::f1, dqm::impl::MonitorElement::Fill(), firstAverageFit_, fitNLumi_, fitPVNLumi_, fitResults, frun, dqm::impl::MonitorElement::getMean(), dqm::impl::MonitorElement::getMeanError(), dqm::impl::MonitorElement::getRMS(), dqm::impl::MonitorElement::getRMSError(), dqm::impl::MonitorElement::getTH1(), dqm::impl::MonitorElement::getTH1F(), dqm::impl::MonitorElement::getTProfile(), h_d0_phi0, h_nTrk_lumi, h_nVtx_lumi, h_nVtx_lumi_all, h_nVtx_st, h_PVx, h_PVy, h_PVz, h_sigmaX0, h_sigmaY0, h_sigmaZ0, h_trk_z0, h_vx_dz, h_vx_vy, h_vy_dz, h_x0, h_y0, h_z0, hs, mps_fire::i, cuy::ii, edm::Service< T >::isAvailable(), dqmiolumiharvest::j, cond::persistency::Logger::logError(), cond::service::PoolDBOutputService::logger(), cond::persistency::Logger::logInfo(), genParticles_cff::map, mapLSBSTrkSize, mapLSCF, mapLSPVStoreSize, mapNPV, mapPVx, mapPVy, mapPVz, SiStripPI::mean, min_Ntrks_, minNrVertices_, dqmiodumpmetadata::n, nFitElements_, nFits_, submitPVValidationJobs::now, nthBSTrk_, onlineDbService_, onlineMode_, preBS, processed_, pvResults, recordName_, refBStime, refPVtime, refTime, reportSummary, reportSummary_, reportSummaryContents, reportSummaryMap, dqm::impl::MonitorElement::Reset(), resetFitNLumi_, resetHistos_, resetPVNLumi_, scrollTH1(), dqm::impl::MonitorElement::setAxisTitle(), BeamSpotObjects::SetBeamWidthX(), BeamSpotObjects::SetBeamWidthXError(), BeamSpotObjects::SetBeamWidthY(), BeamSpotObjects::SetBeamWidthYError(), BeamSpotObjects::SetBetaStar(), dqm::impl::MonitorElement::setBinContent(), BeamSpotObjects::SetCovariance(), BeamSpotOnlineObjects::SetCreationTime(), BeamSpotObjects::Setdxdz(), BeamSpotObjects::Setdydz(), BeamSpotObjects::SetEmittanceX(), BeamSpotObjects::SetEmittanceY(), BeamSpotOnlineObjects::SetLastAnalyzedFill(), BeamSpotOnlineObjects::SetLastAnalyzedLumi(), BeamSpotOnlineObjects::SetLastAnalyzedRun(), BeamSpotOnlineObjects::SetNumPVs(), BeamSpotOnlineObjects::SetNumTracks(), BeamSpotObjects::SetPosition(), BeamSpotObjects::SetSigmaZ(), BeamSpotObjects::SetType(), dqm::impl::MonitorElement::ShiftFillLast(), mathSSE::sqrt(), StartAverage_, startTime, summaryContent_, summarySum_, testScroll(), theBeamFitter, tmpTime, ApeEstimator_cff::width, cond::service::OnlineDBOutputService::writeForNextLumisection(), and HLTMuonOfflineAnalyzer_cfi::z0.

Referenced by beginLuminosityBlock().

formatFitTime()

void BeamMonitor::formatFitTime ( char *  , const std::time_t &    )

private

Definition at line 43 of file BeamMonitor.cc.

                                                         {
  //constexpr int CET(+1);
  constexpr int CEST(+2);
 
  //tm * ptm;
  //ptm = gmtime ( &t );
  //int year = ptm->tm_year;
 
  //get correct year from ctime
  time_t currentTime;
  struct tm* localTime;
  time(&currentTime);                   // Get the current time
  localTime = localtime(&currentTime);  // Convert the current time to the local time
  int year = localTime->tm_year + 1900;
 
  tm* ptm;
  ptm = gmtime(&t);
 
  //check if year is ok
  if (year <= 37)
    year += 2000;
  if (year >= 70 && year <= 137)
    year += 1900;
 
  if (year < 1995) {
    edm::LogError("BadTimeStamp") << "year reported is " << year << " !!" << std::endl;
    //year = 2015; //overwritten later by BeamFitter.cc for fits but needed here for TH1
    //edm::LogError("BadTimeStamp") << "Resetting to " <<year<<std::endl;
  }
  sprintf(ts,
          "%4d-%02d-%02d %02d:%02d:%02d",
          year,
          ptm->tm_mon + 1,
          ptm->tm_mday,
          (ptm->tm_hour + CEST) % 24,
          ptm->tm_min,
          ptm->tm_sec);
 
#ifdef STRIP_TRAILING_BLANKS_IN_TIMEZONE
  unsigned int b = strlen(ts);
  while (ts[--b] == ' ') {
    ts[b] = 0;
  }
#endif
}

References b, CEST, plotBeamSpotDB::ptm, submitPVValidationJobs::t, and protons_cff::time.

Referenced by bookHistograms(), and scrollTH1().

RestartFitting()

void BeamMonitor::RestartFitting ( )

private

Definition at line 1494 of file BeamMonitor.cc.

                                 {
  if (debug_)
    edm::LogInfo("BeamMonitor")
        << " RestartingFitting:: Restart Beami everything to a fresh start !!! because Gap is > 10 LS" << endl;
  //track based fit reset here
  resetHistos_ = true;
  nthBSTrk_ = 0;
  theBeamFitter->resetTrkVector();
  theBeamFitter->resetLSRange();
  theBeamFitter->resetRefTime();
  theBeamFitter->resetPVFitter();
  theBeamFitter->resetCutFlow();
  beginLumiOfBSFit_ = 0;
  refBStime[0] = 0;
  //pv based fit iis reset here
  h_PVx[0]->Reset();
  h_PVy[0]->Reset();
  h_PVz[0]->Reset();
  beginLumiOfPVFit_ = 0;
  refPVtime[0] = 0;
  //Clear all the Maps here
  mapPVx.clear();
  mapPVy.clear();
  mapPVz.clear();
  mapNPV.clear();
  mapBeginBSLS.clear();
  mapBeginPVLS.clear();
  mapBeginBSTime.clear();
  mapBeginPVTime.clear();
  mapLSBSTrkSize.clear();
  mapLSPVStoreSize.clear();
  mapLSCF.clear();
  countGapLumi_ = 0;
  countLumi_ = 0;
  StartAverage_ = false;
}

References beginLumiOfBSFit_, beginLumiOfPVFit_, countGapLumi_, countLumi_, debug_, h_PVx, h_PVy, h_PVz, mapBeginBSLS, mapBeginBSTime, mapBeginPVLS, mapBeginPVTime, mapLSBSTrkSize, mapLSCF, mapLSPVStoreSize, mapNPV, mapPVx, mapPVy, mapPVz, nthBSTrk_, refBStime, refPVtime, dqm::impl::MonitorElement::Reset(), resetHistos_, StartAverage_, and theBeamFitter.

scrollTH1()

void BeamMonitor::scrollTH1 ( TH1 *  , std::time_t    )

private

Definition at line 1550 of file BeamMonitor.cc.

                                              {
  char offsetTime[64];
  formatFitTime(offsetTime, ref);
  TDatime da(offsetTime);
  if (lastNZbin > 0) {
    double val = h->GetBinContent(lastNZbin);
    double valErr = h->GetBinError(lastNZbin);
    h->Reset();
    h->GetXaxis()->SetTimeOffset(da.Convert(kTRUE));
    int bin = (lastNZbin > buffTime ? buffTime : 1);
    h->SetBinContent(bin, val);
    h->SetBinError(bin, valErr);
  } else {
    h->Reset();
    h->GetXaxis()->SetTimeOffset(da.Convert(kTRUE));
  }
}

References newFWLiteAna::bin, buffTime, Calorimetry_cff::da, formatFitTime(), lastNZbin, and heppy_batch::val.

Referenced by FitAndFill().

testScroll()

bool BeamMonitor::testScroll ( std::time_t &  , std::time_t &    )

private

Definition at line 1570 of file BeamMonitor.cc.

                                                               {
  bool scroll_ = false;
  if (tmpTime_ - refTime_ >= intervalInSec_) {
    scroll_ = true;
    edm::LogInfo("BeamMonitor") << "testScroll::  Reset Time Offset" << std::endl;
    lastNZbin = intervalInSec_;
    for (int bin = intervalInSec_; bin >= 1; bin--) {
      if (hs[k_x0_time]->getBinContent(bin) > 0) {
        lastNZbin = bin;
        break;
      }
    }
    edm::LogInfo("BeamMonitor") << "testScroll::  Last non zero bin = " << lastNZbin << std::endl;
    if (tmpTime_ - refTime_ >= intervalInSec_ + lastNZbin) {
      edm::LogInfo("BeamMonitor") << "testScroll::  Time difference too large since last readout" << std::endl;
      lastNZbin = 0;
      refTime_ = tmpTime_ - buffTime;
    } else {
      edm::LogInfo("BeamMonitor") << "testScroll::  Offset to last record" << std::endl;
      int offset = ((lastNZbin > buffTime) ? (lastNZbin - buffTime) : (lastNZbin - 1));
      refTime_ += offset;
    }
  }
  return scroll_;
}

References newFWLiteAna::bin, buffTime, hs, intervalInSec_, lastNZbin, and hltrates_dqm_sourceclient-live_cfg::offset.

Referenced by FitAndFill().

Member Data Documentation

beginLumiOfBSFit_

int BeamMonitor::beginLumiOfBSFit_

private

Definition at line 93 of file BeamMonitor.h.

Referenced by BeamMonitor(), beginLuminosityBlock(), FitAndFill(), and RestartFitting().

beginLumiOfPVFit_

int BeamMonitor::beginLumiOfPVFit_

private

Definition at line 95 of file BeamMonitor.h.

Referenced by BeamMonitor(), beginLuminosityBlock(), FitAndFill(), and RestartFitting().

bsSrc_

edm::EDGetTokenT<reco::BeamSpot> BeamMonitor::bsSrc_

private

Definition at line 73 of file BeamMonitor.h.

Referenced by analyze(), and BeamMonitor().

countEvt_

int BeamMonitor::countEvt_

private

Definition at line 91 of file BeamMonitor.h.

Referenced by analyze().

countGapLumi_

int BeamMonitor::countGapLumi_

private

Definition at line 114 of file BeamMonitor.h.

Referenced by RestartFitting().

countLumi_

int BeamMonitor::countLumi_

private

Definition at line 92 of file BeamMonitor.h.

Referenced by analyze(), beginLuminosityBlock(), FitAndFill(), and RestartFitting().

cutFlowTable

MonitorElement* BeamMonitor::cutFlowTable

private

Definition at line 171 of file BeamMonitor.h.

Referenced by analyze(), and bookHistograms().

DBloggerReturn_

int BeamMonitor::DBloggerReturn_

private

Definition at line 78 of file BeamMonitor.h.

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

debug_

bool BeamMonitor::debug_

private

Definition at line 85 of file BeamMonitor.h.

Referenced by BeamMonitor(), beginLuminosityBlock(), bookHistograms(), dqmEndRun(), FitAndFill(), and RestartFitting().

deltaSigCut_

double BeamMonitor::deltaSigCut_

private

Definition at line 104 of file BeamMonitor.h.

Referenced by BeamMonitor().

dxBin_

const int BeamMonitor::dxBin_

private

Definition at line 56 of file BeamMonitor.h.

Referenced by bookHistograms().

dxMax_

const double BeamMonitor::dxMax_

private

Definition at line 58 of file BeamMonitor.h.

Referenced by bookHistograms().

dxMin_

const double BeamMonitor::dxMin_

private

Definition at line 57 of file BeamMonitor.h.

Referenced by bookHistograms().

dzBin_

const int BeamMonitor::dzBin_

private

Definition at line 68 of file BeamMonitor.h.

Referenced by bookHistograms().

dzMax_

const double BeamMonitor::dzMax_

private

Definition at line 70 of file BeamMonitor.h.

Referenced by bookHistograms().

dzMin_

const double BeamMonitor::dzMin_

private

Definition at line 69 of file BeamMonitor.h.

Referenced by bookHistograms().

endLumiOfBSFit_

int BeamMonitor::endLumiOfBSFit_

private

Definition at line 94 of file BeamMonitor.h.

Referenced by BeamMonitor(), and FitAndFill().

endLumiOfPVFit_

int BeamMonitor::endLumiOfPVFit_

private

Definition at line 96 of file BeamMonitor.h.

Referenced by BeamMonitor(), and FitAndFill().

firstAverageFit_

int BeamMonitor::firstAverageFit_

private

Definition at line 113 of file BeamMonitor.h.

Referenced by FitAndFill().

fitNLumi_

int BeamMonitor::fitNLumi_

private

Definition at line 80 of file BeamMonitor.h.

Referenced by BeamMonitor(), and FitAndFill().

fitPVNLumi_

int BeamMonitor::fitPVNLumi_

private

Definition at line 81 of file BeamMonitor.h.

Referenced by BeamMonitor(), and FitAndFill().

fitResults

MonitorElement* BeamMonitor::fitResults

private

Definition at line 135 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

frun

int BeamMonitor::frun

private

Definition at line 180 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

ftimestamp

edm::TimeValue_t BeamMonitor::ftimestamp

private

Definition at line 178 of file BeamMonitor.h.

Referenced by bookHistograms().

h_d0_phi0

MonitorElement* BeamMonitor::h_d0_phi0

private

Definition at line 128 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

h_nTrk_lumi

MonitorElement* BeamMonitor::h_nTrk_lumi

private

Definition at line 125 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

h_nVtx

MonitorElement* BeamMonitor::h_nVtx

private

Definition at line 142 of file BeamMonitor.h.

Referenced by analyze(), and bookHistograms().

h_nVtx_lumi

MonitorElement* BeamMonitor::h_nVtx_lumi

private

Definition at line 126 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

h_nVtx_lumi_all

MonitorElement* BeamMonitor::h_nVtx_lumi_all

private

Definition at line 127 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

h_nVtx_st

MonitorElement* BeamMonitor::h_nVtx_st

private

Definition at line 143 of file BeamMonitor.h.

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

h_PVx

MonitorElement* BeamMonitor::h_PVx[2]

private

Definition at line 144 of file BeamMonitor.h.

Referenced by analyze(), bookHistograms(), FitAndFill(), and RestartFitting().

h_PVxz

MonitorElement* BeamMonitor::h_PVxz

private

Definition at line 147 of file BeamMonitor.h.

Referenced by analyze(), and bookHistograms().

h_PVy

MonitorElement* BeamMonitor::h_PVy[2]

private

Definition at line 145 of file BeamMonitor.h.

Referenced by analyze(), bookHistograms(), FitAndFill(), and RestartFitting().

h_PVyz

MonitorElement* BeamMonitor::h_PVyz

private

Definition at line 148 of file BeamMonitor.h.

Referenced by analyze(), and bookHistograms().

h_PVz

MonitorElement* BeamMonitor::h_PVz[2]

private

Definition at line 146 of file BeamMonitor.h.

Referenced by analyze(), bookHistograms(), FitAndFill(), and RestartFitting().

h_sigmaX0

MonitorElement* BeamMonitor::h_sigmaX0

private

Definition at line 139 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

h_sigmaY0

MonitorElement* BeamMonitor::h_sigmaY0

private

Definition at line 140 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

h_sigmaZ0

MonitorElement* BeamMonitor::h_sigmaZ0

private

Definition at line 141 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

h_trk_z0

MonitorElement* BeamMonitor::h_trk_z0

private

Definition at line 129 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

h_trkPt

MonitorElement* BeamMonitor::h_trkPt

private

Definition at line 133 of file BeamMonitor.h.

Referenced by analyze(), and bookHistograms().

h_trkVz

MonitorElement* BeamMonitor::h_trkVz

private

Definition at line 134 of file BeamMonitor.h.

Referenced by analyze(), and bookHistograms().

h_vx_dz

MonitorElement* BeamMonitor::h_vx_dz

private

Definition at line 131 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

h_vx_vy

MonitorElement* BeamMonitor::h_vx_vy

private

Definition at line 130 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

h_vy_dz

MonitorElement* BeamMonitor::h_vy_dz

private

Definition at line 132 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

h_x0

MonitorElement* BeamMonitor::h_x0

private

Definition at line 136 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

h_y0

MonitorElement* BeamMonitor::h_y0

private

Definition at line 137 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

h_z0

MonitorElement* BeamMonitor::h_z0

private

Definition at line 138 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

hltSrc_

edm::EDGetTokenT<edm::TriggerResults> BeamMonitor::hltSrc_

private

Definition at line 76 of file BeamMonitor.h.

Referenced by analyze(), and BeamMonitor().

hs

std::vector<MonitorElement*> BeamMonitor::hs

private

Definition at line 150 of file BeamMonitor.h.

Referenced by bookHistograms(), FitAndFill(), and testScroll().

intervalInSec_

int BeamMonitor::intervalInSec_

private

Definition at line 84 of file BeamMonitor.h.

Referenced by BeamMonitor(), bookHistograms(), and testScroll().

jetTrigger_

std::vector<std::string> BeamMonitor::jetTrigger_

private

Definition at line 87 of file BeamMonitor.h.

Referenced by analyze(), and BeamMonitor().

lastlumi_

int BeamMonitor::lastlumi_

private

Definition at line 97 of file BeamMonitor.h.

Referenced by BeamMonitor(), and beginLuminosityBlock().

lastNZbin

int BeamMonitor::lastNZbin

private

Definition at line 181 of file BeamMonitor.h.

Referenced by scrollTH1(), and testScroll().

mapBeginBSLS

std::map<int, int> BeamMonitor::mapBeginBSLS

private

Definition at line 156 of file BeamMonitor.h.

Referenced by beginLuminosityBlock(), dqmEndRun(), and RestartFitting().

mapBeginBSTime

std::map<int, std::time_t> BeamMonitor::mapBeginBSTime

private

Definition at line 157 of file BeamMonitor.h.

Referenced by beginLuminosityBlock(), dqmEndRun(), and RestartFitting().

mapBeginPVLS

std::map<int, int> BeamMonitor::mapBeginPVLS

private

Definition at line 156 of file BeamMonitor.h.

Referenced by beginLuminosityBlock(), dqmEndRun(), and RestartFitting().

mapBeginPVTime

std::map<int, std::time_t> BeamMonitor::mapBeginPVTime

private

Definition at line 157 of file BeamMonitor.h.

Referenced by beginLuminosityBlock(), dqmEndRun(), and RestartFitting().

mapLSBSTrkSize

std::map<int, std::size_t> BeamMonitor::mapLSBSTrkSize

private

Definition at line 159 of file BeamMonitor.h.

Referenced by dqmEndRun(), FitAndFill(), and RestartFitting().

mapLSCF

std::map<int, TH1F> BeamMonitor::mapLSCF

private

Definition at line 162 of file BeamMonitor.h.

Referenced by dqmEndRun(), FitAndFill(), and RestartFitting().

mapLSPVStoreSize

std::map<int, size_t> BeamMonitor::mapLSPVStoreSize

private

Definition at line 160 of file BeamMonitor.h.

Referenced by dqmEndRun(), FitAndFill(), and RestartFitting().

mapNPV

std::map<int, std::vector<int> > BeamMonitor::mapNPV

private

Definition at line 154 of file BeamMonitor.h.

Referenced by analyze(), dqmEndRun(), FitAndFill(), and RestartFitting().

mapPVx

std::map<int, std::vector<float> > BeamMonitor::mapPVx

private

Definition at line 153 of file BeamMonitor.h.

Referenced by analyze(), dqmEndRun(), FitAndFill(), and RestartFitting().

mapPVy

std::map<int, std::vector<float> > BeamMonitor::mapPVy

private

Definition at line 153 of file BeamMonitor.h.

Referenced by analyze(), dqmEndRun(), FitAndFill(), and RestartFitting().

mapPVz

std::map<int, std::vector<float> > BeamMonitor::mapPVz

private

Definition at line 153 of file BeamMonitor.h.

Referenced by analyze(), dqmEndRun(), FitAndFill(), and RestartFitting().

maxZ_

double BeamMonitor::maxZ_

private

Definition at line 106 of file BeamMonitor.h.

Referenced by BeamMonitor().

min_Ntrks_

unsigned int BeamMonitor::min_Ntrks_

private

Definition at line 105 of file BeamMonitor.h.

Referenced by BeamMonitor(), and FitAndFill().

minNrVertices_

unsigned int BeamMonitor::minNrVertices_

private

Definition at line 107 of file BeamMonitor.h.

Referenced by BeamMonitor(), and FitAndFill().

minVtxNdf_

double BeamMonitor::minVtxNdf_

private

Definition at line 108 of file BeamMonitor.h.

Referenced by analyze(), and BeamMonitor().

minVtxWgt_

double BeamMonitor::minVtxWgt_

private

Definition at line 109 of file BeamMonitor.h.

Referenced by analyze(), and BeamMonitor().

monitorName_

std::string BeamMonitor::monitorName_

private

Definition at line 71 of file BeamMonitor.h.

Referenced by BeamMonitor(), and bookHistograms().

nextlumi_

int BeamMonitor::nextlumi_

private

Definition at line 98 of file BeamMonitor.h.

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

nFitElements_

int BeamMonitor::nFitElements_

private

Definition at line 102 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

nFits_

int BeamMonitor::nFits_

private

Definition at line 103 of file BeamMonitor.h.

Referenced by BeamMonitor(), and FitAndFill().

nthBSTrk_

unsigned int BeamMonitor::nthBSTrk_

private

Definition at line 101 of file BeamMonitor.h.

Referenced by FitAndFill(), and RestartFitting().

onlineDbService_

edm::Service<cond::service::OnlineDBOutputService> BeamMonitor::onlineDbService_

private

Definition at line 77 of file BeamMonitor.h.

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

onlineMode_

bool BeamMonitor::onlineMode_

private

Definition at line 86 of file BeamMonitor.h.

Referenced by analyze(), BeamMonitor(), beginLuminosityBlock(), endLuminosityBlock(), and FitAndFill().

phiBin_

const int BeamMonitor::phiBin_

private

Definition at line 64 of file BeamMonitor.h.

Referenced by bookHistograms().

phiMax_

const double BeamMonitor::phiMax_

private

Definition at line 66 of file BeamMonitor.h.

Referenced by bookHistograms().

phiMin_

const double BeamMonitor::phiMin_

private

Definition at line 65 of file BeamMonitor.h.

Referenced by bookHistograms().

preBS

reco::BeamSpot BeamMonitor::preBS

private

Definition at line 122 of file BeamMonitor.h.

Referenced by FitAndFill().

processed_

bool BeamMonitor::processed_

private

Definition at line 116 of file BeamMonitor.h.

Referenced by analyze(), BeamMonitor(), beginLuminosityBlock(), and FitAndFill().

pvResults

MonitorElement* BeamMonitor::pvResults

private

Definition at line 149 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

pvSrc_

edm::EDGetTokenT<reco::VertexCollection> BeamMonitor::pvSrc_

private

Definition at line 75 of file BeamMonitor.h.

Referenced by analyze(), and BeamMonitor().

recordName_

std::string BeamMonitor::recordName_

private

Definition at line 72 of file BeamMonitor.h.

Referenced by BeamMonitor(), and FitAndFill().

refBS

reco::BeamSpot BeamMonitor::refBS

private

Definition at line 121 of file BeamMonitor.h.

Referenced by analyze().

refBStime

std::time_t BeamMonitor::refBStime[2]

private

Definition at line 99 of file BeamMonitor.h.

Referenced by BeamMonitor(), beginLuminosityBlock(), endLuminosityBlock(), FitAndFill(), and RestartFitting().

refPVtime

std::time_t BeamMonitor::refPVtime[2]

private

Definition at line 100 of file BeamMonitor.h.

Referenced by BeamMonitor(), beginLuminosityBlock(), endLuminosityBlock(), FitAndFill(), and RestartFitting().

refTime

std::time_t BeamMonitor::refTime

private

Definition at line 177 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

reportSummary

MonitorElement* BeamMonitor::reportSummary

private

Definition at line 168 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

reportSummary_

Float_t BeamMonitor::reportSummary_

private

Definition at line 165 of file BeamMonitor.h.

Referenced by FitAndFill().

reportSummaryContents

MonitorElement* BeamMonitor::reportSummaryContents[3]

private

Definition at line 169 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

reportSummaryMap

MonitorElement* BeamMonitor::reportSummaryMap

private

Definition at line 170 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

resetFitNLumi_

int BeamMonitor::resetFitNLumi_

private

Definition at line 82 of file BeamMonitor.h.

Referenced by analyze(), BeamMonitor(), and FitAndFill().

resetHistos_

bool BeamMonitor::resetHistos_

private

Definition at line 111 of file BeamMonitor.h.

Referenced by FitAndFill(), and RestartFitting().

resetPVNLumi_

int BeamMonitor::resetPVNLumi_

private

Definition at line 83 of file BeamMonitor.h.

Referenced by BeamMonitor(), and FitAndFill().

StartAverage_

bool BeamMonitor::StartAverage_

private

Definition at line 112 of file BeamMonitor.h.

Referenced by analyze(), beginLuminosityBlock(), FitAndFill(), and RestartFitting().

startTime

std::time_t BeamMonitor::startTime

private

Definition at line 176 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

summaryContent_

Float_t BeamMonitor::summaryContent_[3]

private

Definition at line 167 of file BeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

summarySum_

Float_t BeamMonitor::summarySum_

private

Definition at line 166 of file BeamMonitor.h.

Referenced by FitAndFill().

theBeamFitter

std::unique_ptr<BeamFitter> BeamMonitor::theBeamFitter

private

Definition at line 89 of file BeamMonitor.h.

Referenced by analyze(), BeamMonitor(), FitAndFill(), and RestartFitting().

tmpTime

std::time_t BeamMonitor::tmpTime

private

Definition at line 175 of file BeamMonitor.h.

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

tracksLabel_

edm::EDGetTokenT<reco::TrackCollection> BeamMonitor::tracksLabel_

private

Definition at line 74 of file BeamMonitor.h.

Referenced by analyze(), and BeamMonitor().

vxBin_

const int BeamMonitor::vxBin_

private

Definition at line 60 of file BeamMonitor.h.

Referenced by bookHistograms().

vxMax_

const double BeamMonitor::vxMax_

private

Definition at line 62 of file BeamMonitor.h.

Referenced by bookHistograms().

vxMin_

const double BeamMonitor::vxMin_

private

Definition at line 61 of file BeamMonitor.h.

Referenced by bookHistograms().