testDatabase.cc

Go to the documentation of this file.

#include <iostream>

using namespace std;
#include "database/MJDatabaseClover4x50x80DetCalib.hh"
#include "materials/MJMaterialElement.hh"
#include "database/MJDatabase.hh"
#include "database/MJDatabaseMaterial.hh"
#include "database/MJDatabaseMEGACrystal.hh"
#include "database/MJDatabaseSEGACrystal.hh"
#include "database/MJDatabaseElement.hh"
#include "database/MJDatabaseCloverCrystal.hh"
#include "database/MJDatabaseCloverDetector.hh"
#include "database/MJDatabaseIsotope.hh"
#include "database/MJDatabaseSimple800g.hh"
#include "materials/MJMaterial.hh"
#include "io/MJLogger.hh"
//      Geant4 package includes
//#include "G4RunManager.hh"
//#include "G4UImanager.hh"
//#include "G4UIterminal.hh"
//#include "G4UItcsh.hh"
//#include "G4Material.hh"
#include "G4Box.hh"
#include "globals.hh"

void g4test();
int main(){

  MJDatabaseCloverCrystal* mjclover = MJDatabase::GetCloverCrystal("73771");
  //    mjclover->print(cout);

  MJLog(trace) << *mjclover;  // << endlog;

  
  MJDatabaseElement* mjele1 = MJDatabase::GetElement("Cadmium");
  MJDatabaseElement* mjele2 = MJDatabase::GetElement("Mercury");
  MJLog(trace) << *mjele1;
  MJLog(trace) << *mjele2;

  MJDatabaseMaterial* mat1 = MJDatabase::GetMaterial("Delrin");
  MJLog(trace) << *mat1;

  //  MJMaterialElement* ele = new MJMaterialElement("Carbon");
  //MJLog(trace) << *ele;


  //MJMaterial* mat2 = new MJMaterial("Delrin");
  //MJLog(trace) << *mat2;

  MJDatabaseCloverDetector * det = MJDatabase::GetCloverDetector("LANL1");
  MJLog(trace) << *det;

  MJDatabaseIsotope* iso = MJDatabase::GetIsotope("Boron11");
  MJLog(trace) << *iso;

  MJDatabaseSimple800g* simple =  MJDatabase::GetSimple800g("01");
  MJLog(trace) << *simple;

  vector<MJDatabaseMaterial> matVec = MJDatabase::GetAllMaterials();
  vector<MJDatabaseMaterial>::iterator iter;
  MJLog(trace) << "------ Printing all materials --------------" << endlog;
  for(iter = matVec.begin(); iter != matVec.end(); iter++)
    {
      MJLog(trace) << *iter;
    }

  vector<string> nameVec = MJDatabase::GetMaterialNames();
  vector<string>::iterator stit;
  MJLog(trace) << "------ Printing all material names --------------" << endlog;
  for(stit = nameVec.begin(); stit != nameVec.end(); stit++)
    {
      MJLog(trace) << *stit << "  "<<endlog;;
      
      MJDatabaseMaterial* matt = MJDatabase::GetMaterial(*stit);
      MJLog(trace) << *matt;
   }

  
  /*
  MJDatabaseClover4x50x80DetCalib* calib1 = 
    MJDatabase::GetClover4x50x80Calib("LANL1", "-infinity", "infinity");
    MJLog(trace) << *calib1;
  MJDatabaseClover4x50x80DetCalib* calib2 = 
    MJDatabase::GetClover4x50x80Calib("LANL1", "Jan-1-2004");
    MJLog(trace) << *calib2;
    MJLog(trace) << " E1(A):"<< calib2->GetE1().GetA()<<endlog;
    MJLog(trace) << " E1(B):"<< calib2->GetE1().GetB()<<endlog;
    MJLog(trace) << " E1(C):"<< calib2->GetE1().GetC()<<endlog;

    MJLog(trace) << " E2(A):"<< calib2->GetE2().GetA()<<endlog;
    MJLog(trace) << " E2(B):"<< calib2->GetE2().GetB()<<endlog;
    MJLog(trace) << " E2(C):"<< calib2->GetE2().GetC()<<endlog;

    MJLog(trace) << " E3(A):"<< calib2->GetE3().GetA()<<endlog;
    MJLog(trace) << " E3(B):"<< calib2->GetE3().GetB()<<endlog;
    MJLog(trace) << " E3(C):"<< calib2->GetE3().GetC()<<endlog;
    MJLog(trace) << " Pr(A):"<< calib2->GetPr().GetA()<<endlog;
    MJLog(trace) << " Pr(B):"<< calib2->GetPr().GetB()<<endlog;
    MJLog(trace) << " Pr(C):"<< calib2->GetPr().GetC()<<endlog;
    MJLog(trace) << " energies:"<< (calib2->GetE2().GetEnergy())[0]<<endlog;
    MJLog(trace) << " energies:"<< (calib2->GetE2().GetEnergy())[1]<<endlog;

  MJDatabaseClover4x50x80DetCalib* calib3 = 
    MJDatabase::GetClover4x50x80Calib("LANL1");
    MJLog(trace) << *calib3;
  */
  return 0;
}
void g4test(){
    G4double  worldSizeX = 4.*m;
  G4double  worldSizeYZ = 4.*m;
    G4Box*  solidWorld = new G4Box( "World", worldSizeX/2, worldSizeYZ/2, worldSizeYZ/2 );

}

---