Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/g3tog4/clGeometry/src/G3toG4DetectorConstruction.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  1 //
  2 // ********************************************************************
  3 // * License and Disclaimer                                           *
  4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.                             *
 10 // *                                                                  *
 11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                                                  *
 18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // ********************************************************************
 25 //
 26 /// \file g3tog4/clGeometry/src/G3toG4DetectorConstruction.cc
 27 /// \brief Implementation of the G3toG4DetectorConstruction class
 28 //
 29 //
 30 //
 31 //--------------------------------------------------------------------------
 32 // G3toG4DetectorConstruction. Most the work is Done in
 33 // G4BuildGeom, which returns a G4LogicalVolume*, a pointer to the
 34 // top-level logiical volume in the detector defined by the call List file
 35 // inFile
 36 //--------------------------------------------------------------------------
 37 
 38 #include "G3toG4DetectorConstruction.hh"
 39 
 40 #include "G4Box.hh"
 41 #include "G4Colour.hh"
 42 #include "G4Material.hh"
 43 #include "G4SystemOfUnits.hh"
 44 #include "G4VisAttributes.hh"
 45 #include "G4ios.hh"
 46 
 47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 48 
 49 G3toG4DetectorConstruction::G3toG4DetectorConstruction(G4String inFile)
 50   : G4VUserDetectorConstruction(), fInFile(inFile)
 51 {
 52   G4cout << "Instantiated G3toG4DetectorConstruction using call list file \"" << fInFile << "\""
 53          << G4endl;
 54 }
 55 
 56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 57 
 58 G3toG4DetectorConstruction::~G3toG4DetectorConstruction()
 59 {
 60   // G4cout << "Deleted G3toG4DetectorConstruction..." << G4endl;
 61 }
 62 
 63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 64 
 65 G4VPhysicalVolume* G3toG4DetectorConstruction::Construct()
 66 {
 67   G4LogicalVolume* lv = G4BuildGeom(fInFile);
 68   // G4LogicalVolume* lv = SimpleConstruct();
 69   if (lv) {
 70     G4VPhysicalVolume* pv = new G4PVPlacement(0, G4ThreeVector(), lv, lv->GetName(), 0, false, 0);
 71     G4cout << "Top-level G3toG4 logical volume " << lv->GetName() << " "
 72            << *(lv->GetVisAttributes()) << G4endl;
 73     return pv;
 74   }
 75 
 76   G4cerr << "creation of logical mother failed !!!" << G4endl;
 77   return 0;
 78 }
 79 
 80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 81 
 82 G4LogicalVolume* G3toG4DetectorConstruction::SimpleConstruct()
 83 {
 84   G4String name, symbol;  // a=mass of a mole;
 85   G4double a, z, density, fractionmass;  // z=mean number of protons;
 86   G4int ncomponents;  // iz=number of protons  in an isotope;
 87   // n=number of nucleons in an isotope;
 88 
 89   a = 14.01 * g / mole;
 90   G4Element* eN = new G4Element(name = "Nitrogen", symbol = "N", z = 7., a);
 91 
 92   a = 16.00 * g / mole;
 93   G4Element* eO = new G4Element(name = "Oxygen", symbol = "O", z = 8., a);
 94 
 95   //
 96   // define a material from elements.   case 2: mixture by fractional mass
 97   //
 98 
 99   density = 1.290 * mg / cm3;
100   G4Material* air = new G4Material(name = "Air", density, ncomponents = 2);
101   air->AddElement(eN, fractionmass = 0.7);
102   air->AddElement(eO, fractionmass = 0.3);
103   G4VSolid* mother = new G4Box("TestMother",  // its name
104                                100 * cm, 100 * cm, 100 * cm);  // its size
105 
106   G4VSolid* daughter = new G4Box("TestDaughter", 50 * cm, 20 * cm, 10 * cm);
107 
108   G4LogicalVolume* logicMother = new G4LogicalVolume(mother,  // its solid
109                                                      air,  // its material
110                                                      "LTestMother");  // its name
111 
112   G4LogicalVolume* logicDaughter = new G4LogicalVolume(daughter,  // its solid
113                                                        air,  // its material
114                                                        "LTestDaughter");
115 
116   new G4PVPlacement(0, G4ThreeVector(), logicDaughter, "PTestDaughter", logicMother, false, 0);
117   //
118   // Visualization attributes
119   //
120 
121   logicMother->SetVisAttributes(G4VisAttributes::GetInvisible());
122   G4VisAttributes* daughterVisAtt = new G4VisAttributes(G4Colour(1.0, 1.0, 1.0));
123   daughterVisAtt->SetVisibility(true);
124   logicDaughter->SetVisAttributes(daughterVisAtt);
125   return logicMother;
126 }
127 
128 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
129