exrdmDetectorConstruction Class Reference

#include <exrdmDetectorConstruction.hh>

Inheritance diagram for exrdmDetectorConstruction:

Public Member Functions
	exrdmDetectorConstruction ()
virtual	~exrdmDetectorConstruction ()
virtual G4VPhysicalVolume *	Construct ()
const G4VPhysicalVolume *	GetDetector ()
G4double	GetDetectoFullLength ()
G4double	GetTargetFullLength ()
G4double	GetWorldFullLength ()
G4double	GetDetectorThickness ()
G4double	GetTargetRadius ()
G4double	GetWorldRadius ()
void	SetTargetMaterial (G4String)
void	SetDetectorMaterial (G4String)
void	SetTargetRadius (G4double value)
void	SetDetectorThickness (G4double value)
void	SetTargetLength (G4double value)
void	SetDetectorLength (G4double value)
Public Member Functions inherited from G4VUserDetectorConstruction
	G4VUserDetectorConstruction ()
virtual	~G4VUserDetectorConstruction ()
virtual void	ConstructSDandField ()
virtual void	CloneSD ()
virtual void	CloneF ()
void	RegisterParallelWorld (G4VUserParallelWorld *)
G4int	ConstructParallelGeometries ()
void	ConstructParallelSD ()
G4int	GetNumberOfParallelWorld () const
G4VUserParallelWorld *	GetParallelWorld (G4int i) const

Additional Inherited Members
Protected Member Functions inherited from G4VUserDetectorConstruction
void	SetSensitiveDetector (const G4String &logVolName, G4VSensitiveDetector *aSD, G4bool multi=false)
void	SetSensitiveDetector (G4LogicalVolume *logVol, G4VSensitiveDetector *aSD)

Detailed Description

Definition at line 48 of file exrdmDetectorConstruction.hh.

Constructor & Destructor Documentation

exrdmDetectorConstruction::exrdmDetectorConstruction

(

)

Definition at line 54 of file exrdmDetectorConstruction.cc.

:G4VUserDetectorConstruction(),
 fSolidWorld(0),  fLogicWorld(0),  fPhysiWorld(0),
 fSolidTarget(0), fLogicTarget(0), fPhysiTarget(0), 
 fSolidDetector(0),fLogicDetector(0),fPhysiDetector(0),
 fDetectorMessenger(0),
 fMaterialsManager(0),
 fDefaultMater(0),fTargetMater(0),fDetectorMater(0),
 fTargetLength (1.*cm), fTargetRadius(0.5*cm),
 fDetectorLength(5.0 * cm), fDetectorThickness(2.0 * cm),
 fWorldLength (std::max(fTargetLength,fDetectorLength)),
 fWorldRadius (fTargetRadius + fDetectorThickness),
 fTargetRegion(0), fDetectorRegion(0)
{
  fDetectorMessenger = new exrdmDetectorMessenger(this);
  DefineMaterials();
}

exrdmDetectorConstruction::~exrdmDetectorConstruction ( )

virtual

Definition at line 74 of file exrdmDetectorConstruction.cc.

{
  delete fDetectorMessenger;             
}

Member Function Documentation

G4VPhysicalVolume * exrdmDetectorConstruction::Construct ( void  )

virtual

Implements G4VUserDetectorConstruction.

Definition at line 103 of file exrdmDetectorConstruction.cc.

References G4Region::AddRootLogicalVolume(), G4UImanager::ApplyCommand(), G4UImanager::GetUIpointer(), G4VisAttributes::Invisible, G4INCL::Math::max(), G4LogicalVolume::SetVisAttributes(), and python.hepunit::twopi.

{
//--------- Definitions of Solids, Logical Volumes, Physical Volumes ---------
  //--------- Sizes of the principal geometrical components (solids)  ---------

  fWorldLength = std::max(fTargetLength,fDetectorLength);
  fWorldRadius = fTargetRadius + fDetectorThickness;
   
  //------------------------------ 
  // World
  //------------------------------ 

 fSolidWorld= new G4Tubs("world",0.,fWorldRadius,fWorldLength/2.,0.,twopi);
 fLogicWorld= new G4LogicalVolume( fSolidWorld, fDefaultMater, "World", 0, 0, 0);
  
  //  Must place the World Physical volume unrotated at (0,0,0).
  // 
  fPhysiWorld = new G4PVPlacement(0,              // no rotation
                                 G4ThreeVector(), // at (0,0,0)
                                 fLogicWorld,     // its logical volume
                                 "World",         // its name
                                 0,               // its mother  volume
                                 false,           // no boolean operations
                                 0);              // no field specific to volume
                                 
  //------------------------------ 
  // Target
  //------------------------------
  
  G4ThreeVector positionTarget = G4ThreeVector(0,0,0);
   
  fSolidTarget = new G4Tubs("target",0.,fTargetRadius,fTargetLength/2.,0.,twopi);
  fLogicTarget = new G4LogicalVolume(fSolidTarget,fTargetMater,"Target",0,0,0);
  fPhysiTarget = new G4PVPlacement(0,              // no rotation
                                  positionTarget,  // at (x,y,z)
                                  fLogicTarget,    // its logical volume                                  
                                  "Target",        // its name
                                  fLogicWorld,     // its mother  volume
                                  false,           // no boolean operations
                                  0);              // no particular field 

  //  G4cout << "Target is a cylinder with rdius of " << targetradius/cm << " cm of " 
  //       << fTargetMater->GetName() << G4endl;

  //------------------------------ 
  // Detector
  //------------------------------
  
  G4ThreeVector positionDetector = G4ThreeVector(0,0,0);
  
  fSolidDetector = new G4Tubs("detector",fTargetRadius,fWorldRadius,
                                                          fDetectorLength/2.,0.,twopi);
  fLogicDetector = new G4LogicalVolume(fSolidDetector ,fDetectorMater,
                                                                      "Detector",0,0,0);
  fPhysiDetector = new G4PVPlacement(0,             // no rotation
                                  positionDetector, // at (x,y,z)
                                  fLogicDetector,   // its logical volume                                  
                                  "Detector",       // its name
                                  fLogicWorld,      // its mother  volume
                                  false,            // no boolean operations
                                  0);               // no particular field 

  //------------------------------------------------ 
  // Sensitive detectors
  //------------------------------------------------ 

  //  G4SDManager* SDman = G4SDManager::GetSDMpointer();

  // G4String detectortargetSDname = "exrdm/DetectorTargetSD";
  // exrdmDetectorSD* aDetectorSD = new exrdmDetectorSD( detectorTargetSDname );
  // SDman->AddNewDetector( aDetectorSD );
  //fLogicTarget->SetSensitiveDetector( aDetectorSD );
  // fLogicDetector->SetSensitiveDetector( aDetectorSD );
  //
  //-------------------------------------------------
  // regions
  //
  //  if(fTargetRegion) delete fTargetRegion;
  // if(fDetectorRegion) delete fDetectorRegion;
  fTargetRegion = new G4Region("Target");
  fDetectorRegion   = new G4Region("Detector");
  fTargetRegion->AddRootLogicalVolume(fLogicTarget);
  fDetectorRegion->AddRootLogicalVolume(fLogicDetector);

  //--------- Visualization attributes -------------------------------
  fLogicWorld->SetVisAttributes(G4VisAttributes::Invisible);
  G4VisAttributes* TargetVisAtt= new G4VisAttributes(G4Colour(1.0,1.0,1.0));
  fLogicTarget ->SetVisAttributes(TargetVisAtt);
  G4VisAttributes* DetectorVisAtt= new G4VisAttributes(G4Colour(1.0,1.0,.0));
  fLogicDetector->SetVisAttributes(DetectorVisAtt);


  //------------ set the incident position ------

 // get the pointer to the User Interface manager 
    
  G4UImanager* UI = G4UImanager::GetUIpointer();  
  //      UI->ApplyCommand("/run/verbose 1");
  //      UI->ApplyCommand("/event/verbose 2");
  //      UI->ApplyCommand("/tracking/verbose 1");  

  G4double zpos = -fWorldLength/2.;
  G4String command = "/gps/pos/centre ";
  std::ostringstream os;
  os << zpos ; 
  G4String xs = os.str();
  UI->ApplyCommand(command+"0. 0. "+xs+" mm");
  UI->ApplyCommand("/gps/pos/type Point");
  command = "/gps/position ";
  //  UI->ApplyCommand(command+"0. 0. "+xs+" mm");
  UI->ApplyCommand("/gps/particle proton");
  UI->ApplyCommand("/gps/direction 0 0 1");
  UI->ApplyCommand("/gps/energy 100 MeV");
  //       
  
  return fPhysiWorld;
}

G4double exrdmDetectorConstruction::GetDetectoFullLength ( )

inline

Definition at line 62 of file exrdmDetectorConstruction.hh.

{return fDetectorLength;};

const G4VPhysicalVolume* exrdmDetectorConstruction::GetDetector ( )

inline

Definition at line 60 of file exrdmDetectorConstruction.hh.

{return fPhysiDetector;};

G4double exrdmDetectorConstruction::GetDetectorThickness ( )

inline

Definition at line 65 of file exrdmDetectorConstruction.hh.

{return fDetectorThickness;};

G4double exrdmDetectorConstruction::GetTargetFullLength ( )

inline

Definition at line 63 of file exrdmDetectorConstruction.hh.

{return fTargetLength;};

G4double exrdmDetectorConstruction::GetTargetRadius ( )

inline

Definition at line 66 of file exrdmDetectorConstruction.hh.

{return fTargetRadius;};

G4double exrdmDetectorConstruction::GetWorldFullLength ( )

inline

Definition at line 64 of file exrdmDetectorConstruction.hh.

{return fWorldLength;}; 

G4double exrdmDetectorConstruction::GetWorldRadius ( )

inline

Definition at line 67 of file exrdmDetectorConstruction.hh.

{return fWorldRadius;}; 

void exrdmDetectorConstruction::SetDetectorLength ( G4double  value )

inline

Definition at line 75 of file exrdmDetectorConstruction.hh.

Referenced by exrdmDetectorMessenger::SetNewValue().

{ fDetectorLength = value;};  

void exrdmDetectorConstruction::SetDetectorMaterial

(

G4String

materialName

)

Definition at line 238 of file exrdmDetectorConstruction.cc.

References python.hepunit::cm, G4cout, G4endl, G4Material::GetMaterial(), and G4LogicalVolume::SetMaterial().

Referenced by exrdmDetectorMessenger::SetNewValue().

{
  // search the material by its name 
  G4Material* pttoMaterial = G4Material::GetMaterial(materialName);  
  if (pttoMaterial)
     {fDetectorMater = pttoMaterial;
      if (fLogicDetector) fLogicDetector->SetMaterial(pttoMaterial); 
      G4cout << "\n----> The Deetctor has been changed to" << fDetectorLength/cm
                     << " cm of "
             << materialName << G4endl;
     }             
}

void exrdmDetectorConstruction::SetDetectorThickness ( G4double  value )

inline

Definition at line 73 of file exrdmDetectorConstruction.hh.

Referenced by exrdmDetectorMessenger::SetNewValue().

{ fDetectorThickness = value;};  

void exrdmDetectorConstruction::SetTargetLength ( G4double  value )

inline

Definition at line 74 of file exrdmDetectorConstruction.hh.

Referenced by exrdmDetectorMessenger::SetNewValue().

{ fTargetLength = value; };

void exrdmDetectorConstruction::SetTargetMaterial

(

G4String

materialName

)

Definition at line 223 of file exrdmDetectorConstruction.cc.

References python.hepunit::cm, G4cout, G4endl, G4Material::GetMaterial(), and G4LogicalVolume::SetMaterial().

Referenced by exrdmDetectorMessenger::SetNewValue().

{
  // search the material by its name 
  G4Material* pttoMaterial = G4Material::GetMaterial(materialName);  
  if (pttoMaterial)
     {fTargetMater = pttoMaterial;
      if (fLogicTarget) fLogicTarget->SetMaterial(pttoMaterial); 
      G4cout << "\n----> The target has been changed to " << fTargetLength/cm
                     << " cm of "
             << materialName << G4endl;
     }             
}

void exrdmDetectorConstruction::SetTargetRadius ( G4double  value )

inline

Definition at line 72 of file exrdmDetectorConstruction.hh.

Referenced by exrdmDetectorMessenger::SetNewValue().

{ fTargetRadius = value; };

---

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

• exrdmDetectorConstruction.hh
• exrdmDetectorConstruction.cc