MJGeometryDetectorConstruction.cc

Go to the documentation of this file.

//---------------------------------------------------------------------------//
//bb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nu//
//                                                                           //
//                         MAJORANA Simulation                               //
//                                                                           //
//      This code implementation is the intellectual property of the         //
//      MAJORANA Collaboration. It is based on Geant4, an intellectual       //
//      property of the RD44 GEANT4 collaboration.                           //
//                                                                           //
//                        *********************                              //
//                                                                           //
//    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.    //
//    By copying, distributing or modifying the Program (or any work based   //
//    on on the Program) you indicate your acceptance of this statement,     //
//    and all its terms.                                                     //
//                                                                           //
//bb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nu//
//---------------------------------------------------------------------------//
//                                                          
// $Id: MJGeometryDetectorConstruction.cc,v 1.5 2004/11/26 09:53:18 pandola Exp $ 
//      
// CLASS IMPLEMENTATION:  MJGeometryDetectorConstruction.cc
//
//---------------------------------------------------------------------------//
// 
//---------------------------------------------------------------------------//
//---------------------------------------------------------------------------//
//

#include <vector>

//  GEANT4 headers
#include "globals.hh"
#include "G4Box.hh"
#include "G4LogicalVolume.hh"
#include "G4Material.hh"
#include "G4VPhysicalVolume.hh"
#include "G4PVPlacement.hh"
#include "G4RunManager.hh"

//  MJ headers
#include "geometry/MJGeometryCloverDetector.hh"
#include "geometry/MJGeometryDetectorMessenger.hh"
#include "geometry/MJGeometryDetector.hh"
#include "geometry/MJGeometryCloverInNaIBarrel.hh"
#include "geometry/MJGeometrySolidBlock.hh"
#include "database/MJDatabase.hh"
#include "database/MJDatabaseMaterial.hh"
#include "materials/MJMaterial.hh"
#include "materials/MJGerdaLocalMaterialTable.hh"
#include "io/MJLogger.hh"

//  Detector headers (also fall under the title "MJ headers")
#include "geometry/MJGeometryDetectorConstruction.hh"
#include "gerdageometry/MaGeGeometryGerda.hh"
#include "gerdageometry/MaGeGeometryGerdaHole.hh"
#include "gerdageometry/MaGeGeometryGerdaCrystal.hh"
#include "gerdageometry/MaGeGeometryGerdaCrystalHole.hh"


//Luciano
#include "G4MaterialTable.hh"

using namespace std;

//---------------------------------------------------------------------------//

MJGeometryDetectorConstruction::MJGeometryDetectorConstruction():
  fMJMaterials(0), fNumberOfMaterials(0), theDetector(0),
  fUseDatabase(true)
{
        detectorMessenger = new MJGeometryDetectorMessenger(this);
        localTable = new  MJGerdaLocalMaterialTable;
}

//---------------------------------------------------------------------------//

MJGeometryDetectorConstruction::~MJGeometryDetectorConstruction() 
{
  while(fNumberOfMaterials > 0)
    delete fMJMaterials[--fNumberOfMaterials];
  delete []fMJMaterials;
  delete detectorMessenger;
  delete theDetector;
  delete localTable;
}

//---------------------------------------------------------------------------//

G4VPhysicalVolume *MJGeometryDetectorConstruction::Construct() 
{
  if (fUseDatabase) 
    {
      BuildAllMaterialsInDatabase();
    }
  else 
    {
      localTable->CreateMaterials();
    }
  return ConstructDetector();
}

//---------------------------------------------------------------------------//

G4VPhysicalVolume *MJGeometryDetectorConstruction::ConstructDetector() 
{       
  // RH: Construct Detector here, *after* materials are defined.
  theDetector->ConstructDetector();

        //  Create the physical world
        solidWorld = new G4Box( "solidWorld", 10*m, 10*m, 10*m );
        logicWorld = new G4LogicalVolume( solidWorld, G4Material::GetMaterial("Air"), "World" );
        physiWorld = new G4PVPlacement( 0, G4ThreeVector(), "World", logicWorld, 0, false, 0 );
        
        //the name of the Detector is "Detector" by default
        theDetectorPhysical = new G4PVPlacement( 0, G4ThreeVector(),theDetector->GetDetectorName(), 
                                                 theDetector->GetDetectorLogical(), 
                                                 physiWorld, false, 0 );
        
        //  Always return the physical world
        return physiWorld;
}

//---------------------------------------------------------------------------//

void MJGeometryDetectorConstruction::BuildAllMaterialsInDatabase() 
{
  vector<string> materialNames = MJDatabase::GetMaterialNames();
  vector<string>::iterator namesIterator = materialNames.begin();
  fNumberOfMaterials = materialNames.size();
  if(MJLogger::GetSeverity() == MJLogger::debugging) {
    MJLog(trace) << "Building materials from database: " << endlog;
    while(namesIterator != materialNames.end()) {
      G4cout << *namesIterator << G4endl;
      namesIterator++;
    }      
  }
  namesIterator = materialNames.begin();
  fMJMaterials = new MJMaterial*[fNumberOfMaterials];
  if(!fMJMaterials) {
    MJLog(error) << "Could not allocate memory for materials!" << endlog;
    MJLog(fatal) << endlog;
  }
  G4int i;
  for(i = 0; i < fNumberOfMaterials; i++) { 
    // Horrible hack for now. G4 does not like Hg196. To be fixed...
    // RH. 7/15/2004 
    if(*namesIterator != "Mercury") {
      MJLog(debugging) << "Building material : " << *namesIterator << endlog;
      fMJMaterials[i] = new MJMaterial((*namesIterator).c_str(), true);
    }
    namesIterator++;
  }
  MJLog(routine) << "Built all materials in database : " 
                 << fNumberOfMaterials << " materials built." << endlog;
}

//---------------------------------------------------------------------------//

void MJGeometryDetectorConstruction::SetDetector( G4String detectorType ) 
{
  if(!theDetector) {
        if( detectorType == "clover" )
          theDetector = new MJGeometryCloverDetector( "LANL1" );
        else if (detectorType == "cloverinNaIbarrel")
          theDetector = new MJGeometryCloverInNaIBarrel("0000");
        else if (detectorType == "solidBlock")
          theDetector = new MJGeometrySolidBlock();
        else if (detectorType == "GerdaArray")
          theDetector = new MaGeGeometryGerda();
        else if (detectorType == "GerdaCrystal")
          theDetector = new MaGeGeometryGerdaCrystal();
        else if (detectorType == "GerdaCrystalHole")
          theDetector = new MaGeGeometryGerdaCrystalHole();
        else if (detectorType == "GerdaArrayHole")
          theDetector = new MaGeGeometryGerdaHole();
        else
                G4cout << "Couldn't find detector of type \"" << detectorType 
                       << "\"!" << G4endl;

        MJLog(routine) << "Using the " << detectorType << " detector."<<endlog;
  } else {
    MJLog(error) << "Error: Detector already defined." << endlog;
  }
}

//---------------------------------------------------------------------------//

---