100 |
|
if ( simu->pt() < 26. || simu->pt() > 27. ) return candidate; |
101 |
|
|
102 |
|
*/ |
103 |
< |
// std::cout <<" ------------------ EVENT: " << std::endl; |
103 |
> |
// std::cout <<" ------------------ EVENT: "<<*ev <<" number of digis: "<<vDigi.size()<< std::endl; |
104 |
|
std::vector<Pattern> vpattern(1); |
105 |
|
theEvForPatCounter++; |
106 |
|
static bool skipRpcData = theConfig.getUntrackedParameter<bool>("skipRpcData", false); |
109 |
|
DetId detId( is->first); |
110 |
|
if (skipRpcData && detId.subdetId()==MuonSubdetId::RPC) continue; |
111 |
|
if (skipDtCscData && (detId.subdetId()==MuonSubdetId::DT || detId.subdetId()==MuonSubdetId::CSC) ) continue; |
112 |
+ |
// std::cout << "adding------- "<< is-vDigi.begin()+1 <<" digi det: "<<is->first<<"(rpc:"<<(detId.subdetId()==MuonSubdetId::RPC)<<") data: "<<is->second<< std::endl; |
113 |
|
Pattern::add(vpattern,*is); |
114 |
+ |
// std::cout <<" after vpattern has size: "<<vpattern.size() << std::endl; |
115 |
+ |
if (vpattern.size() > 100) break; |
116 |
|
} |
117 |
|
if (vpattern[0].size() == 0) return candidate; |
118 |
|
// std::cout <<" ------------------ END EVENT, NOW COMPARE, has #patterns: "<<vpattern.size()<<" vpattern[0].size="<<vpattern[0].size() << std::endl; |
130 |
|
GoldenPattern & gp = igps->second; |
131 |
|
GoldenPattern::Result result = gp.compare(pattern); |
132 |
|
// if (!result.hasRpcDet(637602109) && !result.hasRpcDet(637634877) && !result.hasRpcDet(637599914) && !result.hasRpcDet(637632682)) continue; |
133 |
< |
// if (!result.hasRpcDet(igps->first.theDet)) continue; |
134 |
< |
// if (result.nMatchedTot() < 5 )continue; |
135 |
< |
// if (!result) continue; |
133 |
> |
// if (!result.hasRpcDet(igps->first.theDet)) continue; |
134 |
> |
// if (result.nMatchedTot() < 5 )continue; |
135 |
> |
// if (!result) continue; |
136 |
|
// std::cout <<"PATT KEY: "<<igps->first<<" "<<result ; //<<std::endl; |
137 |
|
if (bestMatching < result) { |
138 |
|
bestMatching = result; |