- Timestamp:
- Jun 17, 2009, 10:01:34 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/SmearUtil.cc
r423 r430 1322 1322 1323 1323 //********** returns a segmented value for eta and phi, for calo towers ***** 1324 /*1325 void RESOLution::BinEtaPhi(const float phi, const float eta, float& iPhi, float& iEta){1326 iEta = UNDEFINED;1327 int index= iUNDEFINED;1328 for (unsigned int i=1; i< TOWER_number+1; i++) {1329 if(fabs(eta)>=TOWER_eta_edges[i-1] && fabs(eta)<TOWER_eta_edges[i]) {1330 //if(fabs(eta)>TOWER_eta_edges[i-1] && fabs(eta)<TOWER_eta_edges[i]) {1331 float deta = TOWER_eta_edges[i] - TOWER_eta_edges[i-1];1332 // iEta = (eta>0) ? TOWER_eta_edges[i-1] : -TOWER_eta_edges[i];1333 iEta + deta/2.; // update Delphes > 1.71334 index = i-1;1335 break;1336 }1337 }1338 if(index==UNDEFINED) return;1339 iPhi = UNDEFINED;1340 float dphi = TOWER_dphi[index]*pi/180.;1341 for (unsigned int i=1; i < 360/TOWER_dphi[index]; i++ ) {1342 float low = -pi+(i-1)*dphi;1343 float high= low+dphi;1344 if(phi >= low && phi < high ){1345 //if(phi > low && phi < high ){1346 //iPhi = low;1347 iPhi = low+dphi/2.; // update Delphes > 1.71348 break;1349 }1350 }1351 if (phi > pi-dphi) iPhi = pi-dphi;1352 }1353 */1354 1324 void RESOLution::BinEtaPhi(const float phi, const float eta, float& iPhi, float& iEta){ 1355 1325 iEta = UNDEFINED;
Note:
See TracChangeset
for help on using the changeset viewer.