KVFAZIALNS2016.cpp

//Created by KVClassFactory on Tue Jan 27 11:38:09 2015
//Author: ,,,
 
#include "KVFAZIALNS2016.h"
#include "KVUnits.h"
#include "KVFAZIABlock.h"
#include "TSystem.h"
#include "KVEnv.h"
 
#include <TGeoMatrix.h>
 
ClassImp(KVFAZIALNS2016)
 
// BEGIN_HTML <!--
/* -->
<h2>KVFAZIALNS2016</h2>
<h4>Description of the FAZIA set up</h4>
<!-- */
// --> END_HTML
 
 
 
KVFAZIALNS2016::KVFAZIALNS2016()
   : KVFAZIA()
{
   // Default constructor
   SetTitle(ClassName());
}
 
 
 
 
KVFAZIALNS2016::~KVFAZIALNS2016()
{
   // Destructor
}
 
 
 
 
void KVFAZIALNS2016::GetGeometryParameters()
{
   //Defined number of blocks, the distance from the target and the minimum polar angle
   fNblocks = 5;
   fFDist = 80.0;
   fFThetaMin = 2.1;
}
 
 
 
 
void KVFAZIALNS2016::RutherfordTelescope()
{
   // Telescope for elastic scattering monitoring
   // Two 5mm diameter silicon detectors of 525um thickness
   // placed 2m20 from the target at theta=1.84deg phi=-90deg.
   // distance between centres of detectors = 1mm
 
   KVMaterial silicon("Si");
 
   const double radius = 5 * KVUnits::mm / 2.;
 
   const double thick = 525 * KVUnits::um;
   const double centre_dist = 1 * KVUnits::mm;
   double total_thickness = thick + centre_dist;
 
   TGeoVolume* si_det = gGeoManager->MakeTube("DET_SI", silicon.GetGeoMedium(), 0., radius, thick / 2);
 
   TGeoVolumeAssembly* ruth_tel = gGeoManager->MakeVolumeAssembly("STRUCT_RUTH");
 
   ruth_tel->AddNode(si_det, 1, new TGeoTranslation(0, 0, -centre_dist / 2));
   ruth_tel->AddNode(si_det, 2, new TGeoTranslation(0, 0, centre_dist / 2));
 
   // front entrance of first detector at 2 metres from target
   const double distance = 2.20 * KVUnits::m + 0.5 * total_thickness;
   const double theta = 1.84;
   const double phi = -90;
 
   TGeoRotation rot1, rot2;
   rot2.SetAngles(phi + 90, theta, 0);
   rot1.SetAngles(-90, 0., 0.);
   TGeoTranslation trans(0, 0, distance);
   TGeoHMatrix h = rot2 * trans * rot1;
   gGeoManager->GetTopVolume()->AddNode(ruth_tel, 1, new TGeoHMatrix(h));
}
 
 
 
 
 
void KVFAZIALNS2016::BuildFAZIA()
{
   //Build geometry of FAZIASYM
   //All telescopes are : Si(300µm)-Si(500µm)-CsI(10cm)
   //No attempt has been made to implement real thicknesses
   //
   Info("BuildFAZIA", "Compact geometry, %f cm from target",
        fFDist);
 
   TGeoVolume* top = gGeoManager->GetTopVolume();
 
   Double_t distance_block_cible = fFDist * KVUnits::cm;
   Double_t thick_si1 = 300 * KVUnits::um;
   TGeoTranslation trans;
   trans.SetDz(distance_block_cible + thick_si1 / 2.);
 
   KVFAZIABlock* block = new KVFAZIABlock;
   block->Build();
 
   TGeoRotation rot1, rot2;
   TGeoHMatrix h;
   TGeoHMatrix* ph = 0;
   Double_t theta = 0;
   Double_t phi = 0;
 
   Double_t theta_min = fFThetaMin;//smallest lab polar angle in degrees
   Double_t centre_hole = 2.*tan(theta_min * TMath::DegToRad()) * distance_block_cible;
   Double_t dx = (block->GetTotalSideWithBlindage()) / 2.;
 
   TVector3 centre;
   for (Int_t bb = 0; bb < fNblocks; bb += 1) {
 
      if (bb == 1)        centre.SetXYZ(-1 * (dx - centre_hole / 2), -dx - centre_hole / 2, distance_block_cible);
      else if (bb == 2)   centre.SetXYZ(-1 * (dx + centre_hole / 2), dx - centre_hole / 2, distance_block_cible);
      else if (bb == 3)   centre.SetXYZ(-1 * (-dx + centre_hole / 2), dx + centre_hole / 2, distance_block_cible);
      else if (bb == 0)   centre.SetXYZ(-1 * (-dx - centre_hole / 2), -dx + centre_hole / 2, distance_block_cible);
      else if (bb == 4)   centre.SetXYZ(-1 * (-dx - centre_hole / 2), -3 * dx + centre_hole / 2, distance_block_cible); //centre.SetXYZ(-1 * (dx - centre_hole / 2), -3 * dx - centre_hole / 2, distance_block_cible);
      else {
         KVError::Warning(this, "BuildFAZIA", "Block position definition is done only for %d blocks", fNblocks);
      }
      theta = centre.Theta() * TMath::RadToDeg();
      phi = centre.Phi() * TMath::RadToDeg();
      printf("BLK #%d => theta=%1.2lf - phi=%1.2lf\n", bb, theta, phi);
 
      rot2.SetAngles(phi + 90., theta, 0.);
      rot1.SetAngles(-1.*phi, 0., 0.);
      h = rot2 * trans * rot1;
      ph = new TGeoHMatrix(h);
      top->AddNode(block, bb, ph);
   }
 
   // add telescope for elastic scattering monitoring
//   RutherfordTelescope();
   // Change default geometry import angular range for rutherford telescope
   SetGeometryImportParameters(.25, 1., 1.84);
}
 
 
 
 
void KVFAZIALNS2016::SetNameOfDetectors(KVEnv& env)
{
 
   KVFAZIA::SetNameOfDetectors(env);
   env.SetValue("RUTH_SI_1", "SI1-RUTH");
   env.SetValue("RUTH_SI_2", "SI2-RUTH");
 
}