| 1 |
konec |
1.1 |
#include "UserCode/L1RpcTriggerAnalysis/interface/AnaMuonDistribution.h"
|
| 2 |
|
|
|
| 3 |
konec |
1.4 |
#include "FWCore/ParameterSet/interface/ParameterSet.h"
|
| 4 |
konec |
1.1 |
#include "TProfile.h"
|
| 5 |
|
|
#include "TObjArray.h"
|
| 6 |
|
|
#include "TH2D.h"
|
| 7 |
|
|
#include "TH1D.h"
|
| 8 |
|
|
#include "TGraphErrors.h"
|
| 9 |
|
|
#include "TF1.h"
|
| 10 |
|
|
#include "UserCode/L1RpcTriggerAnalysis/interface/MuonObj.h"
|
| 11 |
|
|
#include "UserCode/L1RpcTriggerAnalysis/interface/Utilities.h"
|
| 12 |
|
|
|
| 13 |
konec |
1.4 |
#include <cmath>
|
| 14 |
|
|
|
| 15 |
|
|
namespace {
|
| 16 |
|
|
TH1D *hMuonPt_DIS, *hMuonEta_DIS, *hMuonPhi_DIS;
|
| 17 |
konec |
1.1 |
TH2D *hMuonPtVsEta_Tk, *hMuonPtVsEta_Ma, *hMuonPtVsEta_Gl;
|
| 18 |
konec |
1.4 |
}
|
| 19 |
|
|
|
| 20 |
|
|
AnaMuonDistribution::AnaMuonDistribution(const edm::ParameterSet& cfg)
|
| 21 |
|
|
: ptMin( cfg.getParameter<double>("ptMin") ),
|
| 22 |
konec |
1.5 |
etaMax (cfg.getParameter<double>("etaMax") ),
|
| 23 |
konec |
1.4 |
minNumberOfMatchedStations( cfg.getParameter<unsigned int>("minNumberOfMatchedStations") ),
|
| 24 |
|
|
minNumberRpcHits( cfg.getParameter<uint>("minNumberRpcHits") ),
|
| 25 |
konec |
1.6 |
minNumberDtCscHits( cfg.getParameter<unsigned int>("minNumberDtCscHits") ),
|
| 26 |
konec |
1.7 |
requireAnyMuon(cfg.getParameter<bool>("requireAnyMuon")),
|
| 27 |
|
|
requireUnique(cfg.getParameter<bool>("requireUnique")),
|
| 28 |
|
|
requireOnlyOne(cfg.getParameter<bool>("requireOnlyOne")),
|
| 29 |
|
|
requireGlobal(cfg.getParameter<bool>("requireGlobal")),
|
| 30 |
|
|
requireInner(cfg.getParameter<bool>("requireInner")),
|
| 31 |
|
|
requireOuter(cfg.getParameter<bool>("requireOuter"))
|
| 32 |
konec |
1.4 |
{ }
|
| 33 |
konec |
1.1 |
|
| 34 |
konec |
1.2 |
void AnaMuonDistribution::init(TObjArray& histos)
|
| 35 |
konec |
1.1 |
{
|
| 36 |
konec |
1.4 |
hMuonPt_DIS = new TH1D("hMuonPt_DIS","All global muons Pt;Glb.muon p_{T} [GeV];Muons / bin",L1PtScale::nPtBins,L1PtScale::ptBins); histos.Add(hMuonPt_DIS);
|
| 37 |
|
|
hMuonEta_DIS = new TH1D("hMuonEta_DIS","All global muons Eta;Glb.muon #eta;Muons / bin",96, -2.4, 2.4); histos.Add(hMuonEta_DIS);
|
| 38 |
|
|
hMuonPhi_DIS = new TH1D("hMuonPhi_DIS","All global muons Phi;Glb.muon #phi [rad];Muons / bin",90,-M_PI,M_PI); histos.Add(hMuonPhi_DIS);
|
| 39 |
konec |
1.1 |
|
| 40 |
|
|
hMuonPtVsEta_Tk = new TH2D("hMuonPtVsEta_Tk","hMuonPtVsEta_Tk", L1RpcEtaScale::nEtaBins, L1RpcEtaScale::etaBins, L1PtScale::nPtBins, L1PtScale::ptBins); histos.Add(hMuonPtVsEta_Tk);
|
| 41 |
|
|
hMuonPtVsEta_Ma = new TH2D("hMuonPtVsEta_Ma","hMuonPtVsEta_Ma", L1RpcEtaScale::nEtaBins, L1RpcEtaScale::etaBins, L1PtScale::nPtBins, L1PtScale::ptBins); histos.Add(hMuonPtVsEta_Ma);
|
| 42 |
|
|
hMuonPtVsEta_Gl = new TH2D("hMuonPtVsEta_Gl","hMuonPtVsEta_Gl", L1RpcEtaScale::nEtaBins, L1RpcEtaScale::etaBins, L1PtScale::nPtBins, L1PtScale::ptBins); histos.Add(hMuonPtVsEta_Gl);
|
| 43 |
|
|
|
| 44 |
|
|
}
|
| 45 |
|
|
|
| 46 |
konec |
1.4 |
bool AnaMuonDistribution::filter(const MuonObj *muon)
|
| 47 |
konec |
1.1 |
{
|
| 48 |
konec |
1.9 |
// std::cout << *muon << std::endl;
|
| 49 |
konec |
1.7 |
if (requireAnyMuon && (muon->nAllMuons==0)) return false;
|
| 50 |
|
|
if (requireOnlyOne && (!muon->nAllMuons==1)) return false;
|
| 51 |
|
|
if (requireUnique && !muon->isUnique) return false;
|
| 52 |
|
|
if (requireGlobal && !muon->isGlobal()) return false;
|
| 53 |
|
|
if (requireInner && !muon->isTracker()) return false;
|
| 54 |
|
|
if (requireOuter && !muon->isOuter()) return false;
|
| 55 |
konec |
1.4 |
if (muon->pt() < ptMin) return false;
|
| 56 |
konec |
1.5 |
if (fabs(muon->eta()) > etaMax) return false;
|
| 57 |
konec |
1.4 |
if (muon->nMatchedStations < minNumberOfMatchedStations) return false;
|
| 58 |
|
|
if (muon->nRPCHits < minNumberRpcHits) return false;
|
| 59 |
|
|
if (muon->nDTHits + muon->nCSCHits < minNumberDtCscHits) return false;
|
| 60 |
|
|
|
| 61 |
|
|
hMuonPt_DIS->Fill(muon->pt());
|
| 62 |
|
|
hMuonEta_DIS->Fill(muon->eta());
|
| 63 |
|
|
hMuonPhi_DIS->Fill(muon->phi());
|
| 64 |
|
|
|
| 65 |
konec |
1.1 |
hMuonPtVsEta_Tk->Fill(muon->eta(), muon->pt());
|
| 66 |
|
|
if (muon->isMatched() ) hMuonPtVsEta_Ma->Fill(muon->eta(), muon->pt());
|
| 67 |
|
|
if (muon->isGlobal()) hMuonPtVsEta_Gl->Fill(muon->eta(), muon->pt());
|
| 68 |
konec |
1.4 |
return true;
|
| 69 |
konec |
1.1 |
}
|
