51 {
52
53
54 MsgStream log(
msgSvc(), name());
56 if(nshower<2) return StatusCode::SUCCESS;
57
58
59
60
61 SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
62 if (!eventHeader) {
63 log << MSG::FATAL << "Could not find Event Header" << endreq;
64 return( StatusCode::FAILURE);
65 }
66
67 SmartDataPtr<MdcDigiCol> mdcDigiCol(eventSvc(),"/Event/Digi/MdcDigiCol");
68 if (!mdcDigiCol) {
69 log << MSG::FATAL << "EmcRec could not find Mdc digi!!" << endreq;
70 return( StatusCode::FAILURE);
71 }
72
73 unsigned int nsec=0;
75 MdcDigiCol::iterator
iter=mdcDigiCol->begin();
76
82
85
88 float phi12=(phi11+phi22-3.1415926)*0.5;
89 float phi21=(phi11+phi22+3.1415926)*0.5;
90 if(phi12<0.) phi12 += 6.2831852;
91 if(phi21>6.2831852) phi21 -= 6.2831852;
92 for(;
iter!= mdcDigiCol->end();
iter++) {
93 id=(*iter)->identify();
96 if(ilayer>=43)
97 log << MSG::ERROR << "MDC(" << ilayer <<","<<iphi <<endreq;
98 float phi=6.2831853*iphi/idmax[ilayer];
99 if(WhetherSector(phi,phi11,phi12)||WhetherSector(phi,phi21,phi22)) nsec++;
100 }
101 log << MSG::INFO << " Number of hits in sector of MDC is " << nsec <<endreq;
102
103
107 if(
m_ef->
getEFVec(37)!=0) log << MSG::WARNING <<
"nsec(0) !=0"<< endreq;
109
111
112 return StatusCode::SUCCESS;
113}
void setValue(float value)
bool getValue(float &value) const
bool addToEFVec(uint32_t val, uint32_t pos)
bool setVecBit(uint32_t val, uint32_t vecpos, uint32_t bbegin, uint32_t bend)
const std::vector< uint32_t > getEFVec() const
bool get(const std::string &name, T &value)
HltStoreSvc * m_HltStoreSvc
static int layer(const Identifier &id)
Values of different levels (failure returns 0)
static int wire(const Identifier &id)