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

#include <vector>

#include "globals.hh"
#include "G4Element.hh"
#include "G4Isotope.hh"

#include "database/MJDatabase.hh"
#include "database/MJDatabaseElement.hh"
#include "io/MJLogger.hh"
#include "materials/MJMaterialIsotope.hh"

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

#include "materials/MJMaterialElement.hh"

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

MJMaterialElement::MJMaterialElement():
 fExists(false), fG4Element(0)
{
  fMJIsotopes.clear();
}

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

MJMaterialElement::MJMaterialElement(G4String name, G4bool registerG4):
  fName(name), fG4Element(0)
{
  fMJIsotopes.clear();
  CreateFromDatabase();
  if(registerG4) RegisterWithG4();
}

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

MJMaterialElement::MJMaterialElement(const MJMaterialElement & other)
{
  fExists = other.fExists;
  fName = other.fName;
  fSymbol = other.fSymbol;
  fNoOfIsotopes = other.fNoOfIsotopes;
  fZ = other.fZ;
  fG4Element = other.fG4Element;
  fIsoName = other.fIsoName;
  fIsoSymb = other.fIsoSymb;
  fIsoAbundance = other.fIsoAbundance;
  fMJIsotopes = other.fMJIsotopes;
}

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

MJMaterialElement::~MJMaterialElement()
{;}

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

void MJMaterialElement::CreateFromDatabase()
{
  MJDatabaseElement* theDBElement = MJDatabase::GetElement(fName);
  if(theDBElement) { 
    fName = theDBElement->GetName();
    fSymbol = theDBElement->GetSymbol();
    fNoOfIsotopes = theDBElement->GetNoOfIsotopes();
    fZ = theDBElement->GetZ();
    fIsoName = theDBElement->GetIsoNames();
    fIsoAbundance = theDBElement->GetIsoAbundance();
    fExists = true;
    MJLog(trace) << fName << " MJ element created from database." << endlog;
  } else {
    fExists = false;
    MJLog(error) << fName << " element not in database." << endlog;
  }
}

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

void MJMaterialElement::RegisterWithG4()
{
  if(fExists) {
    if((fG4Element = G4Element::GetElement(fName))) {
      MJLog(debugging)<<fName<<" element already registered with G4."<<endlog;
    } else {
      // Register new element with G4
      fG4Element = new G4Element(fName, fSymbol, fNoOfIsotopes);

      // Setup loop to add isotopes.
      G4Isotope *g4iso;
      vector<string>::iterator componentIterator = fIsoName.begin();
      vector<G4double>::iterator abundanceIterator = fIsoAbundance.begin();

      // Define and add required isotopes.
      while(componentIterator != fIsoName.end() && 
            abundanceIterator != fIsoAbundance.end()) {
        if(!(g4iso = G4Isotope::GetIsotope(*componentIterator))) {
          MJMaterialIsotope *mjiso = new 
            MJMaterialIsotope(*componentIterator, true);
          fMJIsotopes.push_back(*mjiso);
          if(!(g4iso = G4Isotope::GetIsotope(*componentIterator))){
            MJLog(error) << "Could not register " << *componentIterator 
                         <<" with G4."<<endlog;
            MJLog(fatal) << endlog;
          }
        } // if
        fG4Element->AddIsotope(g4iso, *abundanceIterator * 100.0 * perCent);
        componentIterator++;
        abundanceIterator++;
      } // while
      MJLog(trace) << "Element " << fName << " created in G4." << endlog;
    } // if
  } else {
    MJLog(error) << "Cannot register element " << fName << " with G4 before "
                 << "defining it." << endlog;
  }  
}

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

---