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

#include "database/MJDatabase.hh"
using namespace std;
//---------------------------------------------------------------------------//

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

MJDatabaseSimple800g::MJDatabaseSimple800g()
{

}

MJDatabaseSimple800g::MJDatabaseSimple800g(const MJDatabaseSimple800g & other)
  :fCrystalMat(other.fCrystalMat),
   fInnerContactMat(other.fInnerContactMat)
{
  fSerialNo = other.fSerialNo;
  fHeight = other.fHeight;
  fRadius = other.fRadius;
  fCoreRadius = other.fCoreRadius;
  fCornerR = other.fCornerR;
  fCoreDepth = other.fCoreDepth;
  fInnerContactW = other.fInnerContactW;
  fInnerDeadThickness = other.fInnerDeadThickness;
  fOuterDeadThickness = other.fOuterDeadThickness;
  fDiodType = other.fDiodType;
  fDepletionVolt = other.fDepletionVolt;
  fResolution = other.fResolution;
  fEfficiencies= other.fEfficiencies;
  fInnerBkgModel = other.fInnerBkgModel;
  fBkgModel = other.fBkgModel;
}

MJDatabaseSimple800g::~MJDatabaseSimple800g()
{
}
MJDatabaseSimple800g*
MJDatabaseSimple800g::createCrystal(PGresult * res){

  MJDatabaseUtil dbutil;
  for (int i = 0; i < PQntuples(res); i++)
    {
      fSerialNo = PQgetvalue(res, i, 0);
      string matName =  PQgetvalue(res, i, 1);
      fHeight = dbutil.c_Float(PQgetvalue(res, i, 2));
      fRadius = dbutil.c_Float(PQgetvalue(res, i, 3));

      fBkgModel = PQgetvalue(res, i, 4);
      fCoreRadius = dbutil.c_Float(PQgetvalue(res, i, 5));
      fCoreDepth = dbutil.c_Float(PQgetvalue(res, i, 6));
      fCornerR = dbutil.c_Float(PQgetvalue(res, i, 7));

      fInnerContactW = dbutil.c_Float(PQgetvalue(res, i, 8));
 
      string innerMat = PQgetvalue(res, i, 9);
      fInnerDeadThickness = dbutil.c_Float(PQgetvalue(res, i, 10));
      fOuterDeadThickness = dbutil.c_Float(PQgetvalue(res, i, 11));
      fInnerBkgModel = PQgetvalue(res, i, 12);
      fDiodType = *(char*)PQgetvalue(res, i, 13);
      fDepletionVolt = dbutil.c_Float(PQgetvalue(res, i, 14));

      //get the matrials
      fCrystalMat = *(MJDatabase::GetMaterial(matName));
      fInnerContactMat = *(MJDatabase::GetMaterial(innerMat));

      //get resolution at different energies
      int len = PQgetlength(res,i, 15);
      float * ener ;
      if( len > 0) dbutil.fillArray(ener, PQgetvalue(res, i, 15));
      len = PQgetlength(res,i, 16);
      float * resol;
      if( len > 0) dbutil.fillArray(resol, PQgetvalue(res, i, 16));
      //create the vector of resolutions
      int num =dbutil.getNoEle(PQgetvalue(res, i, 15));  
      for (int j=0; j< num; j++){
        fResolution.push_back(CrystalRes(ener[j], resol[j]));
      }

      //get efficiencies at different energies

      len = PQgetlength(res,i, 17);
      if( len > 0) dbutil.fillArray(ener, PQgetvalue(res, i, 17));
      len = PQgetlength(res,i, 18);
      float * eff;
      if( len > 0) dbutil.fillArray(eff, PQgetvalue(res, i, 18));
      //create the vector of efficiencies
      num =dbutil.getNoEle(PQgetvalue(res, i, 17));
      for (int j=0; j< num; j++){
        fEfficiencies.push_back(CrystalEff(ener[j], eff[j]));
      }

      fcomment = PQgetvalue(res, i, 19);


    }
  return this;
} 
void 
MJDatabaseSimple800g::print(ostream &os){

  os << '\n' 
     << "******* Simple 800g Crystal: " << GetCrystalSerialNo() 
     << " ******" << '\n'
     << "Radius       : " << GetCrystalRadius() << '\n'
     << "Height       : " << GetCrystalHeight() <<'\n'
     << "Core Radius  : " << GetCoreRadius() << '\n'
     << "Corner Radius: " << GetCornerRadius() << '\n'
     << "Core Depth   : " << GetCoreDepth() << '\n';

  vector<CrystalRes>::iterator itr;
  os<< "Resolution   : " ;
  for(itr = fResolution.begin(); itr != fResolution.end(); itr++)
    os <<  itr->resolution() << "(" << itr->energy() << ") "; 
  os << '\n';

  vector<CrystalEff>::iterator it;
  os<< "Efficiencies   : " ;
  for(it = fEfficiencies.begin(); it != fEfficiencies.end(); it++)
    os <<  it->efficiency () << "(" << it->energy() << ") "; 
  os <<endl;

  fCrystalMat.print(os);
  fInnerContactMat.print(os);

}
ostream& 
operator << (ostream& os, MJDatabaseSimple800g& crys){
 
  os << '\n' 
     << "******* Simple 800g Crystal: " << crys.GetCrystalSerialNo() 
     << " ******" << '\n'
     << "Radius       : " << crys.GetCrystalRadius() << '\n'
     << "Height       : " << crys.GetCrystalHeight() <<'\n'
     << "Core Radius  : " << crys.GetCoreRadius() << '\n'
     << "Corner Radius: " << crys.GetCornerRadius() << '\n'
     << "Core Depth   : " << crys.GetCoreDepth() << '\n';

  vector<CrystalRes>::iterator itr;
  os<< "Resolution   : " ;
  for(itr = crys.fResolution.begin(); itr != crys.fResolution.end(); itr++)
    os <<  itr->resolution() << "(" << itr->energy() << ") "; 
  os << '\n';

  vector<CrystalEff>::iterator it;
  os<< "Efficiencies   : " ;
  for(it = crys.fEfficiencies.begin(); it != crys.fEfficiencies.end(); it++)
    os <<  it->efficiency () << "(" << it->energy() << ") "; 
  os <<endl;

  os << crys.fCrystalMat;
  os << crys.fInnerContactMat;
 
  return os;
}

---