#include <Doxymodules_optical.h>

Inheritance diagram for LXeMainVolume:

Public Member Functions
	LXeMainVolume (G4RotationMatrix pRot, const G4ThreeVector &tlate, G4LogicalVolume pMotherLogical, G4bool pMany, G4int pCopyNo, LXeDetectorConstruction *c)

G4LogicalVolume *	GetLogPhotoCath ()

G4LogicalVolume *	GetLogScint ()

std::vector< G4ThreeVector >	GetPmtPositions ()

Private Member Functions
void	VisAttributes ()

void	SurfaceProperties ()

void	PlacePMTs (G4LogicalVolume pmt_Log, G4RotationMatrix rot, G4double &a, G4double &b, G4double da, G4double db, G4double amin, G4double bmin, G4int na, G4int nb, G4double &x, G4double &y, G4double &z, G4int &k)

void	CopyValues ()

Private Attributes
LXeDetectorConstruction *	fConstructor = nullptr

G4double	fScint_x = 0.

G4double	fScint_y = 0.

G4double	fScint_z = 0.

G4double	fD_mtl = 0.

G4int	fNx = 0

G4int	fNy = 0

G4int	fNz = 0

G4double	fOuterRadius_pmt = 0.

G4bool	fSphereOn = false

G4double	fRefl = 0.

G4Box *	fScint_box = nullptr

G4Box *	fHousing_box = nullptr

G4Tubs *	fPmt = nullptr

G4Tubs *	fPhotocath = nullptr

G4Sphere *	fSphere = nullptr

G4LogicalVolume *	fScint_log = nullptr

G4LogicalVolume *	fHousing_log = nullptr

G4LogicalVolume *	fPmt_log = nullptr

G4LogicalVolume *	fPhotocath_log = nullptr

G4LogicalVolume *	fSphere_log = nullptr

std::vector< G4ThreeVector >	fPmtPositions

Detailed Description

Definition at line 54 of file Doxymodules_optical.h.

Constructor & Destructor Documentation

◆ LXeMainVolume()

LXeMainVolume::LXeMainVolume	(	G4RotationMatrix *	pRot,
		const G4ThreeVector &	tlate,
		G4LogicalVolume *	pMotherLogical,
		G4bool	pMany,
		G4int	pCopyNo,
		LXeDetectorConstruction *	c
	)

Definition at line 49 of file LXeMainVolume.cc.

  : G4PVPlacement(pRot, tlate,
                  // Temp logical volume must be created here
                  new G4LogicalVolume(new G4Box("temp", 1, 1, 1),
                                    G4Material::GetMaterial("Vacuum"), "temp"),
                  "housing", pMotherLogical, pMany, pCopyNo)
  , fConstructor(c)
{
  CopyValues();
 
  G4double housing_x = fScint_x + 2. * fD_mtl;
  G4double housing_y = fScint_y + 2. * fD_mtl;
  G4double housing_z = fScint_z + 2. * fD_mtl;
 
  //*************************** housing and scintillator
  fScint_box =
    new G4Box("scint_box", fScint_x / 2., fScint_y / 2., fScint_z / 2.);
  fHousing_box =
    new G4Box("housing_box", housing_x / 2., housing_y / 2., housing_z / 2.);
 
  fScint_log   = new G4LogicalVolume(fScint_box, G4Material::GetMaterial("LXe"),
                                   "scint_log");
  fHousing_log = new G4LogicalVolume(
    fHousing_box, G4Material::GetMaterial("Al"), "housing_log");
 
  new G4PVPlacement(nullptr, G4ThreeVector(), fScint_log, "scintillator",
                    fHousing_log, false, 0);
 
  //*************** Miscellaneous sphere to demonstrate skin surfaces
  fSphere = new G4Sphere("sphere", 0., 2. * cm, 0. * deg, 360. * deg, 0. * deg,
                         360. * deg);
  fSphere_log =
    new G4LogicalVolume(fSphere, G4Material::GetMaterial("Al"), "sphere_log");
  if(fSphereOn)
    new G4PVPlacement(nullptr, G4ThreeVector(5. * cm, 5. * cm, 5. * cm),
                      fSphere_log, "sphere", fScint_log, false, 0);
 
  //****************** Build PMTs
  G4double innerRadius_pmt   = 0.;
  G4double height_pmt        = fD_mtl / 2.;
  G4double startAngle_pmt    = 0.;
  G4double spanningAngle_pmt = 360. * deg;
 
  fPmt = new G4Tubs("pmt_tube", innerRadius_pmt, fOuterRadius_pmt, height_pmt,
                    startAngle_pmt, spanningAngle_pmt);
 
  // the "photocathode" is a metal slab at the back of the glass that
  // is only a very rough approximation of the real thing since it only
  // absorbs or detects the photons based on the efficiency set below
  fPhotocath = new G4Tubs("photocath_tube", innerRadius_pmt, fOuterRadius_pmt,
                          height_pmt / 2., startAngle_pmt, spanningAngle_pmt);
 
  fPmt_log =
    new G4LogicalVolume(fPmt, G4Material::GetMaterial("Glass"), "pmt_log");
  fPhotocath_log = new G4LogicalVolume(
    fPhotocath, G4Material::GetMaterial("Al"), "photocath_log");
 
  new G4PVPlacement(nullptr, G4ThreeVector(0., 0., -height_pmt / 2.),
                    fPhotocath_log, "photocath", fPmt_log, false, 0);
 
  //***********Arrange pmts around the outside of housing**********
 
  G4double dx = fScint_x / fNx;
  G4double dy = fScint_y / fNy;
  G4double dz = fScint_z / fNz;
 
  G4double x, y, z;
  G4double xmin = -fScint_x / 2. - dx / 2.;
  G4double ymin = -fScint_y / 2. - dy / 2.;
  G4double zmin = -fScint_z / 2. - dz / 2.;
  G4int k       = 0;
 
  z = -fScint_z / 2. - height_pmt;  // front
  PlacePMTs(fPmt_log, nullptr, x, y, dx, dy, xmin, ymin, fNx, fNy, x, y, z, k);
 
  auto rm_z = new G4RotationMatrix();
  rm_z->rotateY(180. * deg);
  z = fScint_z / 2. + height_pmt;  // back
  PlacePMTs(fPmt_log, rm_z, x, y, dx, dy, xmin, ymin, fNx, fNy, x, y, z, k);
 
  auto rm_y1 = new G4RotationMatrix();
  rm_y1->rotateY(-90. * deg);
  x = -fScint_x / 2. - height_pmt;  // left
  PlacePMTs(fPmt_log, rm_y1, y, z, dy, dz, ymin, zmin, fNy, fNz, x, y, z, k);
 
  auto rm_y2 = new G4RotationMatrix();
  rm_y2->rotateY(90. * deg);
  x = fScint_x / 2. + height_pmt;  // right
  PlacePMTs(fPmt_log, rm_y2, y, z, dy, dz, ymin, zmin, fNy, fNz, x, y, z, k);
 
  auto rm_x1 = new G4RotationMatrix();
  rm_x1->rotateX(90. * deg);
  y = -fScint_y / 2. - height_pmt;  // bottom
  PlacePMTs(fPmt_log, rm_x1, x, z, dx, dz, xmin, zmin, fNx, fNz, x, y, z, k);
 
  auto rm_x2 = new G4RotationMatrix();
  rm_x2->rotateX(-90. * deg);
  y = fScint_y / 2. + height_pmt;  // top
  PlacePMTs(fPmt_log, rm_x2, x, z, dx, dz, xmin, zmin, fNx, fNz, x, y, z, k);
 
  VisAttributes();
  SurfaceProperties();
 
  SetLogicalVolume(fHousing_log);
}

Member Function Documentation

◆ GetLogPhotoCath()

G4LogicalVolume * LXeMainVolume::GetLogPhotoCath ( )

inline

Definition at line 49 of file LXeMainVolume.hh.

49{ return fPhotocath_log; }

◆ GetLogScint()

G4LogicalVolume * LXeMainVolume::GetLogScint ( )

inline

Definition at line 50 of file LXeMainVolume.hh.

50{ return fScint_log; }

◆ GetPmtPositions()

std::vector< G4ThreeVector > LXeMainVolume::GetPmtPositions ( )

inline

Definition at line 52 of file LXeMainVolume.hh.

52{ return fPmtPositions; }

LXeMainVolume::fPmtPositions

std::vector< G4ThreeVector > fPmtPositions

Definition LXeMainVolume.hh:95

◆ VisAttributes()

void LXeMainVolume::VisAttributes ( )

private

Definition at line 213 of file LXeMainVolume.cc.

{
  auto housing_va = new G4VisAttributes(G4Colour(0.8, 0.8, 0.8));
  fHousing_log->SetVisAttributes(housing_va);
 
  auto sphere_va = new G4VisAttributes();
  sphere_va->SetForceSolid(true);
  fSphere_log->SetVisAttributes(sphere_va);
}

◆ SurfaceProperties()

void LXeMainVolume::SurfaceProperties ( )

private

Definition at line 225 of file LXeMainVolume.cc.

{
  std::vector<G4double> ephoton = { 7.0 * eV, 7.14 * eV };
 
  //**Scintillator housing properties
  std::vector<G4double> reflectivity     = { fRefl, fRefl };
  std::vector<G4double> efficiency       = { 0.0, 0.0 };
  auto scintHsngPT = new G4MaterialPropertiesTable();
  scintHsngPT->AddProperty("REFLECTIVITY", ephoton, reflectivity);
  scintHsngPT->AddProperty("EFFICIENCY", ephoton, efficiency);
  auto OpScintHousingSurface =
    new G4OpticalSurface("HousingSurface", unified, polished, dielectric_metal);
  OpScintHousingSurface->SetMaterialPropertiesTable(scintHsngPT);
 
  //**Sphere surface properties
  std::vector<G4double> sphereReflectivity = { 1.0, 1.0 };
  std::vector<G4double> sphereEfficiency   = { 0.0, 0.0 };
  auto spherePT = new G4MaterialPropertiesTable();
  spherePT->AddProperty("REFLECTIVITY", ephoton, sphereReflectivity);
  spherePT->AddProperty("EFFICIENCY", ephoton, sphereEfficiency);
  auto OpSphereSurface =
    new G4OpticalSurface("SphereSurface", unified, polished, dielectric_metal);
  OpSphereSurface->SetMaterialPropertiesTable(spherePT);
 
  //**Photocathode surface properties
  std::vector<G4double> photocath_EFF     = { 1., 1. };
  std::vector<G4double> photocath_ReR     = { 1.92, 1.92 };
  std::vector<G4double> photocath_ImR     = { 1.69, 1.69 };
  auto photocath_mt = new G4MaterialPropertiesTable();
  photocath_mt->AddProperty("EFFICIENCY", ephoton, photocath_EFF);
  photocath_mt->AddProperty("REALRINDEX", ephoton, photocath_ReR);
  photocath_mt->AddProperty("IMAGINARYRINDEX", ephoton, photocath_ImR);
  auto photocath_opsurf = new G4OpticalSurface(
    "photocath_opsurf", glisur, polished, dielectric_metal);
  photocath_opsurf->SetMaterialPropertiesTable(photocath_mt);
 
  //**Create logical skin surfaces
  new G4LogicalSkinSurface("photocath_surf", fHousing_log,
                           OpScintHousingSurface);
  new G4LogicalSkinSurface("sphere_surface", fSphere_log, OpSphereSurface);
  new G4LogicalSkinSurface("photocath_surf", fPhotocath_log, photocath_opsurf);
}

◆ PlacePMTs()

void LXeMainVolume::PlacePMTs	(	G4LogicalVolume *	pmt_Log,
		G4RotationMatrix *	rot,
		G4double &	a,
		G4double &	b,
		G4double	da,
		G4double	db,
		G4double	amin,
		G4double	bmin,
		G4int	na,
		G4int	nb,
		G4double &	x,
		G4double &	y,
		G4double &	z,
		G4int &	k
	)

private

Definition at line 176 of file LXeMainVolume.cc.

{
  /*  PlacePMTs : a different way to parameterize placement that does not depend
   * on calculating the position from the copy number
   *
   *  pmt_log = logical volume for pmts to be placed
   *  rot = rotation matrix to apply
   *  a,b = coordinates to vary(ie. if varying in the xy plane then pass x,y)
   *  da,db = value to increment a,b by
   *  amin,bmin = start values for a,b
   *  na,nb = number of repitions in a and b
   *  x,y,z = just pass x,y, and z by reference (the same ones passed for a,b)
   *  k = copy number to start with
   *  sd = sensitive detector for pmts
   */
  a = amin;
  for(G4int j = 1; j <= na; ++j)
  {
    a += da;
    b = bmin;
    for(G4int i = 1; i <= nb; ++i)
    {
      b += db;
      new G4PVPlacement(rot, G4ThreeVector(x, y, z), pmt_log, "pmt",
                        fHousing_log, false, k);
      fPmtPositions.push_back(G4ThreeVector(x, y, z));
      ++k;
    }
  }
}

◆ CopyValues()

void LXeMainVolume::CopyValues ( )

private

Definition at line 160 of file LXeMainVolume.cc.

{
  fScint_x         = fConstructor->GetScintX();
  fScint_y         = fConstructor->GetScintY();
  fScint_z         = fConstructor->GetScintZ();
  fD_mtl           = fConstructor->GetHousingThickness();
  fNx              = fConstructor->GetNX();
  fNy              = fConstructor->GetNY();
  fNz              = fConstructor->GetNZ();
  fOuterRadius_pmt = fConstructor->GetPMTRadius();
  fSphereOn        = fConstructor->GetSphereOn();
  fRefl            = fConstructor->GetHousingReflectivity();
}

Member Data Documentation

◆ fConstructor

LXeDetectorConstruction* LXeMainVolume::fConstructor = nullptr

private

Definition at line 65 of file LXeMainVolume.hh.

◆ fScint_x

G4double LXeMainVolume::fScint_x = 0.

private

Definition at line 67 of file LXeMainVolume.hh.

◆ fScint_y

G4double LXeMainVolume::fScint_y = 0.

private

Definition at line 68 of file LXeMainVolume.hh.

◆ fScint_z

G4double LXeMainVolume::fScint_z = 0.

private

Definition at line 69 of file LXeMainVolume.hh.

◆ fD_mtl

G4double LXeMainVolume::fD_mtl = 0.

private

Definition at line 70 of file LXeMainVolume.hh.

◆ fNx

G4int LXeMainVolume::fNx = 0

private

Definition at line 71 of file LXeMainVolume.hh.

◆ fNy

G4int LXeMainVolume::fNy = 0

private

Definition at line 72 of file LXeMainVolume.hh.

◆ fNz

G4int LXeMainVolume::fNz = 0

private

Definition at line 73 of file LXeMainVolume.hh.

◆ fOuterRadius_pmt

G4double LXeMainVolume::fOuterRadius_pmt = 0.

private

Definition at line 74 of file LXeMainVolume.hh.

◆ fSphereOn

G4bool LXeMainVolume::fSphereOn = false

private

Definition at line 75 of file LXeMainVolume.hh.

◆ fRefl

G4double LXeMainVolume::fRefl = 0.

private

Definition at line 76 of file LXeMainVolume.hh.

◆ fScint_box

G4Box* LXeMainVolume::fScint_box = nullptr

private

Definition at line 80 of file LXeMainVolume.hh.

◆ fHousing_box

G4Box* LXeMainVolume::fHousing_box = nullptr

private

Definition at line 81 of file LXeMainVolume.hh.

◆ fPmt

G4Tubs* LXeMainVolume::fPmt = nullptr

private

Definition at line 82 of file LXeMainVolume.hh.

◆ fPhotocath

G4Tubs* LXeMainVolume::fPhotocath = nullptr

private

Definition at line 83 of file LXeMainVolume.hh.

◆ fSphere

G4Sphere* LXeMainVolume::fSphere = nullptr

private

Definition at line 84 of file LXeMainVolume.hh.

◆ fScint_log

G4LogicalVolume* LXeMainVolume::fScint_log = nullptr

private

Definition at line 88 of file LXeMainVolume.hh.

◆ fHousing_log

G4LogicalVolume* LXeMainVolume::fHousing_log = nullptr

private

Definition at line 89 of file LXeMainVolume.hh.

◆ fPmt_log

G4LogicalVolume* LXeMainVolume::fPmt_log = nullptr

private

Definition at line 90 of file LXeMainVolume.hh.

◆ fPhotocath_log

G4LogicalVolume* LXeMainVolume::fPhotocath_log = nullptr

private

Definition at line 91 of file LXeMainVolume.hh.

◆ fSphere_log

G4LogicalVolume* LXeMainVolume::fSphere_log = nullptr

private

Definition at line 92 of file LXeMainVolume.hh.

◆ fPmtPositions

std::vector<G4ThreeVector> LXeMainVolume::fPmtPositions

private

Definition at line 95 of file LXeMainVolume.hh.

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

.doxygen/Doxymodules_optical.h
extended/optical/LXe/include/LXeMainVolume.hh
extended/optical/LXe/src/LXeMainVolume.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ LXeMainVolume()

Member Function Documentation

◆ GetLogPhotoCath()

◆ GetLogScint()

◆ GetPmtPositions()

◆ VisAttributes()

◆ SurfaceProperties()

◆ PlacePMTs()

◆ CopyValues()

Member Data Documentation

◆ fConstructor

◆ fScint_x

◆ fScint_y

◆ fScint_z

◆ fD_mtl

◆ fNx

◆ fNy

◆ fNz

◆ fOuterRadius_pmt

◆ fSphereOn

◆ fRefl

◆ fScint_box

◆ fHousing_box

◆ fPmt

◆ fPhotocath

◆ fSphere

◆ fScint_log

◆ fHousing_log

◆ fPmt_log

◆ fPhotocath_log

◆ fSphere_log

◆ fPmtPositions