ATLAS Offline Software
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Protected Types | Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes | List of all members
BarrelFastSimDedicatedSD Class Reference

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.
G4bool ProcessHits (G4Step *a_step, G4TouchableHistory *) override
 Main processing method.
G4bool SimpleHit (const LArG4Identifier &lar_id, G4double time, G4double energy)
 First method translates to this - also for fast sims.
void EndOfAthenaEvent (LArHitContainer *hitContnainer)
 End of athena event processing.
void setupHelpers (const LArEM_ID *EM, const LArFCAL_ID *FCAL, const LArHEC_ID *HEC)
 Sets the ID helper pointers.
Identifier ConvertID (const LArG4Identifier &a_ident) const
 Helper function for making "real" identifiers from LArG4Identifiers.

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.

Protected Attributes

ILArCalculatorSvc * m_calculator
 Member variable - the calculator we'll use.
G4int m_numberInvalidHits
 Count the number of invalid hits.
LArG4SimpleSD::LArHitTimeBins m_timeBinType
G4float m_timeBinWidth
 Width of the time bins for summing hits - for the uniform binning.
timeBins_t m_timeBins
 The map of hit sets binned in time.
const LArEM_ID * m_larEmID
 Pointers to the identifier helpers.
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.

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.

51 {
52 HitTimeBinDefault = 0,
53 HitTimeBinUniform = 1
54 };

Constructor & Destructor Documentation

◆ BarrelFastSimDedicatedSD()

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

Definition at line 23 of file BarrelFastSimDedicatedSD.cxx.

24 : IFastSimDedicatedSD("BarrelFastSimDedicatedSD", detStore)
25 , m_embManager(nullptr)
26 , m_accordionDetails(nullptr)
27 , m_absorberSections(nullptr)
28{
29 if (verbose) { G4cout << GetName() << "::initialize()" << G4endl; }
30 if ( detStore->retrieve( m_embManager ).isFailure() ){
31 throw std::runtime_error("Could not retrieve EMB manager");
32 }
33 m_accordionDetails=m_embManager->getAccordionDetails();
34 m_absorberSections=m_accordionDetails->getAbsorberSections();
35}
BarrelFastSimDedicatedSD::m_embManager
const EMBDetectorManager * m_embManager
Definition BarrelFastSimDedicatedSD.h:41
BarrelFastSimDedicatedSD::m_absorberSections
const GeoStraightAccSection * m_absorberSections
Definition BarrelFastSimDedicatedSD.h:43
BarrelFastSimDedicatedSD::m_accordionDetails
const EMBAccordionDetails * m_accordionDetails
Definition BarrelFastSimDedicatedSD.h:42
IFastSimDedicatedSD::IFastSimDedicatedSD
IFastSimDedicatedSD(const std::string &name, StoreGateSvc *detStore)
Simple constructor and destructor.
Definition IFastSimDedicatedSD.h:18
verbose
bool verbose
Definition hcg.cxx:73

◆ ~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.

38 {
39
40 // Fill the identifier.
41 Point3D<double> globalPosition=spot.GetPosition();
42
43 static const Transform3D xfNeg = Amg::EigenTransformToCLHEP(m_embManager->getDetectorRegion(0,0,0)->getAbsoluteTransform().inverse());
44 static const Transform3D xfPos = Amg::EigenTransformToCLHEP(m_embManager->getDetectorRegion(1,0,0)->getAbsoluteTransform().inverse());
45
46 Point3D<double> localPosition = globalPosition.z()<0 ? xfNeg*globalPosition : xfPos*globalPosition;
47 int zIndex = globalPosition.z() <0 ? 0:1;
48 double eta = localPosition.getEta();
49 double phi = localPosition.getPhi();
50 double r = localPosition.perp();
51
52 bool implementWaves=true;
53
54 EMBDetectorManager::DetectorRegionConstIterator e;
55 for (e=m_embManager->beginDetectorRegion();e!=m_embManager->endDetectorRegion(); ++e) {
56 const EMBDetectorRegion *region = *e;
57 if (region->getEndcapIndex()!=zIndex) continue;
58 const EMBDetDescr *regionDescriptor=region->getDescriptor();
59 const CellBinning & etaBinning=regionDescriptor->getEtaBinning();
60 if (eta>etaBinning.getStart() && eta<etaBinning.getEnd()) {
61 unsigned int etaIndex = int((eta - etaBinning.getStart())/etaBinning.getDelta()) + etaBinning.getFirstDivisionNumber();
62 unsigned int phiIndex = 0;
63 EMBCellConstLink cellPtr = region->getEMBCell(etaIndex,phiIndex);
64 double rmax = cellPtr->getRLocal(EMBCell::BACK);
65 double rmin = cellPtr->getRLocal(EMBCell::FRONT);
66 if (r>rmin && r<rmax) {
67 const CellBinning & phiBinning=regionDescriptor->getPhiBinning();
68 double minPhi=std::min(phiBinning.getStart(), phiBinning.getEnd());
69 double maxPhi=std::max(phiBinning.getStart(), phiBinning.getEnd());
70
71 while (phi<minPhi) phi += 2*M_PI;
72 while (phi>maxPhi) phi -= 2*M_PI;
73 unsigned int regionIndex = region->getRegionIndex();
74 unsigned int samplingIndex = region->getSamplingIndex();
75 unsigned int phiIndex = int((phi - phiBinning.getStart())/phiBinning.getDelta()) + phiBinning.getFirstDivisionNumber();
76
77 // Now get the nominal cell and see if the indices need to be corrected for the accordion waves:
78 if (samplingIndex!=0) {
79 EMBCellConstLink cellPtr=region->getEMBCell(etaIndex,phiIndex);
80
81 double delta=-2.0*M_PI/1024/2;
82
83 double phiLocalUpper=cellPtr->getPhiLocalUpper()-delta;
84 double phiLocalLower=cellPtr->getPhiLocalLower()-delta;
85
86 while (phiLocalLower<0) phiLocalLower += 2.0*M_PI;
87 while (phiLocalUpper<0) phiLocalUpper += 2.0*M_PI;
88 while (phiLocalLower>2*M_PI) phiLocalLower -= 2.0*M_PI;
89 while (phiLocalUpper>2*M_PI) phiLocalUpper -= 2.0*M_PI;
90
91 int accordionIndexUpper=int(phiLocalUpper/(2*M_PI)*1024+0.25); if (accordionIndexUpper==1024) accordionIndexUpper=0;
92 int accordionIndexLower=int(phiLocalLower/(2*M_PI)*1024+0.25); if (accordionIndexLower==1024) accordionIndexLower=0;
93
94 Point3D<double> A0 =
95 ( Point3D<double>(m_absorberSections->XCent(accordionIndexLower,1),m_absorberSections->YCent(accordionIndexLower,1),0)+
96 Point3D<double>(m_absorberSections->XCent(accordionIndexUpper,1),m_absorberSections->YCent(accordionIndexUpper,1),0))/2.0;
97
98 Point3D<double> A1=
99 (Point3D<double>(m_absorberSections->XCent(accordionIndexLower,12),m_absorberSections->YCent(accordionIndexLower,12),0)+
100 Point3D<double>(m_absorberSections->XCent(accordionIndexUpper,12),m_absorberSections->YCent(accordionIndexUpper,12),0))/2.0-A0;
101
102 Point3D<double> P0=Point3D<double>(localPosition.x(), localPosition.y(),0)-A0;
103 int stackIndex=int(A1.dot(P0)/A1.mag2()*22.0 + 3.0)/2 ;
104
105 if (stackIndex<0 || stackIndex>13) {
106 G4cout << "Warning, bad stack index " << stackIndex << ' ' << rmin << ' ' << r << ' ' << A0.perp() << ' ' << A0 << ' ' << A1 << ' ' << P0 << G4endl;
107 if (implementWaves) return;
108 } else {
109 double xcent = m_absorberSections->XCent (accordionIndexLower,stackIndex);
110 double ycent = m_absorberSections->YCent (accordionIndexLower,stackIndex);
111 double cosU = m_absorberSections->Cosu (accordionIndexLower,stackIndex);
112 double sinU = m_absorberSections->Sinu (accordionIndexLower,stackIndex);
113 double halfLength = m_absorberSections->HalfLength(accordionIndexLower,stackIndex);
114
115 Point3D<double> u(xcent-halfLength*cosU, ycent-halfLength*sinU,0);
116 Point3D<double> v(xcent+halfLength*cosU, ycent+halfLength*sinU,0);
117 Point3D<double> x=v-u,y=Hep3Vector(localPosition.x(),localPosition.y(),0)-u;
118 if ((zIndex==0 && x.cross(y).z()>0) || (zIndex==1&& x.cross(y).z()<0)) {
119 if (implementWaves) {
120 if (phiIndex==0) phiIndex=phiBinning.getNumDivisions()-1;
121 else phiIndex--;
122 }
123 } else {
124 double xcent = m_absorberSections->XCent (accordionIndexUpper,stackIndex);
125 double ycent = m_absorberSections->YCent (accordionIndexUpper,stackIndex);
126 double cosU = m_absorberSections->Cosu (accordionIndexUpper,stackIndex);
127 double sinU = m_absorberSections->Sinu (accordionIndexUpper,stackIndex);
128 double halfLength = m_absorberSections->HalfLength(accordionIndexUpper,stackIndex);
129
130 Point3D<double> u(xcent-halfLength*cosU, ycent-halfLength*sinU,0);
131 Point3D<double> v(xcent+halfLength*cosU, ycent+halfLength*sinU,0);
132 Point3D<double> x=v-u,y=Hep3Vector(localPosition.x(),localPosition.y(),0)-u;
133 if ((zIndex==0 && x.cross(y).z()<0) || (zIndex==1 && x.cross(y).z()>0)) {
134 if (implementWaves) {
135 if (phiIndex==phiBinning.getNumDivisions()-1) phiIndex=0;
136 else phiIndex++;
137 }
138 }
139 }
140 }
141 }
142
143 //static LArG4Identifier id;
144 m_larID.clear();
145 m_larID << 4 // LArCalorimeter
146 << 1 // LArEM
147 << ((zIndex==0) ? -1:1)
148 << samplingIndex
149 << regionIndex
150 << etaIndex
151 << phiIndex;
152 // call process to add this to the collection
153 SimpleHit(m_larID, spot.GetTime(), spot.GetEnergy()*weight);
154
155 return;
156 }
157 }
158 }
159
160}
M_PI
#define M_PI
Definition ActiveFraction.h:11
eta
Scalar eta() const
pseudorapidity method
Definition AmgMatrixBasePlugin.h:83
phi
Scalar phi() const
phi method
Definition AmgMatrixBasePlugin.h:67
Transform3D
Eigen::Affine3d Transform3D
Definition GeoPrimitives.h:46
y
#define y
x
#define x
BarrelFastSimDedicatedSD::m_larID
LArG4Identifier m_larID
My LAr identifier.
Definition BarrelFastSimDedicatedSD.h:46
EMBCell::BACK
@ BACK
Definition EMBCell.h:33
EMBCell::FRONT
@ FRONT
Definition EMBCell.h:33
EMBDetDescr::getPhiBinning
const CellBinning & getPhiBinning() const
The Binning in Phi.
Definition EMBDetDescr.h:177
EMBDetDescr::getEtaBinning
const CellBinning & getEtaBinning() const
The Binning in Eta.
Definition EMBDetDescr.h:184
EMBDetectorManager::DetectorRegionConstIterator
std::vector< constEMBDetectorRegion * >::const_iterator DetectorRegionConstIterator
Definition EMBDetectorManager.h:35
EMBDetectorRegion::getDescriptor
const EMBDetDescr * getDescriptor() const
Returns the Descriptor for this region.
Definition EMBDetectorRegion.h:112
EMBDetectorRegion::getEMBCell
EMBCellConstLink getEMBCell(unsigned int ieta, unsigned int iphi) const
Access to Cells.
Definition EMBDetectorRegion.cxx:22
EMBDetectorRegion::getSamplingIndex
unsigned int getSamplingIndex() const
Returns the Sampling Layer Index.
Definition EMBDetectorRegion.h:119
EMBDetectorRegion::getEndcapIndex
EMBDetectorRegion::DetectorSide getEndcapIndex() const
The endcap index.
Definition EMBDetectorRegion.h:163
EMBDetectorRegion::getRegionIndex
unsigned int getRegionIndex() const
Returns the Region Index.
Definition EMBDetectorRegion.h:126
EnergySpot::GetPosition
G4ThreeVector GetPosition() const
Definition EnergySpot.h:40
EnergySpot::GetEnergy
G4double GetEnergy() const
Definition EnergySpot.h:37
EnergySpot::GetTime
G4double GetTime() const
Definition EnergySpot.h:43
LArG4SimpleSD::SimpleHit
G4bool SimpleHit(const LArG4Identifier &lar_id, G4double time, G4double energy)
First method translates to this - also for fast sims.
r
int r
Definition globals.cxx:22
Amg::EigenTransformToCLHEP
HepGeom::Transform3D EigenTransformToCLHEP(const Amg::Transform3D &eigenTransf)
Converts an Eigen-based Amg::Transform3D into a CLHEP-based HepGeom::Transform3D.
Definition CLHEPtoEigenConverter.h:120
Trk::u
@ u
Enums for curvilinear frames.
Definition ParamDefs.h:77
eflowRec::phiIndex
unsigned int phiIndex(float phi, float binsize)
calculate phi index for a given phi
Definition EtaPhiLUT.cxx:23

◆ 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.

79 {
80 m_larEmID = EM;
81 m_larFcalID = FCAL;
82 m_larHecID = HEC;
83 }
LArG4SimpleSD::m_larHecID
const LArHEC_ID * m_larHecID
Definition LArG4SimpleSD.h:145
LArG4SimpleSD::m_larFcalID
const LArFCAL_ID * m_larFcalID
Definition LArG4SimpleSD.h:144
LArG4SimpleSD::m_larEmID
const LArEM_ID * m_larEmID
Pointers to the identifier helpers.
Definition LArG4SimpleSD.h:143

◆ 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:
Generated on for ATLAS Offline Software by doxygen 1.14.0