BesSCMParameter.cc

Go to the documentation of this file.

// BOOST --- BESIII Object_Oriented Simulation Tool               //
// Description:
// Author: Caogf
// Created: June, 2007
// Modified:
// Comment:
//---------------------------------------------------------------------//
//$ID: BesSCMParameter.cc
 
#include "BesSCMParameter.hh"
#include "SimUtil/ReadBoostRoot.hh"
#include <fstream>
#include <strstream>
using namespace std;
 
BesSCMParameter::BesSCMParameter() {}
BesSCMParameter::~BesSCMParameter() {}
void BesSCMParameter::ReadData() {
  G4String GeometryPath = getenv( "BESSIMROOT" );
  if ( !GeometryPath )
  {
    G4cout << "BOOST environment not set!" << G4endl;
    exit( -1 );
  }
  GeometryPath += "/dat/BesSCM.txt";
 
  ifstream fin;
  fin.open( GeometryPath );
 
  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 ) >> innerCryo[0] >> innerCryo[1] >> innerCryo[2]; break;
    case 2: istrstream( p ) >> l1Adia[0] >> l1Adia[1] >> l1Adia[2]; break;
    case 3: istrstream( p ) >> innerShield[0] >> innerShield[1] >> innerShield[2]; break;
    case 4: istrstream( p ) >> l2Adia[0] >> l2Adia[1] >> l2Adia[2]; break;
    case 5: istrstream( p ) >> outerShield[0] >> outerShield[1] >> outerShield[2]; break;
    case 6: istrstream( p ) >> l3Adia[0] >> l3Adia[1] >> l3Adia[2]; break;
    case 7:
      istrstream( p ) >> endShield[0] >> endShield[1] >> endShield[2] >> endShield[3] >>
          endShield[4];
      break;
    case 8: istrstream( p ) >> l1insu[0] >> l1insu[1] >> l1insu[2]; break;
    case 9: istrstream( p ) >> coil[0] >> coil[1] >> coil[2]; break;
    case 10: istrstream( p ) >> l2insu[0] >> l2insu[1] >> l2insu[2]; break;
    case 11: istrstream( p ) >> supp[0] >> supp[1] >> supp[2]; break;
    case 12: istrstream( p ) >> l4Adia[0] >> l4Adia[1] >> l4Adia[2]; break;
    case 13:
      istrstream( p ) >> endCoil[0] >> endCoil[1] >> endCoil[2] >> endCoil[3] >> endCoil[4];
      break;
    case 14: istrstream( p ) >> outerCryo[0] >> outerCryo[1] >> outerCryo[2]; break;
    case 15:
      istrstream( p ) >> endCryo[0] >> endCryo[1] >> endCryo[2] >> endCryo[3] >> endCryo[4];
      break;
    case 16: istrstream( p ) >> rein[0] >> rein[1] >> rein[2] >> rein[3] >> rein[4]; break;
    case 17:
      istrstream( p ) >> pipe1[0] >> pipe1[1] >> pipe1[2] >> pipe1[3] >> pipe1[4] >> pipe1[5];
      break;
    case 18:
      istrstream( p ) >> pipe2[0] >> pipe2[1] >> pipe2[2] >> pipe2[3] >> pipe2[4] >> pipe2[5];
      break;
    case 19:
      istrstream( p ) >> pipe3[0] >> pipe3[1] >> pipe3[2] >> pipe3[3] >> pipe3[4] >> pipe3[5];
      break;
    case 20:
      istrstream( p ) >> hole[0] >> hole[1] >> hole[2] >> hole[3] >> hole[4] >> hole[5] >>
          hole[6] >> hole[7] >> hole[8] >> hole[9] >> hole[10] >> hole[11] >> hole[12];
    }
  }
}