MuonKinVsEtaAnalyzer Class Reference

#include <MuonKinVsEtaAnalyzer.h>

Inheritance diagram for MuonKinVsEtaAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &)
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	MuonKinVsEtaAnalyzer (const edm::ParameterSet &pSet)
	Constructor. More...
	~MuonKinVsEtaAnalyzer ()
	Destructor. More...
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
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Attributes
std::vector< MonitorElement * >	chi2GlbTrack
std::vector< MonitorElement * >	chi2HighPtTrack
std::vector< MonitorElement * >	chi2LooseTrack
std::vector< MonitorElement * >	chi2MediumTrack
std::vector< MonitorElement * >	chi2probGlbTrack
std::vector< MonitorElement * >	chi2probHighPtTrack
std::vector< MonitorElement * >	chi2probLooseTrack
std::vector< MonitorElement * >	chi2probMediumTrack
std::vector< MonitorElement * >	chi2probSoftTrack
std::vector< MonitorElement * >	chi2probStaTrack
std::vector< MonitorElement * >	chi2probTightTrack
std::vector< MonitorElement * >	chi2probTrack
std::vector< MonitorElement * >	chi2SoftTrack
std::vector< MonitorElement * >	chi2StaTrack
std::vector< MonitorElement * >	chi2TightTrack
std::vector< MonitorElement * >	chi2Track
int	chiBin
double	chiMax
double	chiMin
double	chiprobMax
double	chiprobMin
int	etaBin
double	etaBMax
double	etaBMin
double	EtaCutMax
double	EtaCutMin
double	etaECMax
double	etaECMin
std::vector< MonitorElement * >	etaGlbTrack
std::vector< MonitorElement * >	etaHighPtTrack
std::vector< MonitorElement * >	etaLooseTrack
double	etaMax
std::vector< MonitorElement * >	etaMediumTrack
double	etaMin
double	etaOvlpMax
double	etaOvlpMin
std::vector< MonitorElement * >	etaSoftTrack
std::vector< MonitorElement * >	etaStaTrack
std::vector< MonitorElement * >	etaTightTrack
std::vector< MonitorElement * >	etaTrack
std::string	metname
edm::ParameterSet	parameters
int	pBin
std::vector< MonitorElement * >	pGlbTrack
int	phiBin
std::vector< MonitorElement * >	pHighPtTrack
std::vector< MonitorElement * >	phiGlbTrack
std::vector< MonitorElement * >	phiHighPtTrack
std::vector< MonitorElement * >	phiLooseTrack
double	phiMax
std::vector< MonitorElement * >	phiMediumTrack
double	phiMin
std::vector< MonitorElement * >	phiSoftTrack
std::vector< MonitorElement * >	phiStaTrack
std::vector< MonitorElement * >	phiTightTrack
std::vector< MonitorElement * >	phiTrack
std::vector< MonitorElement * >	pLooseTrack
double	pMax
std::vector< MonitorElement * >	pMediumTrack
double	pMin
std::vector< MonitorElement * >	pSoftTrack
std::vector< MonitorElement * >	pStaTrack
int	ptBin
std::vector< MonitorElement * >	ptGlbTrack
std::vector< MonitorElement * >	ptHighPtTrack
std::vector< MonitorElement * >	pTightTrack
std::vector< MonitorElement * >	ptLooseTrack
double	ptMax
std::vector< MonitorElement * >	ptMediumTrack
double	ptMin
std::vector< MonitorElement * >	pTrack
std::vector< MonitorElement * >	ptSoftTrack
std::vector< MonitorElement * >	ptStaTrack
std::vector< MonitorElement * >	ptTightTrack
std::vector< MonitorElement * >	ptTrack
edm::EDGetTokenT< reco::BeamSpot >	theBeamSpotLabel_
DQMStore *	theDbe
edm::EDGetTokenT < reco::MuonCollection >	theMuonCollectionLabel_
MuonServiceProxy *	theService
edm::EDGetTokenT < reco::VertexCollection >	theVertexLabel_

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
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
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::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

Definition at line 34 of file MuonKinVsEtaAnalyzer.h.

Constructor & Destructor Documentation

MuonKinVsEtaAnalyzer::MuonKinVsEtaAnalyzer

(

const edm::ParameterSet &

pSet

)

Constructor.

Definition at line 20 of file MuonKinVsEtaAnalyzer.cc.

References diMuonHistograms_cfi::etaBin, HLT_25ns10e33_v2_cff::etaMax, HLT_25ns10e33_v2_cff::etaMin, LogTrace, metname, MuonServiceProxy_cff::MuonServiceProxy, cppFunctionSkipper::operator, HLT_25ns10e33_v2_cff::parameters, HLT_25ns10e33_v2_cff::pMin, and PtMinSelector_cfg::ptMin.

                                                                      {
  LogTrace(metname)<<"[MuonKinVsEtaAnalyzer] Parameters initialization";
  
  parameters = pSet;

  // the services
  theService = new MuonServiceProxy(parameters.getParameter<ParameterSet>("ServiceParameters"));
  theDbe = edm::Service<DQMStore>().operator->();

  theMuonCollectionLabel_  = consumes<reco::MuonCollection>(parameters.getParameter<InputTag>("MuonCollection"));
  theVertexLabel_          = consumes<reco::VertexCollection>(parameters.getParameter<edm::InputTag>("VertexLabel"));
  theBeamSpotLabel_        = mayConsume<reco::BeamSpot>      (parameters.getParameter<edm::InputTag>("BeamSpotLabel"));

  // Parameters
  etaBin = parameters.getParameter<int>("etaBin");
  etaMin = parameters.getParameter<double>("etaMin");
  etaMax = parameters.getParameter<double>("etaMax");

  phiBin = parameters.getParameter<int>("phiBin");
  phiMin = parameters.getParameter<double>("phiMin");
  phiMax = parameters.getParameter<double>("phiMax");

  pBin = parameters.getParameter<int>("pBin");
  pMin = parameters.getParameter<double>("pMin");
  pMax = parameters.getParameter<double>("pMax");

  ptBin = parameters.getParameter<int>("ptBin");
  ptMin = parameters.getParameter<double>("ptMin");
  ptMax = parameters.getParameter<double>("ptMax");

  chiBin = parameters.getParameter<int>("chiBin");
  chiMin = parameters.getParameter<double>("chiMin");
  chiMax = parameters.getParameter<double>("chiMax");
  chiprobMin = parameters.getParameter<double>("chiprobMin");
  chiprobMax = parameters.getParameter<double>("chiprobMax");

  etaBMin = parameters.getParameter<double>("etaBMin");
  etaBMax = parameters.getParameter<double>("etaBMax");
  etaECMin = parameters.getParameter<double>("etaECMin");
  etaECMax = parameters.getParameter<double>("etaECMax");
  etaOvlpMin = parameters.getParameter<double>("etaOvlpMin");
  etaOvlpMax = parameters.getParameter<double>("etaOvlpMax");

}

MuonKinVsEtaAnalyzer::~MuonKinVsEtaAnalyzer

(

)

Destructor.

Definition at line 64 of file MuonKinVsEtaAnalyzer.cc.

                                            { 
  delete theService;
}

Member Function Documentation

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

virtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 141 of file MuonKinVsEtaAnalyzer.cc.

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), reco::Muon::combinedMuon(), gather_cfg::cout, relativeConstraints::error, edm::Event::getByToken(), reco::Muon::isGlobalMuon(), muon::isHighPtMuon(), muon::isLooseMuon(), muon::isMediumMuon(), muon::isSoftMuon(), reco::Muon::isStandAloneMuon(), muon::isTightMuon(), reco::Muon::isTrackerMuon(), edm::HandleBase::isValid(), LogTrace, metname, patZpeak::muons, reco::BeamSpot::position(), position, reco::BeamSpot::sigmaZ(), reco::Muon::standAloneMuon(), and reco::Muon::track().

                                                                                     {
  
  LogTrace(metname)<<"[MuonKinVsEtaAnalyzer] Analyze the mu in different eta regions";
  theService->update(iSetup);
  
  edm::Handle<reco::MuonCollection> 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);
  // ==========================================================
  
#ifdef DEBUG
  cout << "[MuonKinVsEtaAnalyzer]: Analyze the mu in different eta regions" << endl;
#endif
  if (!muons.isValid()) return;
  
  for (reco::MuonCollection::const_iterator muonIt = muons->begin(); muonIt!=muons->end(); ++muonIt){
    reco::Muon recoMu = *muonIt;
    
    for(unsigned int iEtaRegion=0;iEtaRegion<4;iEtaRegion++){
      if (iEtaRegion==0) {EtaCutMin= etaBMin; EtaCutMax=etaBMax;}
      if (iEtaRegion==1) {EtaCutMin= etaECMin; EtaCutMax=etaECMax;}
      if (iEtaRegion==2) {EtaCutMin= etaOvlpMin; EtaCutMax=etaOvlpMax;} 
      if (iEtaRegion==3) {EtaCutMin= etaBMin; EtaCutMax=etaECMax;}
    
      if(recoMu.isGlobalMuon()) {
#ifdef DEBUG
    cout << "[MuonKinVsEtaAnalyzer]: The mu is global... Filling the histos" << endl;
#endif
    LogTrace(metname)<<"[MuonKinVsEtaAnalyzer] The mu is global - filling the histos";
    reco::TrackRef recoCombinedGlbTrack = recoMu.combinedMuon();
    // get the track combinig the information from both the glb fit"
    if(fabs(recoCombinedGlbTrack->eta())>EtaCutMin && fabs(recoCombinedGlbTrack->eta())<EtaCutMax){
      etaGlbTrack[iEtaRegion]->Fill(recoCombinedGlbTrack->eta());
      phiGlbTrack[iEtaRegion]->Fill(recoCombinedGlbTrack->phi());
      pGlbTrack[iEtaRegion]->Fill(recoCombinedGlbTrack->p());
      ptGlbTrack[iEtaRegion]->Fill(recoCombinedGlbTrack->pt());
      chi2GlbTrack[iEtaRegion]->Fill(recoCombinedGlbTrack->normalizedChi2());
      chi2probGlbTrack[iEtaRegion]->Fill(TMath::Prob(recoCombinedGlbTrack->normalizedChi2(),recoCombinedGlbTrack->ndof()));
    }
      }

      if(recoMu.isTrackerMuon()) {
#ifdef DEBUG
    cout << "[MuonKinVsEtaAnalyzer]: The mu is tracker... Filling the histos" << endl;
#endif
    LogTrace(metname)<<"[MuonKinVsEtaAnalyzer] The mu is tracker - filling the histos";
    // get the track using only the tracker data
    reco::TrackRef recoTrack = recoMu.track();
    if(fabs(recoTrack->eta())>EtaCutMin && fabs(recoTrack->eta())<EtaCutMax){
      etaTrack[iEtaRegion]->Fill(recoTrack->eta());
      phiTrack[iEtaRegion]->Fill(recoTrack->phi());
      pTrack[iEtaRegion]->Fill(recoTrack->p());
      ptTrack[iEtaRegion]->Fill(recoTrack->pt());
      chi2Track[iEtaRegion]->Fill(recoTrack->normalizedChi2());
      chi2probTrack[iEtaRegion]->Fill(TMath::Prob(recoTrack->normalizedChi2(),recoTrack->ndof()));
    }
      }

      if(recoMu.isStandAloneMuon()) {
#ifdef DEBUG
    cout << "[MuonKinVsEtaAnalyzer]: The mu is standlone... Filling the histos" << endl;
#endif
    LogTrace(metname)<<"[MuonKinVsEtaAnalyzer] The mu is standalone - filling the histos";
    // get the track using only the mu spectrometer data
    reco::TrackRef recoStaTrack = recoMu.standAloneMuon();
    if(fabs(recoStaTrack->eta())>EtaCutMin && fabs(recoStaTrack->eta())<EtaCutMax){
      etaStaTrack[iEtaRegion]->Fill(recoStaTrack->eta());
      phiStaTrack[iEtaRegion]->Fill(recoStaTrack->phi());
      pStaTrack[iEtaRegion]->Fill(recoStaTrack->p());
      ptStaTrack[iEtaRegion]->Fill(recoStaTrack->pt());
      chi2StaTrack[iEtaRegion]->Fill(recoStaTrack->normalizedChi2());
      chi2probStaTrack[iEtaRegion]->Fill(TMath::Prob(recoStaTrack->normalizedChi2(),recoStaTrack->ndof()));
    }
      }

      if ( muon::isTightMuon(recoMu, vtx) ) {
#ifdef DEBUG
    cout << "[MuonKinVsEtaAnalyzer]: The mu is tight... Filling the histos" << endl;
#endif
    LogTrace(metname)<<"[MuonKinVsEtaAnalyzer] The mu is Tight - filling the histos";
    reco::TrackRef recoTightTrack = recoMu.combinedMuon();
    if(fabs(recoTightTrack->eta())>EtaCutMin && fabs(recoTightTrack->eta())<EtaCutMax){
      etaTightTrack[iEtaRegion]->Fill(recoTightTrack->eta());
      phiTightTrack[iEtaRegion]->Fill(recoTightTrack->phi());
      pTightTrack[iEtaRegion]->Fill(recoTightTrack->p());
      ptTightTrack[iEtaRegion]->Fill(recoTightTrack->pt());
      chi2TightTrack[iEtaRegion]->Fill(recoTightTrack->normalizedChi2());
      chi2probTightTrack[iEtaRegion]->Fill(TMath::Prob(recoTightTrack->normalizedChi2(),recoTightTrack->ndof()));
    }
      }


      if ( muon::isLooseMuon(recoMu) ) {
#ifdef DEBUG
    cout << "[MuonKinVsEtaAnalyzer]: The mu is Loose... Filling the histos" << endl;
#endif
    LogTrace(metname)<<"[MuonKinVsEtaAnalyzer] The mu is Loose - filling the histos";
    reco::TrackRef recoLooseTrack; 

    if ( recoMu.isGlobalMuon()) recoLooseTrack = recoMu.combinedMuon();
    else recoLooseTrack = recoMu.track();

    if(fabs(recoLooseTrack->eta())>EtaCutMin && fabs(recoLooseTrack->eta())<EtaCutMax){
      etaLooseTrack[iEtaRegion]->Fill(recoLooseTrack->eta());
      phiLooseTrack[iEtaRegion]->Fill(recoLooseTrack->phi());
      pLooseTrack[iEtaRegion]->Fill(recoLooseTrack->p());
      ptLooseTrack[iEtaRegion]->Fill(recoLooseTrack->pt());
      chi2LooseTrack[iEtaRegion]->Fill(recoLooseTrack->normalizedChi2());
      chi2probLooseTrack[iEtaRegion]->Fill(TMath::Prob(recoLooseTrack->normalizedChi2(),recoLooseTrack->ndof()));
    }
      }

      if ( muon::isMediumMuon(recoMu) ) {
#ifdef DEBUG
    cout << "[MuonKinVsEtaAnalyzer]: The mu is Medium... Filling the histos" << endl;
#endif
    LogTrace(metname)<<"[MuonKinVsEtaAnalyzer] The mu is Medium - filling the histos";
    reco::TrackRef recoMediumTrack; 

    if ( recoMu.isGlobalMuon()) recoMediumTrack = recoMu.combinedMuon();
    else recoMediumTrack = recoMu.track();

    if(fabs(recoMediumTrack->eta())>EtaCutMin && fabs(recoMediumTrack->eta())<EtaCutMax){
      etaMediumTrack[iEtaRegion]->Fill(recoMediumTrack->eta());
      phiMediumTrack[iEtaRegion]->Fill(recoMediumTrack->phi());
      pMediumTrack[iEtaRegion]->Fill(recoMediumTrack->p());
      ptMediumTrack[iEtaRegion]->Fill(recoMediumTrack->pt());
      chi2MediumTrack[iEtaRegion]->Fill(recoMediumTrack->normalizedChi2());
      chi2probMediumTrack[iEtaRegion]->Fill(TMath::Prob(recoMediumTrack->normalizedChi2(),recoMediumTrack->ndof()));
    }
      }

      if ( muon::isSoftMuon(recoMu, vtx) ) {
#ifdef DEBUG
    cout << "[MuonKinVsEtaAnalyzer]: The mu is Soft... Filling the histos" << endl;
#endif
    LogTrace(metname)<<"[MuonKinVsEtaAnalyzer] The mu is Soft - filling the histos";
    reco::TrackRef recoSoftTrack = recoMu.track();
    if(fabs(recoSoftTrack->eta())>EtaCutMin && fabs(recoSoftTrack->eta())<EtaCutMax){
      etaSoftTrack[iEtaRegion]->Fill(recoSoftTrack->eta());
      phiSoftTrack[iEtaRegion]->Fill(recoSoftTrack->phi());
      pSoftTrack[iEtaRegion]->Fill(recoSoftTrack->p());
      ptSoftTrack[iEtaRegion]->Fill(recoSoftTrack->pt());
      chi2SoftTrack[iEtaRegion]->Fill(recoSoftTrack->normalizedChi2());
      chi2probSoftTrack[iEtaRegion]->Fill(TMath::Prob(recoSoftTrack->normalizedChi2(),recoSoftTrack->ndof()));
    }
      }

      if ( muon::isHighPtMuon(recoMu, vtx) ) {
#ifdef DEBUG
    cout << "[MuonKinVsEtaAnalyzer]: The mu is HighPt... Filling the histos" << endl;
#endif
    LogTrace(metname)<<"[MuonKinVsEtaAnalyzer] The mu is HightPt - filling the histos";
    reco::TrackRef recoHighPtTrack = recoMu.combinedMuon();
    if(fabs(recoHighPtTrack->eta())>EtaCutMin && fabs(recoHighPtTrack->eta())<EtaCutMax){
      etaHighPtTrack[iEtaRegion]->Fill(recoHighPtTrack->eta());
      phiHighPtTrack[iEtaRegion]->Fill(recoHighPtTrack->phi());
      pHighPtTrack[iEtaRegion]->Fill(recoHighPtTrack->p());
      ptHighPtTrack[iEtaRegion]->Fill(recoHighPtTrack->pt());
      chi2HighPtTrack[iEtaRegion]->Fill(recoHighPtTrack->normalizedChi2());
      chi2probHighPtTrack[iEtaRegion]->Fill(TMath::Prob(recoHighPtTrack->normalizedChi2(),recoHighPtTrack->ndof()));
    }
      }
    } //end iEtaRegions
  } //end recoMu iteration
}

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 68 of file MuonKinVsEtaAnalyzer.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::cd(), diMuonHistograms_cfi::etaBin, HLT_25ns10e33_v2_cff::etaMax, HLT_25ns10e33_v2_cff::etaMin, HLT_25ns10e33_v2_cff::pMin, PtMinSelector_cfg::ptMin, DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                             {
  ibooker.cd();
  ibooker.setCurrentFolder("Muons/MuonKinVsEtaAnalyzer");


  std::string EtaName;
  for(unsigned int iEtaRegion=0;iEtaRegion<4;iEtaRegion++){
    if (iEtaRegion==0) EtaName = "Barrel";   
    if (iEtaRegion==1) EtaName = "EndCap";   
    if (iEtaRegion==2) EtaName = "Overlap";
    if (iEtaRegion==3) EtaName = "";

    // monitoring of eta parameter
    etaGlbTrack.push_back(ibooker.book1D("GlbMuon_eta_"+EtaName, "#eta_{GLB} "+EtaName, etaBin, etaMin, etaMax));
    etaTrack.push_back(ibooker.book1D("TkMuon_eta_"+EtaName, "#eta_{TK} "+EtaName, etaBin, etaMin, etaMax));
    etaStaTrack.push_back(ibooker.book1D("StaMuon_eta_"+EtaName, "#eta_{STA} "+EtaName, etaBin, etaMin, etaMax));
    etaTightTrack.push_back(ibooker.book1D("TightMuon_eta_"+EtaName, "#eta_{Tight} "+EtaName, etaBin, etaMin, etaMax));
    etaLooseTrack.push_back(ibooker.book1D("LooseMuon_eta_"+EtaName, "#eta_{Loose} "+EtaName, etaBin, etaMin, etaMax));
    etaMediumTrack.push_back(ibooker.book1D("MediumMuon_eta_"+EtaName, "#eta_{Medium} "+EtaName, etaBin, etaMin, etaMax));
    etaSoftTrack.push_back(ibooker.book1D("SoftMuon_eta_"+EtaName, "#eta_{Soft} "+EtaName, etaBin, etaMin, etaMax));
    etaHighPtTrack.push_back(ibooker.book1D("HighPtMuon_eta_"+EtaName, "#eta_{HighPt} "+EtaName, etaBin, etaMin, etaMax));
    
    // monitoring of phi paramater
    phiGlbTrack.push_back(ibooker.book1D("GlbMuon_phi_"+EtaName, "#phi_{GLB} "+EtaName+ "(rad)", phiBin, phiMin, phiMax));
    phiTrack.push_back(ibooker.book1D("TkMuon_phi_"+EtaName, "#phi_{TK}" +EtaName +"(rad)", phiBin, phiMin, phiMax));
    phiStaTrack.push_back(ibooker.book1D("StaMuon_phi_"+EtaName, "#phi_{STA}"+EtaName+" (rad)", phiBin, phiMin, phiMax));
    phiTightTrack.push_back(ibooker.book1D("TightMuon_phi_"+EtaName, "#phi_{Tight}_"+EtaName, phiBin, phiMin, phiMax));
    phiLooseTrack.push_back(ibooker.book1D("LooseMuon_phi_"+EtaName, "#phi_{Loose}_"+EtaName, phiBin, phiMin, phiMax));
    phiMediumTrack.push_back(ibooker.book1D("MediumMuon_phi_"+EtaName, "#phi_{Medium}_"+EtaName, phiBin, phiMin, phiMax));
    phiSoftTrack.push_back(ibooker.book1D("SoftMuon_phi_"+EtaName, "#phi_{Soft}_"+EtaName, phiBin, phiMin, phiMax));
    phiHighPtTrack.push_back(ibooker.book1D("HighPtMuon_phi_"+EtaName, "#phi_{HighPt}_"+EtaName, phiBin, phiMin, phiMax));
    
    // monitoring of the momentum
    pGlbTrack.push_back(ibooker.book1D("GlbMuon_p_"+EtaName, "p_{GLB} "+EtaName, pBin, pMin, pMax));
    pTrack.push_back(ibooker.book1D("TkMuon_p"+EtaName, "p_{TK} "+EtaName, pBin, pMin, pMax));
    pStaTrack.push_back(ibooker.book1D("StaMuon_p"+EtaName, "p_{STA} "+EtaName, pBin, pMin, pMax));
    pTightTrack.push_back(ibooker.book1D("TightMuon_p_"+EtaName, "p_{Tight} "+EtaName, pBin, pMin, pMax));
    pLooseTrack.push_back(ibooker.book1D("LooseMuon_p_"+EtaName, "p_{Loose} "+EtaName, pBin, pMin, pMax));
    pMediumTrack.push_back(ibooker.book1D("MediumMuon_p_"+EtaName, "p_{Medium} "+EtaName, pBin, pMin, pMax));
    pSoftTrack.push_back(ibooker.book1D("SoftMuon_p_"+EtaName, "p_{Soft} "+EtaName, pBin, pMin, pMax));
    pHighPtTrack.push_back(ibooker.book1D("HighPtMuon_p_"+EtaName, "p_{HighPt} "+EtaName, pBin, pMin, pMax));

    // monitoring of the transverse momentum
    ptGlbTrack.push_back(ibooker.book1D("GlbMuon_pt_" +EtaName, "pt_{GLB} "+EtaName, ptBin, ptMin, ptMax));
    ptTrack.push_back(ibooker.book1D("TkMuon_pt_"+EtaName, "pt_{TK} "+EtaName, ptBin, ptMin, ptMax));
    ptStaTrack.push_back(ibooker.book1D("StaMuon_pt_"+EtaName, "pt_{STA} "+EtaName, ptBin, ptMin, pMax));
    ptTightTrack.push_back(ibooker.book1D("TightMuon_pt_"+EtaName, "pt_{Tight} "+EtaName, ptBin, ptMin, ptMax));
    ptLooseTrack.push_back(ibooker.book1D("LooseMuon_pt_"+EtaName, "pt_{Loose} "+EtaName, ptBin, ptMin, ptMax));
    ptMediumTrack.push_back(ibooker.book1D("MediumMuon_pt_"+EtaName, "pt_{Medium} "+EtaName, ptBin, ptMin, ptMax));
    ptSoftTrack.push_back(ibooker.book1D("SoftMuon_pt_"+EtaName, "pt_{Soft} "+EtaName, ptBin, ptMin, ptMax));
    ptHighPtTrack.push_back(ibooker.book1D("HighPtMuon_pt_"+EtaName, "pt_{HighPt} "+EtaName, ptBin, ptMin, ptMax));

    // monitoring chi2 and Prob.Chi2
    chi2GlbTrack.push_back(ibooker.book1D("GlbMuon_chi2_"+EtaName, "#chi^{2}_{GLB} " + EtaName, chiBin, chiMin, chiMax));
    chi2probGlbTrack.push_back(ibooker.book1D("GlbMuon_chi2prob_"+EtaName, "#chi^{2}_{GLB} prob." + EtaName, chiBin, chiprobMin, chiprobMax));
    chi2Track.push_back(ibooker.book1D("TkMuon_chi2_"+EtaName, "#chi^{2}_{TK} " + EtaName, chiBin, chiMin, chiMax));
    chi2probTrack.push_back(ibooker.book1D("TkMuon_chi2prob_"+EtaName, "#chi^{2}_{TK} prob." + EtaName, chiBin, chiprobMin, chiprobMax));
    chi2StaTrack.push_back(ibooker.book1D("StaMuon_chi2_"+EtaName, "#chi^{2}_{STA} " + EtaName, chiBin, chiMin, chiMax));
    chi2probStaTrack.push_back(ibooker.book1D("StaMuon_chi2prob_"+EtaName, "#chi^{2}_{STA} prob." + EtaName, chiBin, chiprobMin, chiprobMax));
    chi2TightTrack.push_back(ibooker.book1D("TightMuon_chi2_"+EtaName, "#chi^{2}_{Tight} " + EtaName, chiBin, chiMin, chiMax));
    chi2probTightTrack.push_back(ibooker.book1D("TightMuon_chi2prob_"+EtaName, "#chi^{2}_{Tight} prob." + EtaName, chiBin, chiprobMin, chiprobMax));
    chi2LooseTrack.push_back(ibooker.book1D("LooseMuon_chi2_"+EtaName, "#chi^{2}_{Loose} " + EtaName, chiBin, chiMin, chiMax));
    chi2MediumTrack.push_back(ibooker.book1D("MediumMuon_chi2_"+EtaName, "#chi^{2}_{Medium} " + EtaName, chiBin, chiMin, chiMax));
    chi2probLooseTrack.push_back(ibooker.book1D("LooseMuon_chi2prob_"+EtaName, "#chi^{2}_{Loose} prob." + EtaName, chiBin, chiprobMin, chiprobMax));
    chi2probMediumTrack.push_back(ibooker.book1D("MediumMuon_chi2prob_"+EtaName, "#chi^{2}_{Medium} prob." + EtaName, chiBin, chiprobMin, chiprobMax));
    chi2SoftTrack.push_back(ibooker.book1D("SoftMuon_chi2_"+EtaName, "#chi^{2}_{Soft} " + EtaName, chiBin, chiMin, chiMax));
    chi2probSoftTrack.push_back(ibooker.book1D("SoftMuon_chi2prob_"+EtaName, "#chi^{2}_{Soft} prob." + EtaName, chiBin, chiprobMin, chiprobMax));
    chi2HighPtTrack.push_back(ibooker.book1D("HighPtMuon_chi2_"+EtaName, "#chi^{2}_{HighPt} " + EtaName, chiBin, chiMin, chiMax));
    chi2probHighPtTrack.push_back(ibooker.book1D("HighPtMuon_chi2prob_"+EtaName, "#chi^{2}_{HighPt} prob." + EtaName, chiBin, chiprobMin, chiprobMax));
  }
}

Member Data Documentation

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2GlbTrack

private

Definition at line 101 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2HighPtTrack

private

Definition at line 157 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2LooseTrack

private

Definition at line 133 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2MediumTrack

private

Definition at line 141 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2probGlbTrack

private

Definition at line 102 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2probHighPtTrack

private

Definition at line 158 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2probLooseTrack

private

Definition at line 134 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2probMediumTrack

private

Definition at line 142 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2probSoftTrack

private

Definition at line 150 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2probStaTrack

private

Definition at line 118 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2probTightTrack

private

Definition at line 126 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2probTrack

private

Definition at line 110 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2SoftTrack

private

Definition at line 149 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2StaTrack

private

Definition at line 117 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2TightTrack

private

Definition at line 125 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::chi2Track

private

Definition at line 109 of file MuonKinVsEtaAnalyzer.h.

int MuonKinVsEtaAnalyzer::chiBin

private

Definition at line 78 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::chiMax

private

Definition at line 80 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::chiMin

private

Definition at line 79 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::chiprobMax

private

Definition at line 83 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::chiprobMin

private

Definition at line 82 of file MuonKinVsEtaAnalyzer.h.

int MuonKinVsEtaAnalyzer::etaBin

private

Definition at line 70 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::etaBMax

private

Definition at line 89 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::etaBMin

private

Definition at line 88 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::EtaCutMax

private

Definition at line 87 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::EtaCutMin

private

Definition at line 86 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::etaECMax

private

Definition at line 91 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::etaECMin

private

Definition at line 90 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::etaGlbTrack

private

Definition at line 97 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::etaHighPtTrack

private

Definition at line 153 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::etaLooseTrack

private

Definition at line 129 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::etaMax

private

Definition at line 72 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::etaMediumTrack

private

Definition at line 137 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::etaMin

private

Definition at line 71 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::etaOvlpMax

private

Definition at line 93 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::etaOvlpMin

private

Definition at line 92 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::etaSoftTrack

private

Definition at line 145 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::etaStaTrack

private

Definition at line 113 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::etaTightTrack

private

Definition at line 121 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::etaTrack

private

Definition at line 105 of file MuonKinVsEtaAnalyzer.h.

std::string MuonKinVsEtaAnalyzer::metname

private

Definition at line 54 of file MuonKinVsEtaAnalyzer.h.

edm::ParameterSet MuonKinVsEtaAnalyzer::parameters

private

Definition at line 51 of file MuonKinVsEtaAnalyzer.h.

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

int MuonKinVsEtaAnalyzer::pBin

private

Definition at line 62 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::pGlbTrack

private

Definition at line 99 of file MuonKinVsEtaAnalyzer.h.

int MuonKinVsEtaAnalyzer::phiBin

private

Definition at line 74 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::pHighPtTrack

private

Definition at line 155 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::phiGlbTrack

private

Definition at line 98 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::phiHighPtTrack

private

Definition at line 154 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::phiLooseTrack

private

Definition at line 130 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::phiMax

private

Definition at line 76 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::phiMediumTrack

private

Definition at line 138 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::phiMin

private

Definition at line 75 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::phiSoftTrack

private

Definition at line 146 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::phiStaTrack

private

Definition at line 114 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::phiTightTrack

private

Definition at line 122 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::phiTrack

private

Definition at line 106 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::pLooseTrack

private

Definition at line 131 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::pMax

private

Definition at line 64 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::pMediumTrack

private

Definition at line 139 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::pMin

private

Definition at line 63 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::pSoftTrack

private

Definition at line 147 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::pStaTrack

private

Definition at line 115 of file MuonKinVsEtaAnalyzer.h.

int MuonKinVsEtaAnalyzer::ptBin

private

Definition at line 66 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::ptGlbTrack

private

Definition at line 100 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::ptHighPtTrack

private

Definition at line 156 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::pTightTrack

private

Definition at line 123 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::ptLooseTrack

private

Definition at line 132 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::ptMax

private

Definition at line 68 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::ptMediumTrack

private

Definition at line 140 of file MuonKinVsEtaAnalyzer.h.

double MuonKinVsEtaAnalyzer::ptMin

private

Definition at line 67 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::pTrack

private

Definition at line 107 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::ptSoftTrack

private

Definition at line 148 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::ptStaTrack

private

Definition at line 116 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::ptTightTrack

private

Definition at line 124 of file MuonKinVsEtaAnalyzer.h.

std::vector<MonitorElement*> MuonKinVsEtaAnalyzer::ptTrack

private

Definition at line 108 of file MuonKinVsEtaAnalyzer.h.

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

private

Definition at line 59 of file MuonKinVsEtaAnalyzer.h.

DQMStore* MuonKinVsEtaAnalyzer::theDbe

private

Definition at line 50 of file MuonKinVsEtaAnalyzer.h.

edm::EDGetTokenT<reco::MuonCollection> MuonKinVsEtaAnalyzer::theMuonCollectionLabel_

private

Definition at line 57 of file MuonKinVsEtaAnalyzer.h.

MuonServiceProxy* MuonKinVsEtaAnalyzer::theService

private

Definition at line 49 of file MuonKinVsEtaAnalyzer.h.

edm::EDGetTokenT<reco::VertexCollection> MuonKinVsEtaAnalyzer::theVertexLabel_

private

Definition at line 58 of file MuonKinVsEtaAnalyzer.h.