CMS 3D CMS Logo
List of all members | Public Member Functions | Protected Member Functions | Private Types | Private Attributes
AlcaBeamMonitor Class Reference

#include <AlcaBeamMonitor.h>

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

Public Member Functions

 AlcaBeamMonitor (const edm::ParameterSet &)
 
 ~AlcaBeamMonitor () override
 
- Public Member Functions inherited from DQMOneLumiEDAnalyzer<>
void beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
 DQMOneLumiEDAnalyzer ()
 
void endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
 
void endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
 
bool getCanSaveByLumi () override
 
- Public Member Functions inherited from DQMOneEDAnalyzer< edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, Args... >
void accumulate (edm::Event const &event, edm::EventSetup const &setup) final
 
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
 
- Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... >
 EDProducer ()=default
 
SerialTaskQueue * globalLuminosityBlocksQueue () final
 
SerialTaskQueue * globalRunsQueue () final
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDProducerBase () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
 ~ProducerBase () noexcept(false) override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Protected Member Functions

void analyze (const edm::Event &iEvent, const edm::EventSetup &iSetup) override
 
void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 
void dqmBeginLuminosityBlock (const edm::LuminosityBlock &iLumi, const edm::EventSetup &iSetup) override
 
void dqmEndLuminosityBlock (const edm::LuminosityBlock &iLumi, const edm::EventSetup &iSetup) override
 
- Protected Member Functions inherited from DQMOneLumiEDAnalyzer<>
virtual void dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
 
virtual void dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
 
- Protected Member Functions inherited from DQMOneEDAnalyzer< edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, Args... >
virtual void analyze (edm::Event const &, edm::EventSetup const &)
 
virtual void dqmBeginRun (edm::Run const &, edm::EventSetup const &)
 
virtual void dqmEndRun (edm::Run const &, edm::EventSetup const &)
 
- Protected Member Functions inherited from edm::ProducerBase
ProducesCollector producesCollector ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename 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)
 

Private Types

typedef std::map< std::string, reco::BeamSpot > BeamSpotContainer
 
typedef std::map< std::string, std::map< std::string, std::map< std::string, MonitorElement * > > > HistosContainer
 
typedef std::map< std::string, std::map< std::string, std::map< std::string, int > > > PositionContainer
 

Private Attributes

edm::InputTag beamSpotLabel_
 
BeamSpotContainer beamSpotsMap_
 
MonitorElement * hD0Phi0_
 
MonitorElement * hDxyBS_
 
std::multimap< std::string, std::string > histoByCategoryNames_
 
HistosContainer histosMap_
 
std::string monitorName_
 
int numberOfValuesToSave_
 
edm::ParameterSet parameters_
 
PositionContainer positionsMap_
 
edm::EDGetTokenT< reco::VertexCollection > primaryVertexLabel_
 
edm::EDGetTokenT< reco::BeamSpot > scalerLabel_
 
BeamFitter * theBeamFitter_
 
PVFitter * thePVFitter_
 
MonitorElement * theValuesContainer_
 
edm::EDGetTokenT< reco::TrackCollection > trackLabel_
 
std::vector< std::string > varNamesV_
 
std::vector< reco::VertexCollection > vertices_
 

Additional Inherited Members

- Public Types inherited from DQMOneEDAnalyzer< edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, Args... >
typedef dqm::reco::DQMStore DQMStore
 
typedef dqm::reco::MonitorElement MonitorElement
 
- Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase ModuleType
 
- Public Types inherited from edm::ProducerBase
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Attributes inherited from DQMOneLumiEDAnalyzer<>
edm::EDPutTokenT< DQMToken > lumiToken_
 
- Protected Attributes inherited from DQMOneEDAnalyzer< edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, Args... >
edm::EDPutTokenT< DQMToken > runToken_
 

Detailed Description

Definition at line 27 of file AlcaBeamMonitor.h.

Member Typedef Documentation

typedef std::map<std::string, reco::BeamSpot> AlcaBeamMonitor::BeamSpotContainer
private

Definition at line 41 of file AlcaBeamMonitor.h.

typedef std::map<std::string, std::map<std::string, std::map<std::string, MonitorElement*> > > AlcaBeamMonitor::HistosContainer
private

Definition at line 43 of file AlcaBeamMonitor.h.

typedef std::map<std::string, std::map<std::string, std::map<std::string, int> > > AlcaBeamMonitor::PositionContainer
private

Definition at line 45 of file AlcaBeamMonitor.h.

Constructor & Destructor Documentation

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

Definition at line 32 of file AlcaBeamMonitor.cc.

References AlcaBeamMonitor_cfi::BeamFitter, edm::EDConsumerBase::consumesCollector(), histoByCategoryNames_, histosMap_, monitorName_, numberOfValuesToSave_, parameters_, positionsMap_, AlcaBeamMonitor_cfi::PVFitter, BeamFitter::resetLSRange(), BeamFitter::resetPVFitter(), BeamFitter::resetRefTime(), BeamFitter::resetTrkVector(), theBeamFitter_, thePVFitter_, and varNamesV_.

33  : parameters_(ps),
34  monitorName_(parameters_.getUntrackedParameter<string>("MonitorName", "YourSubsystemName")),
35  primaryVertexLabel_(
36  consumes<VertexCollection>(parameters_.getUntrackedParameter<InputTag>("PrimaryVertexLabel"))),
37  trackLabel_(consumes<reco::TrackCollection>(parameters_.getUntrackedParameter<InputTag>("TrackLabel"))),
38  scalerLabel_(consumes<BeamSpot>(parameters_.getUntrackedParameter<InputTag>("ScalerLabel"))),
39  beamSpotLabel_(parameters_.getUntrackedParameter<InputTag>("BeamSpotLabel")),
40  numberOfValuesToSave_(0) {
41  if (!monitorName_.empty())
42  monitorName_ = monitorName_ + "/";
43 
44  theBeamFitter_ = new BeamFitter(parameters_, consumesCollector());
45  theBeamFitter_->resetTrkVector();
46  theBeamFitter_->resetLSRange();
47  theBeamFitter_->resetRefTime();
48  theBeamFitter_->resetPVFitter();
49 
50  thePVFitter_ = new PVFitter(parameters_, consumesCollector());
51 
52  varNamesV_.push_back("x");
53  varNamesV_.push_back("y");
54  varNamesV_.push_back("z");
55  varNamesV_.push_back("sigmaX");
56  varNamesV_.push_back("sigmaY");
57  varNamesV_.push_back("sigmaZ");
58 
59  histoByCategoryNames_.insert(pair<string, string>("run", "Coordinate"));
60  histoByCategoryNames_.insert(pair<string, string>("run", "PrimaryVertex fit-DataBase"));
61  histoByCategoryNames_.insert(pair<string, string>("run", "PrimaryVertex fit-BeamFit"));
62  histoByCategoryNames_.insert(pair<string, string>("run", "PrimaryVertex fit-Scalers"));
63  histoByCategoryNames_.insert(pair<string, string>("run", "PrimaryVertex-DataBase"));
64  histoByCategoryNames_.insert(pair<string, string>("run", "PrimaryVertex-BeamFit"));
65  histoByCategoryNames_.insert(pair<string, string>("run", "PrimaryVertex-Scalers"));
66 
67  histoByCategoryNames_.insert(pair<string, string>("lumi", "Lumibased BeamSpotFit"));
68  histoByCategoryNames_.insert(pair<string, string>("lumi", "Lumibased PrimaryVertex"));
69  histoByCategoryNames_.insert(pair<string, string>("lumi", "Lumibased DataBase"));
70  histoByCategoryNames_.insert(pair<string, string>("lumi", "Lumibased Scalers"));
71  histoByCategoryNames_.insert(pair<string, string>("lumi", "Lumibased PrimaryVertex-DataBase fit"));
72  histoByCategoryNames_.insert(pair<string, string>("lumi", "Lumibased PrimaryVertex-Scalers fit"));
73  histoByCategoryNames_.insert(pair<string, string>("validation", "Lumibased Scalers-DataBase fit"));
74  histoByCategoryNames_.insert(pair<string, string>("validation", "Lumibased PrimaryVertex-DataBase"));
75  histoByCategoryNames_.insert(pair<string, string>("validation", "Lumibased PrimaryVertex-Scalers"));
76 
77  for (vector<string>::iterator itV = varNamesV_.begin(); itV != varNamesV_.end(); itV++) {
78  for (multimap<string, string>::iterator itM = histoByCategoryNames_.begin(); itM != histoByCategoryNames_.end();
79  itM++) {
80  if (itM->first == "run") {
81  histosMap_[*itV][itM->first][itM->second] = nullptr;
82  } else {
83  positionsMap_[*itV][itM->first][itM->second] = 3 * numberOfValuesToSave_; //value, error, ok
84  ++numberOfValuesToSave_;
85  }
86  }
87  }
88 
89  // beamSpotsMap_["BF"] = map<LuminosityBlockNumber_t,BeamSpot>();//For each lumi the beamfitter will have a result
90  // beamSpotsMap_["PV"] = map<LuminosityBlockNumber_t,BeamSpot>();//For each lumi the PVfitter will have a result
91  // beamSpotsMap_["DB"] = map<LuminosityBlockNumber_t,BeamSpot>();//For each lumi we take the values that are stored in the database, already collapsed then
92  // beamSpotsMap_["SC"] = map<LuminosityBlockNumber_t,BeamSpot>();//For each lumi we take the beamspot value in the file that is the same as the scaler for the alca reco stream
93  // beamSpotsMap_["BF"] = 0;//For each lumi the beamfitter will have a result
94  // beamSpotsMap_["PV"] = 0;//For each lumi the PVfitter will have a result
95  // beamSpotsMap_["DB"] = 0;//For each lumi we take the values that are stored in the database, already collapsed then
96  // beamSpotsMap_["SC"] = 0;//For each lumi we take the beamspot value in the file that is the same as the scaler for the alca reco stream
97 }
AlcaBeamMonitor_cfi.PVFitter
PVFitter
Definition: AlcaBeamMonitor_cfi.py:39
AlcaBeamMonitor_cfi.BeamFitter
BeamFitter
Definition: AlcaBeamMonitor_cfi.py:11
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
AlcaBeamMonitor::trackLabel_
edm::EDGetTokenT< reco::TrackCollection > trackLabel_
Definition: AlcaBeamMonitor.h:51
AlcaBeamMonitor::histoByCategoryNames_
std::multimap< std::string, std::string > histoByCategoryNames_
Definition: AlcaBeamMonitor.h:70
AlcaBeamMonitor::primaryVertexLabel_
edm::EDGetTokenT< reco::VertexCollection > primaryVertexLabel_
Definition: AlcaBeamMonitor.h:50
BeamFitter::resetTrkVector
void resetTrkVector()
Definition: BeamFitter.h:55
AlcaBeamMonitor::positionsMap_
PositionContainer positionsMap_
Definition: AlcaBeamMonitor.h:68
AlcaBeamMonitor::beamSpotLabel_
edm::InputTag beamSpotLabel_
Definition: AlcaBeamMonitor.h:53
AlcaBeamMonitor::varNamesV_
std::vector< std::string > varNamesV_
Definition: AlcaBeamMonitor.h:69
AlcaBeamMonitor::scalerLabel_
edm::EDGetTokenT< reco::BeamSpot > scalerLabel_
Definition: AlcaBeamMonitor.h:52
edm::EDConsumerBase::consumesCollector
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
Definition: EDConsumerBase.cc:71
AlcaBeamMonitor::theBeamFitter_
BeamFitter * theBeamFitter_
Definition: AlcaBeamMonitor.h:57
BeamFitter::resetLSRange
void resetLSRange()
Definition: BeamFitter.h:57
AlcaBeamMonitor::thePVFitter_
PVFitter * thePVFitter_
Definition: AlcaBeamMonitor.h:58
AlcaBeamMonitor::histosMap_
HistosContainer histosMap_
Definition: AlcaBeamMonitor.h:67
AlcaBeamMonitor::monitorName_
std::string monitorName_
Definition: AlcaBeamMonitor.h:49
BeamFitter::resetPVFitter
void resetPVFitter()
Definition: BeamFitter.h:69
AlcaBeamMonitor::parameters_
edm::ParameterSet parameters_
Definition: AlcaBeamMonitor.h:48
BeamFitter::resetRefTime
void resetRefTime()
Definition: BeamFitter.h:58
AlcaBeamMonitor::~AlcaBeamMonitor ( )
override

Definition at line 99 of file AlcaBeamMonitor.cc.

References theBeamFitter_, and thePVFitter_.

99  {
100  if (theBeamFitter_ != nullptr) {
101  delete theBeamFitter_;
102  }
103 
104  if (thePVFitter_ != nullptr) {
105  delete thePVFitter_;
106  }
107 }
AlcaBeamMonitor::theBeamFitter_
BeamFitter * theBeamFitter_
Definition: AlcaBeamMonitor.h:57
AlcaBeamMonitor::thePVFitter_
PVFitter * thePVFitter_
Definition: AlcaBeamMonitor.h:58

Member Function Documentation

void AlcaBeamMonitor::analyze ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
overrideprotected

Definition at line 252 of file AlcaBeamMonitor.cc.

References beamSpotsMap_, cppFunctionSkipper::exception, reco::BeamSpot::Fake, dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), hD0Phi0_, hDxyBS_, primaryVertexLabel_, edm::Handle< T >::product(), L1TEGammaOffline_cfi::PVCollection, BeamFitter::readEvent(), PVFitter::readEvent(), scalerLabel_, theBeamFitter_, thePVFitter_, HLT_2018_cff::track, reco::BeamSpot::Tracker, trackLabel_, PDWG_EXOHSCP_cff::tracks, vertices_, and cms::Exception::what().

252  {
253  //------ BeamFitter
254  theBeamFitter_->readEvent(iEvent);
255  //------ PVFitter
256  thePVFitter_->readEvent(iEvent);
257 
258  if (beamSpotsMap_.find("DB") != beamSpotsMap_.end()) {
259  //------ Tracks
260  Handle<reco::TrackCollection> TrackCollection;
261  iEvent.getByToken(trackLabel_, TrackCollection);
262  const reco::TrackCollection* tracks = TrackCollection.product();
263  for (reco::TrackCollection::const_iterator track = tracks->begin(); track != tracks->end(); ++track) {
264  hD0Phi0_->Fill(track->phi(), -1 * track->dxy());
265  hDxyBS_->Fill(-1 * track->dxy(beamSpotsMap_["DB"].position()));
266  }
267  }
268 
269  //------ Primary Vertices
270  Handle<VertexCollection> PVCollection;
271  if (iEvent.getByToken(primaryVertexLabel_, PVCollection)) {
272  vertices_.push_back(*PVCollection.product());
273  }
274 
275  if (beamSpotsMap_.find("SC") == beamSpotsMap_.end()) {
276  //BeamSpot from file for this stream is = to the scalar BeamSpot
277  Handle<BeamSpot> recoBeamSpotHandle;
278  try {
279  iEvent.getByToken(scalerLabel_, recoBeamSpotHandle);
280  } catch (cms::Exception& exception) {
281  LogInfo("AlcaBeamMonitor") << exception.what();
282  return;
283  }
284  beamSpotsMap_["SC"] = *recoBeamSpotHandle;
285  if (beamSpotsMap_["SC"].BeamWidthX() != 0) {
286  beamSpotsMap_["SC"].setType(reco::BeamSpot::Tracker);
287  } else {
288  beamSpotsMap_["SC"].setType(reco::BeamSpot::Fake);
289  }
290  }
291 }
AlcaBeamMonitor::hD0Phi0_
MonitorElement * hD0Phi0_
Definition: AlcaBeamMonitor.h:61
AlcaBeamMonitor::trackLabel_
edm::EDGetTokenT< reco::TrackCollection > trackLabel_
Definition: AlcaBeamMonitor.h:51
AlcaBeamMonitor::hDxyBS_
MonitorElement * hDxyBS_
Definition: AlcaBeamMonitor.h:62
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:525
L1TEGammaOffline_cfi.PVCollection
PVCollection
Definition: L1TEGammaOffline_cfi.py:34
reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:14
AlcaBeamMonitor::primaryVertexLabel_
edm::EDGetTokenT< reco::VertexCollection > primaryVertexLabel_
Definition: AlcaBeamMonitor.h:50
cms::Exception::what
char const * what() const override
Definition: Exception.cc:103
BeamFitter::readEvent
void readEvent(const edm::Event &iEvent)
Definition: BeamFitter.cc:224
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:242
PDWG_EXOHSCP_cff.tracks
tracks
Definition: PDWG_EXOHSCP_cff.py:28
PVFitter::readEvent
void readEvent(const edm::Event &iEvent)
Definition: PVFitter.cc:93
AlcaBeamMonitor::scalerLabel_
edm::EDGetTokenT< reco::BeamSpot > scalerLabel_
Definition: AlcaBeamMonitor.h:52
AlcaBeamMonitor::beamSpotsMap_
BeamSpotContainer beamSpotsMap_
Definition: AlcaBeamMonitor.h:66
AlcaBeamMonitor::theBeamFitter_
BeamFitter * theBeamFitter_
Definition: AlcaBeamMonitor.h:57
AlcaBeamMonitor::thePVFitter_
PVFitter * thePVFitter_
Definition: AlcaBeamMonitor.h:58
edm::Handle::product
T const * product() const
Definition: Handle.h:69
HLT_2018_cff.track
track
Definition: HLT_2018_cff.py:10339
AlcaBeamMonitor::vertices_
std::vector< reco::VertexCollection > vertices_
Definition: AlcaBeamMonitor.h:71
void AlcaBeamMonitor::bookHistograms ( DQMStore::IBooker &  ibooker,
edm::Run const &  iRun,
edm::EventSetup const &  iSetup 
)
overrideprotectedvirtual

Implements DQMOneEDAnalyzer< edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, Args... >.

Definition at line 110 of file AlcaBeamMonitor.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::bookProfile(), hD0Phi0_, hDxyBS_, histosMap_, monitorName_, Skims_PA_cff::name, numberOfValuesToSave_, dqm::impl::MonitorElement::setAxisTitle(), dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), dqm::impl::MonitorElement::setLumiFlag(), AlCaHLTBitMon_QueryRunRegistry::string, theValuesContainer_, and overlapproblemtsosanalyzer_cfi::title.

110  {
111  string name;
112  string title;
113  ibooker.setCurrentFolder(monitorName_ + "Debug");
114  for (HistosContainer::iterator itM = histosMap_.begin(); itM != histosMap_.end(); itM++) {
115  for (map<string, MonitorElement*>::iterator itMM = itM->second["run"].begin(); itMM != itM->second["run"].end();
116  itMM++) {
117  name = string("h") + itM->first + itMM->first;
118  title = itM->first + "_{0} " + itMM->first;
119  if (itM->first == "x" || itM->first == "y") {
120  if (itMM->first == "Coordinate") {
121  itMM->second = ibooker.book1D(name, title, 1001, -0.2525, 0.2525);
122  } else if (itMM->first == "PrimaryVertex fit-DataBase" || itMM->first == "PrimaryVertex fit-BeamFit" ||
123  itMM->first == "PrimaryVertex fit-Scalers" || itMM->first == "PrimaryVertex-DataBase" ||
124  itMM->first == "PrimaryVertex-BeamFit" || itMM->first == "PrimaryVertex-Scalers") {
125  itMM->second = ibooker.book1D(name, title, 1001, -0.02525, 0.02525);
126  } else {
127  //assert(0);
128  }
129  } else if (itM->first == "z") {
130  if (itMM->first == "Coordinate") {
131  itMM->second = ibooker.book1D(name, title, 101, -5.05, 5.05);
132  } else if (itMM->first == "PrimaryVertex fit-DataBase" || itMM->first == "PrimaryVertex fit-BeamFit" ||
133  itMM->first == "PrimaryVertex fit-Scalers") {
134  itMM->second = ibooker.book1D(name, title, 101, -0.505, 0.505);
135  } else if (itMM->first == "PrimaryVertex-DataBase" || itMM->first == "PrimaryVertex-BeamFit" ||
136  itMM->first == "PrimaryVertex-Scalers") {
137  itMM->second = ibooker.book1D(name, title, 1001, -5.005, 5.005);
138  } else {
139  //assert(0);
140  }
141  } else if (itM->first == "sigmaX" || itM->first == "sigmaY") {
142  if (itMM->first == "Coordinate") {
143  itMM->second = ibooker.book1D(name, title, 100, 0, 0.015);
144  } else if (itMM->first == "PrimaryVertex fit-DataBase" || itMM->first == "PrimaryVertex fit-BeamFit" ||
145  itMM->first == "PrimaryVertex fit-Scalers" || itMM->first == "PrimaryVertex-DataBase" ||
146  itMM->first == "PrimaryVertex-BeamFit" || itMM->first == "PrimaryVertex-Scalers") {
147  itMM->second = nullptr;
148  } else {
149  //assert(0);
150  }
151  } else if (itM->first == "sigmaZ") {
152  if (itMM->first == "Coordinate") {
153  itMM->second = ibooker.book1D(name, title, 110, 0, 11);
154  } else if (itMM->first == "PrimaryVertex fit-DataBase" || itMM->first == "PrimaryVertex fit-BeamFit" ||
155  itMM->first == "PrimaryVertex fit-Scalers" || itMM->first == "PrimaryVertex-DataBase" ||
156  itMM->first == "PrimaryVertex-BeamFit" || itMM->first == "PrimaryVertex-Scalers") {
157  itMM->second = ibooker.book1D(name, title, 101, -5.05, 5.05);
158  } else {
159  //assert(0);
160  }
161  } else {
162  //assert(0);
163  }
164  if (itMM->second != nullptr) {
165  if (itMM->first == "Coordinate") {
166  itMM->second->setAxisTitle(itM->first + "_{0} (cm)", 1);
167  } else if (itMM->first == "PrimaryVertex fit-DataBase" || itMM->first == "PrimaryVertex fit-BeamFit" ||
168  itMM->first == "PrimaryVertex fit-Scalers" || itMM->first == "PrimaryVertex-DataBase" ||
169  itMM->first == "PrimaryVertex-BeamFit" || itMM->first == "PrimaryVertex-Scalers") {
170  itMM->second->setAxisTitle(itMM->first + " " + itM->first + "_{0} (cm)", 1);
171  }
172  itMM->second->setAxisTitle("Entries", 2);
173  }
174  }
175  }
176  ibooker.setCurrentFolder(monitorName_ + "Service");
177  theValuesContainer_ = ibooker.bookProfile("hHistoLumiValues",
178  "Histo Lumi Values",
179  3 * numberOfValuesToSave_,
180  0.,
181  3 * numberOfValuesToSave_,
182  100.,
183  -100.,
184  9000.,
185  " ");
186  theValuesContainer_->setLumiFlag();
187 
188  // create and cd into new folder
189  ibooker.setCurrentFolder(monitorName_ + "Validation");
190  //Book histograms
191  hD0Phi0_ = ibooker.bookProfile("hD0Phi0", "d_{0} vs. #phi_{0} (All Tracks)", 63, -3.15, 3.15, 100, -0.5, 0.5, "");
192  hD0Phi0_->setAxisTitle("#phi_{0} (rad)", 1);
193  hD0Phi0_->setAxisTitle("d_{0} (cm)", 2);
194 
195  ibooker.setCurrentFolder(monitorName_ + "Debug");
196  hDxyBS_ = ibooker.book1D("hDxyBS", "dxy_{0} w.r.t. Beam spot (All Tracks)", 100, -0.1, 0.1);
197  hDxyBS_->setAxisTitle("dxy_{0} w.r.t. Beam spot (cm)", 1);
198 }
AlcaBeamMonitor::hD0Phi0_
MonitorElement * hD0Phi0_
Definition: AlcaBeamMonitor.h:61
AlcaBeamMonitor::hDxyBS_
MonitorElement * hDxyBS_
Definition: AlcaBeamMonitor.h:62
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
dqm::impl::MonitorElement::setLumiFlag
void setLumiFlag()
this ME is meant to be stored for each luminosity section
Definition: MonitorElement.h:225
overlapproblemtsosanalyzer_cfi.title
title
Definition: overlapproblemtsosanalyzer_cfi.py:7
AlcaBeamMonitor::histosMap_
HistosContainer histosMap_
Definition: AlcaBeamMonitor.h:67
AlcaBeamMonitor::monitorName_
std::string monitorName_
Definition: AlcaBeamMonitor.h:49
AlcaBeamMonitor::theValuesContainer_
MonitorElement * theValuesContainer_
Definition: AlcaBeamMonitor.h:63
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
dqm::impl::MonitorElement::setAxisTitle
virtual void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:813
void AlcaBeamMonitor::dqmBeginLuminosityBlock ( const edm::LuminosityBlock &  iLumi,
const edm::EventSetup &  iSetup 
)
overrideprotected

Definition at line 201 of file AlcaBeamMonitor.cc.

References align::BeamSpot, beamSpotsMap_, cppFunctionSkipper::exception, reco::BeamSpot::Fake, edm::EventSetup::get(), BeamSpotObjects::GetBeamType(), BeamSpotObjects::GetBeamWidthX(), BeamSpotObjects::GetBeamWidthY(), BeamSpotObjects::GetBetaStar(), BeamSpotObjects::GetCovariance(), BeamSpotObjects::Getdxdz(), BeamSpotObjects::Getdydz(), BeamSpotObjects::GetEmittanceX(), BeamSpotObjects::GetEmittanceY(), BeamSpotObjects::GetSigmaZ(), BeamSpotObjects::GetX(), BeamSpotObjects::GetY(), BeamSpotObjects::GetZ(), mps_fire::i, edm::LuminosityBlockBase::id(), edm::ESHandleBase::isValid(), dqmiolumiharvest::j, edm::LuminosityBlockID::luminosityBlock(), makeMuonMisalignmentScenario::matrix, edm::ESHandle< T >::product(), dqm::impl::MonitorElement::Reset(), reco::BeamSpot::setBeamWidthY(), theValuesContainer_, reco::BeamSpot::Tracker, vertices_, and cms::Exception::what().

201  {
202  // Always create a beamspot group for each lumi weather we have results or not! Each Beamspot will be of unknown type!
203 
204  vertices_.clear();
205  theValuesContainer_->Reset();
206  beamSpotsMap_.clear();
207 
208  //Read BeamSpot from DB
209  ESHandle<BeamSpotObjects> bsDBHandle;
210  try {
211  iSetup.get<BeamSpotObjectsRcd>().get(bsDBHandle);
212  } catch (cms::Exception& exception) {
213  LogInfo("AlcaBeamMonitor") << exception.what();
214  return;
215  }
216  if (bsDBHandle.isValid()) { // check the product
217  const BeamSpotObjects* spotDB = bsDBHandle.product();
218 
219  // translate from BeamSpotObjects to reco::BeamSpot
220  BeamSpot::Point apoint(spotDB->GetX(), spotDB->GetY(), spotDB->GetZ());
221 
222  BeamSpot::CovarianceMatrix matrix;
223  for (int i = 0; i < 7; ++i) {
224  for (int j = 0; j < 7; ++j) {
225  matrix(i, j) = spotDB->GetCovariance(i, j);
226  }
227  }
228 
229  beamSpotsMap_["DB"] =
230  BeamSpot(apoint, spotDB->GetSigmaZ(), spotDB->Getdxdz(), spotDB->Getdydz(), spotDB->GetBeamWidthX(), matrix);
231 
232  BeamSpot* aSpot = &(beamSpotsMap_["DB"]);
233 
234  aSpot->setBeamWidthY(spotDB->GetBeamWidthY());
235  aSpot->setEmittanceX(spotDB->GetEmittanceX());
236  aSpot->setEmittanceY(spotDB->GetEmittanceY());
237  aSpot->setbetaStar(spotDB->GetBetaStar());
238 
239  if (spotDB->GetBeamType() == 2) {
240  aSpot->setType(reco::BeamSpot::Tracker);
241  } else {
242  aSpot->setType(reco::BeamSpot::Fake);
243  }
244  //LogInfo("AlcaBeamMonitor")
245  // << *aSpot << std::endl;
246  } else {
247  LogInfo("AlcaBeamMonitor") << "Database BeamSpot is not valid at lumi: " << iLumi.id().luminosityBlock();
248  }
249 }
BeamSpotObjects::Getdydz
double Getdydz() const
get dydz slope, crossing angle in YZ
Definition: BeamSpotObjects.h:81
edm::LuminosityBlockBase::id
LuminosityBlockID id() const
Definition: LuminosityBlockBase.h:44
reco::BeamSpot::CovarianceMatrix
math::Error< dimension >::type CovarianceMatrix
Definition: BeamSpot.h:29
mps_fire.i
i
Definition: mps_fire.py:341
BeamSpotObjects::GetY
double GetY() const
get Y beam position
Definition: BeamSpotObjects.h:69
BeamSpotObjects::GetSigmaZ
double GetSigmaZ() const
get sigma Z, RMS bunch length
Definition: BeamSpotObjects.h:73
BeamSpotObjects::GetBeamWidthX
double GetBeamWidthX() const
get average transverse beam width
Definition: BeamSpotObjects.h:75
reco::BeamSpot::Point
math::XYZPoint Point
point in the space
Definition: BeamSpot.h:27
cms::Exception::what
char const * what() const override
Definition: Exception.cc:103
BeamSpotObjects::GetBeamWidthY
double GetBeamWidthY() const
get average transverse beam width
Definition: BeamSpotObjects.h:77
BeamSpotObjects::GetEmittanceX
double GetEmittanceX() const
get emittance
Definition: BeamSpotObjects.h:103
dqm::impl::MonitorElement::Reset
virtual void Reset()
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:417
reco::BeamSpot::setBeamWidthY
void setBeamWidthY(double v)
Definition: BeamSpot.h:105
BeamSpotObjects::GetBeamType
int GetBeamType() const
get beam type
Definition: BeamSpotObjects.h:101
AlcaBeamMonitor::beamSpotsMap_
BeamSpotContainer beamSpotsMap_
Definition: AlcaBeamMonitor.h:66
BeamSpotObjects::GetZ
double GetZ() const
get Z beam position
Definition: BeamSpotObjects.h:71
BeamSpotObjects::Getdxdz
double Getdxdz() const
get dxdz slope, crossing angle in XZ
Definition: BeamSpotObjects.h:79
BeamSpotObjects::GetX
double GetX() const
get X beam position
Definition: BeamSpotObjects.h:67
edm::LuminosityBlockID::luminosityBlock
LuminosityBlockNumber_t luminosityBlock() const
Definition: LuminosityBlockID.h:42
BeamSpotObjects::GetBetaStar
double GetBetaStar() const
get beta star
Definition: BeamSpotObjects.h:107
BeamSpotObjects::GetCovariance
double GetCovariance(int i, int j) const
get i,j element of the full covariance matrix 7x7
Definition: BeamSpotObjects.h:83
edm::EventSetup::get
T get() const
Definition: EventSetup.h:73
BeamSpotObjects::GetEmittanceY
double GetEmittanceY() const
get emittance
Definition: BeamSpotObjects.h:105
edm::ESHandleBase::isValid
bool isValid() const
Definition: ESHandle.h:44
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
AlcaBeamMonitor::theValuesContainer_
MonitorElement * theValuesContainer_
Definition: AlcaBeamMonitor.h:63
AlcaBeamMonitor::vertices_
std::vector< reco::VertexCollection > vertices_
Definition: AlcaBeamMonitor.h:71
void AlcaBeamMonitor::dqmEndLuminosityBlock ( const edm::LuminosityBlock &  iLumi,
const edm::EventSetup &  iSetup 
)
overrideprotected

Definition at line 294 of file AlcaBeamMonitor.cc.

References beamSpotsMap_, DEFINE_FWK_MODULE, dqm::impl::MonitorElement::Fill(), dqmdumpme::first, BeamFitter::getBeamSpot(), PVFitter::getBeamSpot(), dqm::impl::MonitorElement::getTProfile(), timingPdfMaker::histo, histoByCategoryNames_, histosMap_, position, positionsMap_, funct::pow(), MetAnalyzer::pv(), PVFitter::resetAll(), BeamFitter::resetLSRange(), BeamFitter::resetPVFitter(), BeamFitter::resetRefTime(), BeamFitter::resetTrkVector(), PVFitter::runFitter(), BeamFitter::runPVandTrkFitter(), edm::second(), mathSSE::sqrt(), theBeamFitter_, thePVFitter_, theValuesContainer_, align::Tracker, varNamesV_, and vertices_.

294  {
295  if (theBeamFitter_->runPVandTrkFitter()) {
296  beamSpotsMap_["BF"] = theBeamFitter_->getBeamSpot();
297  }
298  theBeamFitter_->resetTrkVector();
299  theBeamFitter_->resetLSRange();
300  theBeamFitter_->resetRefTime();
301  theBeamFitter_->resetPVFitter();
302 
303  if (thePVFitter_->runFitter()) {
304  beamSpotsMap_["PV"] = thePVFitter_->getBeamSpot();
305  }
306  thePVFitter_->resetAll();
307 
308  // "PV,BF..." Value,Error
309  map<std::string, pair<double, double> > resultsMap;
310  vector<pair<double, double> > vertexResults;
311  MonitorElement* histo = nullptr;
312  int position = 0;
313  for (vector<string>::iterator itV = varNamesV_.begin(); itV != varNamesV_.end(); itV++) {
314  resultsMap.clear();
315  for (BeamSpotContainer::iterator itBS = beamSpotsMap_.begin(); itBS != beamSpotsMap_.end(); itBS++) {
316  if (itBS->second.type() == BeamSpot::Tracker) {
317  if (*itV == "x") {
318  resultsMap[itBS->first] = pair<double, double>(itBS->second.x0(), itBS->second.x0Error());
319  } else if (*itV == "y") {
320  resultsMap[itBS->first] = pair<double, double>(itBS->second.y0(), itBS->second.y0Error());
321  } else if (*itV == "z") {
322  resultsMap[itBS->first] = pair<double, double>(itBS->second.z0(), itBS->second.z0Error());
323  } else if (*itV == "sigmaX") {
324  resultsMap[itBS->first] = pair<double, double>(itBS->second.BeamWidthX(), itBS->second.BeamWidthXError());
325  } else if (*itV == "sigmaY") {
326  resultsMap[itBS->first] = pair<double, double>(itBS->second.BeamWidthY(), itBS->second.BeamWidthYError());
327  } else if (*itV == "sigmaZ") {
328  resultsMap[itBS->first] = pair<double, double>(itBS->second.sigmaZ(), itBS->second.sigmaZ0Error());
329  } else {
330  LogInfo("AlcaBeamMonitor") << "The histosMap_ has been built with the name " << *itV
331  << " that I can't recognize!";
332  //assert(0);
333  }
334  }
335  }
336  vertexResults.clear();
337  for (vector<VertexCollection>::iterator itPV = vertices_.begin(); itPV != vertices_.end(); itPV++) {
338  if (!itPV->empty()) {
339  for (VertexCollection::const_iterator pv = itPV->begin(); pv != itPV->end(); pv++) {
340  if (pv->isFake() || pv->tracksSize() < 10)
341  continue;
342  if (*itV == "x") {
343  vertexResults.push_back(pair<double, double>(pv->x(), pv->xError()));
344  } else if (*itV == "y") {
345  vertexResults.push_back(pair<double, double>(pv->y(), pv->yError()));
346  } else if (*itV == "z") {
347  vertexResults.push_back(pair<double, double>(pv->z(), pv->zError()));
348  } else if (*itV != "sigmaX" && *itV != "sigmaY" && *itV != "sigmaZ") {
349  LogInfo("AlcaBeamMonitor") << "The histosMap_ has been built with the name " << *itV
350  << " that I can't recognize!";
351  //assert(0);
352  }
353  }
354  }
355  }
356  /*
357  histoByCategoryNames_.insert( pair<string,string>("run", "Coordinate"));
358  histoByCategoryNames_.insert( pair<string,string>("run", "PrimaryVertex fit-DataBase"));
359  histoByCategoryNames_.insert( pair<string,string>("run", "PrimaryVertex fit-BeamFit"));
360  histoByCategoryNames_.insert( pair<string,string>("run", "PrimaryVertex fit-Scalers"));
361  histoByCategoryNames_.insert( pair<string,string>("run", "PrimaryVertex-DataBase"));
362  histoByCategoryNames_.insert( pair<string,string>("run", "PrimaryVertex-BeamFit"));
363  histoByCategoryNames_.insert( pair<string,string>("run", "PrimaryVertex-Scalers"));
364 
365  histoByCategoryNames_.insert( pair<string,string>("lumi", "Lumibased BeamSpotFit"));
366  histoByCategoryNames_.insert( pair<string,string>("lumi", "Lumibased PrimaryVertex"));
367  histoByCategoryNames_.insert( pair<string,string>("lumi", "Lumibased DataBase"));
368  histoByCategoryNames_.insert( pair<string,string>("lumi", "Lumibased Scalers"));
369  histoByCategoryNames_.insert( pair<string,string>("lumi", "Lumibased PrimaryVertex-DataBase fit"));
370  histoByCategoryNames_.insert( pair<string,string>("lumi", "Lumibased PrimaryVertex-Scalers fit"));
371  histoByCategoryNames_.insert( pair<string,string>("validation", "Lumibased Scalers-DataBase fit"));
372  histoByCategoryNames_.insert( pair<string,string>("validation", "Lumibased PrimaryVertex-DataBase"));
373  histoByCategoryNames_.insert( pair<string,string>("validation", "Lumibased PrimaryVertex-Scalers"));
374 */
375  for (multimap<string, string>::iterator itM = histoByCategoryNames_.begin(); itM != histoByCategoryNames_.end();
376  itM++) {
377  if (itM->first == "run" && (histo = histosMap_[*itV][itM->first][itM->second]) == nullptr) {
378  continue;
379  } else if (itM->first != "run") {
380  position = positionsMap_[*itV][itM->first][itM->second];
381  }
382  if (itM->second == "Coordinate") {
383  if (beamSpotsMap_.find("DB") != beamSpotsMap_.end()) {
384  histo->Fill(resultsMap["DB"].first);
385  }
386  } else if (itM->second == "PrimaryVertex fit-DataBase") {
387  if (resultsMap.find("PV") != resultsMap.end() && resultsMap.find("DB") != resultsMap.end()) {
388  histo->Fill(resultsMap["PV"].first - resultsMap["DB"].first);
389  }
390  } else if (itM->second == "PrimaryVertex fit-BeamFit") {
391  if (resultsMap.find("PV") != resultsMap.end() && resultsMap.find("BF") != resultsMap.end()) {
392  histo->Fill(resultsMap["PV"].first - resultsMap["BF"].first);
393  }
394  } else if (itM->second == "PrimaryVertex fit-Scalers") {
395  if (resultsMap.find("PV") != resultsMap.end() && resultsMap.find("SC") != resultsMap.end()) {
396  histo->Fill(resultsMap["PV"].first - resultsMap["SC"].first);
397  }
398  } else if (itM->second == "PrimaryVertex-DataBase") {
399  if (resultsMap.find("PV") != resultsMap.end() && resultsMap.find("DB") != resultsMap.end()) {
400  for (vector<pair<double, double> >::iterator itPV = vertexResults.begin(); itPV != vertexResults.end();
401  itPV++) {
402  histo->Fill(itPV->first - resultsMap["DB"].first);
403  }
404  }
405  } else if (itM->second == "PrimaryVertex-BeamFit") {
406  if (resultsMap.find("PV") != resultsMap.end() && resultsMap.find("BF") != resultsMap.end()) {
407  for (vector<pair<double, double> >::iterator itPV = vertexResults.begin(); itPV != vertexResults.end();
408  itPV++) {
409  histo->Fill(itPV->first - resultsMap["BF"].first);
410  }
411  }
412  } else if (itM->second == "PrimaryVertex-Scalers") {
413  if (resultsMap.find("PV") != resultsMap.end() && resultsMap.find("SC") != resultsMap.end()) {
414  for (vector<pair<double, double> >::iterator itPV = vertexResults.begin(); itPV != vertexResults.end();
415  itPV++) {
416  histo->Fill(itPV->first - resultsMap["SC"].first);
417  }
418  }
419  } else if (itM->second == "Lumibased BeamSpotFit") {
420  if (resultsMap.find("BF") != resultsMap.end()) {
421  theValuesContainer_->Fill(position, resultsMap["BF"].first); //Value
422  theValuesContainer_->Fill(position + 1, resultsMap["BF"].second); //Error
423  theValuesContainer_->Fill(position + 2, 1); //ok
424  }
425  } else if (itM->second == "Lumibased PrimaryVertex") {
426  if (resultsMap.find("PV") != resultsMap.end()) {
427  theValuesContainer_->Fill(position, resultsMap["PV"].first); //Value
428  theValuesContainer_->Fill(position + 1, resultsMap["PV"].second); //Error
429  theValuesContainer_->Fill(position + 2, 1); //ok
430  }
431  } else if (itM->second == "Lumibased DataBase") {
432  if (resultsMap.find("DB") != resultsMap.end()) {
433  theValuesContainer_->Fill(position, resultsMap["DB"].first); //Value
434  theValuesContainer_->Fill(position + 1, resultsMap["DB"].second); //Error
435  theValuesContainer_->Fill(position + 2, 1); //ok
436  }
437  } else if (itM->second == "Lumibased Scalers") {
438  if (resultsMap.find("SC") != resultsMap.end()) {
439  theValuesContainer_->Fill(position, resultsMap["SC"].first); //Value
440  theValuesContainer_->Fill(position + 1, resultsMap["SC"].second); //Error
441  theValuesContainer_->Fill(position + 2, 1); //ok
442  }
443  } else if (itM->second == "Lumibased PrimaryVertex-DataBase fit") {
444  if (resultsMap.find("PV") != resultsMap.end() && resultsMap.find("DB") != resultsMap.end()) {
445  theValuesContainer_->Fill(position, resultsMap["PV"].first - resultsMap["DB"].first); //Value
446  theValuesContainer_->Fill(
447  position + 1,
448  std::sqrt(std::pow(resultsMap["PV"].second, 2) + std::pow(resultsMap["DB"].second, 2))); //Error
449  theValuesContainer_->Fill(position + 2, 1); //ok
450  }
451  } else if (itM->second == "Lumibased PrimaryVertex-Scalers fit") {
452  if (resultsMap.find("PV") != resultsMap.end() && resultsMap.find("SC") != resultsMap.end()) {
453  theValuesContainer_->Fill(position, resultsMap["PV"].first - resultsMap["SC"].first); //Value
454  theValuesContainer_->Fill(
455  position + 1,
456  std::sqrt(std::pow(resultsMap["PV"].second, 2) + std::pow(resultsMap["SC"].second, 2))); //Error
457  theValuesContainer_->Fill(position + 2, 1); //ok
458  }
459  } else if (itM->second == "Lumibased Scalers-DataBase fit") {
460  if (resultsMap.find("SC") != resultsMap.end() && resultsMap.find("DB") != resultsMap.end()) {
461  theValuesContainer_->Fill(position, resultsMap["SC"].first - resultsMap["DB"].first); //Value
462  theValuesContainer_->Fill(
463  position + 1,
464  std::sqrt(std::pow(resultsMap["SC"].second, 2) + std::pow(resultsMap["DB"].second, 2))); //Error
465  theValuesContainer_->Fill(position + 2, 1); //ok
466  }
467  } else if (itM->second == "Lumibased PrimaryVertex-DataBase") {
468  if (resultsMap.find("DB") != resultsMap.end() && !vertexResults.empty()) {
469  for (vector<pair<double, double> >::iterator itPV = vertexResults.begin(); itPV != vertexResults.end();
470  itPV++) {
471  theValuesContainer_->Fill(position, (*itPV).first - resultsMap["DB"].first); //Value
472  }
473  /*
474  double error = 0;
475  if(vertexResults.size() != 0){
476  for(vector<pair<double,double> >::iterator itPV=vertexResults.begin(); itPV!=vertexResults.end(); itPV++){
477  error += std::pow((*itPV).first-resultsMap["DB"].first-theValuesContainer_->getTProfile()->GetBinContent(position+1),2.);
478  }
479  error = std::sqrt(error)/vertexResults.size();
480  }
481 // theValuesContainer_->Fill(position+1,std::sqrt(std::pow((*itPV).second,2)+std::pow(resultsMap["DB"].second,2)));//Error
482  theValuesContainer_->Fill(position+1,error);//Error
483 */
484  theValuesContainer_->Fill(position + 1,
485  theValuesContainer_->getTProfile()->GetBinError(position + 1)); //Error
486  theValuesContainer_->Fill(position + 2, 1); //ok
487  }
488  } else if (itM->second == "Lumibased PrimaryVertex-Scalers") {
489  if (resultsMap.find("SC") != resultsMap.end() && !vertexResults.empty()) {
490  for (vector<pair<double, double> >::iterator itPV = vertexResults.begin(); itPV != vertexResults.end();
491  itPV++) {
492  theValuesContainer_->Fill(position, (*itPV).first - resultsMap["SC"].first); //Value
493  }
494  /*
495  double error = 0;
496  if(vertexResults.size() != 0){
497  for(vector<pair<double,double> >::iterator itPV=vertexResults.begin(); itPV!=vertexResults.end(); itPV++){
498  error += std::pow((*itPV).first-resultsMap["SC"].first-theValuesContainer_->getTProfile()->GetBinContent(position+1),2.);
499  }
500  error = std::sqrt(error)/vertexResults.size();
501  }
502 // theValuesContainer_->Fill(position+1,std::sqrt(std::pow((*itPV).second,2)+std::pow(resultsMap["SC"].second,2)));//Error
503  theValuesContainer_->Fill(position+1,error);//Error
504 */
505  theValuesContainer_->Fill(position + 1,
506  theValuesContainer_->getTProfile()->GetBinError(position + 1)); //Error
507  theValuesContainer_->Fill(position + 2, 1); //ok
508  }
509  }
510  // else if(itM->second == "Lumibased Scalers-DataBase"){
511  // if(resultsMap.find("SC") != resultsMap.end() && resultsMap.find("DB") != resultsMap.end()){
512  // itHHH->second->Fill(bin,resultsMap["SC"].first-resultsMap["DB"].first);
513  // }
514  // }
515  else {
516  LogInfo("AlcaBeamMonitor") << "The histosMap_ have a histogram named " << itM->second
517  << " that I can't recognize in this loop!";
518  //assert(0);
519  }
520  }
521  }
522 }
BeamFitter::runPVandTrkFitter
bool runPVandTrkFitter()
Definition: BeamFitter.cc:412
AlcaBeamMonitor::histoByCategoryNames_
std::multimap< std::string, std::string > histoByCategoryNames_
Definition: AlcaBeamMonitor.h:70
BeamFitter::resetTrkVector
void resetTrkVector()
Definition: BeamFitter.h:55
AlcaBeamMonitor::positionsMap_
PositionContainer positionsMap_
Definition: AlcaBeamMonitor.h:68
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:242
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:215
AlcaBeamMonitor::varNamesV_
std::vector< std::string > varNamesV_
Definition: AlcaBeamMonitor.h:69
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
AlcaBeamMonitor::beamSpotsMap_
BeamSpotContainer beamSpotsMap_
Definition: AlcaBeamMonitor.h:66
AlcaBeamMonitor::theBeamFitter_
BeamFitter * theBeamFitter_
Definition: AlcaBeamMonitor.h:57
MetAnalyzer.pv
def pv(vc)
Definition: MetAnalyzer.py:7
BeamFitter::resetLSRange
void resetLSRange()
Definition: BeamFitter.h:57
AlcaBeamMonitor::thePVFitter_
PVFitter * thePVFitter_
Definition: AlcaBeamMonitor.h:58
dqm::impl::MonitorElement::getTProfile
virtual TProfile * getTProfile()
Definition: MonitorElement.cc:1199
PVFitter::runFitter
bool runFitter()
Definition: PVFitter.cc:304
AlcaBeamMonitor::histosMap_
HistosContainer histosMap_
Definition: AlcaBeamMonitor.h:67
PVFitter::resetAll
void resetAll()
Definition: PVFitter.h:80
BeamFitter::resetPVFitter
void resetPVFitter()
Definition: BeamFitter.h:69
position
static int position[264][3]
Definition: ReadPGInfo.cc:289
PVFitter::getBeamSpot
reco::BeamSpot getBeamSpot()
Definition: PVFitter.h:91
AlcaBeamMonitor::theValuesContainer_
MonitorElement * theValuesContainer_
Definition: AlcaBeamMonitor.h:63
BeamFitter::getBeamSpot
reco::BeamSpot getBeamSpot()
Definition: BeamFitter.h:90
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:30
AlcaBeamMonitor::vertices_
std::vector< reco::VertexCollection > vertices_
Definition: AlcaBeamMonitor.h:71
BeamFitter::resetRefTime
void resetRefTime()
Definition: BeamFitter.h:58

Member Data Documentation

edm::InputTag AlcaBeamMonitor::beamSpotLabel_
private

Definition at line 53 of file AlcaBeamMonitor.h.

BeamSpotContainer AlcaBeamMonitor::beamSpotsMap_
private

Definition at line 66 of file AlcaBeamMonitor.h.

Referenced by analyze(), dqmBeginLuminosityBlock(), and dqmEndLuminosityBlock().

MonitorElement* AlcaBeamMonitor::hD0Phi0_
private

Definition at line 61 of file AlcaBeamMonitor.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* AlcaBeamMonitor::hDxyBS_
private

Definition at line 62 of file AlcaBeamMonitor.h.

Referenced by analyze(), and bookHistograms().

std::multimap<std::string, std::string> AlcaBeamMonitor::histoByCategoryNames_
private

Definition at line 70 of file AlcaBeamMonitor.h.

Referenced by AlcaBeamMonitor(), and dqmEndLuminosityBlock().

HistosContainer AlcaBeamMonitor::histosMap_
private

Definition at line 67 of file AlcaBeamMonitor.h.

Referenced by AlcaBeamMonitor(), bookHistograms(), and dqmEndLuminosityBlock().

std::string AlcaBeamMonitor::monitorName_
private

Definition at line 49 of file AlcaBeamMonitor.h.

Referenced by AlcaBeamMonitor(), and bookHistograms().

int AlcaBeamMonitor::numberOfValuesToSave_
private

Definition at line 56 of file AlcaBeamMonitor.h.

Referenced by AlcaBeamMonitor(), and bookHistograms().

edm::ParameterSet AlcaBeamMonitor::parameters_
private

Definition at line 48 of file AlcaBeamMonitor.h.

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

PositionContainer AlcaBeamMonitor::positionsMap_
private

Definition at line 68 of file AlcaBeamMonitor.h.

Referenced by AlcaBeamMonitor(), and dqmEndLuminosityBlock().

edm::EDGetTokenT<reco::VertexCollection> AlcaBeamMonitor::primaryVertexLabel_
private

Definition at line 50 of file AlcaBeamMonitor.h.

Referenced by analyze().

edm::EDGetTokenT<reco::BeamSpot> AlcaBeamMonitor::scalerLabel_
private

Definition at line 52 of file AlcaBeamMonitor.h.

Referenced by analyze().

BeamFitter* AlcaBeamMonitor::theBeamFitter_
private

Definition at line 57 of file AlcaBeamMonitor.h.

Referenced by AlcaBeamMonitor(), analyze(), dqmEndLuminosityBlock(), and ~AlcaBeamMonitor().

PVFitter* AlcaBeamMonitor::thePVFitter_
private

Definition at line 58 of file AlcaBeamMonitor.h.

Referenced by AlcaBeamMonitor(), analyze(), dqmEndLuminosityBlock(), and ~AlcaBeamMonitor().

MonitorElement* AlcaBeamMonitor::theValuesContainer_
private

Definition at line 63 of file AlcaBeamMonitor.h.

Referenced by bookHistograms(), dqmBeginLuminosityBlock(), and dqmEndLuminosityBlock().

edm::EDGetTokenT<reco::TrackCollection> AlcaBeamMonitor::trackLabel_
private

Definition at line 51 of file AlcaBeamMonitor.h.

Referenced by analyze().

std::vector<std::string> AlcaBeamMonitor::varNamesV_
private

Definition at line 69 of file AlcaBeamMonitor.h.

Referenced by AlcaBeamMonitor(), and dqmEndLuminosityBlock().

std::vector<reco::VertexCollection> AlcaBeamMonitor::vertices_
private

Definition at line 71 of file AlcaBeamMonitor.h.

Referenced by analyze(), dqmBeginLuminosityBlock(), and dqmEndLuminosityBlock().

Generated for CMSSW Reference Manual by doxygen 1.8.11