#include "GaudiKernel/SmartDataPtr.h"#include "EvtRecEvent/EvtRecEvent.h"#include "EvtRecEvent/EvtRecTrack.h"Go to the source code of this file.
Functions | |
| double | calcPhotonDang (EvtRecTrack *recTrk, EvtRecEvent *recEvt, EvtRecTrackCol *recTrkCol) |
|
||||||||||||||||
|
find the nearest charged track 00011 {
00012
00013 RecEmcShower *emcTrk = recTrk->emcShower();
00014
00016 Hep3Vector emcpos(emcTrk->x(), emcTrk->y(), emcTrk->z());
00017
00018 double dang = 200.;
00019 for (int j = 0; j < recEvt->totalCharged(); j++)
00020 {
00021 EvtRecTrackIterator jtTrk = recTrkCol->begin() + j;
00022 if ( !(*jtTrk)->isExtTrackValid() ) continue;
00023 RecExtTrack* extTrk = (*jtTrk)->extTrack();
00024 if ( extTrk->emcVolumeNumber() == -1 ) continue;
00025 Hep3Vector extpos = extTrk->emcPosition();
00026 double angd = extpos.angle(emcpos);
00027 if ( angd < dang ) dang = angd;
00028 }
00029
00030 if ( dang < 200 ) dang = dang * 180 / (CLHEP::pi);
00031
00032 return dang;
00033 }
|
1.3.9.1