IEnvelopeDefSvc.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
*/

// IEnvelopeDefSvc.h, (c) ATLAS Detector software
#ifndef IENVELOPEDEFSVC_H
#define IENVELOPEDEFSVC_H
 
// Gaudi
#include "GaudiKernel/IInterface.h"
 
// Atlas Regions
#include "AtlasDetDescr/AtlasRegion.h"
 
// SubDetectorEnvelopes includes
#include "SubDetectorEnvelopes/RZPair.h"
 
#include "CxxUtils/checker_macros.h"
 
class IEnvelopeDefSvc : virtual public IInterface {
 
  public:
    // interface definition
    DeclareInterfaceID(IEnvelopeDefSvc, 1, 0);

    virtual ~IEnvelopeDefSvc() {};

    virtual const RZPairVector &getRZBoundary( AtlasDetDescr::AtlasRegion region ) const = 0;

    virtual const RZPairVector &getRPositiveZBoundary( AtlasDetDescr::AtlasRegion region ) const = 0;

    const RZPairVector &getBeamPipeRZBoundary() const { return getRZBoundary(AtlasDetDescr::fAtlasForward); }
    const RZPairVector &getInDetRZBoundary()    const { return getRZBoundary(AtlasDetDescr::fAtlasID);      }
    const RZPairVector &getCaloRZBoundary()     const { return getRZBoundary(AtlasDetDescr::fAtlasCalo);    }
    const RZPairVector &getMuonRZBoundary()     const { return getRZBoundary(AtlasDetDescr::fAtlasMS);      }
    const RZPairVector &getCavernRZBoundary()   const { return getRZBoundary(AtlasDetDescr::fAtlasCavern);  }
 
  protected:
    inline void mirrorRZ( const RZPairVector &srcRZ, RZPairVector &dstRZ ) const ;
 
};

inline void IEnvelopeDefSvc::mirrorRZ( const RZPairVector &srcRZ, RZPairVector &dstRZ ) const
{
  int numPosPairs = srcRZ.size();
  // the mirrored envelope will have exactly twice as many entries as the given srcRZ
  dstRZ.resize(2*numPosPairs);
 
  // loop over all positive (r,z) pairs
  for ( int curNum = 0; curNum<numPosPairs; curNum++) {
    double curR = srcRZ[curNum].first;
    double curZ = srcRZ[curNum].second;
 
    // debugging output:
    //std::cout << "Envelope: pos=" << curNum << " r=" << curR << " z="<< curZ << std::endl;
 
    // fill in the z<0 side
    dstRZ[numPosPairs-curNum-1] = RZPair(curR, -curZ);
    // fill in the z>0 side
    dstRZ[numPosPairs+curNum]   = RZPair(curR, curZ);
  }
 
  return;
}
 
#endif // IENVELOPEDEFSVC_H