KVFAZIA3DBlock.cpp

//Created by KVClassFactory on Thu May 30 10:13:17 2013
//Author: Eric Bonnet
 
#include "KVFAZIA3DBlock.h"
#include "KVUnits.h"
#include "KVMaterial.h"
#include "TMath.h"
//#include "KVFAZIA.h"
#include "TGeoManager.h"
#include "TGeoMedium.h"
#include "TGeoVolume.h"
#include "TGeoMatrix.h"
#include "KVDataSet.h"
 
ClassImp(KVFAZIA3DBlock)
 
// BEGIN_HTML <!--
/* -->
<h2>KVFAZIA3DBlock</h2>
<h4>Block of 16 telescopes FAZIA</h4>
<!-- */
// --> END_HTML
 
 
 
KVFAZIA3DBlock::KVFAZIA3DBlock(int block_id) : KVFAZIABlock()
{
   idb = block_id;
   // Default constructor
}
 
 
 
 
void KVFAZIA3DBlock::Build()
{
   quartet_position.ReadFile(Form("%s/%s/FAZIA3D.Quartet.Position.env", gDataSet->GetDATADIRFilePath(), gDataSet->GetName()), kEnvLocal);
 
   // Set default dimensions for all parts of geometry
   frame_side = 4.415;
   centre_cross_thick = 0.12;
   side_thick = 0.0475;
   si_holder_side = 2.1;
   si_centre = 0.5 * centre_cross_thick + 0.5 * si_holder_side;
   front_frame_depth = 0.3;
   frame_depth_300 = 0.273;
   frame_depth_500 = 0.24;
   frame_edge_depth = 0.16;
   si300_position = -0.5 * frame_depth_300 + frame_edge_depth;
   si500_position = -0.5 * frame_depth_500 + frame_edge_depth;
   full_width = front_frame_depth + frame_depth_300 + frame_depth_500;
   dist_si2_csi = 0.;// should be ~1mm between si2 ergal edge and csi
   csi_depth = 10;
   full_length_quartet = full_width + dist_si2_csi + csi_depth;
   thick_si1 = 300 * KVUnits::um;
   thick_si2 = 500 * KVUnits::um;
   dist_csi_sarco = 4.25;
   sarco_length = 49.7;
   sarco_width = 9.6;
   sarco_height = 10.4;
   sarco_epaisseur = 0.4;
   sarco_uplift = 0.5; //how much the sarco sides overlap the copper plate
 
   SetMedium(gGeoManager->GetMedium("Vacuum"));//to avoid warnings about STRUCT_BLOCK has dummy medium
 
   TGeoVolume* quartet = MakeQuartet();
 
   for (Int_t nq = 1; nq <= 4; nq += 1) {
      KVString qn = Form("B%03d-Q%d", idb, nq);
 
      // initial translation to set the Si1 entrance window at Z=0
      TGeoTranslation trans0;
      trans0.SetDz(4.9615 + front_frame_depth);
 
      // tree rotations around random axes to set the proper quartet orientation
      double psi = quartet_position.GetValue(Form("%s.Psi", qn.Data()), 0.);
      double theta = quartet_position.GetValue(Form("%s.Theta", qn.Data()), 0.);
      double phi = quartet_position.GetValue(Form("%s.Phi", qn.Data()), 0.);
 
      TGeoRotation rot1;
      rot1.RotateZ(psi);
 
      TGeoRotation rot2;
      rot2.RotateY(theta);
 
      TGeoRotation rot3;
      rot3.RotateZ(phi);
 
      // translation to set the quartet center in place
      TGeoTranslation trans;
      trans.SetDx(quartet_position.GetValue(Form("%s.X", qn.Data()), 0.));
      trans.SetDy(quartet_position.GetValue(Form("%s.Y", qn.Data()), 0.));
      trans.SetDz(quartet_position.GetValue(Form("%s.Z", qn.Data()), 0.));
 
      TGeoHMatrix h = trans * rot3 * rot2 * rot1 * trans0;
      TGeoHMatrix* ph = new TGeoHMatrix(h);
 
      AddNode(quartet, nq, ph);
   }
 
}