LArStraightAbsorbers Class Reference

#include <LArStraightAbsorbers.h>

Collaboration diagram for LArStraightAbsorbers:

Public Member Functions
double	XCentAbs (int stackid, int cellid) const
double	YCentAbs (int stackid, int cellid) const
double	HalfLength (int stackid, int cellid) const
double	Cosu (int stackid, int cellid) const
double	Sinu (int stackid, int cellid) const

Static Public Member Functions
static const LArStraightAbsorbers *	GetInstance (const std::string &strDetector="")

Private Member Functions
	LArStraightAbsorbers (const std::string &strDetector="")
void	initXYCentAbs (const PhysicalVolumeAccessor &theAbsorbers, int stackid, int cellid)
void	initHalfLength (const PhysicalVolumeAccessor &theAbsorbers, int stackid, int cellid)
double	SlantAbs (const PhysicalVolumeAccessor &theAbsorbers, int stackid, int cellid) const

Private Attributes
double	m_xcent [1024][14] {}
double	m_ycent [1024][14] {}
double	m_cosu [1024][14] {}
double	m_sinu [1024][14] {}
double	m_halflength [1024][14] {}
int	m_parity

Detailed Description

Definition at line 11 of file LArStraightAbsorbers.h.

Constructor & Destructor Documentation

LArStraightAbsorbers()

LArStraightAbsorbers::LArStraightAbsorbers ( const std::string &  strDetector = "" )

private

Definition at line 17 of file LArStraightAbsorbers.cxx.

{
  const LArVG4DetectorParameters* parameters = LArVG4DetectorParameters::GetInstance();
  if (parameters->GetValue("LArEMBPhiAtCurvature",0)>0.)  m_parity=0;  // first wave goes up
  else                                                    m_parity=1;  // first wave goes down
 
  static const std::string prefix = strDetector.empty() ? "" : strDetector+"::";
  static const PhysicalVolumeAccessor pva(prefix+"LAr::EMB::STAC",
                                          prefix+"LAr::EMB::ThinAbs::Straight");
 
  for (int stackid=0; stackid<14; stackid++) {
    for (int cellid=0; cellid<1024; cellid++) {
      initXYCentAbs(pva, stackid,cellid);
      initHalfLength(pva, stackid,cellid);
      const double slant = SlantAbs(pva, stackid,cellid);
      sincos(slant, &m_sinu[cellid][stackid], &m_cosu[cellid][stackid]);
    }
  }
}

Member Function Documentation

Cosu()

double LArStraightAbsorbers::Cosu ( int  stackid, int  cellid  ) const

inline

Definition at line 28 of file LArStraightAbsorbers.h.

{ return m_cosu[cellid][stackid]; }

GetInstance()

const LArStraightAbsorbers * LArStraightAbsorbers::GetInstance ( const std::string &  strDetector = "" )

static

Definition at line 11 of file LArStraightAbsorbers.cxx.

{
  static const LArStraightAbsorbers instance(strDetector);
  return &instance;
}

HalfLength()

double LArStraightAbsorbers::HalfLength ( int  stackid, int  cellid  ) const

inline

Definition at line 27 of file LArStraightAbsorbers.h.

{ return m_halflength[cellid][stackid]; }

initHalfLength()

void LArStraightAbsorbers::initHalfLength ( const PhysicalVolumeAccessor &  theAbsorbers, int  stackid, int  cellid  )

private

Definition at line 78 of file LArStraightAbsorbers.cxx.

{
  double l = 0.;
  const int id=cellid+stackid*10000;
  const G4VPhysicalVolume *pv=theAbsorbers.GetPhysicalVolume(id);
  if (pv) {
    const G4LogicalVolume* lv = pv->GetLogicalVolume();
    const G4Trap* trap = (G4Trap*) lv->GetSolid();
    const G4VPhysicalVolume *pv2=theAbsorbers.GetPhysicalVolume(1000000+id);
    if (!pv2) l = trap->GetYHalfLength1();
    else {
      const G4LogicalVolume* lv2 = pv2->GetLogicalVolume();
      const G4Trap* trap2 = (G4Trap*) lv2->GetSolid();
      l = trap->GetYHalfLength1()+trap2->GetYHalfLength1();
    }
  }
  m_halflength[cellid][stackid] = l;
}

initXYCentAbs()

void LArStraightAbsorbers::initXYCentAbs ( const PhysicalVolumeAccessor &  theAbsorbers, int  stackid, int  cellid  )

private

Definition at line 37 of file LArStraightAbsorbers.cxx.

{
  const int id=cellid+stackid*10000;
  const G4VPhysicalVolume *pv=theAbsorbers.GetPhysicalVolume(id);
  if (!pv) {
    m_xcent[cellid][stackid] = 0;
    m_ycent[cellid][stackid] = 0;
    return;
  }
  const G4ThreeVector& tv=pv->GetTranslation();
  const G4VPhysicalVolume *pv2=theAbsorbers.GetPhysicalVolume(1000000+id);
  if (!pv2) {
    m_xcent[cellid][stackid] = tv.x();
    m_ycent[cellid][stackid] = tv.y();
  }
  else {
    const G4ThreeVector& tv2=pv2->GetTranslation();
    const G4LogicalVolume* lv = pv->GetLogicalVolume();
    const G4Trap* trap = (G4Trap*) lv->GetSolid();
    const G4LogicalVolume* lv2 = pv2->GetLogicalVolume();
    const G4Trap* trap2 = (G4Trap*) lv2->GetSolid();
    const double xl1=trap->GetYHalfLength1();
    const double xl2=trap2->GetYHalfLength1();
    m_xcent[cellid][stackid] = (tv.x()*xl1+tv2.x()*xl2)/(xl1+xl2);
    m_ycent[cellid][stackid] = (tv.y()*xl1+tv2.y()*xl2)/(xl1+xl2);
  }
}

Sinu()

double LArStraightAbsorbers::Sinu ( int  stackid, int  cellid  ) const

inline

Definition at line 29 of file LArStraightAbsorbers.h.

{ return m_sinu[cellid][stackid]; }

SlantAbs()

double LArStraightAbsorbers::SlantAbs ( const PhysicalVolumeAccessor &  theAbsorbers, int  stackid, int  cellid  ) const

private

Definition at line 65 of file LArStraightAbsorbers.cxx.

{
  // both stackid and cellid start from 0 in the following code
  const int id=cellid+stackid*10000;
  const G4VPhysicalVolume *pv=theAbsorbers.GetPhysicalVolume(id);
  if (!pv) return 0.;
  const G4RotationMatrix *rm=pv->GetRotation();
  double Slant = (stackid%2 ==m_parity) ? 180*CLHEP::deg-(rm->thetaY()):(rm->thetaY())-180*CLHEP::deg;
  if((stackid%2 == m_parity) && (rm->phiY() > 0)) Slant = 360.*CLHEP::deg - Slant;
  if((stackid%2 == (1-m_parity)) && (rm->phiY() < 0)) Slant = - Slant;
  return Slant;
}

XCentAbs()

double LArStraightAbsorbers::XCentAbs ( int  stackid, int  cellid  ) const

inline

Definition at line 25 of file LArStraightAbsorbers.h.

{ return m_xcent[cellid][stackid]; }

YCentAbs()

double LArStraightAbsorbers::YCentAbs ( int  stackid, int  cellid  ) const

inline

Definition at line 26 of file LArStraightAbsorbers.h.

{ return m_ycent[cellid][stackid]; }

Member Data Documentation

m_cosu

double LArStraightAbsorbers::m_cosu[1024][14] {}

private

Definition at line 19 of file LArStraightAbsorbers.h.

m_halflength

double LArStraightAbsorbers::m_halflength[1024][14] {}

private

Definition at line 21 of file LArStraightAbsorbers.h.

m_parity

int LArStraightAbsorbers::m_parity

private

Definition at line 22 of file LArStraightAbsorbers.h.

m_sinu

double LArStraightAbsorbers::m_sinu[1024][14] {}

private

Definition at line 20 of file LArStraightAbsorbers.h.

m_xcent

double LArStraightAbsorbers::m_xcent[1024][14] {}

private

Definition at line 17 of file LArStraightAbsorbers.h.

m_ycent

double LArStraightAbsorbers::m_ycent[1024][14] {}

private

Definition at line 18 of file LArStraightAbsorbers.h.

---

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

• LArStraightAbsorbers.h
• LArStraightAbsorbers.cxx