|
|
|
#include <DiMuonHistograms.h>
Class DiMuonHistograms
DQM monitoring for dimuon mass
Author: S.Folgueras, U. Oviedo
Definition at line 26 of file DiMuonHistograms.h.
| DiMuonHistograms::DiMuonHistograms | ( | const edm::ParameterSet & | pset, |
| MuonServiceProxy * | theService | ||
| ) |
Definition at line 37 of file DiMuonHistograms.cc.
References parameters.
:MuonAnalyzerBase(theService){ parameters = pSet; }
| DiMuonHistograms::~DiMuonHistograms | ( | ) | [virtual] |
Definition at line 41 of file DiMuonHistograms.cc.
{ }
| void DiMuonHistograms::analyze | ( | const edm::Event & | event, |
| const edm::EventSetup & | eventSetup | ||
| ) |
Definition at line 86 of file DiMuonHistograms.cc.
References SiPixelRawToDigiRegional_cfi::beamSpot, DeDxDiscriminatorTools::charge(), etaBMax, etaBMin, EtaCutMax, EtaCutMin, etaECMax, etaECMin, edm::Event::getByLabel(), GlbGlbMuon_HM, GlbGlbMuon_LM, HighMassMin, edm::HandleBase::isValid(), LogTrace, LowMassMax, metname, patZpeak::muons, reco::BeamSpot::position(), SoftSoftMuon, StaTrkMuon_HM, StaTrkMuon_LM, theMuonCollectionLabel, TightTightMuon, TrkTrkMuon_HM, and TrkTrkMuon_LM.
{
LogTrace(metname)<<"[DiMuonHistograms] Analyze the mu in different eta regions";
edm::Handle<reco::MuonCollection> muons;
iEvent.getByLabel(theMuonCollectionLabel, muons);
reco::BeamSpot beamSpot;
Handle<reco::BeamSpot> beamSpotHandle;
iEvent.getByLabel("offlineBeamSpot", beamSpotHandle);
beamSpot = *beamSpotHandle;
if(!muons.isValid()) return;
// Loop on muon collection
TLorentzVector Mu1, Mu2;
float charge = 99.;
float InvMass = -99.;
for (reco::MuonCollection::const_iterator recoMu1 = muons->begin(); recoMu1!=muons->end(); ++recoMu1) {
LogTrace(metname)<<"[DiMuonHistograms] loop over 1st muon"<<endl;
// Loop on second muons to fill invariant mass plots
for (reco::MuonCollection::const_iterator recoMu2 = recoMu1; recoMu2!=muons->end(); ++recoMu2){
LogTrace(metname)<<"[DiMuonHistograms] loop over 2nd muon"<<endl;
if (recoMu1==recoMu2) continue;
// Global-Global Muon
if (recoMu1->isGlobalMuon() && recoMu2->isGlobalMuon()) {
LogTrace(metname)<<"[DiMuonHistograms] Glb-Glb pair"<<endl;
reco::TrackRef recoCombinedGlbTrack1 = recoMu1->combinedMuon();
reco::TrackRef recoCombinedGlbTrack2 = recoMu2->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 (iEtaRegion==0) {EtaCutMin= 0.; EtaCutMax=2.4; }
if (iEtaRegion==1) {EtaCutMin= etaBMin; EtaCutMax=etaBMax; }
if (iEtaRegion==2) {EtaCutMin= etaECMin; EtaCutMax=etaECMax; }
if(fabs(recoCombinedGlbTrack1->eta())>EtaCutMin && fabs(recoCombinedGlbTrack1->eta())<EtaCutMax &&
fabs(recoCombinedGlbTrack2->eta())>EtaCutMin && fabs(recoCombinedGlbTrack2->eta())<EtaCutMax){
if (InvMass < LowMassMax) GlbGlbMuon_LM[iEtaRegion]->Fill(InvMass);
if (InvMass > HighMassMin) GlbGlbMuon_HM[iEtaRegion]->Fill(InvMass);
}
}
}
// Also Tight-Tight Muon Selection
if (recoMu1->isGlobalMuon() && recoMu1->isTrackerMuon() && recoMu1->combinedMuon()->normalizedChi2()<10.
&& recoMu1->combinedMuon()->hitPattern().numberOfValidMuonHits()>0 && fabs(recoMu1->combinedMuon()->dxy(beamSpot.position()))<0.2
&& recoMu1->combinedMuon()->hitPattern().numberOfValidPixelHits()>0 && recoMu1->numberOfMatches() > 1
&& recoMu2->isGlobalMuon() && recoMu2->isTrackerMuon() && recoMu2->combinedMuon()->normalizedChi2()<10.
&& recoMu2->combinedMuon()->hitPattern().numberOfValidMuonHits()>0 && fabs(recoMu2->combinedMuon()->dxy(beamSpot.position()))<0.2
&& recoMu2->combinedMuon()->hitPattern().numberOfValidPixelHits()>0 && recoMu2->numberOfMatches() > 1) {
LogTrace(metname)<<"[DiMuonHistograms] Tight-Tight pair"<<endl;
for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
if (iEtaRegion==0) {EtaCutMin= 0.; EtaCutMax=2.4; }
if (iEtaRegion==1) {EtaCutMin= etaBMin; EtaCutMax=etaBMax; }
if (iEtaRegion==2) {EtaCutMin= etaECMin; EtaCutMax=etaECMax; }
if(fabs(recoCombinedGlbTrack1->eta())>EtaCutMin && fabs(recoCombinedGlbTrack1->eta())<EtaCutMax &&
fabs(recoCombinedGlbTrack2->eta())>EtaCutMin && fabs(recoCombinedGlbTrack2->eta())<EtaCutMax){
if (InvMass > 55. && InvMass < 125.) TightTightMuon[iEtaRegion]->Fill(InvMass);
}
}
}
}
// Now check for STA-TRK
if (recoMu2->isStandAloneMuon() && recoMu1->isTrackerMuon()) {
LogTrace(metname)<<"[DiMuonHistograms] STA-Trk pair"<<endl;
reco::TrackRef recoStaTrack = recoMu2->standAloneMuon();
reco::TrackRef recoTrack = recoMu1->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 (iEtaRegion==0) {EtaCutMin= 0.; EtaCutMax=2.4; }
if (iEtaRegion==1) {EtaCutMin= etaBMin; EtaCutMax=etaBMax; }
if (iEtaRegion==2) {EtaCutMin= etaECMin; EtaCutMax=etaECMax; }
if(fabs(recoStaTrack->eta())>EtaCutMin && fabs(recoStaTrack->eta())<EtaCutMax &&
fabs(recoTrack->eta())>EtaCutMin && fabs(recoTrack->eta())<EtaCutMax){
if (InvMass < LowMassMax) StaTrkMuon_LM[iEtaRegion]->Fill(InvMass);
if (InvMass > HighMassMin) StaTrkMuon_HM[iEtaRegion]->Fill(InvMass);
}
}
}
}
if (recoMu1->isStandAloneMuon() && recoMu2->isTrackerMuon()) {
LogTrace(metname)<<"[DiMuonHistograms] STA-Trk pair"<<endl;
reco::TrackRef recoStaTrack = recoMu1->standAloneMuon();
reco::TrackRef recoTrack = recoMu2->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 (iEtaRegion==0) {EtaCutMin= 0.; EtaCutMax=2.4; }
if (iEtaRegion==1) {EtaCutMin= etaBMin; EtaCutMax=etaBMax; }
if (iEtaRegion==2) {EtaCutMin= etaECMin; EtaCutMax=etaECMax; }
if(fabs(recoStaTrack->eta())>EtaCutMin && fabs(recoStaTrack->eta())<EtaCutMax &&
fabs(recoTrack->eta())>EtaCutMin && fabs(recoTrack->eta())<EtaCutMax){
if (InvMass < LowMassMax) StaTrkMuon_LM[iEtaRegion]->Fill(InvMass);
if (InvMass > HighMassMin) StaTrkMuon_HM[iEtaRegion]->Fill(InvMass);
}
}
}
}
// TRK-TRK dimuon
if (recoMu1->isTrackerMuon() && recoMu2->isTrackerMuon()) {
LogTrace(metname)<<"[DiMuonHistograms] Trk-Trk dimuon pair"<<endl;
reco::TrackRef recoTrack2 = recoMu2->track();
reco::TrackRef recoTrack1 = recoMu1->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 (iEtaRegion==0) {EtaCutMin= 0.; EtaCutMax=2.4; }
if (iEtaRegion==1) {EtaCutMin= etaBMin; EtaCutMax=etaBMax; }
if (iEtaRegion==2) {EtaCutMin= etaECMin; EtaCutMax=etaECMax; }
if(fabs(recoTrack1->eta())>EtaCutMin && fabs(recoTrack1->eta())<EtaCutMax &&
fabs(recoTrack2->eta())>EtaCutMin && fabs(recoTrack2->eta())<EtaCutMax){
if (InvMass < LowMassMax) TrkTrkMuon_LM[iEtaRegion]->Fill(InvMass);
if (InvMass > HighMassMin) TrkTrkMuon_HM[iEtaRegion]->Fill(InvMass);
}
}
}
if (recoMu1->isTrackerMuon() && recoMu2->isTrackerMuon() &&
recoMu1->innerTrack()->found() > 11 && recoMu2->innerTrack()->found() &&
recoMu1->innerTrack()->chi2()/recoMu1->innerTrack()->ndof() < 4.0 &&
recoMu2->innerTrack()->chi2()/recoMu2->innerTrack()->ndof() < 4.0 &&
recoMu1->numberOfMatches() > 0 && recoMu2->numberOfMatches() > 0 &&
recoMu1->innerTrack()->hitPattern().pixelLayersWithMeasurement() > 1 &&
recoMu2->innerTrack()->hitPattern().pixelLayersWithMeasurement() > 1 &&
fabs(recoMu1->innerTrack()->dxy()) < 3.0 && fabs(recoMu1->innerTrack()->dxy()) < 3.0 &&
fabs(recoMu1->innerTrack()->dz()) < 15.0 && fabs(recoMu1->innerTrack()->dz()) < 15.0){
if (charge < 0) {
InvMass = (Mu1+Mu2).M();
for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
if (iEtaRegion==0) {EtaCutMin= 0.; EtaCutMax=2.4; }
if (iEtaRegion==1) {EtaCutMin= etaBMin; EtaCutMax=etaBMax; }
if (iEtaRegion==2) {EtaCutMin= etaECMin; EtaCutMax=etaECMax; }
if(fabs(recoTrack1->eta())>EtaCutMin && fabs(recoTrack1->eta())<EtaCutMax &&
fabs(recoTrack2->eta())>EtaCutMin && fabs(recoTrack2->eta())<EtaCutMax){
SoftSoftMuon[iEtaRegion]->Fill(InvMass);
}
}
}
}
}
} //muon2
} //Muon1
}
| void DiMuonHistograms::beginJob | ( | DQMStore * | dbe | ) | [virtual] |
Inizialize parameters for histo binning.
Implements MuonAnalyzerBase.
Definition at line 43 of file DiMuonHistograms.cc.
References DQMStore::book1D(), etaBBin, etaBin, etaBMax, etaBMin, etaEBin, etaECMax, etaECMin, EtaName, edm::ParameterSet::getParameter(), GlbGlbMuon_HM, GlbGlbMuon_LM, HighMassMax, HighMassMin, LogTrace, LowMassMax, LowMassMin, metname, parameters, DQMStore::setCurrentFolder(), SoftSoftMuon, StaTrkMuon_HM, StaTrkMuon_LM, theMuonCollectionLabel, TightTightMuon, TrkTrkMuon_HM, and TrkTrkMuon_LM.
{
metname = "DiMuonhistograms";
LogTrace(metname)<<"[DiMuonHistograms] Parameters initialization";
dbe->setCurrentFolder("Muons/DiMuonHistograms");
theMuonCollectionLabel = parameters.getParameter<edm::InputTag>("MuonCollection");
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");
int nBin = 0;
for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
if (iEtaRegion==0) { EtaName = ""; nBin = etaBin;}
if (iEtaRegion==1) { EtaName = "_Barrel"; nBin = etaBBin;}
if (iEtaRegion==2) { EtaName = "_EndCap"; nBin = etaEBin;}
GlbGlbMuon_LM.push_back(dbe->book1D("GlbGlbMuon_LM"+EtaName,"InvMass_{GLB,GLB}"+EtaName,nBin, LowMassMin, LowMassMax));
TrkTrkMuon_LM.push_back(dbe->book1D("TrkTrkMuon_LM"+EtaName,"InvMass_{TRK,TRK}"+EtaName,nBin, LowMassMin, LowMassMax));
StaTrkMuon_LM.push_back(dbe->book1D("StaTrkMuon_LM"+EtaName,"InvMass_{STA,TRK}"+EtaName,nBin, LowMassMin, LowMassMax));
GlbGlbMuon_HM.push_back(dbe->book1D("GlbGlbMuon_HM"+EtaName,"InvMass_{GLB,GLB}"+EtaName,nBin, HighMassMin, HighMassMax));
TrkTrkMuon_HM.push_back(dbe->book1D("TrkTrkMuon_HM"+EtaName,"InvMass_{TRK,TRK}"+EtaName,nBin, HighMassMin, HighMassMax));
StaTrkMuon_HM.push_back(dbe->book1D("StaTrkMuon_HM"+EtaName,"InvMass_{STA,TRK}"+EtaName,nBin, HighMassMin, HighMassMax));
// arround the Z peak
TightTightMuon.push_back(dbe->book1D("TightTightMuon"+EtaName,"InvMass_{Tight,Tight}"+EtaName,nBin, 55.0, 125.0));
// low-mass resonances
SoftSoftMuon.push_back(dbe->book1D("SoftSoftMuon"+EtaName,"InvMass_{Soft,Soft}"+EtaName,nBin, 5.0, 55.0));
}
}
int DiMuonHistograms::etaBBin [protected] |
Definition at line 51 of file DiMuonHistograms.h.
Referenced by beginJob().
int DiMuonHistograms::etaBin [protected] |
Definition at line 50 of file DiMuonHistograms.h.
Referenced by beginJob().
double DiMuonHistograms::etaBMax [protected] |
Definition at line 61 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
double DiMuonHistograms::etaBMin [protected] |
Definition at line 60 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
double DiMuonHistograms::EtaCutMax [protected] |
Definition at line 59 of file DiMuonHistograms.h.
Referenced by analyze().
double DiMuonHistograms::EtaCutMin [protected] |
Definition at line 58 of file DiMuonHistograms.h.
Referenced by analyze().
int DiMuonHistograms::etaEBin [protected] |
Definition at line 52 of file DiMuonHistograms.h.
Referenced by beginJob().
double DiMuonHistograms::etaECMax [protected] |
Definition at line 63 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
double DiMuonHistograms::etaECMin [protected] |
Definition at line 62 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
std::string DiMuonHistograms::EtaName [protected] |
Definition at line 56 of file DiMuonHistograms.h.
Referenced by beginJob().
int DiMuonHistograms::etaOvlpBin [protected] |
Definition at line 53 of file DiMuonHistograms.h.
std::vector<MonitorElement*> DiMuonHistograms::GlbGlbMuon_HM [protected] |
Definition at line 72 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> DiMuonHistograms::GlbGlbMuon_LM [protected] |
Definition at line 71 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
double DiMuonHistograms::HighMassMax [protected] |
Definition at line 69 of file DiMuonHistograms.h.
Referenced by beginJob().
double DiMuonHistograms::HighMassMin [protected] |
Definition at line 68 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
double DiMuonHistograms::LowMassMax [protected] |
Definition at line 67 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
double DiMuonHistograms::LowMassMin [protected] |
Definition at line 66 of file DiMuonHistograms.h.
Referenced by beginJob().
std::string DiMuonHistograms::metname [protected] |
Definition at line 43 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
edm::ParameterSet DiMuonHistograms::parameters [protected] |
Definition at line 40 of file DiMuonHistograms.h.
Referenced by beginJob(), and DiMuonHistograms().
std::vector<MonitorElement*> DiMuonHistograms::SoftSoftMuon [protected] |
Definition at line 79 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> DiMuonHistograms::StaTrkMuon_HM [protected] |
Definition at line 74 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> DiMuonHistograms::StaTrkMuon_LM [protected] |
Definition at line 73 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
Definition at line 47 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
edm::InputTag DiMuonHistograms::theSTACollectionLabel [protected] |
Definition at line 46 of file DiMuonHistograms.h.
std::vector<MonitorElement*> DiMuonHistograms::TightTightMuon [protected] |
Definition at line 78 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> DiMuonHistograms::TrkTrkMuon_HM [protected] |
Definition at line 76 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> DiMuonHistograms::TrkTrkMuon_LM [protected] |
Definition at line 75 of file DiMuonHistograms.h.
Referenced by analyze(), and beginJob().
1.7.3