| 93 |
|
|
| 94 |
|
void L1RpcTreeMaker::analyze(const edm::Event &ev, const edm::EventSetup &es) |
| 95 |
|
{ |
| 96 |
– |
|
| 97 |
– |
|
| 96 |
|
// |
| 97 |
|
// check reference muon |
| 98 |
|
// |
| 162 |
|
// hits and detectors compatible with muon track |
| 163 |
|
// |
| 164 |
|
if ( muon->pt() > 10.) { |
| 167 |
– |
/* |
| 168 |
– |
std::vector<DetHitCompatibleCollector::DetCluDigi> dhcm = theDetHitCollector.compatibleHits( theMuon, ev, es); |
| 169 |
– |
detsHitsCompatibleWithMuon.clear(); |
| 170 |
– |
for (std::vector<DetHitCompatibleCollector::DetCluDigi>::const_iterator it= dhcm.begin(); it < dhcm.end(); ++it) { |
| 171 |
– |
DetCluDigiObj obj; |
| 172 |
– |
obj.det = (*it).first; |
| 173 |
– |
obj.clusterSize = (*it).second.first; |
| 174 |
– |
obj.nDigis = (*it).second.second; |
| 175 |
– |
detsHitsCompatibleWithMuon.push_back(obj); |
| 176 |
– |
} |
| 177 |
– |
*/ |
| 165 |
|
detsHitsCompatibleWithMuon = theDetHitCollector.compatibleHits( theMuon, ev, es); |
| 166 |
|
detsCrossedByMuon = theDetHitCollector.compatibleDets( theMuon, ev, es, false); |
| 167 |
|
detsCrossedByMuonDeepInside = theDetHitCollector.compatibleDets( theMuon, ev, es, true); |
| 168 |
|
if (theConfig.getParameter<bool>("checkDestSIMU")) detsSIMU = theDetHitCollector.compatibleSIMU( theMuon, ev, es); |
| 182 |
– |
// for (uint32_t i=0; i< nDigisCompDets.size(); i++) { |
| 183 |
– |
// if (clSizeCompDets[i] > nDigisCompDets[i]) std::cout <<" PROBLEM, event: "<<theCounter<<" cl:"<<clSizeCompDets[i]<<" nDigis:"<<nDigisCompDets[i]<<std::endl; |
| 184 |
– |
// } |
| 169 |
|
} |
| 170 |
|
|
| 171 |
|
|
| 180 |
|
// |
| 181 |
|
// fill L1 RPCemu |
| 182 |
|
// |
| 183 |
< |
if (theConfig.exists("l1RpcEmu") ) { |
| 183 |
> |
if (theConfig.exists("l1RpcEmu")) { |
| 184 |
|
TrackToL1ObjMatcher matcher(theConfig.getParameter<edm::ParameterSet>("matcherPSet")); |
| 185 |
|
L1ObjMakerRpcEmu l1RpcsFromEmu( theConfig.getParameter<edm::InputTag>("l1RpcEmu"), ev); |
| 186 |
|
std::vector<L1Obj> l1RpcsEmu = l1RpcsFromEmu(); |
