|
|
|
#include <MuonRecoAnalyzer.h>
Definition at line 28 of file MuonRecoAnalyzer.h.
| MuonRecoAnalyzer::MuonRecoAnalyzer | ( | const edm::ParameterSet & | pSet, |
| MuonServiceProxy * | theService | ||
| ) |
Constructor.
Definition at line 30 of file MuonRecoAnalyzer.cc.
References parameters.
:MuonAnalyzerBase(theService) { parameters = pSet; }
| MuonRecoAnalyzer::~MuonRecoAnalyzer | ( | ) | [virtual] |
| void MuonRecoAnalyzer::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| const reco::Muon & | recoMu | ||
| ) |
Get the analysis.
Definition at line 326 of file MuonRecoAnalyzer.cc.
References chi2OvDFGlbTrack, chi2OvDFStaTrack, chi2OvDFTrack, reco::Muon::combinedMuon(), etaEfficiency, etaGlbTrack, etaPull, etaResolution, etaStaTrack, etaTrack, MonitorElement::Fill(), HcalObjRepresent::Fill(), GetRes(), reco::Muon::isCaloMuon(), reco::Muon::isGlobalMuon(), reco::Muon::isStandAloneMuon(), reco::Muon::isTrackerMuon(), LogTrace, metname, DetId::Muon, muReco, muVStkSytemRotation, oneOverpPull, oneOverpResolution, oneOverptPull, oneOverptResolution, pGlbTrack, phiEfficiency, phiGlbTrack, phiPull, phiResolution, phiStaTrack, phiTrack, probchi2GlbTrack, probchi2StaTrack, probchi2Track, pStaTrack, ptGlbTrack, pTrack, ptStaTrack, ptTrack, qGlbTrack, qOverpPull, qOverpResolution, qOverptPull, qOverptResolution, qStaTrack, qTrack, rhAnalysis, reco::Muon::standAloneMuon(), thetaGlbTrack, thetaPull, thetaResolution, thetaStaTrack, thetaTrack, reco::Muon::track(), and DetId::Tracker.
{
LogTrace(metname)<<"[MuonRecoAnalyzer] Analyze the mu";
float res=0, pull=0;
if(recoMu.isGlobalMuon()) {
LogTrace(metname)<<"[MuonRecoAnalyzer] The mu is global - filling the histos";
if(recoMu.isTrackerMuon() && recoMu.isStandAloneMuon())
muReco->Fill(1);
if(!(recoMu.isTrackerMuon()) && recoMu.isStandAloneMuon())
muReco->Fill(2);
if(!recoMu.isStandAloneMuon())
LogTrace(metname)<<"[MuonRecoAnalyzer] ERROR: the mu is global but not standalone!";
// get the track combinig the information from both the Tracker and the Spectrometer
reco::TrackRef recoCombinedGlbTrack = recoMu.combinedMuon();
// get the track using only the tracker data
reco::TrackRef recoTkGlbTrack = recoMu.track();
// get the track using only the mu spectrometer data
reco::TrackRef recoStaGlbTrack = recoMu.standAloneMuon();
etaGlbTrack[0]->Fill(recoCombinedGlbTrack->eta());
etaGlbTrack[1]->Fill(recoTkGlbTrack->eta());
etaGlbTrack[2]->Fill(recoStaGlbTrack->eta());
GetRes(recoTkGlbTrack, recoCombinedGlbTrack, "eta", res, pull);
etaResolution[0]->Fill(res);
GetRes(recoCombinedGlbTrack, recoStaGlbTrack, "eta", res, pull);
etaResolution[1]->Fill(res);
GetRes(recoTkGlbTrack, recoStaGlbTrack, "eta", res, pull);
etaResolution[2]->Fill(res);
etaPull->Fill(pull);
etaResolution[3]->Fill(recoCombinedGlbTrack->eta(), recoTkGlbTrack->eta()-recoCombinedGlbTrack->eta());
etaResolution[4]->Fill(recoCombinedGlbTrack->eta(), -recoStaGlbTrack->eta()+recoCombinedGlbTrack->eta());
etaResolution[5]->Fill(recoCombinedGlbTrack->eta(), recoTkGlbTrack->eta()-recoStaGlbTrack->eta());
thetaGlbTrack[0]->Fill(recoCombinedGlbTrack->theta());
thetaGlbTrack[1]->Fill(recoTkGlbTrack->theta());
thetaGlbTrack[2]->Fill(recoStaGlbTrack->theta());
GetRes(recoTkGlbTrack, recoCombinedGlbTrack, "theta", res, pull);
thetaResolution[0]->Fill(res);
GetRes(recoCombinedGlbTrack, recoStaGlbTrack, "theta", res, pull);
thetaResolution[1]->Fill(res);
GetRes(recoTkGlbTrack, recoStaGlbTrack, "theta", res, pull);
thetaResolution[2]->Fill(res);
thetaPull->Fill(pull);
thetaResolution[3]->Fill(recoCombinedGlbTrack->theta(), recoTkGlbTrack->theta()-recoCombinedGlbTrack->theta());
thetaResolution[4]->Fill(recoCombinedGlbTrack->theta(), -recoStaGlbTrack->theta()+recoCombinedGlbTrack->theta());
thetaResolution[5]->Fill(recoCombinedGlbTrack->theta(), recoTkGlbTrack->theta()-recoStaGlbTrack->theta());
phiGlbTrack[0]->Fill(recoCombinedGlbTrack->phi());
phiGlbTrack[1]->Fill(recoTkGlbTrack->phi());
phiGlbTrack[2]->Fill(recoStaGlbTrack->phi());
GetRes(recoTkGlbTrack, recoCombinedGlbTrack, "phi", res, pull);
phiResolution[0]->Fill(res);
GetRes(recoCombinedGlbTrack, recoStaGlbTrack, "phi", res, pull);
phiResolution[1]->Fill(res);
GetRes(recoTkGlbTrack, recoStaGlbTrack, "phi", res, pull);
phiResolution[2]->Fill(res);
phiPull->Fill(pull);
phiResolution[3]->Fill(recoCombinedGlbTrack->phi(), recoTkGlbTrack->phi()-recoCombinedGlbTrack->phi());
phiResolution[4]->Fill(recoCombinedGlbTrack->phi(), -recoStaGlbTrack->phi()+recoCombinedGlbTrack->phi());
phiResolution[5]->Fill(recoCombinedGlbTrack->phi(), recoTkGlbTrack->phi()-recoStaGlbTrack->phi());
chi2OvDFGlbTrack[0]->Fill(recoCombinedGlbTrack->normalizedChi2());
chi2OvDFGlbTrack[1]->Fill(recoTkGlbTrack->normalizedChi2());
chi2OvDFGlbTrack[2]->Fill(recoStaGlbTrack->normalizedChi2());
//-------------------------
// double probchi = TMath::Prob(recoCombinedGlbTrack->normalizedChi2(),recoCombinedGlbTrack->ndof());
// cout << "rellenando histos."<<endl;
probchi2GlbTrack[0]->Fill(TMath::Prob(recoCombinedGlbTrack->chi2(),recoCombinedGlbTrack->ndof()));
probchi2GlbTrack[1]->Fill(TMath::Prob(recoTkGlbTrack->chi2(),recoTkGlbTrack->ndof()));
probchi2GlbTrack[2]->Fill(TMath::Prob(recoStaGlbTrack->chi2(),recoStaGlbTrack->ndof()));
// cout << "rellenados histos."<<endl;
//-------------------------
pGlbTrack[0]->Fill(recoCombinedGlbTrack->p());
pGlbTrack[1]->Fill(recoTkGlbTrack->p());
pGlbTrack[2]->Fill(recoStaGlbTrack->p());
ptGlbTrack[0]->Fill(recoCombinedGlbTrack->pt());
ptGlbTrack[1]->Fill(recoTkGlbTrack->pt());
ptGlbTrack[2]->Fill(recoStaGlbTrack->pt());
qGlbTrack[0]->Fill(recoCombinedGlbTrack->charge());
qGlbTrack[1]->Fill(recoTkGlbTrack->charge());
qGlbTrack[2]->Fill(recoStaGlbTrack->charge());
if(recoCombinedGlbTrack->charge()==recoStaGlbTrack->charge()) qGlbTrack[3]->Fill(1);
else qGlbTrack[3]->Fill(2);
if(recoCombinedGlbTrack->charge()==recoTkGlbTrack->charge()) qGlbTrack[3]->Fill(3);
else qGlbTrack[3]->Fill(4);
if(recoStaGlbTrack->charge()==recoTkGlbTrack->charge()) qGlbTrack[3]->Fill(5);
else qGlbTrack[3]->Fill(6);
if(recoCombinedGlbTrack->charge()!=recoStaGlbTrack->charge() && recoCombinedGlbTrack->charge()!=recoTkGlbTrack->charge()) qGlbTrack[3]->Fill(7);
if(recoCombinedGlbTrack->charge()==recoStaGlbTrack->charge() && recoCombinedGlbTrack->charge()==recoTkGlbTrack->charge()) qGlbTrack[3]->Fill(8);
GetRes(recoTkGlbTrack, recoCombinedGlbTrack, "qOverp", res, pull);
qOverpResolution[0]->Fill(res);
GetRes(recoCombinedGlbTrack, recoStaGlbTrack, "qOverp", res, pull);
qOverpResolution[1]->Fill(res);
GetRes(recoTkGlbTrack, recoStaGlbTrack, "qOverp", res, pull);
qOverpResolution[2]->Fill(res);
qOverpPull->Fill(pull);
GetRes(recoTkGlbTrack, recoCombinedGlbTrack, "oneOverp", res, pull);
oneOverpResolution[0]->Fill(res);
GetRes(recoCombinedGlbTrack, recoStaGlbTrack, "oneOverp", res, pull);
oneOverpResolution[1]->Fill(res);
GetRes(recoTkGlbTrack, recoStaGlbTrack, "oneOverp", res, pull);
oneOverpResolution[2]->Fill(res);
oneOverpPull->Fill(pull);
GetRes(recoTkGlbTrack, recoCombinedGlbTrack, "qOverpt", res, pull);
qOverptResolution[0]->Fill(res);
GetRes(recoCombinedGlbTrack, recoStaGlbTrack, "qOverpt", res, pull);
qOverptResolution[1]->Fill(res);
GetRes(recoTkGlbTrack, recoStaGlbTrack, "qOverpt", res, pull);
qOverptResolution[2]->Fill(res);
qOverptPull->Fill(pull);
GetRes(recoTkGlbTrack, recoCombinedGlbTrack, "oneOverpt", res, pull);
oneOverptResolution[0]->Fill(res);
GetRes(recoCombinedGlbTrack, recoStaGlbTrack, "oneOverpt", res, pull);
oneOverptResolution[1]->Fill(res);
GetRes(recoTkGlbTrack, recoStaGlbTrack, "oneOverpt", res, pull);
oneOverptResolution[2]->Fill(res);
oneOverptPull->Fill(pull);
oneOverptResolution[3]->Fill(recoCombinedGlbTrack->eta(),(1/recoTkGlbTrack->pt())-(1/recoCombinedGlbTrack->pt()));
oneOverptResolution[4]->Fill(recoCombinedGlbTrack->eta(),-(1/recoStaGlbTrack->pt())+(1/recoCombinedGlbTrack->pt()));
oneOverptResolution[5]->Fill(recoCombinedGlbTrack->eta(),(1/recoTkGlbTrack->pt())-(1/recoStaGlbTrack->pt()));
oneOverptResolution[6]->Fill(recoCombinedGlbTrack->phi(),(1/recoTkGlbTrack->pt())-(1/recoCombinedGlbTrack->pt()));
oneOverptResolution[7]->Fill(recoCombinedGlbTrack->phi(),-(1/recoStaGlbTrack->pt())+(1/recoCombinedGlbTrack->pt()));
oneOverptResolution[8]->Fill(recoCombinedGlbTrack->phi(),(1/recoTkGlbTrack->pt())-(1/recoStaGlbTrack->pt()));
oneOverptResolution[9]->Fill(recoCombinedGlbTrack->pt(),(1/recoTkGlbTrack->pt())-(1/recoCombinedGlbTrack->pt()));
oneOverptResolution[10]->Fill(recoCombinedGlbTrack->pt(),-(1/recoStaGlbTrack->pt())+(1/recoCombinedGlbTrack->pt()));
oneOverptResolution[11]->Fill(recoCombinedGlbTrack->pt(),(1/recoTkGlbTrack->pt())-(1/recoStaGlbTrack->pt()));
// valid hits Glb track
double rhGlb = recoCombinedGlbTrack->found();
// valid hits Glb track from Tracker
double rhGlb_StaProvenance=0;
// valid hits Glb track from Sta system
double rhGlb_TkProvenance=0;
for (trackingRecHit_iterator recHit = recoCombinedGlbTrack->recHitsBegin();
recHit!=recoCombinedGlbTrack->recHitsEnd(); ++recHit){
if((*recHit)->isValid()){
DetId id = (*recHit)->geographicalId();
if (id.det() == DetId::Muon)
rhGlb_StaProvenance++;
if (id.det() == DetId::Tracker)
rhGlb_TkProvenance++;
}
}
// valid hits Sta track associated to Glb track
double rhStaGlb = recoStaGlbTrack->recHitsSize();
// valid hits Traker track associated to Glb track
double rhTkGlb = recoTkGlbTrack->found();
// invalid hits Traker track associated to Glb track
double rhTkGlb_notValid = recoTkGlbTrack->lost();
// fill the histos
rhAnalysis[0]->Fill(rhGlb_StaProvenance/rhGlb);
rhAnalysis[1]->Fill(rhGlb_TkProvenance/rhGlb);
rhAnalysis[2]->Fill(rhGlb_StaProvenance/rhStaGlb);
rhAnalysis[3]->Fill(rhGlb_TkProvenance/rhTkGlb);
rhAnalysis[4]->Fill(rhGlb/(rhStaGlb+rhTkGlb));
rhAnalysis[5]->Fill(rhTkGlb_notValid/rhGlb);
// aligment plots (mu system w.r.t. tracker rotation)
if(recoCombinedGlbTrack->charge()>0)
muVStkSytemRotation[0]->Fill(recoCombinedGlbTrack->pt(),recoTkGlbTrack->pt()/recoCombinedGlbTrack->pt());
else
muVStkSytemRotation[1]->Fill(recoCombinedGlbTrack->pt(),recoTkGlbTrack->pt()/recoCombinedGlbTrack->pt());
}
if(recoMu.isTrackerMuon() && !(recoMu.isGlobalMuon())) {
LogTrace(metname)<<"[MuonRecoAnalyzer] The mu is tracker only - filling the histos";
if(recoMu.isStandAloneMuon())
muReco->Fill(3);
if(!(recoMu.isStandAloneMuon()))
muReco->Fill(4);
// get the track using only the tracker data
reco::TrackRef recoTrack = recoMu.track();
etaTrack->Fill(recoTrack->eta());
thetaTrack->Fill(recoTrack->theta());
phiTrack->Fill(recoTrack->phi());
chi2OvDFTrack->Fill(recoTrack->normalizedChi2());
probchi2Track->Fill(TMath::Prob(recoTrack->chi2(),recoTrack->ndof()));
pTrack->Fill(recoTrack->p());
ptTrack->Fill(recoTrack->pt());
qTrack->Fill(recoTrack->charge());
}
if(recoMu.isStandAloneMuon() && !(recoMu.isGlobalMuon())) {
LogTrace(metname)<<"[MuonRecoAnalyzer] The mu is STA only - filling the histos";
if(!(recoMu.isTrackerMuon()))
muReco->Fill(5);
// get the track using only the mu spectrometer data
reco::TrackRef recoStaTrack = recoMu.standAloneMuon();
etaStaTrack->Fill(recoStaTrack->eta());
thetaStaTrack->Fill(recoStaTrack->theta());
phiStaTrack->Fill(recoStaTrack->phi());
chi2OvDFStaTrack->Fill(recoStaTrack->normalizedChi2());
probchi2StaTrack->Fill(TMath::Prob(recoStaTrack->chi2(),recoStaTrack->ndof()));
pStaTrack->Fill(recoStaTrack->p());
ptStaTrack->Fill(recoStaTrack->pt());
qStaTrack->Fill(recoStaTrack->charge());
}
if(recoMu.isCaloMuon() && !(recoMu.isGlobalMuon()) && !(recoMu.isTrackerMuon()) && !(recoMu.isStandAloneMuon()))
muReco->Fill(6);
//efficiency plots
// get the track using only the mu spectrometer data
reco::TrackRef recoStaGlbTrack = recoMu.standAloneMuon();
if(recoMu.isStandAloneMuon()){
etaEfficiency[0]->Fill(recoStaGlbTrack->eta());
phiEfficiency[0]->Fill(recoStaGlbTrack->phi());
}
if(recoMu.isStandAloneMuon() && recoMu.isGlobalMuon()){
etaEfficiency[1]->Fill(recoStaGlbTrack->eta());
phiEfficiency[1]->Fill(recoStaGlbTrack->phi());
}
}
| void MuonRecoAnalyzer::beginJob | ( | DQMStore * | dbe | ) | [virtual] |
Inizialize parameters for histo binning.
Implements MuonAnalyzerBase.
Definition at line 38 of file MuonRecoAnalyzer.cc.
References DQMStore::book1D(), DQMStore::book2D(), chi2Bin, chi2Max, chi2Min, chi2OvDFGlbTrack, chi2OvDFStaTrack, chi2OvDFTrack, etaBin, etaEfficiency, etaGlbTrack, etaMax, etaMin, etaPull, etaResolution, etaStaTrack, etaTrack, edm::ParameterSet::getParameter(), LogTrace, metname, muReco, muVStkSytemRotation, oneOverpPull, oneOverpResolution, oneOverptPull, oneOverptResolution, parameters, pBin, pGlbTrack, phiBin, phiEfficiency, phiGlbTrack, phiMax, phiMin, phiPull, phiResolution, phiStaTrack, phiTrack, pMax, pMin, pResBin, pResMax, pResMin, probchi2GlbTrack, probchi2StaTrack, probchi2Track, pStaTrack, ptBin, ptGlbTrack, ptMax, ptMin, pTrack, ptStaTrack, ptTrack, qGlbTrack, qOverpPull, qOverpResolution, qOverptPull, qOverptResolution, qStaTrack, qTrack, rhAnalysis, rhBin, rhMax, rhMin, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), thetaBin, thetaGlbTrack, thetaMax, thetaMin, thetaPull, thetaResolution, thetaStaTrack, and thetaTrack.
{
metname = "muRecoAnalyzer";
LogTrace(metname)<<"[MuonRecoAnalyzer] Parameters initialization";
dbe->setCurrentFolder("Muons/MuonRecoAnalyzer");
muReco = dbe->book1D("muReco", "muon reconstructed tracks", 6, 1, 7);
muReco->setBinLabel(1,"glb+tk+sta");
muReco->setBinLabel(2,"glb+sta");
muReco->setBinLabel(3,"tk+sta");
muReco->setBinLabel(4,"tk");
muReco->setBinLabel(5,"sta");
muReco->setBinLabel(6,"calo");
int binFactor = 4;
// monitoring of eta parameter
etaBin = parameters.getParameter<int>("etaBin");
etaMin = parameters.getParameter<double>("etaMin");
etaMax = parameters.getParameter<double>("etaMax");
std::string histname = "GlbMuon_";
etaGlbTrack.push_back(dbe->book1D(histname+"Glb_eta", "#eta_{GLB}", etaBin, etaMin, etaMax));
etaGlbTrack.push_back(dbe->book1D(histname+"Tk_eta", "#eta_{TKfromGLB}", etaBin, etaMin, etaMax));
etaGlbTrack.push_back(dbe->book1D(histname+"Sta_eta", "#eta_{STAfromGLB}", etaBin, etaMin, etaMax));
etaResolution.push_back(dbe->book1D("Res_TkGlb_eta", "#eta_{TKfromGLB} - #eta_{GLB}", etaBin*binFactor, etaMin/3000, etaMax/3000));
etaResolution.push_back(dbe->book1D("Res_GlbSta_eta", "#eta_{GLB} - #eta_{STAfromGLB}", etaBin*binFactor, etaMin/100, etaMax/100));
etaResolution.push_back(dbe->book1D("Res_TkSta_eta", "#eta_{TKfromGLB} - #eta_{STAfromGLB}", etaBin*binFactor, etaMin/100, etaMax/100));
etaResolution.push_back(dbe->book2D("ResVsEta_TkGlb_eta", "(#eta_{TKfromGLB} - #eta_{GLB}) vs #eta_{GLB}", etaBin, etaMin, etaMax, etaBin*binFactor, etaMin/3000, etaMax/3000));
etaResolution.push_back(dbe->book2D("ResVsEta_GlbSta_eta", "(#eta_{GLB} - #eta_{STAfromGLB}) vs #eta_{GLB}", etaBin, etaMin, etaMax, etaBin*binFactor, etaMin/100, etaMax/100));
etaResolution.push_back(dbe->book2D("ResVsEta_TkSta_eta", "(#eta_{TKfromGLB} - #eta_{STAfromGLB}) vs #eta_{TKfromGLB}", etaBin, etaMin, etaMax, etaBin*binFactor, etaMin/100, etaMax/100));
etaPull = dbe->book1D("Pull_TkSta_eta", "#eta_{TKfromGLB} - #eta_{GLB} / error", 100,-10,10);
etaTrack = dbe->book1D("TkMuon_eta", "#eta_{TK}", etaBin, etaMin, etaMax);
etaStaTrack = dbe->book1D("StaMuon_eta", "#eta_{STA}", etaBin, etaMin, etaMax);
etaEfficiency.push_back(dbe->book1D("StaEta", "#eta_{STAfromGLB}", etaBin, etaMin, etaMax));
etaEfficiency.push_back(dbe->book1D("StaEta_ifCombinedAlso", "#eta_{STAfromGLB} if isGlb=true", etaBin, etaMin, etaMax));
// monitoring of theta parameter
thetaBin = parameters.getParameter<int>("thetaBin");
thetaMin = parameters.getParameter<double>("thetaMin");
thetaMax = parameters.getParameter<double>("thetaMax");
thetaGlbTrack.push_back(dbe->book1D(histname+"Glb_theta", "#theta_{GLB}", thetaBin, thetaMin, thetaMax));
thetaGlbTrack[0]->setAxisTitle("rad");
thetaGlbTrack.push_back(dbe->book1D(histname+"Tk_theta", "#theta_{TKfromGLB}", thetaBin, thetaMin, thetaMax));
thetaGlbTrack[1]->setAxisTitle("rad");
thetaGlbTrack.push_back(dbe->book1D(histname+"Sta_theta", "#theta_{STAfromGLB}", thetaBin, thetaMin, thetaMax));
thetaGlbTrack[2]->setAxisTitle("rad");
thetaResolution.push_back(dbe->book1D("Res_TkGlb_theta", "#theta_{TKfromGLB} - #theta_{GLB}", thetaBin*binFactor, -(thetaMax/3000), thetaMax/3000));
thetaResolution[0]->setAxisTitle("rad");
thetaResolution.push_back(dbe->book1D("Res_GlbSta_theta", "#theta_{GLB} - #theta_{STAfromGLB}", thetaBin*binFactor,-(thetaMax/100), thetaMax/100));
thetaResolution[1]->setAxisTitle("rad");
thetaResolution.push_back(dbe->book1D("Res_TkSta_theta", "#theta_{TKfromGLB} - #theta_{STAfromGLB}", thetaBin*binFactor, -(thetaMax/100), thetaMax/100));
thetaResolution[2]->setAxisTitle("rad");
thetaResolution.push_back(dbe->book2D("ResVsTheta_TkGlb_theta", "(#theta_{TKfromGLB} - #theta_{GLB}) vs #theta_{GLB}", thetaBin, thetaMin, thetaMax, thetaBin*binFactor, -(thetaMax/3000), thetaMax/3000));
thetaResolution[3]->setAxisTitle("rad",1);
thetaResolution[3]->setAxisTitle("rad",2);
thetaResolution.push_back(dbe->book2D("ResVsTheta_GlbSta_theta", "(#theta_{GLB} - #theta_{STAfromGLB}) vs #theta_{GLB}", thetaBin, thetaMin, thetaMax, thetaBin*binFactor, -(thetaMax/100), thetaMax/100));
thetaResolution[4]->setAxisTitle("rad",1);
thetaResolution[4]->setAxisTitle("rad",2);
thetaResolution.push_back(dbe->book2D("ResVsTheta_TkSta_theta", "(#theta_{TKfromGLB} - #theta_{STAfromGLB}) vs #theta_{TKfromGLB}", thetaBin, thetaMin, thetaMax, thetaBin*binFactor, -(thetaMax/100), thetaMax/100));
thetaResolution[5]->setAxisTitle("rad",1);
thetaResolution[5]->setAxisTitle("rad",2);
thetaPull = dbe->book1D("Pull_TkSta_theta", "#theta_{TKfromGLB} - #theta_{STAfromGLB} / error", 100,-10,10);
thetaTrack = dbe->book1D("TkMuon_theta", "#theta_{TK}", thetaBin, thetaMin, thetaMax);
thetaTrack->setAxisTitle("rad");
thetaStaTrack = dbe->book1D("StaMuon_theta", "#theta_{STA}", thetaBin, thetaMin, thetaMax);
thetaStaTrack->setAxisTitle("rad");
// monitoring of phi paramater
phiBin = parameters.getParameter<int>("phiBin");
phiMin = parameters.getParameter<double>("phiMin");
phiMax = parameters.getParameter<double>("phiMax");
phiGlbTrack.push_back(dbe->book1D(histname+"Glb_phi", "#phi_{GLB}", phiBin, phiMin, phiMax));
phiGlbTrack[0]->setAxisTitle("rad");
phiGlbTrack.push_back(dbe->book1D(histname+"Tk_phi", "#phi_{TKfromGLB}", phiBin, phiMin, phiMax));
phiGlbTrack[1]->setAxisTitle("rad");
phiGlbTrack.push_back(dbe->book1D(histname+"Sta_phi", "#phi_{STAfromGLB}", phiBin, phiMin, phiMax));
phiGlbTrack[2]->setAxisTitle("rad");
phiResolution.push_back(dbe->book1D("Res_TkGlb_phi", "#phi_{TKfromGLB} - #phi_{GLB}", phiBin*binFactor, phiMin/3000, phiMax/3000));
phiResolution[0]->setAxisTitle("rad");
phiResolution.push_back(dbe->book1D("Res_GlbSta_phi", "#phi_{GLB} - #phi_{STAfromGLB}", phiBin*binFactor, phiMin/100, phiMax/100));
phiResolution[1]->setAxisTitle("rad");
phiResolution.push_back(dbe->book1D("Res_TkSta_phi", "#phi_{TKfromGLB} - #phi_{STAfromGLB}", phiBin*binFactor, phiMin/100, phiMax/100));
phiResolution[2]->setAxisTitle("rad");
phiResolution.push_back(dbe->book2D("ResVsPhi_TkGlb_phi", "(#phi_{TKfromGLB} - #phi_{GLB}) vs #phi_GLB", phiBin, phiMin, phiMax, phiBin*binFactor, phiMin/3000, phiMax/3000));
phiResolution[3]->setAxisTitle("rad",1);
phiResolution[3]->setAxisTitle("rad",2);
phiResolution.push_back(dbe->book2D("ResVsPhi_GlbSta_phi", "(#phi_{GLB} - #phi_{STAfromGLB}) vs #phi_{GLB}", phiBin, phiMin, phiMax, phiBin*binFactor, phiMin/100, phiMax/100));
phiResolution[4]->setAxisTitle("rad",1);
phiResolution[4]->setAxisTitle("rad",2);
phiResolution.push_back(dbe->book2D("ResVsPhi_TkSta_phi", "(#phi_{TKfromGLB} - #phi_{STAfromGLB}) vs #phi_{TKfromGLB}", phiBin, phiMin, phiMax, phiBin*binFactor, phiMin/100, phiMax/100));
phiResolution[5]->setAxisTitle("rad",1);
phiResolution[5]->setAxisTitle("rad",2);
phiPull = dbe->book1D("Pull_TkSta_phi", "#phi_{TKfromGLB} - #phi_{STAfromGLB} / error", 100,-10,10);
phiTrack = dbe->book1D("TkMuon_phi", "#phi_{TK}", phiBin, phiMin, phiMax);
phiTrack->setAxisTitle("rad");
phiStaTrack = dbe->book1D("StaMuon_phi", "#phi_{STA}", phiBin, phiMin, phiMax);
phiStaTrack->setAxisTitle("rad");
phiEfficiency.push_back(dbe->book1D("StaPhi", "#phi_{STAfromGLB}", phiBin, phiMin, phiMax));
phiEfficiency[0]->setAxisTitle("rad");
phiEfficiency.push_back(dbe->book1D("StaPhi_ifCombinedAlso", "#phi_{STAfromGLB} if the isGlb=true", phiBin, phiMin, phiMax));
phiEfficiency[1]->setAxisTitle("rad");
// monitoring of the chi2 parameter
chi2Bin = parameters.getParameter<int>("chi2Bin");
chi2Min = parameters.getParameter<double>("chi2Min");
chi2Max = parameters.getParameter<double>("chi2Max");
chi2OvDFGlbTrack.push_back(dbe->book1D(histname+"Glb_chi2OverDf", "#chi_{2}OverDF_{GLB}", chi2Bin, chi2Min, chi2Max));
chi2OvDFGlbTrack.push_back(dbe->book1D(histname+"Tk_chi2OverDf", "#chi_{2}OverDF_{TKfromGLB}", phiBin, chi2Min, chi2Max));
chi2OvDFGlbTrack.push_back(dbe->book1D(histname+"Sta_chi2OverDf", "#chi_{2}OverDF_{STAfromGLB}", chi2Bin, chi2Min, chi2Max));
chi2OvDFTrack = dbe->book1D("TkMuon_chi2OverDf", "#chi_{2}OverDF_{TK}", chi2Bin, chi2Min, chi2Max);
chi2OvDFStaTrack = dbe->book1D("StaMuon_chi2OverDf", "#chi_{2}OverDF_{STA}", chi2Bin, chi2Min, chi2Max);
//--------------------------
probchi2GlbTrack.push_back(dbe->book1D(histname+"Glb_probchi", "Prob #chi_{GLB}", 120, chi2Min, 1.20));
probchi2GlbTrack.push_back(dbe->book1D(histname+"Tk_probchi", "Prob #chi_{TKfromGLB}", 120, chi2Min, 1.20));
probchi2GlbTrack.push_back(dbe->book1D(histname+"Sta_probchi", "Prob #chi_{STAfromGLB}", 120, chi2Min, 1.20));
probchi2Track=dbe->book1D("TkMuon_probchi", "Prob #chi_{TK}", 120, chi2Min, 1.20);
probchi2StaTrack=dbe->book1D("StaMuon_probchi", "Prob #chi_{STA}", 120, chi2Min, 1.20);
// monitoring of the momentum
pBin = parameters.getParameter<int>("pBin");
pMin = parameters.getParameter<double>("pMin");
pMax = parameters.getParameter<double>("pMax");
pGlbTrack.push_back(dbe->book1D(histname+"Glb_p", "p_{GLB}", pBin, pMin, pMax));
pGlbTrack[0]->setAxisTitle("GeV");
pGlbTrack.push_back(dbe->book1D(histname+"Tk_p", "p_{TKfromGLB}", pBin, pMin, pMax));
pGlbTrack[1]->setAxisTitle("GeV");
pGlbTrack.push_back(dbe->book1D(histname+"Sta_p", "p_{STAfromGLB}", pBin, pMin, pMax));
pGlbTrack[2]->setAxisTitle("GeV");
pTrack = dbe->book1D("TkMuon_p", "p_{TK}", pBin, pMin, pMax);
pTrack->setAxisTitle("GeV");
pStaTrack = dbe->book1D("StaMuon_p", "p_{STA}", pBin, pMin, pMax);
pStaTrack->setAxisTitle("GeV");
// monitoring of the transverse momentum
ptBin = parameters.getParameter<int>("ptBin");
ptMin = parameters.getParameter<double>("ptMin");
ptMax = parameters.getParameter<double>("ptMax");
ptGlbTrack.push_back(dbe->book1D(histname+"Glb_pt", "pt_{GLB}", ptBin, ptMin, ptMax));
ptGlbTrack[0]->setAxisTitle("GeV");
ptGlbTrack.push_back(dbe->book1D(histname+"Tk_pt", "pt_{TKfromGLB}", ptBin, ptMin, ptMax));
ptGlbTrack[1]->setAxisTitle("GeV");
ptGlbTrack.push_back(dbe->book1D(histname+"Sta_pt", "pt_{STAfromGLB}", ptBin, ptMin, ptMax));
ptGlbTrack[2]->setAxisTitle("GeV");
ptTrack = dbe->book1D("TkMuon_pt", "pt_{TK}", ptBin, ptMin, ptMax);
ptTrack->setAxisTitle("GeV");
ptStaTrack = dbe->book1D("StaMuon_pt", "pt_{STA}", ptBin, ptMin, pMax);
ptStaTrack->setAxisTitle("GeV");
// monitoring of the muon charge
qGlbTrack.push_back(dbe->book1D(histname+"Glb_q", "q_{GLB}", 5, -2.5, 2.5));
qGlbTrack.push_back(dbe->book1D(histname+"Tk_q", "q_{TKfromGLB}", 5, -2.5, 2.5));
qGlbTrack.push_back(dbe->book1D(histname+"Sta_q", "q_{STAformGLB}", 5, -2.5, 2.5));
qGlbTrack.push_back(dbe->book1D(histname+"qComparison", "comparison between q_{GLB} and q_{TKfromGLB}, q_{STAfromGLB}", 8, 0.5, 8.5));
qGlbTrack[3]->setBinLabel(1,"qGlb=qSta");
qGlbTrack[3]->setBinLabel(2,"qGlb!=qSta");
qGlbTrack[3]->setBinLabel(3,"qGlb=qTk");
qGlbTrack[3]->setBinLabel(4,"qGlb!=qTk");
qGlbTrack[3]->setBinLabel(5,"qSta=qTk");
qGlbTrack[3]->setBinLabel(6,"qSta!=qTk");
qGlbTrack[3]->setBinLabel(7,"qGlb!=qSta,qGlb!=Tk");
qGlbTrack[3]->setBinLabel(8,"qGlb=qSta,qGlb=Tk");
qTrack = dbe->book1D("TkMuon_q", "q_{TK}", 5, -2.5, 2.5);
qStaTrack = dbe->book1D("StaMuon_q", "q_{STA}", 5, -2.5, 2.5);
// monitoring of the momentum resolution
pResBin = parameters.getParameter<int>("pResBin");
pResMin = parameters.getParameter<double>("pResMin");
pResMax = parameters.getParameter<double>("pResMax");
qOverpResolution.push_back(dbe->book1D("Res_TkGlb_qOverp", "(q/p)_{TKfromGLB} - (q/p)_{GLB}", pResBin*binFactor*2, pResMin/10, pResMax/10));
qOverpResolution[0]->setAxisTitle("GeV^{-1}");
qOverpResolution.push_back(dbe->book1D("Res_GlbSta_qOverp", "(q/p)_{GLB} - (q/p)_{STAfromGLB}", pResBin*binFactor, pResMin, pResMax));
qOverpResolution[1]->setAxisTitle("GeV^{-1}");
qOverpResolution.push_back(dbe->book1D("Res_TkSta_qOverp", "(q/p)_{TKfromGLB} - (q/p)_{STAfromGLB}", pResBin*binFactor, pResMin, pResMax));
qOverpResolution[2]->setAxisTitle("GeV^{-1}");
qOverpPull = dbe->book1D("Pull_TkSta_qOverp", "(q/p)_{TKfromGLB} - (q/p)_{STAfromGLB} / error", 100,-10,10);
oneOverpResolution.push_back(dbe->book1D("Res_TkGlb_oneOverp", "(1/p)_{TKfromGLB} - (1/p)_{GLB}", pResBin*binFactor*2, pResMin/10, pResMax/10));
oneOverpResolution[0]->setAxisTitle("GeV^{-1}");
oneOverpResolution.push_back(dbe->book1D("Res_GlbSta_oneOverp", "(1/p)_{GLB} - (1/p)_{STAfromGLB}", pResBin*binFactor, pResMin, pResMax));
oneOverpResolution[1]->setAxisTitle("GeV^{-1}");
oneOverpResolution.push_back(dbe->book1D("Res_TkSta_oneOverp", "(q/p)_{TKfromGLB} - (q/p)_{STAfromGLB}", pResBin*binFactor, pResMin, pResMax));
oneOverpResolution[2]->setAxisTitle("GeV^{-1}");
oneOverpPull = dbe->book1D("Pull_TkSta_oneOverp", "(1/p)_{TKfromGLB} - (1/p)_{STAfromGLB} / error", 100,-10,10);
qOverptResolution.push_back(dbe->book1D("Res_TkGlb_qOverpt", "(q/p_{t})_{TKfromGLB} - (q/p_{t})_{GLB}", pResBin*binFactor*2, pResMin/10, pResMax/10));
qOverptResolution[0]->setAxisTitle("GeV^{-1}");
qOverptResolution.push_back(dbe->book1D("Res_GlbSta_qOverpt", "(q/p_{t})_{GLB} - (q/p_{t})_{STAfromGLB}", pResBin*binFactor, pResMin, pResMax));
qOverptResolution[1]->setAxisTitle("GeV^{-1}");
qOverptResolution.push_back(dbe->book1D("Res_TkSta_qOverpt", "(q/p_{t})_{TKfromGLB} - (q/p_{t})_{STAfromGLB}", pResBin*binFactor, pResMin, pResMax));
qOverptResolution[2]->setAxisTitle("GeV^{-1}");
qOverptPull = dbe->book1D("Pull_TkSta_qOverpt", "(q/pt)_{TKfromGLB} - (q/pt)_{STAfromGLB} / error", 100,-10,10);
oneOverptResolution.push_back(dbe->book1D("Res_TkGlb_oneOverpt", "(1/p_{t})_{TKfromGLB} - (1/p_{t})_{GLB}", pResBin*binFactor*2, pResMin/10, pResMax/10));
oneOverptResolution[0]->setAxisTitle("GeV^{-1}");
oneOverptResolution.push_back(dbe->book1D("Res_GlbSta_oneOverpt", "(1/p_{t})_{GLB} - (1/p_{t})_{STAfromGLB}", pResBin*binFactor, pResMin, pResMax));
oneOverptResolution[1]->setAxisTitle("GeV^{-1}");
oneOverptResolution.push_back(dbe->book1D("Res_TkSta_oneOverpt", "(1/p_{t})_{TKfromGLB} - (1/p_{t})_{STAfromGLB}", pResBin*binFactor, pResMin, pResMax));
oneOverptResolution[2]->setAxisTitle("GeV^{-1}");
oneOverptResolution.push_back(dbe->book2D("ResVsEta_TkGlb_oneOverpt", "(#eta_{TKfromGLB} - #eta_{GLB}) vs (1/p_{t})_{GLB}", etaBin, etaMin, etaMax, pResBin*binFactor*2, pResMin/10, pResMax/10));
oneOverptResolution[3]->setAxisTitle("GeV^{-1}",2);
oneOverptResolution.push_back(dbe->book2D("ResVsEta_GlbSta_oneOverpt", "(#eta_{GLB} - #eta_{STAfromGLB} vs (1/p_{t})_{GLB}", etaBin, etaMin, etaMax, pResBin*binFactor, pResMin, pResMax));
oneOverptResolution[4]->setAxisTitle("GeV^{-1}",2);
oneOverptResolution.push_back(dbe->book2D("ResVsEta_TkSta_oneOverpt", "(#eta_{TKfromGLB} - #eta_{STAfromGLB}) vs (1/p_{t})_{TKfromGLB}", etaBin, etaMin, etaMax, pResBin*binFactor, pResMin, pResMax));
oneOverptResolution[5]->setAxisTitle("GeV^{-1}",2);
oneOverptResolution.push_back(dbe->book2D("ResVsPhi_TkGlb_oneOverpt", "(#phi_{TKfromGLB} - #phi_{GLB}) vs (1/p_{t})_{GLB}", phiBin, phiMin, phiMax, pResBin*binFactor*2, pResMin/10, pResMax/10));
oneOverptResolution[6]->setAxisTitle("rad",1);
oneOverptResolution[6]->setAxisTitle("GeV^{-1}",2);
oneOverptResolution.push_back(dbe->book2D("ResVsPhi_GlbSta_oneOverpt", "(#phi_{GLB} - #phi_{STAfromGLB} vs (1/p_{t})_{GLB}", phiBin, phiMin, phiMax, pResBin*binFactor, pResMin, pResMax));
oneOverptResolution[7]->setAxisTitle("rad",1);
oneOverptResolution[7]->setAxisTitle("GeV^{-1}",2);
oneOverptResolution.push_back(dbe->book2D("ResVsPhi_TkSta_oneOverpt", "(#phi_{TKfromGLB} - #phi_{STAfromGLB}) vs (1/p_{t})_{TKfromGLB}", phiBin, phiMin, phiMax, pResBin*binFactor, pResMin, pResMax));
oneOverptResolution[8]->setAxisTitle("rad",1);
oneOverptResolution[8]->setAxisTitle("GeV^{-1}",2);
oneOverptResolution.push_back(dbe->book2D("ResVsPt_TkGlb_oneOverpt", "((1/p_{t})_{TKfromGLB} - (1/p_{t})_{GLB}) vs (1/p_{t})_{GLB}", ptBin/5, ptMin, ptMax/100, pResBin*binFactor*2, pResMin/10, pResMax/10));
oneOverptResolution[9]->setAxisTitle("GeV^{-1}",1);
oneOverptResolution[9]->setAxisTitle("GeV^{-1}",2);
oneOverptResolution.push_back(dbe->book2D("ResVsPt_GlbSta_oneOverpt", "((1/p_{t})_{GLB} - (1/p_{t})_{STAfromGLB} vs (1/p_{t})_{GLB}", ptBin/5, ptMin, ptMax/100, pResBin*binFactor, pResMin, pResMax));
oneOverptResolution[10]->setAxisTitle("GeV^{-1}",1);
oneOverptResolution[10]->setAxisTitle("GeV^{-1}",2);
oneOverptResolution.push_back(dbe->book2D("ResVsPt_TkSta_oneOverpt", "((1/p_{t})_{TKfromGLB} - (1/p_{t})_{STAfromGLB}) vs (1/p_{t})_{TKfromGLB}", ptBin/5, ptMin, ptMax/100, pResBin*binFactor, pResMin, pResMax));
oneOverptResolution[11]->setAxisTitle("GeV^{-1}",1);
oneOverptResolution[11]->setAxisTitle("GeV^{-1}",2);
oneOverptPull = dbe->book1D("Pull_TkSta_oneOverpt", "(1/pt)_{TKfromGLB} - (1/pt)_{STAfromGLB} / error", 100,-10,10);
// monitoring of the recHits provenance
rhBin=parameters.getParameter<int>("rhBin");
rhMin=parameters.getParameter<double>("rhMin");
rhMax=parameters.getParameter<double>("rhMax");
rhAnalysis.push_back(dbe->book1D("StaRh_Frac_inGlb", "recHits_{STAinGLB} / recHits_{GLB}", rhBin, rhMin, rhMax));
rhAnalysis.push_back(dbe->book1D("TkRh_Frac_inGlb", "recHits_{TKinGLB} / recHits_{GLB}", rhBin, rhMin, rhMax));
rhAnalysis.push_back(dbe->book1D("StaRh_inGlb_Div_RhAssoSta", "recHits_{STAinGLB} / recHits_{STAfromGLB}", rhBin, rhMin, rhMax));
rhAnalysis.push_back(dbe->book1D("TkRh_inGlb_Div_RhAssoTk", "recHits_{TKinGLB} / recHits_{TKfromGLB}", rhBin, rhMin, rhMax));
rhAnalysis.push_back(dbe->book1D("GlbRh_Div_RhAssoStaTk", "recHits_{GLB} / (recHits_{TKfromGLB}+recHits_{STAfromGLB})", rhBin, rhMin, rhMax));
rhAnalysis.push_back(dbe->book1D("invalidRh_Frac_inTk", "Invalid recHits / rechits_{GLB}", rhBin, rhMin, rhMax));
// monitoring of the muon system rotation w.r.t. tracker
muVStkSytemRotation.push_back(dbe->book2D("muVStkSytemRotation_posMu", "pT_{TK} / pT_{GLB} vs pT_{GLB} for #mu^{+}", 50,0,200,100,0.8,1.2));
muVStkSytemRotation.push_back(dbe->book2D("muVStkSytemRotation_negMu", "pT_{TK} / pT_{GLB} vs pT_{GLB} for #mu^{-}", 50,0,200,100,0.8,1.2));
}
| void MuonRecoAnalyzer::GetRes | ( | reco::TrackRef | t1, |
| reco::TrackRef | t2, | ||
| std::string | par, | ||
| float & | res, | ||
| float & | pull | ||
| ) |
Definition at line 283 of file MuonRecoAnalyzer.cc.
References p1, p2, and mathSSE::sqrt().
Referenced by analyze().
{
float p1=0, p2=0, p1e=1, p2e=1;
if(par == "eta") {
p1 = t1->eta(); p1e = t1->etaError();
p2 = t2->eta(); p2e = t2->etaError();
}
else if(par == "theta") {
p1 = t1->theta(); p1e = t1->thetaError();
p2 = t2->theta(); p2e = t2->thetaError();
}
else if(par == "phi") {
p1 = t1->phi(); p1e = t1->phiError();
p2 = t2->phi(); p2e = t2->phiError();
}
else if(par == "qOverp") {
p1 = t1->charge()/t1->p(); p1e = t1->qoverpError();
p2 = t2->charge()/t2->p(); p2e = t2->qoverpError();
}
else if(par == "oneOverp") {
p1 = 1./t1->p(); p1e = t1->qoverpError();
p2 = 1./t2->p(); p2e = t2->qoverpError();
}
else if(par == "qOverpt") {
p1 = t1->charge()/t1->pt(); p1e = t1->ptError()*p1*p1;
p2 = t2->charge()/t2->pt(); p2e = t2->ptError()*p2*p2;
}
else if(par == "oneOverpt") {
p1 = 1./t1->pt(); p1e = t1->ptError()*p1*p1;
p2 = 1./t2->pt(); p2e = t2->ptError()*p2*p2;
}
res = p1 - p2;
if(p1e!=0 || p2e!=0) pull = res / sqrt(p1e*p1e + p2e*p2e);
else pull = -99;
return;
}
int MuonRecoAnalyzer::chi2Bin [private] |
Definition at line 68 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
double MuonRecoAnalyzer::chi2Max [private] |
Definition at line 70 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
double MuonRecoAnalyzer::chi2Min [private] |
Definition at line 69 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::chi2OvDFGlbTrack [private] |
Definition at line 97 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::chi2OvDFStaTrack [private] |
Definition at line 130 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::chi2OvDFTrack [private] |
Definition at line 121 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
int MuonRecoAnalyzer::etaBin [private] |
Definition at line 56 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::etaEfficiency [private] |
Definition at line 136 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::etaGlbTrack [private] |
Definition at line 91 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
double MuonRecoAnalyzer::etaMax [private] |
Definition at line 58 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
double MuonRecoAnalyzer::etaMin [private] |
Definition at line 57 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
MonitorElement* MuonRecoAnalyzer::etaPull [private] |
Definition at line 109 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::etaResolution [private] |
Definition at line 92 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::etaStaTrack [private] |
Definition at line 127 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::etaTrack [private] |
Definition at line 118 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
std::string MuonRecoAnalyzer::metname [private] |
Definition at line 51 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::muReco [private] |
Definition at line 89 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::muVStkSytemRotation [private] |
Definition at line 107 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::oneOverpPull [private] |
Definition at line 114 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::oneOverpResolution [private] |
Definition at line 104 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::oneOverptPull [private] |
Definition at line 115 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::oneOverptResolution [private] |
Definition at line 105 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
Definition at line 49 of file MuonRecoAnalyzer.h.
Referenced by beginJob(), and MuonRecoAnalyzer().
int MuonRecoAnalyzer::pBin [private] |
Definition at line 72 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::pGlbTrack [private] |
Definition at line 99 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
int MuonRecoAnalyzer::phiBin [private] |
Definition at line 64 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::phiEfficiency [private] |
Definition at line 137 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::phiGlbTrack [private] |
Definition at line 95 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
double MuonRecoAnalyzer::phiMax [private] |
Definition at line 66 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
double MuonRecoAnalyzer::phiMin [private] |
Definition at line 65 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
MonitorElement* MuonRecoAnalyzer::phiPull [private] |
Definition at line 111 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::phiResolution [private] |
Definition at line 96 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::phiStaTrack [private] |
Definition at line 129 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::phiTrack [private] |
Definition at line 120 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
double MuonRecoAnalyzer::pMax [private] |
Definition at line 74 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
double MuonRecoAnalyzer::pMin [private] |
Definition at line 73 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
int MuonRecoAnalyzer::pResBin [private] |
Definition at line 80 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
double MuonRecoAnalyzer::pResMax [private] |
Definition at line 82 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
double MuonRecoAnalyzer::pResMin [private] |
Definition at line 81 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::probchi2GlbTrack [private] |
Definition at line 98 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::probchi2StaTrack [private] |
Definition at line 131 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::probchi2Track [private] |
Definition at line 122 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::pStaTrack [private] |
Definition at line 132 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
int MuonRecoAnalyzer::ptBin [private] |
Definition at line 76 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::ptGlbTrack [private] |
Definition at line 100 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
double MuonRecoAnalyzer::ptMax [private] |
Definition at line 78 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
double MuonRecoAnalyzer::ptMin [private] |
Definition at line 77 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
MonitorElement* MuonRecoAnalyzer::pTrack [private] |
Definition at line 123 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::ptStaTrack [private] |
Definition at line 133 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::ptTrack [private] |
Definition at line 124 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::qGlbTrack [private] |
Definition at line 101 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::qOverpPull [private] |
Definition at line 112 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::qOverpResolution [private] |
Definition at line 102 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::qOverptPull [private] |
Definition at line 113 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::qOverptResolution [private] |
Definition at line 103 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::qStaTrack [private] |
Definition at line 134 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::qTrack [private] |
Definition at line 125 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::rhAnalysis [private] |
Definition at line 106 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
int MuonRecoAnalyzer::rhBin [private] |
Definition at line 84 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
double MuonRecoAnalyzer::rhMax [private] |
Definition at line 86 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
double MuonRecoAnalyzer::rhMin [private] |
Definition at line 85 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
Definition at line 53 of file MuonRecoAnalyzer.h.
int MuonRecoAnalyzer::thetaBin [private] |
Definition at line 60 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::thetaGlbTrack [private] |
Definition at line 93 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
double MuonRecoAnalyzer::thetaMax [private] |
Definition at line 62 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
double MuonRecoAnalyzer::thetaMin [private] |
Definition at line 61 of file MuonRecoAnalyzer.h.
Referenced by beginJob().
MonitorElement* MuonRecoAnalyzer::thetaPull [private] |
Definition at line 110 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
std::vector<MonitorElement*> MuonRecoAnalyzer::thetaResolution [private] |
Definition at line 94 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::thetaStaTrack [private] |
Definition at line 128 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
MonitorElement* MuonRecoAnalyzer::thetaTrack [private] |
Definition at line 119 of file MuonRecoAnalyzer.h.
Referenced by analyze(), and beginJob().
1.7.3