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

#include "io/MJLogger.hh"

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

#include "database/MJDatabaseMaterial.hh"      //Present MJ Class Headers 
#include "database/MJDatabaseUtil.hh"
//---------------------------------------------------------------------------//

MJDatabaseMaterial::MJDatabaseMaterial()
{
}
MJDatabaseMaterial::MJDatabaseMaterial(string name, double den, 
                                       double temp, double pres, int nc,
                                       string defbkg, string com)
  :fName(name),
   fDensity (den),
   fTemperature(temp),
   fPressure (pres),
   fNumOfComponents (nc),
   fDefaultBKG (defbkg),
   fComment (com)
{
}

MJDatabaseMaterial::MJDatabaseMaterial(const MJDatabaseMaterial & other)
{

  fName = other.fName;
  fDensity = other.fDensity;
  fTemperature = other.fTemperature;
  fPressure = other.fPressure;
  fNumOfComponents = other.fNumOfComponents;
  fDefaultBKG = other.fDefaultBKG;
  fComment = other.fComment;
  fComponents = other.fComponents;
  fNumAtoms = other.fNumAtoms;
  fFracMass = other.fFracMass;
}

MJDatabaseMaterial::~MJDatabaseMaterial()
{
}
MJDatabaseMaterial*
MJDatabaseMaterial::createMaterial(PGresult * res, int n){

  MJDatabaseUtil dbutil;
  // for (int i = 0; i < n; i++){
    fName = PQgetvalue(res, n, 0);
    fDensity = dbutil.c_Float(PQgetvalue(res, n, 1));
    fTemperature = dbutil.c_Float(PQgetvalue(res, n, 2));
    fPressure = dbutil.c_Float(PQgetvalue(res, n, 3));
    
    fNumOfComponents = dbutil.c_Int(PQgetvalue(res, n, 4));
    
    fDefaultBKG = PQgetvalue(res, n, 6);
    
    //get the time stamp
    
    //get the comments
   for (int i = 0; i < fNumOfComponents; i++){
   
    fComment = PQgetvalue(res, i+n, 8);
    fComponents.push_back(PQgetvalue(res, i+n, 9));
    
    fNumAtoms.push_back(dbutil.c_Int(PQgetvalue(res, i+n, 10)));
    fFracMass.push_back(dbutil.c_Float(PQgetvalue(res, i+n, 11)));
  }
  return this;
}
void 
MJDatabaseMaterial::print(ostream & os){
  
 os << '\n' 
    << "******* Material: " << GetName() << " ******" << '\n'
    << "Density: " << GetDensity() << '\n'
    << "Temperature: " << GetTemperature() <<'\n'
    << "Pressure:    " << GetPressure() << '\n'
    << "Number of components: " << GetNoComponents() << '\n'
    <<endlog;
 
  
}
ostream& 
operator << (ostream& os, MJDatabaseMaterial & mat){
  os << '\n' 
     << "******* Material: " << mat.GetName() << " ******" << '\n'
     << "Density     : " << mat.GetDensity() << '\n'
     << "Temperature : " << mat.GetTemperature() <<'\n'
     << "Pressure    : " << mat.GetPressure() << '\n'
     << "Number of components: " << mat.GetNoComponents() << '\n';
  vector<string>::iterator itr1;
  os << "Components  : { ";
  for (itr1 = mat.fComponents.begin(); itr1 != mat.fComponents.end(); itr1++)
    os << *itr1 << " ";
  os << "}\n";

  vector<int>::iterator itr2;
  os << "No of atoms  : { ";
  for (itr2 = mat.fNumAtoms.begin(); itr2 != mat.fNumAtoms.end(); itr2++)
    os << *itr2 << " ";
  os << "}\n";

  vector<double>::iterator itr3;
  os << "Fractional mass: { ";
  for (itr3 = mat.fFracMass.begin(); itr3 != mat.fFracMass.end(); itr3++)
    os << *itr3 << " ";
  os << "}\n";


  os<<endl;
 
  return os;
}

---