Skip to content

Commit

Permalink
add MuonMonitor
Browse files Browse the repository at this point in the history
  • Loading branch information
Alicia Calderon committed Oct 23, 2019
1 parent 1e43325 commit 812b438
Show file tree
Hide file tree
Showing 5 changed files with 531 additions and 0 deletions.
21 changes: 21 additions & 0 deletions DQM/MuonMonitor/BuildFile.xml
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,21 @@
<use name="Geometry/CSCGeometry"/>
<use name="FWCore/Framework"/>
<use name="DQMServices/Core"/>
<use name="DQMServices/Components"/>
<use name="FWCore/PluginManager"/>
<use name="DataFormats/MuonReco"/>
<use name="DataFormats/L1GlobalMuonTrigger"/>
<use name="DataFormats/FEDRawData"/>
<use name="RecoMuon/TrackingTools"/>
<use name="TrackingTools/TransientTrack"/>
<use name="DataFormats/CSCDigi"/>
<use name="DataFormats/CSCRecHit"/>
<use name="DataFormats/DTRecHit"/>
<use name="DataFormats/RPCDigi"/>
<use name="DataFormats/PatCandidates"/>
<use name="CondFormats/DTObjects"/>
<use name="CondFormats/CSCObjects"/>
<use name="CondFormats/DataRecord"/>
<use name="EventFilter/CSCRawToDigi"/>
<use name="CommonTools/TriggerUtils"/>
<flags EDM_PLUGIN="1"/>
111 changes: 111 additions & 0 deletions DQM/MuonMonitor/interface/CosmicMuonRecoAnalyzer.h
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,111 @@
#ifndef CosmicMuonRecoAnalyzer_H
#define CosmicMuonRecoAnalyzer_H

/** \class CosmicMuRecoAnalyzer
*
* DQM monitoring source for muon reco track
*
* \author A. Calderon IFCA-CSIC-Unv.Cantabria
*/

#include <memory>
#include <fstream>
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ServiceRegistry/interface/Service.h"

#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"

#include "RecoMuon/TrackingTools/interface/MuonServiceProxy.h"

#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include "DataFormats/Scalers/interface/DcsStatus.h"

class CosmicMuonRecoAnalyzer : public DQMEDAnalyzer {
public:
/// Constructor
CosmicMuonRecoAnalyzer(const edm::ParameterSet&);

/// Destructor
~CosmicMuonRecoAnalyzer() override;

/// Inizialize parameters for histo binning
void analyze(const edm::Event&, const edm::EventSetup&) override;
void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;

private:
// ----------member data ---------------------------
MuonServiceProxy* theService;
edm::ParameterSet parameters;

edm::EDGetTokenT<edm::View<reco::Track> > theMuonCollectionLabel_;

//histo binning parameters

int hitsBin;
int hitsMin;
int hitsMax;

int nTrkBin;
int nTrkMax;
int nTrkMin;

int etaBin;
double etaMin;
double etaMax;

int thetaBin;
double thetaMin;
double thetaMax;

int phiBin;
double phiMin;
double phiMax;

int chi2Bin;
double chi2Min;
double chi2Max;

int pBin;
double pMin;
double pMax;

int ptBin;
double ptMin;
double ptMax;

int pResBin;
double pResMin;
double pResMax;

// sta muon
MonitorElement* nTracksSta;

MonitorElement* etaStaTrack;
MonitorElement* thetaStaTrack;
MonitorElement* phiStaTrack;
MonitorElement* chi2OvDFStaTrack;
MonitorElement* probchi2StaTrack;
MonitorElement* pStaTrack;
MonitorElement* qOverPStaTrack;
MonitorElement* ptStaTrack;
MonitorElement* qStaTrack;
MonitorElement* nValidHitsStaTrack;

MonitorElement* qOverPStaTrack_p;
MonitorElement* phiVsetaStaTrack;
MonitorElement* nValidHitsStaTrack_eta;
MonitorElement* nValidHitsStaTrack_phi;

std::string theFolder;
};
#endif
26 changes: 26 additions & 0 deletions DQM/MuonMonitor/python/muonCosmicAnalyzer_cff.py
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,26 @@
import FWCore.ParameterSet.Config as cms


from RecoMuon.TrackingTools.MuonServiceProxy_cff import *
from DQMServices.Core.DQMEDAnalyzer import DQMEDAnalyzer
dqmInfoMuons = DQMEDAnalyzer('DQMEventInfo',
subSystemFolder = cms.untracked.string('Muons')
)


from DQMOffline.Muon.segmentTrackAnalyzer_cfi import *
from DQMOffline.Muon.muonSeedsAnalyzer_cfi import *
from DQM.MuonMonitor.muonRecoAnalyzer_cfi import *

#muonCosmicGlbSegmentAnalyzer = glbMuonSegmentAnalyzer.clone()
muonCosmicStaSegmentAnalyzer = staMuonSegmentAnalyzer.clone()
muonCosmicSeedAnalyzer = muonSeedsAnalyzer.clone()


#muonCosmicGlbSegmentAnalyzer.MuTrackCollection = cms.InputTag('globalCosmicMuons')
muonCosmicSeedAnalyzer.SeedCollection = cms.InputTag('CosmicMuonSeed')
muonCosmicStaSegmentAnalyzer.MuTrackCollection = cms.InputTag('cosmicMuons')



muonCosmicAnalyzer = cms.Sequence(muonCosmicStaSegmentAnalyzer * muonCosmicSeedAnalyzer * muonRecoAnalyzer)
178 changes: 178 additions & 0 deletions DQM/MuonMonitor/src/CosmicMuonRecoAnalyzer.cc
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,178 @@
#include "DQM/MuonMonitor/interface/CosmicMuonRecoAnalyzer.h"

#include "FWCore/PluginManager/interface/ModuleDef.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/MuonReco/interface/MuonEnergy.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackBase.h"

#include "DataFormats/TrackReco/interface/TrackFwd.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <string>
#include "TMath.h"
#include "TList.h"

using namespace std;
using namespace edm;

CosmicMuonRecoAnalyzer::CosmicMuonRecoAnalyzer(const edm::ParameterSet& pSet) {
parameters = pSet;

// the services:
theService = new MuonServiceProxy(parameters.getParameter<ParameterSet>("ServiceParameters"));
theMuonCollectionLabel_ = consumes<edm::View<reco::Track> >(parameters.getParameter<edm::InputTag>("MuonCollection"));

nTrkBin = parameters.getParameter<int>("nTrkBin");
nTrkMax = parameters.getParameter<int>("nTrkMax");
nTrkMin = parameters.getParameter<int>("nTrkMin");
ptBin = parameters.getParameter<int>("ptBin");
ptMin = parameters.getParameter<double>("ptMin");
ptMax = parameters.getParameter<double>("ptMax");
pBin = parameters.getParameter<int>("pBin");
pMin = parameters.getParameter<double>("pMin");
pMax = parameters.getParameter<double>("pMax");
chi2Bin = parameters.getParameter<int>("chi2Bin");
chi2Min = parameters.getParameter<double>("chi2Min");
chi2Max = parameters.getParameter<double>("chi2Max");
phiBin = parameters.getParameter<int>("phiBin");
phiMin = parameters.getParameter<double>("phiMin");
phiMax = parameters.getParameter<double>("phiMax");
thetaBin = parameters.getParameter<int>("thetaBin");
thetaMin = parameters.getParameter<double>("thetaMin");
thetaMax = parameters.getParameter<double>("thetaMax");
etaBin = parameters.getParameter<int>("etaBin");
etaMin = parameters.getParameter<double>("etaMin");
etaMax = parameters.getParameter<double>("etaMax");
pResBin = parameters.getParameter<int>("pResBin");
pResMin = parameters.getParameter<double>("pResMin");
pResMax = parameters.getParameter<double>("pResMax");

hitsBin = parameters.getParameter<int>("hitsBin");
hitsMin = parameters.getParameter<int>("hitsMin");
hitsMax = parameters.getParameter<int>("hitsMax");

theFolder = parameters.getParameter<string>("folder");
}

CosmicMuonRecoAnalyzer::~CosmicMuonRecoAnalyzer() { delete theService; }

void CosmicMuonRecoAnalyzer::bookHistograms(DQMStore::IBooker& ibooker,
edm::Run const& /*iRun*/,
edm::EventSetup const& /* iSetup */) {
ibooker.cd();
ibooker.setCurrentFolder(theFolder);

/////////////////////////////////////////////////////
// monitoring of eta parameter
/////////////////////////////////////////////////////
std::string histname = "cosmicMuon";

nTracksSta = ibooker.book1D(histname + "_traks", "#tracks", nTrkBin, nTrkMin, nTrkMax);

etaStaTrack = ibooker.book1D(histname + "_eta", "#eta_{STA}", etaBin, etaMin, etaMax);

//////////////////////////////////////////////////////
// monitoring of theta parameter
/////////////////////////////////////////////////////

thetaStaTrack = ibooker.book1D(histname + "_theta", "#theta_{STA}", thetaBin, thetaMin, thetaMax);
thetaStaTrack->setAxisTitle("rad");

// monitoring of phi paramater

phiStaTrack = ibooker.book1D(histname + "_phi", "#phi_{STA}", phiBin, phiMin, phiMax);
phiStaTrack->setAxisTitle("rad");

// monitoring of the chi2 parameter

chi2OvDFStaTrack = ibooker.book1D(histname + "_chi2OverDf", "#chi_{2}OverDF_{STA}", chi2Bin, chi2Min, chi2Max);

//--------------------------

probchi2StaTrack = ibooker.book1D(histname + "_probchi", "Prob #chi_{STA}", 120, chi2Min, 1.20);

// monitoring of the momentum

pStaTrack = ibooker.book1D(histname + "_p", "p_{STA}", pBin, pMin, pMax);
pStaTrack->setAxisTitle("GeV");

qOverPStaTrack = ibooker.book1D(histname + "_qoverp", "qoverp_{STA}", pResBin, pResMin, pResMax);
qOverPStaTrack->setAxisTitle("1/GeV");

qOverPStaTrack_p =
ibooker.book2D(histname + "_qoverp_p", "qoverp_p_{STA}", pBin, pMin, pMax, pResBin, pResMin, pResMax);

// monitoring of the transverse momentum

ptStaTrack = ibooker.book1D(histname + "_pt", "pt_{STA}", ptBin, ptMin, pMax);
ptStaTrack->setAxisTitle("GeV");

// monitoring of the muon charge

qStaTrack = ibooker.book1D(histname + "_q", "q_{STA}", 5, -2.5, 2.5);

//////////////////////////////////////////////////////////////
// monitoring of the phi-eta

phiVsetaStaTrack = ibooker.book2D(
histname + "_phiVseta", "#phi vs #eta (STA)", etaBin / 2, etaMin, etaMax, phiBin / 2, phiMin, phiMax);
phiVsetaStaTrack->setAxisTitle("eta", 1);
phiVsetaStaTrack->setAxisTitle("phi", 2);

// monitoring the hits
nValidHitsStaTrack = ibooker.book1D(histname + "_nValidHits", "#valid hits (STA)", hitsBin, hitsMin, hitsMax);

nValidHitsStaTrack_eta = ibooker.book2D(
histname + "_nValidHits_eta", "#valid hits vs eta (STA)", etaBin, etaMin, etaMax, hitsBin, hitsMin, hitsMax);

nValidHitsStaTrack_phi = ibooker.book2D(
histname + "_nValidHits_phi", "#valid hits vs phi (STA)", phiBin, phiMin, phiMax, hitsBin, hitsMin, hitsMax);
}

void CosmicMuonRecoAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
LogTrace(metname) << "[MuonRecoAnalyzer] Analyze the mu";
theService->update(iSetup);

// Take the muon container
edm::Handle<edm::View<reco::Track> > muons;
iEvent.getByToken(theMuonCollectionLabel_, muons);

int nTracks_ = 0;

if (!muons.isValid())
return;

for (edm::View<reco::Track>::const_iterator muon = muons->begin(); muon != muons->end(); ++muon) {
nTracks_++;

//Needed for MVA soft muon

// get the track using only the mu spectrometer data

etaStaTrack->Fill(muon->eta());
thetaStaTrack->Fill(muon->theta());
phiStaTrack->Fill(muon->phi());
chi2OvDFStaTrack->Fill(muon->normalizedChi2());
probchi2StaTrack->Fill(TMath::Prob(muon->chi2(), muon->ndof()));
pStaTrack->Fill(muon->p());
ptStaTrack->Fill(muon->pt());
qStaTrack->Fill(muon->charge());
qOverPStaTrack->Fill(muon->qoverp());
qOverPStaTrack_p->Fill(muon->p(), muon->qoverp());
phiVsetaStaTrack->Fill(muon->eta(), muon->phi());

nValidHitsStaTrack->Fill(muon->numberOfValidHits());
nValidHitsStaTrack_eta->Fill(muon->eta(), muon->numberOfValidHits());
nValidHitsStaTrack_phi->Fill(muon->phi(), muon->numberOfValidHits());
}

nTracksSta->Fill(nTracks_);
}

DEFINE_FWK_MODULE(CosmicMuonRecoAnalyzer);
Loading

0 comments on commit 812b438

Please sign in to comment.