#include <FakeBeamMonitor.h>

Inheritance diagram for FakeBeamMonitor:

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

	~FakeBeamMonitor () override

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

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

	DQMOneEDAnalyzer ()

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

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args...>
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector < edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

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)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::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	dqmBeginRun (edm::Run const &, edm::EventSetup const &) 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 edm::ProducerBase
ProducesCollector	producesCollector ()

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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

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

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

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

void	consumesMany (const TypeToGet &id)

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

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

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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...

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

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

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

void	resetItemsToGetFrom (BranchType iType)

Private Member Functions
void	FitAndFill (const edm::LuminosityBlock &lumiSeg, int &, int &, int &)

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

std::string	getGMTstring (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_

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

std::vector< MonitorElement * >	hs

int	intervalInSec_

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

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_

TRandom3 *	rndm_

bool	StartAverage_

std::time_t	startTime

Float_t	summaryContent_ [3]

Float_t	summarySum_

std::time_t	tmpTime

bool	useLockRecords_

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
Francesco Brivio - Milano-Bicocca : Same behaviour of BeamMonitor except the beamspot is randomly generated and not fitted from tracks/PVs

Definition at line 34 of file FakeBeamMonitor.h.

Constructor & Destructor Documentation

FakeBeamMonitor::FakeBeamMonitor ( const edm::ParameterSet & ps )

Definition at line 104 of file FakeBeamMonitor.cc.

References beginLumiOfBSFit_, beginLumiOfPVFit_, debug_, deltaSigCut_, endLumiOfBSFit_, endLumiOfPVFit_, fitNLumi_, fitPVNLumi_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), intervalInSec_, lastlumi_, maxZ_, min_Ntrks_, minNrVertices_, minVtxNdf_, minVtxWgt_, monitorName_, nextlumi_, nFits_, onlineMode_, processed_, recordName_, refBStime, refPVtime, resetFitNLumi_, resetPVNLumi_, rndm_, and useLockRecords_.

     : 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");
   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");
   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");
   useLockRecords_ = ps.getUntrackedParameter<bool>("useLockRecords");
 
   if (!monitorName_.empty())
     monitorName_ = monitorName_ + "/";
 
   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;
 
   rndm_ = new TRandom3(0);
 }

FakeBeamMonitor::~FakeBeamMonitor ( )

override

Definition at line 102 of file FakeBeamMonitor.cc.

102 { delete rndm_; };

FakeBeamMonitor::rndm_

TRandom3 * rndm_

Definition: FakeBeamMonitor.h:181

Member Function Documentation

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

overrideprotectedvirtual

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

Definition at line 655 of file FakeBeamMonitor.cc.

References countEvt_, countLumi_, edm::EventBase::luminosityBlock(), mapNPV, mapPVx, mapPVy, mapPVz, nextlumi_, onlineMode_, processed_, resetFitNLumi_, and StartAverage_.

                                                                            {
   const int nthlumi = iEvent.luminosityBlock();
 
   if (onlineMode_ && (nthlumi < nextlumi_)) {
     edm::LogInfo("FakeBeamMonitor") << "analyze::  Spilt event from previous lumi section!" << std::endl;
     return;
   }
   if (onlineMode_ && (nthlumi > nextlumi_)) {
     edm::LogInfo("FakeBeamMonitor") << "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) {
   for (int tmp_idx = 0; tmp_idx < 10; tmp_idx++) {
     //--- 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(tmp_idx);
     mapPVy[countLumi_].push_back(tmp_idx);
     mapPVz[countLumi_].push_back(tmp_idx);
 
     //      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;
 }

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

overrideprotected

Definition at line 525 of file FakeBeamMonitor.cc.

References beginLumiOfBSFit_, beginLumiOfPVFit_, edm::LuminosityBlockBase::beginTime(), countLumi_, DBloggerReturn_, debug_, FitAndFill(), edm::LuminosityBlock::getRun(), 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, edm::RunBase::run(), cond::persistency::Logger::start(), StartAverage_, and edm::Timestamp::value().

                                                                                                     {
   // start DB logger
   DBloggerReturn_ = 0;
   if (onlineDbService_.isAvailable()) {
     onlineDbService_->logger().start();
     onlineDbService_->logger().logInfo() << "FakeBeamMonitor::beginLuminosityBlock - LS: " << lumiSeg.luminosityBlock()
                                          << " - Run: " << lumiSeg.getRun().run();
   }
 
   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("FakeBeamMonitor") << " 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("FakeBeamMonitor") << " 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("FakeBeamMonitor") << "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("FakeBeamMonitor") << " 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("FakeBeamMonitor") << " beginLuminosityBlock::  Begin of Lumi: " << nthlumi << endl;
 }

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

overrideprotectedvirtual

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

Definition at line 213 of file FakeBeamMonitor.cc.

References cms::cuda::assert(), edm::RunBase::beginTime(), dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookFloat(), dqm::implementation::IBooker::bookProfile(), cutFlowTable, plotBeamSpotDB::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, hs, mps_fire::i, intervalInSec_, label, monitorName_, dqmiodumpmetadata::n, nFitElements_, AlCaHLTBitMon_ParallelJobs::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_, and vxMin_.

                                                                                                               {
   frun = iRun.run();
   ftimestamp = iRun.beginTime().value();
   tmpTime = ftimestamp >> 32;
   startTime = refTime = tmpTime;
   char eventTime[64];
   formatFitTime(eventTime, tmpTime);
   edm::LogInfo("FakeBeamMonitor") << "TimeOffset = " << eventTime << std::endl;
   TDatime da(eventTime);
   if (debug_) {
     edm::LogInfo("FakeBeamMonitor") << "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.clear();
   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("FakeBeamMonitor") << "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.);
   }
 }

void FakeBeamMonitor::dqmBeginRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

overrideprotectedvirtual

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

Definition at line 207 of file FakeBeamMonitor.cc.

References edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::lockRecords(), onlineDbService_, and useLockRecords_.

                                                                      {
   if (useLockRecords_ && onlineDbService_.isAvailable()) {
     onlineDbService_->lockRecords();
   }
 }

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

overrideprotectedvirtual

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

Definition at line 1582 of file FakeBeamMonitor.cc.

References debug_, edm::Service< T >::isAvailable(), mapBeginBSLS, mapBeginBSTime, mapBeginPVLS, mapBeginPVTime, mapLSBSTrkSize, mapLSCF, mapLSPVStoreSize, mapNPV, mapPVx, mapPVy, mapPVz, onlineDbService_, cond::service::PoolDBOutputService::releaseLocks(), and useLockRecords_.

                                                                        {
   if (debug_)
     edm::LogInfo("FakeBeamMonitor") << "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();
   if (useLockRecords_ && onlineDbService_.isAvailable()) {
     onlineDbService_->releaseLocks();
   }
 }

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

overrideprotected

Definition at line 794 of file FakeBeamMonitor.cc.

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().

                                                                                                  {
   int nthlumi = lumiSeg.id().luminosityBlock();
   edm::LogInfo("FakeBeamMonitor") << "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() << "FakeBeamMonitor::endLuminosityBlock";
     onlineDbService_->logger().end(DBloggerReturn_);
   }
 }

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

private

Definition at line 813 of file FakeBeamMonitor.cc.

Referenced by beginLuminosityBlock().

                                                                                                            {
   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("FakeBeamMonitor") << "FitAndFill::  Lumi of the current fit: " << currentlumi << endl;
   lastlumi = currentlumi;
   endLumiOfBSFit_ = currentlumi;
   endLumiOfPVFit_ = currentlumi;
 
   //---------Fix for Runninv average-------------
   mapLSPVStoreSize[countLumi_] = 10;  //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_);
   DipPVInfo_.push_back(rndm_->Gaus(1000., 100.));  // Events used
   DipPVInfo_.push_back(rndm_->Gaus(100., 10.));    // Mean PV
   DipPVInfo_.push_back(rndm_->Gaus(10., 5.));      // Mean PV err
   DipPVInfo_.push_back(rndm_->Gaus(10., 5.));      // Rms PV
   DipPVInfo_.push_back(rndm_->Gaus(5., 3.));       // Rms PV err
   DipPVInfo_.push_back(rndm_->Gaus(100., 10.));    // Max PVs
   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("FakeBeamMonitor") << "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("FakeBeamMonitor") << "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("FakeBeamMonitor") << "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("FakeBeamMonitor") << "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("FakeBeamMonitor") << "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("FakeBeamMonitor") << " 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("FakeBeamMonitor") << " The Current Recored for this fit is  =" << nthBSTrk_ << endl;
 
   //  if (countFitting)
   //    edm::LogInfo("FakeBeamMonitor") << "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("FakeBeamMonitor") << "FitAndFill:: [DebugTime] refBStime[0] = " << refBStime[0]
                                   << "; address =  " << &refBStime[0] << std::endl;
   edm::LogInfo("FakeBeamMonitor") << "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();
   std::pair<int, int> fitLS(beginLumiOfBSFit_, endLumiOfBSFit_);
   //    edm::LogInfo("BeamMonitor") << "FitAndFill::  [BeamFitter] Do BeamSpot Fit for LS = " << fitLS.first << " to "
   //                                << fitLS.second << std::endl;
   edm::LogInfo("FakeBeamMonitor") << "FitAndFill::  [FakeBeamMonitor] 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();
 
   //Create fake BS here
   reco::BeamSpot::CovarianceMatrix matrix;
   for (int j = 0; j < 7; ++j) {
     for (int k = j; k < 7; ++k) {
       matrix(j, k) = 0;
     }
   }
 
   // random values for fake BeamSpot
   float tmp_BSx = rndm_->Gaus(0.1, 0.1);            // [cm]
   float tmp_BSy = rndm_->Gaus(0.1, 0.1);            // [cm]
   float tmp_BSz = rndm_->Gaus(0.1, 0.1);            // [cm]
   float tmp_BSwidthX = rndm_->Gaus(0.001, 0.0005);  // [cm]
   float tmp_BSwidthY = rndm_->Gaus(0.001, 0.0005);  // [cm]
   float tmp_BSwidthZ = rndm_->Gaus(3.5, 0.5);       // [cm]
 
   reco::BeamSpot bs(reco::BeamSpot::Point(tmp_BSx, tmp_BSy, tmp_BSz),
                     tmp_BSwidthZ,
                     0,
                     0,
                     tmp_BSwidthX,
                     matrix,
                     reco::BeamSpot::Tracker);
   bs.setBeamWidthY(tmp_BSwidthY);
 
   if (bs.type() > 0)  // with good beamwidth fit
     preBS = bs;       // cache good fit results
 
   edm::LogInfo("FakeBeamMonitor") << "\n RESULTS OF DEFAULT FIT:" << endl;
   edm::LogInfo("FakeBeamMonitor") << bs << endl;
   edm::LogInfo("FakeBeamMonitor") << "[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(beginLumiOfBSFit_, endLumiOfBSFit_);  //= 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;
   BSOnline.SetLastAnalyzedLumi(LSRange.second);
   BSOnline.SetLastAnalyzedRun(frun);
   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(50);
   BSOnline.SetNumPVs(10);
   BSOnline.SetUsedEvents((int)DipPVInfo_[0]);
   BSOnline.SetMeanPV(DipPVInfo_[1]);
   BSOnline.SetMeanErrorPV(DipPVInfo_[2]);
   BSOnline.SetRmsPV(DipPVInfo_[3]);
   BSOnline.SetRmsErrorPV(DipPVInfo_[4]);
   BSOnline.SetMaxPVs((int)DipPVInfo_[5]);
   auto creationTime =
       std::chrono::duration_cast<std::chrono::microseconds>(std::chrono::system_clock::now().time_since_epoch()).count();
   BSOnline.SetCreationTime(creationTime);
 
   // use fake timestamps from 1970.01.01 00:00:00 to 1970.01.01 00:00:01 GMT
   std::pair<time_t, time_t> timeForDIP = std::make_pair(0, 1);
   BSOnline.SetStartTimeStamp(timeForDIP.first);
   BSOnline.SetStartTime(getGMTstring(timeForDIP.first));
   BSOnline.SetEndTimeStamp(timeForDIP.second);
   BSOnline.SetEndTime(getGMTstring(timeForDIP.second));
 
   edm::LogInfo("FakeBeamMonitor") << "FitAndFill::[PayloadCreation] BeamSpotOnline object created: \n" << std::endl;
   edm::LogInfo("FakeBeamMonitor") << BSOnline << std::endl;
 
   // Create the payload for BeamSpotOnlineObjects object
   if (onlineDbService_.isAvailable()) {
     edm::LogInfo("FakeBeamMonitor") << "FitAndFill::[PayloadCreation] onlineDbService available \n" << std::endl;
     onlineDbService_->logger().logInfo() << "FakeBeamMonitor::FitAndFill - Lumi of the current fit: " << currentlumi;
     onlineDbService_->logger().logInfo() << "FakeBeamMonitor::FitAndFill - Do PV Fitting for LS = " << beginLumiOfPVFit_
                                          << " to " << endLumiOfPVFit_;
     onlineDbService_->logger().logInfo() << "FakeBeamMonitor::FitAndFill - [BeamFitter] Do BeamSpot Fit for LS = "
                                          << LSRange.first << " to " << LSRange.second;
     onlineDbService_->logger().logInfo() << "FakeBeamMonitor::FitAndFill - [FakeBeamMonitor] Do BeamSpot Fit for LS = "
                                          << beginLumiOfBSFit_ << " to " << endLumiOfBSFit_;
     onlineDbService_->logger().logInfo() << "FakeBeamMonitor - RESULTS OF DEFAULT FIT:";
     onlineDbService_->logger().logInfo() << "\n" << bs;
     onlineDbService_->logger().logInfo()
         << "FakeBeamMonitor::FitAndFill - [PayloadCreation] BeamSpotOnline object created:";
     onlineDbService_->logger().logInfo() << "\n" << BSOnline;
     onlineDbService_->logger().logInfo() << "FakeBeamMonitor - Additional parameters for DIP:";
     onlineDbService_->logger().logInfo() << "Events used in the fit: " << BSOnline.GetUsedEvents();
     onlineDbService_->logger().logInfo() << "Mean PV               : " << BSOnline.GetMeanPV();
     onlineDbService_->logger().logInfo() << "Mean PV Error         : " << BSOnline.GetMeanErrorPV();
     onlineDbService_->logger().logInfo() << "Rms PV                : " << BSOnline.GetRmsPV();
     onlineDbService_->logger().logInfo() << "Rms PV Error          : " << BSOnline.GetRmsErrorPV();
     onlineDbService_->logger().logInfo() << "Max PVs               : " << BSOnline.GetMaxPVs();
     onlineDbService_->logger().logInfo() << "StartTime             : " << BSOnline.GetStartTime();
     onlineDbService_->logger().logInfo() << "StartTimeStamp        : " << BSOnline.GetStartTimeStamp();
     onlineDbService_->logger().logInfo() << "EndTime               : " << BSOnline.GetEndTime();
     onlineDbService_->logger().logInfo() << "EndTimeStamp          : " << BSOnline.GetEndTimeStamp();
     onlineDbService_->logger().logInfo() << "FakeBeamMonitor::FitAndFill - [PayloadCreation] onlineDbService available";
     onlineDbService_->logger().logInfo() << "FakeBeamMonitor::FitAndFill - [PayloadCreation] SetCreationTime: "
                                          << creationTime << " [epoch in microseconds]";
     try {
       onlineDbService_->writeIOVForNextLumisection(BSOnline, recordName_);
       onlineDbService_->logger().logInfo()
           << "FakeBeamMonitor::FitAndFill - [PayloadCreation] writeIOVForNextLumisection executed correctly";
     } catch (const std::exception& e) {
       onlineDbService_->logger().logError() << "FakeBeamMonitor - Error writing record: " << recordName_
                                             << " for Run: " << frun << " - Lumi: " << LSRange.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("FakeBeamMonitor") << "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;
   //
   //      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;
   //
   //    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("FakeBeamMonitor") << "FitAndFill:: The flag is ON for running average Beam Spot fit" << endl;
     StartAverage_ = true;
     firstAverageFit_++;
     resetHistos_ = true;
     nthBSTrk_ = 0;
     beginLumiOfBSFit_ = 0;
     refBStime[0] = 0;
   }
 }

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

private

Definition at line 44 of file FakeBeamMonitor.cc.

References b, CEST, and plotBeamSpotDB::ptm.

Referenced by bookHistograms(), and scrollTH1().

                                                              {
   //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
 }

std::string FakeBeamMonitor::getGMTstring ( const std::time_t & )

private

Definition at line 92 of file FakeBeamMonitor.cc.

References AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by FitAndFill().

                                                                    {
   char buff[32];
   std::strftime(buff, sizeof(buff), "%Y.%m.%d %H:%M:%S GMT", gmtime(&timeToConvert));
   std::string timeStr(buff);
   return timeStr;
 }

void FakeBeamMonitor::RestartFitting ( )

private

Definition at line 1544 of file FakeBeamMonitor.cc.

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_, and StartAverage_.

                                      {  //
   if (debug_)
     edm::LogInfo("FakeBeamMonitor")
         << " 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;
 }

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

private

Definition at line 1603 of file FakeBeamMonitor.cc.

References newFWLiteAna::bin, buffTime, plotBeamSpotDB::da, formatFitTime(), lastNZbin, and hgcalPerformanceValidation::val.

Referenced by FitAndFill().

                                                   {
   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));
   }
 }

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

private

Definition at line 1623 of file FakeBeamMonitor.cc.

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

Referenced by FitAndFill().

                                                                    {
   bool scroll_ = false;
   if (tmpTime_ - refTime_ >= intervalInSec_) {
     scroll_ = true;
     edm::LogInfo("FakeBeamMonitor") << "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("FakeBeamMonitor") << "testScroll::  Last non zero bin = " << lastNZbin << std::endl;
     if (tmpTime_ - refTime_ >= intervalInSec_ + lastNZbin) {
       edm::LogInfo("FakeBeamMonitor") << "testScroll::  Time difference too large since last readout" << std::endl;
       lastNZbin = 0;
       refTime_ = tmpTime_ - buffTime;
     } else {
       edm::LogInfo("FakeBeamMonitor") << "testScroll::  Offset to last record" << std::endl;
       int offset = ((lastNZbin > buffTime) ? (lastNZbin - buffTime) : (lastNZbin - 1));
       refTime_ += offset;
     }
   }
   return scroll_;
 }

Member Data Documentation

int FakeBeamMonitor::beginLumiOfBSFit_

private

Definition at line 91 of file FakeBeamMonitor.h.

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

int FakeBeamMonitor::beginLumiOfPVFit_

private

Definition at line 93 of file FakeBeamMonitor.h.

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

int FakeBeamMonitor::countEvt_

private

Definition at line 89 of file FakeBeamMonitor.h.

Referenced by analyze().

int FakeBeamMonitor::countGapLumi_

private

Definition at line 112 of file FakeBeamMonitor.h.

Referenced by RestartFitting().

int FakeBeamMonitor::countLumi_

private

Definition at line 90 of file FakeBeamMonitor.h.

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

MonitorElement* FakeBeamMonitor::cutFlowTable

private

Definition at line 169 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

int FakeBeamMonitor::DBloggerReturn_

private

Definition at line 79 of file FakeBeamMonitor.h.

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

bool FakeBeamMonitor::debug_

private

Definition at line 86 of file FakeBeamMonitor.h.

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

double FakeBeamMonitor::deltaSigCut_

private

Definition at line 102 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor().

const int FakeBeamMonitor::dxBin_

private

Definition at line 61 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

const double FakeBeamMonitor::dxMax_

private

Definition at line 63 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

const double FakeBeamMonitor::dxMin_

private

Definition at line 62 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

const int FakeBeamMonitor::dzBin_

private

Definition at line 73 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

const double FakeBeamMonitor::dzMax_

private

Definition at line 75 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

const double FakeBeamMonitor::dzMin_

private

Definition at line 74 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

int FakeBeamMonitor::endLumiOfBSFit_

private

Definition at line 92 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

int FakeBeamMonitor::endLumiOfPVFit_

private

Definition at line 94 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

int FakeBeamMonitor::firstAverageFit_

private

Definition at line 111 of file FakeBeamMonitor.h.

Referenced by FitAndFill().

int FakeBeamMonitor::fitNLumi_

private

Definition at line 81 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

int FakeBeamMonitor::fitPVNLumi_

private

Definition at line 82 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

MonitorElement* FakeBeamMonitor::fitResults

private

Definition at line 133 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

int FakeBeamMonitor::frun

private

Definition at line 178 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

edm::TimeValue_t FakeBeamMonitor::ftimestamp

private

Definition at line 176 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

MonitorElement* FakeBeamMonitor::h_d0_phi0

private

Definition at line 126 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::h_nTrk_lumi

private

Definition at line 123 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::h_nVtx

private

Definition at line 140 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

MonitorElement* FakeBeamMonitor::h_nVtx_lumi

private

Definition at line 124 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::h_nVtx_lumi_all

private

Definition at line 125 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

MonitorElement* FakeBeamMonitor::h_nVtx_st

private

Definition at line 141 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::h_PVx[2]

private

Definition at line 142 of file FakeBeamMonitor.h.

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

MonitorElement* FakeBeamMonitor::h_PVxz

private

Definition at line 145 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

MonitorElement* FakeBeamMonitor::h_PVy[2]

private

Definition at line 143 of file FakeBeamMonitor.h.

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

MonitorElement* FakeBeamMonitor::h_PVyz

private

Definition at line 146 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

MonitorElement* FakeBeamMonitor::h_PVz[2]

private

Definition at line 144 of file FakeBeamMonitor.h.

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

MonitorElement* FakeBeamMonitor::h_sigmaX0

private

Definition at line 137 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::h_sigmaY0

private

Definition at line 138 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::h_sigmaZ0

private

Definition at line 139 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::h_trk_z0

private

Definition at line 127 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::h_trkPt

private

Definition at line 131 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

MonitorElement* FakeBeamMonitor::h_trkVz

private

Definition at line 132 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

MonitorElement* FakeBeamMonitor::h_vx_dz

private

Definition at line 129 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::h_vx_vy

private

Definition at line 128 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::h_vy_dz

private

Definition at line 130 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::h_x0

private

Definition at line 134 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::h_y0

private

Definition at line 135 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::h_z0

private

Definition at line 136 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

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

private

Definition at line 148 of file FakeBeamMonitor.h.

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

int FakeBeamMonitor::intervalInSec_

private

Definition at line 85 of file FakeBeamMonitor.h.

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

int FakeBeamMonitor::lastlumi_

private

Definition at line 95 of file FakeBeamMonitor.h.

Referenced by beginLuminosityBlock(), and FakeBeamMonitor().

int FakeBeamMonitor::lastNZbin

private

Definition at line 179 of file FakeBeamMonitor.h.

Referenced by scrollTH1(), and testScroll().

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

private

Definition at line 154 of file FakeBeamMonitor.h.

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

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

private

Definition at line 155 of file FakeBeamMonitor.h.

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

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

private

Definition at line 154 of file FakeBeamMonitor.h.

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

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

private

Definition at line 155 of file FakeBeamMonitor.h.

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

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

private

Definition at line 157 of file FakeBeamMonitor.h.

Referenced by dqmEndRun(), and RestartFitting().

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

private

Definition at line 160 of file FakeBeamMonitor.h.

Referenced by dqmEndRun(), and RestartFitting().

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

private

Definition at line 158 of file FakeBeamMonitor.h.

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

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

private

Definition at line 152 of file FakeBeamMonitor.h.

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

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

private

Definition at line 151 of file FakeBeamMonitor.h.

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

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

private

Definition at line 151 of file FakeBeamMonitor.h.

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

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

private

Definition at line 151 of file FakeBeamMonitor.h.

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

double FakeBeamMonitor::maxZ_

private

Definition at line 104 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor().

unsigned int FakeBeamMonitor::min_Ntrks_

private

Definition at line 103 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

unsigned int FakeBeamMonitor::minNrVertices_

private

Definition at line 105 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

double FakeBeamMonitor::minVtxNdf_

private

Definition at line 106 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor().

double FakeBeamMonitor::minVtxWgt_

private

Definition at line 107 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor().

std::string FakeBeamMonitor::monitorName_

private

Definition at line 76 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FakeBeamMonitor().

int FakeBeamMonitor::nextlumi_

private

Definition at line 96 of file FakeBeamMonitor.h.

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

int FakeBeamMonitor::nFitElements_

private

Definition at line 100 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

int FakeBeamMonitor::nFits_

private

Definition at line 101 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

unsigned int FakeBeamMonitor::nthBSTrk_

private

Definition at line 99 of file FakeBeamMonitor.h.

Referenced by FitAndFill(), and RestartFitting().

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

private

Definition at line 78 of file FakeBeamMonitor.h.

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

bool FakeBeamMonitor::onlineMode_

private

Definition at line 87 of file FakeBeamMonitor.h.

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

const int FakeBeamMonitor::phiBin_

private

Definition at line 69 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

const double FakeBeamMonitor::phiMax_

private

Definition at line 71 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

const double FakeBeamMonitor::phiMin_

private

Definition at line 70 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

reco::BeamSpot FakeBeamMonitor::preBS

private

Definition at line 120 of file FakeBeamMonitor.h.

Referenced by FitAndFill().

bool FakeBeamMonitor::processed_

private

Definition at line 114 of file FakeBeamMonitor.h.

Referenced by analyze(), beginLuminosityBlock(), and FakeBeamMonitor().

MonitorElement* FakeBeamMonitor::pvResults

private

Definition at line 147 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

std::string FakeBeamMonitor::recordName_

private

Definition at line 77 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

reco::BeamSpot FakeBeamMonitor::refBS

private

Definition at line 119 of file FakeBeamMonitor.h.

std::time_t FakeBeamMonitor::refBStime[2]

private

Definition at line 97 of file FakeBeamMonitor.h.

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

std::time_t FakeBeamMonitor::refPVtime[2]

private

Definition at line 98 of file FakeBeamMonitor.h.

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

std::time_t FakeBeamMonitor::refTime

private

Definition at line 175 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::reportSummary

private

Definition at line 166 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

Float_t FakeBeamMonitor::reportSummary_

private

Definition at line 163 of file FakeBeamMonitor.h.

Referenced by FitAndFill().

MonitorElement* FakeBeamMonitor::reportSummaryContents[3]

private

Definition at line 167 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

MonitorElement* FakeBeamMonitor::reportSummaryMap

private

Definition at line 168 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

int FakeBeamMonitor::resetFitNLumi_

private

Definition at line 83 of file FakeBeamMonitor.h.

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

bool FakeBeamMonitor::resetHistos_

private

Definition at line 109 of file FakeBeamMonitor.h.

Referenced by FitAndFill(), and RestartFitting().

int FakeBeamMonitor::resetPVNLumi_

private

Definition at line 84 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

TRandom3* FakeBeamMonitor::rndm_

private

Definition at line 181 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

bool FakeBeamMonitor::StartAverage_

private

Definition at line 110 of file FakeBeamMonitor.h.

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

std::time_t FakeBeamMonitor::startTime

private

Definition at line 174 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

Float_t FakeBeamMonitor::summaryContent_[3]

private

Definition at line 165 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

Float_t FakeBeamMonitor::summarySum_

private

Definition at line 164 of file FakeBeamMonitor.h.

Referenced by FitAndFill().

std::time_t FakeBeamMonitor::tmpTime

private

Definition at line 173 of file FakeBeamMonitor.h.

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

bool FakeBeamMonitor::useLockRecords_

private

Definition at line 115 of file FakeBeamMonitor.h.

Referenced by dqmBeginRun(), dqmEndRun(), and FakeBeamMonitor().

const int FakeBeamMonitor::vxBin_

private

Definition at line 65 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

const double FakeBeamMonitor::vxMax_

private

Definition at line 67 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

const double FakeBeamMonitor::vxMin_

private

Definition at line 66 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation