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 "TrkExtAlg/ExtBesDetectorConstruction.h"
 
#include "G4Colour.hh"
#include "globals.hh"
#include <cstdlib>
#include <fstream>
#include "G4String.hh"
 
#include "GDMLProcessor.hh"
#include "GDMLExpressionEvaluator.hh"
 
// Added here just to help resolve properly dependencies
#include "G4BooleanSolid.hh"
#include "G4Geo/BesG4Geo.h"
 
 
ExtBesDetectorConstruction::ExtBesDetectorConstruction(const bool bFieldOn, int tofversion)
{
    if(bFieldOn) besField = new ExtBesMagneticField();
    fWorld=0;
    m_tofversion=tofversion;
}
 
 
ExtBesDetectorConstruction::~ExtBesDetectorConstruction()
{
    if(besField)
    delete besField;
}
 
 
G4VPhysicalVolume* ExtBesDetectorConstruction::Construct()
{
    BesG4Geo* aBesG4Geo = new BesG4Geo();
 
    fWorld=(G4VPhysicalVolume *)GDMLProcessor::GetInstance()->GetWorldVolume();  
    if( fWorld == 0 ) {
              G4Exception("World volume not set properly check your setup selection criteria or GDML input!");
        }
 
    delete aBesG4Geo;
                
    G4LogicalVolume* logicalBes = const_cast<G4LogicalVolume*>( GDMLProcessor::GetInstance()->GetLogicalVolume("logicalWorld" ) );
    Mdc.Construct(logicalBes);
 
    if(m_tofversion==2)
      { 
        std::cout << "ExBesDetectorConstruction:     You are using the old Tofdetector " << std::endl;
        Tof.Construct(logicalBes);
      }
    else if(m_tofversion==4)
      {
        std::cout << "ExBesDetectorConstruction:     You are using the new (MRPC) Tofdetector " << std::endl;
            Tof.Construct_mrpc(logicalBes);
          }
    else
      {
        std::cout << "ExBesDetectorConstruction:     TofDetector is not well defined! \n  " << std::endl;
        std::cout << "                               Use in your jobOption script " << std::endl;
        std::cout << "                               TrkExtAlg.Tof = 4; for the new Tof (MRPC Endcaps) " << std::endl;
        std::cout << "                               TrkExtAlg.Tof = 2; for the old Tof. " << std::endl;
      }
 
    
        Emc.Construct(logicalBes);
        Muc.Construct(logicalBes);
    
    return fWorld;
}