MJMaterial.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: MJMaterial.cc,v 1.2 2004/11/09 13:42:39 xliu Exp $ 
//      
// CLASS IMPLEMENTATION:  MJMaterial.cc
//
//---------------------------------------------------------------------------//
// 
//---------------------------------------------------------------------------//
//---------------------------------------------------------------------------//
//

#include "globals.hh"
#include "G4ios.hh"
#include "G4Material.hh"               
#include "G4Element.hh"
  
#include "database/MJDatabase.hh"
#include "materials/MJMaterialElement.hh"
#include "io/MJLogger.hh"

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

#include "materials/MJMaterial.hh"

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

MJMaterial::MJMaterial():
 fExists(false), fG4Material(0)
{
  fMJElements.clear();
}

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

MJMaterial::MJMaterial(G4String name, G4bool registerG4):
  fMaterialName(name), fG4Material(0)
{
  fMJElements.clear();
  CreateFromDatabase();
  if(registerG4) RegisterWithG4();
}

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

MJMaterial::~MJMaterial() {;}

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

void MJMaterial::CreateFromDatabase() 
{
    // Get Material information from database
    MJDatabaseMaterial *theDBMaterial = MJDatabase::GetMaterial(fMaterialName);
    if(theDBMaterial) {
      fDensity = theDBMaterial->GetDensity() * g / cm3;
      fTemperature = theDBMaterial->GetTemperature() * kelvin;
      fPressure = theDBMaterial->GetPressure() * pascal;
      fNumComponents = theDBMaterial->GetNoComponents();
      fNumOfAtoms = theDBMaterial->GetNoAtoms();
      fComponents = theDBMaterial->GetComponents();
      fFracMass = theDBMaterial->GetFracMass(); // CAUTION ... UNDEFINED UNITS.
      fDefaultBKG = theDBMaterial->GetBackground();
      fMaterialComment = theDBMaterial->GetComment();
      fExists = true;
      MJLog(trace) << "Material "<<fMaterialName<<" created from database."
                   << endlog;
    } else {
      fExists = false;
      MJLog(fatal) << "Material "<<fMaterialName<<" not found in database."
                   << endlog;
    }
 }

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

void MJMaterial::RegisterWithG4()
{
  // Note: This method can currently only build a material from elements, 
  // not other materials. 
  // This will be fixed in the next version.
  //
  if(fExists) {
    if((fG4Material = G4Material::GetMaterial(fMaterialName))) {
      MJLog(debugging)<<fMaterialName<<" already registered with G4."<<endlog;
    } else {
      // Register new materials with G4
      G4State state = kStateUndefined;
      fG4Material = new G4Material(fMaterialName, fDensity, fNumComponents,
                                   state, fTemperature, fPressure);

      // Setup loop to add elements.
      vector<string>::iterator  componentsIterator = fComponents.begin();
      vector<G4int>::iterator  NatomsIterator = fNumOfAtoms.begin();
      G4Element *g4ele;
      G4String componentName;

      // Define and add required elements
      while (componentsIterator != fComponents.end() && 
             NatomsIterator != fNumOfAtoms.end()) {
        componentName = *componentsIterator;
        if(!(g4ele = G4Element::GetElement(componentName))) {
          MJMaterialElement *mjElement =
            new MJMaterialElement(componentName, true);
          fMJElements.push_back(*mjElement);
          if(!(g4ele = G4Element::GetElement(componentName))) {
            MJLog(fatal) << "Could not register " << *componentName 
                         << " with G4." << endlog;
          }
        }
        fG4Material->AddElement(g4ele, *NatomsIterator);
        componentsIterator++; 
        NatomsIterator++;
      }
      MJLog(trace) << "Material " << fMaterialName << " created in G4 :"
      << "Density : " << fDensity << " " << fNumComponents << endlog;
    }
  } else {
    MJLog(error) << "Cannot register material " << fMaterialName << " with G4"
                 << " before defining it." << endlog;
  }
}
//---------------------------------------------------------------------------//
//---------------------------------------------------------------------------//

---