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

#include "io/MJLogger.hh"
#include "database/MJDatabase.hh"
#include "database/MJDatabasePsqlImpl.hh"
#include "database/MJDatabaseUtil.hh"
static void
exit_nicely(PGconn *conn)
{
        PQfinish(conn);
        exit(1);
}
//---------------------------------------------------------------------------//

#include "database/MJDatabaseCloverDetector.hh"      //Present MJ Class Headers 

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

MJDatabaseCloverDetector::MJDatabaseCloverDetector()
{

}

MJDatabaseCloverDetector::MJDatabaseCloverDetector(string serialNo)
  :fSerialNo( serialNo )
{
  //construct material here
}
MJDatabaseCloverDetector::MJDatabaseCloverDetector(const MJDatabaseCloverDetector & other)
{

}
MJDatabaseCloverDetector::~MJDatabaseCloverDetector()
{
}
MJDatabaseCloverDetector* 
MJDatabaseCloverDetector::createDetector(PGresult *res, PGconn * conn)
{
  MJDatabaseUtil dbutil;
  for (int i = 0; i < PQntuples(res); i++)
    {
      fSerialNo = PQgetvalue(res, i, 0);
      fSpacerWidth =dbutil.c_Float(PQgetvalue(res, i, 1));
      string spacerMat = PQgetvalue(res, i, 2);
      fSpacerBGModel = PQgetvalue(res, i, 3);
      fOpVoltage =dbutil.c_Float( PQgetvalue(res, i, 4));
      fCryoEndGap = dbutil.c_Float(PQgetvalue(res, i, 5));
      fCryoThinWallThickness = dbutil.c_Float(PQgetvalue(res, i, 6));
      fCryoThinWallLength = dbutil.c_Float(PQgetvalue(res, i, 7));
      fCryoThickWallThickness = dbutil.c_Float(PQgetvalue(res, i, 8));
      fCryoLength = dbutil.c_Float(PQgetvalue(res, i, 9));
      fCryoWidth = dbutil.c_Float(PQgetvalue(res, i, 10));
      fCryoCornerRadius = dbutil.c_Float(PQgetvalue(res, i, 11));
      string cryoMat = PQgetvalue(res, i, 12);
      fCryoBGModel = PQgetvalue(res, i, 13);
      fColdplateThickness =dbutil.c_Float(PQgetvalue(res, i, 14));
      fColdplateGap =dbutil.c_Float( PQgetvalue(res, i, 15));
      string coldPMat = PQgetvalue(res, i, 16);
      fColdplateBGModel = PQgetvalue(res, i, 17);

      //create the material objects
      fSpacerMaterial = *(MJDatabase::GetMaterial(spacerMat));
      fCryoMaterial = *(MJDatabase::GetMaterial(cryoMat));
      fColdplateMaterial = *(MJDatabase::GetMaterial(coldPMat));

      int len, num;
      //get Pl resolution at different energies
      len = PQgetlength(res,i, 18);
      float * ener ;
      if( len > 0) dbutil.fillArray(ener, PQgetvalue(res, i, 18));
      len = PQgetlength(res,i, 19);
      float * resol;
      if( len > 0)  dbutil.fillArray(resol, PQgetvalue(res, i, 19));
      //create the vector of resolutions
      num =dbutil.c_Int(PQgetvalue(res, i, 18)+2*sizeof(int));  
      for (int j=0; j< num; j++){
        fResPl.push_back(DetectorRes(ener[j], resol[j]));
      }
     //get Pm resolution at different energies 
      len = PQgetlength(res,i, 20);
      if( len > 0)  dbutil.fillArray(ener, PQgetvalue(res, i, 20));
      len = PQgetlength(res,i, 21);
      if( len > 0)  dbutil.fillArray(resol, PQgetvalue(res, i, 21));
      //create the vector of resolutions
      num =dbutil.c_Int(PQgetvalue(res, i, 20)+2*sizeof(int));  
      for (int j=0; j< num; j++){
        fResPm.push_back(DetectorRes(ener[j], resol[j]));
      }
     //get Pr resolution at different energies
      len = PQgetlength(res,i, 22);
      if( len > 0)  dbutil.fillArray(ener, PQgetvalue(res, i, 22));
      len = PQgetlength(res,i, 23);
      if( len > 0)  dbutil.fillArray(resol, PQgetvalue(res, i, 23));
      //create the vector of resolutions
      num =dbutil.c_Int(PQgetvalue(res, i, 22)+2*sizeof(int));  
      for (int j=0; j< num; j++){
        fResPr.push_back(DetectorRes(ener[j], resol[j]));
      }
      //get the parts
      getParts(conn);
    }
  return this;
}
void
MJDatabaseCloverDetector::getParts(PGconn * conn){
  PGresult *res;
  string cmd = "DECLARE Aportal CURSOR FOR select * from CloverPart_tot" ;
  cmd += " where parentid = '" + fSerialNo + "'";
  res = PQexec(conn,  cmd.c_str());
  if (PQresultStatus(res) != PGRES_COMMAND_OK)
    {
      MJLog(error) << "DECLARE CURSOR failed: " << PQerrorMessage(conn) << endlog;
      PQclear(res);
      exit_nicely(conn);
    }
  PQclear(res);
  res = PQexec(conn, "FETCH ALL in Aportal");

  //check to see if we found the object in database
  if( PQntuples(res) == 0 ) {
    MJLog(error) << "Could not find " << fSerialNo << " in the database " << endlog;
    return ;
  }

  for (int i = 0; i < PQntuples(res); i++){
    fCrystalSerialNo[i] = PQgetvalue(res,i,0);
    fCrystals.push_back(*(MJDatabase::GetCloverCrystal(fCrystalSerialNo[i])));
  }
}
ostream & operator << (ostream& os, MJDatabaseCloverDetector& det)
{
  os << '\n'
     <<"****** Clover Detector: " << det.GetDetectorSerialNo() << '\n'
     <<"Cryo material      : " << det.GetCryoMaterialName() << '\n'
     <<"Coldplate material : " << det.GetColdplateMaterialName() << '\n'
     <<"Spacer material    : " << det.GetSpacerMaterialName() << '\n'
     <<"Crystal serial no. : { " ;
  for (int i=0; i<4; i++)
    os << det.fCrystalSerialNo[i] << " ";
  os << "}" << '\n';

  return os;

}

---