#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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

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_

MuonServiceProxy *	theService

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_

Detailed Description

Definition at line 29 of file MuonTiming.h.

Member Enumeration Documentation

enum MuonTiming::eta_

private

Enumerator
overlap
barrel
endcap

Definition at line 70 of file MuonTiming.h.

70 {overlap, barrel, endcap};

MuonTiming::endcap

Definition: MuonTiming.h:70

MuonTiming::barrel

Definition: MuonTiming.h:70

MuonTiming::overlap

Definition: MuonTiming.h:70

enum MuonTiming::object_

private

Enumerator
sta
glb

Definition at line 71 of file MuonTiming.h.

71 {sta, glb};

MuonTiming::sta

Definition: MuonTiming.h:71

MuonTiming::glb

Definition: MuonTiming.h:71

Constructor & Destructor Documentation

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

Constructor.

Definition at line 20 of file MuonTiming.cc.

References edm::ParameterSet::getParameter(), and MuonServiceProxy_cff::MuonServiceProxy.

                                                   {
   const edm::ParameterSet& parameters = pSet;
 
   // Input booleans
   
   // the services:
   theService = new MuonServiceProxy(parameters.getParameter<ParameterSet>("ServiceParameters"));
   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 62 of file MuonTiming.cc.

                         {
   delete theService;
 }

Member Function Documentation

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

override

Inizialize parameters for histo binning.

Definition at line 130 of file MuonTiming.cc.

References Reference_intrackfit_cff::barrel, makeMuonMisalignmentScenario::endcap, HcalObjRepresent::Fill(), edm::Event::getByToken(), glb, edm::HandleBase::isValid(), LogTrace, metsig::muon, extraflags_cff::muons, reco::MuonTime::nDof, muon::overlap(), protons_cff::time, reco::MuonTime::timeAtIpInOut, reco::MuonTime::timeAtIpInOutErr, HiIsolationCommonParameters_cff::track, and TriggerAnalyzer::veto.

                                                                            {
 
   LogTrace(metname_)<<"[MuonTiming] Analyze the mu";
   theService->update(iSetup);
  
   // 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 &
	)

override

Definition at line 66 of file MuonTiming.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::cd(), MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), and DQMStore::IBooker::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 66 of file MuonTiming.h.

double MuonTiming::etaBarrelMin_

private

Definition at line 66 of file MuonTiming.h.

double MuonTiming::etaEndcapMax_

private

Definition at line 66 of file MuonTiming.h.

double MuonTiming::etaEndcapMin_

private

Definition at line 66 of file MuonTiming.h.

double MuonTiming::etamax_

private

Definition at line 65 of file MuonTiming.h.

double MuonTiming::etamin_

private

Definition at line 65 of file MuonTiming.h.

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

private

Definition at line 69 of file MuonTiming.h.

int MuonTiming::etanbins_

private

Definition at line 58 of file MuonTiming.h.

double MuonTiming::etaOverlapMax_

private

Definition at line 66 of file MuonTiming.h.

double MuonTiming::etaOverlapMin_

private

Definition at line 66 of file MuonTiming.h.

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

private

Definition at line 89 of file MuonTiming.h.

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

private

Definition at line 90 of file MuonTiming.h.

std::string MuonTiming::metname_

private

Definition at line 48 of file MuonTiming.h.

double MuonTiming::ndofmax_

private

Definition at line 63 of file MuonTiming.h.

double MuonTiming::ndofmin_

private

Definition at line 63 of file MuonTiming.h.

int MuonTiming::ndofnbins_

private

Definition at line 56 of file MuonTiming.h.

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

private

Definition at line 69 of file MuonTiming.h.

double MuonTiming::ptmax_

private

Definition at line 64 of file MuonTiming.h.

double MuonTiming::ptmin_

private

Definition at line 64 of file MuonTiming.h.

int MuonTiming::ptnbins_

private

Definition at line 57 of file MuonTiming.h.

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

private

Definition at line 91 of file MuonTiming.h.

double MuonTiming::terrmax_

private

Definition at line 61 of file MuonTiming.h.

double MuonTiming::terrmaxrpc_

private

Definition at line 62 of file MuonTiming.h.

double MuonTiming::terrmin_

private

Definition at line 61 of file MuonTiming.h.

double MuonTiming::terrminrpc_

private

Definition at line 62 of file MuonTiming.h.

int MuonTiming::terrnbins_

private

Definition at line 54 of file MuonTiming.h.

int MuonTiming::terrnbinsrpc_

private

Definition at line 55 of file MuonTiming.h.

std::string MuonTiming::theFolder_

private

Definition at line 68 of file MuonTiming.h.

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

private

Definition at line 46 of file MuonTiming.h.

MuonServiceProxy* MuonTiming::theService

private

Definition at line 44 of file MuonTiming.h.

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

private

Definition at line 85 of file MuonTiming.h.

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

private

Definition at line 87 of file MuonTiming.h.

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

private

Definition at line 88 of file MuonTiming.h.

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

private

Definition at line 86 of file MuonTiming.h.

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

private

Definition at line 84 of file MuonTiming.h.

double MuonTiming::tmax_

private

Definition at line 59 of file MuonTiming.h.

double MuonTiming::tmaxrpc_

private

Definition at line 60 of file MuonTiming.h.

double MuonTiming::tmin_

private

Definition at line 59 of file MuonTiming.h.

double MuonTiming::tminrpc_

private

Definition at line 60 of file MuonTiming.h.

int MuonTiming::tnbins_

private

Definition at line 52 of file MuonTiming.h.

int MuonTiming::tnbinsrpc_

private

Definition at line 53 of file MuonTiming.h.

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

private

Definition at line 92 of file MuonTiming.h.

Public Member Functions

Private Types

Private Attributes

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation