6 |
|
#include "TGraphErrors.h" |
7 |
|
#include "UserCode/L1RpcTriggerAnalysis/interface/RPCDetIdUtil.h" |
8 |
|
#include "UserCode/L1RpcTriggerAnalysis/interface/MuonObj.h" |
9 |
+ |
#include "UserCode/L1RpcTriggerAnalysis/interface/DetCluDigiObj.h" |
10 |
|
#include <sstream> |
11 |
|
#include <bitset> |
12 |
|
#include <algorithm> |
20 |
|
|
21 |
|
void AnaDet::run( |
22 |
|
const MuonObj* muon, |
23 |
< |
const std::vector<uint32_t> & detsHitsCompatibleWithMuon, |
23 |
> |
const std::vector<DetCluDigiObj> & detsHitsCompatibleWithMuon, |
24 |
|
const std::vector<uint32_t> & detsCrossedByMuon, |
25 |
|
const std::vector<uint32_t> & detsCrossedByMuonDeepInside) |
26 |
|
{ |
27 |
|
typedef std::vector<uint32_t>::const_iterator IDET; |
28 |
+ |
typedef std::vector<DetCluDigiObj>::const_iterator IDCD; |
29 |
|
if( muon->pt() < 10.) return; |
30 |
|
// if (fabs(muon->eta()) > 1.14 && fabs(muon->eta()) < 1.24) debug = true; |
31 |
|
|
32 |
|
DetEfficiencyManager aManDIS; |
33 |
< |
for( IDET it= detsHitsCompatibleWithMuon.begin(); it < detsHitsCompatibleWithMuon.end(); it++) aManDIS.addDetHit( *it ); |
33 |
> |
for( IDCD it= detsHitsCompatibleWithMuon.begin(); it < detsHitsCompatibleWithMuon.end(); it++) aManDIS.addDetHit( it->det ); |
34 |
|
for( IDET it= detsCrossedByMuonDeepInside.begin(); it < detsCrossedByMuonDeepInside.end(); it++) aManDIS.addDetMuon( *it ); |
35 |
|
const std::vector<DetEfficiency> & detEfficDIS = aManDIS.dets(); |
36 |
|
|
55 |
|
// |
56 |
|
DetEfficiencyManager aMan; |
57 |
|
if (debug) std::cout <<"detsCrossedByMuon: "<<detsCrossedByMuon.size()<<" detsHitsCompatibleWithMuon: "<<detsHitsCompatibleWithMuon.size()<<std::endl; |
58 |
< |
for( IDET it= detsHitsCompatibleWithMuon.begin(); it < detsHitsCompatibleWithMuon.end(); it++) aMan.addDetHit( *it ); |
58 |
> |
for( IDCD it= detsHitsCompatibleWithMuon.begin(); it < detsHitsCompatibleWithMuon.end(); it++) aMan.addDetHit( it->det ); |
59 |
|
for( IDET it= detsCrossedByMuon.begin(); it < detsCrossedByMuon.end(); it++) aMan.addDetMuon( *it ); |
60 |
|
const std::vector<DetEfficiency> & detEffic= aMan.dets(); |
61 |
|
// |