| 1 |
#include "UserCode/L1RpcTriggerAnalysis/interface/SynchroSelectorL1Muon.h"
|
| 2 |
|
| 3 |
|
| 4 |
#include "FWCore/Framework/interface/Event.h"
|
| 5 |
#include "FWCore/Framework/interface/EventSetup.h"
|
| 6 |
#include "FWCore/Framework/interface/ESHandle.h"
|
| 7 |
#include "FWCore/Utilities/interface/InputTag.h"
|
| 8 |
|
| 9 |
#include "DataFormats/L1GlobalMuonTrigger/interface/L1MuRegionalCand.h"
|
| 10 |
#include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTCand.h"
|
| 11 |
#include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTExtendedCand.h"
|
| 12 |
#include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTReadoutCollection.h"
|
| 13 |
|
| 14 |
#include "Geometry/CommonDetUnit/interface/GeomDet.h"
|
| 15 |
#include "Geometry/CommonDetUnit/interface/GlobalTrackingGeometry.h"
|
| 16 |
#include "Geometry/Records/interface/GlobalTrackingGeometryRecord.h"
|
| 17 |
|
| 18 |
#include "DataFormats/MuonDetId/interface/RPCDetId.h"
|
| 19 |
#include "DataFormats/Math/interface/deltaPhi.h"
|
| 20 |
|
| 21 |
#include "TObjArray.h"
|
| 22 |
#include "TH1F.h"
|
| 23 |
#include <cmath>
|
| 24 |
using namespace edm;
|
| 25 |
using namespace std;
|
| 26 |
|
| 27 |
|
| 28 |
SynchroSelectorL1Muon::SynchroSelectorL1Muon(const edm::ParameterSet& cfg, TObjArray& histos)
|
| 29 |
: SynchroSelector(cfg)
|
| 30 |
{
|
| 31 |
hDeltaEta = new TH1F("hDeltaEta","hDeltaEta",100.,-1.,1.); histos.Add(hDeltaEta);
|
| 32 |
hDeltaPhi = new TH1F("hDeltaPhi","hDeltaPhi",100.,-1.,1.); histos.Add(hDeltaPhi);
|
| 33 |
}
|
| 34 |
|
| 35 |
void SynchroSelectorL1Muon::update(const edm::Event&ev, const edm::EventSetup& es)
|
| 36 |
{
|
| 37 |
theL1Muons.clear();
|
| 38 |
edm::Handle<L1MuGMTReadoutCollection> pCollection;
|
| 39 |
ev.getByLabel(theConfig.getParameter<edm::InputTag>("l1MuReadout"),pCollection);
|
| 40 |
|
| 41 |
L1MuGMTReadoutCollection const* gmtrc = pCollection.product();
|
| 42 |
if (!gmtrc) return;
|
| 43 |
|
| 44 |
vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
|
| 45 |
|
| 46 |
vector<string> systems = theConfig.getParameter<vector<string> >("systems");
|
| 47 |
|
| 48 |
for (vector<string>::const_iterator ic = systems.begin(); ic != systems.end(); ++ic) {
|
| 49 |
for( vector<L1MuGMTReadoutRecord>::const_iterator it = gmt_records.begin() ; it != gmt_records.end() ; it++ ) {
|
| 50 |
vector<L1MuRegionalCand> cands;
|
| 51 |
if ( (*ic) == "rpcb" ) cands = it->getBrlRPCCands();
|
| 52 |
else if ( (*ic) == "rpcf" ) cands = it->getFwdRPCCands();
|
| 53 |
else if ( (*ic) == "csc" ) cands = it->getCSCCands();
|
| 54 |
else if ( (*ic) == "dt" ) cands = it->getDTBXCands();
|
| 55 |
for (vector<L1MuRegionalCand>::const_iterator il=cands.begin(); il < cands.end(); ++il) {
|
| 56 |
if (!il->empty() && ( (il->bx() == 0) || !theConfig.getParameter<bool>("checkBX0") ) )theL1Muons.push_back(*il);
|
| 57 |
}
|
| 58 |
}
|
| 59 |
}
|
| 60 |
}
|
| 61 |
|
| 62 |
|
| 63 |
bool SynchroSelectorL1Muon::takeIt(const RPCDetId & det, const edm::Event&ev, const edm::EventSetup& es)
|
| 64 |
{
|
| 65 |
bool result = false;
|
| 66 |
edm::ESHandle<GlobalTrackingGeometry> globalGeometry;
|
| 67 |
es.get<GlobalTrackingGeometryRecord>().get(globalGeometry);
|
| 68 |
|
| 69 |
|
| 70 |
GlobalPoint detPos = globalGeometry->idToDet(det)->position();
|
| 71 |
float detEta = detPos.eta();
|
| 72 |
float detPhi = detPos.phi();
|
| 73 |
for(vector<L1MuRegionalCand>::const_iterator it =theL1Muons.begin(); it != theL1Muons.end(); ++it) {
|
| 74 |
if (it->empty()) continue;
|
| 75 |
float l1Eta = it->etaValue();
|
| 76 |
float l1Phi = it->phiValue();
|
| 77 |
bool matchedDeltaPhi = fabs(deltaPhi(l1Phi,detPhi)) < theConfig.getParameter<double>("maxDeltaPhi");
|
| 78 |
bool matchedDeltaEta = fabs(l1Eta-detEta) < theConfig.getParameter<double>("maxDeltaEta");
|
| 79 |
if (matchedDeltaPhi) hDeltaEta->Fill(l1Eta-detEta);
|
| 80 |
if (matchedDeltaEta) hDeltaPhi->Fill(deltaPhi(l1Phi,detPhi));
|
| 81 |
if ( !matchedDeltaPhi ) continue;
|
| 82 |
if ( !matchedDeltaEta ) continue;
|
| 83 |
|
| 84 |
// std::cout <<" compatible det, (eta,phi,r,z)= ("<<detEta<<", "<<detPhi<<", "<<detPos.perp()<<", "<<detPos.z()
|
| 85 |
// <<") L1 (eta,phi): ("<< l1Eta<<", "<<l1Phi<<")"<<std::endl;
|
| 86 |
|
| 87 |
result = true;
|
| 88 |
}
|
| 89 |
|
| 90 |
return result;
|
| 91 |
}
|
| 92 |
|
