#include <BeamMonitorBx.h>

Inheritance diagram for BeamMonitorBx:

Public Types
typedef std::map< BxNum, reco::BeamSpot >	BeamSpotMapBx

typedef int	BxNum

typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

	~BeamMonitorBx () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

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

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

void	beginJob () override

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

void	beginRun (const edm::Run &r, const edm::EventSetup &c) override

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

void	endRun (const edm::Run &r, const edm::EventSetup &c) override

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

template<BranchType B>
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)

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

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

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

Private Member Functions
void	BookTables (int, std::map< std::string, std::string > &, std::string)

void	BookTrendHistos (bool, int, std::map< std::string, std::string > &, std::string, const TString &, const TString &)

void	FillTables (int, int, std::map< std::string, std::string > &, reco::BeamSpot &, std::string)

void	FillTrendHistos (int, int, std::map< std::string, std::string > &, reco::BeamSpot &, const TString &)

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

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

void	weight (BeamSpotMapBx &, const BeamSpotMapBx &)

void	weight (double &mean, double &meanError, const double &val, const double &valError)

Private Attributes
int	beginLumiOfBSFit_

edm::InputTag	bsSrc_

unsigned int	countBx_

int	countEvt_

int	countGoodFit_

int	countLumi_

DQMStore *	dbe_

bool	debug_

int	endLumiOfBSFit_

BeamSpotMapBx	fbspotMap

int	firstlumi_

int	fitNLumi_

edm::TimeValue_t	ftimestamp

std::map< TString, MonitorElement * >	hs

std::map< TString, MonitorElement * >	hst

int	lastlumi_

std::string	monitorName_

int	nextlumi_

edm::ParameterSet	parameters_

bool	processed_

std::time_t	refBStime [2]

std::time_t	refTime

int	resetFitNLumi_

bool	resetHistos_

std::time_t	startTime

BeamFitter *	theBeamFitter

std::time_t	tmpTime

std::map< std::string, std::string >	varMap

std::map< std::string, std::string >	varMap1

Additional Inherited Members
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Detailed Description

Author
Geng-yuan Jeng/UC Riverside Francisco Yumiceva/FNAL

Definition at line 25 of file BeamMonitorBx.h.

Member Typedef Documentation

◆ BeamSpotMapBx

typedef std::map<BxNum, reco::BeamSpot> BeamMonitorBx::BeamSpotMapBx

Definition at line 34 of file BeamMonitorBx.h.

◆ BxNum

typedef int BeamMonitorBx::BxNum

Definition at line 33 of file BeamMonitorBx.h.

◆ DQMStore

typedef dqm::legacy::DQMStore BeamMonitorBx::DQMStore

Definition at line 28 of file BeamMonitorBx.h.

◆ MonitorElement

typedef dqm::legacy::MonitorElement BeamMonitorBx::MonitorElement

Definition at line 27 of file BeamMonitorBx.h.

Constructor & Destructor Documentation

◆ BeamMonitorBx()

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

Definition at line 54 of file BeamMonitorBx.cc.

                                                    : countBx_(0), countEvt_(0), countLumi_(0), resetHistos_(false) {
   parameters_ = ps;
   monitorName_ = parameters_.getUntrackedParameter<string>("monitorName", "YourSubsystemName");
   bsSrc_ = parameters_.getUntrackedParameter<InputTag>("beamSpot");
   fitNLumi_ = parameters_.getUntrackedParameter<int>("fitEveryNLumi", -1);
   resetFitNLumi_ = parameters_.getUntrackedParameter<int>("resetEveryNLumi", -1);
  
   dbe_ = Service<DQMStore>().operator->();
  
   if (!monitorName_.empty())
     monitorName_ = monitorName_ + "/";
  
   theBeamFitter = new BeamFitter(parameters_, consumesCollector());
   theBeamFitter->resetTrkVector();
   theBeamFitter->resetLSRange();
   theBeamFitter->resetRefTime();
   theBeamFitter->resetPVFitter();
  
   if (fitNLumi_ <= 0)
     fitNLumi_ = 1;
   beginLumiOfBSFit_ = endLumiOfBSFit_ = 0;
   refBStime[0] = refBStime[1] = 0;
   lastlumi_ = 0;
   nextlumi_ = 0;
   firstlumi_ = 0;
   processed_ = false;
   countGoodFit_ = 0;
 }

References AlcaBeamMonitor_cfi::BeamFitter, beginLumiOfBSFit_, bsSrc_, edm::EDConsumerBase::consumesCollector(), countGoodFit_, dbe_, endLumiOfBSFit_, firstlumi_, fitNLumi_, edm::ParameterSet::getUntrackedParameter(), lastlumi_, monitorName_, nextlumi_, parameters_, processed_, refBStime, resetFitNLumi_, BeamFitter::resetLSRange(), BeamFitter::resetPVFitter(), BeamFitter::resetRefTime(), BeamFitter::resetTrkVector(), and theBeamFitter.

◆ ~BeamMonitorBx()

BeamMonitorBx::~BeamMonitorBx ( )

override

Definition at line 83 of file BeamMonitorBx.cc.

83 { delete theBeamFitter; }

References theBeamFitter.

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Implements edm::EDAnalyzer.

Definition at line 158 of file BeamMonitorBx.cc.

                                                                          {
   bool readEvent_ = true;
   const int nthlumi = iEvent.luminosityBlock();
   if (nthlumi < nextlumi_) {
     readEvent_ = false;
     edm::LogWarning("BX|BeamMonitorBx") << "Spilt event from previous lumi section!" << endl;
   }
   if (nthlumi > nextlumi_) {
     readEvent_ = false;
     edm::LogWarning("BX|BeamMonitorBx") << "Spilt event from next lumi section!!!" << endl;
   }
  
   if (readEvent_) {
     countEvt_++;
     theBeamFitter->readEvent(iEvent);
     processed_ = true;
   }
 }

References countEvt_, iEvent, nextlumi_, processed_, BeamFitter::readEvent(), and theBeamFitter.

◆ beginJob()

void BeamMonitorBx::beginJob ( void )

overrideprotectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 86 of file BeamMonitorBx.cc.

                              {
   varMap["x0_bx"] = "X_{0}";
   varMap["y0_bx"] = "Y_{0}";
   varMap["z0_bx"] = "Z_{0}";
   varMap["sigmaX_bx"] = "#sigma_{X}";
   varMap["sigmaY_bx"] = "#sigma_{Y}";
   varMap["sigmaZ_bx"] = "#sigma_{Z}";
  
   varMap1["x"] = "X_{0} [cm]";
   varMap1["y"] = "Y_{0} [cm]";
   varMap1["z"] = "Z_{0} [cm]";
   varMap1["sigmaX"] = "#sigma_{X} [cm]";
   varMap1["sigmaY"] = "#sigma_{Y} [cm]";
   varMap1["sigmaZ"] = "#sigma_{Z} [cm]";
  
   // create and cd into new folder
   dbe_->setCurrentFolder(monitorName_ + "FitBx");
   // Results of good fit:
   BookTables(1, varMap, "");
   //if (resetFitNLumi_ > 0) BookTables(1,varMap,"all");
  
   // create and cd into new folders
   for (std::map<std::string, std::string>::const_iterator varName = varMap1.begin(); varName != varMap1.end();
        ++varName) {
     string subDir_ = "FitBx";
     subDir_ += "/";
     subDir_ += "All_";
     subDir_ += varName->first;
     dbe_->setCurrentFolder(monitorName_ + subDir_);
   }
   dbe_->setCurrentFolder(monitorName_ + "FitBx/All_nPVs");
 }

References BookTables(), dbe_, monitorName_, dqm::implementation::DQMStore::setCurrentFolder(), varMap, and varMap1.

◆ beginLuminosityBlock()

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

overrideprotectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 127 of file BeamMonitorBx.cc.

                                                                                                   {
   int nthlumi = lumiSeg.luminosityBlock();
   const edm::TimeValue_t fbegintimestamp = lumiSeg.beginTime().value();
   const std::time_t ftmptime = fbegintimestamp >> 32;
  
   if (countLumi_ == 0) {
     beginLumiOfBSFit_ = nthlumi;
     refBStime[0] = ftmptime;
   }
   if (beginLumiOfBSFit_ == 0)
     beginLumiOfBSFit_ = nextlumi_;
  
   if (nthlumi < nextlumi_)
     return;
  
   if (nthlumi > nextlumi_) {
     if (countLumi_ != 0 && processed_) {
       FitAndFill(lumiSeg, lastlumi_, nextlumi_, nthlumi);
     }
     nextlumi_ = nthlumi;
     edm::LogInfo("LS|BX|BeamMonitorBx") << "Next Lumi to Fit: " << nextlumi_ << endl;
     if (refBStime[0] == 0)
       refBStime[0] = ftmptime;
   }
   countLumi_++;
   if (processed_)
     processed_ = false;
   edm::LogInfo("LS|BX|BeamMonitorBx") << "Begin of Lumi: " << nthlumi << endl;
 }

References beginLumiOfBSFit_, edm::LuminosityBlockBase::beginTime(), countLumi_, FitAndFill(), lastlumi_, edm::LuminosityBlockBase::luminosityBlock(), nextlumi_, processed_, refBStime, and edm::Timestamp::value().

◆ beginRun()

void BeamMonitorBx::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

overrideprotectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 120 of file BeamMonitorBx.cc.

                                                                        {
   ftimestamp = r.beginTime().value();
   tmpTime = ftimestamp >> 32;
   startTime = refTime = tmpTime;
 }

References ftimestamp, alignCSCRings::r, refTime, startTime, and tmpTime.

◆ BookTables()

void BeamMonitorBx::BookTables	(	int	,
		std::map< std::string, std::string > &	,
		std::string
	)

private

Definition at line 190 of file BeamMonitorBx.cc.

                                                                                  {
   // to rebin histograms when number of bx increases
   dbe_->cd(monitorName_ + "FitBx");
  
   for (std::map<std::string, std::string>::const_iterator varName = vMap.begin(); varName != vMap.end(); ++varName) {
     string tmpName = varName->first;
     if (!suffix_.empty()) {
       tmpName += "_";
       tmpName += suffix_;
     }
  
     hs[tmpName] = dbe_->book2D(tmpName, varName->second, 3, 0, 3, nBx, 0, nBx);
     hs[tmpName]->setBinLabel(1, "bx", 1);
     hs[tmpName]->setBinLabel(2, varName->second, 1);
     hs[tmpName]->setBinLabel(3, "#Delta " + varName->second, 1);
     for (int i = 0; i < nBx; i++) {
       hs[tmpName]->setBinLabel(i + 1, " ", 2);
     }
     hs[tmpName]->getTH1()->SetOption("text");
     hs[tmpName]->Reset();
   }
 }

References dqm::implementation::IBooker::book2D(), dqm::implementation::DQMStore::cd(), dbe_, hs, mps_fire::i, and monitorName_.

Referenced by beginJob(), and FitAndFill().

◆ BookTrendHistos()

void BeamMonitorBx::BookTrendHistos	(	bool	,
		int	,
		std::map< std::string, std::string > &	,
		std::string	,
		const TString &	,
		const TString &
	)

private

Definition at line 214 of file BeamMonitorBx.cc.

                                                                                                                      {
   int nType_ = 2;
   if (plotPV)
     nType_ = 4;
   std::ostringstream ss;
   std::ostringstream ss1;
   ss << setfill('0') << setw(5) << nBx;
   ss1 << nBx;
  
   for (int i = 0; i < nType_; i++) {
     for (std::map<std::string, std::string>::const_iterator varName = vMap.begin(); varName != vMap.end(); ++varName) {
       string tmpDir_ = subDir_ + "/All_" + varName->first;
       dbe_->cd(monitorName_ + tmpDir_);
       TString histTitle(varName->first);
       TString tmpName;
       if (prefix_ != "")
         tmpName = prefix_ + "_" + varName->first;
       if (suffix_ != "")
         tmpName = tmpName + "_" + suffix_;
       tmpName = tmpName + "_" + ss.str();
  
       TString histName(tmpName);
       string ytitle(varName->second);
       string xtitle("");
       string options("E1");
       bool createHisto = true;
       switch (i) {
         case 1:  // BS vs time
           histName.Insert(histName.Index("_bx_", 4), "_time");
           xtitle = "Time [UTC]  [Bx# " + ss1.str() + "]";
           if (ytitle.find("sigma") == string::npos)
             histTitle += " coordinate of beam spot vs time (Fit)";
           else
             histTitle = histTitle.Insert(5, " ") + " of beam spot vs time (Fit)";
           break;
         case 2:  // PV +/- sigmaPV vs lumi
           if (ytitle.find("sigma") == string::npos) {
             histName.Insert(0, "PV");
             histName.Insert(histName.Index("_bx_", 4), "_lumi");
             histTitle.Insert(0, "Avg. ");
             histTitle += " position of primary vtx vs lumi";
             xtitle = "Lumisection  [Bx# " + ss1.str() + "]";
             ytitle.insert(0, "PV");
             ytitle += " #pm #sigma_{PV";
             ytitle += varName->first;
             ytitle += "} (cm)";
           } else
             createHisto = false;
           break;
         case 3:  // PV +/- sigmaPV vs time
           if (ytitle.find("sigma") == string::npos) {
             histName.Insert(0, "PV");
             histName.Insert(histName.Index("_bx_", 4), "_time");
             histTitle.Insert(0, "Avg. ");
             histTitle += " position of primary vtx vs time";
             xtitle = "Time [UTC]  [Bx# " + ss1.str() + "]";
             ytitle.insert(0, "PV");
             ytitle += " #pm #sigma_{PV";
             ytitle += varName->first;
             ytitle += "} (cm)";
           } else
             createHisto = false;
           break;
         default:  // BS vs lumi
           histName.Insert(histName.Index("_bx_", 4), "_lumi");
           xtitle = "Lumisection  [Bx# " + ss1.str() + "]";
           if (ytitle.find("sigma") == string::npos)
             histTitle += " coordinate of beam spot vs lumi (Fit)";
           else
             histTitle = histTitle.Insert(5, " ") + " of beam spot vs lumi (Fit)";
           break;
       }
       // check if already exist
       if (dbe_->get(monitorName_ + tmpDir_ + "/" + string(histName.Data())))
         continue;
  
       if (createHisto) {
         edm::LogInfo("BX|BeamMonitorBx") << "histName = " << histName << "; histTitle = " << histTitle << std::endl;
         hst[histName] = dbe_->book1D(histName, histTitle, 40, 0.5, 40.5);
         hst[histName]->getTH1()->SetCanExtend(TH1::kAllAxes);
         hst[histName]->setAxisTitle(xtitle, 1);
         hst[histName]->setAxisTitle(ytitle, 2);
         hst[histName]->getTH1()->SetOption("E1");
         if (histName.Contains("time")) {
           hst[histName]->getTH1()->SetBins(3600, 0.5, 3600 + 0.5);
           hst[histName]->setAxisTimeDisplay(1);
           hst[histName]->setAxisTimeFormat("%H:%M:%S", 1);
  
           char eventTime[64];
           formatFitTime(eventTime, startTime);
           TDatime da(eventTime);
           if (debug_) {
             edm::LogInfo("BX|BeamMonitorBx") << "TimeOffset = ";
             da.Print();
           }
           hst[histName]->getTH1()->GetXaxis()->SetTimeOffset(da.Convert(kTRUE));
         }
       }
     }  //End of variable loop
   }    // End of type loop (lumi, time)
  
   // num of PVs(#Bx) per LS
   dbe_->cd(monitorName_ + subDir_ + "/All_nPVs");
   TString histName = "Trending_nPVs_lumi_bx_" + ss.str();
   string xtitle = "Lumisection  [Bx# " + ss1.str() + "]";
  
   hst[histName] = dbe_->book1D(histName, "Number of Good Reconstructed Vertices", 40, 0.5, 40.5);
   hst[histName]->setAxisTitle(xtitle, 1);
   hst[histName]->getTH1()->SetCanExtend(TH1::kAllAxes);
   hst[histName]->getTH1()->SetOption("E1");
 }

References dqm::implementation::IBooker::book1D(), dqm::implementation::DQMStore::cd(), Calorimetry_cff::da, dbe_, debug_, formatFitTime(), dqm::implementation::IGetter::get(), B2GTnPMonitor_cfi::histTitle, hst, mps_fire::i, monitorName_, AlcaSiPixelAliHarvester0T_cff::options, contentValuesCheck::ss, startTime, hgcalPlots::xtitle, and hgcalPlots::ytitle.

Referenced by FitAndFill().

◆ endLuminosityBlock()

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

overrideprotectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 178 of file BeamMonitorBx.cc.

                                                                                                {
   int nthlumi = lumiSeg.id().luminosityBlock();
   edm::LogInfo("LS|BX|BeamMonitorBx") << "Lumi of the last event before endLuminosityBlock: " << nthlumi << endl;
  
   if (nthlumi < nextlumi_)
     return;
   const edm::TimeValue_t fendtimestamp = lumiSeg.endTime().value();
   const std::time_t fendtime = fendtimestamp >> 32;
   tmpTime = refBStime[1] = fendtime;
 }

References edm::LuminosityBlockBase::endTime(), edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), nextlumi_, refBStime, tmpTime, and edm::Timestamp::value().

◆ endRun()

void BeamMonitorBx::endRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

overrideprotectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 566 of file BeamMonitorBx.cc.

566 {}

◆ FillTables()

void BeamMonitorBx::FillTables	(	int	,
		int	,
		std::map< std::string, std::string > &	,
		reco::BeamSpot &	,
		std::string
	)

private

Definition at line 450 of file BeamMonitorBx.cc.

                                                                                                                    {
   map<string, pair<double, double> > val_;
   val_["x0_bx"] = pair<double, double>(bs_.x0(), bs_.x0Error());
   val_["y0_bx"] = pair<double, double>(bs_.y0(), bs_.y0Error());
   val_["z0_bx"] = pair<double, double>(bs_.z0(), bs_.z0Error());
   val_["sigmaX_bx"] = pair<double, double>(bs_.BeamWidthX(), bs_.BeamWidthXError());
   val_["sigmaY_bx"] = pair<double, double>(bs_.BeamWidthY(), bs_.BeamWidthYError());
   val_["sigmaZ_bx"] = pair<double, double>(bs_.sigmaZ(), bs_.sigmaZ0Error());
  
   for (std::map<std::string, std::string>::const_iterator varName = vMap.begin(); varName != vMap.end(); ++varName) {
     TString tmpName = varName->first;
     if (!suffix_.empty())
       tmpName += TString("_" + suffix_);
     hs[tmpName]->setBinContent(1, nthbin_, nthbx);
     hs[tmpName]->setBinContent(2, nthbin_, val_[varName->first].first);
     hs[tmpName]->setBinContent(3, nthbin_, val_[varName->first].second);
   }
 }

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthXError(), reco::BeamSpot::BeamWidthY(), reco::BeamSpot::BeamWidthYError(), hs, reco::BeamSpot::sigmaZ(), reco::BeamSpot::sigmaZ0Error(), reco::BeamSpot::x0(), reco::BeamSpot::x0Error(), reco::BeamSpot::y0(), reco::BeamSpot::y0Error(), reco::BeamSpot::z0(), and reco::BeamSpot::z0Error().

Referenced by FitAndFill().

◆ FillTrendHistos()

void BeamMonitorBx::FillTrendHistos	(	int	,
		int	,
		std::map< std::string, std::string > &	,
		reco::BeamSpot &	,
		const TString &
	)

private

Definition at line 470 of file BeamMonitorBx.cc.

                                                                                                 {
   map<TString, pair<double, double> > val_;
   val_[TString(prefix_ + "_x")] = pair<double, double>(bs_.x0(), bs_.x0Error());
   val_[TString(prefix_ + "_y")] = pair<double, double>(bs_.y0(), bs_.y0Error());
   val_[TString(prefix_ + "_z")] = pair<double, double>(bs_.z0(), bs_.z0Error());
   val_[TString(prefix_ + "_sigmaX")] = pair<double, double>(bs_.BeamWidthX(), bs_.BeamWidthXError());
   val_[TString(prefix_ + "_sigmaY")] = pair<double, double>(bs_.BeamWidthY(), bs_.BeamWidthYError());
   val_[TString(prefix_ + "_sigmaZ")] = pair<double, double>(bs_.sigmaZ(), bs_.sigmaZ0Error());
  
   std::ostringstream ss;
   ss << setfill('0') << setw(5) << nthBx;
   int ntbin_ = tmpTime - startTime;
   for (map<TString, MonitorElement*>::iterator itHst = hst.begin(); itHst != hst.end(); ++itHst) {
     if (!(itHst->first.Contains(ss.str())))
       continue;
     if (itHst->first.Contains("nPVs"))
       continue;
     edm::LogInfo("BX|BeamMonitorBx") << "Filling histogram: " << itHst->first << endl;
     if (itHst->first.Contains("time")) {
       int idx = itHst->first.Index("_time", 5);
       itHst->second->setBinContent(ntbin_, val_[itHst->first(0, idx)].first);
       itHst->second->setBinError(ntbin_, val_[itHst->first(0, idx)].second);
     }
     if (itHst->first.Contains("lumi")) {
       int idx = itHst->first.Index("_lumi", 5);
       itHst->second->setBinContent(endLumiOfBSFit_, val_[itHst->first(0, idx)].first);
       itHst->second->setBinError(endLumiOfBSFit_, val_[itHst->first(0, idx)].second);
     }
   }
   TString histName = "Trending_nPVs_lumi_bx_" + ss.str();
   if (hst[histName])
     hst[histName]->setBinContent(endLumiOfBSFit_, nPVs_);
 }

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthXError(), reco::BeamSpot::BeamWidthY(), reco::BeamSpot::BeamWidthYError(), endLumiOfBSFit_, hst, training_settings::idx, reco::BeamSpot::sigmaZ(), reco::BeamSpot::sigmaZ0Error(), contentValuesCheck::ss, startTime, tmpTime, reco::BeamSpot::x0(), reco::BeamSpot::x0Error(), reco::BeamSpot::y0(), reco::BeamSpot::y0Error(), reco::BeamSpot::z0(), and reco::BeamSpot::z0Error().

Referenced by FitAndFill().

◆ FitAndFill()

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

private

Definition at line 328 of file BeamMonitorBx.cc.

                                                                                                          {
   if (nthlumi <= nextlumi)
     return;
  
   int currentlumi = nextlumi;
   edm::LogInfo("LS|BX|BeamMonitorBx") << "Lumi of the current fit: " << currentlumi << endl;
   lastlumi = currentlumi;
   endLumiOfBSFit_ = currentlumi;
  
   edm::LogInfo("BX|BeamMonitorBx") << "Time lapsed = " << tmpTime - refTime << std::endl;
  
   if (resetHistos_) {
     edm::LogInfo("BX|BeamMonitorBx") << "Resetting Histograms" << endl;
     theBeamFitter->resetCutFlow();
     resetHistos_ = false;
   }
  
   if (fitNLumi_ > 0)
     if (currentlumi % fitNLumi_ != 0)
       return;
  
   std::pair<int, int> fitLS = theBeamFitter->getFitLSRange();
   edm::LogInfo("LS|BX|BeamMonitorBx") << " [Fitter] Do BeamSpot Fit for LS = " << fitLS.first << " to " << fitLS.second
                                       << endl;
   edm::LogInfo("LS|BX|BeamMonitorBx") << " [BX] Do BeamSpot Fit for LS = " << beginLumiOfBSFit_ << " to "
                                       << endLumiOfBSFit_ << endl;
  
   edm::LogInfo("BX|BeamMonitorBx") << " [BxDebugTime] refBStime[0] = " << refBStime[0]
                                    << "; address =  " << &refBStime[0] << std::endl;
   edm::LogInfo("BX|BeamMonitorBx") << " [BxDebugTime] refBStime[1] = " << refBStime[1]
                                    << "; address =  " << &refBStime[1] << std::endl;
  
   if (theBeamFitter->runPVandTrkFitter()) {
     countGoodFit_++;
     edm::LogInfo("BX|BeamMonitorBx") << "Number of good fit = " << countGoodFit_ << endl;
     BeamSpotMapBx bsmap = theBeamFitter->getBeamSpotMap();
     std::map<int, int> npvsmap = theBeamFitter->getNPVsperBX();
     edm::LogInfo("BX|BeamMonitorBx") << "Number of bx = " << bsmap.size() << endl;
     if (bsmap.empty())
       return;
     if (countBx_ < bsmap.size()) {
       countBx_ = bsmap.size();
       BookTables(countBx_, varMap, "");
       BookTables(countBx_, varMap, "all");
       for (BeamSpotMapBx::const_iterator abspot = bsmap.begin(); abspot != bsmap.end(); ++abspot) {
         int bx = abspot->first;
         BookTrendHistos(false, bx, varMap1, "FitBx", "Trending", "bx");
       }
     }
  
     std::pair<int, int> LSRange = theBeamFitter->getFitLSRange();
     char tmpTitle[50];
     sprintf(tmpTitle, "%s %i %s %i %s", " [cm] (LS: ", LSRange.first, " to ", LSRange.second, ")");
     for (std::map<std::string, std::string>::const_iterator varName = varMap.begin(); varName != varMap.end();
          ++varName) {
       hs[varName->first]->setTitle(varName->second + " " + tmpTitle);
       hs[varName->first]->Reset();
     }
  
     if (countGoodFit_ == 1)
       firstlumi_ = LSRange.first;
  
     if (resetFitNLumi_ > 0) {
       char tmpTitle1[60];
       if (countGoodFit_ > 1)
         snprintf(tmpTitle1,
                  sizeof(tmpTitle1),
                  "%s %i %s %i %s",
                  " [cm] (LS: ",
                  firstlumi_,
                  " to ",
                  LSRange.second,
                  ") [weighted average]");
       else
         snprintf(tmpTitle1, sizeof(tmpTitle1), "%s", "Need at least two fits to calculate weighted average");
       for (std::map<std::string, std::string>::const_iterator varName = varMap.begin(); varName != varMap.end();
            ++varName) {
         TString tmpName = varName->first + "_all";
         hs[tmpName]->setTitle(varName->second + " " + tmpTitle1);
         hs[tmpName]->Reset();
       }
     }
  
     int nthBin = countBx_;
     for (BeamSpotMapBx::const_iterator abspot = bsmap.begin(); abspot != bsmap.end(); ++abspot, nthBin--) {
       reco::BeamSpot bs = abspot->second;
       int bx = abspot->first;
       int nPVs = npvsmap.find(bx)->second;
       edm::LogInfo("BeamMonitorBx") << "Number of PVs of bx#[" << bx << "] = " << nPVs << endl;
       // Tables
       FillTables(bx, nthBin, varMap, bs, "");
       // Histograms
       FillTrendHistos(bx, nPVs, varMap1, bs, "Trending");
     }
     // Calculate weighted beam spots
     weight(fbspotMap, bsmap);
     // Fill the results
     nthBin = countBx_;
     if (resetFitNLumi_ > 0 && countGoodFit_ > 1) {
       for (BeamSpotMapBx::const_iterator abspot = fbspotMap.begin(); abspot != fbspotMap.end(); ++abspot, nthBin--) {
         reco::BeamSpot bs = abspot->second;
         int bx = abspot->first;
         FillTables(bx, nthBin, varMap, bs, "all");
       }
     }
   }
   //   else
   //     edm::LogInfo("BeamMonitorBx") << "Bad Fit!!!" << endl;
  
   if (resetFitNLumi_ > 0 && currentlumi % resetFitNLumi_ == 0) {
     edm::LogInfo("LS|BX|BeamMonitorBx") << "Reset track collection for beam fit!!!" << endl;
     resetHistos_ = true;
     theBeamFitter->resetTrkVector();
     theBeamFitter->resetLSRange();
     theBeamFitter->resetRefTime();
     theBeamFitter->resetPVFitter();
     beginLumiOfBSFit_ = 0;
     refBStime[0] = 0;
   }
 }

References beginLumiOfBSFit_, BookTables(), BookTrendHistos(), cms::cuda::bs, l1GtPatternGenerator_cfi::bx, countBx_, countGoodFit_, endLumiOfBSFit_, fbspotMap, FillTables(), FillTrendHistos(), firstlumi_, fitNLumi_, BeamFitter::getBeamSpotMap(), BeamFitter::getFitLSRange(), BeamFitter::getNPVsperBX(), hs, refBStime, refTime, BeamFitter::resetCutFlow(), resetFitNLumi_, resetHistos_, BeamFitter::resetLSRange(), BeamFitter::resetPVFitter(), BeamFitter::resetRefTime(), BeamFitter::resetTrkVector(), BeamFitter::runPVandTrkFitter(), theBeamFitter, tmpTime, varMap, varMap1, and weight().

Referenced by beginLuminosityBlock().

◆ formatFitTime()

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

private

Definition at line 28 of file BeamMonitorBx.cc.

                                                            {
 #define CET (+1)
 #define CEST (+2)
  
   tm* ptm;
   ptm = gmtime(&t);
   sprintf(ts,
           "%4d-%02d-%02d %02d:%02d:%02d",
           ptm->tm_year,
           ptm->tm_mon + 1,
           ptm->tm_mday,
           (ptm->tm_hour + CEST) % 24,
           ptm->tm_min,
           ptm->tm_sec);
  
 #ifdef STRIP_TRAILING_BLANKS_IN_TIMEZONE
   unsigned int b = strlen(ts);
   while (ts[--b] == ' ') {
     ts[b] = 0;
   }
 #endif
 }

References b, CEST, plotBeamSpotDB::ptm, and OrderedSet::t.

Referenced by BookTrendHistos().

◆ weight() [1/2]

void BeamMonitorBx::weight	(	BeamSpotMapBx &	weightedMap_,
		const BeamSpotMapBx &	newMap_
	)

private

Definition at line 506 of file BeamMonitorBx.cc.

                                                                                     {
   for (BeamSpotMapBx::const_iterator it = newMap_.begin(); it != newMap_.end(); ++it) {
     if (weightedMap_.find(it->first) == weightedMap_.end() || (it->second.type() != 2)) {
       weightedMap_[it->first] = it->second;
       continue;
     }
  
     BeamSpot obs = weightedMap_[it->first];
     double val_[8] = {
         obs.x0(), obs.y0(), obs.z0(), obs.sigmaZ(), obs.dxdz(), obs.dydz(), obs.BeamWidthX(), obs.BeamWidthY()};
     double valErr_[8] = {obs.x0Error(),
                          obs.y0Error(),
                          obs.z0Error(),
                          obs.sigmaZ0Error(),
                          obs.dxdzError(),
                          obs.dydzError(),
                          obs.BeamWidthXError(),
                          obs.BeamWidthYError()};
  
     reco::BeamSpot::BeamType type = reco::BeamSpot::Unknown;
     weight(val_[0], valErr_[0], it->second.x0(), it->second.x0Error());
     weight(val_[1], valErr_[1], it->second.y0(), it->second.y0Error());
     weight(val_[2], valErr_[2], it->second.z0(), it->second.z0Error());
     weight(val_[3], valErr_[3], it->second.sigmaZ(), it->second.sigmaZ0Error());
     weight(val_[4], valErr_[4], it->second.dxdz(), it->second.dxdzError());
     weight(val_[5], valErr_[5], it->second.dydz(), it->second.dydzError());
     weight(val_[6], valErr_[6], it->second.BeamWidthX(), it->second.BeamWidthXError());
     weight(val_[7], valErr_[7], it->second.BeamWidthY(), it->second.BeamWidthYError());
     if (it->second.type() == reco::BeamSpot::Tracker) {
       type = reco::BeamSpot::Tracker;
     }
  
     reco::BeamSpot::Point bsPosition(val_[0], val_[1], val_[2]);
     reco::BeamSpot::CovarianceMatrix error;
     for (int i = 0; i < 7; ++i) {
       error(i, i) = valErr_[i] * valErr_[i];
     }
     reco::BeamSpot weightedBeamSpot(bsPosition, val_[3], val_[4], val_[5], val_[6], error, type);
     weightedBeamSpot.setBeamWidthY(val_[7]);
     LogInfo("BX|BeamMonitorBx") << weightedBeamSpot << endl;
     weightedMap_.erase(it->first);
     weightedMap_[it->first] = weightedBeamSpot;
   }
 }

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthXError(), reco::BeamSpot::BeamWidthY(), reco::BeamSpot::BeamWidthYError(), reco::BeamSpot::dxdz(), reco::BeamSpot::dxdzError(), reco::BeamSpot::dydz(), reco::BeamSpot::dydzError(), relativeConstraints::error, mps_fire::i, reco::BeamSpot::setBeamWidthY(), reco::BeamSpot::sigmaZ(), reco::BeamSpot::sigmaZ0Error(), reco::BeamSpot::Tracker, reco::BeamSpot::Unknown, reco::BeamSpot::x0(), reco::BeamSpot::x0Error(), reco::BeamSpot::y0(), reco::BeamSpot::y0Error(), reco::BeamSpot::z0(), and reco::BeamSpot::z0Error().

Referenced by FitAndFill().

◆ weight() [2/2]

void BeamMonitorBx::weight	(	double &	mean,
		double &	meanError,
		const double &	val,
		const double &	valError
	)

private

Definition at line 552 of file BeamMonitorBx.cc.

                                                                                                      {
   double tmpError = 0;
   if (meanError < 1e-8) {
     tmpError = 1 / (valError * valError);
     mean = val * tmpError;
   } else {
     tmpError = 1 / (meanError * meanError) + 1 / (valError * valError);
     mean = mean / (meanError * meanError) + val / (valError * valError);
   }
   mean = mean / tmpError;
   meanError = std::sqrt(1 / tmpError);
 }

References MillePedeFileConverter_cfg::e, SiStripPI::mean, mathSSE::sqrt(), and heppy_batch::val.

Member Data Documentation

◆ beginLumiOfBSFit_

int BeamMonitorBx::beginLumiOfBSFit_

private

Definition at line 75 of file BeamMonitorBx.h.

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

◆ bsSrc_

edm::InputTag BeamMonitorBx::bsSrc_

private

Definition at line 63 of file BeamMonitorBx.h.

Referenced by BeamMonitorBx().

◆ countBx_

unsigned int BeamMonitorBx::countBx_

private

Definition at line 72 of file BeamMonitorBx.h.

Referenced by FitAndFill().

◆ countEvt_

int BeamMonitorBx::countEvt_

private

Definition at line 73 of file BeamMonitorBx.h.

Referenced by analyze().

◆ countGoodFit_

int BeamMonitorBx::countGoodFit_

private

Definition at line 80 of file BeamMonitorBx.h.

Referenced by BeamMonitorBx(), and FitAndFill().

◆ countLumi_

int BeamMonitorBx::countLumi_

private

Definition at line 74 of file BeamMonitorBx.h.

Referenced by beginLuminosityBlock().

◆ dbe_

DQMStore* BeamMonitorBx::dbe_

private

Definition at line 69 of file BeamMonitorBx.h.

Referenced by BeamMonitorBx(), beginJob(), BookTables(), and BookTrendHistos().

◆ debug_

bool BeamMonitorBx::debug_

private

Definition at line 67 of file BeamMonitorBx.h.

Referenced by BookTrendHistos().

◆ endLumiOfBSFit_

int BeamMonitorBx::endLumiOfBSFit_

private

Definition at line 76 of file BeamMonitorBx.h.

Referenced by BeamMonitorBx(), FillTrendHistos(), and FitAndFill().

◆ fbspotMap

BeamSpotMapBx BeamMonitorBx::fbspotMap

private

Definition at line 86 of file BeamMonitorBx.h.

Referenced by FitAndFill().

◆ firstlumi_

int BeamMonitorBx::firstlumi_

private

Definition at line 79 of file BeamMonitorBx.h.

Referenced by BeamMonitorBx(), and FitAndFill().

◆ fitNLumi_

int BeamMonitorBx::fitNLumi_

private

Definition at line 65 of file BeamMonitorBx.h.

Referenced by BeamMonitorBx(), and FitAndFill().

◆ ftimestamp

edm::TimeValue_t BeamMonitorBx::ftimestamp

private

Definition at line 100 of file BeamMonitorBx.h.

Referenced by beginRun().

◆ hs

std::map<TString, MonitorElement*> BeamMonitorBx::hs

private

Definition at line 90 of file BeamMonitorBx.h.

Referenced by BookTables(), FillTables(), and FitAndFill().

◆ hst

std::map<TString, MonitorElement*> BeamMonitorBx::hst

private

Definition at line 91 of file BeamMonitorBx.h.

Referenced by BookTrendHistos(), and FillTrendHistos().

◆ lastlumi_

int BeamMonitorBx::lastlumi_

private

Definition at line 77 of file BeamMonitorBx.h.

Referenced by BeamMonitorBx(), and beginLuminosityBlock().

◆ monitorName_

std::string BeamMonitorBx::monitorName_

private

Definition at line 62 of file BeamMonitorBx.h.

Referenced by BeamMonitorBx(), beginJob(), BookTables(), and BookTrendHistos().

◆ nextlumi_

int BeamMonitorBx::nextlumi_

private

Definition at line 78 of file BeamMonitorBx.h.

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

◆ parameters_

edm::ParameterSet BeamMonitorBx::parameters_

private

Definition at line 61 of file BeamMonitorBx.h.

Referenced by BeamMonitorBx(), Mixins._TypedParameterizable::clone(), Types.PSet::clone(), Types.EDAlias::clone(), and Mixins._TypedParameterizable::copy().

◆ processed_

bool BeamMonitorBx::processed_

private

Definition at line 84 of file BeamMonitorBx.h.

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

◆ refBStime

std::time_t BeamMonitorBx::refBStime[2]

private

Definition at line 81 of file BeamMonitorBx.h.

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

◆ refTime

std::time_t BeamMonitorBx::refTime

private

Definition at line 98 of file BeamMonitorBx.h.

Referenced by beginRun(), and FitAndFill().

◆ resetFitNLumi_

int BeamMonitorBx::resetFitNLumi_

private

Definition at line 66 of file BeamMonitorBx.h.

Referenced by BeamMonitorBx(), and FitAndFill().

◆ resetHistos_

bool BeamMonitorBx::resetHistos_

private

Definition at line 83 of file BeamMonitorBx.h.

Referenced by FitAndFill().

◆ startTime

std::time_t BeamMonitorBx::startTime

private

Definition at line 99 of file BeamMonitorBx.h.

Referenced by beginRun(), BookTrendHistos(), and FillTrendHistos().

◆ theBeamFitter

BeamFitter* BeamMonitorBx::theBeamFitter

private

Definition at line 70 of file BeamMonitorBx.h.

Referenced by analyze(), BeamMonitorBx(), FitAndFill(), and ~BeamMonitorBx().

◆ tmpTime

std::time_t BeamMonitorBx::tmpTime

private

Definition at line 97 of file BeamMonitorBx.h.

Referenced by beginRun(), endLuminosityBlock(), FillTrendHistos(), and FitAndFill().

◆ varMap

std::map<std::string, std::string> BeamMonitorBx::varMap

private

Definition at line 87 of file BeamMonitorBx.h.

Referenced by beginJob(), and FitAndFill().

◆ varMap1

std::map<std::string, std::string> BeamMonitorBx::varMap1

private

Definition at line 88 of file BeamMonitorBx.h.

Referenced by beginJob(), and FitAndFill().

Public Types

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ BeamSpotMapBx

◆ BxNum

◆ DQMStore

◆ MonitorElement

Constructor & Destructor Documentation

◆ BeamMonitorBx()

◆ ~BeamMonitorBx()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ beginLuminosityBlock()

◆ beginRun()

◆ BookTables()

◆ BookTrendHistos()

◆ endLuminosityBlock()

◆ endRun()

◆ FillTables()

◆ FillTrendHistos()

◆ FitAndFill()

◆ formatFitTime()

◆ weight() [1/2]

◆ weight() [2/2]

Member Data Documentation

◆ beginLumiOfBSFit_

◆ bsSrc_

◆ countBx_

◆ countEvt_

◆ countGoodFit_

◆ countLumi_

◆ dbe_

◆ debug_

◆ endLumiOfBSFit_

◆ fbspotMap

◆ firstlumi_

◆ fitNLumi_

◆ ftimestamp

◆ hs

◆ hst

◆ lastlumi_

◆ monitorName_

◆ nextlumi_

◆ parameters_

◆ processed_

◆ refBStime

◆ refTime

◆ resetFitNLumi_

◆ resetHistos_

◆ startTime

◆ theBeamFitter

◆ tmpTime

◆ varMap

◆ varMap1