#include <LArBarrelGeometry.h>

Inheritance diagram for LArG4::Barrel::Geometry:

Collaboration diagram for LArG4::Barrel::Geometry:

Public Member Functions
	Geometry (const std::string &name, ISvcLocator *pSvcLocator)

virtual	~Geometry ()=default

virtual StatusCode	initialize () override final

virtual StatusCode	finalize () override final

virtual void	initializeForSDCreation () override final

virtual LArG4Identifier	CalculateIdentifier (const G4Step *) const override final

virtual void	findCell (CalcData &currentCellData, const double &x, const double &y, const double &z, const double &r, const double &eta, const double &phi, const bool detail) const override final

Protected Member Functions
	Geometry ()

Private Member Functions
LArG4Identifier	CalculateECAMIdentifier (const G4Step *, const G4int indEcam, const bool inSTAC=true, int zside=1) const

bool	CheckLArIdentifier (int sampling, int region, int eta, int phi) const

bool	CheckDMIdentifier (int type, int sampling, int region, int eta, int phi) const

double	Distance_Ele (const double &x, const double &y, const int &PhiC, int &Num_Straight, const int &Num_Coude, double &xl) const

double	Distance_Abs (const double &x, const double &y, const int &nabs, const int &Num_Straight, const int &Num_Coude) const

G4int	SampSeg (G4double, G4double, G4double, G4int &, G4int &, G4int &, G4int &, G4int &) const

void	GetRphi ()
	phi vs r of first absorber in nominal geometry More...

G4double	Phi0 (G4double) const

G4int	PhiGap (const double &, const double &, const double &) const

Private Attributes
Gaudi::Property< std::string >	m_detectorName {this, "DetectorName", "LArMgr"}

Gaudi::Property< bool >	m_testbeam {this, "TestBeam", false}

G4String	m_ecamName

double	m_rMinAccordion {0.}

double	m_rMaxAccordion {0.}

double	m_zMinBarrel {0.}

double	m_zMaxBarrel {0.}

double	m_zMaxBarrelDMMargin {0.}

double	m_etaMaxBarrel {0.}

int	m_NCellTot {0}

int	m_NCellMax {0}

int	m_Nbrt {0}

int	m_Nbrt1 {0}

double	m_gam0 {0.}

double	m_rint_eleFib {0.}

double *	m_rc {nullptr}

double *	m_phic {nullptr}

double *	m_xc {nullptr}

double *	m_yc {nullptr}

double *	m_delta {nullptr}

int	m_parity {0}

const LArCoudeElectrodes *	m_coudeelec {nullptr}

const LArCoudeAbsorbers *	m_coudeabs {nullptr}

const LArStraightElectrodes *	m_electrode {nullptr}

const LArStraightAbsorbers *	m_absorber {nullptr}

bool	m_iflSAG {false}

G4int	m_NRphi {0}

G4double	m_Rmin {0.}

G4double	m_Rmax {0.}

G4double	m_Rphi [5000] = {0}

G4double	m_dR {0.}

Detailed Description

Definition at line 36 of file LArBarrelGeometry.h.

Constructor & Destructor Documentation

◆ Geometry() [1/2]

LArG4::Barrel::Geometry::Geometry	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Definition at line 75 of file LArBarrelGeometry.cxx.

82 { m_parity=1; }

◆ ~Geometry()

virtual LArG4::Barrel::Geometry::~Geometry ( )

virtualdefault

◆ Geometry() [2/2]

LArG4::Barrel::Geometry::Geometry ( )

protected

Member Function Documentation

◆ CalculateECAMIdentifier()

LArG4Identifier LArG4::Barrel::Geometry::CalculateECAMIdentifier	(	const G4Step *	a_step,
		const G4int	indEcam,
		const bool	inSTAC = `true`,
		int	zside = `1`
	)		const

private

Definition at line 1035 of file LArBarrelGeometry.cxx.

                                                                                   : midinLocal.phi();
       const G4double radius2Zpos = xZpos*xZpos + yZpos*yZpos;
       const G4double radiusZpos = sqrt(radius2Zpos);
  
       CalcData currentCellData;
       if (m_testbeam) {
         currentCellData.zSide = 1;
       }
       else {
         currentCellData.zSide = zside;
       }
  
       // Check if the hit is in the fiducial range and in the STAC volume
       //  if yes this is active or inactive material
  
       if (inSTAC && radiusZpos >=m_rMinAccordion && radiusZpos <= m_rMaxAccordion &&
           zZpos <= m_zMaxBarrel && zZpos >= m_zMinBarrel && etaZpos <= m_etaMaxBarrel) {
  
 #ifdef  DEBUGHITS
         ATH_MSG_VERBOSE("This hit is in the STAC volume !!!!! ");
 #endif
  
         //   DETERMINATION of currentCellData.cellID,
         //   currentCellData.zSide, currentCellData.sampling,
         //   currentCellData.phiBin, currentCellData.etaBin,
         //   m_stackNumID
         const bool MapDetail(false);
         this->findCell( currentCellData, xZpos, yZpos, zZpos, radiusZpos, etaZpos, phiZpos, MapDetail );
  
         // adjust phi in the negative half barrel frame
  
         if( currentCellData.zSide == -1 )
           {
             if( currentCellData.sampling == 1 && currentCellData.region ==0 )
               {
                 currentCellData.phiBin = 31 - currentCellData.phiBin;
                 if(currentCellData.phiBin < 0 ) currentCellData.phiBin += 64;
               }
             if( currentCellData.sampling == 1 && currentCellData.region ==1 )
               {
                 currentCellData.phiBin = 127 - currentCellData.phiBin;
                 if(currentCellData.phiBin < 0 ) currentCellData.phiBin += 256;
               }
             if( currentCellData.sampling >= 2 )
               {
                 currentCellData.phiBin = 127 - currentCellData.phiBin;
                 if(currentCellData.phiBin < 0 ) currentCellData.phiBin += 256;
               }
           }
  
         // there are few hundred microns between the 4mm nominal beginning of the active region
         //  and the real beginning of the first strip at eta=0.025/8
         //  to avoid inactive energy with strip=0 assign this to strip=1
         if (currentCellData.sampling==1 && currentCellData.region==0 && currentCellData.etaBin==0) {
           currentCellData.etaBin=1;
         }
  
         result << 4          // LArCalorimeter
                << 1          // LArEM
                << currentCellData.zSide
                << currentCellData.sampling
                << currentCellData.region
                << currentCellData.etaBin
                << currentCellData.phiBin;
  
 #ifdef  DEBUGHITS
         ATH_MSG_VERBOSE("Here the identifier for the barrel ACTIVE ----> ");
         ATH_MSG_VERBOSE("eta in local frame --> " << etaZpos);
         ATH_MSG_VERBOSE("currentCellData.zSide  ----> " << currentCellData.zSide);
         ATH_MSG_VERBOSE("currentCellData.sampling  ----> " << currentCellData.sampling);
         ATH_MSG_VERBOSE("currentCellData.region  ----> " <<  currentCellData.region);
         ATH_MSG_VERBOSE("currentCellData.etaBin  ----> " << currentCellData.etaBin);
         ATH_MSG_VERBOSE("currentCellData.phiBin  ----> " << currentCellData.phiBin);
         ATH_MSG_VERBOSE("And also etafirst ----> " << thisStepPoint->GetPosition().pseudoRapidity());
 #endif
  
         //    if (!Geometry::CheckLArIdentifier(currentCellData.sampling,currentCellData.region,currentCellData.etaBin,currentCellData.phiBin)) {
         //      ATH_MSG_ERROR(" **  Bad LAr identifier " << currentCellData.sampling << " " << currentCellData.region << " "
         //                << currentCellData.etaBin << " " << currentCellData.phiBin);
         //      ATH_MSG_ERROR(" x,y,z,eta,phi " <<  xZpos << " " << yZpos << " " << zZpos
         //                << " " << radiusZpos << " " << etaZpos << " " << phiZpos);
         //    }
  
  
       }
       // hits in dead material part
       else {
  
         G4int sampling=0;
         G4int region=0;
         const G4int numDeadPhiBins = 64;
         double abs_eta = std::fabs(etaZpos);
         const double DM1EtaWidth = 0.1 ;
         const double DM1PhiWidth = 2.*M_PI / numDeadPhiBins ;
         currentCellData.etaBin = (G4int) ( abs_eta * (1./DM1EtaWidth) ) ;
         currentCellData.phiBin = (G4int) (phiZpos/ DM1PhiWidth );
         G4int type=1;
         // protect against rounding error for phi ~2pi
         if (currentCellData.phiBin==numDeadPhiBins) currentCellData.phiBin=currentCellData.phiBin-1;
  
         // adjust phi for negative half barrel
  
         if ( currentCellData.zSide == -1 ) {
           currentCellData.phiBin = 31 - currentCellData.phiBin;
           if (currentCellData.phiBin < 0 ) currentCellData.phiBin +=64 ;
         }
  
         // material in front of the active accordion
         if ( radiusZpos < m_rMinAccordion ) {
           sampling =1 ;
           region = 3 ;
           if (currentCellData.etaBin > 14) currentCellData.etaBin=14;
  
 #ifdef  DEBUGHITS
           ATH_MSG_VERBOSE("This hit is in the ECAM volume in front of the accordion (DEAD MATERIAL) !!!!! ");
 #endif
  
         } else if (radiusZpos >= m_rMaxAccordion){  // material behind the active accordion
           sampling = 2;
  
           if (abs_eta < 1.0 ) {
             region = 0 ;
 #ifdef  DEBUGHITS
             ATH_MSG_VERBOSE("This hit is in the ECAM volume behind accordion (DEAD MATERIAL 0)  !!!!! ");
 #endif
           } else if ( abs_eta >= 1.0 && abs_eta < 1.5) {
             region = 2;
             currentCellData.etaBin = currentCellData.etaBin - 10;    // to have etabin between 0 and 4
 #ifdef  DEBUGHITS
             ATH_MSG_VERBOSE("This hit is in the ECAM volume behind accordion (DEAD MATERIAL 2) !!!!! ");
 #endif
           } else {
             ATH_MSG_ERROR(" LArBarrelGeometry: hit behind accordion at eta>1.5 !!! ");
             region = 2;
             currentCellData.etaBin = 4;
           }
  
         } else if (zZpos <= m_zMinBarrel) {   // inactive matter between two EMB halves
           type=2;
           region=0;
           const G4int phisave=currentCellData.phiBin;
           const G4bool MapDetail(false);
           this->findCell( currentCellData, xZpos, yZpos, zZpos, radiusZpos, etaZpos, phiZpos, MapDetail );
           sampling = currentCellData.sampling; // sampling as in normal definition
           currentCellData.etaBin=0;
           currentCellData.phiBin=phisave;
  
         } else if (zZpos >= m_zMaxBarrel - m_zMaxBarrelDMMargin || abs_eta >= 1.40) { // inactive matter between EMB and scintillator including some margin for error
           if (abs_eta >1.0 && abs_eta < 1.5) {
             sampling=2;
             region=2;
             currentCellData.etaBin = currentCellData.etaBin - 10;    // to have etabin between 0 and 4
           } else if (abs_eta < 1.6) {
             sampling=1;
             region=4;
             currentCellData.etaBin=0;
           } else {
             ATH_MSG_ERROR(" LArBarrelGeometry: hit at eta>1.6 !!! ");
             sampling=1;
             region=4;
             currentCellData.etaBin=0;
           }
         } else {
           if (!m_testbeam) {
             const G4double thisStepEnergyDeposit = a_step->GetTotalEnergyDeposit() * a_step->GetTrack()->GetWeight();
             std::ostringstream dmLog;
             dmLog << "LArBarrelGeometry: cannot find region for DM hit..." << std::endl;
             dmLog << "LArBarrelGeometry: cannot find region for DM hit..." << std::endl;
             dmLog << "m_zMinBarrel: " << m_zMinBarrel << std::endl;
             dmLog << "m_zMaxBarrel: " << m_zMaxBarrel << std::endl;
             dmLog << "m_rMinAccordion: " << m_rMinAccordion << std::endl;
             dmLog << "m_rMaxAccordion: " << m_rMaxAccordion << std::endl;
             dmLog << "r,z,eta,phi " << radiusZpos << " " << zZpos << " " << etaZpos << " " << phiZpos << std::endl;
             dmLog << "x,y,z (Atlas) " << p.x() << " " << p.y() << " " << p.z() << std::endl;
             dmLog << " inSTAC " << inSTAC << std::endl;
             dmLog << " eDeposited " << thisStepEnergyDeposit << std::endl;
             const G4VPhysicalVolume* vol = thisStepPoint->GetPhysicalVolume();
             const G4String volName = vol->GetName();
             dmLog << " volName " << volName << std::endl;
             const G4int ndep = g4navigation->GetDepth();
             for (G4int ii=0;ii<=ndep;ii++) {
               const G4VPhysicalVolume* v1 = g4navigation->GetVolume(ii);
               const G4String vname = v1->GetName();
               dmLog << "vname " << vname << std::endl;
             }
             if (thisStepEnergyDeposit > 1.*CLHEP::MeV) {
               ATH_MSG_ERROR(dmLog.str());
             } else {
               ATH_MSG_WARNING(dmLog.str());
             }
           }
           // in test beam case, we can get there for leakage on the side (in phi) of the module
           // in this case, we attribute an identifier like inactive material
           else
             {
               G4bool MapDetail=false;
               this->findCell( currentCellData, xZpos, yZpos, zZpos, radiusZpos, etaZpos, phiZpos, MapDetail );
               //         ATH_MSG_ERROR(" Lateral lakage r,eta,phi " << radiusZpos << " " << etaZpos << " "
               //                   << phiZpos << "  sampling/region/eta/phi " << currentCellData.sampling << " " <<
               //              currentCellData.region << " " << currentCellData.etaBin << " " << currentCellData.phiBin);
               // protect against small space between z=4m and real beginning of ieta=1 in strips
               if (currentCellData.sampling==1 && currentCellData.region==0 && currentCellData.etaBin==0) {
                 currentCellData.etaBin=1;
                 //            ATH_MSG_ERROR("S1R0 etabin 0 found  r,z,phi local " << radiusZpos << " "
                 //                   << " " << zZpos << " " << phiZpos);
               }
               result << 4          // LArCalorimeter
                      << 1          // LArEM
                      << currentCellData.zSide
                      << currentCellData.sampling
                      << currentCellData.region
                      << currentCellData.etaBin
                      << currentCellData.phiBin;
               return result;
             }
         }
  
         result << 10             // LArCalorimeter
                << currentCellData.zSide * 4     // LArBarrel
                << type
                << sampling
                << region
                << currentCellData.etaBin
                << currentCellData.phiBin;
  
 #ifdef  DEBUGHITS
         ATH_MSG_VERBOSE("Here the identifier for the barrel DEAD materials ---->");
         ATH_MSG_VERBOSE("Type     ----> " << type);
         ATH_MSG_VERBOSE("Sampling ----> " << sampling);
         ATH_MSG_VERBOSE("Region   ----> " << region);
         ATH_MSG_VERBOSE("zSide  ----> "   << currentCellData.zSide*4);
         ATH_MSG_VERBOSE("etaBin   ----> " << currentCellData.etaBin);
         ATH_MSG_VERBOSE("phiBin   ----> " << currentCellData.phiBin);
 #endif
  
         //    if (!Geometry::CheckDMIdentifier(type,sampling,region,currentCellData.etaBin,currentCellData.phiBin)) {
         //        ATH_MSG_ERROR(" **  Bad DM identifier " << type << " " << sampling << " " << region << " "
         //                  << currentCellData.etaBin << " " << currentCellData.phiBin);
         //        ATH_MSG_ERROR("x,y,z,r,eta,phi" << xZpos << " " << yZpos << " " << zZpos <<
         //         " " << radiusZpos << " " << etaZpos << " " << phiZpos);
         //    }
  
       }
  
       return result;
  
     }
  
     bool Geometry::CheckDMIdentifier(int type, int sampling, int region, int eta, int phi) const
     {
  
       if (type <1 || type > 2) return false;
       if (type==1) {
         if (sampling<1 || sampling>2) return false;
         if (sampling==1) {
           if (region!=3 && region !=4) return false;
           if (phi<0 || phi>63) return false;
           if (region==3) {
             if (eta<0 || eta>14) return false;
           }
           if (region==4) {
             if (eta !=0) return false;
           }
         }
         if (sampling==2) {
           if (region !=0 && region !=2) return false;
           if (phi<0 || phi>63) return false;
           if (region==0){
             if (eta<0 || eta>9) return false;
           }
           if (region==2) {
             if (eta<0 || eta>4) return false;
           }
         }
       }
       if (type==2) {
         if (sampling<1 || sampling >3) return false;
         if (region !=0) return false;
         if (eta!=0) return false;

◆ CalculateIdentifier()

LArG4Identifier LArG4::Barrel::Geometry::CalculateIdentifier ( const G4Step * a_step ) const

finaloverridevirtual

Definition at line 955 of file LArBarrelGeometry.cxx.

◆ CheckDMIdentifier()

bool LArG4::Barrel::Geometry::CheckDMIdentifier	(	int	type,
		int	sampling,
		int	region,
		int	eta,
		int	phi
	)		const

private

Definition at line 1320 of file LArBarrelGeometry.cxx.

     {
       if (sampling<0 || sampling >3) return false;
       if (sampling==0) {
         if (region!=0) return false;
         if (eta<0 || eta>60) return false;
         if (phi<0 || phi>63) return false;
       }
       if (sampling==1) {
         if (region<0 || region >1) return false;
         if (region==0) {
           if (eta<1 || eta>447) return false;
           if (phi<0 || phi>63) return false;
         }
         if (region==1) {
           if (eta<0 || eta>2) return false;
           if (phi<0 || phi>255) return false;
         }
       }
       if (sampling==2) {
         if (region<0 || region >1) return false;
         if (region==0) {
           if (eta<0 || eta>55) return false;
           if (phi<0 || phi>255) return false;
         }
         if (region==1) {
           if (eta!=0) return false;
           if (phi<0 || phi>255) return false;
         }
       }

◆ CheckLArIdentifier()

bool LArG4::Barrel::Geometry::CheckLArIdentifier	(	int	sampling,
		int	region,
		int	eta,
		int	phi
	)		const

private

Definition at line 1357 of file LArBarrelGeometry.cxx.

◆ Distance_Abs()

double LArG4::Barrel::Geometry::Distance_Abs	(	const double &	x,
		const double &	y,
		const int &	nabs,
		const int &	Num_Straight,
		const int &	Num_Coude
	)		const

private

Definition at line 259 of file LArBarrelGeometry.cxx.

                                { m_coudeabs->XCentCoude(Num_Coude, PhiCell), m_coudeabs->YCentCoude(Num_Coude, PhiCell) };
         // c_Hit Vector components and its length
         dx = xhit - Xc[0];
         dy = yhit - Xc[1];
         const double dr = sqrt( dx*dx + dy*dy);
         return (Num_Coude%2 == m_parity) ? (m_rint_eleFib-dr) : (dr - m_rint_eleFib);
       }      // end of Fold Regions
     } // end of the function Distance_Abs
  
  
     //=============================================================================
     //  Function SampSeg
     //
     // eta-sampling segmentation of barrel calorimeter   GU, January 2005
     //  input values: eta,radius in half-barrel frame
     //
     //  return value of function:  true=active area, false=inactive area
     //  return arguments: iregion,isampling,ieta
     //                             take into account detailed electrode drawing
     //                             with readout strips
     //                            isamp2,ieta2 do not take into account
     //                             readout strips and can be used to access current
     //                             maps.
     //
     //  iregion=0 (eta<1.4)  or 1 (eta=1.4-1.475)
     //  for region 0: isampling = 1 (strips), 2 (middle), 3 (back)
     //  for region 1: isampling = 1 or isampling = 2
     //  ieta= eta cell number
     //    region0,samp1:  ieta=1->448 (strip 0 does not exist)
     //    region0,samp2:  ieta=0->55
     //    region0,samp3:  ieta=0->26  (max eta 1.325)
     //    region1,samp1:  ieta=0->2 (deta=0.025)
     //    region1,samp2:  ieta=0 (only 1 cell)
  
     G4int Geometry::SampSeg(G4double eta, G4double radius, G4double z,
                             G4int& iregion, G4int& isampling, G4int& ieta,
                             G4int& isamp2, G4int& ieta2) const
     {
       // Helper struct to hold pre-calculated geometry information
       struct Geo {

◆ Distance_Ele()

double LArG4::Barrel::Geometry::Distance_Ele	(	const double &	x,
		const double &	y,
		const int &	PhiC,
		int &	Num_Straight,
		const int &	Num_Coude,
		double &	xl
	)		const

private

Definition at line 185 of file LArBarrelGeometry.cxx.

                            { m_electrode->XCentEle(Num_Straight, PhiCell), m_electrode->YCentEle(Num_Straight, PhiCell) };
       // m_Hit Vector components
       double dx = xhit - Xm[0];
       double dy = yhit - Xm[1];
  
       // First compute algebric distance m_hit (2D) the 2D_projection of the
       // m_Hit Vector on this electrode neutral fiber.
       const double hit = dx*u[0] + dy*u[1];
  
       //
       // Flat of Fold Region ?
       //
       const G4double Half_Elec(m_electrode->HalfLength(Num_Straight,PhiCell));
  
       if(std::fabs(hit) < Half_Elec) {
         // Flat Region
         xl=hit/Half_Elec;
         return u[0]*dy - u[1]*dx;
       }
       else {
         // Fold region
         // c_Hit Vector components and its length
         dx = xhit - Xc[0];
         dy = yhit - Xc[1];
         const double dr = sqrt( dx*dx + dy*dy);
         if (Num_Coude==Num_Straight) { xl=-1.; }
         else xl=+1;
         return (Num_Coude%2 == m_parity) ? (m_rint_eleFib-dr) : (dr - m_rint_eleFib);
       }      // end of Fold Regions
     } // end of the function Distance_Ele
  
  
     //======================================================================================
     //  Algebric distance to absorber
     //
     // inputs: xhit,yhit = x,y positions in local half barrel
     //         PhiCell      = absorber number in phi (0 to 1023 for Atlas case)
     //         Num_Straight = number (0 to 13) of the straight section
     //         Num_Coude    = number (0 to 14) of closest fold
     //
     // output:  Function value = algebric distance to electrode
  
     double Geometry::Distance_Abs(const double & xhit,
                                   const double &yhit, const int &PhiCell, const int &Num_Straight,
                                   const int &Num_Coude) const
     {
       //
       // FrameWork is consistent with the one used to PhiCell determination
       // e.g. it assumes HERE to be the LOCAL one of "stac_phys1",
       // (mother of ACCordion volumes) from which Z> 0. and Z < 0. half_barrel
       // parts are then defined.
       //
       // One needs POINTERS to Electrode neutral fibers
       // either for straight parts or for folds
       //
       // u unit 2D_Vector along straight part of the electrode neutral fiber
       const G4double u[2] = { m_absorber->Cosu(Num_Straight, PhiCell), m_absorber->Sinu(Num_Straight, PhiCell) };
       // Middle m_coordinates of this straight part of the electrode neutral fiber
       const G4double Xm[2] = { m_absorber->XCentAbs(Num_Straight, PhiCell), m_absorber->YCentAbs(Num_Straight, PhiCell) };
       // m_Hit Vector components
       double dx = xhit - Xm[0];  double dy = yhit - Xm[1];
  

◆ finalize()

StatusCode LArG4::Barrel::Geometry::finalize ( )

finaloverridevirtual

Definition at line 156 of file LArBarrelGeometry.cxx.

     {
       //
       // FrameWork is consistent with the one used to PhiCell determination
       // e.g. it assumes HERE to be the LOCAL one of "stac_phys1",
       // (mother of ACCordion volumes) from which Z> 0. and Z < 0. half_barrel
       // parts are then defined.
       //
       // One needs POINTERS to Electrode neutral fibers
       // either for straight parts or for folds
       //

◆ findCell()

void LArG4::Barrel::Geometry::findCell	(	CalcData &	currentCellData,
		const double &	x,
		const double &	y,
		const double &	z,
		const double &	r,
		const double &	eta,
		const double &	phi,
		const bool	detail
	)		const

finaloverridevirtual

Definition at line 612 of file LArBarrelGeometry.cxx.

                                                                                                           {
         currentCellData.nfold +=1;
       }
       if (currentCellData.nfold <0 || currentCellData.nfold >= m_Nbrt1) {
         ATH_MSG_ERROR("Invalid fold number " << currentCellData.nfold);
         return;
       }
  
  
 #ifdef DEBUGHITS
       ATH_MSG_VERBOSE("  BarrelGeometry: radius,eta,phi " << aRadius << " " << anEta << " ");
       ATH_MSG_VERBOSE("  Straight/Fold numbers " << currentCellData.nstraight << " " << currentCellData.nfold);
 #endif
  
       // eta and longitudinal segmentations
       G4int ireg,isamp,ieta,isamp2,ieta2;
       currentCellData.cellID = this->SampSeg(anEta,aRadius,zPosition,ireg,isamp,ieta,isamp2,ieta2);
  
       currentCellData.etaBin   = ieta;
       currentCellData.sampling = isamp;
       currentCellData.region   = ireg;
       currentCellData.etaMap   = ieta2;
       currentCellData.sampMap  = isamp2;
  
       // compute electrode number in phi
       int phicell = this->PhiGap(aRadius,xPosition,yPosition);
       if (phicell<0) phicell=0;
       // for test beam, some protection
       if (m_NCellTot !=1024) {
         if (phicell>=m_NCellTot) {
           if (phicell<512) { phicell=m_NCellTot-1; }
           else { phicell=0; }
           currentCellData.cellID=0;
         }
       }
  
 #ifdef DEBUGHITS
       ATH_MSG_VERBOSE(" phigap " << phicell);
 #endif
  
       // compute readout cell number
       int sampling_phi_nGaps=4;
       if (currentCellData.region==0 && currentCellData.sampling==1) { sampling_phi_nGaps=16; }
  
       if (currentCellData.cellID==0) {
         currentCellData.phiBin = (G4int) ( phicell/sampling_phi_nGaps );
         currentCellData.distElec=9999.;
         return;
       }
  
       // compute distance to electrode
       G4double xl;
       G4int nstr = currentCellData.nstraight;
       const G4double distElec = this->Distance_Ele(xPosition,yPosition,phicell,nstr,currentCellData.nfold,xl);
  
 #ifdef DEBUGHITS
       ATH_MSG_VERBOSE(" distElec " << distElec);
 #endif
  
       // if distance is < 2.5mm we are sure to be in the correct gap
       if (std::fabs(distElec) > 2.5) {
         // try +-2 electrode in phi to get minimum distance
         double dElecMin=distElec;
         double xlmin=xl;
         int phicellmin=phicell;
         for (int ii=-2;ii<3;ii++) {
           if (ii==0) { continue; }
           int phicellnew = phicell+ii;
           //  for test beam no phi wrapping
           if (m_NCellTot != 1024 && ( phicellnew<0 || phicellnew >= m_NCellTot)) { continue; }
           if (phicellnew < 0) { phicellnew += m_NCellTot; }
           if (phicellnew >= m_NCellTot) { phicellnew -= m_NCellTot; }
           double xln;
           int nstr2=currentCellData.nstraight;
           double dElec = Distance_Ele(xPosition,yPosition,phicellnew,nstr2,currentCellData.nfold,xln);
           if (std::fabs(dElec)<std::fabs(dElecMin)) {
             phicellmin=phicellnew;
             xlmin=xln;
             dElecMin = dElec;
             nstr=nstr2;
           }
         }
         currentCellData.phiGap = phicellmin;
         currentCellData.distElec = dElecMin;
         currentCellData.xl = xlmin;
         currentCellData.nstraight = nstr;
       }       // end distance >2.5mm
       else {
         currentCellData.phiGap=phicell;
         currentCellData.distElec=distElec;
         currentCellData.xl=xl;
         currentCellData.nstraight=nstr;
       }
  
 #ifdef DEBUGHITS
       ATH_MSG_VERBOSE("  final phiGap,distElec,xl " << currentCellData.phiGap << " " << currentCellData.distElec << " "
                       << currentCellData.xl);
 #endif
  
       // compute distance to absorber
  
       G4int nabs;
       if (currentCellData.distElec<0) nabs=currentCellData.phiGap;
       else nabs=currentCellData.phiGap+1;
       if (nabs >= m_NCellMax) nabs -= m_NCellMax;
       currentCellData.distAbs = Distance_Abs(xPosition,yPosition,nabs,currentCellData.nstraight,currentCellData.nfold);
 #ifdef DEBUGHITS
       ATH_MSG_VERBOSE("  nabs,distAbs " << nabs << " " << currentCellData.distAbs);
 #endif
  
       // in some rare cases near fold, the closest distance could give the wrong gap
       //  in such case, the signs of distAbs and distElec are not opposite as they should
       if ((currentCellData.distAbs>0. && currentCellData.distElec>0) ||
           (currentCellData.distAbs<0. && currentCellData.distElec<0) ) {
         //    ATH_MSG_VERBOSE("distElec and distAbs same sign " << currentCellData.distElec << " " << currentCellData.distAbs);
         //    ATH_MSG_VERBOSE(" currentCellData.phiGap " << currentCellData.phiGap);
         if (std::fabs(currentCellData.distElec)>std::fabs(currentCellData.distAbs)) {
           if (currentCellData.distAbs>0) { currentCellData.phiGap += 1; }
           if (currentCellData.distAbs<0) { currentCellData.phiGap -= 1; }
           if (m_NCellTot != 1024) {
             if (currentCellData.phiGap <0) { currentCellData.phiGap=0; }
             if (currentCellData.phiGap >= m_NCellTot) { currentCellData.phiGap = m_NCellTot-1; }
           }
           else {
             if (currentCellData.phiGap < 0) { currentCellData.phiGap += m_NCellTot; }
             if (currentCellData.phiGap >= m_NCellTot) { currentCellData.phiGap -= m_NCellTot; }
           }
           currentCellData.distElec = Distance_Ele(xPosition,yPosition,currentCellData.phiGap,currentCellData.nstraight,currentCellData.nfold,currentCellData.xl);
           //        ATH_MSG_VERBOSE(" new phiGap,distElec " << currentCellData.phiGap << " " << currentCellData.distElec);
         }
       }
  
       currentCellData.phiBin = (G4int) ( currentCellData.phiGap/sampling_phi_nGaps );
  
       if (MapDetail) {
         // compute x0,y0 coordinates in local electrode frame, using closest fold
         // as reference
         const G4double alpha = m_coudeelec->PhiRot(currentCellData.nfold,currentCellData.phiGap);
         const G4double dx=xPosition-m_coudeelec->XCentCoude(currentCellData.nfold,currentCellData.phiGap);
         const G4double dy=yPosition-m_coudeelec->YCentCoude(currentCellData.nfold,currentCellData.phiGap);
         const G4double dx1=dx*cos(alpha)-dy*sin(alpha);
         const G4double dy1=dx*sin(alpha)+dy*cos(alpha);
         currentCellData.x0 = dx1 + m_xc[currentCellData.nfold];
         currentCellData.y0 = dy1 + m_yc[currentCellData.nfold];
         if (m_parity==1) { currentCellData.y0 = -1*currentCellData.y0; }
       }
  
  
     } // end of findCell method
  
     // =============================================================================
     //  initialize phi0 vs radius of first absorber (for gam=0)
     void Geometry::GetRphi()
     {
       const G4double dl=0.001;
       const G4double inv_dl = 1. / dl;
       G4double cenx[15],ceny[15];
       //G4double xl,xl2;
       std::vector<G4double> sum1(5000);
       std::vector<G4double> sumx(5000);
       //xl=0;
       //xl2=0.;
       m_NRphi=5000;
       m_Rmin=1500.;
       m_dR=0.10;
       m_Rmax=0.;
  
       const G4double rint= m_rint_eleFib;
       const G4double inv_rint = 1. / rint;
       const G4double dt=dl * inv_rint;
       const G4double inv_dt = 1. / dt;
  
       for (G4int i=0;i<m_NRphi;i++) {
         sum1[i]=0.;
         sumx[i]=0.;
       }
       for (G4int i=0;i<15;i++) {
         cenx[i]=m_rc[i]*cos(m_phic[i]);
         ceny[i]=m_rc[i]*sin(m_phic[i]);
       }
  

◆ GetRphi()

void LArG4::Barrel::Geometry::GetRphi ( )

private

phi vs r of first absorber in nominal geometry

Initialize r-phi reference map (called from constructor)

Definition at line 796 of file LArBarrelGeometry.cxx.

                   {
           // first fold goes up
           if (m_parity==0) {
             phi0=-CLHEP::pi/2.;
             phi1=-m_delta[0];
           }
           // first fold goes down
           else {
             phi0=m_delta[0];
             phi1=CLHEP::pi/2;
           }
         }
         else if (i==14) {
           if (m_parity==0) {
             phi0=-CLHEP::pi+m_delta[13];
             phi1=-CLHEP::pi/2.;
           }
           else {
             phi0=CLHEP::pi/2;
             phi1=CLHEP::pi - m_delta[13];
           }
         }
         else {
           if (i%2==(1-m_parity)) {
             phi0=m_delta[i];
             phi1=CLHEP::pi-m_delta[i-1];
           }
           else {
             phi0=-CLHEP::pi+m_delta[i-1];
             phi1=-m_delta[i];
           }
         }
         //xl2+=rint*std::fabs(phi1-phi0);
         const G4int nstep=int((phi1-phi0)*inv_dt)+1;
         for (int ii=0;ii<nstep;ii++) {
           //xl+=dl;
           const G4double phi=phi0+dt*((G4double)ii);
           const G4double x=cenx[i]+rint*cos(phi);
           const G4double y=ceny[i]+rint*sin(phi);
           const G4double radius=sqrt(x*x+y*y);
           if (radius>m_Rmax) { m_Rmax=radius; }
           const G4double phid=atan(y/x);
           const G4int ir=((int) ((radius-m_Rmin)/m_dR) );
           if (ir>=0 && ir < m_NRphi) {
             sum1[ir]+=1.;
             sumx[ir]+=phid;
           }
         }
  
         // straight section
         if (i<14) {
           const G4double dx=cenx[i+1]-cenx[i];
           const G4double dy=ceny[i+1]-ceny[i];
           const G4double along=std::sqrt(dx*dx+dy*dy-4.*rint*rint);
           const G4double x0=0.5*(cenx[i+1]+cenx[i]);
           const G4double y0=0.5*(ceny[i+1]+ceny[i]);
           const G4double phi = (i%2==m_parity) ? CLHEP::pi/2-m_delta[i] : -CLHEP::pi/2.+m_delta[i];
           const G4double x1=x0-0.5*along*cos(phi);
           const G4double y1=y0-0.5*along*sin(phi);
           //xl2+=along;
           const int nstep=int(along*inv_dl)+1;
           for (int ii=0;ii<nstep;ii++) {
             //xl+=dl;
             const G4double x=x1+dl*((G4double)ii)*cos(phi);
             const G4double y=y1+dl*((G4double)ii)*sin(phi);
             const G4double radius=sqrt(x*x+y*y);
             if (radius>m_Rmax) { m_Rmax=radius; }
             const G4double phid=atan(y/x);
             const G4int ir=((int) ((radius-m_Rmin)/m_dR) );
             if (ir>=0 && ir < m_NRphi) {
               sum1[ir]+=1.;
               sumx[ir]+=phid;
             }
           }
         }
       }
       //  ATH_MSG_VERBOSE("total electrode length " << xl << " " << xl2);
       //  ATH_MSG_VERBOSE("rmax in accordion " << m_Rmax);
       for (int i=0; i<m_NRphi; i++) {
         if (sum1[i]>0) {
           m_Rphi[i]=sumx[i]/sum1[i];
           // Not used:
           //G4double radius = m_Rmin + ((G4double(i))+0.5)*m_dR;
           //ATH_MSG_VERBOSE(" GUTEST  r,phi0 " << radius << " " << m_Rphi[i]);
         }
         else { m_Rphi[i]=0.; }
       }
     }
  
     // ======================================================================================
     // phi of first absorber as function of radius for nominal accordion geometry
     //  (before sagging)
     G4double Geometry::Phi0(G4double radius) const
     {
       // TODO This function could be simplified.
       G4int ir;
       if (radius < m_Rmin) { ir=0; }
       else {
         if (radius > m_Rmax) radius=m_Rmax-0.0001;
         ir=((int) ((radius-m_Rmin)/m_dR) );
       }
       return m_Rphi[ir];
     }
  
     // ======================================================================================
     // compute number (0 to 1023) of closest electrode according to nominal
     // accordion geometry
     G4int Geometry::PhiGap(const double & radius, const double & xhit, const double &yhit) const
     {
       const G4double phi0=Phi0(radius)+m_gam0;   // from -pi to pi
       const G4double phi_hit=atan2(yhit,xhit);  // from -pi to pi
       G4double dphi=phi_hit-phi0;
       // bring back to 0-2pi
       if (dphi<0) dphi=dphi+2*M_PI;
       if (dphi>=2*M_PI) dphi=dphi-2*M_PI;
       dphi=dphi/(2*M_PI)*1024;
       const G4int ngap=((int) dphi);
 #ifdef DEBUGHITS
       ATH_MSG_VERBOSE(" phi0 " << phi0 << " dphi, ngap " << dphi << " " << ngap);
 #endif

◆ initialize()

StatusCode LArG4::Barrel::Geometry::initialize ( )

finaloverridevirtual

Definition at line 82 of file LArBarrelGeometry.cxx.

                         { m_parity=1; }
       //
       m_rMinAccordion  =   parameters->GetValue("LArEMBRadiusInnerAccordion");
       m_rMaxAccordion  =   parameters->GetValue("LArEMBFiducialRmax");
       m_etaMaxBarrel   =   parameters->GetValue("LArEMBMaxEtaAcceptance");
       m_zMinBarrel     =   parameters->GetValue("LArEMBfiducialMothZmin");
       m_zMaxBarrel     =   parameters->GetValue("LArEMBfiducialMothZmax");
       m_zMaxBarrelDMMargin = 10.0; // 10 mm margin
       // === GU 11/06/2003   total number of cells in phi
       // to distinguish 1 module (testbeam case) from full Atlas
       m_NCellTot = (int) (parameters->GetValue("LArEMBnoOFPhysPhiCell"));
       // total number of cells in phi to distinguish 1 module (testbeam case) from full Atlas
       m_testbeam=false;
       if (m_NCellTot != 1024) {
         m_testbeam=true;
       }
       m_NCellMax=1024;
       // ===
  
       // Initialize r-phi reference map
       this->GetRphi();
  
       if (m_detectorName.empty()) m_ecamName  = "LAr::EMB::ECAM";
       else                        m_ecamName  = m_detectorName+"::LAr::EMB::ECAM";
  
  
       return StatusCode::SUCCESS;
     }
  
     // ====================================================================================
  
     void Geometry::initializeForSDCreation()
     {
       // get pointers to access G4 geometry
       m_electrode = LArStraightElectrodes::GetInstance(m_detectorName);
       m_absorber  = LArStraightAbsorbers::GetInstance(m_detectorName);
       m_coudeelec = LArCoudeElectrodes::GetInstance(m_detectorName);
       m_coudeabs  = LArCoudeAbsorbers::GetInstance(m_detectorName);
     }
  
     // ====================================================================================
  
     StatusCode Geometry::finalize()
     {
       if (m_rc) delete [] m_rc;
       if (m_phic) delete [] m_phic;
       if (m_delta) delete [] m_delta;
       if (m_xc) delete [] m_xc;
       if (m_yc) delete [] m_yc;
  
       return StatusCode::SUCCESS;
     }
  
     //======================================================================================
     //
     // Here INTRINSIC Distance_to_electrode determination (G.P.)
     //
     // This retuns an ALGEBRICDistEle value, the distance from electrode
     //neutral fiber TOWARDS the Sub_Step in LAr (measured on a local perpendicular
     //vector unit oriented upwards i.e. following increasing Phi values).

◆ initializeForSDCreation()

void LArG4::Barrel::Geometry::initializeForSDCreation ( )

finaloverridevirtual

Definition at line 145 of file LArBarrelGeometry.cxx.

◆ Phi0()

G4double LArG4::Barrel::Geometry::Phi0 ( G4double radius ) const

private

Definition at line 921 of file LArBarrelGeometry.cxx.

     {
  
       // The default result is a blank identifier.
       LArG4Identifier result = LArG4Identifier();
  
       // Get all the required information from the current step
       const G4NavigationHistory* g4navigation = a_step->GetPreStepPoint()->GetTouchable()->GetHistory();
       const G4int ndep = g4navigation->GetDepth();

◆ PhiGap()

G4int LArG4::Barrel::Geometry::PhiGap	(	const double &	radius,
		const double &	xhit,
		const double &	yhit
	)		const

private

Definition at line 936 of file LArBarrelGeometry.cxx.

                                      {
         const G4String& vname = g4navigation->GetVolume(ii)->GetName();
         // FIXME Need to find a way to avoid these string-comparisons
         if ( indECAM<0 && vname == m_ecamName ) indECAM=ii;
         if ( !inSTAC && vname.find("STAC") !=std::string::npos) inSTAC=true;
         if ( vname.find("NegPhysical") != std::string::npos) zside=-1;
       }
       if (indECAM>=0)
         result = this->CalculateECAMIdentifier( a_step , indECAM, inSTAC, zside) ;
       else
         ATH_MSG_ERROR("LArBarrel::Geometry::CalculateIdentifier  ECAM volume not found in hierarchy");
  
       return result;
     }

◆ SampSeg()

G4int LArG4::Barrel::Geometry::SampSeg	(	G4double	eta,
		G4double	radius,
		G4double	z,
		G4int &	iregion,
		G4int &	isampling,
		G4int &	ieta,
		G4int &	isamp2,
		G4int &	ieta2
	)		const

private

Definition at line 326 of file LArBarrelGeometry.cxx.

                            {
  
         // get radius for end of strips
         istrip=(int) (aeta/deta*8.);
         if (istrip<0 || istrip >=448) {
           ATH_MSG_ERROR(" Problem aeta,istrip " << aeta << " " << istrip);
           return 0;
         }
         r12=Rmax1[istrip];
  
         // get radius for end of middle
         imid = (int) (aeta/deta);
         if (imid <0 || imid >=56) {
           ATH_MSG_ERROR(" Problem aeta,imid " << aeta << " " << imid);
           return 0;
         }
         r23=Rmax2[imid];
  
         iregion=0;
  
         // strips
         if (radius <= r12) {
           isampling=1;
           ieta=istrip;
           if (ieta==0) iactive=0;
           isamp2=1;
           ieta2=istrip;
         }
  
         // region between strips and middle => not active, same identifier as strips
         else if (radius < (r12+Dr_s12)) {
           isampling=1;
           ieta=istrip;
           iactive=0;
           isamp2=1;
           ieta2=istrip;
         }
  
         else {
  
           // eta<1.325, we can be in the back
           if (aeta<Eta_max_s3) {
             // radius<r23 we are in the middle
             if (radius <= r23) {
               isampling=2;
               ieta=imid;
               isamp2=2;
               ieta2=imid;
             }
             // for radius >r23 we have to take care of the readout strips at high z
             //   and attribute some of the energy to other cells
             else {     // radius>r23
               if (z>g.zmax1) {
                 isampling=2;
                 ieta=55;
               }
               else if (z>g.zmax2) {
                 isampling=2;
                 ieta=54;
               }
               else if (z>g.zmax3) {
                 isampling=2;
                 ieta=53;
               }
               else if (z>g.zmax4) {
                 isampling=3;
                 ieta=26;
               }
               else if (aeta<1.3 && z>g.zmax5) {
                 isampling=2;
                 ieta=52;
               }
               else if (aeta<1.3 && z>g.zmax6) {
                 isampling=2;
                 ieta=51;
               }
               else if (radius>g.rmax4 && z<g.zmax5 && aeta>1.2) {
                 if (radius>g.rmax1) {
                   isampling=2;
                   ieta=51;
                 }
                 else if(radius>g.rmax2) {
                   isampling=3;
                   ieta=25;
                 }
                 else if (radius>g.rmax3) {
                   if (z<g.zmax7) {
                     isampling=2;
                     ieta=50;
                   }
                   else {
                     isampling=3;
                     ieta=25;
                   }
                 }
                 else {
                   if (aeta<1.25) {
                     isampling=2;
                     ieta=49;
                   }
                   else {
                     isampling=3;
                     ieta=25;
                   }
                 }
               }
               // normal back cell
               else {
                 isampling=3;
                 ieta=imid/2;
                 isamp2=3;
                 ieta2=ieta;
               }
               isamp2=3;
               ieta2=imid/2;
             }   // end radius>r23
             // put into middle energy deposited along readout strips across the back
             if (isampling==3 && z<g.zmax4 && (radius<g.rmax4 || aeta<1.2) ) {
               const double etastr = (imid%2==0) ? 0.025*imid : 0.025*(imid+1);
               const double delta=radius*(sinh(etastr)-sinh(aeta))/cosh(etastr);
               double deltastr;
               if (aeta<0.475) { deltastr=1.5;}
               else if (aeta<0.80) { deltastr=2.75;}
               else if (aeta<0.85) { deltastr=1.5;}
               else if (aeta<1.1)  { deltastr=2.75;}
               else { deltastr=3.25;}
  
               if (std::fabs(delta)<deltastr) {
                 isampling=2;
                 ieta=imid;
               }
             }     // end if sampling==3
           }     // end if eta<1.325
           else {
             isampling=2;
             ieta=imid;
             if (z>g.zmax1) {
               ieta=55;
             }
             else if (z>g.zmax2 && aeta<1.375) {
               ieta=54;
             }
             isamp2=2;
             ieta2=imid;
           }       // end eta>1.352
         }          // end radius selection
       }     // end eta1.4
  
       // eta between 1.4 and 1.475
  
       if (aeta>=Eta_max_s1 && aeta<g.Eta_max) {
         r12 = Rmax1[447];   // radius for end of sampling 1
         r23=g.Z_max_acc/sinh(aeta);   // radius of end of sampling 2, bounded by high z end
  
         const double zmax = g.Z_max_lowr + g.dzdr*(radius-g.R_min_acc);
  
         iregion=1;
         if (radius <=r12) {
           isampling=1;
           ieta=int((aeta-Eta_max_s1)/deta);
           if (z>zmax) { iactive=0; }
         }
         else if (radius < (r12+Dr_s12)) {
           isampling=1;
           ieta=int((aeta-Eta_max_s1)/deta);
           iactive=0;
         }
         else if (radius <= r23) {
           isampling=2;
           ieta=0;
           if (z>zmax) { iactive=0; }
         }
         else {
           isampling=2;
           ieta=0;
           iactive=0;
         }
         isamp2=isampling;
         ieta2=ieta;
       }
       // eta above 1.475, not fiducial region, but still returns something
       //  for calibration hits
       if (aeta>g.Eta_max) {
         iregion=1;
         r12 = Rmax1[447];
         if (radius <=r12) {
           isampling=1;
           ieta=2;
         }
         else {
           isampling=2;
           ieta=0;
         }
         isamp2=isampling;
         ieta2=ieta;
         iactive=0;
       }
  
       // cross-check of active region
       if (z>g.Z_max_acc || radius>g.R_max_acc || radius<g.R_min_acc || aeta > g.Eta_max) iactive=0;
  
       return iactive;
     }
     // =======================================================================
     //   function findCell
     //
     // compute cell in EM accordion for hit at position x,y,z,radius,eta,phi
     //  given in LOCAL half barrel coordinate system (Stac Geant volume)
     // It has already been checked that the hit is in the accordion sensitive volume
     //
  
     void Geometry::findCell(CalcData & currentCellData,
                             const double &xPosition,
                             const double &yPosition,
                             const double &zPosition,
                             const double &aRadius,
                             const double &anEta,
                             const double &/*aPhi*/,
                             const bool MapDetail) const
     {
  
       currentCellData.cellID = 0;
  
       if (aRadius < m_rc[0] || aRadius >= m_rc[m_Nbrt1-1]) {
 #ifdef DEBUGHITS
         ATH_MSG_VERBOSE(" Outside Accordion " << aRadius << " " << m_rc[0] << " " << m_rc[m_Nbrt1-1]);
 #endif
         return;    // outside accordion
       }
  
       // set the straight section number
       currentCellData.nstraight=0;
       for (int i=1;i<m_Nbrt1;i++) {
         if (m_rc[i] > aRadius) { break; }
         currentCellData.nstraight++;

Member Data Documentation

◆ m_absorber

const LArStraightAbsorbers* LArG4::Barrel::Geometry::m_absorber {nullptr}

private

Definition at line 102 of file LArBarrelGeometry.h.

◆ m_coudeabs

const LArCoudeAbsorbers* LArG4::Barrel::Geometry::m_coudeabs {nullptr}

private

Definition at line 100 of file LArBarrelGeometry.h.

◆ m_coudeelec

const LArCoudeElectrodes* LArG4::Barrel::Geometry::m_coudeelec {nullptr}

private

Definition at line 99 of file LArBarrelGeometry.h.

◆ m_delta

double* LArG4::Barrel::Geometry::m_delta {nullptr}

private

Definition at line 95 of file LArBarrelGeometry.h.

◆ m_detectorName

Gaudi::Property<std::string> LArG4::Barrel::Geometry::m_detectorName {this, "DetectorName", "LArMgr"}

private

Definition at line 64 of file LArBarrelGeometry.h.

◆ m_dR

G4double LArG4::Barrel::Geometry::m_dR {0.}

private

Definition at line 111 of file LArBarrelGeometry.h.

◆ m_ecamName

G4String LArG4::Barrel::Geometry::m_ecamName

private

Definition at line 69 of file LArBarrelGeometry.h.

◆ m_electrode

const LArStraightElectrodes* LArG4::Barrel::Geometry::m_electrode {nullptr}

private

Definition at line 101 of file LArBarrelGeometry.h.

◆ m_etaMaxBarrel

double LArG4::Barrel::Geometry::m_etaMaxBarrel {0.}

private

Definition at line 77 of file LArBarrelGeometry.h.

◆ m_gam0

double LArG4::Barrel::Geometry::m_gam0 {0.}

private

Definition at line 88 of file LArBarrelGeometry.h.

◆ m_iflSAG

bool LArG4::Barrel::Geometry::m_iflSAG {false}

private

Definition at line 104 of file LArBarrelGeometry.h.

◆ m_Nbrt

int LArG4::Barrel::Geometry::m_Nbrt {0}

private

Definition at line 84 of file LArBarrelGeometry.h.

◆ m_Nbrt1

int LArG4::Barrel::Geometry::m_Nbrt1 {0}

private

Definition at line 85 of file LArBarrelGeometry.h.

◆ m_NCellMax

int LArG4::Barrel::Geometry::m_NCellMax {0}

private

Definition at line 81 of file LArBarrelGeometry.h.

◆ m_NCellTot

int LArG4::Barrel::Geometry::m_NCellTot {0}

private

Definition at line 80 of file LArBarrelGeometry.h.

◆ m_NRphi

G4int LArG4::Barrel::Geometry::m_NRphi {0}

private

Definition at line 107 of file LArBarrelGeometry.h.

◆ m_parity

int LArG4::Barrel::Geometry::m_parity {0}

private

Definition at line 96 of file LArBarrelGeometry.h.

◆ m_phic

double* LArG4::Barrel::Geometry::m_phic {nullptr}

private

Definition at line 92 of file LArBarrelGeometry.h.

◆ m_rc

double* LArG4::Barrel::Geometry::m_rc {nullptr}

private

Definition at line 91 of file LArBarrelGeometry.h.

◆ m_rint_eleFib

double LArG4::Barrel::Geometry::m_rint_eleFib {0.}

private

Definition at line 89 of file LArBarrelGeometry.h.

◆ m_Rmax

G4double LArG4::Barrel::Geometry::m_Rmax {0.}

private

Definition at line 109 of file LArBarrelGeometry.h.

◆ m_rMaxAccordion

double LArG4::Barrel::Geometry::m_rMaxAccordion {0.}

private

Definition at line 73 of file LArBarrelGeometry.h.

◆ m_Rmin

G4double LArG4::Barrel::Geometry::m_Rmin {0.}

private

Definition at line 108 of file LArBarrelGeometry.h.

◆ m_rMinAccordion

double LArG4::Barrel::Geometry::m_rMinAccordion {0.}

private

Definition at line 72 of file LArBarrelGeometry.h.

◆ m_Rphi

G4double LArG4::Barrel::Geometry::m_Rphi[5000] = {0}

private

Definition at line 110 of file LArBarrelGeometry.h.

◆ m_testbeam

Gaudi::Property<bool> LArG4::Barrel::Geometry::m_testbeam {this, "TestBeam", false}

private

Definition at line 67 of file LArBarrelGeometry.h.

◆ m_xc

double* LArG4::Barrel::Geometry::m_xc {nullptr}

private

Definition at line 93 of file LArBarrelGeometry.h.

◆ m_yc

double* LArG4::Barrel::Geometry::m_yc {nullptr}

private

Definition at line 94 of file LArBarrelGeometry.h.

◆ m_zMaxBarrel

double LArG4::Barrel::Geometry::m_zMaxBarrel {0.}

private

Definition at line 75 of file LArBarrelGeometry.h.

◆ m_zMaxBarrelDMMargin

double LArG4::Barrel::Geometry::m_zMaxBarrelDMMargin {0.}

private

Definition at line 76 of file LArBarrelGeometry.h.

◆ m_zMinBarrel

double LArG4::Barrel::Geometry::m_zMinBarrel {0.}

private

Definition at line 74 of file LArBarrelGeometry.h.

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

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Geometry() [1/2]

◆ ~Geometry()

◆ Geometry() [2/2]

Member Function Documentation

◆ CalculateECAMIdentifier()

◆ CalculateIdentifier()

◆ CheckDMIdentifier()

◆ CheckLArIdentifier()

◆ Distance_Abs()

◆ Distance_Ele()

◆ finalize()

◆ findCell()

◆ GetRphi()

◆ initialize()

◆ initializeForSDCreation()

◆ Phi0()

◆ PhiGap()

◆ SampSeg()

Member Data Documentation

◆ m_absorber

◆ m_coudeabs

◆ m_coudeelec

◆ m_delta

◆ m_detectorName

◆ m_dR

◆ m_ecamName

◆ m_electrode

◆ m_etaMaxBarrel

◆ m_gam0

◆ m_iflSAG

◆ m_Nbrt

◆ m_Nbrt1

◆ m_NCellMax

◆ m_NCellTot

◆ m_NRphi

◆ m_parity

◆ m_phic

◆ m_rc

◆ m_rint_eleFib

◆ m_Rmax

◆ m_rMaxAccordion

◆ m_Rmin

◆ m_rMinAccordion

◆ m_Rphi

◆ m_testbeam

◆ m_xc

◆ m_yc

◆ m_zMaxBarrel

◆ m_zMaxBarrelDMMargin

◆ m_zMinBarrel