#include <DiMuonHistograms.h>

Inheritance diagram for DiMuonHistograms:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

	DiMuonHistograms (const edm::ParameterSet &pset)

	~DiMuonHistograms () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

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

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

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

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

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

	~DQMEDAnalyzer () override=default

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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

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)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
int	etaBBin

int	etaBin

double	etaBMax

double	etaBMin

double	EtaCutMax

double	EtaCutMin

int	etaEBin

double	etaECMax

double	etaECMin

std::string	EtaName [3]

int	etaOvlpBin

std::vector< MonitorElement * >	GlbGlbMuon_HM

std::vector< MonitorElement * >	GlbGlbMuon_LM

double	HighMassMax

double	HighMassMin

std::vector< MonitorElement * >	LooseLooseMuon

std::vector< MonitorElement * >	LooseLooseMuonBadFrac

double	LowMassMax

double	LowMassMin

std::vector< MonitorElement * >	MediumMediumMuon

std::vector< MonitorElement * >	MediumMediumMuonBadFrac

std::string	metname

int	nGlbGlb

int	nLooseLoose

int	nMediumMedium

int	nTightTight

edm::ParameterSet	parameters

std::vector< MonitorElement * >	SoftSoftMuon

std::vector< MonitorElement * >	SoftSoftMuonBadFrac

std::vector< MonitorElement * >	StaTrkMuon_HM

std::vector< MonitorElement * >	StaTrkMuon_LM

MonitorElement *	test

edm::EDGetTokenT< reco::BeamSpot >	theBeamSpotLabel_

std::string	theFolder

edm::EDGetTokenT< edm::View< reco::Muon > >	theMuonCollectionLabel_

MuonServiceProxy *	theService

edm::EDGetTokenT< reco::VertexCollection >	theVertexLabel_

std::vector< MonitorElement * >	TightTightMuon

std::vector< MonitorElement * >	TightTightMuonBadFrac

std::vector< MonitorElement * >	TrkTrkMuon_HM

std::vector< MonitorElement * >	TrkTrkMuon_LM

Additional Inherited Members
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 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)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

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 conversionPostprocessing_cfi::etaBin.

                                                              {
   
   // initialise parameters:
   parameters = pSet;
 
   // counter
   nTightTight = 0;
   nMediumMedium = 0;
   nLooseLoose = 0;
   nGlbGlb     = 0;
 
   // declare consumes:
   theMuonCollectionLabel_ = consumes<edm::View<reco::Muon> >  (parameters.getParameter<edm::InputTag>("MuonCollection"));
   theVertexLabel_          = consumes<reco::VertexCollection>(parameters.getParameter<edm::InputTag>("VertexLabel"));
 
   theBeamSpotLabel_        = mayConsume<reco::BeamSpot>      (parameters.getParameter<edm::InputTag>("BeamSpotLabel"));
 
   etaBin = parameters.getParameter<int>("etaBin");
   etaBBin = parameters.getParameter<int>("etaBBin");
   etaEBin = parameters.getParameter<int>("etaEBin");
  
   etaBMin = parameters.getParameter<double>("etaBMin");
   etaBMax = parameters.getParameter<double>("etaBMax");
   etaECMin = parameters.getParameter<double>("etaECMin");
   etaECMax = parameters.getParameter<double>("etaECMax");
   
   LowMassMin = parameters.getParameter<double>("LowMassMin");  
   LowMassMax = parameters.getParameter<double>("LowMassMax");
   HighMassMin = parameters.getParameter<double>("HighMassMin");
   HighMassMax = parameters.getParameter<double>("HighMassMax");
 
   theFolder = parameters.getParameter<string>("folder");
 }

DiMuonHistograms::~DiMuonHistograms ( )

override

Definition at line 66 of file DiMuonHistograms.cc.

66 { }

Member Function Documentation

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

override

Definition at line 104 of file DiMuonHistograms.cc.

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), ALCARECOTkAlJpsiMuMu_cff::charge, relativeConstraints::error, edm::Event::getByToken(), muon::isLooseMuon(), muon::isMediumMuon(), muon::isSoftMuon(), muon::isTightMuon(), edm::HandleBase::isValid(), LogTrace, metname, extraflags_cff::muons, reco::BeamSpot::position(), position, reco::BeamSpot::sigmaZ(), and badGlobalMuonTaggersAOD_cff::vtx.

                                                                                  {
 
   LogTrace(metname)<<"[DiMuonHistograms] Analyze the mu in different eta regions";
   edm::Handle<edm::View<reco::Muon> > muons; 
   iEvent.getByToken(theMuonCollectionLabel_, muons);
 
   // =================================================================================
   // Look for the Primary Vertex (and use the BeamSpot instead, if you can't find it):
   reco::Vertex::Point posVtx;
   reco::Vertex::Error errVtx;
   unsigned int theIndexOfThePrimaryVertex = 999.;
   
   edm::Handle<reco::VertexCollection> vertex;
   iEvent.getByToken(theVertexLabel_, vertex);
   if (vertex.isValid()){
     for (unsigned int ind=0; ind<vertex->size(); ++ind) {
       if ( (*vertex)[ind].isValid() && !((*vertex)[ind].isFake()) ) {
     theIndexOfThePrimaryVertex = ind;
     break;
       }
     }
   }
 
   if (theIndexOfThePrimaryVertex<100) {
     posVtx = ((*vertex)[theIndexOfThePrimaryVertex]).position();
     errVtx = ((*vertex)[theIndexOfThePrimaryVertex]).error();
   }   
   else {
     LogInfo("RecoMuonValidator") << "reco::PrimaryVertex not found, use BeamSpot position instead\n";
     
     edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
     iEvent.getByToken(theBeamSpotLabel_,recoBeamSpotHandle);
     reco::BeamSpot bs = *recoBeamSpotHandle;
     
     posVtx = bs.position();
     errVtx(0,0) = bs.BeamWidthX();
     errVtx(1,1) = bs.BeamWidthY();
     errVtx(2,2) = bs.sigmaZ();
   }
   
   const reco::Vertex vtx(posVtx,errVtx);
   
   if(!muons.isValid()) return;
 
   // Loop on muon collection
   TLorentzVector Mu1, Mu2;
   float charge = 99.;
   float InvMass = -99.;
 
   //Eta regions
   double EtaCutMin[] = {0, etaBMin, etaECMin};
   double EtaCutMax[] = {2.4, etaBMax, etaECMax};
 
   for (edm::View<reco::Muon>::const_iterator muon1 = muons->begin(); muon1 != muons->end(); ++muon1){
     LogTrace(metname)<<"[DiMuonHistograms] loop over 1st muon"<<endl;
 
     // Loop on second muons to fill invariant mass plots
     for (edm::View<reco::Muon>::const_iterator muon2 = muon1; muon2 != muons->end(); ++muon2){
       LogTrace(metname)<<"[DiMuonHistograms] loop over 2nd muon"<<endl;
       if (muon1==muon2) continue;
       
       // Global-Global Muon
       if (muon1->isGlobalMuon() && muon2->isGlobalMuon()) {
     LogTrace(metname)<<"[DiMuonHistograms] Glb-Glb pair"<<endl;
     reco::TrackRef recoCombinedGlbTrack1 = muon1->combinedMuon();
     reco::TrackRef recoCombinedGlbTrack2 = muon2->combinedMuon();
     Mu1.SetPxPyPzE(recoCombinedGlbTrack1->px(), recoCombinedGlbTrack1->py(),recoCombinedGlbTrack1->pz(), recoCombinedGlbTrack1->p());
     Mu2.SetPxPyPzE(recoCombinedGlbTrack2->px(), recoCombinedGlbTrack2->py(),recoCombinedGlbTrack2->pz(), recoCombinedGlbTrack2->p());
     
     charge  = recoCombinedGlbTrack1->charge()*recoCombinedGlbTrack2->charge();
     if (charge < 0) {
       InvMass = (Mu1+Mu2).M();
       for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
         if(fabs(recoCombinedGlbTrack1->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack1->eta())<EtaCutMax[iEtaRegion] && 
            fabs(recoCombinedGlbTrack2->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack2->eta())<EtaCutMax[iEtaRegion]){
           if (InvMass < LowMassMax)  GlbGlbMuon_LM[iEtaRegion]->Fill(InvMass);
           if (InvMass > HighMassMin) GlbGlbMuon_HM[iEtaRegion]->Fill(InvMass);
         }
       }
     }
         // Also Tight-Tight Muon Selection
     if ( muon::isTightMuon(*muon1, vtx)  && 
          muon::isTightMuon(*muon2, vtx) ) { 
       test->Fill(InvMass);
       LogTrace(metname)<<"[DiMuonHistograms] Tight-Tight pair"<<endl;
       if (charge < 0){
         for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
           if(fabs(recoCombinedGlbTrack1->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack1->eta())<EtaCutMax[iEtaRegion] && 
          fabs(recoCombinedGlbTrack2->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack2->eta())<EtaCutMax[iEtaRegion]){
         if (InvMass > 55. && InvMass < 125.) {
             TightTightMuon[iEtaRegion]->Fill(InvMass);
             TightTightMuonBadFrac[iEtaRegion]->Fill(muon1->innerTrack()->lost()/muon1->innerTrack()->found());
         }
           }
         }
       }
     }
     // Also Medium-Medium Muon Selection
     if ( muon::isMediumMuon(*muon1)  && 
          muon::isMediumMuon(*muon2) ) { 
       test->Fill(InvMass);
       LogTrace(metname)<<"[DiMuonHistograms] Medium-Medium pair"<<endl;
       if (charge < 0){
         for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
           if(fabs(recoCombinedGlbTrack1->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack1->eta())<EtaCutMax[iEtaRegion] && 
          fabs(recoCombinedGlbTrack2->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack2->eta())<EtaCutMax[iEtaRegion]){
         if (InvMass > 55. && InvMass < 125.) {
             MediumMediumMuon[iEtaRegion]->Fill(InvMass);
             MediumMediumMuonBadFrac[iEtaRegion]->Fill(muon1->innerTrack()->lost()/muon1->innerTrack()->found());
         }
           }
         }
       }
     }
     // Also Loose-Loose Muon Selection
     if ( muon::isLooseMuon(*muon1)  && 
          muon::isLooseMuon(*muon2) ) { 
       test->Fill(InvMass);
       LogTrace(metname)<<"[DiMuonHistograms] Loose-Loose pair"<<endl;
       if (charge < 0){
         for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
           if(fabs(recoCombinedGlbTrack1->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack1->eta())<EtaCutMax[iEtaRegion] && 
          fabs(recoCombinedGlbTrack2->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack2->eta())<EtaCutMax[iEtaRegion]){
         if (InvMass > 55. && InvMass < 125.) {
             LooseLooseMuon[iEtaRegion]->Fill(InvMass);
             LooseLooseMuonBadFrac[iEtaRegion]->Fill(muon1->innerTrack()->lost()/muon1->innerTrack()->found());
         }
           }
         }
       }
     }   
       }
     
       // Now check for STA-TRK 
       if (muon2->isStandAloneMuon() && muon1->isTrackerMuon()) {
     LogTrace(metname)<<"[DiMuonHistograms] STA-Trk pair"<<endl;
     reco::TrackRef recoStaTrack = muon2->standAloneMuon();
     reco::TrackRef recoTrack    = muon1->track();
     Mu2.SetPxPyPzE(recoStaTrack->px(), recoStaTrack->py(),recoStaTrack->pz(), recoStaTrack->p());
     Mu1.SetPxPyPzE(recoTrack->px(), recoTrack->py(),recoTrack->pz(), recoTrack->p());
 
         charge  = recoStaTrack->charge()*recoTrack->charge();
         if (charge < 0) {
           InvMass = (Mu1+Mu2).M();
       for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
         if(fabs(recoStaTrack->eta())>EtaCutMin[iEtaRegion] && fabs(recoStaTrack->eta())<EtaCutMax[iEtaRegion] &&
            fabs(recoTrack->eta())>EtaCutMin[iEtaRegion] && fabs(recoTrack->eta())<EtaCutMax[iEtaRegion]){
           if (InvMass < LowMassMax)  StaTrkMuon_LM[iEtaRegion]->Fill(InvMass);
           if (InvMass > HighMassMin) StaTrkMuon_HM[iEtaRegion]->Fill(InvMass);
         }
       }
     }
       }
       if (muon1->isStandAloneMuon() && muon2->isTrackerMuon()) {
     LogTrace(metname)<<"[DiMuonHistograms] STA-Trk pair"<<endl;
     reco::TrackRef recoStaTrack = muon1->standAloneMuon();
     reco::TrackRef recoTrack    = muon2->track();
     Mu1.SetPxPyPzE(recoStaTrack->px(), recoStaTrack->py(),recoStaTrack->pz(), recoStaTrack->p());
     Mu2.SetPxPyPzE(recoTrack->px(), recoTrack->py(),recoTrack->pz(), recoTrack->p());
 
         charge  = recoStaTrack->charge()*recoTrack->charge();
         if (charge < 0) {
           InvMass = (Mu1+Mu2).M();
       for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
         if(fabs(recoStaTrack->eta())>EtaCutMin[iEtaRegion] && fabs(recoStaTrack->eta())<EtaCutMax[iEtaRegion] &&
            fabs(recoTrack->eta())>EtaCutMin[iEtaRegion] && fabs(recoTrack->eta())<EtaCutMax[iEtaRegion]){
           if (InvMass < LowMassMax)  StaTrkMuon_LM[iEtaRegion]->Fill(InvMass);
           if (InvMass > HighMassMin) StaTrkMuon_HM[iEtaRegion]->Fill(InvMass);
         }
       }
     }
       }
 
       // TRK-TRK dimuon 
       if (muon1->isTrackerMuon() && muon2->isTrackerMuon()) {
     LogTrace(metname)<<"[DiMuonHistograms] Trk-Trk dimuon pair"<<endl;
     reco::TrackRef recoTrack2 = muon2->track();
     reco::TrackRef recoTrack1 = muon1->track();
     Mu2.SetPxPyPzE(recoTrack2->px(), recoTrack2->py(),recoTrack2->pz(), recoTrack2->p());
     Mu1.SetPxPyPzE(recoTrack1->px(), recoTrack1->py(),recoTrack1->pz(), recoTrack1->p());
     
     charge  = recoTrack1->charge()*recoTrack2->charge();
         if (charge < 0) {
           InvMass = (Mu1+Mu2).M();
       for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
         if(fabs(recoTrack1->eta())>EtaCutMin[iEtaRegion] && fabs(recoTrack1->eta())<EtaCutMax[iEtaRegion] &&
            fabs(recoTrack2->eta())>EtaCutMin[iEtaRegion] && fabs(recoTrack2->eta())<EtaCutMax[iEtaRegion]){
           if (InvMass < LowMassMax)  TrkTrkMuon_LM[iEtaRegion]->Fill(InvMass);
           if (InvMass > HighMassMin) TrkTrkMuon_HM[iEtaRegion]->Fill(InvMass);
         }
       }
     }
 
 
     LogTrace(metname)<<"[DiMuonHistograms] Soft-Soft pair"<<endl;
 
     if (muon::isSoftMuon(*muon1, vtx)  && 
         muon::isSoftMuon(*muon2, vtx) ) { 
       
       if (charge < 0) {
         InvMass = (Mu1+Mu2).M();
         for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
           
           if(fabs(recoTrack1->eta())>EtaCutMin[iEtaRegion] && fabs(recoTrack1->eta())<EtaCutMax[iEtaRegion] &&
          fabs(recoTrack2->eta())>EtaCutMin[iEtaRegion] && fabs(recoTrack2->eta())<EtaCutMax[iEtaRegion]){
             SoftSoftMuon[iEtaRegion]->Fill(InvMass); 
             SoftSoftMuonBadFrac[iEtaRegion]->Fill(muon1->innerTrack()->lost()/muon1->innerTrack()->found());
           }
         }
       }
     } 
       }
     } //muon2
   } //Muon1
    
 
 }

void DiMuonHistograms::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 68 of file DiMuonHistograms.cc.

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

                                               {
   
   ibooker.cd();
   ibooker.setCurrentFolder(theFolder);  
 
   int nBin[3] = {etaBin,etaBBin,etaEBin};
   EtaName[0]  = ""; EtaName[1] = "_Barrel"; EtaName[2] = "_EndCap";
   test = ibooker.book1D("test","InvMass_{Tight,Tight}",100, 0., 200.);
   for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
     
     GlbGlbMuon_LM.push_back(ibooker.book1D("GlbGlbMuon_LM"+EtaName[iEtaRegion],"InvMass_{GLB,GLB}"+EtaName[iEtaRegion],nBin[iEtaRegion], LowMassMin, LowMassMax));
     TrkTrkMuon_LM.push_back(ibooker.book1D("TrkTrkMuon_LM"+EtaName[iEtaRegion],"InvMass_{TRK,TRK}"+EtaName[iEtaRegion],nBin[iEtaRegion], LowMassMin, LowMassMax));
     StaTrkMuon_LM.push_back(ibooker.book1D("StaTrkMuon_LM"+EtaName[iEtaRegion],"InvMass_{STA,TRK}"+EtaName[iEtaRegion],nBin[iEtaRegion], LowMassMin, LowMassMax));
     
     GlbGlbMuon_HM.push_back(ibooker.book1D("GlbGlbMuon_HM"+EtaName[iEtaRegion],"InvMass_{GLB,GLB}"+EtaName[iEtaRegion],nBin[iEtaRegion], HighMassMin, HighMassMax));
     TrkTrkMuon_HM.push_back(ibooker.book1D("TrkTrkMuon_HM"+EtaName[iEtaRegion],"InvMass_{TRK,TRK}"+EtaName[iEtaRegion],nBin[iEtaRegion], HighMassMin, HighMassMax));
     StaTrkMuon_HM.push_back(ibooker.book1D("StaTrkMuon_HM"+EtaName[iEtaRegion],"InvMass_{STA,TRK}"+EtaName[iEtaRegion],nBin[iEtaRegion], HighMassMin, HighMassMax));
     
     // arround the Z peak
     TightTightMuon.push_back(ibooker.book1D("TightTightMuon"+EtaName[iEtaRegion],"InvMass_{Tight,Tight}"+EtaName[iEtaRegion], nBin[iEtaRegion], LowMassMin, LowMassMax));
     MediumMediumMuon.push_back(ibooker.book1D("MediumMediumMuon"+EtaName[iEtaRegion],"InvMass_{Medium,Medium}"+EtaName[iEtaRegion],nBin[iEtaRegion], LowMassMin, LowMassMax));
     LooseLooseMuon.push_back(ibooker.book1D("LooseLooseMuon"+EtaName[iEtaRegion],"InvMass_{Loose,Loose}"+EtaName[iEtaRegion], nBin[iEtaRegion], LowMassMin, LowMassMax));
     //Fraction of bad hits in the tracker track to the total
     TightTightMuonBadFrac.push_back(ibooker.book1D("TightTightMuonBadFrac"+EtaName[iEtaRegion],"BadFrac_{Tight,Tight}"+EtaName[iEtaRegion], 10, 0, 0.4));
     MediumMediumMuonBadFrac.push_back(ibooker.book1D("MediumMediumMuonBadFrac"+EtaName[iEtaRegion],"BadFrac_{Medium,Medium}"+EtaName[iEtaRegion], 10, 0, 0.4));
     LooseLooseMuonBadFrac.push_back(ibooker.book1D("LooseLooseMuonBadFrac"+EtaName[iEtaRegion],"BadFrac_{Loose,Loose}"+EtaName[iEtaRegion], 10, 0, 0.4));
 
     // low-mass resonances
     SoftSoftMuon.push_back(ibooker.book1D("SoftSoftMuon"+EtaName[iEtaRegion],"InvMass_{Soft,Soft}"+EtaName[iEtaRegion],nBin[iEtaRegion], 0.0, 55.0));
     SoftSoftMuonBadFrac.push_back(ibooker.book1D("SoftSoftMuonBadFrac"+EtaName[iEtaRegion],"BadFrac_{Soft,Soft}"+EtaName[iEtaRegion], 10, 0, 0.4));
 
   }
 }

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[3]

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 76 of file DiMuonHistograms.h.

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

private

Definition at line 75 of file DiMuonHistograms.h.

double DiMuonHistograms::HighMassMax

private

Definition at line 73 of file DiMuonHistograms.h.

double DiMuonHistograms::HighMassMin

private

Definition at line 72 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::LooseLooseMuon

private

Definition at line 82 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::LooseLooseMuonBadFrac

private

Definition at line 86 of file DiMuonHistograms.h.

double DiMuonHistograms::LowMassMax

private

Definition at line 71 of file DiMuonHistograms.h.

double DiMuonHistograms::LowMassMin

private

Definition at line 70 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::MediumMediumMuon

private

Definition at line 83 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::MediumMediumMuonBadFrac

private

Definition at line 87 of file DiMuonHistograms.h.

std::string DiMuonHistograms::metname

private

Definition at line 50 of file DiMuonHistograms.h.

int DiMuonHistograms::nGlbGlb

private

Definition at line 103 of file DiMuonHistograms.h.

int DiMuonHistograms::nLooseLoose

private

Definition at line 102 of file DiMuonHistograms.h.

int DiMuonHistograms::nMediumMedium

private

Definition at line 101 of file DiMuonHistograms.h.

int DiMuonHistograms::nTightTight

private

Definition at line 100 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 85 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::SoftSoftMuonBadFrac

private

Definition at line 89 of file DiMuonHistograms.h.

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

private

Definition at line 78 of file DiMuonHistograms.h.

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

private

Definition at line 77 of file DiMuonHistograms.h.

MonitorElement* DiMuonHistograms::test

private

Definition at line 91 of file DiMuonHistograms.h.

Referenced by edmIntegrityCheck.IntegrityCheck::report(), and edmIntegrityCheck.IntegrityCheck::structured().

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

private

Definition at line 96 of file DiMuonHistograms.h.

std::string DiMuonHistograms::theFolder

private

Definition at line 98 of file DiMuonHistograms.h.

edm::EDGetTokenT<edm::View<reco::Muon> > DiMuonHistograms::theMuonCollectionLabel_

private

Definition at line 94 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 95 of file DiMuonHistograms.h.

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

private

Definition at line 84 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::TightTightMuonBadFrac

private

Definition at line 88 of file DiMuonHistograms.h.

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

private

Definition at line 80 of file DiMuonHistograms.h.

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

private

Definition at line 79 of file DiMuonHistograms.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation