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


//#include "MJStuff.hh"                   //additional MJ Headers

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

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

MJDatabaseClover4x50x80DetCalib::MJDatabaseClover4x50x80DetCalib()
{

}
MJDatabaseClover4x50x80DetCalib::MJDatabaseClover4x50x80DetCalib(const MJDatabaseClover4x50x80DetCalib & other)
{

}

MJDatabaseClover4x50x80DetCalib::~MJDatabaseClover4x50x80DetCalib()
{
}
MJDatabaseClover4x50x80DetCalib*
MJDatabaseClover4x50x80DetCalib::createCalib(PGresult* res)
{
  MJDatabaseUtil dbutil;
  for(int i=0; i<PQntuples(res); i++)
    {
      string name = PQgetvalue(res, i, 0);
      fDetName = PQgetvalue(res, i, 1);
      float Acoef = dbutil.c_Float(PQgetvalue(res, i, 2));
      float Bcoef = dbutil.c_Float(PQgetvalue(res, i, 3));
      float Ccoef = dbutil.c_Float(PQgetvalue(res, i, 4));
      float stdvA = dbutil.c_Float(PQgetvalue(res, i, 5));
      float stdvB = dbutil.c_Float(PQgetvalue(res, i, 6));
      
      //get list of energies at which measurements were made
      int len = PQgetlength(res,i, 7);
      float *ener;
      if( len > 0) dbutil.fillArray(ener, PQgetvalue(res, i, 7));
      float * fwhm;
      len = PQgetlength(res,i, 8);
      if( len > 0) dbutil.fillArray(fwhm, PQgetvalue(res, i, 8));
      int* channel;
      len = PQgetlength(res,i, 9);
      if( len > 0) dbutil.fillArray(channel, PQgetvalue(res, i, 9));

      //this need fixing. this returns long long and should be 
      //converted inot the time 
      //char* bt = dbutil.c_Time(PQgetvalue(res, i, 10));
      //char* et = dbutil.c_Time(PQgetvalue(res, i, 11));

      char* bt = PQgetvalue(res, i, 10);
      char* et = PQgetvalue(res, i, 11);

      //create a calibration object and add it to the vector
      fCalib.insert(make_pair(name, CloverDet4x50x80Calib
                              ((double)Acoef, 
                               (double)Bcoef, (double)Ccoef,
                               (double)stdvA, (double)stdvB, bt, et,
                               (double*)ener, (double*)fwhm, channel)));
      fName.push_back(name);


    }
  return this;
}
ostream& operator << (ostream& os, const CloverDet4x50x80Calib& calib)
{
  os << calib.GetA() << "  " <<calib.GetB() << "  " << calib.GetC() 
     << "  " << calib.GetSDEVA()
     << "  "<< calib.GetSDEVB() << "  " << calib.GetBeginTime() 
     << "  " << calib.GetEndTime() 
     << "  " << calib.fenergy[0] << " , " << calib.fenergy[1] 
     << " , " << calib.fenergy[2] 
     << "  " << calib.fFWHM[0] << " , " << calib.fFWHM[1] 
     << " , " << calib.fFWHM[2] 
     << "  " << calib.fchannelNo[0] << " , " << calib.fchannelNo[1] 
     << " , " << calib.fchannelNo[2]
     << "    " << calib.GetBeginTime() 
     << "    " << calib.GetEndTime() 
    ;
  return os;
}

ostream & operator << (ostream& os, MJDatabaseClover4x50x80DetCalib& calib)
{
  os << '\n'
     <<"****** Clover Detector calibration for: " << calib.GetDetName() << '\n'
     <<" E1   :\n " << calib.GetE1()<< '\n'     
     <<" E2   :\n " << calib.GetE2() << '\n'
     <<" E3   :\n " << calib.GetE3()<< '\n'
     <<" E4   :\n " << calib.GetE4()<< '\n'
     <<" Pl   :\n " << calib.GetPl()<< '\n'
     <<" Pm   :\n " << calib.GetPm()<< '\n'
     <<" Pr   :\n " << calib.GetPr()<< '\n'
     <<'\n';

  return os;
}

---