#include <MuonTiming.h>

Inheritance diagram for MuonTiming:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	Inizialize parameters for histo binning. More...

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

	MuonTiming (const edm::ParameterSet &)
	Constructor. More...

	~MuonTiming () override
	Destructor. More...

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

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

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

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer ()

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

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Types
enum	eta_ { overlap, barrel, endcap }

enum	object_ { sta, glb }

Private Attributes
double	etaBarrelMax_

double	etaBarrelMin_

double	etaEndcapMax_

double	etaEndcapMin_

double	etamax_

double	etamin_

std::vector< std::string >	EtaName_

int	etanbins_

double	etaOverlapMax_

double	etaOverlapMin_

std::array< MonitorElement *, 1 >	etaptVeto_

std::array< MonitorElement *, 1 >	etaVeto_

std::string	metname_

double	ndofmax_

double	ndofmin_

int	ndofnbins_

std::vector< std::string >	ObjectName_

double	ptmax_

double	ptmin_

int	ptnbins_

std::array< MonitorElement *, 1 >	ptVeto_

double	terrmax_

double	terrmaxrpc_

double	terrmin_

double	terrminrpc_

int	terrnbins_

int	terrnbinsrpc_

std::string	theFolder_

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

std::array< std::array < MonitorElement *, 1 >, 3 >	timeAtIpInOut_

std::array< std::array < MonitorElement *, 1 >, 3 >	timeAtIpInOutErr_

std::array< std::array < MonitorElement *, 1 >, 3 >	timeAtIpInOutErrRPC_

std::array< std::array < MonitorElement *, 1 >, 3 >	timeAtIpInOutRPC_

std::array< std::array < MonitorElement *, 1 >, 3 >	timeNDof_

double	tmax_

double	tmaxrpc_

double	tmin_

double	tminrpc_

int	tnbins_

int	tnbinsrpc_

std::array< MonitorElement *, 1 >	yields_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr < DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 26 of file MuonTiming.h.

Member Enumeration Documentation

enum MuonTiming::eta_

private

Enumerator
overlap
barrel
endcap

Definition at line 64 of file MuonTiming.h.

64 { overlap, barrel, endcap };

MuonTiming::endcap

Definition: MuonTiming.h:64

MuonTiming::barrel

Definition: MuonTiming.h:64

MuonTiming::overlap

Definition: MuonTiming.h:64

enum MuonTiming::object_

private

Enumerator
sta
glb

Definition at line 65 of file MuonTiming.h.

65 { sta, glb };

MuonTiming::sta

Definition: MuonTiming.h:65

MuonTiming::glb

Definition: MuonTiming.h:65

Constructor & Destructor Documentation

MuonTiming::MuonTiming ( const edm::ParameterSet & pSet )

Constructor.

Definition at line 18 of file MuonTiming.cc.

References edm::ParameterSet::getParameter().

                                                   {
   const edm::ParameterSet& parameters = pSet;
 
   // Input booleans
 
   // the services:
   theMuonCollectionLabel_ = consumes<edm::View<reco::Muon> >(parameters.getParameter<edm::InputTag>("MuonCollection"));
 
   tnbins_ = parameters.getParameter<int>("tnbins");
   tnbinsrpc_ = parameters.getParameter<int>("tnbinsrpc");
   terrnbins_ = parameters.getParameter<int>("terrnbins");
   terrnbinsrpc_ = parameters.getParameter<int>("terrnbinsrpc");
   ndofnbins_ = parameters.getParameter<int>("ndofnbins");
   ptnbins_ = parameters.getParameter<int>("ptnbins");
   etanbins_ = parameters.getParameter<int>("etanbins");
   tmax_ = parameters.getParameter<double>("tmax");
   tmaxrpc_ = parameters.getParameter<double>("tmaxrpc");
   terrmax_ = parameters.getParameter<double>("terrmax");
   terrmaxrpc_ = parameters.getParameter<double>("terrmaxrpc");
   ndofmax_ = parameters.getParameter<double>("ndofmax");
   ptmax_ = parameters.getParameter<double>("ptmax");
   etamax_ = parameters.getParameter<double>("etamax");
   tmin_ = parameters.getParameter<double>("tmin");
   tminrpc_ = parameters.getParameter<double>("tminrpc");
   terrmin_ = parameters.getParameter<double>("terrmin");
   terrminrpc_ = parameters.getParameter<double>("terrminrpc");
   ndofmin_ = parameters.getParameter<double>("ndofmin");
   ptmin_ = parameters.getParameter<double>("ptmin");
   etamin_ = parameters.getParameter<double>("etamin");
 
   etaBarrelMin_ = parameters.getParameter<double>("etaBarrelMin");
   etaBarrelMax_ = parameters.getParameter<double>("etaBarrelMax");
   etaEndcapMin_ = parameters.getParameter<double>("etaEndcapMin");
   etaEndcapMax_ = parameters.getParameter<double>("etaEndcapMax");
   etaOverlapMin_ = parameters.getParameter<double>("etaOverlapMin");
   etaOverlapMax_ = parameters.getParameter<double>("etaOverlapMax");
 
   theFolder_ = parameters.getParameter<string>("folder");
 }

MuonTiming::~MuonTiming ( )

override

Destructor.

Definition at line 58 of file MuonTiming.cc.

58 {}

Member Function Documentation

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

overridevirtual

Inizialize parameters for histo binning.

Reimplemented from DQMEDAnalyzer.

Definition at line 153 of file MuonTiming.cc.

References Reference_intrackfit_cff::barrel, Reference_intrackfit_cff::endcap, HcalObjRepresent::Fill(), edm::Event::getByToken(), glb, edm::HandleBase::isValid(), LogTrace, HLT_FULL_cff::muon, patZpeak::muons, reco::MuonTime::nDof, muon::overlap(), reco::MuonTime::timeAtIpInOut, reco::MuonTime::timeAtIpInOutErr, HLT_FULL_cff::track, and relval_cleanedupgrade::veto.

                                                                             {
   LogTrace(metname_) << "[MuonTiming] Analyze the mu";
 
   // Take the muon container
   edm::Handle<edm::View<reco::Muon> > muons;
   iEvent.getByToken(theMuonCollectionLabel_, muons);
 
   if (!muons.isValid())
     return;
 
   for (edm::View<reco::Muon>::const_iterator muon = muons->begin(); muon != muons->end(); ++muon) {
     const reco::MuonTime time = muon->time();
     const reco::MuonTime rpcTime = muon->rpcTime();
     //Only creating so far the timing information for STA muons
     if (!muon->isStandAloneMuon() || muon->isGlobalMuon())
       continue;
     reco::TrackRef track;
     //Select whether it's a global or standalone muon
     object_ theObject = sta;
     if (muon->isGlobalMuon()) {
       track = muon->combinedMuon();
       theObject = glb;
     } else {
       track = muon->standAloneMuon();
       theObject = sta;
     }
 
     //These definitions have been taken from Piotr Traczyk
     bool cmbok = (time.nDof > 7);
     bool rpcok = (rpcTime.nDof > 1 && rpcTime.timeAtIpInOutErr == 0);
     bool veto = false;
     if (rpcok) {
       if ((fabs(rpcTime.timeAtIpInOut) > 10) && !(cmbok && fabs(time.timeAtIpInOut) < 10))
         veto = true;
       else if (cmbok && (time.timeAtIpInOut > 20 || time.timeAtIpInOut < -45))
         veto = true;
     }
 
     //std::cout << time.timeAtIpInOut << std::endl;
     //Filling the yields histogram
     if (muon->isTimeValid())
       yields_[theObject]->Fill(1);
     else
       yields_[theObject]->Fill(0);
 
     if (cmbok)
       yields_[theObject]->Fill(3);
     else
       yields_[theObject]->Fill(2);
 
     if (rpcok)
       yields_[theObject]->Fill(5);
     else
       yields_[theObject]->Fill(4);
 
     if (cmbok && !rpcok)
       yields_[theObject]->Fill(6);
     if (!cmbok && rpcok)
       yields_[theObject]->Fill(7);
 
     if (veto)
       yields_[theObject]->Fill(8);
     else
       yields_[theObject]->Fill(9);
 
     //Starting now with the pt and eta for vetoed and not vetoed muons
     if (veto) {
       etaptVeto_[theObject]->Fill(track->pt(), track->eta());
       etaVeto_[theObject]->Fill(track->eta());
       ptVeto_[theObject]->Fill(track->pt());
     }
 
     //Check the eta region of the muon
     eta_ theEta = barrel;
     if (fabs(track->eta()) >= etaBarrelMin_ && fabs(track->eta()) <= etaBarrelMax_)
       theEta = barrel;
     if (fabs(track->eta()) >= etaOverlapMin_ && fabs(track->eta()) <= etaOverlapMax_)
       theEta = overlap;
     if (fabs(track->eta()) >= etaEndcapMin_ && fabs(track->eta()) <= etaEndcapMax_)
       theEta = endcap;
     timeNDof_[theEta][theObject]->Fill(time.nDof);
     timeAtIpInOut_[theEta][theObject]->Fill(time.timeAtIpInOut);
     timeAtIpInOutRPC_[theEta][theObject]->Fill(rpcTime.timeAtIpInOut);
     timeAtIpInOutErr_[theEta][theObject]->Fill(time.timeAtIpInOutErr);
     timeAtIpInOutErrRPC_[theEta][theObject]->Fill(rpcTime.timeAtIpInOutErr);
   }
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 60 of file MuonTiming.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::NavigatorBase::cd(), dqm::impl::MonitorElement::setAxisTitle(), dqm::impl::MonitorElement::setBinLabel(), and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                         {
   ibooker.cd();
   ibooker.setCurrentFolder(theFolder_);
 
   EtaName_.push_back("_Overlap");
   EtaName_.push_back("_Barrel");
   EtaName_.push_back("_Endcap");
   ObjectName_.push_back("Sta_");
   ObjectName_.push_back("Glb_");
 
   for (unsigned int iEtaRegion = 0; iEtaRegion < 3; iEtaRegion++) {
     /*std::array<MonitorElement*, 1> timeNDofv_;
     std::array<MonitorElement*, 1> timeAtIpInOutv_;
     std::array<MonitorElement*, 1> timeAtIpInOutRPCv_;
     std::array<MonitorElement*, 1> timeAtIpInOutErrv_;
     std::array<MonitorElement*, 1> timeAtIpInOutErrRPCv_;*/
     //Only creating so far the timing information for STA muons, however the code can be extended to also Glb by just setting the limit of this loop to 2
     for (unsigned int iObjectName = 0; iObjectName < 1; iObjectName++) {
       timeNDof_[iEtaRegion][iObjectName] = ibooker.book1D(
           ObjectName_[iObjectName] + "timenDOF" + EtaName_[iEtaRegion], "muon time ndof", ndofnbins_, 0, ndofmax_);
       timeAtIpInOut_[iEtaRegion][iObjectName] = ibooker.book1D(
           ObjectName_[iObjectName] + "timeAtIpInOut" + EtaName_[iEtaRegion], "muon time", tnbins_, tmin_, tmax_);
       timeAtIpInOutRPC_[iEtaRegion][iObjectName] =
           ibooker.book1D(ObjectName_[iObjectName] + "timeAtIpInOutRPC" + EtaName_[iEtaRegion],
                          "muon rpc time",
                          tnbinsrpc_,
                          tminrpc_,
                          tmaxrpc_);
       timeAtIpInOutErr_[iEtaRegion][iObjectName] =
           ibooker.book1D(ObjectName_[iObjectName] + "timeAtIpInOutErr" + EtaName_[iEtaRegion],
                          "muon time error",
                          terrnbins_,
                          terrmin_,
                          terrmax_);
       timeAtIpInOutErrRPC_[iEtaRegion][iObjectName] =
           ibooker.book1D(ObjectName_[iObjectName] + "timeAtIpInOutRPCErr" + EtaName_[iEtaRegion],
                          "muon rpc time error",
                          terrnbinsrpc_,
                          terrminrpc_,
                          terrmaxrpc_);
       timeNDof_[iEtaRegion][iObjectName]->setAxisTitle("Time nDof");
       timeAtIpInOut_[iEtaRegion][iObjectName]->setAxisTitle("Combined time [ns]");
       timeAtIpInOutErr_[iEtaRegion][iObjectName]->setAxisTitle("Combined time Error [ns]");
       timeAtIpInOutRPC_[iEtaRegion][iObjectName]->setAxisTitle("RPC time [ns]");
       timeAtIpInOutErrRPC_[iEtaRegion][iObjectName]->setAxisTitle("RPC time Error [ns]");
     }
     /*
     timeNDof_[iEtaregion] = timeNDofv_);
     timeAtIpInOut_.push_back(timeAtIpInOutv_);
     timeAtIpInOutRPC_.push_back(timeAtIpInOutRPCv_);
     timeAtIpInOutErr_.push_back(timeAtIpInOutErrv_);
     timeAtIpInOutErrRPC_.push_back(timeAtIpInOutErrRPCv_);
     */
   }
 
   //Only creating so far the timing information for STA muons, however the code can be extended to also Glb by just setting the limit of this loop to 2
   for (unsigned int iObjectName = 0; iObjectName < 1; iObjectName++) {
     etaptVeto_[iObjectName] = ibooker.book2D(ObjectName_[iObjectName] + "etapt",
                                              "Eta and Pt distribution for muons not passing the veto",
                                              ptnbins_,
                                              ptmin_,
                                              ptmax_,
                                              etanbins_,
                                              etamin_,
                                              etamax_);
     etaVeto_[iObjectName] = ibooker.book1D(ObjectName_[iObjectName] + "eta",
                                            "Eta distribution for muons not passing the veto",
                                            etanbins_,
                                            etamin_,
                                            etamax_);
     ptVeto_[iObjectName] = ibooker.book1D(
         ObjectName_[iObjectName] + "pt", "Pt distribution for muons not passing the veto", ptnbins_, ptmin_, ptmax_);
     yields_[iObjectName] = ibooker.book1D(
         ObjectName_[iObjectName] + "yields", "Number of muons passing/not passing the different conditions", 10, 0, 10);
     yields_[iObjectName]->setBinLabel(1, "Not valid time");
     yields_[iObjectName]->setBinLabel(2, "Valid time");
     yields_[iObjectName]->setBinLabel(3, "Not Combined time");
     yields_[iObjectName]->setBinLabel(4, "Combined time");
     yields_[iObjectName]->setBinLabel(5, "Not RPC time");
     yields_[iObjectName]->setBinLabel(6, "RPC time");
     yields_[iObjectName]->setBinLabel(7, "Combined not RPC");
     yields_[iObjectName]->setBinLabel(8, "RPC not Combined");
     yields_[iObjectName]->setBinLabel(9, "Not passing veto");
     yields_[iObjectName]->setBinLabel(10, "Passing veto");
     etaptVeto_[iObjectName]->setAxisTitle("p_{T} [GeV]");
     etaptVeto_[iObjectName]->setAxisTitle("#eta#", 2);
     ptVeto_[iObjectName]->setAxisTitle("p_{T} [GeV]");
     etaVeto_[iObjectName]->setAxisTitle("#eta");
   }
 }

Member Data Documentation

double MuonTiming::etaBarrelMax_

private

Definition at line 60 of file MuonTiming.h.

double MuonTiming::etaBarrelMin_

private

Definition at line 60 of file MuonTiming.h.

double MuonTiming::etaEndcapMax_

private

Definition at line 60 of file MuonTiming.h.

double MuonTiming::etaEndcapMin_

private

Definition at line 60 of file MuonTiming.h.

double MuonTiming::etamax_

private

Definition at line 59 of file MuonTiming.h.

double MuonTiming::etamin_

private

Definition at line 59 of file MuonTiming.h.

std::vector<std::string> MuonTiming::EtaName_

private

Definition at line 63 of file MuonTiming.h.

int MuonTiming::etanbins_

private

Definition at line 52 of file MuonTiming.h.

double MuonTiming::etaOverlapMax_

private

Definition at line 60 of file MuonTiming.h.

double MuonTiming::etaOverlapMin_

private

Definition at line 60 of file MuonTiming.h.

std::array<MonitorElement *, 1> MuonTiming::etaptVeto_

private

Definition at line 83 of file MuonTiming.h.

std::array<MonitorElement *, 1> MuonTiming::etaVeto_

private

Definition at line 84 of file MuonTiming.h.

std::string MuonTiming::metname_

private

Definition at line 43 of file MuonTiming.h.

double MuonTiming::ndofmax_

private

Definition at line 57 of file MuonTiming.h.

double MuonTiming::ndofmin_

private

Definition at line 57 of file MuonTiming.h.

int MuonTiming::ndofnbins_

private

Definition at line 50 of file MuonTiming.h.

std::vector<std::string> MuonTiming::ObjectName_

private

Definition at line 63 of file MuonTiming.h.

double MuonTiming::ptmax_

private

Definition at line 58 of file MuonTiming.h.

double MuonTiming::ptmin_

private

Definition at line 58 of file MuonTiming.h.

int MuonTiming::ptnbins_

private

Definition at line 51 of file MuonTiming.h.

std::array<MonitorElement *, 1> MuonTiming::ptVeto_

private

Definition at line 85 of file MuonTiming.h.

double MuonTiming::terrmax_

private

Definition at line 55 of file MuonTiming.h.

double MuonTiming::terrmaxrpc_

private

Definition at line 56 of file MuonTiming.h.

double MuonTiming::terrmin_

private

Definition at line 55 of file MuonTiming.h.

double MuonTiming::terrminrpc_

private

Definition at line 56 of file MuonTiming.h.

int MuonTiming::terrnbins_

private

Definition at line 48 of file MuonTiming.h.

int MuonTiming::terrnbinsrpc_

private

Definition at line 49 of file MuonTiming.h.

std::string MuonTiming::theFolder_

private

Definition at line 62 of file MuonTiming.h.

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

private

Definition at line 41 of file MuonTiming.h.

std::array<std::array<MonitorElement *, 1>, 3> MuonTiming::timeAtIpInOut_

private

Definition at line 79 of file MuonTiming.h.

std::array<std::array<MonitorElement *, 1>, 3> MuonTiming::timeAtIpInOutErr_

private

Definition at line 81 of file MuonTiming.h.

std::array<std::array<MonitorElement *, 1>, 3> MuonTiming::timeAtIpInOutErrRPC_

private

Definition at line 82 of file MuonTiming.h.

std::array<std::array<MonitorElement *, 1>, 3> MuonTiming::timeAtIpInOutRPC_

private

Definition at line 80 of file MuonTiming.h.

std::array<std::array<MonitorElement *, 1>, 3> MuonTiming::timeNDof_

private

Definition at line 78 of file MuonTiming.h.

double MuonTiming::tmax_

private

Definition at line 53 of file MuonTiming.h.

double MuonTiming::tmaxrpc_

private

Definition at line 54 of file MuonTiming.h.

double MuonTiming::tmin_

private

Definition at line 53 of file MuonTiming.h.

double MuonTiming::tminrpc_

private

Definition at line 54 of file MuonTiming.h.

int MuonTiming::tnbins_

private

Definition at line 46 of file MuonTiming.h.

int MuonTiming::tnbinsrpc_

private

Definition at line 47 of file MuonTiming.h.

std::array<MonitorElement *, 1> MuonTiming::yields_

private

Definition at line 86 of file MuonTiming.h.

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation