WLSDetectorConstruction Class Reference

#include <Doxymodules_optical.h>

Inheritance diagram for WLSDetectorConstruction:

Public Member Functions
	WLSDetectorConstruction ()
	~WLSDetectorConstruction () override
G4VPhysicalVolume *	Construct () override
G4VPhysicalVolume *	ConstructDetector ()
void	ConstructSDandField () override
void	SetWorldMaterial (G4String)
void	SetWLSFiberMaterial (G4String)
void	SetCoupleMaterial (G4String)
void	SetPhotonDetGeometry (G4String)
void	SetNumberOfCladding (G4int)
void	SetWLSLength (G4double)
void	SetWLSRadius (G4double)
void	SetClad1Radius (G4double)
void	SetClad2Radius (G4double)
void	SetPhotonDetHalfLength (G4double)
void	SetGap (G4double)
void	SetPhotonDetAlignment (G4double)
void	SetXYRatio (G4double)
void	SetSurfaceRoughness (G4double)
void	SetMirrorReflectivity (G4double)
void	SetMirrorPolish (G4double)
void	SetPhotonDetReflectivity (G4double)
void	SetPhotonDetPolish (G4double)
void	SetMirror (G4bool)
void	SetBarLength (G4double)
void	SetBarBase (G4double)
void	SetHoleRadius (G4double)
void	SetCoatingThickness (G4double)
void	SetCoatingRadius (G4double)
G4double	GetWLSFiberLength ()
G4double	GetWLSFiberEnd ()
G4double	GetWLSFiberRMax ()
G4double	GetSurfaceRoughness ()
G4bool	IsPerfectFiber ()
G4double	GetBarLength ()
G4double	GetBarBase ()
G4double	GetHoleRadius ()
G4double	GetHoleLength ()
G4double	GetFiberRadius ()
G4double	GetCoatingThickness ()
G4double	GetCoatingRadius ()
G4Material *	FindMaterial (G4String)

Private Member Functions
void	UpdateGeometryParameters ()

Private Attributes
std::vector< G4VisAttributes * >	fVisAttributes
WLSMaterials *	fMaterials = nullptr
G4LogicalVolume *	fLogicHole = nullptr
G4LogicalVolume *	fLogicWorld = nullptr
G4VPhysicalVolume *	fPhysiWorld = nullptr
G4VPhysicalVolume *	fPhysiHole = nullptr
G4double	fWorldSizeX = -1.
G4double	fWorldSizeY = -1.
G4double	fWorldSizeZ = -1.
G4double	fWLSfiberRX = -1.
G4double	fWLSfiberRY = 0.5 * CLHEP::mm
G4double	fWLSfiberZ = 1. * CLHEP::m
G4double	fClad1RX = -1.
G4double	fClad1RY = -1.
G4double	fClad1Z = -1.
G4double	fClad2RX = -1.
G4double	fClad2RY = -1.
G4double	fClad2Z = -1.
G4double	fClrfiberHalfL = -1.
G4double	fClrfiberZ = -1.
G4double	fCoupleRX = -1.
G4double	fCoupleRY = -1.
G4double	fCoupleZ = -1.
G4double	fMirrorRmax = -1.
G4double	fMirrorZ = 0.1 * CLHEP::mm
G4bool	fMirrorToggle = true
G4String	fMPPCShape = "Circle"
G4double	fMPPCHalfL = -1.
G4double	fMPPCZ = 0.05 * CLHEP::mm
G4double	fMPPCDist = 0.
G4double	fMPPCTheta = 0
G4double	fWLSfiberOrigin = 0.
G4double	fCoupleOrigin = 0.
G4double	fMirrorOrigin = 0.
G4double	fMPPCOriginX = 0.
G4double	fMPPCOriginZ = 0.
G4int	fNumOfCladLayers = 0
G4double	fMirrorPolish = 1.
G4double	fMirrorReflectivity = 1.
G4double	fMPPCPolish = 1.
G4double	fMPPCReflectivity = 0.
G4double	fExtrusionPolish = 1.
G4double	fExtrusionReflectivity = 1.
G4double	fSurfaceRoughness = 1.
G4double	fXYRatio = 1.
G4double	fBarLength = 1. * CLHEP::m
G4double	fBarBase = 9.6 * CLHEP::mm
G4double	fHoleRadius = 0.9 * CLHEP::mm
G4double	fHoleLength = -1.
G4double	fCoatingThickness = 0.25 * CLHEP::mm
G4double	fCoatingRadius = 1.875 * CLHEP::mm
WLSDetectorMessenger *	fDetectorMessenger = nullptr
G4Cache< WLSPhotonDetSD * >	fmppcSD

Detailed Description

Definition at line 85 of file Doxymodules_optical.h.

Constructor & Destructor Documentation

WLSDetectorConstruction()

WLSDetectorConstruction::WLSDetectorConstruction

(

)

Definition at line 65 of file WLSDetectorConstruction.cc.

  : fVisAttributes()
{
  fDetectorMessenger = new WLSDetectorMessenger(this);
 
  fMPPCHalfL = fWLSfiberRY;
  fClrfiberZ = fMPPCZ + 10. * nm;
  fHoleLength       = fBarLength;
}

~WLSDetectorConstruction()

WLSDetectorConstruction::~WLSDetectorConstruction ( )

override

Definition at line 77 of file WLSDetectorConstruction.cc.

{
  delete fDetectorMessenger;
  delete fMaterials;
  for (auto visAttributes: fVisAttributes)
  {
    delete visAttributes;
  }
}

Member Function Documentation

Construct()

G4VPhysicalVolume * WLSDetectorConstruction::Construct ( )

override

Definition at line 89 of file WLSDetectorConstruction.cc.

{
  if(fPhysiWorld)
  {
    G4GeometryManager::GetInstance()->OpenGeometry();
    G4PhysicalVolumeStore::GetInstance()->Clean();
    G4LogicalVolumeStore::GetInstance()->Clean();
    G4SolidStore::GetInstance()->Clean();
    G4LogicalSkinSurface::CleanSurfaceTable();
    G4LogicalBorderSurface::CleanSurfaceTable();
  }
 
  fMaterials = WLSMaterials::GetInstance();
  UpdateGeometryParameters();
 
  return ConstructDetector();
}

ConstructDetector()

G4VPhysicalVolume * WLSDetectorConstruction::ConstructDetector

(

)

Definition at line 109 of file WLSDetectorConstruction.cc.

{
 
  auto air = FindMaterial("G4_AIR");
  //G4cout << "\nMaterial Properties Table for G4_AIR:" << G4endl;
  //air->GetMaterialPropertiesTable()->DumpTable();
 
  //--------------------------------------------------
  // World
  //--------------------------------------------------
 
  G4VSolid* solidWorld =
    new G4Box("World", fWorldSizeX, fWorldSizeY, fWorldSizeZ);
 
  fLogicWorld =
    new G4LogicalVolume(solidWorld, air, "World");
 
  fPhysiWorld =
    new G4PVPlacement(nullptr, G4ThreeVector(), fLogicWorld, "World", nullptr,
                      false, 0);
 
  //--------------------------------------------------
  // Extrusion
  //--------------------------------------------------
 
  auto coating = FindMaterial("Coating");
 
  G4VSolid* solidExtrusion = new G4Box("Extrusion", GetBarBase() / 2.,
                                       GetBarBase() / 2., GetBarLength() / 2.);
 
  auto logicExtrusion =
    new G4LogicalVolume(solidExtrusion, coating, "Extrusion");
 
  auto TiO2Surface = new G4OpticalSurface(
    "TiO2Surface", glisur, ground, dielectric_metal, fExtrusionPolish);
 
  auto TiO2SurfaceProperty = new G4MaterialPropertiesTable();
 
  std::vector<G4double> p_TiO2 = { 2.00 * eV, 3.47 * eV };
 
  std::vector<G4double> refl_TiO2 = { fExtrusionReflectivity,
                                      fExtrusionReflectivity };
  std::vector<G4double> effi_TiO2 = { 0., 0. };
 
  TiO2SurfaceProperty->AddProperty("REFLECTIVITY", p_TiO2, refl_TiO2);
  TiO2SurfaceProperty->AddProperty("EFFICIENCY", p_TiO2, effi_TiO2);
 
  TiO2Surface->SetMaterialPropertiesTable(TiO2SurfaceProperty);
 
  new G4PVPlacement(nullptr, G4ThreeVector(), logicExtrusion, "Extrusion",
                    fLogicWorld, false, 0);
 
  new G4LogicalSkinSurface("TiO2Surface", logicExtrusion, TiO2Surface);
 
  //--------------------------------------------------
  // Scintillator
  //--------------------------------------------------
 
  auto polystyrene = FindMaterial("Polystyrene");
  //G4cout << "\nMaterial Properties Table for Polystyrene:" << G4endl;
  //polystyrene->GetMaterialPropertiesTable()->DumpTable();
 
  G4VSolid* solidScintillator =
    new G4Box("Scintillator",
              GetBarBase() / 2. - GetCoatingThickness() - GetCoatingRadius(),
              GetBarBase() / 2. - GetCoatingThickness() - GetCoatingRadius(),
              GetBarLength() / 2.);
 
  auto logicScintillator = new G4LogicalVolume(
    solidScintillator, polystyrene, "Scintillator");
 
  new G4PVPlacement(nullptr, G4ThreeVector(), logicScintillator, "Scintillator",
                    logicExtrusion, false, 0);
  
  G4LogicalVolume* logicScintSide = nullptr;
  G4LogicalVolume* logicScintCrnr = nullptr;
  if(GetCoatingRadius() > 0.)
  {
    G4VSolid* solidScintside =
      new G4Box("SideOfBar",
                GetBarBase() / 2. - GetCoatingThickness() - GetCoatingRadius(),
                GetCoatingRadius() / 2., GetBarLength() / 2.);
 
    G4VSolid* solidScintcrnr =
      new G4Tubs("CrnrOfBar", 0.0 * cm, GetCoatingRadius(), GetBarLength() / 2.,
                 0. * deg, 90. * deg);
 
    logicScintSide = new G4LogicalVolume(
      solidScintside, polystyrene, "SideOfBar");
 
    logicScintCrnr = new G4LogicalVolume(
      solidScintcrnr, polystyrene, "CrnrOfBar");
 
    G4double pos =
      GetBarBase() / 2. - GetCoatingThickness() - GetCoatingRadius() / 2.;
 
    new G4PVPlacement(nullptr, G4ThreeVector(0., -pos, 0.), logicScintSide,
                      "SideOfBar", logicExtrusion, false, 0);
 
    new G4PVPlacement(nullptr, G4ThreeVector(0., pos, 0.), logicScintSide,
                      "SideOfBar", logicExtrusion, false, 1);
 
    auto rot1 = new G4RotationMatrix();
    rot1->rotateZ(-90.*deg);
 
    new G4PVPlacement(rot1, G4ThreeVector(pos, 0., 0.), logicScintSide,
                      "SideOfBar", logicExtrusion, false, 2);
 
    new G4PVPlacement(rot1, G4ThreeVector(-pos, 0., 0.), logicScintSide,
                      "SideOfBar", logicExtrusion, false, 3);
 
    pos = GetBarBase() / 2. - GetCoatingThickness() - GetCoatingRadius();
 
    new G4PVPlacement(nullptr, G4ThreeVector(pos, pos, 0.), logicScintCrnr,
                      "CrnrOfBar", logicExtrusion, false, 0);
 
    new G4PVPlacement(rot1, G4ThreeVector(-pos, pos, 0.), logicScintCrnr,
                      "CrnrOfBar", logicExtrusion, false, 1);
 
    auto rot2 = new G4RotationMatrix();
    rot2->rotateZ(-180.*deg);
 
    new G4PVPlacement(rot2, G4ThreeVector(-pos, -pos, 0.), logicScintCrnr,
                      "CrnrOfBar", logicExtrusion, false, 2);
 
    auto rot3 = new G4RotationMatrix();
    rot3->rotateZ(-270.*deg);
 
    new G4PVPlacement(rot3, G4ThreeVector(pos, -pos, 0.), logicScintCrnr,
                      "CrnrOfBar", logicExtrusion, false, 3);
  }
 
  if(GetFiberRadius() < GetHoleRadius())
  {
    G4VSolid* solidHole = new G4Tubs(
      "Hole", 0., GetHoleRadius(), GetHoleLength() / 2., 0. * deg, 360. * deg);
 
    fLogicHole = new G4LogicalVolume(solidHole, air, "Hole");
 
    fPhysiHole = new G4PVPlacement(nullptr, G4ThreeVector(), fLogicHole, "Hole",
                                   logicScintillator, false, 0);
  }
 
  //--------------------------------------------------
  // Fiber
  //--------------------------------------------------
 
  if(!(fLogicHole) || !(fPhysiHole))
  {
    G4ExceptionDescription ed;
    ed << "The Fiber Hole has not been constructed";
    G4Exception("WLSDetectorConstruction", "wls001", FatalException, ed);
  }
 
  // Pointers to the most recently constructed volume
  G4LogicalVolume* logicPlacement   = fLogicHole;
  G4VPhysicalVolume* physiPlacement = fPhysiHole;
 
  //--------------------------------------------------
  // Fiber Construction
  //--------------------------------------------------
 
  // Boundary Surface Properties
  G4OpticalSurface* opSurface = nullptr;
 
  if(fSurfaceRoughness < 1.)
    opSurface = new G4OpticalSurface("RoughSurface", glisur, ground,
                                     dielectric_dielectric, fSurfaceRoughness);
  
  G4LogicalVolume* logicWLSfiber = nullptr;
  G4LogicalVolume* logicClad1    = nullptr;
  G4LogicalVolume* logicClad2    = nullptr;
  G4VPhysicalVolume* physiClad1  = nullptr;
  G4VPhysicalVolume* physiClad2  = nullptr;
 
  auto fpethylene = FindMaterial("FPethylene");
  auto pethylene = FindMaterial("Pethylene");
  auto pmma = FindMaterial("PMMA");
 
  // Determine the number of cladding layers to be built
  switch(fNumOfCladLayers)
  {
    case 2:
 
      //--------------------------------------------------
      // Cladding 2
      //--------------------------------------------------
 
      //G4cout << "\nMaterial Properties Table for fPethylene:" << G4endl;
      //fpethylene->GetMaterialPropertiesTable()->DumpTable();
 
      G4VSolid* solidClad2;
 
      if(fXYRatio == 1.)
        solidClad2 = new G4Tubs("Clad2", 0., fClad2RX, fClad2Z, 0., twopi);
      else
        solidClad2 = new G4EllipticalTube("Clad2", fClad2RX, fClad2RY, fClad2Z);
 
      logicClad2 =
        new G4LogicalVolume(solidClad2, fpethylene, "Clad2");
 
      physiClad2 =
        new G4PVPlacement(nullptr, G4ThreeVector(0.0, 0.0, fWLSfiberOrigin),
                          logicClad2, "Clad2", logicPlacement, false, 0);
 
      // Place the rough surface only if needed
      if(opSurface)
      {
        new G4LogicalBorderSurface("surfaceClad2Out", physiClad2,
                                   physiPlacement, opSurface);
        new G4LogicalBorderSurface("surfaceClad2In", physiPlacement, physiClad2,
                                   opSurface);
      }
 
      logicPlacement = logicClad2;
      physiPlacement = physiClad2;
      [[fallthrough]];
 
    case 1:
 
      //--------------------------------------------------
      // Cladding 1
      //--------------------------------------------------
 
      //G4cout << "\nMaterial Properties Table for Pethylene:" << G4endl;
      //pethylene->GetMaterialPropertiesTable()->DumpTable();
 
      G4VSolid* solidClad1;
 
      if(fXYRatio == 1.)
        solidClad1 = new G4Tubs("Clad1", 0., fClad1RX, fClad1Z, 0., twopi);
      else
        solidClad1 = new G4EllipticalTube("Clad1", fClad1RX, fClad1RY, fClad1Z);
 
      logicClad1 =
        new G4LogicalVolume(solidClad1, pethylene, "Clad1");
 
      physiClad1 =
        new G4PVPlacement(nullptr, G4ThreeVector(0., 0., fWLSfiberOrigin),
                          logicClad1, "Clad1", logicPlacement, false, 0);
 
      // Place the rough surface only if needed
      if(opSurface)
      {
        new G4LogicalBorderSurface("surfaceClad1Out", physiClad1,
                                   physiPlacement, opSurface);
 
        new G4LogicalBorderSurface("surfaceClad1In", physiPlacement, physiClad1,
                                   opSurface);
      }
 
      logicPlacement = logicClad1;
      physiPlacement = physiClad1;
      [[fallthrough]];
 
    default:
 
      //--------------------------------------------------
      // WLS Fiber
      //--------------------------------------------------
 
      //G4cout << "\nMaterial Properties Table for PMMA:" << G4endl;
      //pmma->GetMaterialPropertiesTable()->DumpTable();
 
      G4VSolid* solidWLSfiber;
 
      if(fXYRatio == 1.)
      {
        solidWLSfiber =
          new G4Tubs("WLSFiber", 0., fWLSfiberRX, fWLSfiberZ, 0., twopi);
      }
      else
      {
        solidWLSfiber = new G4EllipticalTube("WLSFiber", fWLSfiberRX,
                                             fWLSfiberRY, fWLSfiberZ);
      }
 
      logicWLSfiber =
        new G4LogicalVolume(solidWLSfiber, pmma, "WLSFiber");
 
      logicWLSfiber->SetUserLimits(
        new G4UserLimits(DBL_MAX, DBL_MAX, 10. * ms));
 
      G4VPhysicalVolume* physiWLSfiber =
        new G4PVPlacement(nullptr, G4ThreeVector(0., 0., fWLSfiberOrigin),
                          logicWLSfiber, "WLSFiber", logicPlacement, false, 0);
 
      // Place the rough surface only if needed
      if(opSurface)
      {
        new G4LogicalBorderSurface("surfaceWLSOut", physiWLSfiber,
                                   physiPlacement, opSurface);
 
        new G4LogicalBorderSurface("surfaceWLSIn", physiPlacement,
                                   physiWLSfiber, opSurface);
      }
  }
 
  //--------------------------------------------------
  // Mirror for reflection at one of the end
  //--------------------------------------------------
 
  // Place the mirror only if the user wants the mirror
  G4LogicalVolume* logicMirror = nullptr;
 
  auto aluminum = FindMaterial("G4_Al");
 
  if(fMirrorToggle)
  {
    G4VSolid* solidMirror =
      new G4Box("Mirror", fMirrorRmax, fMirrorRmax, fMirrorZ);
 
    logicMirror =
      new G4LogicalVolume(solidMirror, aluminum, "Mirror");
 
    auto mirrorSurface = new G4OpticalSurface(
      "MirrorSurface", glisur, ground, dielectric_metal, fMirrorPolish);
 
    auto mirrorSurfaceProperty = new G4MaterialPropertiesTable();
 
    std::vector<G4double> p_mirror    = { 2.00 * eV, 3.47 * eV };
    std::vector<G4double> refl_mirror = { fMirrorReflectivity,
                                          fMirrorReflectivity };
    std::vector<G4double> effi_mirror = { 0., 0. };
 
    mirrorSurfaceProperty->AddProperty("REFLECTIVITY", p_mirror, refl_mirror);
    mirrorSurfaceProperty->AddProperty("EFFICIENCY", p_mirror, effi_mirror);
 
    mirrorSurface->SetMaterialPropertiesTable(mirrorSurfaceProperty);
 
    new G4PVPlacement(nullptr, G4ThreeVector(0., 0., fMirrorOrigin),
                      logicMirror, "Mirror", fLogicWorld, false, 0);
 
    new G4LogicalSkinSurface("MirrorSurface", logicMirror, mirrorSurface);
  }
 
  //--------------------------------------------------
  // Coupling at the read-out end
  //--------------------------------------------------
 
  // Clear Fiber (Coupling Layer)
  G4VSolid* solidCouple = new G4Box("Couple", fCoupleRX, fCoupleRY, fCoupleZ);
 
  auto logicCouple = new G4LogicalVolume(solidCouple, air, "Couple");
 
  new G4PVPlacement(nullptr, G4ThreeVector(0., 0., fCoupleOrigin), logicCouple,
                    "Couple", fLogicWorld, false, 0);
 
  //--------------------------------------------------
  // A logical layer in front of PhotonDet
  //--------------------------------------------------
 
  // Purpose: Preventing direct dielectric to metal contact
 
  // Check for valid placement of PhotonDet
  if(fMPPCTheta > std::atan(fMPPCDist / fMPPCHalfL))
  {
    fMPPCTheta   = 0.;
    fMPPCOriginX = std::sin(fMPPCTheta) * (fMPPCDist + fClrfiberZ);
    fMPPCOriginZ = -fCoupleZ + std::cos(fMPPCTheta) * (fMPPCDist + fClrfiberZ);
    G4ExceptionDescription ed;
    ed << "Invalid alignment.  Alignment reset to 0.";
    G4Exception("WLSDetectorConstruction", "wls002", JustWarning, ed);
  }
 
  // Clear Fiber (Coupling Layer)
  G4VSolid* solidClrfiber;
 
  if(fMPPCShape == "Square")
  {
    solidClrfiber =
      new G4Box("ClearFiber", fClrfiberHalfL, fClrfiberHalfL, fClrfiberZ);
  }
  else
  {
    solidClrfiber =
      new G4Tubs("ClearFiber", 0., fClrfiberHalfL, fClrfiberZ, 0., twopi);
  }
 
  auto logicClrfiber = new G4LogicalVolume(solidClrfiber, air, "ClearFiber");
 
  new G4PVPlacement(new G4RotationMatrix(CLHEP::HepRotationY(-fMPPCTheta)),
                    G4ThreeVector(fMPPCOriginX, 0.0, fMPPCOriginZ),
                    logicClrfiber, "ClearFiber", logicCouple, false, 0);
 
  //--------------------------------------------------
  // PhotonDet (Sensitive Detector)
  //--------------------------------------------------
 
  // Physical Construction
  G4VSolid* solidPhotonDet = nullptr;
 
  if(fMPPCShape == "Square")
    solidPhotonDet = new G4Box("PhotonDet", fMPPCHalfL, fMPPCHalfL, fMPPCZ);
  else
    solidPhotonDet = new G4Tubs("PhotonDet", 0., fMPPCHalfL, fMPPCZ, 0., twopi);
 
  auto logicPhotonDet =
    new G4LogicalVolume(solidPhotonDet, aluminum, "PhotonDet_LV");
 
  new G4PVPlacement(nullptr, G4ThreeVector(0., 0., 0.), logicPhotonDet,
                    "PhotonDet", logicClrfiber, false, 0);
 
  // PhotonDet Surface Properties
  auto photonDetSurface = new G4OpticalSurface(
    "PhotonDetSurface", glisur, ground, dielectric_metal, fMPPCPolish);
 
  auto photonDetSurfaceProperty = new G4MaterialPropertiesTable();
 
  std::vector<G4double> p_mppc    = { 2.00 * eV, 3.47 * eV };
  std::vector<G4double> refl_mppc = { fMPPCReflectivity, fMPPCReflectivity };
  std::vector<G4double> effi_mppc = { 1., 1. };
 
  photonDetSurfaceProperty->AddProperty("REFLECTIVITY", p_mppc, refl_mppc);
  photonDetSurfaceProperty->AddProperty("EFFICIENCY", p_mppc, effi_mppc);
 
  photonDetSurface->SetMaterialPropertiesTable(photonDetSurfaceProperty);
 
  new G4LogicalSkinSurface("PhotonDetSurface", logicPhotonDet,
                           photonDetSurface);
 
  // visualization attributes -------------------------------------------------
 
  auto visAttributes = new G4VisAttributes(G4Colour(1.0,1.0,1.0));
  visAttributes->SetVisibility(false);
  fLogicWorld->SetVisAttributes(visAttributes);
  fVisAttributes.push_back(visAttributes);
 
  visAttributes = new G4VisAttributes(G4Colour(0.2,0.2,0.2,0.5));
  visAttributes->SetVisibility(true);
  logicExtrusion->SetVisAttributes(visAttributes);
  fVisAttributes.push_back(visAttributes);
 
  visAttributes = new G4VisAttributes(G4Colour(0.0,0.0,1.0,0.9));
  visAttributes->SetVisibility(true);
  logicScintillator->SetVisAttributes(visAttributes);
  fVisAttributes.push_back(visAttributes);
 
  visAttributes = new G4VisAttributes(G4Colour(0.0,0.8,0.2,0.2));
  visAttributes->SetVisibility(true);
  logicScintSide->SetVisAttributes(visAttributes);
  fVisAttributes.push_back(visAttributes);
 
  visAttributes = new G4VisAttributes(G4Colour(0.0,0.8,0.2,0.2));
  visAttributes->SetVisibility(true);
  logicScintCrnr->SetVisAttributes(visAttributes);
  fVisAttributes.push_back(visAttributes);
 
  visAttributes = new G4VisAttributes(G4Colour(0.4,0.0,0.0,0.5));
  visAttributes->SetVisibility(true);
  fLogicHole->SetVisAttributes(visAttributes);
  fVisAttributes.push_back(visAttributes);
 
  if(logicClad1 != nullptr)
  {
    visAttributes = new G4VisAttributes(G4Colour(0.0,0.8,0.5,0.5));
    visAttributes->SetVisibility(true);
    logicClad1->SetVisAttributes(visAttributes);
    fVisAttributes.push_back(visAttributes);
  }
 
  if(logicClad2 != nullptr)
  {
    visAttributes = new G4VisAttributes(G4Colour(0.0,0.5,0.8,0.5));
    visAttributes->SetVisibility(true);
    logicClad2->SetVisAttributes(visAttributes);
    fVisAttributes.push_back(visAttributes);
  }
 
  visAttributes = new G4VisAttributes(G4Colour(0.8,0.8,1.0));
  visAttributes->SetVisibility(true);
  logicWLSfiber->SetVisAttributes(visAttributes);
  fVisAttributes.push_back(visAttributes);
 
  if(fMirrorToggle == true)
  {
    visAttributes = new G4VisAttributes(G4Colour(0.3,0.3,1.0,0.3));
    visAttributes->SetVisibility(true);
    logicMirror->SetVisAttributes(visAttributes);
    fVisAttributes.push_back(visAttributes);
  }
 
  visAttributes = new G4VisAttributes(G4Colour(0.0,0.0,0.5,0.5));
  visAttributes->SetVisibility(true);
  logicCouple->SetVisAttributes(visAttributes);
  fVisAttributes.push_back(visAttributes);
 
  visAttributes = new G4VisAttributes(G4Colour(0.3,0.3,0.3,0.5));
  visAttributes->SetVisibility(true);
  logicClrfiber->SetVisAttributes(visAttributes);
  fVisAttributes.push_back(visAttributes);
 
  visAttributes = new G4VisAttributes(G4Colour(1.0,1.0,1.0,0.8));
  visAttributes->SetVisibility(true);
  logicPhotonDet->SetVisAttributes(visAttributes);
  fVisAttributes.push_back(visAttributes);
 
  return fPhysiWorld;
}

ConstructSDandField()

void WLSDetectorConstruction::ConstructSDandField ( )

override

Definition at line 613 of file WLSDetectorConstruction.cc.

{
  if(!fmppcSD.Get())
  {
    G4String mppcSDName    = "WLS/PhotonDet";
    auto mppcSD = new WLSPhotonDetSD(mppcSDName);
    G4SDManager::GetSDMpointer()->AddNewDetector(mppcSD);
    fmppcSD.Put(mppcSD);
  }
  SetSensitiveDetector("PhotonDet_LV", fmppcSD.Get(), true);
}

SetWorldMaterial()

void WLSDetectorConstruction::SetWorldMaterial

(

G4String

)

SetWLSFiberMaterial()

void WLSDetectorConstruction::SetWLSFiberMaterial

(

G4String

)

SetCoupleMaterial()

void WLSDetectorConstruction::SetCoupleMaterial

(

G4String

)

SetPhotonDetGeometry()

void WLSDetectorConstruction::SetPhotonDetGeometry

(

G4String

shape

)

Definition at line 659 of file WLSDetectorConstruction.cc.

{
  if(shape == "Circle" || shape == "Square")
    fMPPCShape = shape;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetNumberOfCladding()

void WLSDetectorConstruction::SetNumberOfCladding

(

G4int

num

)

Definition at line 670 of file WLSDetectorConstruction.cc.

{
  fNumOfCladLayers = num;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetWLSLength()

void WLSDetectorConstruction::SetWLSLength

(

G4double

length

)

Definition at line 680 of file WLSDetectorConstruction.cc.

{
  fWLSfiberZ = length;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetWLSRadius()

void WLSDetectorConstruction::SetWLSRadius

(

G4double

radius

)

Definition at line 689 of file WLSDetectorConstruction.cc.

{
  fWLSfiberRY = radius;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetClad1Radius()

void WLSDetectorConstruction::SetClad1Radius

(

G4double

radius

)

Definition at line 698 of file WLSDetectorConstruction.cc.

{
  fClad1RY = radius;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetClad2Radius()

void WLSDetectorConstruction::SetClad2Radius

(

G4double

radius

)

Definition at line 707 of file WLSDetectorConstruction.cc.

{
  fClad2RY = radius;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetPhotonDetHalfLength()

void WLSDetectorConstruction::SetPhotonDetHalfLength

(

G4double

halfL

)

Definition at line 716 of file WLSDetectorConstruction.cc.

{
  fMPPCHalfL = halfL;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetGap()

void WLSDetectorConstruction::SetGap

(

G4double

gap

)

Definition at line 726 of file WLSDetectorConstruction.cc.

{
  fMPPCDist = gap;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetPhotonDetAlignment()

void WLSDetectorConstruction::SetPhotonDetAlignment

(

G4double

theta

)

Definition at line 735 of file WLSDetectorConstruction.cc.

{
  fMPPCTheta = theta;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetXYRatio()

void WLSDetectorConstruction::SetXYRatio

(

G4double

r

)

Definition at line 807 of file WLSDetectorConstruction.cc.

{
  fXYRatio = r;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetSurfaceRoughness()

void WLSDetectorConstruction::SetSurfaceRoughness

(

G4double

roughness

)

Definition at line 747 of file WLSDetectorConstruction.cc.

{
  fSurfaceRoughness = roughness;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetMirrorReflectivity()

void WLSDetectorConstruction::SetMirrorReflectivity

(

G4double

reflectivity

)

Definition at line 767 of file WLSDetectorConstruction.cc.

{
  fMirrorReflectivity = reflectivity;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetMirrorPolish()

void WLSDetectorConstruction::SetMirrorPolish

(

G4double

polish

)

Definition at line 757 of file WLSDetectorConstruction.cc.

{
  fMirrorPolish = polish;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetPhotonDetReflectivity()

void WLSDetectorConstruction::SetPhotonDetReflectivity

(

G4double

reflectivity

)

Definition at line 787 of file WLSDetectorConstruction.cc.

{
  fMPPCReflectivity = reflectivity;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetPhotonDetPolish()

void WLSDetectorConstruction::SetPhotonDetPolish

(

G4double

polish

)

Definition at line 777 of file WLSDetectorConstruction.cc.

{
  fMPPCPolish = polish;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetMirror()

void WLSDetectorConstruction::SetMirror

(

G4bool

flag

)

Definition at line 797 of file WLSDetectorConstruction.cc.

{
  fMirrorToggle = flag;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetBarLength()

void WLSDetectorConstruction::SetBarLength

(

G4double

length

)

Definition at line 817 of file WLSDetectorConstruction.cc.

{
  fBarLength = length;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetBarBase()

void WLSDetectorConstruction::SetBarBase

(

G4double

side

)

Definition at line 826 of file WLSDetectorConstruction.cc.

{
  fBarBase = side;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetHoleRadius()

void WLSDetectorConstruction::SetHoleRadius

(

G4double

radius

)

Definition at line 835 of file WLSDetectorConstruction.cc.

{
  fHoleRadius = radius;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetCoatingThickness()

void WLSDetectorConstruction::SetCoatingThickness

(

G4double

thick

)

Definition at line 844 of file WLSDetectorConstruction.cc.

{
  fCoatingThickness = thick;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

SetCoatingRadius()

void WLSDetectorConstruction::SetCoatingRadius

(

G4double

radius

)

Definition at line 853 of file WLSDetectorConstruction.cc.

{
  fCoatingRadius = radius;
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

GetWLSFiberLength()

G4double WLSDetectorConstruction::GetWLSFiberLength

(

)

Definition at line 862 of file WLSDetectorConstruction.cc.

{ return fWLSfiberZ; }

GetWLSFiberEnd()

G4double WLSDetectorConstruction::GetWLSFiberEnd

(

)

Definition at line 897 of file WLSDetectorConstruction.cc.

{
  return fWLSfiberOrigin + fWLSfiberZ;
}

GetWLSFiberRMax()

G4double WLSDetectorConstruction::GetWLSFiberRMax

(

)

Definition at line 904 of file WLSDetectorConstruction.cc.

{
  if(fNumOfCladLayers == 2)
    return fClad2RY;
  if(fNumOfCladLayers == 1)
    return fClad1RY;
  return fWLSfiberRY;
}

GetSurfaceRoughness()

G4double WLSDetectorConstruction::GetSurfaceRoughness

(

)

Definition at line 915 of file WLSDetectorConstruction.cc.

{
  return fSurfaceRoughness;
}

IsPerfectFiber()

G4bool WLSDetectorConstruction::IsPerfectFiber

(

)

Definition at line 923 of file WLSDetectorConstruction.cc.

{
  return fSurfaceRoughness == 1. && fXYRatio == 1. &&
         (!fMirrorToggle || (fMirrorPolish == 1. && fMirrorReflectivity == 1.));
}

GetBarLength()

G4double WLSDetectorConstruction::GetBarLength

(

)

Definition at line 866 of file WLSDetectorConstruction.cc.

{ return fBarLength; }

GetBarBase()

G4double WLSDetectorConstruction::GetBarBase

(

)

Definition at line 870 of file WLSDetectorConstruction.cc.

{ return fBarBase; }

GetHoleRadius()

G4double WLSDetectorConstruction::GetHoleRadius

(

)

Definition at line 874 of file WLSDetectorConstruction.cc.

{ return fHoleRadius; }

GetHoleLength()

G4double WLSDetectorConstruction::GetHoleLength

(

)

Definition at line 878 of file WLSDetectorConstruction.cc.

{ return fHoleLength; }

GetFiberRadius()

G4double WLSDetectorConstruction::GetFiberRadius

(

)

Definition at line 882 of file WLSDetectorConstruction.cc.

{ return GetWLSFiberRMax(); }

GetCoatingThickness()

G4double WLSDetectorConstruction::GetCoatingThickness

(

)

Definition at line 886 of file WLSDetectorConstruction.cc.

{
  return fCoatingThickness;
}

GetCoatingRadius()

G4double WLSDetectorConstruction::GetCoatingRadius

(

)

Definition at line 893 of file WLSDetectorConstruction.cc.

{ return fCoatingRadius; }

FindMaterial()

G4Material * WLSDetectorConstruction::FindMaterial

(

G4String

name

)

Definition at line 931 of file WLSDetectorConstruction.cc.

{
  G4Material* material = G4Material::GetMaterial(name, true);
  return material;
}

UpdateGeometryParameters()

void WLSDetectorConstruction::UpdateGeometryParameters ( )

private

Definition at line 627 of file WLSDetectorConstruction.cc.

{
  fWLSfiberRX = fXYRatio * fWLSfiberRY;
 
  fClad1RX = fWLSfiberRX + 0.03 * fWLSfiberRX;
  fClad1RY = fWLSfiberRY + 0.03 * fWLSfiberRY;
  fClad1Z  = fWLSfiberZ;
 
  fClad2RX = fClad1RX + 0.03 * fWLSfiberRX;
  fClad2RY = fClad1RY + 0.03 * fWLSfiberRY;
  fClad2Z  = fWLSfiberZ;
 
  fWorldSizeX = fClad2RX + fMPPCDist + fMPPCHalfL + 1. * cm;
  fWorldSizeY = fClad2RY + fMPPCDist + fMPPCHalfL + 1. * cm;
  fWorldSizeZ = fWLSfiberZ + fMPPCDist + fMPPCHalfL + 1. * cm;
 
  fCoupleRX = fWorldSizeX;
  fCoupleRY = fWorldSizeY;
  fCoupleZ  = (fWorldSizeZ - fWLSfiberZ) / 2.;
 
  fClrfiberHalfL = fMPPCHalfL;
 
  fMirrorRmax = fClad2RY;
 
  fCoupleOrigin = fWLSfiberOrigin + fWLSfiberZ + fCoupleZ;
  fMirrorOrigin = fWLSfiberOrigin - fWLSfiberZ - fMirrorZ;
  fMPPCOriginX  = std::sin(fMPPCTheta) * (fMPPCDist + fClrfiberZ);
  fMPPCOriginZ  = -fCoupleZ + std::cos(fMPPCTheta) * (fMPPCDist + fClrfiberZ);
}

Member Data Documentation

fVisAttributes

std::vector<G4VisAttributes*> WLSDetectorConstruction::fVisAttributes

private

Definition at line 120 of file WLSDetectorConstruction.hh.

fMaterials

WLSMaterials* WLSDetectorConstruction::fMaterials = nullptr

private

Definition at line 122 of file WLSDetectorConstruction.hh.

fLogicHole

G4LogicalVolume* WLSDetectorConstruction::fLogicHole = nullptr

private

Definition at line 124 of file WLSDetectorConstruction.hh.

fLogicWorld

G4LogicalVolume* WLSDetectorConstruction::fLogicWorld = nullptr

private

Definition at line 125 of file WLSDetectorConstruction.hh.

fPhysiWorld

G4VPhysicalVolume* WLSDetectorConstruction::fPhysiWorld = nullptr

private

Definition at line 127 of file WLSDetectorConstruction.hh.

fPhysiHole

G4VPhysicalVolume* WLSDetectorConstruction::fPhysiHole = nullptr

private

Definition at line 128 of file WLSDetectorConstruction.hh.

fWorldSizeX

G4double WLSDetectorConstruction::fWorldSizeX = -1.

private

Definition at line 130 of file WLSDetectorConstruction.hh.

fWorldSizeY

G4double WLSDetectorConstruction::fWorldSizeY = -1.

private

Definition at line 131 of file WLSDetectorConstruction.hh.

fWorldSizeZ

G4double WLSDetectorConstruction::fWorldSizeZ = -1.

private

Definition at line 132 of file WLSDetectorConstruction.hh.

fWLSfiberRX

G4double WLSDetectorConstruction::fWLSfiberRX = -1.

private

Definition at line 134 of file WLSDetectorConstruction.hh.

fWLSfiberRY

G4double WLSDetectorConstruction::fWLSfiberRY = 0.5 * CLHEP::mm

private

Definition at line 135 of file WLSDetectorConstruction.hh.

fWLSfiberZ

G4double WLSDetectorConstruction::fWLSfiberZ = 1. * CLHEP::m

private

Definition at line 136 of file WLSDetectorConstruction.hh.

fClad1RX

G4double WLSDetectorConstruction::fClad1RX = -1.

private

Definition at line 138 of file WLSDetectorConstruction.hh.

fClad1RY

G4double WLSDetectorConstruction::fClad1RY = -1.

private

Definition at line 139 of file WLSDetectorConstruction.hh.

fClad1Z

G4double WLSDetectorConstruction::fClad1Z = -1.

private

Definition at line 140 of file WLSDetectorConstruction.hh.

fClad2RX

G4double WLSDetectorConstruction::fClad2RX = -1.

private

Definition at line 142 of file WLSDetectorConstruction.hh.

fClad2RY

G4double WLSDetectorConstruction::fClad2RY = -1.

private

Definition at line 143 of file WLSDetectorConstruction.hh.

fClad2Z

G4double WLSDetectorConstruction::fClad2Z = -1.

private

Definition at line 144 of file WLSDetectorConstruction.hh.

fClrfiberHalfL

G4double WLSDetectorConstruction::fClrfiberHalfL = -1.

private

Definition at line 146 of file WLSDetectorConstruction.hh.

fClrfiberZ

G4double WLSDetectorConstruction::fClrfiberZ = -1.

private

Definition at line 147 of file WLSDetectorConstruction.hh.

fCoupleRX

G4double WLSDetectorConstruction::fCoupleRX = -1.

private

Definition at line 149 of file WLSDetectorConstruction.hh.

fCoupleRY

G4double WLSDetectorConstruction::fCoupleRY = -1.

private

Definition at line 150 of file WLSDetectorConstruction.hh.

fCoupleZ

G4double WLSDetectorConstruction::fCoupleZ = -1.

private

Definition at line 151 of file WLSDetectorConstruction.hh.

fMirrorRmax

G4double WLSDetectorConstruction::fMirrorRmax = -1.

private

Definition at line 153 of file WLSDetectorConstruction.hh.

fMirrorZ

G4double WLSDetectorConstruction::fMirrorZ = 0.1 * CLHEP::mm

private

Definition at line 154 of file WLSDetectorConstruction.hh.

fMirrorToggle

G4bool WLSDetectorConstruction::fMirrorToggle = true

private

Definition at line 155 of file WLSDetectorConstruction.hh.

fMPPCShape

G4String WLSDetectorConstruction::fMPPCShape = "Circle"

private

Definition at line 157 of file WLSDetectorConstruction.hh.

fMPPCHalfL

G4double WLSDetectorConstruction::fMPPCHalfL = -1.

private

Definition at line 158 of file WLSDetectorConstruction.hh.

fMPPCZ

G4double WLSDetectorConstruction::fMPPCZ = 0.05 * CLHEP::mm

private

Definition at line 159 of file WLSDetectorConstruction.hh.

fMPPCDist

G4double WLSDetectorConstruction::fMPPCDist = 0.

private

Definition at line 160 of file WLSDetectorConstruction.hh.

fMPPCTheta

G4double WLSDetectorConstruction::fMPPCTheta = 0

private

Definition at line 161 of file WLSDetectorConstruction.hh.

fWLSfiberOrigin

G4double WLSDetectorConstruction::fWLSfiberOrigin = 0.

private

Definition at line 163 of file WLSDetectorConstruction.hh.

fCoupleOrigin

G4double WLSDetectorConstruction::fCoupleOrigin = 0.

private

Definition at line 164 of file WLSDetectorConstruction.hh.

fMirrorOrigin

G4double WLSDetectorConstruction::fMirrorOrigin = 0.

private

Definition at line 165 of file WLSDetectorConstruction.hh.

fMPPCOriginX

G4double WLSDetectorConstruction::fMPPCOriginX = 0.

private

Definition at line 166 of file WLSDetectorConstruction.hh.

fMPPCOriginZ

G4double WLSDetectorConstruction::fMPPCOriginZ = 0.

private

Definition at line 167 of file WLSDetectorConstruction.hh.

fNumOfCladLayers

G4int WLSDetectorConstruction::fNumOfCladLayers = 0

private

Definition at line 169 of file WLSDetectorConstruction.hh.

fMirrorPolish

G4double WLSDetectorConstruction::fMirrorPolish = 1.

private

Definition at line 171 of file WLSDetectorConstruction.hh.

fMirrorReflectivity

G4double WLSDetectorConstruction::fMirrorReflectivity = 1.

private

Definition at line 172 of file WLSDetectorConstruction.hh.

fMPPCPolish

G4double WLSDetectorConstruction::fMPPCPolish = 1.

private

Definition at line 173 of file WLSDetectorConstruction.hh.

fMPPCReflectivity

G4double WLSDetectorConstruction::fMPPCReflectivity = 0.

private

Definition at line 174 of file WLSDetectorConstruction.hh.

fExtrusionPolish

G4double WLSDetectorConstruction::fExtrusionPolish = 1.

private

Definition at line 175 of file WLSDetectorConstruction.hh.

fExtrusionReflectivity

G4double WLSDetectorConstruction::fExtrusionReflectivity = 1.

private

Definition at line 176 of file WLSDetectorConstruction.hh.

fSurfaceRoughness

G4double WLSDetectorConstruction::fSurfaceRoughness = 1.

private

Definition at line 177 of file WLSDetectorConstruction.hh.

fXYRatio

G4double WLSDetectorConstruction::fXYRatio = 1.

private

Definition at line 178 of file WLSDetectorConstruction.hh.

fBarLength

G4double WLSDetectorConstruction::fBarLength = 1. * CLHEP::m

private

Definition at line 180 of file WLSDetectorConstruction.hh.

fBarBase

G4double WLSDetectorConstruction::fBarBase = 9.6 * CLHEP::mm

private

Definition at line 181 of file WLSDetectorConstruction.hh.

fHoleRadius

G4double WLSDetectorConstruction::fHoleRadius = 0.9 * CLHEP::mm

private

Definition at line 182 of file WLSDetectorConstruction.hh.

fHoleLength

G4double WLSDetectorConstruction::fHoleLength = -1.

private

Definition at line 183 of file WLSDetectorConstruction.hh.

fCoatingThickness

G4double WLSDetectorConstruction::fCoatingThickness = 0.25 * CLHEP::mm

private

Definition at line 184 of file WLSDetectorConstruction.hh.

fCoatingRadius

G4double WLSDetectorConstruction::fCoatingRadius = 1.875 * CLHEP::mm

private

Definition at line 185 of file WLSDetectorConstruction.hh.

fDetectorMessenger

WLSDetectorMessenger* WLSDetectorConstruction::fDetectorMessenger = nullptr

private

Definition at line 189 of file WLSDetectorConstruction.hh.

fmppcSD

G4Cache<WLSPhotonDetSD*> WLSDetectorConstruction::fmppcSD

private

Definition at line 190 of file WLSDetectorConstruction.hh.

---

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

• .doxygen/Doxymodules_optical.h
• extended/optical/wls/include/WLSDetectorConstruction.hh
• extended/optical/wls/src/WLSDetectorConstruction.cc