NEEDS DOCUMENTATION. More...

#include <BarrelFastSimDedicatedSD.h>

Inheritance diagram for BarrelFastSimDedicatedSD:

Collaboration diagram for BarrelFastSimDedicatedSD:

Public Types
enum	LArHitTimeBins { HitTimeBinDefault = 0, HitTimeBinUniform = 1 }

Public Member Functions
	BarrelFastSimDedicatedSD (StoreGateSvc *, bool verbose)

	~BarrelFastSimDedicatedSD ()

void	ProcessSpot (const EnergySpot &spot, double weight) override final
	ProcessHitsMethod. More...

G4bool	ProcessHits (G4Step a_step, G4TouchableHistory ) override
	Main processing method. More...

G4bool	SimpleHit (const LArG4Identifier &lar_id, G4double time, G4double energy)
	First method translates to this - also for fast sims. More...

void	EndOfAthenaEvent (LArHitContainer *hitContnainer)
	End of athena event processing. More...

void	setupHelpers (const LArEM_ID EM, const LArFCAL_ID FCAL, const LArHEC_ID *HEC)
	Sets the ID helper pointers. More...

Identifier	ConvertID (const LArG4Identifier &a_ident) const
	Helper function for making "real" identifiers from LArG4Identifiers. More...

Protected Types
typedef std::set< LArHit *, LessHit >	hits_t

typedef std::map< G4int, hits_t * >	timeBins_t

Protected Member Functions
G4int	getTimeBin (G4double time) const
	Helper method for time info. More...

Protected Attributes
ILArCalculatorSvc *	m_calculator
	Member variable - the calculator we'll use. More...

G4int	m_numberInvalidHits
	Count the number of invalid hits. More...

LArG4SimpleSD::LArHitTimeBins	m_timeBinType

G4float	m_timeBinWidth
	Width of the time bins for summing hits - for the uniform binning. More...

timeBins_t	m_timeBins
	The map of hit sets binned in time. More...

const LArEM_ID *	m_larEmID
	Pointers to the identifier helpers. More...

const LArFCAL_ID *	m_larFcalID

const LArHEC_ID *	m_larHecID

Private Member Functions
	FRIEND_TEST (LArG4SimpleSDtest, ProcessHits)

	FRIEND_TEST (LArG4SimpleSDtest, EndOfAthenaEvent)

	FRIEND_TEST (LArG4SimpleSDtest, setupHelpers)

	FRIEND_TEST (LArG4SimpleSDtest, getTimeBin)

	FRIEND_TEST (LArG4SimpleSDtest, SimpleHit)

	FRIEND_TEST (LArG4SimpleSDtest, ConvertID)

	FRIEND_TEST (SDWrappertest, Initialize)

	FRIEND_TEST (SDWrappertest, EndOfAthenaEvent)

	FRIEND_TEST (SDWrappertest, addSD)

	FRIEND_TEST (SDWrappertest, addFastSimSD)

	FRIEND_TEST (SDWrappertest, ProcessHits)

Private Attributes
const EMBDetectorManager *	m_embManager

const EMBAccordionDetails *	m_accordionDetails

const GeoStraightAccSection *	m_absorberSections

LArG4Identifier	m_larID
	My LAr identifier. More...

Detailed Description

NEEDS DOCUMENTATION.

Definition at line 25 of file BarrelFastSimDedicatedSD.h.

Member Typedef Documentation

◆ hits_t

typedef std::set< LArHit*, LessHit > LArG4SimpleSD::hits_t

protectedinherited

Definition at line 113 of file LArG4SimpleSD.h.

◆ timeBins_t

typedef std::map< G4int, hits_t* > LArG4SimpleSD::timeBins_t

protectedinherited

Definition at line 119 of file LArG4SimpleSD.h.

Member Enumeration Documentation

◆ LArHitTimeBins

enum LArG4SimpleSD::LArHitTimeBins

inherited

Enumerator
HitTimeBinDefault
HitTimeBinUniform

Definition at line 50 of file LArG4SimpleSD.h.

   {
       HitTimeBinDefault = 0,
       HitTimeBinUniform = 1
   };

Constructor & Destructor Documentation

◆ BarrelFastSimDedicatedSD()

BarrelFastSimDedicatedSD::BarrelFastSimDedicatedSD	(	StoreGateSvc *	detStore,
		bool	verbose
	)

Definition at line 23 of file BarrelFastSimDedicatedSD.cxx.

   : IFastSimDedicatedSD("BarrelFastSimDedicatedSD", detStore)
   , m_embManager(nullptr)
   , m_accordionDetails(nullptr)
   , m_absorberSections(nullptr)
 {
   if (verbose) { G4cout << GetName() << "::initialize()" << G4endl; }
   if ( detStore->retrieve( m_embManager ).isFailure()  ){
     throw std::runtime_error("Could not retrieve EMB manager");
   }
   m_accordionDetails=m_embManager->getAccordionDetails();
   m_absorberSections=m_accordionDetails->getAbsorberSections();
 }

◆ ~BarrelFastSimDedicatedSD()

BarrelFastSimDedicatedSD::~BarrelFastSimDedicatedSD ( )

inline

Definition at line 34 of file BarrelFastSimDedicatedSD.h.

34 {}

Member Function Documentation

◆ ConvertID()

Identifier LArG4SimpleSD::ConvertID ( const LArG4Identifier & a_ident ) const

inherited

Helper function for making "real" identifiers from LArG4Identifiers.

◆ EndOfAthenaEvent()

void LArG4SimpleSD::EndOfAthenaEvent ( LArHitContainer * hitContnainer )

inherited

End of athena event processing.

◆ FRIEND_TEST() [1/11]

LArG4SimpleSD::FRIEND_TEST	(	LArG4SimpleSDtest	,
		ConvertID
	)

privateinherited

◆ FRIEND_TEST() [2/11]

LArG4SimpleSD::FRIEND_TEST	(	LArG4SimpleSDtest	,
		EndOfAthenaEvent
	)

privateinherited

◆ FRIEND_TEST() [3/11]

LArG4SimpleSD::FRIEND_TEST	(	LArG4SimpleSDtest	,
		getTimeBin
	)

privateinherited

◆ FRIEND_TEST() [4/11]

LArG4SimpleSD::FRIEND_TEST	(	LArG4SimpleSDtest	,
		ProcessHits
	)

privateinherited

◆ FRIEND_TEST() [5/11]

LArG4SimpleSD::FRIEND_TEST	(	LArG4SimpleSDtest	,
		setupHelpers
	)

privateinherited

◆ FRIEND_TEST() [6/11]

LArG4SimpleSD::FRIEND_TEST	(	LArG4SimpleSDtest	,
		SimpleHit
	)

privateinherited

◆ FRIEND_TEST() [7/11]

LArG4SimpleSD::FRIEND_TEST	(	SDWrappertest	,
		addFastSimSD
	)

privateinherited

◆ FRIEND_TEST() [8/11]

LArG4SimpleSD::FRIEND_TEST	(	SDWrappertest	,
		addSD
	)

privateinherited

◆ FRIEND_TEST() [9/11]

LArG4SimpleSD::FRIEND_TEST	(	SDWrappertest	,
		EndOfAthenaEvent
	)

privateinherited

◆ FRIEND_TEST() [10/11]

LArG4SimpleSD::FRIEND_TEST	(	SDWrappertest	,
		Initialize
	)

privateinherited

◆ FRIEND_TEST() [11/11]

LArG4SimpleSD::FRIEND_TEST	(	SDWrappertest	,
		ProcessHits
	)

privateinherited

◆ getTimeBin()

G4int LArG4SimpleSD::getTimeBin ( G4double time ) const

protectedinherited

Helper method for time info.

◆ ProcessHits()

G4bool LArG4SimpleSD::ProcessHits	(	G4Step *	a_step,
		G4TouchableHistory *
	)

overrideinherited

Main processing method.

◆ ProcessSpot()

void BarrelFastSimDedicatedSD::ProcessSpot	(	const EnergySpot &	spot,
		double	weight
	)

finaloverridevirtual

ProcessHitsMethod.

Process a single energy spot from a frozen shower. The appropriate region of the sensitive detector is calculated and a LArIdentifier is constructed

Implements IFastSimDedicatedSD.

Definition at line 38 of file BarrelFastSimDedicatedSD.cxx.

                                                                                  {
  
   // Fill the identifier.
   Point3D<double> globalPosition=spot.GetPosition();
  
   static const Transform3D xfNeg = Amg::EigenTransformToCLHEP(m_embManager->getDetectorRegion(0,0,0)->getAbsoluteTransform().inverse());
   static const Transform3D xfPos = Amg::EigenTransformToCLHEP(m_embManager->getDetectorRegion(1,0,0)->getAbsoluteTransform().inverse());
  
   Point3D<double> localPosition = globalPosition.z()<0 ? xfNeg*globalPosition : xfPos*globalPosition;
   int zIndex = globalPosition.z() <0 ? 0:1;
   double eta = localPosition.getEta();
   double phi = localPosition.getPhi();
   double r   = localPosition.perp();
  
   bool implementWaves=true;
  
   EMBDetectorManager::DetectorRegionConstIterator e;
   for (e=m_embManager->beginDetectorRegion();e!=m_embManager->endDetectorRegion();  ++e) {
     const EMBDetectorRegion *region = *e;
     if (region->getEndcapIndex()!=zIndex) continue;
     const EMBDetDescr       *regionDescriptor=region->getDescriptor();
     const CellBinning       & etaBinning=regionDescriptor->getEtaBinning();
     if (eta>etaBinning.getStart() && eta<etaBinning.getEnd()) {
       unsigned int etaIndex = int((eta - etaBinning.getStart())/etaBinning.getDelta()) + etaBinning.getFirstDivisionNumber();
       unsigned int phiIndex = 0;
       EMBCellConstLink cellPtr = region->getEMBCell(etaIndex,phiIndex);
       double rmax = cellPtr->getRLocal(EMBCell::BACK);
       double rmin = cellPtr->getRLocal(EMBCell::FRONT);
       if (r>rmin && r<rmax) {
         const CellBinning       & phiBinning=regionDescriptor->getPhiBinning();
         double minPhi=std::min(phiBinning.getStart(), phiBinning.getEnd());
         double maxPhi=std::max(phiBinning.getStart(), phiBinning.getEnd());
  
         while (phi<minPhi)   phi += 2*M_PI;
         while (phi>maxPhi)   phi -= 2*M_PI;
         unsigned int regionIndex = region->getRegionIndex();
         unsigned int samplingIndex = region->getSamplingIndex();
         unsigned int phiIndex      =  int((phi - phiBinning.getStart())/phiBinning.getDelta()) + phiBinning.getFirstDivisionNumber();
  
         // Now get the nominal cell and see if the indices need to be corrected for the accordion waves:
         if (samplingIndex!=0) {
           EMBCellConstLink cellPtr=region->getEMBCell(etaIndex,phiIndex);
  
           double delta=-2.0*M_PI/1024/2;
  
           double phiLocalUpper=cellPtr->getPhiLocalUpper()-delta;
           double phiLocalLower=cellPtr->getPhiLocalLower()-delta;
  
           while (phiLocalLower<0)      phiLocalLower += 2.0*M_PI;
           while (phiLocalUpper<0)      phiLocalUpper += 2.0*M_PI;
           while (phiLocalLower>2*M_PI) phiLocalLower -= 2.0*M_PI;
           while (phiLocalUpper>2*M_PI) phiLocalUpper -= 2.0*M_PI;
  
           int accordionIndexUpper=int(phiLocalUpper/(2*M_PI)*1024+0.25); if (accordionIndexUpper==1024) accordionIndexUpper=0;
           int accordionIndexLower=int(phiLocalLower/(2*M_PI)*1024+0.25); if (accordionIndexLower==1024) accordionIndexLower=0;
  
           Point3D<double> A0     =
             ( Point3D<double>(m_absorberSections->XCent(accordionIndexLower,1),m_absorberSections->YCent(accordionIndexLower,1),0)+
               Point3D<double>(m_absorberSections->XCent(accordionIndexUpper,1),m_absorberSections->YCent(accordionIndexUpper,1),0))/2.0;
  
           Point3D<double> A1=
             (Point3D<double>(m_absorberSections->XCent(accordionIndexLower,12),m_absorberSections->YCent(accordionIndexLower,12),0)+
              Point3D<double>(m_absorberSections->XCent(accordionIndexUpper,12),m_absorberSections->YCent(accordionIndexUpper,12),0))/2.0-A0;
  
           Point3D<double> P0=Point3D<double>(localPosition.x(), localPosition.y(),0)-A0;
           int stackIndex=int(A1.dot(P0)/A1.mag2()*22.0 + 3.0)/2 ;
  
           if (stackIndex<0 || stackIndex>13) {
             G4cout << "Warning, bad stack index " << stackIndex << ' ' << rmin << ' ' << r << ' ' << A0.perp() << ' ' << A0 <<  ' ' << A1 << ' ' << P0 << G4endl;
             if (implementWaves) return;
           } else {
             double xcent      = m_absorberSections->XCent     (accordionIndexLower,stackIndex);
             double ycent      = m_absorberSections->YCent     (accordionIndexLower,stackIndex);
             double cosU       = m_absorberSections->Cosu      (accordionIndexLower,stackIndex);
             double sinU       = m_absorberSections->Sinu      (accordionIndexLower,stackIndex);
             double halfLength = m_absorberSections->HalfLength(accordionIndexLower,stackIndex);
  
             Point3D<double> u(xcent-halfLength*cosU, ycent-halfLength*sinU,0);
             Point3D<double> v(xcent+halfLength*cosU, ycent+halfLength*sinU,0);
             Point3D<double> x=v-u,y=Hep3Vector(localPosition.x(),localPosition.y(),0)-u;
             if ((zIndex==0 && x.cross(y).z()>0) || (zIndex==1&& x.cross(y).z()<0)) {
               if (implementWaves) {
                 if (phiIndex==0) phiIndex=phiBinning.getNumDivisions()-1;
                 else phiIndex--;
               }
             } else {
               double xcent      = m_absorberSections->XCent     (accordionIndexUpper,stackIndex);
               double ycent      = m_absorberSections->YCent     (accordionIndexUpper,stackIndex);
               double cosU       = m_absorberSections->Cosu      (accordionIndexUpper,stackIndex);
               double sinU       = m_absorberSections->Sinu      (accordionIndexUpper,stackIndex);
               double halfLength = m_absorberSections->HalfLength(accordionIndexUpper,stackIndex);
  
               Point3D<double> u(xcent-halfLength*cosU, ycent-halfLength*sinU,0);
               Point3D<double> v(xcent+halfLength*cosU, ycent+halfLength*sinU,0);
               Point3D<double> x=v-u,y=Hep3Vector(localPosition.x(),localPosition.y(),0)-u;
               if ((zIndex==0 && x.cross(y).z()<0) || (zIndex==1 && x.cross(y).z()>0)) {
                 if (implementWaves) {
                   if (phiIndex==phiBinning.getNumDivisions()-1) phiIndex=0;
                   else phiIndex++;
                 }
               }
             }
           }
         }
  
         //static LArG4Identifier id;
         m_larID.clear();
         m_larID << 4          // LArCalorimeter
                 << 1          // LArEM
                 << ((zIndex==0) ? -1:1)
                 << samplingIndex
                 << regionIndex
                 << etaIndex
                 << phiIndex;
         // call process to add this to the collection 
         SimpleHit(m_larID, spot.GetTime(), spot.GetEnergy()*weight);
  
         return;
       }
     }
   }
  
 }

◆ setupHelpers()

void LArG4SimpleSD::setupHelpers	(	const LArEM_ID *	EM,
		const LArFCAL_ID *	FCAL,
		const LArHEC_ID *	HEC
	)

inlineinherited

Sets the ID helper pointers.

Definition at line 77 of file LArG4SimpleSD.h.

                                             {
     m_larEmID = EM;
     m_larFcalID = FCAL;
     m_larHecID = HEC;
   }

◆ SimpleHit()

G4bool LArG4SimpleSD::SimpleHit	(	const LArG4Identifier &	lar_id,
		G4double	time,
		G4double	energy
	)

inherited

First method translates to this - also for fast sims.

Member Data Documentation

◆ m_absorberSections

const GeoStraightAccSection* BarrelFastSimDedicatedSD::m_absorberSections

private

Definition at line 43 of file BarrelFastSimDedicatedSD.h.

◆ m_accordionDetails

const EMBAccordionDetails* BarrelFastSimDedicatedSD::m_accordionDetails

private

Definition at line 42 of file BarrelFastSimDedicatedSD.h.

◆ m_calculator

ILArCalculatorSvc* LArG4SimpleSD::m_calculator

protectedinherited

Member variable - the calculator we'll use.

Definition at line 93 of file LArG4SimpleSD.h.

◆ m_embManager

const EMBDetectorManager* BarrelFastSimDedicatedSD::m_embManager

private

Definition at line 41 of file BarrelFastSimDedicatedSD.h.

◆ m_larEmID

const LArEM_ID* LArG4SimpleSD::m_larEmID

protectedinherited

Pointers to the identifier helpers.

Definition at line 143 of file LArG4SimpleSD.h.

◆ m_larFcalID

const LArFCAL_ID* LArG4SimpleSD::m_larFcalID

protectedinherited

Definition at line 144 of file LArG4SimpleSD.h.

◆ m_larHecID

const LArHEC_ID* LArG4SimpleSD::m_larHecID

protectedinherited

Definition at line 145 of file LArG4SimpleSD.h.

◆ m_larID

LArG4Identifier BarrelFastSimDedicatedSD::m_larID

private

My LAr identifier.

Definition at line 46 of file BarrelFastSimDedicatedSD.h.

◆ m_numberInvalidHits

G4int LArG4SimpleSD::m_numberInvalidHits

protectedinherited

Count the number of invalid hits.

Definition at line 96 of file LArG4SimpleSD.h.

◆ m_timeBins

timeBins_t LArG4SimpleSD::m_timeBins

protectedinherited

The map of hit sets binned in time.

Definition at line 140 of file LArG4SimpleSD.h.

◆ m_timeBinType

LArG4SimpleSD::LArHitTimeBins LArG4SimpleSD::m_timeBinType

protectedinherited

Definition at line 134 of file LArG4SimpleSD.h.

◆ m_timeBinWidth

G4float LArG4SimpleSD::m_timeBinWidth

protectedinherited

Width of the time bins for summing hits - for the uniform binning.

Definition at line 137 of file LArG4SimpleSD.h.

The documentation for this class was generated from the following files:

Public Types

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ hits_t

◆ timeBins_t

Member Enumeration Documentation

◆ LArHitTimeBins

Constructor & Destructor Documentation

◆ BarrelFastSimDedicatedSD()

◆ ~BarrelFastSimDedicatedSD()

Member Function Documentation

◆ ConvertID()

◆ EndOfAthenaEvent()

◆ FRIEND_TEST() [1/11]

◆ FRIEND_TEST() [2/11]

◆ FRIEND_TEST() [3/11]

◆ FRIEND_TEST() [4/11]

◆ FRIEND_TEST() [5/11]

◆ FRIEND_TEST() [6/11]

◆ FRIEND_TEST() [7/11]

◆ FRIEND_TEST() [8/11]

◆ FRIEND_TEST() [9/11]

◆ FRIEND_TEST() [10/11]

◆ FRIEND_TEST() [11/11]

◆ getTimeBin()

◆ ProcessHits()

◆ ProcessSpot()

◆ setupHelpers()

◆ SimpleHit()

Member Data Documentation

◆ m_absorberSections

◆ m_accordionDetails

◆ m_calculator

◆ m_embManager

◆ m_larEmID

◆ m_larFcalID

◆ m_larHecID

◆ m_larID

◆ m_numberInvalidHits

◆ m_timeBins

◆ m_timeBinType

◆ m_timeBinWidth