BesEmcParameter.cc

Go to the documentation of this file.

//---------------------------------------------------------------------------//
//      BOOST --- BESIII Object_Oriented Simulation Tool                     //
//---------------------------------------------------------------------------//
/// Description:
// Author: He Miao
// Created: Nov,12 2004
// Modified:
// Comment:
//---------------------------------------------------------------------------//
//$ID: BesEmcParameter.hh
 
#include "EmcSim/BesEmcParameter.hh"
#include <fstream>
#include <strstream>
// #include "ReadBoostRoot.hh"
#include "EmcGeneralClass/EmcStructure.h"
using namespace std;
 
// Initialize static data member
BesEmcParameter* BesEmcParameter::fpInstance = 0;
 
BesEmcParameter::BesEmcParameter() {}
 
BesEmcParameter::~BesEmcParameter() {}
 
// static method
// Access to an instance
BesEmcParameter& BesEmcParameter::GetInstance() {
  if ( !Exist() )
  {
    fpInstance = new BesEmcParameter;
    fpInstance->ReadData();
  }
  return *fpInstance;
}
 
bool BesEmcParameter::Exist() { return fpInstance != 0; }
 
void BesEmcParameter::Kill() {
  if ( Exist() )
  {
    delete fpInstance;
    fpInstance = 0;
  }
}
 
void BesEmcParameter::ReadData() {
  G4String ParaPath = getenv( "EMCSIMROOT" );
  if ( !ParaPath )
  {
    G4cout << "BOOST environment not set!" << G4endl;
    exit( -1 );
  }
 
  G4String ParaPath1 = ParaPath;
 
  ParaPath += "/dat/BesEmc.txt";
  ifstream fin;
  fin.open( ParaPath );
  assert( fin );
 
  const int maxCharOfOneLine = 255;
  char      temp[maxCharOfOneLine], *p;
  int       lineNo = 0, inputNo = 0;
 
  while ( fin.peek() != EOF )
  {
    fin.getline( temp, maxCharOfOneLine );
    p = temp;
    lineNo++;
    while ( *p != '\0' )
    {
      if ( *p == '#' )
      {
        *p = '\0'; // delete the comments.
        break;
      }
      p++;
    }
    p = temp; // reset the pointer to the beginning of the string.
    while ( *p == ' ' || *p == '\t' ) p++;
    if ( *p == '\0' ) continue;
    inputNo++;
    switch ( inputNo )
    {
    case 1:
      istrstream( p ) >> array_size >> m_tau >> m_highRange >> m_midRange >> m_lowRange >>
          m_sampleTime >> m_bitNb >> m_photonsPerMeV >> m_nonuniformity >> m_peakTime >>
          m_timeOffset;
      break;
    case 2:
      istrstream( p ) >> WorldRmin1 >> WorldRmax1 >> WorldRmin2 >> WorldRmax2 >> WorldDz >>
          WorldZPosition >> CrystalLength >> CrystalLength1;
      break;
    case 3:
      istrstream( p ) >> cryNumInOneLayer[0] >> cryNumInOneLayer[1] >> cryNumInOneLayer[2] >>
          cryNumInOneLayer[3] >> cryNumInOneLayer[4] >> cryNumInOneLayer[5];
      break;
    case 4:
      istrstream( p ) >> pentaInOneSector[0] >> pentaInOneSector[1] >> pentaInOneSector[2] >>
          pentaInOneSector[3] >> pentaInOneSector[4];
      break;
    case 5: istrstream( p ) >> fTyvekThickness >> fAlThickness >> fMylarThickness; break;
    case 6:
      istrstream( p ) >> BSCRmin >> BSCDz >> BSCRmin1 >> BSCRmax1 >> BSCRmin2 >> BSCRmax2 >>
          BSCDz1;
      break;
    case 7: istrstream( p ) >> BSCAngleRotat >> BSCNbPhi >> BSCNbTheta; break;
    case 8:
      istrstream( p ) >> BSCYFront0 >> BSCYFront >> BSCYFront1 >> BSCPosition0 >> BSCPosition1;
      break;
    case 9:
      istrstream( p ) >> TaperRingDz >> TaperRingThickness1 >> TaperRingThickness2 >>
          TaperRingThickness3 >> TaperRingTheta >> TaperRingInnerLength >>
          TaperRingOuterLength;
    case 10:
      istrstream( p ) >> rearBoxLength >> rearBoxDz >> HangingPlateDz >> OCGirderAngle >>
          rearCasingThickness;
    case 11:
      istrstream( p ) >> orgGlassLengthX >> orgGlassLengthY >> orgGlassLengthZ >> PDLengthX >>
          PDLengthY >> PDLengthZ >> AlPlateDz >> PABoxDz >> PABoxThickness;
    case 12:
      istrstream( p ) >> cableDr >> waterPipeDr >> waterPipeThickness >> SPBarThickness >>
          SPBarThickness1 >> SPBarwidth >> EndRingDz >> EndRingDr >> EndRingRmin;
    default:;
    }
  }
  fin.close();
 
  ifstream fin1;
  ParaPath1 += "/dat/EmcLightOutput.dat";
  fin1.open( ParaPath1 );
  assert( fin1 );
  for ( G4int i = 0; i < 6240; i++ ) { fin1 >> lightOutput[i]; }
  fin1.close();
}
 
G4double BesEmcParameter::GetLightOutput( G4int partId, G4int theta, G4int phi ) {
  EmcStructure struc;
  struc.setEmcStruc();
  G4int index = struc.getGeomIndex( partId, theta, phi );
  if ( index >= 0 && index < 6240 ) { return lightOutput[index]; }
  else { return 1.; }
}