PisaMAG.C

Go to the documentation of this file.

#include "PisaMAG.h"

#include <cassert>
#include <cstdlib>

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

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

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

//_____________________________________________________________________________
void
PisaMAG::buildCM(PHCompositeNode* topNode)
{
  PisaPhnxPar* parameters = PisaHelper::getClass<PisaPhnxPar>(topNode, "PHNXPAR");
  assert(parameters != 0);

  const PisaPhnxSubPar* mag_par = parameters->get
    ("mag_par");
  assert(mag_par != 0);

  TGeoManager* gm = PisaHelper::getGeoManager();

  const char* MAG[] = { "MAGD","MAGU" };
  const char* AMG[] = { "AMGD","AMGU" };
  const int sign[] = { +1, -1 };

  for ( size_t i = 0; i < 2; ++i )
    {
      // Loop on north/south.

      int nz = 14; // number of z planes for the following PCON

      TGeoPcon* mag_shape = new TGeoPcon(0.0,      // phi min
                                         360.0,    // dphi
                                         nz);
      
      for ( int iz = 0; iz < nz; ++iz )
        {
          int index = (nz-1)*i + sign[i]*iz;

          mag_shape->DefineSection
            (iz,
             sign[i]*mag_par->getDouble("z_cent",index),
             mag_par->getDouble("rmin_cent",index),
             mag_par->getDouble("rmax_cent",index)
             );
        }
      
      TGeoMedium* medium = gm->GetMedium(mag_par->getInt("nmed_cent"));
      assert(medium != 0);
      
      TGeoVolume* magd = new TGeoVolume(MAG[i], mag_shape, medium);
      
      int nzextra = 3;
      
      mag_shape = new TGeoPcon(0.0, 360.0, nzextra);
      
      for ( int iz = 0; iz < nzextra; ++iz )
        {
          int index = (nzextra-1)*i + sign[i]*iz + nz-1;
          mag_shape->DefineSection(iz,
                                   sign[i]*mag_par->getDouble("z_cent",index),
                                   mag_par->getDouble("rmin_cent",index),
                                   mag_par->getDouble("rmax_cent",index));
        }
      
      TGeoVolume* amgd = new TGeoVolume(AMG[i], mag_shape, medium);
      
      // The inner and outer coils
      
      medium = gm->GetMedium(mag_par->getInt("cencoil_med", 0));
      assert(medium != 0);
      
      double z1 = mag_par->getDouble("cencoil_z", 0);
      double z2 = mag_par->getDouble("cencoil_z", 1);
      
      TGeoVolume* m1co = gm->MakeTube("M1CO", medium,
                                      mag_par->getDouble("cencoil_rin", 0),
                                      mag_par->getDouble("cencoil_rout", 0),
                                      0.5 * (z2 - z1));
      
      magd->AddNode(m1co, 1, new TGeoTranslation(0.0, 0.0, 
                                                 sign[i]*0.5*(z1 + z2)));
      
      
      medium = gm->GetMedium(mag_par->getInt("cencoil_med", 1));
      assert(medium != 0);
      
      z1 = mag_par->getDouble("cencoil_z", 2);
      z2 = mag_par->getDouble("cencoil_z", 3);
      double zsize = 0.5 * (z2 - z1) / 6.0;
      
      TGeoVolume* m2co = gm->MakeTube("M2CO", medium,
                                      mag_par->getDouble("cencoil_rin", 1),
                                      mag_par->getDouble("cencoil_rout", 1),
                                      zsize);
      
      for ( int j = 1; j <= 6; ++j )
        {
          double zcoil = sign[i]*(0.5 * (z2 + z1) + (j - 3) * zsize * 2.0) 
            + 0.1;
          magd->AddNode(m2co,j, new TGeoTranslation(0.0, 0.0, zcoil));
        }
      
      // Attach the volumes to the HALL
      
      TGeoVolume* hall = gm->GetVolume("HALL");
      assert(hall != 0);
      
      hall->AddNode(amgd, 1);
      hall->AddNode(magd, 1);
      
    }

  // FIXME: use the colors of the phnx.par file, if available

 //  gm->GetVolume("MAGD")->SetLineColor(1);
//   gm->GetVolume("AMGD")->SetLineColor(2);
//   gm->GetVolume("M1CO")->SetLineColor(3);
//   gm->GetVolume("M2CO")->SetLineColor(4);
}

//_____________________________________________________________________________
void
PisaMAG::buildYokes(PHCompositeNode* topNode)
{
  PisaPhnxPar* parameters = PisaHelper::getClass<PisaPhnxPar>(topNode, "PHNXPAR");
  assert(parameters != 0);

  const PisaPhnxSubPar* mag_par = parameters->get
    ("mag_par");
  assert(mag_par != 0);

  TGeoManager* gm = PisaHelper::getGeoManager();

  int cyoke_flag = mag_par->getInt("cyoke_flag");

  //         cyoke_flag = 0   means no yoke volumes installed
  //         cyoke_flag = 1   means only upper (+Y) return yoke is installed
  //         cyoke_flag = 2   means both upper and lower return yokes installed
  //         cyoke_flag = -1  means only lower (-Y) return yoke is installed
  

  // CYOKE is the "horizontal" cross piece from Z = -242 to +242 cm,
  //          159.85 full width along X, and 100 cm full height along Y
  //
  // DYOKE is the "vertical" piece on the North and South ends; it goes from
  //          Z = 135 to 242 cm (North Side), and from 217 to 295 cm vertical;
  //          The X full width is the same as the CYOKE.
  //          This model is slightly inaccurate since there is a base X width
  //          of 180.0 full width extending from Y = 217.0 to 234.0 cm, and
  //          the DYOKE volume overlaps with the central steel shape.
  //

  TGeoMedium* medium = gm->GetMedium(mag_par->getInt("nmed_cent"));
  assert(medium != 0);

  if ( cyoke_flag == 1 || cyoke_flag == 2 )
    {
      // Upper Parts of Return Yoke
      TGeoVolume* crtu = gm->MakeBox("CRTU", medium,
                                     mag_par->getDouble("cyoke_dx"),
                                     mag_par->getDouble("cyoke_dy"),
                                     mag_par->getDouble("cyoke_dz"));

      TGeoVolume* crru = gm->MakeBox("CRRU", medium,
                                     mag_par->getDouble("dyoke_dx"),
                                     mag_par->getDouble("dyoke_dy"),
                                     mag_par->getDouble("dyoke_dz"));

      TGeoVolume* crlu = gm->MakeBox("CRLU", medium,
                                     mag_par->getDouble("dyoke_dx"),
                                     mag_par->getDouble("dyoke_dy"),
                                     mag_par->getDouble("dyoke_dz"));

      TGeoVolume* hall = gm->GetVolume("HALL");
      assert(hall != 0);

      hall->AddNode(crtu, 0, new TGeoTranslation
                    (0.0, mag_par->getDouble("cyoke_yp"), 0.0));

      hall->AddNode(crru, 0, new TGeoTranslation
                    (0.0, mag_par->getDouble("dyoke_yp"),
                     mag_par->getDouble("dyoke_zp")));

      hall->AddNode(crlu, 0, new TGeoTranslation
                    (0.0, mag_par->getDouble("dyoke_yp"),
                     -mag_par->getDouble("dyoke_zp")));
    }

  if ( cyoke_flag == -1 || cyoke_flag == 2 )
    {
      // Lower Parts of Return Yoke
      TGeoVolume* crtd = gm->MakeBox("CRTD", medium,
                                     mag_par->getDouble("cyoke_dx"),
                                     mag_par->getDouble("cyoke_dy"),
                                     mag_par->getDouble("cyoke_dz"));

      TGeoVolume* crrd = gm->MakeBox("CRRD", medium,
                                     mag_par->getDouble("dyoke_dx"),
                                     mag_par->getDouble("dyoke_dy"),
                                     mag_par->getDouble("dyoke_dz"));

      TGeoVolume* crld = gm->MakeBox("CRLD", medium,
                                     mag_par->getDouble("dyoke_dx"),
                                     mag_par->getDouble("dyoke_dy"),
                                     mag_par->getDouble("dyoke_dz"));

      TGeoVolume* hall = gm->GetVolume("HALL");
      assert(hall != 0);

      hall->AddNode(crtd, 0, new TGeoTranslation
                    (0.0, -mag_par->getDouble("cyoke_yp"), 0.0));

      hall->AddNode(crrd, 0, new TGeoTranslation
                    (0.0, -mag_par->getDouble("dyoke_yp"),
                     mag_par->getDouble("dyoke_zp")));

      hall->AddNode(crld, 0, new TGeoTranslation
                    (0.0, -mag_par->getDouble("dyoke_yp"),
                     -mag_par->getDouble("dyoke_zp")));
    }
}

//_____________________________________________________________________________
void
PisaMAG::createGeometry(PHCompositeNode* topNode)
{
  // A few checks to start with.

  assert(gMC->IsRootGeometrySupported() == true); // insure we deal with VMC

  TGeoManager* gm = PisaHelper::getGeoManager();

  TGeoVolume* hall = gm->GetVolume("HALL");

  if (!hall)
    {
      std::cerr << PHWHERE << "HALL Volume does not exist ! Aborting now"
                << std::endl;
      exit(1);
    }

  buildCM(topNode);

  buildYokes(topNode);
}

//_____________________________________________________________________________
void
PisaMAG::createMaterials(PHCompositeNode*)
{
  TGeoMaterial* iron = new TGeoMaterial("IRON", 55.85, 26, 7.87);
  iron->SetUniqueID(10);

  TGeoMaterial* copper = new TGeoMaterial("COPPER", 63.54, 29, 8.96);
  copper->SetUniqueID(11);

  new TGeoMedium("YOKE", 5, iron->GetUniqueID(),
                 0, 1, 5, 0.3, 0.1000000E+11, 0.25,
                 0.1, 0.2813469E-01);

  new TGeoMedium("CU", 15, copper->GetUniqueID(), 0, 1, 5, 1,
                 0.1000000E+11, 0.25, 0.1, 0.2813862E-01);
}