59 return StatusCode::SUCCESS;
78 return StatusCode::SUCCESS;
89 std::vector<CaloAffectedRegionInfo>::const_iterator reg1 = affRegVec->begin();
90 std::vector<CaloAffectedRegionInfo>::const_iterator reg2 = affRegVec->end();
91 for (;reg1 != reg2; ++reg1) {
94 if (problem == 0)
continue;
134 ATH_MSG_DEBUG(
"EMBAPS " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
139 ATH_MSG_DEBUG(
"EMBCPS " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
141 phi+=(2*TMath::Pi()/32);
142 if(fabs(phi)<1.e-4) phi=0.;
147 while(phi<TMath::Pi()){
151 ATH_MSG_DEBUG(
"EMBAPS " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
156 ATH_MSG_DEBUG(
"EMBCPS " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
158 phi+=(2*TMath::Pi()/32);
159 if(fabs(phi)<1.e-4) phi=0.;
166 ATH_MSG_DEBUG(
"EMBAPS " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
171 ATH_MSG_DEBUG(
"EMBCPS " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
173 phi+=(2*TMath::Pi()/32);
174 if(fabs(phi)<1.e-4) phi=0.;
193 ATH_MSG_DEBUG(
"EMBA " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
198 ATH_MSG_DEBUG(
"EMBC " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
200 phi+=(2*TMath::Pi()/128);
201 if(fabs(phi)<1.e-4) phi=0.;
206 while(phi<TMath::Pi()){
210 ATH_MSG_DEBUG(
"EMBA " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
215 ATH_MSG_DEBUG(
"EMBC " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
217 phi+=(2*TMath::Pi()/128);
218 if(fabs(phi)<1.e-4) phi=0.;
225 ATH_MSG_DEBUG(
"EMBA " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
230 ATH_MSG_DEBUG(
"EMBC " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
232 phi+=(2*TMath::Pi()/128);
233 if(fabs(phi)<1.e-4) phi=0.;
252 ATH_MSG_DEBUG(
"EMECAPS " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
257 ATH_MSG_DEBUG(
"EMECCPS " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
259 phi+=(2*TMath::Pi()/64);
260 if(fabs(phi)<1.e-4) phi=0.;
265 while(phi<TMath::Pi()){
269 ATH_MSG_DEBUG(
"EMECAPS " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
274 ATH_MSG_DEBUG(
"EMECCPS " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
276 phi+=(2*TMath::Pi()/64);
277 if(fabs(phi)<1.e-4) phi=0.01;
279 phi = -TMath::Pi()+0.01;
284 ATH_MSG_DEBUG(
"EMECAPS " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
289 ATH_MSG_DEBUG(
"EMECCPS " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
291 phi+=(2*TMath::Pi()/64);
292 if(fabs(phi)<1.e-4) phi=0.;
295 if(fabs(eta)<2.5)eta+=0.2;
296 if(fabs(eta)>=2.5)eta+=0.3;
312 ATH_MSG_DEBUG(
"EMECA " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
317 ATH_MSG_DEBUG(
"EMECC " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
319 phi+=(2*TMath::Pi()/256);
320 if(fabs(phi)<1.e-4) phi=0.;
325 while(phi<TMath::Pi()){
329 ATH_MSG_DEBUG(
"EMECA " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
334 ATH_MSG_DEBUG(
"EMECC " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
336 phi+=(2*TMath::Pi()/256);
337 if(fabs(phi)<1.e-4) phi=0.;
339 phi = -TMath::Pi()+0.01;
344 ATH_MSG_DEBUG(
"EMECA " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
349 ATH_MSG_DEBUG(
"EMECC " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
351 phi+=(2*TMath::Pi()/256);
352 if(fabs(phi)<1.e-4) phi=0.;
355 if(fabs(eta)>=2.8) eta+=0.4;
356 else if(fabs(eta)>=2.5)eta+=0.3;
362 if(layermin>=8 && layermin<=11){
388 fill(hec_layer_name,etaPOS,phi,problem);
389 ATH_MSG_DEBUG(
"HECA " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
393 fill(hec_layer_name,etaNEG,phi,problem);
394 ATH_MSG_DEBUG(
"HECC " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << region );
396 phi+=(2*TMath::Pi()/32);
397 if(fabs(phi)<1.e-4) phi=0.;
402 while(phi<TMath::Pi()){
405 fill(hec_layer_name,etaPOS,phi,problem);
406 ATH_MSG_DEBUG(
"HECA " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
410 fill(hec_layer_name,etaNEG,phi,problem);
411 ATH_MSG_DEBUG(
"HECC " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << region );
413 phi+=(2*TMath::Pi()/32);
414 if(fabs(phi)<1.e-4) phi=0.;
420 fill(hec_layer_name,etaPOS,phi,problem);
421 ATH_MSG_DEBUG(
"HECA " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
425 fill(hec_layer_name,etaNEG,phi,problem);
426 ATH_MSG_DEBUG(
"HECC " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << region );
428 phi+=(2*TMath::Pi()/32);
429 if(fabs(phi)<1.e-4) phi=0.;
435 ATH_MSG_DEBUG(
"layermin=" << layermin <<
" layermin-21=" << layermin-21);
436 if(layermin>=21 && layermin<=23){
439 switch(layermin-21) {
459 fill(fcal_layer_name,etaPOS,phi,problem);
460 ATH_MSG_DEBUG(
"FCALA " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
464 fill(fcal_layer_name,etaNEG,phi,problem);
465 ATH_MSG_DEBUG(
"FCALC " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
467 phi+=(2*TMath::Pi()/16);
468 if(fabs(phi)<1.e-4) phi=0.;
473 while(phi<TMath::Pi()){
476 fill(fcal_layer_name,etaPOS,phi,problem);
477 ATH_MSG_DEBUG(
"FCALA " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
481 fill(fcal_layer_name,etaNEG,phi,problem);
482 ATH_MSG_DEBUG(
"FCALC " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
484 phi+=(2*TMath::Pi()/16);
485 if(fabs(phi)<1.e-4) phi=0.;
487 phi = -TMath::Pi()+0.1;
491 fill(fcal_layer_name,etaPOS,phi,problem);
492 ATH_MSG_DEBUG(
"FCALA " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
496 fill(fcal_layer_name,etaNEG,phi,problem);
497 ATH_MSG_DEBUG(
"FCALC " << eta <<
" " << phi <<
" " << problem <<
" " << layermin <<
" " << layermax <<
" " << phimin <<
" " << phimax <<
" " <<
etamin <<
" " << etamax <<
" " << region );
499 phi+=(2*TMath::Pi()/16);
500 if(fabs(phi)<1.e-4) phi=0.;
506 return StatusCode::SUCCESS;