1162 |
int nangbin = TrkParams::GetPFA_NbinsAngle(); |
int nangbin = TrkParams::GetPFA_NbinsAngle(); |
1163 |
int istart = (int)(nangbin/2); |
int istart = (int)(nangbin/2); |
1164 |
int idbin = -1; |
int idbin = -1; |
1165 |
|
// cout << angle <<" ??" << endl; |
1166 |
if( angle > pfa_.angR[istart] ){ |
if( angle > pfa_.angR[istart] ){ |
1167 |
idbin = nangbin-1; |
idbin = nangbin-1; |
1168 |
for(int id=istart+1; id<nangbin; id++ ){ |
for(int id=istart+1; id<nangbin; id++ ){ |
1169 |
|
// cout << id << endl; |
1170 |
if( pfa_.angR[id] >= angle ){ |
if( pfa_.angR[id] >= angle ){ |
1171 |
idbin = id; |
idbin = id; |
1172 |
break; |
break; |
1175 |
}else{ |
}else{ |
1176 |
idbin = 0; |
idbin = 0; |
1177 |
for(int id=istart; id>=0; id-- ){ |
for(int id=istart; id>=0; id-- ){ |
1178 |
|
// cout << id << endl; |
1179 |
if( pfa_.angL[id] < angle ){ |
if( pfa_.angL[id] < angle ){ |
1180 |
idbin = id; |
idbin = id; |
1181 |
break; |
break; |
1183 |
} |
} |
1184 |
|
|
1185 |
} |
} |
1186 |
|
// cout << ">>> "<<idbin<<" "<<pfa_.angL[idbin]<<" "<<pfa_.angR[idbin]<<endl; |
1187 |
return idbin; |
return idbin; |
1188 |
} |
} |
1189 |
|
|