|
ATLAS Offline Software
|
Go to the documentation of this file.
17 #ifndef TRIGT2CALO_T2GEOMETRY
18 #define TRIGT2CALO_T2GEOMETRY
43 int EtaPhiRange(
const int nCaloType,
const int nCaloSamp,
44 const double& energyEta,
const double& energyPhi)
override;
54 const double& etaCell,
const double& phiCell)
override;
63 const double& etaCell,
const double& phiCell)
override;
72 const double& etaCell,
const double& phiCell)
override;
155 const int& nStripPhi,
const double& energyEta,
156 const double& energyPhi,
const int index_size)
override;
203 const int& nStripEta,
const int& nStripPhi,
204 const double& energyEta,
const double& energyPhi,
const int index_size){
206 double aeta = fabs(energyEta);
209 double CellEtaMin=0.;
210 double CellEtaMax=0.;
211 double CellPhiMin=0.;
212 double CellPhiMax=0.;
214 CellEtaMin = aeta - dgra*nStripEta*0.5;
216 if ( CellEtaMin < 0. ) {
218 CellEtaMin = fabs(CellEtaMin );
221 " aeta " << aeta <<
" dgra " << dgra <<
" netaregions " << netaregions
222 <<
" First CellEtaMin " << CellEtaMin );
240 for (irmin=0;irmin<netaregions;irmin++) {
247 if(CellEtaMin - trunc(CellEtaMin) > 0.7 ) CellEtaMin++;
248 CellEtaMin = dgra1*trunc(CellEtaMin) +
m_etareg[
layer][irmin];
250 CellEtaMin = CellEtaMin*EtaMinSign;
251 CellEtaMax = aeta + dgra*nStripEta * 0.5;
254 for (irmax=irmin;irmax<netaregions;irmax++) {
259 if (CellEtaMax-trunc(CellEtaMax) >= 0.3) CellEtaMax++;
271 CellEtaMin=-CellEtaMax;
284 CellPhiMin = energyPhi - dgrap*nStripPhi*0.5;
287 for (irmin=0;irmin<nphiregions;irmin++) {
293 CellPhiMin = dgra1*rint(CellPhiMin/dgra1);
294 CellPhiMax = CellPhiMin + dgrap*nStripPhi;
299 if (CellPhiMax >
M_PI) CellPhiMax -= 2*
M_PI;
300 if (CellPhiMin <= -
M_PI) CellPhiMin += 2*
M_PI;
301 if (CellPhiMin >
M_PI) CellPhiMin -= 2*
M_PI;
302 if (CellPhiMax <= -
M_PI) CellPhiMax += 2*
M_PI;
311 else if (index_size == 1)
318 else if (index_size == 2)
331 const double& energyEta,
const double& energyPhi){
334 if (nCaloType < 0 || nCaloType > 1) {
338 if (nCaloSamp < 0 || nCaloSamp > 3 || (nCaloSamp == 3 && nCaloType == 1)) {
342 int layer = nCaloType*4+nCaloSamp;
394 const int nCaloSamp,
const double& etaCell,
const double& phiCell){
397 if (nCaloType < 0 || nCaloType > 1) {
401 if (nCaloSamp < 0 || nCaloSamp > 3 || (nCaloSamp == 3 && nCaloType == 1)) {
406 int layer = nCaloType*4+nCaloSamp;
422 if ( IetaPass == 1 && IphiPass == 1 )
return true;
427 const int nCaloSamp,
const double& etaCell,
const double& phiCell){
429 if (nCaloType < 0 || nCaloType > 1) {
433 if (nCaloSamp < 0 || nCaloSamp > 3 || (nCaloSamp == 3 && nCaloType == 1)) {
437 int layer = nCaloType*4+nCaloSamp;
453 if ( IetaPass == 1 && IphiPass == 1 )
return true;
458 const int nCaloSamp,
const double& etaCell,
const double& phiCell){
460 if (nCaloType < 0 || nCaloType > 1) {
464 if (nCaloSamp < 0 || nCaloSamp > 3 || (nCaloSamp == 3 && nCaloType == 1)) {
468 int layer = nCaloType*4+nCaloSamp;
485 if ( IetaPass == 1 && IphiPass == 1 )
return true;
::StatusCode StatusCode
StatusCode definition for legacy code.