GB03DetectorConstruction Class Reference

#include <Doxymodules_biasing.h>

Inheritance diagram for GB03DetectorConstruction:

Public Member Functions
	GB03DetectorConstruction ()
virtual	~GB03DetectorConstruction ()
virtual G4VPhysicalVolume *	Construct ()
virtual void	ConstructSDandField ()
void	PrintCalorParameters () const
void	SetAbsorberMaterial (G4String materialChoice)
G4String	GetAbsorberMaterial () const
void	SetGapMaterial (G4String materialChoice)
G4String	GetGapMaterial () const
void	SetNumberOfLayers (G4int nl)
G4int	GetVerboseLevel () const
void	SetVerboseLevel (G4int val)

Static Public Member Functions
static G4int	GetNumberOfLayers ()

Private Member Functions
void	DefineMaterials ()
void	SetupGeometry ()
void	SetupDetectors ()
void	SetupBiasing ()

Private Attributes
G4double	fTotalThickness
G4double	fLayerThickness
	total thinkness of one calorimeter
G4bool	fConstructed
	= fTotalThickness / fNumberOfLayers
G4String	fCalName
G4Material *	fWorldMaterial
G4Material *	fAbsorberMaterial
G4Material *	fGapMaterial
G4Box *	fLayerSolid
G4Box *	fGapSolid
G4LogicalVolume *	fWorldLogical
G4LogicalVolume *	fCalorLogical
G4LogicalVolume *	fLayerLogical
G4LogicalVolume *	fGapLogical
G4VPhysicalVolume *	fWorldPhysical
G4VPhysicalVolume *	fCalorPhysical
G4PVReplica *	fLayerPhysical
G4VPhysicalVolume *	fGapPhysical
GB03DetectorMessenger *	fDetectorMessenger
G4int	fVerboseLevel

Static Private Attributes
static G4int	fNumberOfLayers = 40
static G4ThreadLocal G4bool	fConstructedSDandField = false

Detailed Description

Definition at line 96 of file Doxymodules_biasing.h.

Constructor & Destructor Documentation

GB03DetectorConstruction()

GB03DetectorConstruction::GB03DetectorConstruction

(

)

Definition at line 63 of file GB03DetectorConstruction.cc.

: G4VUserDetectorConstruction(),
  fTotalThickness (2.0*m), fLayerThickness(0.),
  fConstructed(false), 
  fWorldMaterial(0), fAbsorberMaterial(0), fGapMaterial(0),
  fLayerSolid(0), fGapSolid(0),
  fWorldLogical(0), fCalorLogical(0), fLayerLogical(0), fGapLogical(0),
  fWorldPhysical(0), fCalorPhysical(0), fLayerPhysical(0), fGapPhysical(0),
  fDetectorMessenger(0), fVerboseLevel(1)
{
  fLayerThickness = fTotalThickness / fNumberOfLayers;
  fCalName = "Calor";
  fDetectorMessenger = new GB03DetectorMessenger(this);
}

~GB03DetectorConstruction()

GB03DetectorConstruction::~GB03DetectorConstruction ( )

virtual

Definition at line 80 of file GB03DetectorConstruction.cc.

{ delete fDetectorMessenger;}

Member Function Documentation

Construct()

G4VPhysicalVolume * GB03DetectorConstruction::Construct ( )

virtual

Definition at line 85 of file GB03DetectorConstruction.cc.

{
  if(!fConstructed)
  {
    fConstructed = true;
    DefineMaterials();
    SetupGeometry();
  }
  if (GetVerboseLevel()>0)
  { PrintCalorParameters(); }
 
  return fWorldPhysical;
}

ConstructSDandField()

void GB03DetectorConstruction::ConstructSDandField ( )

virtual

Definition at line 101 of file GB03DetectorConstruction.cc.

{
  if(!fConstructedSDandField)
  {
    fConstructedSDandField = true;
    SetupDetectors();
    SetupBiasing();
  }
}

PrintCalorParameters()

void GB03DetectorConstruction::PrintCalorParameters

(

)

const

Definition at line 321 of file GB03DetectorConstruction.cc.

{
  G4cout 
    << "--------------------------------------------------------" << G4endl;
  G4cout 
    << " Absorber is made of " << fAbsorberMaterial->GetName() << G4endl
    << " Gap is made of " << fGapMaterial->GetName() << G4endl
    << "--------------------------------------------------------" << G4endl;
}

SetAbsorberMaterial()

void GB03DetectorConstruction::SetAbsorberMaterial

(

G4String

materialChoice

)

Definition at line 333 of file GB03DetectorConstruction.cc.

{
  // search the material by its name   
  G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice);     
  if(pttoMaterial)
  {
    fAbsorberMaterial = pttoMaterial;
    if(fConstructed) 
    {
      fCalorLogical->SetMaterial(fAbsorberMaterial);
      fLayerLogical->SetMaterial(fAbsorberMaterial);
    }
    G4RunManager::GetRunManager()->GeometryHasBeenModified();
    if (GetVerboseLevel()>1) {
      PrintCalorParameters();
    }
  }
  else
  { 
    G4cerr 
      << materialChoice << " is not defined. - Command is ignored." << G4endl; 
  }
}

GetAbsorberMaterial()

G4String GB03DetectorConstruction::GetAbsorberMaterial

(

)

const

Definition at line 359 of file GB03DetectorConstruction.cc.

{ return fAbsorberMaterial->GetName(); }

SetGapMaterial()

void GB03DetectorConstruction::SetGapMaterial

(

G4String

materialChoice

)

Definition at line 364 of file GB03DetectorConstruction.cc.

{
  // search the material by its name 
  G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice);  
  if(pttoMaterial)
  {
    fGapMaterial = pttoMaterial;
    if(fConstructed) 
    { fGapLogical->SetMaterial(fGapMaterial); }
    G4RunManager::GetRunManager()->GeometryHasBeenModified();
    if (GetVerboseLevel()>1) {
      PrintCalorParameters();
    }
  }
  else
  { 
    G4cerr 
      << materialChoice << " is not defined. - Command is ignored." << G4endl; 
  }
}

GetGapMaterial()

G4String GB03DetectorConstruction::GetGapMaterial

(

)

const

Definition at line 387 of file GB03DetectorConstruction.cc.

{ return fGapMaterial->GetName(); }

SetNumberOfLayers()

void GB03DetectorConstruction::SetNumberOfLayers

(

G4int

nl

)

Definition at line 392 of file GB03DetectorConstruction.cc.

{
  fNumberOfLayers = nl;
  fLayerThickness = fTotalThickness/fNumberOfLayers;
  if(!fConstructed) return;
 
  fLayerSolid->SetZHalfLength(fLayerThickness/2.);
  fGapSolid->SetZHalfLength(fLayerThickness/4.);
 
  fCalorLogical->RemoveDaughter(fLayerPhysical);
  delete fLayerPhysical;
  fLayerPhysical 
      = new G4PVReplica(fCalName+"_Layer",fLayerLogical,fCalorLogical,
                        kZAxis,fNumberOfLayers,fLayerThickness);
  fGapPhysical->SetTranslation(G4ThreeVector(0.,0.,fLayerThickness/4.));
  
  G4RunManager::GetRunManager()->GeometryHasBeenModified();
}

GetNumberOfLayers()

static G4int GB03DetectorConstruction::GetNumberOfLayers ( )

inlinestatic

Definition at line 58 of file GB03DetectorConstruction.hh.

{ return fNumberOfLayers; }

GetVerboseLevel()

G4int GB03DetectorConstruction::GetVerboseLevel ( ) const

inline

Definition at line 60 of file GB03DetectorConstruction.hh.

{ return  fVerboseLevel; }

SetVerboseLevel()

void GB03DetectorConstruction::SetVerboseLevel ( G4int  val )

inline

Definition at line 61 of file GB03DetectorConstruction.hh.

{ fVerboseLevel = val; }

DefineMaterials()

void GB03DetectorConstruction::DefineMaterials ( )

private

Definition at line 113 of file GB03DetectorConstruction.cc.

{ 
  G4String name, symbol;             //a=mass of a mole;
  G4double a, z, density;            //z=mean number of protons;  
  G4int iz;                          //iz=number of protons  in an isotope; 
  G4int n;                           // n=number of nucleons in an isotope;
 
  G4int ncomponents, natoms;
  G4double abundance, fractionmass;
  G4double temperature, pressure;
 
  //
  // define Elements
  //
 
  a = 1.01*g/mole;
  G4Element* H  = new G4Element(name="Hydrogen",symbol="H" , z= 1., a);
 
  a = 12.01*g/mole;
  G4Element* C  = new G4Element(name="Carbon"  ,symbol="C" , z= 6., a);
 
  a = 14.01*g/mole;
  G4Element* N  = new G4Element(name="Nitrogen",symbol="N" , z= 7., a);
 
  a = 16.00*g/mole;
  G4Element* O  = new G4Element(name="Oxygen"  ,symbol="O" , z= 8., a);
 
  //
  // define an Element from isotopes, by relative abundance 
  //
 
  G4Isotope* U5 = new G4Isotope(name="U235", iz=92, n=235, a=235.01*g/mole);
  G4Isotope* U8 = new G4Isotope(name="U238", iz=92, n=238, a=238.03*g/mole);
 
  G4Element* U  = new G4Element(name="enriched Uranium",symbol="U",ncomponents=2);
  U->AddIsotope(U5, abundance= 90.*perCent);
  U->AddIsotope(U8, abundance= 10.*perCent);
 
  //
  // define simple materials
  //
 
  new G4Material(name="Aluminium", z=13., a=26.98*g/mole, density=2.700*g/cm3);
  new G4Material(name="Silicon", z=14., a= 28.09*g/mole, density= 2.33*g/cm3);
  new G4Material(name="Iron", z=26., a=55.85*g/mole, density=7.87*g/cm3);
  new G4Material(name="ArgonGas",z=18., a= 39.95*g/mole, density=1.782*mg/cm3);
  new G4Material(name="He", z=2., a=4.0*g/mole, density=0.1786e-03*g/cm3);
 
  density = 1.390*g/cm3;
  a = 39.95*g/mole;
  new G4Material(name="liquidArgon", z=18., a, density);
 
  density = 11.35*g/cm3;
  a = 207.19*g/mole;
  G4Material* Pb = new G4Material(name="Lead"     , z=82., a, density);
 
  //
  // define a material from elements.   case 1: chemical molecule
  //
 
  density = 1.000*g/cm3;
  G4Material* H2O = new G4Material(name="Water", density, ncomponents=2);
  H2O->AddElement(H, natoms=2);
  H2O->AddElement(O, natoms=1);
 
  density = 1.032*g/cm3;
  G4Material* Sci = new G4Material(name="Scintillator", density, ncomponents=2);
  Sci->AddElement(C, natoms=9);
  Sci->AddElement(H, natoms=10);
 
  //
  // define a material from elements.   case 2: mixture by fractional mass
  //
 
  density = 1.290*mg/cm3;
  G4Material* Air = new G4Material(name="Air"  , density, ncomponents=2);
  Air->AddElement(N, fractionmass=0.7);
  Air->AddElement(O, fractionmass=0.3);
 
  //
  // examples of vacuum
  //
 
  density     = universe_mean_density;
  pressure    = 3.e-18*pascal;
  temperature = 2.73*kelvin;
  G4Material* Vacuum = new G4Material(name="Galactic", z=1., a=1.01*g/mole,
                                    density,kStateGas,temperature,pressure);
 
  if (GetVerboseLevel()>1) {
    G4cout << *(G4Material::GetMaterialTable()) << G4endl;
  }
 
  //default materials of the calorimeter
  fWorldMaterial    = Vacuum;
  fAbsorberMaterial = Pb;
  fGapMaterial      = Sci;
}

SetupGeometry()

void GB03DetectorConstruction::SetupGeometry ( )

private

Definition at line 214 of file GB03DetectorConstruction.cc.

{
  //     
  // World
  //
  G4VSolid* worldSolid = new G4Box("World",2.*m,2.*m,fTotalThickness*2.);
  fWorldLogical = new G4LogicalVolume(worldSolid,fWorldMaterial,"World");
  fWorldPhysical = new G4PVPlacement(0,G4ThreeVector(),fWorldLogical,"World",
                        0,false,0);
  
  //                               
  // Calorimeter
  //  
  G4VSolid* calorSolid = new G4Box("Calor",0.5*m,0.5*m,fTotalThickness/2.);
  fCalorLogical = new G4LogicalVolume(calorSolid,fAbsorberMaterial,fCalName);
  fCalorPhysical = new G4PVPlacement(0, G4ThreeVector(0.,0.,0.),
                       fCalorLogical,fCalName,fWorldLogical,false,0);
 
  //                                 
  // Layers --- as absorbers
  //
  fLayerSolid = new G4Box("Layer",0.5*m,0.5*m,fLayerThickness/2.);
  fLayerLogical
      = new G4LogicalVolume(fLayerSolid,fAbsorberMaterial,fCalName+"_LayerLog");
  fLayerPhysical
      = new G4PVReplica(fCalName+"_Layer",fLayerLogical,fCalorLogical,
                        kZAxis,fNumberOfLayers,fLayerThickness);
   
  //
  // Gap
  //
  fGapSolid = new G4Box("Gap",0.5*m,0.5*m,fLayerThickness/4.);
  fGapLogical = new G4LogicalVolume(fGapSolid,fGapMaterial,fCalName+"_Gap");
  fGapPhysical = new G4PVPlacement(0,G4ThreeVector(0.,0.,fLayerThickness/4.),
                fGapLogical,fCalName+"_gap",fLayerLogical,false,0);
 
  //                                        
  // Visualization attributes
  //
  fWorldLogical->SetVisAttributes(G4VisAttributes::GetInvisible());
  G4VisAttributes* simpleBoxVisAtt= new G4VisAttributes(G4Colour(1.0,1.0,1.0));
  simpleBoxVisAtt->SetVisibility(true);
  fCalorLogical->SetVisAttributes(simpleBoxVisAtt);
  fLayerLogical->SetVisAttributes(simpleBoxVisAtt);
  fGapLogical->SetVisAttributes(simpleBoxVisAtt);
  
}

SetupDetectors()

void GB03DetectorConstruction::SetupDetectors ( )

private

Definition at line 264 of file GB03DetectorConstruction.cc.

{
  G4SDManager::GetSDMpointer()->SetVerboseLevel(1);
  G4String filterName;
 
  G4SDNeutralFilter* neutralFilter
    = new G4SDNeutralFilter(filterName="neutralFilter");
  G4SDChargedFilter* chargedFilter
    = new G4SDChargedFilter(filterName="chargedFilter");
 
  for(G4int j=0;j<2;j++)
  {
    // Loop counter j = 0 : absorber
    //                = 1 : gap
    G4String detName = fCalName;
    if(j==0)
    { detName += "_abs"; }
    else
    { detName += "_gap"; }
    G4MultiFunctionalDetector* det = new G4MultiFunctionalDetector(detName);
    G4SDManager::GetSDMpointer()->AddNewDetector(det);
    // The second argument in each primitive means the "level" of geometrical 
    // hierarchy, the copy number of that level is used as the key of the 
    // G4THitsMap.
    // For absorber (j = 0), the copy number of its own physical volume is used.
    // For gap (j = 1), the copy number of its mother physical volume is used, 
    // since there is only one physical volume of gap is placed with respect 
    // to its mother.
    G4VPrimitiveScorer* primitive;
    primitive = new G4PSEnergyDeposit("eDep",j);
    det->RegisterPrimitive(primitive);
    primitive = new G4PSFlatSurfaceFlux("nNeutral",1,j);
    primitive->SetFilter(neutralFilter);
    det->RegisterPrimitive(primitive);
    primitive = new G4PSFlatSurfaceFlux("nCharged",1,j);
    primitive->SetFilter(chargedFilter);
    det->RegisterPrimitive(primitive);
 
    if(j==0)
    { SetSensitiveDetector(fLayerLogical, det); }
    else
    { SetSensitiveDetector(fGapLogical, det);}
   
  }
  G4SDManager::GetSDMpointer()->SetVerboseLevel(0);
}

SetupBiasing()

void GB03DetectorConstruction::SetupBiasing ( )

private

Definition at line 313 of file GB03DetectorConstruction.cc.

{
  GB03BOptrGeometryBasedBiasing* biasingOperator = new GB03BOptrGeometryBasedBiasing();
  biasingOperator->AttachTo(fLayerLogical);
}

Member Data Documentation

fNumberOfLayers

G4int GB03DetectorConstruction::fNumberOfLayers = 40

staticprivate

Definition at line 70 of file GB03DetectorConstruction.hh.

fTotalThickness

G4double GB03DetectorConstruction::fTotalThickness

private

Definition at line 72 of file GB03DetectorConstruction.hh.

fLayerThickness

G4double GB03DetectorConstruction::fLayerThickness

private

total thinkness of one calorimeter

Definition at line 73 of file GB03DetectorConstruction.hh.

fConstructed

G4bool GB03DetectorConstruction::fConstructed

private

= fTotalThickness / fNumberOfLayers

Definition at line 75 of file GB03DetectorConstruction.hh.

fConstructedSDandField

G4ThreadLocal G4bool GB03DetectorConstruction::fConstructedSDandField = false

staticprivate

Definition at line 76 of file GB03DetectorConstruction.hh.

fCalName

G4String GB03DetectorConstruction::fCalName

private

Definition at line 78 of file GB03DetectorConstruction.hh.

fWorldMaterial

G4Material* GB03DetectorConstruction::fWorldMaterial

private

Definition at line 80 of file GB03DetectorConstruction.hh.

fAbsorberMaterial

G4Material* GB03DetectorConstruction::fAbsorberMaterial

private

Definition at line 81 of file GB03DetectorConstruction.hh.

fGapMaterial

G4Material* GB03DetectorConstruction::fGapMaterial

private

Definition at line 82 of file GB03DetectorConstruction.hh.

fLayerSolid

G4Box* GB03DetectorConstruction::fLayerSolid

private

Definition at line 84 of file GB03DetectorConstruction.hh.

fGapSolid

G4Box* GB03DetectorConstruction::fGapSolid

private

Definition at line 85 of file GB03DetectorConstruction.hh.

fWorldLogical

G4LogicalVolume* GB03DetectorConstruction::fWorldLogical

private

Definition at line 87 of file GB03DetectorConstruction.hh.

fCalorLogical

G4LogicalVolume* GB03DetectorConstruction::fCalorLogical

private

Definition at line 88 of file GB03DetectorConstruction.hh.

fLayerLogical

G4LogicalVolume* GB03DetectorConstruction::fLayerLogical

private

Definition at line 89 of file GB03DetectorConstruction.hh.

fGapLogical

G4LogicalVolume* GB03DetectorConstruction::fGapLogical

private

Definition at line 90 of file GB03DetectorConstruction.hh.

fWorldPhysical

G4VPhysicalVolume* GB03DetectorConstruction::fWorldPhysical

private

Definition at line 92 of file GB03DetectorConstruction.hh.

fCalorPhysical

G4VPhysicalVolume* GB03DetectorConstruction::fCalorPhysical

private

Definition at line 93 of file GB03DetectorConstruction.hh.

fLayerPhysical

G4PVReplica* GB03DetectorConstruction::fLayerPhysical

private

Definition at line 94 of file GB03DetectorConstruction.hh.

fGapPhysical

G4VPhysicalVolume* GB03DetectorConstruction::fGapPhysical

private

Definition at line 95 of file GB03DetectorConstruction.hh.

fDetectorMessenger

GB03DetectorMessenger* GB03DetectorConstruction::fDetectorMessenger

private

Definition at line 97 of file GB03DetectorConstruction.hh.

fVerboseLevel

G4int GB03DetectorConstruction::fVerboseLevel

private

Definition at line 99 of file GB03DetectorConstruction.hh.

---

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

• .doxygen/Doxymodules_biasing.h
• extended/biasing/GB03/include/GB03DetectorConstruction.hh
• extended/biasing/GB03/src/GB03DetectorConstruction.cc