ExtBesDetectorConstruction.cxx

Go to the documentation of this file.

//
// File: ExtBesDetectorConstruction.cxx
// Author: L.L.Wang
// Description: Constuction of BESIII detector using package GdmlToG4.
// History: 2005.8.5 created by L.L.Wang
//
 
#include "ExtBesDetectorConstruction.h"
 
#include "G4Colour.hh"
#include "G4String.hh"
#include "globals.hh"
#include <cstdlib>
#include <fstream>
 
#include "G4Evaluator/GDMLExpressionEvaluator.hh"
#include "G4Processor/GDMLProcessor.hh"
 
// Added here just to help resolve properly dependencies
#include "G4BooleanSolid.hh"
#include "G4Geo/BesG4Geo.h"
 
ExtBesDetectorConstruction::ExtBesDetectorConstruction( const bool bFieldOn, int detVer ) {
  if ( bFieldOn ) besField = new ExtBesMagneticField();
  fWorld   = 0;
  m_detVer = detVer;
}
 
ExtBesDetectorConstruction::~ExtBesDetectorConstruction() {
  if ( besField ) delete besField;
}
 
G4VPhysicalVolume* ExtBesDetectorConstruction::Construct() {
  BesG4Geo* aBesG4Geo = new BesG4Geo();
 
  fWorld = (G4VPhysicalVolume*)GDMLProcessor::GetInstance()->GetWorldVolume();
  if ( fWorld == 0 )
  {
    G4cout
        << "World volume not set properly check your setup selection criteria or GDML input!"
        << G4endl;
    exit( -1 );
  }
 
  delete aBesG4Geo;
 
  G4LogicalVolume* logicalBes = const_cast<G4LogicalVolume*>(
      GDMLProcessor::GetInstance()->GetLogicalVolume( "logicalWorld" ) );
  Mdc.Construct( logicalBes );
 
  if ( m_detVer == 1 || m_detVer == 2 )
  {
    std::cout << "ExBesDetectorConstruction:     You are using the Sintillator Tof Detector!"
              << std::endl;
    Tof.Construct( logicalBes );
  }
  else if ( m_detVer == 3 )
  {
    std::cout
        << "ExBesDetectorConstruction:     You are using the upgraded (MRPC) Tof Detector!"
        << std::endl;
    Tof.ConstructETF( logicalBes );
  }
  else
  {
    std::cout << "ExBesDetectorConstruction:     Detector Version is not well defined! "
              << m_detVer << " \n  " << std::endl;
    std::cout << "                               Use in your jobOption script " << std::endl;
    std::cout << "                               TrkExtAlg.Tof = 1; for the upgraded Tof "
                 "(MRPC Endcaps, double sided readout) "
              << std::endl;
    std::cout << "                               TrkExtAlg.Tof = 0; for the Sintillator Tof. "
              << std::endl;
  }
 
  Emc.Construct( logicalBes );
  Muc.Construct( logicalBes );
 
  return fWorld;
}