CMS 3D CMS Logo

List of all members | Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes
FakeBeamMonitor Class Reference

#include <FakeBeamMonitor.h>

Inheritance diagram for FakeBeamMonitor:
DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks > edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... > edm::one::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

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
 
SerialTaskQueueglobalLuminosityBlocksQueue () final
 
SerialTaskQueueglobalRunsQueue () 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 EDProduceroperator= (const EDProducer &)=delete
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () const final
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDProducerBase () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
std::vector< bool > const & recordProvenanceList () const
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
TypeLabelList const & typeLabelList () const
 used by the fwk to register the list of products of this module More...
 
 ~ProducerBase () noexcept(false) override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESResolverIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESResolverIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

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
template<Transition Tr = Transition::Event>
auto produces (std::string instanceName) noexcept
 declare what type of product will make and with which optional label More...
 
template<Transition B>
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
template<BranchType B>
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<class ProductType >
BranchAliasSetterT< ProductType > produces ()
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces ()
 
template<class ProductType >
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces ()
 
template<Transition Tr = Transition::Event>
auto produces () noexcept
 
ProducesCollector producesCollector ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< Bconsumes (edm::InputTag tag) noexcept
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
constexpr auto esConsumes ()
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag)
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
void resetItemsToGetFrom (BranchType iType)
 

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_
 
MonitorElementcutFlowTable
 
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_
 
MonitorElementfitResults
 
int frun
 
edm::TimeValue_t ftimestamp
 
MonitorElementh_d0_phi0
 
MonitorElementh_nTrk_lumi
 
MonitorElementh_nVtx
 
MonitorElementh_nVtx_lumi
 
MonitorElementh_nVtx_lumi_all
 
MonitorElementh_nVtx_st
 
MonitorElementh_PVx [2]
 
MonitorElementh_PVxz
 
MonitorElementh_PVy [2]
 
MonitorElementh_PVyz
 
MonitorElementh_PVz [2]
 
MonitorElementh_sigmaX0
 
MonitorElementh_sigmaY0
 
MonitorElementh_sigmaZ0
 
MonitorElementh_trk_z0
 
MonitorElementh_trkPt
 
MonitorElementh_trkVz
 
MonitorElementh_vx_dz
 
MonitorElementh_vx_vy
 
MonitorElementh_vy_dz
 
MonitorElementh_x0
 
MonitorElementh_y0
 
MonitorElementh_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::OnlineDBOutputServiceonlineDbService_
 
bool onlineMode_
 
const int phiBin_
 
const double phiMax_
 
const double phiMin_
 
reco::BeamSpot preBS
 
bool processed_
 
MonitorElementpvResults
 
std::string recordName_
 
reco::BeamSpot refBS
 
std::time_t refBStime [2]
 
std::time_t refPVtime [2]
 
std::time_t refTime
 
MonitorElementreportSummary
 
Float_t reportSummary_
 
MonitorElementreportSummaryContents [3]
 
MonitorElementreportSummaryMap
 
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
template<typename T >
using BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >
 
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > >
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- 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< DQMTokenrunToken_
 

Detailed Description

Definition at line 34 of file FakeBeamMonitor.h.

Constructor & Destructor Documentation

◆ FakeBeamMonitor()

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

105  : dxBin_(ps.getParameter<int>("dxBin")),
106  dxMin_(ps.getParameter<double>("dxMin")),
107  dxMax_(ps.getParameter<double>("dxMax")),
108 
109  vxBin_(ps.getParameter<int>("vxBin")),
110  vxMin_(ps.getParameter<double>("vxMin")),
111  vxMax_(ps.getParameter<double>("vxMax")),
112 
113  phiBin_(ps.getParameter<int>("phiBin")),
114  phiMin_(ps.getParameter<double>("phiMin")),
115  phiMax_(ps.getParameter<double>("phiMax")),
116 
117  dzBin_(ps.getParameter<int>("dzBin")),
118  dzMin_(ps.getParameter<double>("dzMin")),
119  dzMax_(ps.getParameter<double>("dzMax")),
120 
121  countEvt_(0),
122  countLumi_(0),
123  nthBSTrk_(0),
124  nFitElements_(3),
125  resetHistos_(false),
126  StartAverage_(false),
127  firstAverageFit_(0),
128  countGapLumi_(0) {
129  monitorName_ = ps.getUntrackedParameter<string>("monitorName", "YourSubsystemName");
130  recordName_ = ps.getUntrackedParameter<string>("recordName");
131  intervalInSec_ = ps.getUntrackedParameter<int>("timeInterval", 920); //40 LS X 23"
132  fitNLumi_ = ps.getUntrackedParameter<int>("fitEveryNLumi", -1);
133  resetFitNLumi_ = ps.getUntrackedParameter<int>("resetEveryNLumi", -1);
134  fitPVNLumi_ = ps.getUntrackedParameter<int>("fitPVEveryNLumi", -1);
135  resetPVNLumi_ = ps.getUntrackedParameter<int>("resetPVEveryNLumi", -1);
136  deltaSigCut_ = ps.getUntrackedParameter<double>("deltaSignificanceCut", 15);
137  debug_ = ps.getUntrackedParameter<bool>("Debug");
138  onlineMode_ = ps.getUntrackedParameter<bool>("OnlineMode");
139  min_Ntrks_ = ps.getParameter<ParameterSet>("BeamFitter").getUntrackedParameter<int>("MinimumInputTracks");
140  maxZ_ = ps.getParameter<ParameterSet>("BeamFitter").getUntrackedParameter<double>("MaximumZ");
141  minNrVertices_ = ps.getParameter<ParameterSet>("PVFitter").getUntrackedParameter<unsigned int>("minNrVerticesForFit");
142  minVtxNdf_ = ps.getParameter<ParameterSet>("PVFitter").getUntrackedParameter<double>("minVertexNdf");
143  minVtxWgt_ = ps.getParameter<ParameterSet>("PVFitter").getUntrackedParameter<double>("minVertexMeanWeight");
144  useLockRecords_ = ps.getUntrackedParameter<bool>("useLockRecords");
145 
146  if (!monitorName_.empty())
147  monitorName_ = monitorName_ + "/";
148 
149  if (fitNLumi_ <= 0)
150  fitNLumi_ = 1;
152  refBStime[0] = refBStime[1] = refPVtime[0] = refPVtime[1] = 0;
153  maxZ_ = std::fabs(maxZ_);
154  lastlumi_ = 0;
155  nextlumi_ = 0;
156  processed_ = false;
157 
158  rndm_ = new TRandom3(0);
159 }
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
const double phiMax_
const double vxMax_
const double phiMin_
const double dxMin_
unsigned int nthBSTrk_
std::time_t refBStime[2]
const double vxMin_
unsigned int minNrVertices_
T getUntrackedParameter(std::string const &, T const &) const
std::time_t refPVtime[2]
const double dzMin_
std::string monitorName_
const double dxMax_
const double dzMax_
std::string recordName_
unsigned int min_Ntrks_

◆ ~FakeBeamMonitor()

FakeBeamMonitor::~FakeBeamMonitor ( )
override

Definition at line 102 of file FakeBeamMonitor.cc.

102 { delete rndm_; };

Member Function Documentation

◆ analyze()

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_, iEvent, mapNPV, mapPVx, mapPVy, mapPVz, nextlumi_, onlineMode_, processed_, resetFitNLumi_, and StartAverage_.

655  {
656  const int nthlumi = iEvent.luminosityBlock();
657 
658  if (onlineMode_ && (nthlumi < nextlumi_)) {
659  edm::LogInfo("FakeBeamMonitor") << "analyze:: Spilt event from previous lumi section!" << std::endl;
660  return;
661  }
662  if (onlineMode_ && (nthlumi > nextlumi_)) {
663  edm::LogInfo("FakeBeamMonitor") << "analyze:: Spilt event from next lumi section!!!" << std::endl;
664  return;
665  }
666 
667  countEvt_++;
668  // theBeamFitter->readEvent(
669  // iEvent); //Remember when track fitter read the event in the same place the PVFitter read the events !!!!!!!!!
670 
671  // Handle<reco::BeamSpot> recoBeamSpotHandle;
672  // iEvent.getByToken(bsSrc_, recoBeamSpotHandle);
673  // refBS = *recoBeamSpotHandle;
674 
675  // //------Cut Flow Table filled every event!--------------------------------------
676  // {
677  // // Make a copy of the cut flow table from the beam fitter.
678  // auto tmphisto = static_cast<TH1F*>(theBeamFitter->getCutFlow());
679  // cutFlowTable->getTH1()->SetBins(
680  // tmphisto->GetNbinsX(), tmphisto->GetXaxis()->GetXmin(), tmphisto->GetXaxis()->GetXmax());
681  // // Update the bin labels
682  // if (countEvt_ == 1) // SetLabel just once
683  // for (int n = 0; n < tmphisto->GetNbinsX(); n++)
684  // cutFlowTable->setBinLabel(n + 1, tmphisto->GetXaxis()->GetBinLabel(n + 1), 1);
685  // cutFlowTable->Reset();
686  // cutFlowTable->getTH1()->Add(tmphisto);
687  // }
688 
689  //----Reco tracks -------------------------------------
690  // Handle<reco::TrackCollection> TrackCollection;
691  // iEvent.getByToken(tracksLabel_, TrackCollection);
692  // const reco::TrackCollection* tracks = TrackCollection.product();
693  // for (reco::TrackCollection::const_iterator track = tracks->begin(); track != tracks->end(); ++track) {
694  // h_trkPt->Fill(track->pt()); //no need to change here for average bs
695  // h_trkVz->Fill(track->vz());
696  // }
697 
698  //-------HLT Trigger --------------------------------
699  // edm::Handle<TriggerResults> triggerResults;
700  // bool JetTrigPass = false;
701  // if (iEvent.getByToken(hltSrc_, triggerResults)) {
702  // const edm::TriggerNames& trigNames = iEvent.triggerNames(*triggerResults);
703  // for (unsigned int i = 0; i < triggerResults->size(); i++) {
704  // const std::string& trigName = trigNames.triggerName(i);
705  //
706  // if (JetTrigPass)
707  // continue;
708  //
709  // for (size_t t = 0; t < jetTrigger_.size(); ++t) {
710  // if (JetTrigPass)
711  // continue;
712  //
713  // string string_search(jetTrigger_[t]);
714  // size_t found = trigName.find(string_search);
715  //
716  // if (found != string::npos) {
717  // int thisTrigger_ = trigNames.triggerIndex(trigName);
718  // if (triggerResults->accept(thisTrigger_))
719  // JetTrigPass = true;
720  // } //if trigger found
721  // } //for(t=0;..)
722  // } //for(i=0; ..)
723  // } //if trigger colleciton exist)
724 
725  //------ Primary Vertices-------
726  // edm::Handle<reco::VertexCollection> PVCollection;
727 
728  // if (iEvent.getByToken(pvSrc_, PVCollection)) {
729  int nPVcount_ST = 0; //For Single Trigger(hence ST)
730 
731  // for (reco::VertexCollection::const_iterator pv = PVCollection->begin(); pv != PVCollection->end(); ++pv) {
732  for (int tmp_idx = 0; tmp_idx < 10; tmp_idx++) {
733  //--- vertex selection
734  // if (pv->isFake() || pv->tracksSize() == 0)
735  // continue;
736  //if (JetTrigPass)
737  nPVcount_ST++; //non-fake pv with a specific trigger
738 
739  // if (pv->ndof() < minVtxNdf_ || (pv->ndof() + 3.) / pv->tracksSize() < 2 * minVtxWgt_)
740  // continue;
741 
742  //Fill this map to store xyx for pv so that later we can remove the first one for run aver
743  mapPVx[countLumi_].push_back(tmp_idx);
744  mapPVy[countLumi_].push_back(tmp_idx);
745  mapPVz[countLumi_].push_back(tmp_idx);
746 
747  // if (!StartAverage_) { //for first N LS
748  // h_PVx[0]->Fill(pv->x());
749  // h_PVy[0]->Fill(pv->y());
750  // h_PVz[0]->Fill(pv->z());
751  // h_PVxz->Fill(pv->z(), pv->x());
752  // h_PVyz->Fill(pv->z(), pv->y());
753  // } //for first N LiS
754  // else {
755  // h_PVxz->Fill(pv->z(), pv->x());
756  // h_PVyz->Fill(pv->z(), pv->y());
757  // }
758 
759  } //loop over pvs
760 
761  mapNPV[countLumi_].push_back((nPVcount_ST));
762 
763  // if (!StartAverage_) {
764  // h_nVtx_st->Fill(nPVcount_ST * 1.);
765  // }
766 
767  // } //if pv collection is availaable
768 
769  if (StartAverage_) {
770  map<int, std::vector<float> >::iterator itpvx = mapPVx.begin();
771  map<int, std::vector<float> >::iterator itpvy = mapPVy.begin();
772  map<int, std::vector<float> >::iterator itpvz = mapPVz.begin();
773 
774  map<int, std::vector<int> >::iterator itbspvinfo = mapNPV.begin();
775 
776  if ((int)mapPVx.size() > resetFitNLumi_) { //sometimes the events is not there but LS is there!
777  mapPVx.erase(itpvx);
778  mapPVy.erase(itpvy);
779  mapPVz.erase(itpvz);
780  mapNPV.erase(itbspvinfo);
781  } //loop over Last N lumi collected
782 
783  } //StartAverage==true
784 
785  processed_ = true;
786 }
std::map< int, std::vector< int > > mapNPV
int iEvent
Definition: GenABIO.cc:224
std::map< int, std::vector< float > > mapPVy
std::map< int, std::vector< float > > mapPVz
std::map< int, std::vector< float > > mapPVx
Log< level::Info, false > LogInfo

◆ beginLuminosityBlock()

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

525  {
526  // start DB logger
527  DBloggerReturn_ = 0;
530  onlineDbService_->logger().logInfo() << "FakeBeamMonitor::beginLuminosityBlock - LS: " << lumiSeg.luminosityBlock()
531  << " - Run: " << lumiSeg.getRun().run();
532  }
533 
534  int nthlumi = lumiSeg.luminosityBlock();
535  const edm::TimeValue_t fbegintimestamp = lumiSeg.beginTime().value();
536  const std::time_t ftmptime = fbegintimestamp >> 32;
537 
538  if (countLumi_ == 0 && (!processed_)) {
540  refBStime[0] = refPVtime[0] = ftmptime;
541  mapBeginBSLS[countLumi_] = nthlumi;
542  mapBeginPVLS[countLumi_] = nthlumi;
543  mapBeginBSTime[countLumi_] = ftmptime;
544  mapBeginPVTime[countLumi_] = ftmptime;
545  } //for the first record
546 
547  if (nthlumi > nextlumi_) {
548  if (processed_) {
549  countLumi_++;
550  //store here them will need when we remove the first one of Last N LS
551  mapBeginBSLS[countLumi_] = nthlumi;
552  mapBeginPVLS[countLumi_] = nthlumi;
553  mapBeginBSTime[countLumi_] = ftmptime;
554  mapBeginPVTime[countLumi_] = ftmptime;
555  } //processed passed but not the first lumi
556  if ((!processed_) && countLumi_ != 0) {
557  mapBeginBSLS[countLumi_] = nthlumi;
558  mapBeginPVLS[countLumi_] = nthlumi;
559  mapBeginBSTime[countLumi_] = ftmptime;
560  mapBeginPVTime[countLumi_] = ftmptime;
561  } //processed fails for last lumi
562  } //nthLumi > nextlumi
563 
564  if (StartAverage_) {
565  //Just Make sure it get rest
566  refBStime[0] = 0;
567  refPVtime[0] = 0;
568  beginLumiOfPVFit_ = 0;
569  beginLumiOfBSFit_ = 0;
570 
571  if (debug_)
572  edm::LogInfo("FakeBeamMonitor") << " beginLuminosityBlock: Size of mapBeginBSLS before = "
573  << mapBeginBSLS.size() << endl;
574  if (nthlumi >
575  nextlumi_) { //this make sure that it does not take into account this lumi for fitting and only look forward for new lumi
576  //as countLumi also remains the same so map value get overwritten once return to normal running.
577  //even if few LS are misssing and DQM module do not sees them then it catchs up again
578  map<int, int>::iterator itbs = mapBeginBSLS.begin();
579  map<int, int>::iterator itpv = mapBeginPVLS.begin();
580  map<int, std::time_t>::iterator itbstime = mapBeginBSTime.begin();
581  map<int, std::time_t>::iterator itpvtime = mapBeginPVTime.begin();
582 
583  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
584  mapBeginBSLS.erase(itbs);
585  mapBeginPVLS.erase(itpv);
586  mapBeginBSTime.erase(itbstime);
587  mapBeginPVTime.erase(itpvtime);
588  }
589  /*//not sure if want this or not ??
590  map<int, int>::iterator itgapb=mapBeginBSLS.begin();
591  map<int, int>::iterator itgape=mapBeginBSLS.end(); itgape--;
592  countGapLumi_ = ( (itgape->second) - (itgapb->second) );
593  //if we see Gap more than then 2*resetNFitLumi !!!!!!!
594  //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
595  // so better start as fresh and reset everything like starting in the begining!
596  if(countGapLumi_ >= 2*resetFitNLumi_){RestartFitting(); mapBeginBSLS[countLumi_] = nthlumi;}
597  */
598  }
599 
600  if (debug_)
601  edm::LogInfo("FakeBeamMonitor") << " beginLuminosityBlock:: Size of mapBeginBSLS After = " << mapBeginBSLS.size()
602  << endl;
603 
604  map<int, int>::iterator bbs = mapBeginBSLS.begin();
605  map<int, int>::iterator bpv = mapBeginPVLS.begin();
606  map<int, std::time_t>::iterator bbst = mapBeginBSTime.begin();
607  map<int, std::time_t>::iterator bpvt = mapBeginPVTime.begin();
608 
609  if (beginLumiOfPVFit_ == 0)
610  beginLumiOfPVFit_ = bpv->second; //new begin time after removing the LS
611  if (beginLumiOfBSFit_ == 0)
612  beginLumiOfBSFit_ = bbs->second;
613  if (refBStime[0] == 0)
614  refBStime[0] = bbst->second;
615  if (refPVtime[0] == 0)
616  refPVtime[0] = bpvt->second;
617 
618  } //same logic for average fit as above commented line
619 
620  map<int, std::time_t>::iterator nbbst = mapBeginBSTime.begin();
621  map<int, std::time_t>::iterator nbpvt = mapBeginPVTime.begin();
622 
623  if (onlineMode_ && (nthlumi < nextlumi_))
624  return;
625 
626  if (onlineMode_) {
627  if (nthlumi > nextlumi_) {
628  if (countLumi_ != 0 && processed_)
629  FitAndFill(lumiSeg, lastlumi_, nextlumi_, nthlumi);
630  nextlumi_ = nthlumi;
631  edm::LogInfo("FakeBeamMonitor") << "beginLuminosityBlock:: Next Lumi to Fit: " << nextlumi_ << endl;
632  if ((StartAverage_) && refBStime[0] == 0)
633  refBStime[0] = nbbst->second;
634  if ((StartAverage_) && refPVtime[0] == 0)
635  refPVtime[0] = nbpvt->second;
636  }
637  } else {
638  if (processed_)
639  FitAndFill(lumiSeg, lastlumi_, nextlumi_, nthlumi);
640  nextlumi_ = nthlumi;
641  edm::LogInfo("FakeBeamMonitor") << " beginLuminosityBlock:: Next Lumi to Fit: " << nextlumi_ << endl;
642  if ((StartAverage_) && refBStime[0] == 0)
643  refBStime[0] = nbbst->second;
644  if ((StartAverage_) && refPVtime[0] == 0)
645  refPVtime[0] = nbpvt->second;
646  }
647 
648  //countLumi_++;
649  if (processed_)
650  processed_ = false;
651  edm::LogInfo("FakeBeamMonitor") << " beginLuminosityBlock:: Begin of Lumi: " << nthlumi << endl;
652 }
void FitAndFill(const edm::LuminosityBlock &lumiSeg, int &, int &, int &)
Run const & getRun() const
std::map< int, std::time_t > mapBeginPVTime
std::time_t refBStime[2]
edm::Service< cond::service::OnlineDBOutputService > onlineDbService_
RunNumber_t run() const
Definition: RunBase.h:40
std::time_t refPVtime[2]
cond::persistency::Logger & logger()
unsigned long long TimeValue_t
Definition: Timestamp.h:21
std::map< int, std::time_t > mapBeginBSTime
Log< level::Info, false > LogInfo
std::map< int, int > mapBeginPVLS
TimeValue_t value() const
Definition: Timestamp.h:38
Timestamp const & beginTime() const
std::map< int, int > mapBeginBSLS
bool isAvailable() const
Definition: Service.h:40
EchoedLogStream< edm::LogInfo > logInfo()
Definition: Logger.cc:157
LuminosityBlockNumber_t luminosityBlock() const

◆ bookHistograms()

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

213  {
214  frun = iRun.run();
215  ftimestamp = iRun.beginTime().value();
216  tmpTime = ftimestamp >> 32;
218  char eventTime[64];
219  formatFitTime(eventTime, tmpTime);
220  edm::LogInfo("FakeBeamMonitor") << "TimeOffset = " << eventTime << std::endl;
221  TDatime da(eventTime);
222  if (debug_) {
223  edm::LogInfo("FakeBeamMonitor") << "TimeOffset = ";
224  da.Print();
225  }
226  auto daTime = da.Convert(kTRUE);
227 
228  // book some histograms here
229 
230  // create and cd into new folder
231  iBooker.setCurrentFolder(monitorName_ + "Fit");
232 
233  h_nTrk_lumi = iBooker.book1D("nTrk_lumi", "Num. of selected tracks vs lumi (Fit)", 20, 0.5, 20.5);
234  h_nTrk_lumi->setAxisTitle("Lumisection", 1);
235  h_nTrk_lumi->setAxisTitle("Num of Tracks for Fit", 2);
236 
237  //store vtx vs lumi for monitoring why fits fail
238  h_nVtx_lumi = iBooker.book1D("nVtx_lumi", "Num. of selected Vtx vs lumi (Fit)", 20, 0.5, 20.5);
239  h_nVtx_lumi->setAxisTitle("Lumisection", 1);
240  h_nVtx_lumi->setAxisTitle("Num of Vtx for Fit", 2);
241 
242  h_nVtx_lumi_all = iBooker.book1D("nVtx_lumi_all", "Num. of selected Vtx vs lumi (Fit) all", 20, 0.5, 20.5);
243  h_nVtx_lumi_all->getTH1()->SetCanExtend(TH1::kAllAxes);
244  h_nVtx_lumi_all->setAxisTitle("Lumisection", 1);
245  h_nVtx_lumi_all->setAxisTitle("Num of Vtx for Fit", 2);
246 
247  h_d0_phi0 = iBooker.bookProfile(
248  "d0_phi0", "d_{0} vs. #phi_{0} (Selected Tracks)", phiBin_, phiMin_, phiMax_, dxBin_, dxMin_, dxMax_, "");
249  h_d0_phi0->setAxisTitle("#phi_{0} (rad)", 1);
250  h_d0_phi0->setAxisTitle("d_{0} (cm)", 2);
251 
252  h_vx_vy = iBooker.book2D(
253  "trk_vx_vy", "Vertex (PCA) position of selected tracks", vxBin_, vxMin_, vxMax_, vxBin_, vxMin_, vxMax_);
254  h_vx_vy->setOption("COLZ");
255  // h_vx_vy->getTH1()->SetCanExtend(TH1::kAllAxes);
256  h_vx_vy->setAxisTitle("x coordinate of input track at PCA (cm)", 1);
257  h_vx_vy->setAxisTitle("y coordinate of input track at PCA (cm)", 2);
258 
259  {
260  TDatime* da = new TDatime();
261  gStyle->SetTimeOffset(da->Convert(kTRUE));
262  }
263 
264  const int nvar_ = 6;
265  string coord[nvar_] = {"x", "y", "z", "sigmaX", "sigmaY", "sigmaZ"};
266  string label[nvar_] = {
267  "x_{0} (cm)", "y_{0} (cm)", "z_{0} (cm)", "#sigma_{X_{0}} (cm)", "#sigma_{Y_{0}} (cm)", "#sigma_{Z_{0}} (cm)"};
268 
269  hs.clear();
270  hs.reserve(kNumHists);
271  for (int i = 0; i < 4; i++) {
272  iBooker.setCurrentFolder(monitorName_ + "Fit");
273  for (int ic = 0; ic < nvar_; ++ic) {
274  TString histName(coord[ic]);
275  TString histTitle(coord[ic]);
276  string ytitle(label[ic]);
277  string xtitle("");
278  string options("E1");
279  bool createHisto = true;
280  switch (i) {
281  case 1: // BS vs time
282  histName += "0_time";
283  xtitle = "Time [UTC]";
284  if (ic < 3)
285  histTitle += " coordinate of beam spot vs time (Fit)";
286  else
287  histTitle = histTitle.Insert(5, " ") + " of beam spot vs time (Fit)";
288  break;
289  case 2: // PV vs lumi
290  if (ic < 3) {
291  iBooker.setCurrentFolder(monitorName_ + "PrimaryVertex");
292  histName.Insert(0, "PV");
293  histName += "_lumi";
294  histTitle.Insert(0, "Avg. ");
295  histTitle += " position of primary vtx vs lumi";
296  xtitle = "Lumisection";
297  ytitle.insert(0, "PV");
298  ytitle += " #pm #sigma_{PV";
299  ytitle += coord[ic];
300  ytitle += "} (cm)";
301  } else
302  createHisto = false;
303  break;
304  case 3: // PV vs time
305  if (ic < 3) {
306  iBooker.setCurrentFolder(monitorName_ + "PrimaryVertex");
307  histName.Insert(0, "PV");
308  histName += "_time";
309  histTitle.Insert(0, "Avg. ");
310  histTitle += " position of primary vtx vs time";
311  xtitle = "Time [UTC]";
312  ytitle.insert(0, "PV");
313  ytitle += " #pm #sigma_{PV";
314  ytitle += coord[ic];
315  ytitle += "} (cm)";
316  } else
317  createHisto = false;
318  break;
319  default: // BS vs lumi
320  histName += "0_lumi";
321  xtitle = "Lumisection";
322  if (ic < 3)
323  histTitle += " coordinate of beam spot vs lumi (Fit)";
324  else
325  histTitle = histTitle.Insert(5, " ") + " of beam spot vs lumi (Fit)";
326  break;
327  }
328  if (createHisto) {
329  edm::LogInfo("FakeBeamMonitor") << "hitsName = " << histName << "; histTitle = " << histTitle << std::endl;
330  auto tmpHs = iBooker.book1D(histName, histTitle, 40, 0.5, 40.5);
331  hs.push_back(tmpHs);
332  tmpHs->setAxisTitle(xtitle, 1);
333  tmpHs->setAxisTitle(ytitle, 2);
334  tmpHs->getTH1()->SetOption("E1");
335  if (histName.Contains("time")) {
336  //int nbins = (intervalInSec_/23 > 0 ? intervalInSec_/23 : 40);
337  tmpHs->getTH1()->SetBins(intervalInSec_, 0.5, intervalInSec_ + 0.5);
338  tmpHs->setAxisTimeDisplay(1);
339  tmpHs->setAxisTimeFormat("%H:%M:%S", 1);
340  tmpHs->getTH1()->GetXaxis()->SetTimeOffset(daTime);
341  }
342  histName += "_all";
343  histTitle += " all";
344  tmpHs = iBooker.book1D(histName, histTitle, 40, 0.5, 40.5);
345  hs.push_back(tmpHs);
346  tmpHs->getTH1()->SetCanExtend(TH1::kAllAxes);
347  tmpHs->setAxisTitle(xtitle, 1);
348  tmpHs->setAxisTitle(ytitle, 2);
349  tmpHs->getTH1()->SetOption("E1");
350  if (histName.Contains("time")) {
351  //int nbins = (intervalInSec_/23 > 0 ? intervalInSec_/23 : 40);
352  tmpHs->getTH1()->SetBins(intervalInSec_, 0.5, intervalInSec_ + 0.5);
353  tmpHs->setAxisTimeDisplay(1);
354  tmpHs->setAxisTimeFormat("%H:%M:%S", 1);
355  tmpHs->getTH1()->GetXaxis()->SetTimeOffset(daTime);
356  }
357  }
358  }
359  }
360  assert(hs.size() == kNumHists);
361  assert(0 == strcmp(hs[k_sigmaY0_time]->getTH1()->GetName(), "sigmaY0_time"));
362  assert(0 == strcmp(hs[k_PVz_lumi_all]->getTH1()->GetName(), "PVz_lumi_all"));
363 
364  iBooker.setCurrentFolder(monitorName_ + "Fit");
365 
366  h_trk_z0 = iBooker.book1D("trk_z0", "z_{0} of selected tracks", dzBin_, dzMin_, dzMax_);
367  h_trk_z0->setAxisTitle("z_{0} of selected tracks (cm)", 1);
368 
369  h_vx_dz = iBooker.bookProfile(
370  "vx_dz", "v_{x} vs. dz of selected tracks", dzBin_, dzMin_, dzMax_, dxBin_, dxMin_, dxMax_, "");
371  h_vx_dz->setAxisTitle("dz (cm)", 1);
372  h_vx_dz->setAxisTitle("x coordinate of input track at PCA (cm)", 2);
373 
374  h_vy_dz = iBooker.bookProfile(
375  "vy_dz", "v_{y} vs. dz of selected tracks", dzBin_, dzMin_, dzMax_, dxBin_, dxMin_, dxMax_, "");
376  h_vy_dz->setAxisTitle("dz (cm)", 1);
377  h_vy_dz->setAxisTitle("y coordinate of input track at PCA (cm)", 2);
378 
379  h_x0 = iBooker.book1D("BeamMonitorFeedBack_x0", "x coordinate of beam spot (Fit)", 100, -0.01, 0.01);
380  h_x0->setAxisTitle("x_{0} (cm)", 1);
381  h_x0->getTH1()->SetCanExtend(TH1::kAllAxes);
382 
383  h_y0 = iBooker.book1D("BeamMonitorFeedBack_y0", "y coordinate of beam spot (Fit)", 100, -0.01, 0.01);
384  h_y0->setAxisTitle("y_{0} (cm)", 1);
385  h_y0->getTH1()->SetCanExtend(TH1::kAllAxes);
386 
387  h_z0 = iBooker.book1D("BeamMonitorFeedBack_z0", "z coordinate of beam spot (Fit)", dzBin_, dzMin_, dzMax_);
388  h_z0->setAxisTitle("z_{0} (cm)", 1);
389  h_z0->getTH1()->SetCanExtend(TH1::kAllAxes);
390 
391  h_sigmaX0 = iBooker.book1D("BeamMonitorFeedBack_sigmaX0", "sigma x0 of beam spot (Fit)", 100, 0, 0.05);
392  h_sigmaX0->setAxisTitle("#sigma_{X_{0}} (cm)", 1);
393  h_sigmaX0->getTH1()->SetCanExtend(TH1::kAllAxes);
394 
395  h_sigmaY0 = iBooker.book1D("BeamMonitorFeedBack_sigmaY0", "sigma y0 of beam spot (Fit)", 100, 0, 0.05);
396  h_sigmaY0->setAxisTitle("#sigma_{Y_{0}} (cm)", 1);
397  h_sigmaY0->getTH1()->SetCanExtend(TH1::kAllAxes);
398 
399  h_sigmaZ0 = iBooker.book1D("BeamMonitorFeedBack_sigmaZ0", "sigma z0 of beam spot (Fit)", 100, 0, 10);
400  h_sigmaZ0->setAxisTitle("#sigma_{Z_{0}} (cm)", 1);
401  h_sigmaZ0->getTH1()->SetCanExtend(TH1::kAllAxes);
402 
403  // Histograms of all reco tracks (without cuts):
404  h_trkPt = iBooker.book1D("trkPt", "p_{T} of all reco'd tracks (no selection)", 200, 0., 50.);
405  h_trkPt->setAxisTitle("p_{T} (GeV/c)", 1);
406 
407  h_trkVz = iBooker.book1D("trkVz", "Z coordinate of PCA of all reco'd tracks (no selection)", dzBin_, dzMin_, dzMax_);
408  h_trkVz->setAxisTitle("V_{Z} (cm)", 1);
409 
410  cutFlowTable = iBooker.book1D("cutFlowTable", "Cut flow table of track selection", 9, 0, 9);
411 
412  // Results of previous good fit:
413  fitResults = iBooker.book2D("fitResults", "Results of previous good beam fit", 2, 0, 2, 8, 0, 8);
414  fitResults->setAxisTitle("Fitted Beam Spot (cm)", 1);
415  fitResults->setBinLabel(8, "x_{0}", 2);
416  fitResults->setBinLabel(7, "y_{0}", 2);
417  fitResults->setBinLabel(6, "z_{0}", 2);
418  fitResults->setBinLabel(5, "#sigma_{Z}", 2);
419  fitResults->setBinLabel(4, "#frac{dx}{dz} (rad)", 2);
420  fitResults->setBinLabel(3, "#frac{dy}{dz} (rad)", 2);
421  fitResults->setBinLabel(2, "#sigma_{X}", 2);
422  fitResults->setBinLabel(1, "#sigma_{Y}", 2);
423  fitResults->setBinLabel(1, "Mean", 1);
424  fitResults->setBinLabel(2, "Stat. Error", 1);
425  fitResults->getTH1()->SetOption("text");
426 
427  // Histos of PrimaryVertices:
428  iBooker.setCurrentFolder(monitorName_ + "PrimaryVertex");
429 
430  h_nVtx = iBooker.book1D("vtxNbr", "Reconstructed Vertices(non-fake) in all Event", 60, -0.5, 59.5);
431  h_nVtx->setAxisTitle("Num. of reco. vertices", 1);
432 
433  //For one Trigger only
434  h_nVtx_st = iBooker.book1D("vtxNbr_SelectedTriggers", "Reconstructed Vertices(non-fake) in Events", 60, -0.5, 59.5);
435  //h_nVtx_st->setAxisTitle("Num. of reco. vertices for Un-Prescaled Jet Trigger",1);
436 
437  // Monitor only the PV with highest sum pt of assoc. trks:
438  h_PVx[0] = iBooker.book1D("PVX", "x coordinate of Primary Vtx", 50, -0.01, 0.01);
439  h_PVx[0]->setAxisTitle("PVx (cm)", 1);
440  h_PVx[0]->getTH1()->SetCanExtend(TH1::kAllAxes);
441 
442  h_PVy[0] = iBooker.book1D("PVY", "y coordinate of Primary Vtx", 50, -0.01, 0.01);
443  h_PVy[0]->setAxisTitle("PVy (cm)", 1);
444  h_PVy[0]->getTH1()->SetCanExtend(TH1::kAllAxes);
445 
446  h_PVz[0] = iBooker.book1D("PVZ", "z coordinate of Primary Vtx", dzBin_, dzMin_, dzMax_);
447  h_PVz[0]->setAxisTitle("PVz (cm)", 1);
448 
449  h_PVx[1] = iBooker.book1D("PVXFit", "x coordinate of Primary Vtx (Last Fit)", 50, -0.01, 0.01);
450  h_PVx[1]->setAxisTitle("PVx (cm)", 1);
451  h_PVx[1]->getTH1()->SetCanExtend(TH1::kAllAxes);
452 
453  h_PVy[1] = iBooker.book1D("PVYFit", "y coordinate of Primary Vtx (Last Fit)", 50, -0.01, 0.01);
454  h_PVy[1]->setAxisTitle("PVy (cm)", 1);
455  h_PVy[1]->getTH1()->SetCanExtend(TH1::kAllAxes);
456 
457  h_PVz[1] = iBooker.book1D("PVZFit", "z coordinate of Primary Vtx (Last Fit)", dzBin_, dzMin_, dzMax_);
458  h_PVz[1]->setAxisTitle("PVz (cm)", 1);
459 
460  h_PVxz = iBooker.bookProfile("PVxz", "PVx vs. PVz", dzBin_ / 2, dzMin_, dzMax_, dxBin_ / 2, dxMin_, dxMax_, "");
461  h_PVxz->setAxisTitle("PVz (cm)", 1);
462  h_PVxz->setAxisTitle("PVx (cm)", 2);
463 
464  h_PVyz = iBooker.bookProfile("PVyz", "PVy vs. PVz", dzBin_ / 2, dzMin_, dzMax_, dxBin_ / 2, dxMin_, dxMax_, "");
465  h_PVyz->setAxisTitle("PVz (cm)", 1);
466  h_PVyz->setAxisTitle("PVy (cm)", 2);
467 
468  // Results of previous good fit:
469  pvResults = iBooker.book2D("pvResults", "Results of fitting Primary Vertices", 2, 0, 2, 6, 0, 6);
470  pvResults->setAxisTitle("Fitted Primary Vertex (cm)", 1);
471  pvResults->setBinLabel(6, "PVx", 2);
472  pvResults->setBinLabel(5, "PVy", 2);
473  pvResults->setBinLabel(4, "PVz", 2);
474  pvResults->setBinLabel(3, "#sigma_{X}", 2);
475  pvResults->setBinLabel(2, "#sigma_{Y}", 2);
476  pvResults->setBinLabel(1, "#sigma_{Z}", 2);
477  pvResults->setBinLabel(1, "Mean", 1);
478  pvResults->setBinLabel(2, "Stat. Error", 1);
479  pvResults->getTH1()->SetOption("text");
480 
481  // Summary plots:
482  iBooker.setCurrentFolder(monitorName_ + "EventInfo");
483 
484  reportSummary = iBooker.bookFloat("reportSummary");
485  if (reportSummary)
486  reportSummary->Fill(std::numeric_limits<double>::quiet_NaN());
487 
488  char histo[20];
489  iBooker.setCurrentFolder(monitorName_ + "EventInfo/reportSummaryContents");
490  for (int n = 0; n < nFitElements_; n++) {
491  switch (n) {
492  case 0:
493  sprintf(histo, "x0_status");
494  break;
495  case 1:
496  sprintf(histo, "y0_status");
497  break;
498  case 2:
499  sprintf(histo, "z0_status");
500  break;
501  }
502  reportSummaryContents[n] = iBooker.bookFloat(histo);
503  }
504 
505  for (int i = 0; i < nFitElements_; i++) {
506  summaryContent_[i] = 0.;
507  reportSummaryContents[i]->Fill(std::numeric_limits<double>::quiet_NaN());
508  }
509 
510  iBooker.setCurrentFolder(monitorName_ + "EventInfo");
511 
512  reportSummaryMap = iBooker.book2D("reportSummaryMap", "Beam Spot Summary Map", 1, 0, 1, 3, 0, 3);
514  reportSummaryMap->setAxisTitle("Fitted Beam Spot", 2);
515  reportSummaryMap->setBinLabel(1, " ", 1);
516  reportSummaryMap->setBinLabel(1, "x_{0}", 2);
517  reportSummaryMap->setBinLabel(2, "y_{0}", 2);
518  reportSummaryMap->setBinLabel(3, "z_{0}", 2);
519  for (int i = 0; i < nFitElements_; i++) {
520  reportSummaryMap->setBinContent(1, i + 1, -1.);
521  }
522 }
MonitorElement * pvResults
MonitorElement * h_z0
std::vector< MonitorElement * > hs
const double phiMax_
MonitorElement * h_nTrk_lumi
MonitorElement * h_sigmaY0
const double vxMax_
MonitorElement * cutFlowTable
virtual void setOption(const char *option)
MonitorElement * h_vy_dz
const double phiMin_
const double dxMin_
MonitorElement * reportSummary
MonitorElement * h_nVtx_lumi_all
Float_t summaryContent_[3]
assert(be >=bs)
const double vxMin_
MonitorElement * h_PVz[2]
MonitorElement * h_trkPt
void Fill(long long x)
char const * label
MonitorElement * h_nVtx
MonitorElement * h_trkVz
void formatFitTime(char *, const std::time_t &)
const double dzMin_
std::time_t startTime
MonitorElement * h_sigmaZ0
std::string monitorName_
MonitorElement * h_x0
MonitorElement * h_d0_phi0
virtual void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
const double dxMax_
MonitorElement * h_sigmaX0
edm::TimeValue_t ftimestamp
Log< level::Info, false > LogInfo
virtual void setBinContent(int binx, double content)
set content of bin (1-D)
MonitorElement * h_trk_z0
std::time_t refTime
MonitorElement * reportSummaryContents[3]
MonitorElement * h_nVtx_lumi
MonitorElement * h_vx_vy
MonitorElement * h_vx_dz
const double dzMax_
MonitorElement * h_PVy[2]
MonitorElement * h_y0
MonitorElement * h_PVxz
MonitorElement * h_nVtx_st
std::time_t tmpTime
MonitorElement * reportSummaryMap
MonitorElement * h_PVyz
MonitorElement * h_PVx[2]
MonitorElement * fitResults
virtual void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)

◆ dqmBeginRun()

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

207  {
210  }
211 }
edm::Service< cond::service::OnlineDBOutputService > onlineDbService_
bool isAvailable() const
Definition: Service.h:40

◆ dqmEndRun()

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

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

Definition at line 1552 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_.

1552  {
1553  if (debug_)
1554  edm::LogInfo("FakeBeamMonitor") << "dqmEndRun:: Clearing all the Maps " << endl;
1555  //Clear all the Maps here
1556  mapPVx.clear();
1557  mapPVy.clear();
1558  mapPVz.clear();
1559  mapNPV.clear();
1560  mapBeginBSLS.clear();
1561  mapBeginPVLS.clear();
1562  mapBeginBSTime.clear();
1563  mapBeginPVTime.clear();
1564  mapLSBSTrkSize.clear();
1565  mapLSPVStoreSize.clear();
1566  mapLSCF.clear();
1569  }
1570 }
std::map< int, size_t > mapLSPVStoreSize
std::map< int, std::time_t > mapBeginPVTime
std::map< int, std::size_t > mapLSBSTrkSize
std::map< int, std::vector< int > > mapNPV
edm::Service< cond::service::OnlineDBOutputService > onlineDbService_
std::map< int, std::vector< float > > mapPVy
std::map< int, std::vector< float > > mapPVz
std::map< int, std::vector< float > > mapPVx
std::map< int, std::time_t > mapBeginBSTime
Log< level::Info, false > LogInfo
std::map< int, int > mapBeginPVLS
std::map< int, TH1F > mapLSCF
std::map< int, int > mapBeginBSLS
bool isAvailable() const
Definition: Service.h:40

◆ endLuminosityBlock()

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

Definition at line 789 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().

789  {
790  int nthlumi = lumiSeg.id().luminosityBlock();
791  edm::LogInfo("FakeBeamMonitor") << "endLuminosityBlock:: Lumi of the last event before endLuminosityBlock: "
792  << nthlumi << endl;
793 
794  if (onlineMode_ && nthlumi < nextlumi_)
795  return;
796  const edm::TimeValue_t fendtimestamp = lumiSeg.endTime().value();
797  const std::time_t fendtime = fendtimestamp >> 32;
798  tmpTime = refBStime[1] = refPVtime[1] = fendtime;
799 
800  // end DB logger
802  onlineDbService_->logger().logInfo() << "FakeBeamMonitor::endLuminosityBlock";
804  }
805 }
LuminosityBlockNumber_t luminosityBlock() const
std::time_t refBStime[2]
edm::Service< cond::service::OnlineDBOutputService > onlineDbService_
std::time_t refPVtime[2]
cond::persistency::Logger & logger()
unsigned long long TimeValue_t
Definition: Timestamp.h:21
Log< level::Info, false > LogInfo
LuminosityBlockID id() const
TimeValue_t value() const
Definition: Timestamp.h:38
Timestamp const & endTime() const
std::time_t tmpTime
bool isAvailable() const
Definition: Service.h:40
EchoedLogStream< edm::LogInfo > logInfo()
Definition: Logger.cc:157
void end(int retCode)
Definition: Logger.cc:80

◆ FitAndFill()

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

Definition at line 808 of file FakeBeamMonitor.cc.

References simBeamSpotPI::alpha, reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthXError(), reco::BeamSpot::BeamWidthY(), reco::BeamSpot::BeamWidthYError(), beginLumiOfBSFit_, beginLumiOfPVFit_, cms::cuda::bs, submitPVResolutionJobs::count, countLumi_, beamSpotPI::creationTime, DBloggerReturn_, MillePedeFileConverter_cfg::e, mps_fire::end, endLumiOfBSFit_, endLumiOfPVFit_, BeamSpotOnlineObjects::endTime(), BeamSpotOnlineObjects::endTimeStamp(), cppFunctionSkipper::exception, DeadROC_duringRun::f1, dqm::impl::MonitorElement::Fill(), firstAverageFit_, fitNLumi_, fitPVNLumi_, fitResults, frun, getGMTstring(), dqm::impl::MonitorElement::getTH1(), dqm::impl::MonitorElement::getTH1F(), dqm::impl::MonitorElement::getTProfile(), h_d0_phi0, h_nTrk_lumi, h_nVtx_lumi, h_nVtx_st, h_PVx, h_PVy, h_PVz, h_sigmaX0, h_sigmaY0, h_sigmaZ0, h_trk_z0, h_vx_dz, h_vx_vy, h_vy_dz, h_x0, h_y0, h_z0, hs, mps_fire::i, cuy::ii, edm::Service< T >::isAvailable(), dqmiolumiharvest::j, dqmdumpme::k, cond::persistency::Logger::logError(), cond::service::PoolDBOutputService::logger(), cond::persistency::Logger::logInfo(), genParticles_cff::map, mapLSPVStoreSize, mapNPV, mapPVx, mapPVy, mapPVz, makeMuonMisalignmentScenario::matrix, BeamSpotOnlineObjects::maxPVs(), SiStripPI::mean, BeamSpotOnlineObjects::meanErrorPV(), BeamSpotOnlineObjects::meanPV(), min_Ntrks_, minNrVertices_, dqmiodumpmetadata::n, nFitElements_, nFits_, submitPVValidationJobs::now, nthBSTrk_, onlineDbService_, onlineMode_, preBS, pvResults, recordName_, refBStime, refPVtime, refTime, reportSummary, reportSummary_, reportSummaryContents, reportSummaryMap, dqm::impl::MonitorElement::Reset(), resetFitNLumi_, resetHistos_, resetPVNLumi_, BeamSpotOnlineObjects::rmsErrorPV(), BeamSpotOnlineObjects::rmsPV(), rndm_, scrollTH1(), dqm::impl::MonitorElement::setAxisTitle(), BeamSpotObjects::setBeamWidthX(), BeamSpotObjects::setBeamWidthXError(), BeamSpotObjects::setBeamWidthY(), BeamSpotObjects::setBeamWidthYError(), BeamSpotObjects::setBetaStar(), dqm::impl::MonitorElement::setBinContent(), BeamSpotObjects::setCovariance(), BeamSpotOnlineObjects::setCreationTime(), BeamSpotObjects::setdxdz(), BeamSpotObjects::setdydz(), BeamSpotObjects::setEmittanceX(), BeamSpotObjects::setEmittanceY(), BeamSpotOnlineObjects::setEndTime(), BeamSpotOnlineObjects::setEndTimeStamp(), BeamSpotOnlineObjects::setLastAnalyzedFill(), BeamSpotOnlineObjects::setLastAnalyzedLumi(), BeamSpotOnlineObjects::setLastAnalyzedRun(), BeamSpotOnlineObjects::setMaxPVs(), BeamSpotOnlineObjects::setMeanErrorPV(), BeamSpotOnlineObjects::setMeanPV(), BeamSpotOnlineObjects::setNumPVs(), BeamSpotOnlineObjects::setNumTracks(), BeamSpotObjects::setPosition(), BeamSpotOnlineObjects::setRmsErrorPV(), BeamSpotOnlineObjects::setRmsPV(), BeamSpotObjects::setSigmaZ(), BeamSpotOnlineObjects::setStartTime(), BeamSpotOnlineObjects::setStartTimeStamp(), BeamSpotObjects::setType(), BeamSpotOnlineObjects::setUsedEvents(), dqm::impl::MonitorElement::ShiftFillLast(), mathSSE::sqrt(), StartAverage_, BeamSpotOnlineObjects::startTime(), startTime, BeamSpotOnlineObjects::startTimeStamp(), summaryContent_, summarySum_, testScroll(), tmpTime, reco::BeamSpot::Tracker, BeamSpotOnlineObjects::usedEvents(), ApeEstimator_cff::width, and cond::service::OnlineDBOutputService::writeIOVForNextLumisection().

Referenced by beginLuminosityBlock().

808  {
809  if (onlineMode_ && (nthlumi <= nextlumi))
810  return;
811 
812  //set the correct run number when no event in the LS for fake output
813  // if ((processed_) && theBeamFitter->getRunNumber() != frun)
814  // theBeamFitter->setRun(frun);
815 
816  int currentlumi = nextlumi;
817  edm::LogInfo("FakeBeamMonitor") << "FitAndFill:: Lumi of the current fit: " << currentlumi << endl;
818  lastlumi = currentlumi;
819  endLumiOfBSFit_ = currentlumi;
820  endLumiOfPVFit_ = currentlumi;
821 
822  //---------Fix for Runninv average-------------
823  mapLSPVStoreSize[countLumi_] = 10; //theBeamFitter->getPVvectorSize();
824 
825  // if (StartAverage_) {
826  // std::map<int, std::size_t>::iterator rmLSPVi = mapLSPVStoreSize.begin();
827  // size_t SizeToRemovePV = rmLSPVi->second;
828  // for (std::map<int, std::size_t>::iterator rmLSPVe = mapLSPVStoreSize.end(); ++rmLSPVi != rmLSPVe;)
829  // rmLSPVi->second -= SizeToRemovePV;
830  //
831  // theBeamFitter->resizePVvector(SizeToRemovePV);
832  //
833  // map<int, std::size_t>::iterator tmpItpv = mapLSPVStoreSize.begin();
834  // mapLSPVStoreSize.erase(tmpItpv);
835  // }
836  // if (debug_)
837  // edm::LogInfo("BeamMonitor") << "FitAndFill:: Size of thePVvector After removing the PVs = "
838  // << theBeamFitter->getPVvectorSize() << endl;
839 
840  //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
841  bool resetHistoFlag_ = false;
842  if ((int)mapPVx.size() >= resetFitNLumi_ && (StartAverage_)) {
843  h_PVx[0]->Reset();
844  h_PVy[0]->Reset();
845  h_PVz[0]->Reset();
846  h_nVtx_st->Reset();
847  resetHistoFlag_ = true;
848  }
849 
850  int MaxPVs = 0;
851 
852  std::map<int, std::vector<int> >::iterator mnpv = mapNPV.begin();
853  std::map<int, std::vector<float> >::iterator mpv2 = mapPVy.begin();
854  std::map<int, std::vector<float> >::iterator mpv3 = mapPVz.begin();
855 
856  for (std::map<int, std::vector<float> >::iterator mpv1 = mapPVx.begin(); mpv1 != mapPVx.end();
857  ++mpv1, ++mpv2, ++mpv3, ++mnpv) {
858  std::vector<float>::iterator mpvs2 = (mpv2->second).begin();
859  std::vector<float>::iterator mpvs3 = (mpv3->second).begin();
860  for (std::vector<float>::iterator mpvs1 = (mpv1->second).begin(); mpvs1 != (mpv1->second).end();
861  ++mpvs1, ++mpvs2, ++mpvs3) {
862  if (resetHistoFlag_) {
863  h_PVx[0]->Fill(*mpvs1); //these histogram are reset after StartAverage_ flag is ON
864  h_PVy[0]->Fill(*mpvs2);
865  h_PVz[0]->Fill(*mpvs3);
866  }
867  } //loop over second
868 
869  //Do the same here for nPV distr.
870  for (std::vector<int>::iterator mnpvs = (mnpv->second).begin(); mnpvs != (mnpv->second).end(); ++mnpvs) {
871  if ((*mnpvs > 0) && (resetHistoFlag_))
872  h_nVtx_st->Fill((*mnpvs) * (1.0));
873  if ((*mnpvs) > MaxPVs)
874  MaxPVs = (*mnpvs);
875  } //loop over second of mapNPV
876 
877  } //loop over last N lumis
878 
879  char tmpTitlePV[100];
880  sprintf(tmpTitlePV, "%s %i %s %i", "Num. of reco. vertices for LS: ", beginLumiOfPVFit_, " to ", endLumiOfPVFit_);
881  h_nVtx_st->setAxisTitle(tmpTitlePV, 1);
882 
883  std::vector<float> DipPVInfo_;
884  DipPVInfo_.clear();
885  DipPVInfo_.push_back(rndm_->Gaus(1000., 100.)); // Events used
886  DipPVInfo_.push_back(rndm_->Gaus(100., 10.)); // Mean PV
887  DipPVInfo_.push_back(rndm_->Gaus(10., 5.)); // Mean PV err
888  DipPVInfo_.push_back(rndm_->Gaus(10., 5.)); // Rms PV
889  DipPVInfo_.push_back(rndm_->Gaus(5., 3.)); // Rms PV err
890  DipPVInfo_.push_back(rndm_->Gaus(100., 10.)); // Max PVs
891 
892  if (onlineMode_) { // filling LS gap
893  // FIXME: need to add protection for the case if the gap is at the resetting LS!
894  const int countLS_bs = hs[k_x0_lumi]->getTH1()->GetEntries();
895  const int countLS_pv = hs[k_PVx_lumi]->getTH1()->GetEntries();
896  edm::LogInfo("FakeBeamMonitor") << "FitAndFill:: countLS_bs = " << countLS_bs << " ; countLS_pv = " << countLS_pv
897  << std::endl;
898  int LSgap_bs = currentlumi / fitNLumi_ - countLS_bs;
899  int LSgap_pv = currentlumi / fitPVNLumi_ - countLS_pv;
900  if (currentlumi % fitNLumi_ == 0)
901  LSgap_bs--;
902  if (currentlumi % fitPVNLumi_ == 0)
903  LSgap_pv--;
904  edm::LogInfo("FakeBeamMonitor") << "FitAndFill:: LSgap_bs = " << LSgap_bs << " ; LSgap_pv = " << LSgap_pv
905  << std::endl;
906  // filling previous fits if LS gap ever exists
907  for (int ig = 0; ig < LSgap_bs; ig++) {
908  hs[k_x0_lumi]->ShiftFillLast(0., 0., fitNLumi_); //x0 , x0err, fitNLumi_; see DQMCore....
909  hs[k_y0_lumi]->ShiftFillLast(0., 0., fitNLumi_);
910  hs[k_z0_lumi]->ShiftFillLast(0., 0., fitNLumi_);
911  hs[k_sigmaX0_lumi]->ShiftFillLast(0., 0., fitNLumi_);
912  hs[k_sigmaY0_lumi]->ShiftFillLast(0., 0., fitNLumi_);
913  hs[k_sigmaZ0_lumi]->ShiftFillLast(0., 0., fitNLumi_);
915  }
916  for (int ig = 0; ig < LSgap_pv; ig++) {
917  hs[k_PVx_lumi]->ShiftFillLast(0., 0., fitPVNLumi_);
918  hs[k_PVy_lumi]->ShiftFillLast(0., 0., fitPVNLumi_);
919  hs[k_PVz_lumi]->ShiftFillLast(0., 0., fitPVNLumi_);
920  }
921  const int previousLS = h_nTrk_lumi->getTH1()->GetEntries();
922  for (int i = 1; i < (currentlumi - previousLS);
923  i++) //if (current-previoius)= 1 then never go inside the for loop!!!!!!!!!!!
925  }
926 
927  edm::LogInfo("FakeBeamMonitor") << "FitAndFill:: Time lapsed since last scroll = " << tmpTime - refTime << std::endl;
928 
929  if (testScroll(tmpTime, refTime)) {
930  scrollTH1(hs[k_x0_time]->getTH1(), refTime);
931  scrollTH1(hs[k_y0_time]->getTH1(), refTime);
932  scrollTH1(hs[k_z0_time]->getTH1(), refTime);
933  scrollTH1(hs[k_sigmaX0_time]->getTH1(), refTime);
934  scrollTH1(hs[k_sigmaY0_time]->getTH1(), refTime);
935  scrollTH1(hs[k_sigmaZ0_time]->getTH1(), refTime);
936  scrollTH1(hs[k_PVx_time]->getTH1(), refTime);
937  scrollTH1(hs[k_PVy_time]->getTH1(), refTime);
938  scrollTH1(hs[k_PVz_time]->getTH1(), refTime);
939  }
940 
941  // bool doPVFit = false;
942  //
943  // if (fitPVNLumi_ > 0) {
944  // if (onlineMode_) {
945  // if (currentlumi % fitPVNLumi_ == 0)
946  // doPVFit = true;
947  // } else if (countLumi_ % fitPVNLumi_ == 0)
948  // doPVFit = true;
949  // } else
950  // doPVFit = true;
951  //
952  // if (doPVFit) {
953  edm::LogInfo("FakeBeamMonitor") << "FitAndFill:: Do PV Fitting for LS = " << beginLumiOfPVFit_ << " to "
954  << endLumiOfPVFit_ << std::endl;
955  // Primary Vertex Fit:
956  if (h_PVx[0]->getTH1()->GetEntries() > minNrVertices_) {
957  pvResults->Reset();
958  char tmpTitle[50];
959  sprintf(tmpTitle, "%s %i %s %i", "Fitted Primary Vertex (cm) of LS: ", beginLumiOfPVFit_, " to ", endLumiOfPVFit_);
960  pvResults->setAxisTitle(tmpTitle, 1);
961 
962  std::unique_ptr<TF1> fgaus{new TF1("fgaus", "gaus")};
963  double mean, width, meanErr, widthErr;
964  fgaus->SetLineColor(4);
965  h_PVx[0]->getTH1()->Fit(fgaus.get(), "QLM0");
966  mean = fgaus->GetParameter(1);
967  width = fgaus->GetParameter(2);
968  meanErr = fgaus->GetParError(1);
969  widthErr = fgaus->GetParError(2);
970 
971  hs[k_PVx_lumi]->ShiftFillLast(mean, width, fitPVNLumi_);
972  hs[k_PVx_lumi_all]->setBinContent(currentlumi, mean);
973  hs[k_PVx_lumi_all]->setBinError(currentlumi, width);
974  int nthBin = tmpTime - refTime;
975  if (nthBin < 0)
976  edm::LogInfo("FakeBeamMonitor") << "FitAndFill:: Event time outside current range of time histograms!"
977  << std::endl;
978  if (nthBin > 0) {
979  hs[k_PVx_time]->setBinContent(nthBin, mean);
980  hs[k_PVx_time]->setBinError(nthBin, width);
981  }
982  int jthBin = tmpTime - startTime;
983  if (jthBin > 0) {
984  hs[k_PVx_time_all]->setBinContent(jthBin, mean);
985  hs[k_PVx_time_all]->setBinError(jthBin, width);
986  }
987  pvResults->setBinContent(1, 6, mean);
988  pvResults->setBinContent(1, 3, width);
989  pvResults->setBinContent(2, 6, meanErr);
990  pvResults->setBinContent(2, 3, widthErr);
991 
992  {
993  // snap shot of the fit
994  auto tmphisto = h_PVx[0]->getTH1F();
995  h_PVx[1]->getTH1()->SetBins(
996  tmphisto->GetNbinsX(), tmphisto->GetXaxis()->GetXmin(), tmphisto->GetXaxis()->GetXmax());
997  h_PVx[1]->Reset();
998  h_PVx[1]->getTH1()->Add(tmphisto);
999  h_PVx[1]->getTH1()->Fit(fgaus.get(), "QLM");
1000  }
1001 
1002  h_PVy[0]->getTH1()->Fit(fgaus.get(), "QLM0");
1003  mean = fgaus->GetParameter(1);
1004  width = fgaus->GetParameter(2);
1005  meanErr = fgaus->GetParError(1);
1006  widthErr = fgaus->GetParError(2);
1007  hs[k_PVy_lumi]->ShiftFillLast(mean, width, fitPVNLumi_);
1008  hs[k_PVy_lumi_all]->setBinContent(currentlumi, mean);
1009  hs[k_PVy_lumi_all]->setBinError(currentlumi, width);
1010  if (nthBin > 0) {
1011  hs[k_PVy_time]->setBinContent(nthBin, mean);
1012  hs[k_PVy_time]->setBinError(nthBin, width);
1013  }
1014  if (jthBin > 0) {
1015  hs[k_PVy_time_all]->setBinContent(jthBin, mean);
1016  hs[k_PVy_time_all]->setBinError(jthBin, width);
1017  }
1018  pvResults->setBinContent(1, 5, mean);
1019  pvResults->setBinContent(1, 2, width);
1020  pvResults->setBinContent(2, 5, meanErr);
1021  pvResults->setBinContent(2, 2, widthErr);
1022  // snap shot of the fit
1023  {
1024  auto tmphisto = h_PVy[0]->getTH1F();
1025  h_PVy[1]->getTH1()->SetBins(
1026  tmphisto->GetNbinsX(), tmphisto->GetXaxis()->GetXmin(), tmphisto->GetXaxis()->GetXmax());
1027  h_PVy[1]->Reset();
1028  h_PVy[1]->getTH1()->Add(tmphisto);
1029  h_PVy[1]->getTH1()->Fit(fgaus.get(), "QLM");
1030  }
1031 
1032  h_PVz[0]->getTH1()->Fit(fgaus.get(), "QLM0");
1033  mean = fgaus->GetParameter(1);
1034  width = fgaus->GetParameter(2);
1035  meanErr = fgaus->GetParError(1);
1036  widthErr = fgaus->GetParError(2);
1037  hs[k_PVz_lumi]->ShiftFillLast(mean, width, fitPVNLumi_);
1038  hs[k_PVz_lumi_all]->setBinContent(currentlumi, mean);
1039  hs[k_PVz_lumi_all]->setBinError(currentlumi, width);
1040  if (nthBin > 0) {
1041  hs[k_PVz_time]->setBinContent(nthBin, mean);
1042  hs[k_PVz_time]->setBinError(nthBin, width);
1043  }
1044  if (jthBin > 0) {
1045  hs[k_PVz_time_all]->setBinContent(jthBin, mean);
1046  hs[k_PVz_time_all]->setBinError(jthBin, width);
1047  }
1048  pvResults->setBinContent(1, 4, mean);
1049  pvResults->setBinContent(1, 1, width);
1050  pvResults->setBinContent(2, 4, meanErr);
1051  pvResults->setBinContent(2, 1, widthErr);
1052  {
1053  // snap shot of the fit
1054  auto tmphisto = h_PVz[0]->getTH1F();
1055  h_PVz[1]->getTH1()->SetBins(
1056  tmphisto->GetNbinsX(), tmphisto->GetXaxis()->GetXmin(), tmphisto->GetXaxis()->GetXmax());
1057  h_PVz[1]->Reset();
1058  h_PVz[1]->getTH1()->Add(tmphisto);
1059  h_PVz[1]->getTH1()->Fit(fgaus.get(), "QLM");
1060  }
1061  } //check if found min Vertices
1062  // } //do PVfit
1063 
1064  if ((resetPVNLumi_ > 0 && countLumi_ == resetPVNLumi_) || StartAverage_) {
1065  beginLumiOfPVFit_ = 0;
1066  refPVtime[0] = 0;
1067  }
1068 
1069  //---------Readjustment of theBSvector, RefTime, beginLSofFit---------
1070  // vector<BSTrkParameters> theBSvector1 = theBeamFitter->getBSvector();
1071  // mapLSBSTrkSize[countLumi_] = (theBSvector1.size());
1072  size_t PreviousRecords = 0; //needed to fill nth record of tracks in GUI
1073 
1074  // if (StartAverage_) {
1075  // size_t SizeToRemove = 0;
1076  // std::map<int, std::size_t>::iterator rmls = mapLSBSTrkSize.begin();
1077  // SizeToRemove = rmls->second;
1078  // if (debug_)
1079  // edm::LogInfo("BeamMonitor") << " The size to remove is = " << SizeToRemove << endl;
1080  // int changedAfterThis = 0;
1081  // for (std::map<int, std::size_t>::iterator rmLS = mapLSBSTrkSize.begin(); rmLS != mapLSBSTrkSize.end();
1082  // ++rmLS, ++changedAfterThis) {
1083  // if (changedAfterThis > 0) {
1084  // (rmLS->second) = (rmLS->second) - SizeToRemove;
1085  // if ((mapLSBSTrkSize.size() - (size_t)changedAfterThis) == 2)
1086  // PreviousRecords = (rmLS->second);
1087  // }
1088  // }
1089  //
1090  // theBeamFitter->resizeBSvector(SizeToRemove);
1091  //
1092  // map<int, std::size_t>::iterator tmpIt = mapLSBSTrkSize.begin();
1093  // mapLSBSTrkSize.erase(tmpIt);
1094  //
1095  // std::pair<int, int> checkfitLS = theBeamFitter->getFitLSRange();
1096  // std::pair<time_t, time_t> checkfitTime = theBeamFitter->getRefTime();
1097  // theBeamFitter->setFitLSRange(beginLumiOfBSFit_, checkfitLS.second);
1098  // theBeamFitter->setRefTime(refBStime[0], checkfitTime.second);
1099  // }
1100 
1101  //Fill the track for this fit
1102  // vector<BSTrkParameters> theBSvector = theBeamFitter->getBSvector();
1103  // h_nTrk_lumi->ShiftFillLast(theBSvector.size());
1104  //
1105  // if (debug_)
1106  // edm::LogInfo("BeamMonitor") << "FitAndFill:: Size of theBSViector.size() After =" << theBSvector.size() << endl;
1107 
1108  // bool countFitting = false;
1109  // if (theBSvector.size() >= PreviousRecords && theBSvector.size() >= min_Ntrks_) {
1110  // countFitting = true;
1111  // }
1112 
1113  //---Fix for Cut Flow Table for Running average in a same way//the previous code has problem for resetting!!!
1114  // mapLSCF[countLumi_] = *theBeamFitter->getCutFlow();
1115  // if (StartAverage_ && !mapLSCF.empty()) {
1116  // const TH1F& cutFlowToSubtract = mapLSCF.begin()->second;
1117  // // Subtract the last cut flow from all of the others.
1118  // std::map<int, TH1F>::iterator cf = mapLSCF.begin();
1119  // // Start on second entry
1120  // for (; cf != mapLSCF.end(); ++cf) {
1121  // cf->second.Add(&cutFlowToSubtract, -1);
1122  // }
1123  // theBeamFitter->subtractFromCutFlow(&cutFlowToSubtract);
1124  // // Remove the obsolete lumi section
1125  // mapLSCF.erase(mapLSCF.begin());
1126  // }
1127 
1128  if (resetHistos_) {
1129  h_d0_phi0->Reset();
1130  h_vx_vy->Reset();
1131  h_vx_dz->Reset();
1132  h_vy_dz->Reset();
1133  h_trk_z0->Reset();
1134  resetHistos_ = false;
1135  }
1136 
1137  if (StartAverage_)
1138  nthBSTrk_ = PreviousRecords; //after average proccess is ON//for 2-6 LS fit PreviousRecords is size from 2-5 LS
1139 
1140  edm::LogInfo("FakeBeamMonitor") << " The Previous Recored for this fit is =" << nthBSTrk_ << endl;
1141 
1142  // unsigned int itrk = 0;
1143  // for (vector<BSTrkParameters>::const_iterator BSTrk = theBSvector.begin(); BSTrk != theBSvector.end();
1144  // ++BSTrk, ++itrk) {
1145  // if (itrk >= nthBSTrk_) { //fill for this record only !!
1146  // h_d0_phi0->Fill(BSTrk->phi0(), BSTrk->d0());
1147  // double vx = BSTrk->vx();
1148  // double vy = BSTrk->vy();
1149  // double z0 = BSTrk->z0();
1150  // h_vx_vy->Fill(vx, vy);
1151  // h_vx_dz->Fill(z0, vx);
1152  // h_vy_dz->Fill(z0, vy);
1153  // h_trk_z0->Fill(z0);
1154  // }
1155  // }
1156 
1157  // nthBSTrk_ = theBSvector.size(); // keep track of num of tracks filled so far
1158 
1159  edm::LogInfo("FakeBeamMonitor") << " The Current Recored for this fit is =" << nthBSTrk_ << endl;
1160 
1161  // if (countFitting)
1162  // edm::LogInfo("FakeBeamMonitor") << "FitAndFill:: Num of tracks collected = " << nthBSTrk_ << endl;
1163 
1164  if (fitNLumi_ > 0) {
1165  if (onlineMode_) {
1166  if (currentlumi % fitNLumi_ != 0) {
1167  // for (std::map<TString,MonitorElement*>::iterator itAll = hs.begin();
1168  // itAll != hs.end(); ++itAll) {
1169  // if ((*itAll).first.Contains("all")) {
1170  // (*itAll).second->setBinContent(currentlumi,0.);
1171  // (*itAll).second->setBinError(currentlumi,0.);
1172  // }
1173  // }
1174  return;
1175  }
1176  } else if (countLumi_ % fitNLumi_ != 0)
1177  return;
1178  }
1179 
1180  edm::LogInfo("FakeBeamMonitor") << "FitAndFill:: [DebugTime] refBStime[0] = " << refBStime[0]
1181  << "; address = " << &refBStime[0] << std::endl;
1182  edm::LogInfo("FakeBeamMonitor") << "FitAndFill:: [DebugTime] refBStime[1] = " << refBStime[1]
1183  << "; address = " << &refBStime[1] << std::endl;
1184 
1185  //Fill for all LS even if fit fails
1186  // h_nVtx_lumi->ShiftFillLast((theBeamFitter->getPVvectorSize()), 0., fitNLumi_);
1187  // h_nVtx_lumi_all->setBinContent(currentlumi, (theBeamFitter->getPVvectorSize()));
1188 
1189  // if (countFitting) {
1190  nFits_++;
1191  // std::pair<int, int> fitLS = theBeamFitter->getFitLSRange();
1192  std::pair<int, int> fitLS(beginLumiOfBSFit_, endLumiOfBSFit_);
1193  // edm::LogInfo("BeamMonitor") << "FitAndFill:: [BeamFitter] Do BeamSpot Fit for LS = " << fitLS.first << " to "
1194  // << fitLS.second << std::endl;
1195  edm::LogInfo("FakeBeamMonitor") << "FitAndFill:: [FakeBeamMonitor] Do BeamSpot Fit for LS = " << beginLumiOfBSFit_
1196  << " to " << endLumiOfBSFit_ << std::endl;
1197 
1198  //Now Run the PV and Track Fitter over the collected tracks and pvs
1199  // if (theBeamFitter->runPVandTrkFitter()) {
1200  // reco::BeamSpot bs = theBeamFitter->getBeamSpot();
1201 
1202  //Create fake BS here
1204  for (int j = 0; j < 7; ++j) {
1205  for (int k = j; k < 7; ++k) {
1206  matrix(j, k) = 0;
1207  }
1208  }
1209 
1210  // random values for fake BeamSpot
1211  float tmp_BSx = rndm_->Gaus(0.1, 0.1); // [cm]
1212  float tmp_BSy = rndm_->Gaus(0.1, 0.1); // [cm]
1213  float tmp_BSz = rndm_->Gaus(0.1, 0.1); // [cm]
1214  float tmp_BSwidthX = rndm_->Gaus(0.001, 0.0005); // [cm]
1215  float tmp_BSwidthY = rndm_->Gaus(0.001, 0.0005); // [cm]
1216  float tmp_BSwidthZ = rndm_->Gaus(3.5, 0.5); // [cm]
1217 
1218  reco::BeamSpot bs(reco::BeamSpot::Point(tmp_BSx, tmp_BSy, tmp_BSz),
1219  tmp_BSwidthZ,
1220  0,
1221  0,
1222  tmp_BSwidthX,
1223  matrix,
1225  bs.setBeamWidthY(tmp_BSwidthY);
1226 
1227  if (bs.type() > 0) // with good beamwidth fit
1228  preBS = bs; // cache good fit results
1229 
1230  edm::LogInfo("FakeBeamMonitor") << "\n RESULTS OF DEFAULT FIT:" << endl;
1231  edm::LogInfo("FakeBeamMonitor") << bs << endl;
1232  edm::LogInfo("FakeBeamMonitor") << "[BeamFitter] fitting done \n" << endl;
1233 
1234  hs[k_x0_lumi]->ShiftFillLast(bs.x0(), bs.x0Error(), fitNLumi_);
1235  hs[k_y0_lumi]->ShiftFillLast(bs.y0(), bs.y0Error(), fitNLumi_);
1236  hs[k_z0_lumi]->ShiftFillLast(bs.z0(), bs.z0Error(), fitNLumi_);
1237  hs[k_sigmaX0_lumi]->ShiftFillLast(bs.BeamWidthX(), bs.BeamWidthXError(), fitNLumi_);
1238  hs[k_sigmaY0_lumi]->ShiftFillLast(bs.BeamWidthY(), bs.BeamWidthYError(), fitNLumi_);
1239  hs[k_sigmaZ0_lumi]->ShiftFillLast(bs.sigmaZ(), bs.sigmaZ0Error(), fitNLumi_);
1240  hs[k_x0_lumi_all]->setBinContent(currentlumi, bs.x0());
1241  hs[k_x0_lumi_all]->setBinError(currentlumi, bs.x0Error());
1242  hs[k_y0_lumi_all]->setBinContent(currentlumi, bs.y0());
1243  hs[k_y0_lumi_all]->setBinError(currentlumi, bs.y0Error());
1244  hs[k_z0_lumi_all]->setBinContent(currentlumi, bs.z0());
1245  hs[k_z0_lumi_all]->setBinError(currentlumi, bs.z0Error());
1246  hs[k_sigmaX0_lumi_all]->setBinContent(currentlumi, bs.BeamWidthX());
1247  hs[k_sigmaX0_lumi_all]->setBinError(currentlumi, bs.BeamWidthXError());
1248  hs[k_sigmaY0_lumi_all]->setBinContent(currentlumi, bs.BeamWidthY());
1249  hs[k_sigmaY0_lumi_all]->setBinError(currentlumi, bs.BeamWidthYError());
1250  hs[k_sigmaZ0_lumi_all]->setBinContent(currentlumi, bs.sigmaZ());
1251  hs[k_sigmaZ0_lumi_all]->setBinError(currentlumi, bs.sigmaZ0Error());
1252 
1253  int nthBin = tmpTime - refTime;
1254  if (nthBin > 0) {
1255  hs[k_x0_time]->setBinContent(nthBin, bs.x0());
1256  hs[k_y0_time]->setBinContent(nthBin, bs.y0());
1257  hs[k_z0_time]->setBinContent(nthBin, bs.z0());
1258  hs[k_sigmaX0_time]->setBinContent(nthBin, bs.BeamWidthX());
1259  hs[k_sigmaY0_time]->setBinContent(nthBin, bs.BeamWidthY());
1260  hs[k_sigmaZ0_time]->setBinContent(nthBin, bs.sigmaZ());
1261  hs[k_x0_time]->setBinError(nthBin, bs.x0Error());
1262  hs[k_y0_time]->setBinError(nthBin, bs.y0Error());
1263  hs[k_z0_time]->setBinError(nthBin, bs.z0Error());
1264  hs[k_sigmaX0_time]->setBinError(nthBin, bs.BeamWidthXError());
1265  hs[k_sigmaY0_time]->setBinError(nthBin, bs.BeamWidthYError());
1266  hs[k_sigmaZ0_time]->setBinError(nthBin, bs.sigmaZ0Error());
1267  }
1268 
1269  int jthBin = tmpTime - startTime;
1270  if (jthBin > 0) {
1271  hs[k_x0_time_all]->setBinContent(jthBin, bs.x0());
1272  hs[k_y0_time_all]->setBinContent(jthBin, bs.y0());
1273  hs[k_z0_time_all]->setBinContent(jthBin, bs.z0());
1274  hs[k_sigmaX0_time_all]->setBinContent(jthBin, bs.BeamWidthX());
1275  hs[k_sigmaY0_time_all]->setBinContent(jthBin, bs.BeamWidthY());
1276  hs[k_sigmaZ0_time_all]->setBinContent(jthBin, bs.sigmaZ());
1277  hs[k_x0_time_all]->setBinError(jthBin, bs.x0Error());
1278  hs[k_y0_time_all]->setBinError(jthBin, bs.y0Error());
1279  hs[k_z0_time_all]->setBinError(jthBin, bs.z0Error());
1280  hs[k_sigmaX0_time_all]->setBinError(jthBin, bs.BeamWidthXError());
1281  hs[k_sigmaY0_time_all]->setBinError(jthBin, bs.BeamWidthYError());
1282  hs[k_sigmaZ0_time_all]->setBinError(jthBin, bs.sigmaZ0Error());
1283  }
1284 
1285  h_x0->Fill(bs.x0());
1286  h_y0->Fill(bs.y0());
1287  h_z0->Fill(bs.z0());
1288  if (bs.type() > 0) { // with good beamwidth fit
1289  h_sigmaX0->Fill(bs.BeamWidthX());
1290  h_sigmaY0->Fill(bs.BeamWidthY());
1291  }
1292  h_sigmaZ0->Fill(bs.sigmaZ());
1293 
1294  if (nthBSTrk_ >= 2 * min_Ntrks_) {
1295  double amp = std::sqrt(bs.x0() * bs.x0() + bs.y0() * bs.y0());
1296  double alpha = std::atan2(bs.y0(), bs.x0());
1297  std::unique_ptr<TF1> f1{new TF1("f1", "[0]*sin(x-[1])", -3.14, 3.14)};
1298  f1->SetParameters(amp, alpha);
1299  f1->SetParLimits(0, amp - 0.1, amp + 0.1);
1300  f1->SetParLimits(1, alpha - 0.577, alpha + 0.577);
1301  f1->SetLineColor(4);
1302  h_d0_phi0->getTProfile()->Fit(f1.get(), "QR");
1303 
1304  double mean = bs.z0();
1305  double width = bs.sigmaZ();
1306  std::unique_ptr<TF1> fgaus{new TF1("fgaus", "gaus")};
1307  fgaus->SetParameters(mean, width);
1308  fgaus->SetLineColor(4);
1309  h_trk_z0->getTH1()->Fit(fgaus.get(), "QLRM", "", mean - 3 * width, mean + 3 * width);
1310  }
1311 
1312  fitResults->Reset();
1313  std::pair<int, int> LSRange(beginLumiOfBSFit_, endLumiOfBSFit_); //= theBeamFitter->getFitLSRange();
1314  char tmpTitle[50];
1315  sprintf(tmpTitle, "%s %i %s %i", "Fitted Beam Spot (cm) of LS: ", LSRange.first, " to ", LSRange.second);
1316  fitResults->setAxisTitle(tmpTitle, 1);
1317  fitResults->setBinContent(1, 8, bs.x0());
1318  fitResults->setBinContent(1, 7, bs.y0());
1319  fitResults->setBinContent(1, 6, bs.z0());
1320  fitResults->setBinContent(1, 5, bs.sigmaZ());
1321  fitResults->setBinContent(1, 4, bs.dxdz());
1322  fitResults->setBinContent(1, 3, bs.dydz());
1323  if (bs.type() > 0) { // with good beamwidth fit
1324  fitResults->setBinContent(1, 2, bs.BeamWidthX());
1325  fitResults->setBinContent(1, 1, bs.BeamWidthY());
1326  } else { // fill cached widths
1329  }
1330 
1331  fitResults->setBinContent(2, 8, bs.x0Error());
1332  fitResults->setBinContent(2, 7, bs.y0Error());
1333  fitResults->setBinContent(2, 6, bs.z0Error());
1334  fitResults->setBinContent(2, 5, bs.sigmaZ0Error());
1335  fitResults->setBinContent(2, 4, bs.dxdzError());
1336  fitResults->setBinContent(2, 3, bs.dydzError());
1337  if (bs.type() > 0) { // with good beamwidth fit
1338  fitResults->setBinContent(2, 2, bs.BeamWidthXError());
1339  fitResults->setBinContent(2, 1, bs.BeamWidthYError());
1340  } else { // fill cached width errors
1343  }
1344 
1345  // count good fit
1346  // if (std::fabs(refBS.x0()-bs.x0())/bs.x0Error() < deltaSigCut_) { // disabled temporarily
1347  summaryContent_[0] += 1.;
1348  // }
1349  // if (std::fabs(refBS.y0()-bs.y0())/bs.y0Error() < deltaSigCut_) { // disabled temporarily
1350  summaryContent_[1] += 1.;
1351  // }
1352  // if (std::fabs(refBS.z0()-bs.z0())/bs.z0Error() < deltaSigCut_) { // disabled temporarily
1353  summaryContent_[2] += 1.;
1354  // }
1355 
1356  // Create the BeamSpotOnlineObjects object
1357  BeamSpotOnlineObjects BSOnline;
1358  BSOnline.setLastAnalyzedLumi(LSRange.second);
1359  BSOnline.setLastAnalyzedRun(frun);
1360  BSOnline.setLastAnalyzedFill(0); // To be updated with correct LHC Fill number
1361  BSOnline.setPosition(bs.x0(), bs.y0(), bs.z0());
1362  BSOnline.setSigmaZ(bs.sigmaZ());
1363  BSOnline.setBeamWidthX(bs.BeamWidthX());
1364  BSOnline.setBeamWidthY(bs.BeamWidthY());
1365  BSOnline.setBeamWidthXError(bs.BeamWidthXError());
1366  BSOnline.setBeamWidthYError(bs.BeamWidthYError());
1367  BSOnline.setdxdz(bs.dxdz());
1368  BSOnline.setdydz(bs.dydz());
1369  BSOnline.setType(bs.type());
1370  BSOnline.setEmittanceX(bs.emittanceX());
1371  BSOnline.setEmittanceY(bs.emittanceY());
1372  BSOnline.setBetaStar(bs.betaStar());
1373  for (int i = 0; i < 7; ++i) {
1374  for (int j = 0; j < 7; ++j) {
1375  BSOnline.setCovariance(i, j, bs.covariance(i, j));
1376  }
1377  }
1378  BSOnline.setNumTracks(50);
1379  BSOnline.setNumPVs(10);
1380  BSOnline.setUsedEvents((int)DipPVInfo_[0]);
1381  BSOnline.setMeanPV(DipPVInfo_[1]);
1382  BSOnline.setMeanErrorPV(DipPVInfo_[2]);
1383  BSOnline.setRmsPV(DipPVInfo_[3]);
1384  BSOnline.setRmsErrorPV(DipPVInfo_[4]);
1385  BSOnline.setMaxPVs((int)DipPVInfo_[5]);
1386  auto creationTime =
1387  std::chrono::duration_cast<std::chrono::microseconds>(std::chrono::system_clock::now().time_since_epoch()).count();
1388  BSOnline.setCreationTime(creationTime);
1389 
1390  // use fake timestamps from 1970.01.01 00:00:00 to 1970.01.01 00:00:01 GMT
1391  std::pair<time_t, time_t> timeForDIP = std::make_pair(0, 1);
1392  BSOnline.setStartTimeStamp(timeForDIP.first);
1393  BSOnline.setStartTime(getGMTstring(timeForDIP.first));
1394  BSOnline.setEndTimeStamp(timeForDIP.second);
1395  BSOnline.setEndTime(getGMTstring(timeForDIP.second));
1396 
1397  edm::LogInfo("FakeBeamMonitor") << "FitAndFill::[PayloadCreation] BeamSpotOnline object created: \n" << std::endl;
1398  edm::LogInfo("FakeBeamMonitor") << BSOnline << std::endl;
1399 
1400  // Create the payload for BeamSpotOnlineObjects object
1401  if (onlineDbService_.isAvailable()) {
1402  edm::LogInfo("FakeBeamMonitor") << "FitAndFill::[PayloadCreation] onlineDbService available \n" << std::endl;
1403  onlineDbService_->logger().logInfo() << "FakeBeamMonitor::FitAndFill - Lumi of the current fit: " << currentlumi;
1404  onlineDbService_->logger().logInfo() << "FakeBeamMonitor::FitAndFill - Do PV Fitting for LS = " << beginLumiOfPVFit_
1405  << " to " << endLumiOfPVFit_;
1406  onlineDbService_->logger().logInfo() << "FakeBeamMonitor::FitAndFill - [BeamFitter] Do BeamSpot Fit for LS = "
1407  << LSRange.first << " to " << LSRange.second;
1408  onlineDbService_->logger().logInfo() << "FakeBeamMonitor::FitAndFill - [FakeBeamMonitor] Do BeamSpot Fit for LS = "
1409  << beginLumiOfBSFit_ << " to " << endLumiOfBSFit_;
1410  onlineDbService_->logger().logInfo() << "FakeBeamMonitor - RESULTS OF DEFAULT FIT:";
1411  onlineDbService_->logger().logInfo() << "\n" << bs;
1413  << "FakeBeamMonitor::FitAndFill - [PayloadCreation] BeamSpotOnline object created:";
1414  onlineDbService_->logger().logInfo() << "\n" << BSOnline;
1415  onlineDbService_->logger().logInfo() << "FakeBeamMonitor - Additional parameters for DIP:";
1416  onlineDbService_->logger().logInfo() << "Events used in the fit: " << BSOnline.usedEvents();
1417  onlineDbService_->logger().logInfo() << "Mean PV : " << BSOnline.meanPV();
1418  onlineDbService_->logger().logInfo() << "Mean PV Error : " << BSOnline.meanErrorPV();
1419  onlineDbService_->logger().logInfo() << "Rms PV : " << BSOnline.rmsPV();
1420  onlineDbService_->logger().logInfo() << "Rms PV Error : " << BSOnline.rmsErrorPV();
1421  onlineDbService_->logger().logInfo() << "Max PVs : " << BSOnline.maxPVs();
1422  onlineDbService_->logger().logInfo() << "StartTime : " << BSOnline.startTime();
1423  onlineDbService_->logger().logInfo() << "StartTimeStamp : " << BSOnline.startTimeStamp();
1424  onlineDbService_->logger().logInfo() << "EndTime : " << BSOnline.endTime();
1425  onlineDbService_->logger().logInfo() << "EndTimeStamp : " << BSOnline.endTimeStamp();
1426  onlineDbService_->logger().logInfo() << "FakeBeamMonitor::FitAndFill - [PayloadCreation] onlineDbService available";
1427  onlineDbService_->logger().logInfo() << "FakeBeamMonitor::FitAndFill - [PayloadCreation] SetCreationTime: "
1428  << creationTime << " [epoch in microseconds]";
1429  try {
1432  << "FakeBeamMonitor::FitAndFill - [PayloadCreation] writeIOVForNextLumisection executed correctly";
1433  } catch (const std::exception& e) {
1434  onlineDbService_->logger().logError() << "FakeBeamMonitor - Error writing record: " << recordName_
1435  << " for Run: " << frun << " - Lumi: " << LSRange.second;
1436  onlineDbService_->logger().logError() << "Error is: " << e.what();
1437  onlineDbService_->logger().logError() << "RESULTS OF DEFAULT FIT WAS:";
1438  onlineDbService_->logger().logError() << "\n" << bs;
1439  DBloggerReturn_ = 2;
1440  }
1441  }
1442  edm::LogInfo("FakeBeamMonitor") << "FitAndFill::[PayloadCreation] BeamSpotOnline payload created \n" << std::endl;
1443 
1444  // } //if (theBeamFitter->runPVandTrkFitter())
1445  // else { // beam fit fails
1446  // reco::BeamSpot bs = theBeamFitter->getBeamSpot();
1447  // edm::LogInfo("BeamMonitor") << "FitAndFill:: [BeamMonitor] Beam fit fails!!! \n" << endl;
1448  // edm::LogInfo("BeamMonitor") << "FitAndFill:: [BeamMonitor] Output beam spot for DIP \n" << endl;
1449  // edm::LogInfo("BeamMonitor") << bs << endl;
1450  //
1451  // hs[k_sigmaX0_lumi]->ShiftFillLast(bs.BeamWidthX(), bs.BeamWidthXError(), fitNLumi_);
1452  // hs[k_sigmaY0_lumi]->ShiftFillLast(bs.BeamWidthY(), bs.BeamWidthYError(), fitNLumi_);
1453  // hs[k_sigmaZ0_lumi]->ShiftFillLast(bs.sigmaZ(), bs.sigmaZ0Error(), fitNLumi_);
1454  // hs[k_x0_lumi]->ShiftFillLast(bs.x0(), bs.x0Error(), fitNLumi_);
1455  // hs[k_y0_lumi]->ShiftFillLast(bs.y0(), bs.y0Error(), fitNLumi_);
1456  // hs[k_z0_lumi]->ShiftFillLast(bs.z0(), bs.z0Error(), fitNLumi_);
1457  // } // end of beam fit fails
1458 
1459  // } //-------- end of countFitting------------------------------------------
1460  // else { // no fit
1461  // // Overwrite Fit LS and fit time when no event processed or no track selected
1462  // theBeamFitter->setFitLSRange(beginLumiOfBSFit_, endLumiOfBSFit_);
1463  // theBeamFitter->setRefTime(refBStime[0], refBStime[1]);
1464  // if (theBeamFitter->runPVandTrkFitter()) {
1465  // } // Dump fake beam spot for DIP
1466  // reco::BeamSpot bs = theBeamFitter->getBeamSpot();
1467  // edm::LogInfo("BeamMonitor") << "FitAndFill:: [BeamMonitor] No fitting \n" << endl;
1468  // edm::LogInfo("BeamMonitor") << "FitAndFill:: [BeamMonitor] Output fake beam spot for DIP \n" << endl;
1469  // edm::LogInfo("BeamMonitor") << bs << endl;
1470  //
1471  // hs[k_sigmaX0_lumi]->ShiftFillLast(bs.BeamWidthX(), bs.BeamWidthXError(), fitNLumi_);
1472  // hs[k_sigmaY0_lumi]->ShiftFillLast(bs.BeamWidthY(), bs.BeamWidthYError(), fitNLumi_);
1473  // hs[k_sigmaZ0_lumi]->ShiftFillLast(bs.sigmaZ(), bs.sigmaZ0Error(), fitNLumi_);
1474  // hs[k_x0_lumi]->ShiftFillLast(bs.x0(), bs.x0Error(), fitNLumi_);
1475  // hs[k_y0_lumi]->ShiftFillLast(bs.y0(), bs.y0Error(), fitNLumi_);
1476  // hs[k_z0_lumi]->ShiftFillLast(bs.z0(), bs.z0Error(), fitNLumi_);
1477  // }
1478 
1479  // Fill summary report
1480  // if (countFitting) {
1481  for (int n = 0; n < nFitElements_; n++) {
1483  }
1484 
1485  summarySum_ = 0;
1486  for (int ii = 0; ii < nFitElements_; ii++) {
1488  }
1490  if (reportSummary)
1492 
1493  for (int bi = 0; bi < nFitElements_; bi++) {
1494  reportSummaryMap->setBinContent(1, bi + 1, summaryContent_[bi] / (float)nFits_);
1495  }
1496  // }
1497 
1498  if ((resetFitNLumi_ > 0 &&
1499  ((onlineMode_ &&
1500  countLumi_ == resetFitNLumi_) || //OR it should be currentLumi_ (if in sequence then does not mattar)
1501  (!onlineMode_ && countLumi_ == resetFitNLumi_))) ||
1502  (StartAverage_)) {
1503  edm::LogInfo("FakeBeamMonitor") << "FitAndFill:: The flag is ON for running average Beam Spot fit" << endl;
1504  StartAverage_ = true;
1505  firstAverageFit_++;
1506  resetHistos_ = true;
1507  nthBSTrk_ = 0;
1508  beginLumiOfBSFit_ = 0;
1509  refBStime[0] = 0;
1510  }
1511 }
MonitorElement * pvResults
math::Error< dimension >::type CovarianceMatrix
Definition: BeamSpot.h:29
double BeamWidthX() const
beam width X
Definition: BeamSpot.h:82
void setBeamWidthYError(double val)
set beam width Y error
void setCovariance(int i, int j, double val)
set i,j element of the full covariance matrix 7x7
std::map< int, size_t > mapLSPVStoreSize
void setStartTimeStamp(cond::Time_t val)
MonitorElement * h_z0
std::vector< MonitorElement * > hs
double BeamWidthYError() const
error on beam width Y, assume error in X = Y
Definition: BeamSpot.h:101
MonitorElement * h_nTrk_lumi
void setEmittanceY(double val)
set emittance
MonitorElement * h_sigmaY0
MonitorElement * h_vy_dz
std::string getGMTstring(const std::time_t &)
void setBetaStar(double val)
set beta star
unsigned int nthBSTrk_
bool testScroll(std::time_t &, std::time_t &)
std::string endTime() const
void setEndTimeStamp(cond::Time_t val)
math::XYZPoint Point
point in the space
Definition: BeamSpot.h:27
MonitorElement * reportSummary
Float_t summaryContent_[3]
std::time_t refBStime[2]
std::string startTime() const
std::map< int, std::vector< int > > mapNPV
void setType(int type)
set beam type
MonitorElement * h_PVz[2]
void setEndTime(std::string val)
unsigned int minNrVertices_
edm::Service< cond::service::OnlineDBOutputService > onlineDbService_
void Fill(long long x)
virtual void Reset()
Remove all data from the ME, keept the empty histogram with all its settings.
DQM_DEPRECATED void ShiftFillLast(double y, double ye=0., int32_t xscale=1)
void setEmittanceX(double val)
set emittance
EchoedLogStream< edm::LogError > logError()
Definition: Logger.cc:165
std::time_t refPVtime[2]
cond::persistency::Logger & logger()
T sqrt(T t)
Definition: SSEVec.h:19
std::time_t startTime
MonitorElement * h_sigmaZ0
void setdydz(double val)
set dydz slope, crossing angle in XZ
cond::Time_t endTimeStamp() const
MonitorElement * h_x0
double BeamWidthY() const
beam width Y
Definition: BeamSpot.h:84
std::map< int, std::vector< float > > mapPVy
std::map< int, std::vector< float > > mapPVz
MonitorElement * h_d0_phi0
cond::Time_t writeIOVForNextLumisection(const PayloadType &payload, const std::string &recordName)
std::map< int, std::vector< float > > mapPVx
virtual TProfile * getTProfile()
void setdxdz(double val)
set dxdz slope, crossing angle
ii
Definition: cuy.py:589
MonitorElement * h_sigmaX0
double BeamWidthXError() const
error on beam width X, assume error in X = Y
Definition: BeamSpot.h:99
Log< level::Info, false > LogInfo
cond::Time_t startTimeStamp() const
void scrollTH1(TH1 *, std::time_t)
virtual void setBinContent(int binx, double content)
set content of bin (1-D)
void setBeamWidthXError(double val)
set beam width X error
MonitorElement * h_trk_z0
std::time_t refTime
void setSigmaZ(double val)
set sigma Z, RMS bunch length
MonitorElement * reportSummaryContents[3]
MonitorElement * h_nVtx_lumi
MonitorElement * h_vx_vy
MonitorElement * h_vx_dz
MonitorElement * h_PVy[2]
std::string recordName_
MonitorElement * h_y0
reco::BeamSpot preBS
MonitorElement * h_nVtx_st
void setStartTime(std::string val)
void setPosition(double x, double y, double z)
set XYZ position
void setBeamWidthX(double val)
set average transverse beam width X
std::time_t tmpTime
void setCreationTime(cond::Time_t val)
MonitorElement * reportSummaryMap
bool isAvailable() const
Definition: Service.h:40
MonitorElement * h_PVx[2]
EchoedLogStream< edm::LogInfo > logInfo()
Definition: Logger.cc:157
unsigned int min_Ntrks_
void setBeamWidthY(double val)
set average transverse beam width Y
MonitorElement * fitResults
virtual void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)

◆ formatFitTime()

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

Definition at line 44 of file FakeBeamMonitor.cc.

References b, CEST, plotBeamSpotDB::ptm, submitPVValidationJobs::t, hcalRecHitTable_cff::time, and HLT_2023v12_cff::year.

Referenced by bookHistograms(), and scrollTH1().

44  {
45  //constexpr int CET(+1);
46  constexpr int CEST(+2);
47 
48  //tm * ptm;
49  //ptm = gmtime ( &t );
50  //int year = ptm->tm_year;
51 
52  //get correct year from ctime
53  time_t currentTime;
54  struct tm* localTime;
55  time(&currentTime); // Get the current time
56  localTime = localtime(&currentTime); // Convert the current time to the local time
57  int year = localTime->tm_year + 1900;
58 
59  tm* ptm;
60  ptm = gmtime(&t);
61 
62  //check if year is ok
63  if (year <= 37)
64  year += 2000;
65  if (year >= 70 && year <= 137)
66  year += 1900;
67 
68  if (year < 1995) {
69  edm::LogError("BadTimeStamp") << "year reported is " << year << " !!" << std::endl;
70  //year = 2015; //overwritten later by BeamFitter.cc for fits but needed here for TH1
71  //edm::LogError("BadTimeStamp") << "Resetting to " <<year<<std::endl;
72  }
73  sprintf(ts,
74  "%4d-%02d-%02d %02d:%02d:%02d",
75  year,
76  ptm->tm_mon + 1,
77  ptm->tm_mday,
78  (ptm->tm_hour + CEST) % 24,
79  ptm->tm_min,
80  ptm->tm_sec);
81 
82 #ifdef STRIP_TRAILING_BLANKS_IN_TIMEZONE
83  unsigned int b = strlen(ts);
84  while (ts[--b] == ' ') {
85  ts[b] = 0;
86  }
87 #endif
88 }
#define CEST
Log< level::Error, false > LogError
double b
Definition: hdecay.h:120

◆ getGMTstring()

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

Definition at line 92 of file FakeBeamMonitor.cc.

References AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by FitAndFill().

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

◆ RestartFitting()

void FakeBeamMonitor::RestartFitting ( )
private

Definition at line 1514 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_.

1514  { //
1515  if (debug_)
1516  edm::LogInfo("FakeBeamMonitor")
1517  << " RestartingFitting:: Restart Beami everything to a fresh start !!! because Gap is > 10 LS" << endl;
1518  //track based fit reset here
1519  resetHistos_ = true;
1520  nthBSTrk_ = 0;
1521  // theBeamFitter->resetTrkVector();
1522  // theBeamFitter->resetLSRange();
1523  // theBeamFitter->resetRefTime();
1524  // theBeamFitter->resetPVFitter();
1525  // theBeamFitter->resetCutFlow();
1526  beginLumiOfBSFit_ = 0;
1527  refBStime[0] = 0;
1528  //pv based fit iis reset here
1529  h_PVx[0]->Reset();
1530  h_PVy[0]->Reset();
1531  h_PVz[0]->Reset();
1532  beginLumiOfPVFit_ = 0;
1533  refPVtime[0] = 0;
1534  //Clear all the Maps here
1535  mapPVx.clear();
1536  mapPVy.clear();
1537  mapPVz.clear();
1538  mapNPV.clear();
1539  mapBeginBSLS.clear();
1540  mapBeginPVLS.clear();
1541  mapBeginBSTime.clear();
1542  mapBeginPVTime.clear();
1543  mapLSBSTrkSize.clear();
1544  mapLSPVStoreSize.clear();
1545  mapLSCF.clear();
1546  countGapLumi_ = 0;
1547  countLumi_ = 0;
1548  StartAverage_ = false;
1549 }
std::map< int, size_t > mapLSPVStoreSize
std::map< int, std::time_t > mapBeginPVTime
unsigned int nthBSTrk_
std::map< int, std::size_t > mapLSBSTrkSize
std::time_t refBStime[2]
std::map< int, std::vector< int > > mapNPV
MonitorElement * h_PVz[2]
virtual void Reset()
Remove all data from the ME, keept the empty histogram with all its settings.
std::time_t refPVtime[2]
std::map< int, std::vector< float > > mapPVy
std::map< int, std::vector< float > > mapPVz
std::map< int, std::vector< float > > mapPVx
std::map< int, std::time_t > mapBeginBSTime
Log< level::Info, false > LogInfo
std::map< int, int > mapBeginPVLS
std::map< int, TH1F > mapLSCF
MonitorElement * h_PVy[2]
std::map< int, int > mapBeginBSLS
MonitorElement * h_PVx[2]

◆ scrollTH1()

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

Definition at line 1573 of file FakeBeamMonitor.cc.

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

Referenced by FitAndFill().

1573  {
1574  char offsetTime[64];
1575  formatFitTime(offsetTime, ref);
1576  TDatime da(offsetTime);
1577  if (lastNZbin > 0) {
1578  double val = h->GetBinContent(lastNZbin);
1579  double valErr = h->GetBinError(lastNZbin);
1580  h->Reset();
1581  h->GetXaxis()->SetTimeOffset(da.Convert(kTRUE));
1582  int bin = (lastNZbin > buffTime ? buffTime : 1);
1583  h->SetBinContent(bin, val);
1584  h->SetBinError(bin, valErr);
1585  } else {
1586  h->Reset();
1587  h->GetXaxis()->SetTimeOffset(da.Convert(kTRUE));
1588  }
1589 }
void formatFitTime(char *, const std::time_t &)
The Signals That Services Can Subscribe To This is based on ActivityRegistry h
Helper function to determine trigger accepts.
Definition: Activities.doc:4
static constexpr int buffTime

◆ testScroll()

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

Definition at line 1593 of file FakeBeamMonitor.cc.

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

Referenced by FitAndFill().

1593  {
1594  bool scroll_ = false;
1595  if (tmpTime_ - refTime_ >= intervalInSec_) {
1596  scroll_ = true;
1597  edm::LogInfo("FakeBeamMonitor") << "testScroll:: Reset Time Offset" << std::endl;
1599  for (int bin = intervalInSec_; bin >= 1; bin--) {
1600  if (hs[k_x0_time]->getBinContent(bin) > 0) {
1601  lastNZbin = bin;
1602  break;
1603  }
1604  }
1605  edm::LogInfo("FakeBeamMonitor") << "testScroll:: Last non zero bin = " << lastNZbin << std::endl;
1606  if (tmpTime_ - refTime_ >= intervalInSec_ + lastNZbin) {
1607  edm::LogInfo("FakeBeamMonitor") << "testScroll:: Time difference too large since last readout" << std::endl;
1608  lastNZbin = 0;
1609  refTime_ = tmpTime_ - buffTime;
1610  } else {
1611  edm::LogInfo("FakeBeamMonitor") << "testScroll:: Offset to last record" << std::endl;
1612  int offset = ((lastNZbin > buffTime) ? (lastNZbin - buffTime) : (lastNZbin - 1));
1613  refTime_ += offset;
1614  }
1615  }
1616  return scroll_;
1617 }
std::vector< MonitorElement * > hs
Log< level::Info, false > LogInfo
static constexpr int buffTime

Member Data Documentation

◆ beginLumiOfBSFit_

int FakeBeamMonitor::beginLumiOfBSFit_
private

◆ beginLumiOfPVFit_

int FakeBeamMonitor::beginLumiOfPVFit_
private

◆ countEvt_

int FakeBeamMonitor::countEvt_
private

Definition at line 89 of file FakeBeamMonitor.h.

Referenced by analyze().

◆ countGapLumi_

int FakeBeamMonitor::countGapLumi_
private

Definition at line 112 of file FakeBeamMonitor.h.

Referenced by RestartFitting().

◆ countLumi_

int FakeBeamMonitor::countLumi_
private

Definition at line 90 of file FakeBeamMonitor.h.

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

◆ cutFlowTable

MonitorElement* FakeBeamMonitor::cutFlowTable
private

Definition at line 169 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ DBloggerReturn_

int FakeBeamMonitor::DBloggerReturn_
private

Definition at line 79 of file FakeBeamMonitor.h.

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

◆ debug_

bool FakeBeamMonitor::debug_
private

◆ deltaSigCut_

double FakeBeamMonitor::deltaSigCut_
private

Definition at line 102 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor().

◆ dxBin_

const int FakeBeamMonitor::dxBin_
private

Definition at line 61 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ dxMax_

const double FakeBeamMonitor::dxMax_
private

Definition at line 63 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ dxMin_

const double FakeBeamMonitor::dxMin_
private

Definition at line 62 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ dzBin_

const int FakeBeamMonitor::dzBin_
private

Definition at line 73 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ dzMax_

const double FakeBeamMonitor::dzMax_
private

Definition at line 75 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ dzMin_

const double FakeBeamMonitor::dzMin_
private

Definition at line 74 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ endLumiOfBSFit_

int FakeBeamMonitor::endLumiOfBSFit_
private

Definition at line 92 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

◆ endLumiOfPVFit_

int FakeBeamMonitor::endLumiOfPVFit_
private

Definition at line 94 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

◆ firstAverageFit_

int FakeBeamMonitor::firstAverageFit_
private

Definition at line 111 of file FakeBeamMonitor.h.

Referenced by FitAndFill().

◆ fitNLumi_

int FakeBeamMonitor::fitNLumi_
private

Definition at line 81 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

◆ fitPVNLumi_

int FakeBeamMonitor::fitPVNLumi_
private

Definition at line 82 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

◆ fitResults

MonitorElement* FakeBeamMonitor::fitResults
private

Definition at line 133 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ frun

int FakeBeamMonitor::frun
private

Definition at line 178 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ ftimestamp

edm::TimeValue_t FakeBeamMonitor::ftimestamp
private

Definition at line 176 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ h_d0_phi0

MonitorElement* FakeBeamMonitor::h_d0_phi0
private

Definition at line 126 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ h_nTrk_lumi

MonitorElement* FakeBeamMonitor::h_nTrk_lumi
private

Definition at line 123 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ h_nVtx

MonitorElement* FakeBeamMonitor::h_nVtx
private

Definition at line 140 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ h_nVtx_lumi

MonitorElement* FakeBeamMonitor::h_nVtx_lumi
private

Definition at line 124 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ h_nVtx_lumi_all

MonitorElement* FakeBeamMonitor::h_nVtx_lumi_all
private

Definition at line 125 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ h_nVtx_st

MonitorElement* FakeBeamMonitor::h_nVtx_st
private

Definition at line 141 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ h_PVx

MonitorElement* FakeBeamMonitor::h_PVx[2]
private

Definition at line 142 of file FakeBeamMonitor.h.

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

◆ h_PVxz

MonitorElement* FakeBeamMonitor::h_PVxz
private

Definition at line 145 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ h_PVy

MonitorElement* FakeBeamMonitor::h_PVy[2]
private

Definition at line 143 of file FakeBeamMonitor.h.

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

◆ h_PVyz

MonitorElement* FakeBeamMonitor::h_PVyz
private

Definition at line 146 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ h_PVz

MonitorElement* FakeBeamMonitor::h_PVz[2]
private

Definition at line 144 of file FakeBeamMonitor.h.

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

◆ h_sigmaX0

MonitorElement* FakeBeamMonitor::h_sigmaX0
private

Definition at line 137 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ h_sigmaY0

MonitorElement* FakeBeamMonitor::h_sigmaY0
private

Definition at line 138 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ h_sigmaZ0

MonitorElement* FakeBeamMonitor::h_sigmaZ0
private

Definition at line 139 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ h_trk_z0

MonitorElement* FakeBeamMonitor::h_trk_z0
private

Definition at line 127 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ h_trkPt

MonitorElement* FakeBeamMonitor::h_trkPt
private

Definition at line 131 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ h_trkVz

MonitorElement* FakeBeamMonitor::h_trkVz
private

Definition at line 132 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ h_vx_dz

MonitorElement* FakeBeamMonitor::h_vx_dz
private

Definition at line 129 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ h_vx_vy

MonitorElement* FakeBeamMonitor::h_vx_vy
private

Definition at line 128 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ h_vy_dz

MonitorElement* FakeBeamMonitor::h_vy_dz
private

Definition at line 130 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ h_x0

MonitorElement* FakeBeamMonitor::h_x0
private

Definition at line 134 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ h_y0

MonitorElement* FakeBeamMonitor::h_y0
private

Definition at line 135 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ h_z0

MonitorElement* FakeBeamMonitor::h_z0
private

Definition at line 136 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ hs

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

Definition at line 148 of file FakeBeamMonitor.h.

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

◆ intervalInSec_

int FakeBeamMonitor::intervalInSec_
private

Definition at line 85 of file FakeBeamMonitor.h.

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

◆ lastlumi_

int FakeBeamMonitor::lastlumi_
private

Definition at line 95 of file FakeBeamMonitor.h.

Referenced by beginLuminosityBlock(), and FakeBeamMonitor().

◆ lastNZbin

int FakeBeamMonitor::lastNZbin
private

Definition at line 179 of file FakeBeamMonitor.h.

Referenced by scrollTH1(), and testScroll().

◆ mapBeginBSLS

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

Definition at line 154 of file FakeBeamMonitor.h.

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

◆ mapBeginBSTime

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

Definition at line 155 of file FakeBeamMonitor.h.

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

◆ mapBeginPVLS

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

Definition at line 154 of file FakeBeamMonitor.h.

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

◆ mapBeginPVTime

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

Definition at line 155 of file FakeBeamMonitor.h.

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

◆ mapLSBSTrkSize

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

Definition at line 157 of file FakeBeamMonitor.h.

Referenced by dqmEndRun(), and RestartFitting().

◆ mapLSCF

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

Definition at line 160 of file FakeBeamMonitor.h.

Referenced by dqmEndRun(), and RestartFitting().

◆ mapLSPVStoreSize

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

Definition at line 158 of file FakeBeamMonitor.h.

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

◆ mapNPV

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

Definition at line 152 of file FakeBeamMonitor.h.

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

◆ mapPVx

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

Definition at line 151 of file FakeBeamMonitor.h.

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

◆ mapPVy

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

Definition at line 151 of file FakeBeamMonitor.h.

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

◆ mapPVz

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

Definition at line 151 of file FakeBeamMonitor.h.

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

◆ maxZ_

double FakeBeamMonitor::maxZ_
private

Definition at line 104 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor().

◆ min_Ntrks_

unsigned int FakeBeamMonitor::min_Ntrks_
private

Definition at line 103 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

◆ minNrVertices_

unsigned int FakeBeamMonitor::minNrVertices_
private

Definition at line 105 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

◆ minVtxNdf_

double FakeBeamMonitor::minVtxNdf_
private

Definition at line 106 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor().

◆ minVtxWgt_

double FakeBeamMonitor::minVtxWgt_
private

Definition at line 107 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor().

◆ monitorName_

std::string FakeBeamMonitor::monitorName_
private

Definition at line 76 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FakeBeamMonitor().

◆ nextlumi_

int FakeBeamMonitor::nextlumi_
private

◆ nFitElements_

int FakeBeamMonitor::nFitElements_
private

Definition at line 100 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ nFits_

int FakeBeamMonitor::nFits_
private

Definition at line 101 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

◆ nthBSTrk_

unsigned int FakeBeamMonitor::nthBSTrk_
private

Definition at line 99 of file FakeBeamMonitor.h.

Referenced by FitAndFill(), and RestartFitting().

◆ onlineDbService_

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

◆ onlineMode_

bool FakeBeamMonitor::onlineMode_
private

◆ phiBin_

const int FakeBeamMonitor::phiBin_
private

Definition at line 69 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ phiMax_

const double FakeBeamMonitor::phiMax_
private

Definition at line 71 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ phiMin_

const double FakeBeamMonitor::phiMin_
private

Definition at line 70 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ preBS

reco::BeamSpot FakeBeamMonitor::preBS
private

Definition at line 120 of file FakeBeamMonitor.h.

Referenced by FitAndFill().

◆ processed_

bool FakeBeamMonitor::processed_
private

Definition at line 114 of file FakeBeamMonitor.h.

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

◆ pvResults

MonitorElement* FakeBeamMonitor::pvResults
private

Definition at line 147 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ recordName_

std::string FakeBeamMonitor::recordName_
private

Definition at line 77 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

◆ refBS

reco::BeamSpot FakeBeamMonitor::refBS
private

Definition at line 119 of file FakeBeamMonitor.h.

◆ refBStime

std::time_t FakeBeamMonitor::refBStime[2]
private

◆ refPVtime

std::time_t FakeBeamMonitor::refPVtime[2]
private

◆ refTime

std::time_t FakeBeamMonitor::refTime
private

Definition at line 175 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ reportSummary

MonitorElement* FakeBeamMonitor::reportSummary
private

Definition at line 166 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ reportSummary_

Float_t FakeBeamMonitor::reportSummary_
private

Definition at line 163 of file FakeBeamMonitor.h.

Referenced by FitAndFill().

◆ reportSummaryContents

MonitorElement* FakeBeamMonitor::reportSummaryContents[3]
private

Definition at line 167 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ reportSummaryMap

MonitorElement* FakeBeamMonitor::reportSummaryMap
private

Definition at line 168 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ resetFitNLumi_

int FakeBeamMonitor::resetFitNLumi_
private

Definition at line 83 of file FakeBeamMonitor.h.

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

◆ resetHistos_

bool FakeBeamMonitor::resetHistos_
private

Definition at line 109 of file FakeBeamMonitor.h.

Referenced by FitAndFill(), and RestartFitting().

◆ resetPVNLumi_

int FakeBeamMonitor::resetPVNLumi_
private

Definition at line 84 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

◆ rndm_

TRandom3* FakeBeamMonitor::rndm_
private

Definition at line 181 of file FakeBeamMonitor.h.

Referenced by FakeBeamMonitor(), and FitAndFill().

◆ StartAverage_

bool FakeBeamMonitor::StartAverage_
private

Definition at line 110 of file FakeBeamMonitor.h.

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

◆ startTime

std::time_t FakeBeamMonitor::startTime
private

Definition at line 174 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ summaryContent_

Float_t FakeBeamMonitor::summaryContent_[3]
private

Definition at line 165 of file FakeBeamMonitor.h.

Referenced by bookHistograms(), and FitAndFill().

◆ summarySum_

Float_t FakeBeamMonitor::summarySum_
private

Definition at line 164 of file FakeBeamMonitor.h.

Referenced by FitAndFill().

◆ tmpTime

std::time_t FakeBeamMonitor::tmpTime
private

Definition at line 173 of file FakeBeamMonitor.h.

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

◆ useLockRecords_

bool FakeBeamMonitor::useLockRecords_
private

Definition at line 115 of file FakeBeamMonitor.h.

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

◆ vxBin_

const int FakeBeamMonitor::vxBin_
private

Definition at line 65 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ vxMax_

const double FakeBeamMonitor::vxMax_
private

Definition at line 67 of file FakeBeamMonitor.h.

Referenced by bookHistograms().

◆ vxMin_

const double FakeBeamMonitor::vxMin_
private

Definition at line 66 of file FakeBeamMonitor.h.

Referenced by bookHistograms().