PisaHALL.C

Go to the documentation of this file.

#include "PisaHALL.h"

#include <cassert>
#include <string>
#include <sstream>

#include "TGeoManager.h"
#include "TGeoMaterial.h"
#include "TGeoMatrix.h"
#include "TGeoMedium.h"
#include "TGeoPcon.h"
#include "TVirtualMC.h"

#include "PHCompositeNode.h"
#include "PisaPhnxPar.h"
#include "PisaHelper.h"
#include "PisaServer.h"

//_____________________________________________________________________________
PisaHALL::~PisaHALL()
{
}

//_____________________________________________________________________________
void
PisaHALL::createMaterials(PHCompositeNode*)
{
  TGeoManager* gm = PisaHelper::getGeoManager();

  Double_t a2[2] = { 14.01, 16.00};
  Double_t z2[2] = {  7.0,   8.0};
  Double_t w2[2] = {  0.7,   0.3};
  Double_t density = 1.29e-03;  
  int imatair = 15;

  gm->Mixture("Air", a2, z2, density, 2, w2, imatair);

  int isvol = 0; // not a sensitive volume
  Int_t ifield = 2;          // User defined magnetic field
  Double_t fieldm = 5.;     // Maximum field value (in kiloGauss)
  Double_t epsil  = .001;    // Tracking precision, 
  Double_t stemax = -0.01;   // Maximum displacement for multiple scat 
  Double_t tmaxfd = -20.;    // Maximum angle due to field deflection 
  Double_t deemax = -.3;     // Maximum fractional energy loss, DLS 
  Double_t stmin  = -.8;

  fImedAIRF = 1;
  fImedAIRLOWF = 2;

  // AIR, High Magnetic Field
  gm->Medium("AIRF",fImedAIRF,imatair,isvol,
             ifield, fieldm, tmaxfd, stemax, 
             deemax, epsil, stmin);
                                     
 
  // AIR, Low Magnetic Field
  fieldm = 1.0;
  gm->Medium("AIRLOWF",fImedAIRLOWF,imatair,isvol,
             ifield, fieldm, tmaxfd, stemax, 
             deemax, epsil, stmin);

}

//_____________________________________________________________________________
void
PisaHALL::createGeometry(PHCompositeNode* topNode)
{
  assert(gMC->IsRootGeometrySupported()==true);
                                
  PisaPhnxPar* parameters = PisaHelper::getClass<PisaPhnxPar>(topNode,"PHNXPAR");
  assert(parameters!=0);

  TGeoManager* gm = PisaHelper::getGeoManager();

  TGeoMedium* mairf = gm->GetMedium("AIRF");
  assert(mairf!=0);
  TGeoMedium* mairlowf = gm->GetMedium("AIRLOWF");
  assert(mairlowf!=0);

  const PisaPhnxSubPar* hall_par = parameters->get("hall_par");
  assert(hall_par!=0);
 
  double dx = hall_par->getDouble("dim",0);
  double dy = hall_par->getDouble("dim",1);
  double dz = hall_par->getDouble("dim",2);

  TGeoVolume* hall = gm->MakeBox("HALL",mairf,dx,dy,dz);
  
  PisaServer* ps = PisaServer::instance();

  if ( ps->hasSubsystem("MVD") )
    {
      const PisaPhnxSubPar* vert_par = parameters->get("vert_par");
      assert(vert_par!=0);
      TGeoVolume* vert = gm->MakeTube("VERT",mairf,
                                      vert_par->getDouble("dim",0),
                                      vert_par->getDouble("dim",1),
                                      vert_par->getDouble("dim",2));

      hall->AddNode(vert,0,new TGeoTranslation
                    (vert_par->getDouble("pos",0),
                     vert_par->getDouble("pos",1),
                     vert_par->getDouble("pos",2)
                     ));
    }

  if ( ps->hasSubsystem("DC") && hall_par->getInt("arm_shift")==0 )
    {
      const PisaPhnxSubPar* intr_par = parameters->get("intr_par");
      assert(intr_par!=0);
      TGeoVolume* intr = gm->MakeTube("INTR",mairf,
                                      intr_par->getDouble("dim",0),
                                      intr_par->getDouble("dim",1),
                                      intr_par->getDouble("dim",2));

      hall->AddNode(intr,0);
    }
  
  // The tracking medium, from PC2 and beyond, is AIR, Low Field

  if ( hall_par->getInt("arm_shift")==0 )
    {
      // Outer pseudo-volume (region beyond the RICH)
      // This is the mother volume for the PC2, TRD, PC3, TOF, and EMCAL
      
      const PisaPhnxSubPar* emcl_par = parameters->get("emcl_par");
      assert(emcl_par!=0);
      int nz = emcl_par->getInt("nz"); // number of PCON sections.
      TGeoPcon* emcl_shape = new TGeoPcon(emcl_par->getDouble("phi"),
                                          emcl_par->getDouble("dphi"),
                                          nz);

      for ( int isection = 0; isection < nz; ++isection )
        {
          std::ostringstream tmp;

          tmp << "section" << isection;

          std::string section = tmp.str();

          emcl_shape->DefineSection(isection,
                                     emcl_par->getDouble(section.c_str(),0),
                                     emcl_par->getDouble(section.c_str(),1),
                                     emcl_par->getDouble(section.c_str(),2));
        }

      TGeoVolume* emcl = new TGeoVolume("EMCL",emcl_shape,mairlowf);

      hall->AddNode(emcl,0);
    }
}