CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Attributes
DiMuonHistograms Class Reference

#include <DiMuonHistograms.h>

Inheritance diagram for DiMuonHistograms:
DQMEDAnalyzer edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > > edm::stream::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &)
 
void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 
 DiMuonHistograms (const edm::ParameterSet &pset)
 
virtual ~DiMuonHistograms ()
 
- Public Member Functions inherited from DQMEDAnalyzer
virtual void beginRun (edm::Run const &, edm::EventSetup const &) final
 
virtual void beginStream (edm::StreamID id) final
 
virtual void dqmBeginRun (edm::Run const &, edm::EventSetup const &)
 
 DQMEDAnalyzer (void)
 
virtual void endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
virtual void endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
uint32_t streamId () const
 
- Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
 EDAnalyzer ()=default
 
- Public Member Functions inherited from edm::stream::EDAnalyzerBase
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDAnalyzerBase ()
 
- Public Member Functions inherited from edm::EDConsumerBase
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 		itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Attributes

int etaBBin
 
int etaBin
 
double etaBMax
 
double etaBMin
 
double EtaCutMax
 
double EtaCutMin
 
int etaEBin
 
double etaECMax
 
double etaECMin
 
std::string EtaName
 
int etaOvlpBin
 
std::vector< MonitorElement * > GlbGlbMuon_HM
 
std::vector< MonitorElement * > GlbGlbMuon_LM
 
double HighMassMax
 
double HighMassMin
 
double LowMassMax
 
double LowMassMin
 
std::string metname
 
edm::ParameterSet parameters
 
std::vector< MonitorElement * > SoftSoftMuon
 
std::vector< MonitorElement * > StaTrkMuon_HM
 
std::vector< MonitorElement * > StaTrkMuon_LM
 
edm::EDGetTokenT< reco::BeamSpot > theBeamSpotLabel_
 
edm::EDGetTokenT
< reco::MuonCollection > 		theMuonCollectionLabel_
 
MuonServiceProxy * theService
 
edm::EDGetTokenT
< reco::VertexCollection > 		theVertexLabel_
 
std::vector< MonitorElement * > TightTightMuon
 
std::vector< MonitorElement * > TrkTrkMuon_HM
 
std::vector< MonitorElement * > TrkTrkMuon_LM
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...> CacheTypes
 
typedef CacheTypes::GlobalCache GlobalCache
 
typedef AbilityChecker< T...> HasAbility
 
typedef
CacheTypes::LuminosityBlockCache 		LuminosityBlockCache
 
typedef
LuminosityBlockContextT
< LuminosityBlockCache,
RunCache, GlobalCache > 		LuminosityBlockContext
 
typedef
CacheTypes::LuminosityBlockSummaryCache 		LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache,
GlobalCache > 		RunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 
- Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase ModuleType
 
- Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr
< dqmDetails::NoCache > 		globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
 
static std::shared_ptr
< dqmDetails::NoCache > 		globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
 
static void globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
 
static void globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
 
- Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::stream::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
- 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 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)
 

Detailed Description

Class DiMuonHistograms

DQM monitoring for dimuon mass

Author: S.Folgueras, U. Oviedo

Definition at line 33 of file DiMuonHistograms.h.

Constructor & Destructor Documentation

DiMuonHistograms::DiMuonHistograms ( const edm::ParameterSet &  pset)

Definition at line 32 of file DiMuonHistograms.cc.

References Parameters::parameters.

32  {
33 
34  // initialise parameters:
35  parameters = pSet;
36 
37 
38  // declare consumes:
39  theMuonCollectionLabel_ = consumes<reco::MuonCollection> (parameters.getParameter<edm::InputTag>("MuonCollection"));
40  theBeamSpotLabel_ = mayConsume<reco::BeamSpot> (parameters.getParameter<edm::InputTag>("BeamSpotLabel"));
41  theVertexLabel_ = consumes<reco::VertexCollection>(parameters.getParameter<edm::InputTag>("VertexLabel"));
42 
43  etaBin = parameters.getParameter<int>("etaBin");
44  etaBBin = parameters.getParameter<int>("etaBBin");
45  etaEBin = parameters.getParameter<int>("etaEBin");
46 
47  etaBMin = parameters.getParameter<double>("etaBMin");
48  etaBMax = parameters.getParameter<double>("etaBMax");
49  etaECMin = parameters.getParameter<double>("etaECMin");
50  etaECMax = parameters.getParameter<double>("etaECMax");
51 
52  LowMassMin = parameters.getParameter<double>("LowMassMin");
53  LowMassMax = parameters.getParameter<double>("LowMassMax");
54  HighMassMin = parameters.getParameter<double>("HighMassMin");
55  HighMassMax = parameters.getParameter<double>("HighMassMax");
56 
57 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
DiMuonHistograms::parameters
edm::ParameterSet parameters
Definition: DiMuonHistograms.h:47
DiMuonHistograms::theMuonCollectionLabel_
edm::EDGetTokenT< reco::MuonCollection > theMuonCollectionLabel_
Definition: DiMuonHistograms.h:85
DiMuonHistograms::theBeamSpotLabel_
edm::EDGetTokenT< reco::BeamSpot > theBeamSpotLabel_
Definition: DiMuonHistograms.h:87
DiMuonHistograms::theVertexLabel_
edm::EDGetTokenT< reco::VertexCollection > theVertexLabel_
Definition: DiMuonHistograms.h:86
DiMuonHistograms::~DiMuonHistograms ( )
virtual

Definition at line 59 of file DiMuonHistograms.cc.

59 { }

Member Function Documentation

void DiMuonHistograms::analyze ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
virtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 90 of file DiMuonHistograms.cc.

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), DeDxDiscriminatorTools::charge(), relativeConstraints::error, edm::Event::getByToken(), muon::isSoftMuon(), muon::isTightMuon(), edm::HandleBase::isValid(), LogTrace, metname, patZpeak::muons, reco::BeamSpot::position(), position, and reco::BeamSpot::sigmaZ().

90  {
91 
92  LogTrace(metname)<<"[DiMuonHistograms] Analyze the mu in different eta regions";
93  edm::Handle<reco::MuonCollection> muons;
94  iEvent.getByToken(theMuonCollectionLabel_, muons);
95 
96  // =================================================================================
97  // Look for the Primary Vertex (and use the BeamSpot instead, if you can't find it):
98  reco::Vertex::Point posVtx;
99  reco::Vertex::Error errVtx;
100  unsigned int theIndexOfThePrimaryVertex = 999.;
101 
102  edm::Handle<reco::VertexCollection> vertex;
103  iEvent.getByToken(theVertexLabel_, vertex);
104  if (vertex.isValid()){
105  for (unsigned int ind=0; ind<vertex->size(); ++ind) {
106  if ( (*vertex)[ind].isValid() && !((*vertex)[ind].isFake()) ) {
107  theIndexOfThePrimaryVertex = ind;
108  break;
109  }
110  }
111  }
112 
113  if (theIndexOfThePrimaryVertex<100) {
114  posVtx = ((*vertex)[theIndexOfThePrimaryVertex]).position();
115  errVtx = ((*vertex)[theIndexOfThePrimaryVertex]).error();
116  }
117  else {
118  LogInfo("RecoMuonValidator") << "reco::PrimaryVertex not found, use BeamSpot position instead\n";
119 
120  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
121  iEvent.getByToken(theBeamSpotLabel_,recoBeamSpotHandle);
122  reco::BeamSpot bs = *recoBeamSpotHandle;
123 
124  posVtx = bs.position();
125  errVtx(0,0) = bs.BeamWidthX();
126  errVtx(1,1) = bs.BeamWidthY();
127  errVtx(2,2) = bs.sigmaZ();
128  }
129 
130  const reco::Vertex vtx(posVtx,errVtx);
131 
132  if(!muons.isValid()) return;
133 
134  // Loop on muon collection
135  TLorentzVector Mu1, Mu2;
136  float charge = 99.;
137  float InvMass = -99.;
138 
139  for (reco::MuonCollection::const_iterator recoMu1 = muons->begin(); recoMu1!=muons->end(); ++recoMu1) {
140  LogTrace(metname)<<"[DiMuonHistograms] loop over 1st muon"<<endl;
141 
142  // Loop on second muons to fill invariant mass plots
143  for (reco::MuonCollection::const_iterator recoMu2 = recoMu1; recoMu2!=muons->end(); ++recoMu2){
144  LogTrace(metname)<<"[DiMuonHistograms] loop over 2nd muon"<<endl;
145  if (recoMu1==recoMu2) continue;
146 
147  // Global-Global Muon
148  if (recoMu1->isGlobalMuon() && recoMu2->isGlobalMuon()) {
149  LogTrace(metname)<<"[DiMuonHistograms] Glb-Glb pair"<<endl;
150  reco::TrackRef recoCombinedGlbTrack1 = recoMu1->combinedMuon();
151  reco::TrackRef recoCombinedGlbTrack2 = recoMu2->combinedMuon();
152  Mu1.SetPxPyPzE(recoCombinedGlbTrack1->px(), recoCombinedGlbTrack1->py(),recoCombinedGlbTrack1->pz(), recoCombinedGlbTrack1->p());
153  Mu2.SetPxPyPzE(recoCombinedGlbTrack2->px(), recoCombinedGlbTrack2->py(),recoCombinedGlbTrack2->pz(), recoCombinedGlbTrack2->p());
154 
155  charge = recoCombinedGlbTrack1->charge()*recoCombinedGlbTrack2->charge();
156  if (charge < 0) {
157  InvMass = (Mu1+Mu2).M();
158  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
159  if (iEtaRegion==0) {EtaCutMin= 0.; EtaCutMax=2.4; }
160  if (iEtaRegion==1) {EtaCutMin= etaBMin; EtaCutMax=etaBMax; }
161  if (iEtaRegion==2) {EtaCutMin= etaECMin; EtaCutMax=etaECMax; }
162 
163  if(fabs(recoCombinedGlbTrack1->eta())>EtaCutMin && fabs(recoCombinedGlbTrack1->eta())<EtaCutMax &&
164  fabs(recoCombinedGlbTrack2->eta())>EtaCutMin && fabs(recoCombinedGlbTrack2->eta())<EtaCutMax){
165  if (InvMass < LowMassMax) GlbGlbMuon_LM[iEtaRegion]->Fill(InvMass);
166  if (InvMass > HighMassMin) GlbGlbMuon_HM[iEtaRegion]->Fill(InvMass);
167  }
168  }
169  }
170  // Also Tight-Tight Muon Selection
171 
172  if ( muon::isTightMuon(*recoMu1, vtx) &&
173  muon::isTightMuon(*recoMu2, vtx) ) {
174 
175  LogTrace(metname)<<"[DiMuonHistograms] Tight-Tight pair"<<endl;
176  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
177  if (iEtaRegion==0) {EtaCutMin= 0.; EtaCutMax=2.4; }
178  if (iEtaRegion==1) {EtaCutMin= etaBMin; EtaCutMax=etaBMax; }
179  if (iEtaRegion==2) {EtaCutMin= etaECMin; EtaCutMax=etaECMax; }
180 
181  if(fabs(recoCombinedGlbTrack1->eta())>EtaCutMin && fabs(recoCombinedGlbTrack1->eta())<EtaCutMax &&
182  fabs(recoCombinedGlbTrack2->eta())>EtaCutMin && fabs(recoCombinedGlbTrack2->eta())<EtaCutMax){
183  if (InvMass > 55. && InvMass < 125.) TightTightMuon[iEtaRegion]->Fill(InvMass);
184  }
185  }
186  }
187  }
188 
189  // Now check for STA-TRK
190  if (recoMu2->isStandAloneMuon() && recoMu1->isTrackerMuon()) {
191  LogTrace(metname)<<"[DiMuonHistograms] STA-Trk pair"<<endl;
192  reco::TrackRef recoStaTrack = recoMu2->standAloneMuon();
193  reco::TrackRef recoTrack = recoMu1->track();
194  Mu2.SetPxPyPzE(recoStaTrack->px(), recoStaTrack->py(),recoStaTrack->pz(), recoStaTrack->p());
195  Mu1.SetPxPyPzE(recoTrack->px(), recoTrack->py(),recoTrack->pz(), recoTrack->p());
196 
197  charge = recoStaTrack->charge()*recoTrack->charge();
198  if (charge < 0) {
199  InvMass = (Mu1+Mu2).M();
200  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
201  if (iEtaRegion==0) {EtaCutMin= 0.; EtaCutMax=2.4; }
202  if (iEtaRegion==1) {EtaCutMin= etaBMin; EtaCutMax=etaBMax; }
203  if (iEtaRegion==2) {EtaCutMin= etaECMin; EtaCutMax=etaECMax; }
204 
205  if(fabs(recoStaTrack->eta())>EtaCutMin && fabs(recoStaTrack->eta())<EtaCutMax &&
206  fabs(recoTrack->eta())>EtaCutMin && fabs(recoTrack->eta())<EtaCutMax){
207  if (InvMass < LowMassMax) StaTrkMuon_LM[iEtaRegion]->Fill(InvMass);
208  if (InvMass > HighMassMin) StaTrkMuon_HM[iEtaRegion]->Fill(InvMass);
209  }
210  }
211  }
212  }
213  if (recoMu1->isStandAloneMuon() && recoMu2->isTrackerMuon()) {
214  LogTrace(metname)<<"[DiMuonHistograms] STA-Trk pair"<<endl;
215  reco::TrackRef recoStaTrack = recoMu1->standAloneMuon();
216  reco::TrackRef recoTrack = recoMu2->track();
217  Mu1.SetPxPyPzE(recoStaTrack->px(), recoStaTrack->py(),recoStaTrack->pz(), recoStaTrack->p());
218  Mu2.SetPxPyPzE(recoTrack->px(), recoTrack->py(),recoTrack->pz(), recoTrack->p());
219 
220  charge = recoStaTrack->charge()*recoTrack->charge();
221  if (charge < 0) {
222  InvMass = (Mu1+Mu2).M();
223  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
224  if (iEtaRegion==0) {EtaCutMin= 0.; EtaCutMax=2.4; }
225  if (iEtaRegion==1) {EtaCutMin= etaBMin; EtaCutMax=etaBMax; }
226  if (iEtaRegion==2) {EtaCutMin= etaECMin; EtaCutMax=etaECMax; }
227 
228  if(fabs(recoStaTrack->eta())>EtaCutMin && fabs(recoStaTrack->eta())<EtaCutMax &&
229  fabs(recoTrack->eta())>EtaCutMin && fabs(recoTrack->eta())<EtaCutMax){
230  if (InvMass < LowMassMax) StaTrkMuon_LM[iEtaRegion]->Fill(InvMass);
231  if (InvMass > HighMassMin) StaTrkMuon_HM[iEtaRegion]->Fill(InvMass);
232  }
233  }
234  }
235  }
236 
237  // TRK-TRK dimuon
238  if (recoMu1->isTrackerMuon() && recoMu2->isTrackerMuon()) {
239  LogTrace(metname)<<"[DiMuonHistograms] Trk-Trk dimuon pair"<<endl;
240  reco::TrackRef recoTrack2 = recoMu2->track();
241  reco::TrackRef recoTrack1 = recoMu1->track();
242  Mu2.SetPxPyPzE(recoTrack2->px(), recoTrack2->py(),recoTrack2->pz(), recoTrack2->p());
243  Mu1.SetPxPyPzE(recoTrack1->px(), recoTrack1->py(),recoTrack1->pz(), recoTrack1->p());
244 
245  charge = recoTrack1->charge()*recoTrack2->charge();
246  if (charge < 0) {
247  InvMass = (Mu1+Mu2).M();
248  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
249  if (iEtaRegion==0) {EtaCutMin= 0.; EtaCutMax=2.4; }
250  if (iEtaRegion==1) {EtaCutMin= etaBMin; EtaCutMax=etaBMax; }
251  if (iEtaRegion==2) {EtaCutMin= etaECMin; EtaCutMax=etaECMax; }
252 
253  if(fabs(recoTrack1->eta())>EtaCutMin && fabs(recoTrack1->eta())<EtaCutMax &&
254  fabs(recoTrack2->eta())>EtaCutMin && fabs(recoTrack2->eta())<EtaCutMax){
255  if (InvMass < LowMassMax) TrkTrkMuon_LM[iEtaRegion]->Fill(InvMass);
256  if (InvMass > HighMassMin) TrkTrkMuon_HM[iEtaRegion]->Fill(InvMass);
257  }
258  }
259  }
260 
261 
262  LogTrace(metname)<<"[DiMuonHistograms] Soft-Soft pair"<<endl;
263 
264  if (muon::isSoftMuon(*recoMu1, vtx) &&
265  muon::isSoftMuon(*recoMu2, vtx) ) {
266 
267  if (charge < 0) {
268  InvMass = (Mu1+Mu2).M();
269  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
270  if (iEtaRegion==0) {EtaCutMin= 0.; EtaCutMax=2.4; }
271  if (iEtaRegion==1) {EtaCutMin= etaBMin; EtaCutMax=etaBMax; }
272  if (iEtaRegion==2) {EtaCutMin= etaECMin; EtaCutMax=etaECMax; }
273 
274  if(fabs(recoTrack1->eta())>EtaCutMin && fabs(recoTrack1->eta())<EtaCutMax &&
275  fabs(recoTrack2->eta())>EtaCutMin && fabs(recoTrack2->eta())<EtaCutMax){
276  SoftSoftMuon[iEtaRegion]->Fill(InvMass);
277  }
278  }
279  }
280  }
281  }
282  } //muon2
283  } //Muon1
284 }
DiMuonHistograms::metname
std::string metname
Definition: DiMuonHistograms.h:50
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:434
DiMuonHistograms::GlbGlbMuon_LM
std::vector< MonitorElement * > GlbGlbMuon_LM
Definition: DiMuonHistograms.h:74
reco::Vertex::Error
math::Error< dimension >::type Error
covariance error matrix (3x3)
Definition: Vertex.h:43
DeDxDiscriminatorTools::charge
double charge(const std::vector< uint8_t > &Ampls)
Definition: DeDxDiscriminatorTools.cc:43
muon::isSoftMuon
bool isSoftMuon(const reco::Muon &, const reco::Vertex &)
Definition: MuonSelectors.cc:767
DiMuonHistograms::theMuonCollectionLabel_
edm::EDGetTokenT< reco::MuonCollection > theMuonCollectionLabel_
Definition: DiMuonHistograms.h:85
DiMuonHistograms::theBeamSpotLabel_
edm::EDGetTokenT< reco::BeamSpot > theBeamSpotLabel_
Definition: DiMuonHistograms.h:87
DiMuonHistograms::SoftSoftMuon
std::vector< MonitorElement * > SoftSoftMuon
Definition: DiMuonHistograms.h:82
reco::Vertex::Point
math::XYZPoint Point
point in the space
Definition: Vertex.h:39
reco::BeamSpot::BeamWidthX
double BeamWidthX() const
beam width X
Definition: BeamSpot.h:86
DiMuonHistograms::TrkTrkMuon_LM
std::vector< MonitorElement * > TrkTrkMuon_LM
Definition: DiMuonHistograms.h:78
DiMuonHistograms::StaTrkMuon_LM
std::vector< MonitorElement * > StaTrkMuon_LM
Definition: DiMuonHistograms.h:76
DiMuonHistograms::theVertexLabel_
edm::EDGetTokenT< reco::VertexCollection > theVertexLabel_
Definition: DiMuonHistograms.h:86
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
DiMuonHistograms::TrkTrkMuon_HM
std::vector< MonitorElement * > TrkTrkMuon_HM
Definition: DiMuonHistograms.h:79
reco::BeamSpot::sigmaZ
double sigmaZ() const
sigma z
Definition: BeamSpot.h:80
reco::BeamSpot::BeamWidthY
double BeamWidthY() const
beam width Y
Definition: BeamSpot.h:88
patZpeak.muons
tuple muons
Definition: patZpeak.py:38
position
static int position[264][3]
Definition: ReadPGInfo.cc:509
DiMuonHistograms::StaTrkMuon_HM
std::vector< MonitorElement * > StaTrkMuon_HM
Definition: DiMuonHistograms.h:77
DiMuonHistograms::GlbGlbMuon_HM
std::vector< MonitorElement * > GlbGlbMuon_HM
Definition: DiMuonHistograms.h:75
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:62
DiMuonHistograms::TightTightMuon
std::vector< MonitorElement * > TightTightMuon
Definition: DiMuonHistograms.h:81
muon::isTightMuon
bool isTightMuon(const reco::Muon &, const reco::Vertex &)
Definition: MuonSelectors.cc:745
void DiMuonHistograms::bookHistograms ( DQMStore::IBooker &  ibooker,
edm::Run const &  ,
edm::EventSetup const &   
)
overridevirtual

Implements DQMEDAnalyzer.

Definition at line 61 of file DiMuonHistograms.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::cd(), and DQMStore::IBooker::setCurrentFolder().

63  {
64 
65  ibooker.cd();
66  ibooker.setCurrentFolder("Muons/DiMuonHistograms");
67 
68  int nBin = 0;
69  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
70  if (iEtaRegion==0) { EtaName = ""; nBin = etaBin;}
71  if (iEtaRegion==1) { EtaName = "_Barrel"; nBin = etaBBin;}
72  if (iEtaRegion==2) { EtaName = "_EndCap"; nBin = etaEBin;}
73 
74  GlbGlbMuon_LM.push_back(ibooker.book1D("GlbGlbMuon_LM"+EtaName,"InvMass_{GLB,GLB}"+EtaName,nBin, LowMassMin, LowMassMax));
75  TrkTrkMuon_LM.push_back(ibooker.book1D("TrkTrkMuon_LM"+EtaName,"InvMass_{TRK,TRK}"+EtaName,nBin, LowMassMin, LowMassMax));
76  StaTrkMuon_LM.push_back(ibooker.book1D("StaTrkMuon_LM"+EtaName,"InvMass_{STA,TRK}"+EtaName,nBin, LowMassMin, LowMassMax));
77 
78  GlbGlbMuon_HM.push_back(ibooker.book1D("GlbGlbMuon_HM"+EtaName,"InvMass_{GLB,GLB}"+EtaName,nBin, HighMassMin, HighMassMax));
79  TrkTrkMuon_HM.push_back(ibooker.book1D("TrkTrkMuon_HM"+EtaName,"InvMass_{TRK,TRK}"+EtaName,nBin, HighMassMin, HighMassMax));
80  StaTrkMuon_HM.push_back(ibooker.book1D("StaTrkMuon_HM"+EtaName,"InvMass_{STA,TRK}"+EtaName,nBin, HighMassMin, HighMassMax));
81 
82  // arround the Z peak
83  TightTightMuon.push_back(ibooker.book1D("TightTightMuon"+EtaName,"InvMass_{Tight,Tight}"+EtaName,nBin, 55.0, 125.0));
84 
85  // low-mass resonances
86  SoftSoftMuon.push_back(ibooker.book1D("SoftSoftMuon"+EtaName,"InvMass_{Soft,Soft}"+EtaName,nBin, 5.0, 55.0));
87  }
88 }
DQMStore::IBooker::cd
void cd(void)
Definition: DQMStore.cc:266
DiMuonHistograms::GlbGlbMuon_LM
std::vector< MonitorElement * > GlbGlbMuon_LM
Definition: DiMuonHistograms.h:74
DiMuonHistograms::SoftSoftMuon
std::vector< MonitorElement * > SoftSoftMuon
Definition: DiMuonHistograms.h:82
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:113
DiMuonHistograms::TrkTrkMuon_LM
std::vector< MonitorElement * > TrkTrkMuon_LM
Definition: DiMuonHistograms.h:78
DiMuonHistograms::StaTrkMuon_LM
std::vector< MonitorElement * > StaTrkMuon_LM
Definition: DiMuonHistograms.h:76
DiMuonHistograms::TrkTrkMuon_HM
std::vector< MonitorElement * > TrkTrkMuon_HM
Definition: DiMuonHistograms.h:79
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:274
DiMuonHistograms::StaTrkMuon_HM
std::vector< MonitorElement * > StaTrkMuon_HM
Definition: DiMuonHistograms.h:77
DiMuonHistograms::GlbGlbMuon_HM
std::vector< MonitorElement * > GlbGlbMuon_HM
Definition: DiMuonHistograms.h:75
DiMuonHistograms::TightTightMuon
std::vector< MonitorElement * > TightTightMuon
Definition: DiMuonHistograms.h:81
DiMuonHistograms::EtaName
std::string EtaName
Definition: DiMuonHistograms.h:59

Member Data Documentation

int DiMuonHistograms::etaBBin
private

Definition at line 54 of file DiMuonHistograms.h.

int DiMuonHistograms::etaBin
private

Definition at line 53 of file DiMuonHistograms.h.

double DiMuonHistograms::etaBMax
private

Definition at line 64 of file DiMuonHistograms.h.

double DiMuonHistograms::etaBMin
private

Definition at line 63 of file DiMuonHistograms.h.

double DiMuonHistograms::EtaCutMax
private

Definition at line 62 of file DiMuonHistograms.h.

double DiMuonHistograms::EtaCutMin
private

Definition at line 61 of file DiMuonHistograms.h.

int DiMuonHistograms::etaEBin
private

Definition at line 55 of file DiMuonHistograms.h.

double DiMuonHistograms::etaECMax
private

Definition at line 66 of file DiMuonHistograms.h.

double DiMuonHistograms::etaECMin
private

Definition at line 65 of file DiMuonHistograms.h.

std::string DiMuonHistograms::EtaName
private

Definition at line 59 of file DiMuonHistograms.h.

int DiMuonHistograms::etaOvlpBin
private

Definition at line 56 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::GlbGlbMuon_HM
private

Definition at line 75 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::GlbGlbMuon_LM
private

Definition at line 74 of file DiMuonHistograms.h.

double DiMuonHistograms::HighMassMax
private

Definition at line 72 of file DiMuonHistograms.h.

double DiMuonHistograms::HighMassMin
private

Definition at line 71 of file DiMuonHistograms.h.

double DiMuonHistograms::LowMassMax
private

Definition at line 70 of file DiMuonHistograms.h.

double DiMuonHistograms::LowMassMin
private

Definition at line 69 of file DiMuonHistograms.h.

std::string DiMuonHistograms::metname
private

Definition at line 50 of file DiMuonHistograms.h.

edm::ParameterSet DiMuonHistograms::parameters
private

Definition at line 47 of file DiMuonHistograms.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

std::vector<MonitorElement*> DiMuonHistograms::SoftSoftMuon
private

Definition at line 82 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::StaTrkMuon_HM
private

Definition at line 77 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::StaTrkMuon_LM
private

Definition at line 76 of file DiMuonHistograms.h.

edm::EDGetTokenT<reco::BeamSpot> DiMuonHistograms::theBeamSpotLabel_
private

Definition at line 87 of file DiMuonHistograms.h.

edm::EDGetTokenT<reco::MuonCollection> DiMuonHistograms::theMuonCollectionLabel_
private

Definition at line 85 of file DiMuonHistograms.h.

MuonServiceProxy* DiMuonHistograms::theService
private

Definition at line 46 of file DiMuonHistograms.h.

edm::EDGetTokenT<reco::VertexCollection> DiMuonHistograms::theVertexLabel_
private

Definition at line 86 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::TightTightMuon
private

Definition at line 81 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::TrkTrkMuon_HM
private

Definition at line 79 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::TrkTrkMuon_LM
private

Definition at line 78 of file DiMuonHistograms.h.

Generated for CMSSW Reference Manual by doxygen 1.8.5