// // ******************************************************************** // * License and Disclaimer * // * * // * The Geant4 software is copyright of the Copyright Holders of * // * the Geant4 Collaboration. It is provided under the terms and * // * conditions of the Geant4 Software License, included in the file * // * LICENSE and available at http://cern.ch/geant4/license . These * // * include a list of copyright holders. * // * * // * Neither the authors of this software system, nor their employing * // * institutes,nor the agencies providing financial support for this * // * work make any representation or warranty, express or implied, * // * regarding this software system or assume any liability for its * // * use. Please see the license in the file LICENSE and URL above * // * for the full disclaimer and the limitation of liability. * // * * // * This code implementation is the result of the scientific and * // * technical work of the GEANT4 collaboration. * // * By using, copying, modifying or distributing the software (or * // * any work based on the software) you agree to acknowledge its * // * use in resulting scientific publications, and indicate your * // * acceptance of all terms of the Geant4 Software license. * // ******************************************************************** // // // $Id: ExDivDetectorConstruction.cc,v 1.3 2006/06/29 18:19:53 gunter Exp $ // GEANT4 tag $Name: geant4-09-04-beta-cand-01 $ // //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... #include "ExDivDetectorConstruction.hh" #include "G4Material.hh" #include "G4Box.hh" #include "G4VisAttributes.hh" #include "G4Colour.hh" #include "G4ios.hh" #include "ExDivTesterBox.hh" #include "ExDivTesterTubs.hh" #include "ExDivTesterCons.hh" #include "ExDivTesterTrd.hh" #include "ExDivTesterPara.hh" #include "ExDivTesterPolycone.hh" #include "ExDivTesterPolyhedra.hh" //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... ExDivDetectorConstruction:: ExDivDetectorConstruction( const G4String& solidTypeStr, const G4String& PVTypeStr, const G4String& PosTypeStr, const std::vector& extraPars ) : theSolidTypeStr( solidTypeStr ), thePVTypeStr( PVTypeStr ), thePosTypeStr( PosTypeStr ), theExtraPars( extraPars ) { } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... ExDivDetectorConstruction::~ExDivDetectorConstruction() { } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... G4VPhysicalVolume* ExDivDetectorConstruction::Construct() { theDivTester = CreateSolidTester(theSolidTypeStr, thePVTypeStr, thePosTypeStr, theExtraPars); //- SolidType soltype = getSolidType( theSolidTypeStr ); // PVType pvtype = pvDivision; G4VPhysicalVolume* myTopNode = theDivTester->BuildGeometry( ); theDivTester->GenerateScanPoints(); /* // Repeat tests but with full voxels G4GeometryManager::GetInstance()->OpenGeometry(); G4GeometryManager::GetInstance()->CloseGeometry(true); testG4Navigator1(myTopNode); testG4Navigator2(myTopNode); */ // theDivTester->PrintParentSolid( G4cout ); // theDivTester->PrintChildrenSolids( G4cout ); // G4GeometryManager::GetInstance()->OpenGeometry(); return myTopNode; } //-------------------------------------------------------------------------- ExVDivTester* ExDivDetectorConstruction::CreateSolidTester( const G4String& stype, const G4String& thePVTypeStr, const G4String& thePosTypeStr, std::vector& extraPars ) { PVType pvtype = getPVType( thePVTypeStr ); PlaceType postype = getPosType( thePosTypeStr ); ExVDivTester* theSolidTester = 0; if( stype == "box" ) { theSolidTester = new ExDivTesterBox( pvtype, postype, extraPars ); } else if( stype == "tubs" ) { theSolidTester = new ExDivTesterTubs( pvtype, postype, extraPars ); ExVDivTester::bDivCylindrical = 1; } else if( stype == "cons" ) { theSolidTester = new ExDivTesterCons( pvtype, postype, extraPars ); ExVDivTester::bDivCylindrical = 1; } else if( stype == "trd" ) { theSolidTester = new ExDivTesterTrd( pvtype, postype, extraPars ); } else if( stype == "para" ) { theSolidTester = new ExDivTesterPara( pvtype, postype, extraPars ); } else if( stype == "pcone" ) { theSolidTester = new ExDivTesterPolycone( pvtype, postype, extraPars ); ExVDivTester::bDivCylindrical = 1; } else if( stype == "phedra" ) { theSolidTester = new ExDivTesterPolyhedra( pvtype, postype, extraPars ); } else { G4cout << "ERROR - ExDivDetectorConstruction::CreateSolidTester()" << G4endl << " Only the following solid types are allowed:" << G4endl << " 'tubs', 'cons', 'trd', 'para', 'pcone', 'phedra'." << G4endl << " Wrong identifier: " << stype << G4endl; G4Exception("ExDivDetectorConstruction::CreateSolidTester()", "InvalidSetup", FatalException, "Unknown solid type."); } return theSolidTester; } //-------------------------------------------------------------------------- PVType ExDivDetectorConstruction::getPVType( const G4String& pvt ) { G4cout << pvt << G4endl; PVType vtype = pvPlacement; if( pvt == "division") { vtype = pvDivision; } else if( pvt == "replica") { vtype = pvReplica; } else { G4Exception("ExDivDetectorConstruction::getPVType()", "InvalidArgument", FatalException, "PV type can only be 'division' or 'replica' "); } return vtype; } //-------------------------------------------------------------------------- PlaceType ExDivDetectorConstruction::getPosType( const G4String& pos ) { G4cout << pos << G4endl; PlaceType ptype = pvNormal; if( pos == "normal") { ptype = pvNormal; } else if( pos == "reflected") { ptype = pvReflected; } else { G4Exception("ExDivDetectorConstruction::getPosType()", "InvalidArgument", FatalException, "The positioning type can only be 'normal' or 'reflected' "); } return ptype; }