ATLAS Offline Software
Loading...
Searching...
No Matches
Public Types | Public Member Functions | 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 *, std::string hitCollectionName, bool verbose)
 ~BarrelFastSimDedicatedSD ()
void ProcessSpot (const EnergySpot &spot, double weight) override final
 ProcessHitsMethod.
void Initialize (G4HCofThisEvent *) override
 Register this regular SD with the event-owned builder.
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 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 Member Functions

G4int getTimeBin (G4double time) const
 Helper method for time info.
LArHitContainerBuilder * getHitContainer () const

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.
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, setupHelpers)
 FRIEND_TEST (LArG4SimpleSDtest, getTimeBin)
 FRIEND_TEST (LArG4SimpleSDtest, SimpleHit)
 FRIEND_TEST (LArG4SimpleSDtest, ConvertID)

Private Attributes

const EMBDetectorManager * m_embManager
const EMBAccordionDetails * m_accordionDetails
const GeoStraightAccSection * m_absorberSections
LArG4Identifier m_larID
 My LAr identifier.
std::string m_hitCollectionName
std::string m_hitSourceName

Detailed Description

NEEDS DOCUMENTATION.

Definition at line 25 of file BarrelFastSimDedicatedSD.h.

Member Enumeration Documentation

◆ LArHitTimeBins

enum LArG4SimpleSD::LArHitTimeBins
inherited
Enumerator
HitTimeBinDefault 
HitTimeBinUniform 

Definition at line 46 of file LArG4SimpleSD.h.

47 {
48 HitTimeBinDefault = 0,
49 HitTimeBinUniform = 1
50 };

Constructor & Destructor Documentation

◆ BarrelFastSimDedicatedSD()

BarrelFastSimDedicatedSD::BarrelFastSimDedicatedSD ( StoreGateSvc * detStore,
std::string hitCollectionName,
bool verbose )

Definition at line 25 of file BarrelFastSimDedicatedSD.cxx.

28 : IFastSimDedicatedSD("BarrelFastSimDedicatedSD", detStore,
29 std::move(hitCollectionName))
30 , m_embManager(nullptr)
31 , m_accordionDetails(nullptr)
32 , m_absorberSections(nullptr)
33{
34 if (verbose) { G4cout << GetName() << "::initialize()" << G4endl; }
35 if ( detStore->retrieve( m_embManager ).isFailure() ){
36 throw std::runtime_error("Could not retrieve EMB manager");
37 }
38 m_accordionDetails=m_embManager->getAccordionDetails();
39 m_absorberSections=m_accordionDetails->getAbsorberSections();
40}
BarrelFastSimDedicatedSD::m_embManager
const EMBDetectorManager * m_embManager
Definition BarrelFastSimDedicatedSD.h:42
BarrelFastSimDedicatedSD::m_absorberSections
const GeoStraightAccSection * m_absorberSections
Definition BarrelFastSimDedicatedSD.h:44
BarrelFastSimDedicatedSD::m_accordionDetails
const EMBAccordionDetails * m_accordionDetails
Definition BarrelFastSimDedicatedSD.h:43
IFastSimDedicatedSD::IFastSimDedicatedSD
IFastSimDedicatedSD(const std::string &name, StoreGateSvc *detStore, std::string hitCollectionName)
Simple constructor and destructor.
Definition IFastSimDedicatedSD.h:21
StoreGateSvc::retrieve
StatusCode retrieve(const T *&ptr) const
Retrieve the default object into a const T*.
verbose
bool verbose
Definition hcg.cxx:75

◆ ~BarrelFastSimDedicatedSD()

BarrelFastSimDedicatedSD::~BarrelFastSimDedicatedSD ( )
inline

Definition at line 35 of file BarrelFastSimDedicatedSD.h.

35{}

Member Function Documentation

◆ ConvertID()

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

Helper function for making "real" identifiers from LArG4Identifiers.

◆ FRIEND_TEST() [1/5]

LArG4SimpleSD::FRIEND_TEST ( LArG4SimpleSDtest ,
ConvertID  )
privateinherited

◆ FRIEND_TEST() [2/5]

LArG4SimpleSD::FRIEND_TEST ( LArG4SimpleSDtest ,
getTimeBin  )
privateinherited

◆ FRIEND_TEST() [3/5]

LArG4SimpleSD::FRIEND_TEST ( LArG4SimpleSDtest ,
ProcessHits  )
privateinherited

◆ FRIEND_TEST() [4/5]

LArG4SimpleSD::FRIEND_TEST ( LArG4SimpleSDtest ,
setupHelpers  )
privateinherited

◆ FRIEND_TEST() [5/5]

LArG4SimpleSD::FRIEND_TEST ( LArG4SimpleSDtest ,
SimpleHit  )
privateinherited

◆ getHitContainer()

LArHitContainerBuilder * LArG4SimpleSD::getHitContainer ( ) const
protectedinherited

◆ getTimeBin()

G4int LArG4SimpleSD::getTimeBin ( G4double time) const
protectedinherited

Helper method for time info.

◆ Initialize()

void LArG4SimpleSD::Initialize ( G4HCofThisEvent * )
overrideinherited

Register this regular SD with the event-owned builder.

◆ 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 43 of file BarrelFastSimDedicatedSD.cxx.

43 {
44
45 // Fill the identifier.
46 Point3D<double> globalPosition=spot.GetPosition();
47
48 static const Transform3D xfNeg = Amg::EigenTransformToCLHEP(m_embManager->getDetectorRegion(0,0,0)->getAbsoluteTransform().inverse());
49 static const Transform3D xfPos = Amg::EigenTransformToCLHEP(m_embManager->getDetectorRegion(1,0,0)->getAbsoluteTransform().inverse());
50
51 Point3D<double> localPosition = globalPosition.z()<0 ? xfNeg*globalPosition : xfPos*globalPosition;
52 int zIndex = globalPosition.z() <0 ? 0:1;
53 double eta = localPosition.getEta();
54 double phi = localPosition.getPhi();
55 double r = localPosition.perp();
56
57 bool implementWaves=true;
58
59 EMBDetectorManager::DetectorRegionConstIterator e;
60 for (e=m_embManager->beginDetectorRegion();e!=m_embManager->endDetectorRegion(); ++e) {
61 const EMBDetectorRegion *region = *e;
62 if (region->getEndcapIndex()!=zIndex) continue;
63 const EMBDetDescr *regionDescriptor=region->getDescriptor();
64 const CellBinning & etaBinning=regionDescriptor->getEtaBinning();
65 if (eta>etaBinning.getStart() && eta<etaBinning.getEnd()) {
66 unsigned int etaIndex = int((eta - etaBinning.getStart())/etaBinning.getDelta()) + etaBinning.getFirstDivisionNumber();
67 unsigned int phiIndex = 0;
68 EMBCellConstLink cellPtr = region->getEMBCell(etaIndex,phiIndex);
69 double rmax = cellPtr->getRLocal(EMBCell::BACK);
70 double rmin = cellPtr->getRLocal(EMBCell::FRONT);
71 if (r>rmin && r<rmax) {
72 const CellBinning & phiBinning=regionDescriptor->getPhiBinning();
73 double minPhi=std::min(phiBinning.getStart(), phiBinning.getEnd());
74 double maxPhi=std::max(phiBinning.getStart(), phiBinning.getEnd());
75
76 while (phi<minPhi) phi += 2*M_PI;
77 while (phi>maxPhi) phi -= 2*M_PI;
78 unsigned int regionIndex = region->getRegionIndex();
79 unsigned int samplingIndex = region->getSamplingIndex();
80 unsigned int phiIndex = int((phi - phiBinning.getStart())/phiBinning.getDelta()) + phiBinning.getFirstDivisionNumber();
81
82 // Now get the nominal cell and see if the indices need to be corrected for the accordion waves:
83 if (samplingIndex!=0) {
84 EMBCellConstLink cellPtr=region->getEMBCell(etaIndex,phiIndex);
85
86 double delta=-2.0*M_PI/1024/2;
87
88 double phiLocalUpper=cellPtr->getPhiLocalUpper()-delta;
89 double phiLocalLower=cellPtr->getPhiLocalLower()-delta;
90
91 while (phiLocalLower<0) phiLocalLower += 2.0*M_PI;
92 while (phiLocalUpper<0) phiLocalUpper += 2.0*M_PI;
93 while (phiLocalLower>2*M_PI) phiLocalLower -= 2.0*M_PI;
94 while (phiLocalUpper>2*M_PI) phiLocalUpper -= 2.0*M_PI;
95
96 int accordionIndexUpper=int(phiLocalUpper/(2*M_PI)*1024+0.25); if (accordionIndexUpper==1024) accordionIndexUpper=0;
97 int accordionIndexLower=int(phiLocalLower/(2*M_PI)*1024+0.25); if (accordionIndexLower==1024) accordionIndexLower=0;
98
99 Point3D<double> A0 =
100 ( Point3D<double>(m_absorberSections->XCent(accordionIndexLower,1),m_absorberSections->YCent(accordionIndexLower,1),0)+
101 Point3D<double>(m_absorberSections->XCent(accordionIndexUpper,1),m_absorberSections->YCent(accordionIndexUpper,1),0))/2.0;
102
103 Point3D<double> A1=
104 (Point3D<double>(m_absorberSections->XCent(accordionIndexLower,12),m_absorberSections->YCent(accordionIndexLower,12),0)+
105 Point3D<double>(m_absorberSections->XCent(accordionIndexUpper,12),m_absorberSections->YCent(accordionIndexUpper,12),0))/2.0-A0;
106
107 Point3D<double> P0=Point3D<double>(localPosition.x(), localPosition.y(),0)-A0;
108 int stackIndex=int(A1.dot(P0)/A1.mag2()*22.0 + 3.0)/2 ;
109
110 if (stackIndex<0 || stackIndex>13) {
111 G4cout << "Warning, bad stack index " << stackIndex << ' ' << rmin << ' ' << r << ' ' << A0.perp() << ' ' << A0 << ' ' << A1 << ' ' << P0 << G4endl;
112 if (implementWaves) return;
113 } else {
114 double xcent = m_absorberSections->XCent (accordionIndexLower,stackIndex);
115 double ycent = m_absorberSections->YCent (accordionIndexLower,stackIndex);
116 double cosU = m_absorberSections->Cosu (accordionIndexLower,stackIndex);
117 double sinU = m_absorberSections->Sinu (accordionIndexLower,stackIndex);
118 double halfLength = m_absorberSections->HalfLength(accordionIndexLower,stackIndex);
119
120 Point3D<double> u(xcent-halfLength*cosU, ycent-halfLength*sinU,0);
121 Point3D<double> v(xcent+halfLength*cosU, ycent+halfLength*sinU,0);
122 Point3D<double> x=v-u,y=Hep3Vector(localPosition.x(),localPosition.y(),0)-u;
123 if ((zIndex==0 && x.cross(y).z()>0) || (zIndex==1&& x.cross(y).z()<0)) {
124 if (implementWaves) {
125 if (phiIndex==0) phiIndex=phiBinning.getNumDivisions()-1;
126 else phiIndex--;
127 }
128 } else {
129 double xcent = m_absorberSections->XCent (accordionIndexUpper,stackIndex);
130 double ycent = m_absorberSections->YCent (accordionIndexUpper,stackIndex);
131 double cosU = m_absorberSections->Cosu (accordionIndexUpper,stackIndex);
132 double sinU = m_absorberSections->Sinu (accordionIndexUpper,stackIndex);
133 double halfLength = m_absorberSections->HalfLength(accordionIndexUpper,stackIndex);
134
135 Point3D<double> u(xcent-halfLength*cosU, ycent-halfLength*sinU,0);
136 Point3D<double> v(xcent+halfLength*cosU, ycent+halfLength*sinU,0);
137 Point3D<double> x=v-u,y=Hep3Vector(localPosition.x(),localPosition.y(),0)-u;
138 if ((zIndex==0 && x.cross(y).z()<0) || (zIndex==1 && x.cross(y).z()>0)) {
139 if (implementWaves) {
140 if (phiIndex==phiBinning.getNumDivisions()-1) phiIndex=0;
141 else phiIndex++;
142 }
143 }
144 }
145 }
146 }
147
148 //static LArG4Identifier id;
149 m_larID.clear();
150 m_larID << 4 // LArCalorimeter
151 << 1 // LArEM
152 << ((zIndex==0) ? -1:1)
153 << samplingIndex
154 << regionIndex
155 << etaIndex
156 << phiIndex;
157 // call process to add this to the collection
158 SimpleHit(m_larID, spot.GetTime(), spot.GetEnergy()*weight);
159
160 return;
161 }
162 }
163 }
164
165}
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:47
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 75 of file LArG4SimpleSD.h.

77 {
78 m_larEmID = EM;
79 m_larFcalID = FCAL;
80 m_larHecID = HEC;
81 }
LArG4SimpleSD::m_larHecID
const LArHEC_ID * m_larHecID
Definition LArG4SimpleSD.h:117
LArG4SimpleSD::m_larFcalID
const LArFCAL_ID * m_larFcalID
Definition LArG4SimpleSD.h:116
LArG4SimpleSD::m_larEmID
const LArEM_ID * m_larEmID
Pointers to the identifier helpers.
Definition LArG4SimpleSD.h:115

◆ 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 44 of file BarrelFastSimDedicatedSD.h.

◆ m_accordionDetails

const EMBAccordionDetails* BarrelFastSimDedicatedSD::m_accordionDetails
private

Definition at line 43 of file BarrelFastSimDedicatedSD.h.

◆ m_calculator

ILArCalculatorSvc* LArG4SimpleSD::m_calculator
protectedinherited

Member variable - the calculator we'll use.

Definition at line 91 of file LArG4SimpleSD.h.

◆ m_embManager

const EMBDetectorManager* BarrelFastSimDedicatedSD::m_embManager
private

Definition at line 42 of file BarrelFastSimDedicatedSD.h.

◆ m_hitCollectionName

std::string LArG4SimpleSD::m_hitCollectionName
privateinherited

Definition at line 122 of file LArG4SimpleSD.h.

◆ m_hitSourceName

std::string LArG4SimpleSD::m_hitSourceName
privateinherited

Definition at line 123 of file LArG4SimpleSD.h.

◆ m_larEmID

const LArEM_ID* LArG4SimpleSD::m_larEmID
protectedinherited

Pointers to the identifier helpers.

Definition at line 115 of file LArG4SimpleSD.h.

◆ m_larFcalID

const LArFCAL_ID* LArG4SimpleSD::m_larFcalID
protectedinherited

Definition at line 116 of file LArG4SimpleSD.h.

◆ m_larHecID

const LArHEC_ID* LArG4SimpleSD::m_larHecID
protectedinherited

Definition at line 117 of file LArG4SimpleSD.h.

◆ m_larID

LArG4Identifier BarrelFastSimDedicatedSD::m_larID
private

My LAr identifier.

Definition at line 47 of file BarrelFastSimDedicatedSD.h.

◆ m_numberInvalidHits

G4int LArG4SimpleSD::m_numberInvalidHits
protectedinherited

Count the number of invalid hits.

Definition at line 94 of file LArG4SimpleSD.h.

◆ m_timeBinType

LArG4SimpleSD::LArHitTimeBins LArG4SimpleSD::m_timeBinType
protectedinherited

Definition at line 109 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 112 of file LArG4SimpleSD.h.

The documentation for this class was generated from the following files:
Generated on for ATLAS Offline Software by doxygen 1.16.1