Delete construction.cc

GEANT4/construction.cc

-#include "construction.hh"
-
-MyDetectorConstruction::MyDetectorConstruction()
-{
-    nCols = 10;
-    nRows = 10;
-
-    fMessenger = new G4GenericMessenger(this, "/detector/", "Detector Construction");
-
-    fMessenger->DeclareProperty("nCols", nCols, "Number of cols");
-    fMessenger->DeclareProperty("nRows", nRows, "Number of rows");
-    fMessenger->DeclareProperty("atmosphere", isAtmosphere, "Construct Atmosphere");
-
-    DefineMaterials();
-
-    isAtmosphere = true;
-
-    xWorld = 2300.*m;
-    yWorld = 2300.*m;
-    zWorld = 1500.*m;
-}
-
-MyDetectorConstruction::~MyDetectorConstruction()
-{}
-
-void MyDetectorConstruction::DefineMaterials()
-{
-    G4NistManager *nist = G4NistManager::Instance();
-
-    worldMat = nist->FindOrBuildMaterial("G4_AIR");
-
-    G4double energy[2] = {1.239841939*eV/0.2, 1.239841939*eV/0.9};
-    G4double rindexWorld[2] = {1.0, 1.0};
-
-    G4MaterialPropertiesTable *mptWorld = new G4MaterialPropertiesTable();
-    mptWorld->AddProperty("RINDEX", energy, rindexWorld, 2);
-
-
-    G4double fraction[2] = {1.0, 1.0};
-
-    worldMat->SetMaterialPropertiesTable(mptWorld);
-
-    G4double reflectivity[2] = {1.0, 1.0};
-
-    G4double density0 = 1.29*kg/m3;
-    G4double aN = 14.01*g/mole;
-    G4double aO = 16.00*g/mole;
-
-    G4Element *elN = new G4Element("Nitrogen", "N", 7, aN);
-    G4Element *elO = new G4Element("Oxygen", "O", 8, aO);
-
-    G4double f = 3;
-    G4double R = 8.3144626181532;
-    G4double g0 = 9.81;
-    G4double kappa = (f+2)/f;
-    G4double T = 293.15;
-    G4double M = (0.3*14.01 + 0.7*16.0)/1000.;
-
-    for(G4int i = 0; i < 10; i++)
-    {
-        std::stringstream stri;
-        stri << i;
-        G4double h = 2.3e3/10.*i;
-        G4double density = density0*pow((1-(kappa-1)/kappa*M*g0*h/(R*T)),(1/(kappa-1)));
-        Air[i] = new G4Material("G4_AIR_"+stri.str(), density, 2);
-        Air[i]->AddElement(elN, 70*perCent);
-        Air[i]->AddElement(elO, 30*perCent);
-    }
-}
-
-void MyDetectorConstruction::ConstructAtmosphere()
-{
-    solidAir = new G4Box("solidAir", xWorld, yWorld, zWorld/10.);
-
-    for(G4int i = 0; i < 10; i++)
-    {
-        logicAir[i] = new G4LogicalVolume(solidAir, Air[i], "logicAir");
-
-        physAir[i] = new G4PVPlacement(0, G4ThreeVector(0, 0, zWorld/10.*2*i - zWorld + zWorld/10.), logicAir[i], "physAir", logicWorld, false, i, true);
-    }
-}
-
-G4VPhysicalVolume *MyDetectorConstruction::Construct()
-{
-    solidWorld = new G4Box("solidWorld", xWorld, yWorld, zWorld);
-
-    logicWorld = new G4LogicalVolume(solidWorld, worldMat, "logicWorld");
-
-    physWorld = new G4PVPlacement(0, G4ThreeVector(0., 0., 0.), logicWorld, "physWorld", 0, false, 0, true);
-
-    if(isAtmosphere)
-        ConstructAtmosphere();
-
-    return physWorld;
-}