StRoot: StiSvt/StiSvtDetectorBuilder.cxx Source File

00001 #include <stdexcept>
00002 #include "StDbUtilities/StGlobalCoordinate.hh"
00003 #include "StDbUtilities/StSvtLocalCoordinate.hh"
00004 #include "StDbUtilities/StSvtCoordinateTransform.hh"
00005 #include "StiSvtDetectorBuilder.h" 
00006 #include "StSvtClassLibrary/StSvtConfig.hh"
00007 #include "StSvtClassLibrary/StSvtGeometry.hh"
00008 #include "StSvtClassLibrary/StSvtWaferGeometry.hh"
00009 #include "StSvtDbMaker/StSvtDbMaker.h"
00010 #include "StDbUtilities/St_svtRDOstrippedC.h"
00011 #include "Sti/Base/Factory.h"
00012 #include "Sti/StiPlanarShape.h"
00013 #include "Sti/StiMaterial.h"
00014 #include "Sti/StiPlacement.h"
00015 #include "Sti/StiDetector.h"
00016 #include "Sti/StiToolkit.h"
00017 #include "Sti/StiNeverActiveFunctor.h"
00018 #include "StiIsSvtActiveFunctor.h"
00019 #include "Sti/StiElossCalculator.h"
00020 #include "StDetectorDbMaker/StiSvtHitErrorCalculator.h"
00021 #include <stdio.h>
00022 #include "tables/St_HitError_Table.h"
00023 #include "StiSvtLayerLadder.h"
00024 #include "StThreeVectorD.hh"
00025 /*
00026   Geant names: SVTT the mother of all SVT volumes
00027                  SFMO: is the mother of all Silicon Strip Detector volumes
00028                  SOUM: Outer shielding structure
00029                  SXR[L,1,2]: Circular water feeds
00030                  SCBM: Mother of All Cables
00031                  SALM: aluminum shield mesh
00032                  SOSH: SVT outer shield "
00033                  SISH: SVT inner shield "
00034                  SLY[D,1,2,3,4,5]: layer mother
00035                     SLS[D,1,2,3,4,5]: ladder mother
00036                       SELE: electronics mother volume
00037                       SLDI: a ladder volume
00038                         SPCB: the G10 PCB
00039                         SRHC: roha cell wafer supports
00040                         STAB: the Berrillium tabs and the ends of the wafer carriers
00041                         SBER: the Berillium wafer carrier rails
00042                         STLI: the wafer pack container
00043                           STSI:  a single waver container
00044                             SVTD: an active wafer volume
00045                               SSIR: a non-sensitive up-down border of the wafer
00046                               SSID: a non-sensitive left-right border of the wafer
00047                               STRA: a trapezoid of triangular shape
00048                                 Total epoxy: ~10 pounds = 4.57 kG ==> 4.57 / 216 = 21.16 g/wafer *7 => 148.1 g/big ladder
00049                  SROD: Support rod
00050                  SBSP: Beampipe support mother
00051                  SCON: Support cone mother
00052                  SBWC: water manifold to support cone bracket mother
00053                  SWMM: water manifold mother
00054                  SIES: Volume to hold inner endring screws
00055                  SOES: Volume to hold outer endring screws
00056                  SBRG: Bracket joining the end rungs
00057                  SOER: outer end ring
00058                  SIRT: inner end ring tube piece 
00059                  SIRP: inner end ring polygon piece 
00060 
00061   SVT Layer and Ladder Naming Convention
00062   Hardware       StiSvtBuilder
00063   ----------------------------------
00064   Layer Ladder   Layer Ladder
00065     1     2        0     0
00066     1     4        0     1
00067     1     6        0     2
00068     1     8        0     3
00069     2     1        1     0
00070     2     3        1     1
00071     2     5        1     2
00072     2     7        1     3
00073 
00074     3     2        2     0
00075     3     4        2     1
00076     3     6        2     2
00077     3     8        2     3
00078     3     10       2     4
00079     3     12       2     5
00080     4     1        3     0
00081     4     3        3     1
00082     4     5        3     2
00083     4     7        3     3
00084     4     9        3     4
00085     4     11       3     5
00086 
00087     5     2        4     0
00088     5     4        4     1
00089     5     6        4     2
00090     5     8        4     3
00091     5     10       4     4
00092     5     12       4     5
00093     5     14       4     6
00094     5     16       4     7
00095     6     1        5     0
00096     6     3        5     1
00097     6     5        5     2
00098     6     7        5     3
00099     6     9        5     4
00100     6     11       5     5
00101     6     13       5     6
00102     6     15       5     7
00103  */
00104 static Int_t _debug = 0;
00105 StiSvtDetectorBuilder::StiSvtDetectorBuilder(bool active, const string & inputFile)
00106   : StiDetectorBuilder("Svt",active,inputFile), _siMat(0), _hybridMat(0){}
00107 
00108 StiSvtDetectorBuilder::~StiSvtDetectorBuilder() {}
00109 
00110 void StiSvtDetectorBuilder::buildDetectors(StMaker & source)
00111 {
00112   char name[50];  
00113   int nRows;
00114   cout << "StiSvtDetectorBuilder::buildDetectors() -I- Started" << endl;
00115 
00116   St_DataSet *dataSet = NULL;
00117   dataSet = source.GetDataSet("StSvtConfig");
00118   if (!dataSet) throw runtime_error("StiSvtDetectorBuilder::loadDb() -E- dataSet==0 while getting StSvtConfig");
00119   _config = static_cast<StSvtConfig*>(dataSet->GetObject());
00120   if (!_config) throw runtime_error("StiSvtDetectorBuilder::loadDb() -E- _config==0");
00121 
00122   dataSet = source.GetDataSet("StSvtGeometry");
00123         if (!dataSet)   throw runtime_error("StiSvtDetectorBuilder::loadDb() -E- dataSet==0 while getting StSvtGeometry");
00124   _geometry = static_cast<StSvtGeometry*>(dataSet->GetObject());
00125   if (!_geometry) throw runtime_error("StiSvtDetectorBuilder::loadDb() -E- _geometry==0");
00126   nRows = 2* _config->getNumberOfBarrels();
00127   setNRows(nRows);
00128   
00129   cout << "SVT Number of  Rows : "<<2* _config->getNumberOfBarrels()<<endl
00130                          << "  Layer#  numberOfLadders      Radius" << endl;
00131   cout << "StiSvtDetectorBuilder::buildDetectors() -I- Define Svt Materials" << endl;
00132   if (! _gasMat) 
00133     _gasMat    = add(new StiMaterial("Air",7.3, 14.61, 0.001205, 30420.*0.001205, 7.3*12.e-9));
00134   if (! _siMat)
00135     _siMat     = add(new StiMaterial("Si",     14.,      28.0855,   2.33,     21.82,           14.*12.*1e-9) );
00136   if (! _hybridMat)
00137     _hybridMat = add(new StiMaterial("Hybrid", 14.,      28.0855,   2.33,     21.82,           14.*12.*1e-9) );
00138   cout << "StiSvtDetectorBuilder::buildDetectors() -I- Define Svt Shapes" << endl;
00139 
00140   double ionization = _siMat->getIonization();
00141   StiElossCalculator * siElossCalculator = 
00142     new StiElossCalculator(_siMat->getZOverA(), ionization*ionization, _siMat->getA(), _siMat->getZ(), _siMat->getDensity());
00143   for (int layer=0;layer<nRows;layer++)    {
00144     cout << "  "<<layer<<"     "<<_config->getNumberOfLadders(1+layer/2)/2 << "   " 
00145          << _geometry->getBarrelRadius(layer+1) << endl;
00146     Int_t svtLayer = layer+1;
00147     Int_t svtBarrel = getSvtBarrel(svtLayer);
00148     int nWafers = _config->getNumberOfWafers(svtBarrel);
00149     // Si wafer
00150     sprintf(name, "Svt/Layer_%d/Wafers", layer);
00151     _waferShape[layer] = new StiPlanarShape(name,
00152                                             nWafers*3.15, //_geometry->getWaferLength(),
00153                                             2.*_geometry->getWaferThickness(),
00154                                             3.15);        //_geometry->getWaferWidth() );
00155     add(_waferShape[layer]);
00156     int nSectors = _config->getNumberOfLadders(svtBarrel)/2;
00157     setNSectors(layer,nSectors); 
00158     // calculate generic params for this layer
00159     // number of ladders per layer (not barrel) & phi increment between ladders
00160     //    float fDeltaPhi = M_PI/nSectors;
00161     // width of gap between the edges of 2 adjacent ladders:
00162     //   first, the angle subtended by 1/2 of the ladder
00163     float fLadderRadius  = _geometry->getBarrelRadius(svtLayer);
00164     if (fLadderRadius<=0)       throw runtime_error("StiSvtDetectorBuilder::buildDetectors() - FATAL - fLadderRadius<=0");
00165     StSvtWaferGeometry* waferGeom;
00166     Int_t index1, index2;
00167     StSvtWaferGeometry* waferGeom2;
00168     for(unsigned int ladder = 0; ladder<getNSectors(layer); ladder++)   {
00169       Int_t svtLadder = 2*(ladder+1) - (svtLayer-1)%2;
00170       Int_t wafer = nWafers/2+1;
00171       index1 = _geometry->getWaferIndex(svtBarrel,svtLadder,wafer);
00172       assert (index1 >= 0);
00173       waferGeom = (StSvtWaferGeometry*) _geometry->at(index1);
00174       if (_debug) waferGeom->print();
00175       StThreeVectorD centerVector(waferGeom->x(0), waferGeom->x(1), waferGeom->x(2) );
00176       if ( nWafers%2 == 0) {
00177         index2 = _geometry->getWaferIndex(svtBarrel,svtLadder,wafer-1);
00178         assert(index2 >= 0);
00179         waferGeom2 = (StSvtWaferGeometry*) _geometry->at(index2);
00180         StThreeVectorD centerVector2(waferGeom2->x(0), waferGeom2->x(1), waferGeom2->x(2) );
00181         if (_debug) waferGeom2->print();
00182         centerVector += centerVector2;
00183         centerVector *= 0.5;
00184       }
00185       StThreeVectorD normalVector(waferGeom->n(0), waferGeom->n(1), waferGeom->n(2) );
00186       Double_t prod = centerVector.x()*normalVector.x() + centerVector.y()*normalVector.y();
00187       if (prod < 0) normalVector *= -1;
00188       double phi  = centerVector.phi();
00189       double phiD = normalVector.phi();
00190       double r = centerVector.perp();
00191       cout <<"Det Id = "<<waferGeom->getID()<<"\tcv\t:"<<centerVector<<"\tphi:\t"<<phi<<"\tr:\t"<<r<<"\tz:\t" << centerVector.z() << endl;
00192       StiPlacement *pPlacement = new StiPlacement;
00193       pPlacement->setZcenter(centerVector.z());
00194       pPlacement->setLayerRadius(r); //this is only used for ordering in detector container...
00195       pPlacement->setLayerAngle(phi); //this is only used for ordering in detector container...
00196       pPlacement->setRegion(StiPlacement::kMidRapidity);
00197       pPlacement->setNormalRep(phiD, r*TMath::Cos(phi-phiD), r*TMath::Sin(phi-phiD)); 
00198       sprintf(name, "Svt/Layer_%d/Ladder_%d/Wafers", layer, ladder);
00199       StiDetector *pLadder = _detectorFactory->getInstance();
00200       pLadder->setName(name);
00201       pLadder->setIsOn(true);
00202       pLadder->setIsActive(new StiIsSvtActiveFunctor(svtBarrel,svtLadder, nWafers, _geometry->getWaferWidth(), _geometry->getWaferLength(), _active));
00203       pLadder->setIsContinuousMedium(true);
00204       pLadder->setIsDiscreteScatterer(true);
00205       pLadder->setGas(_gasMat);
00206       pLadder->setMaterial(_siMat);
00207       pLadder->setShape(_waferShape[layer]);
00208       pLadder->setPlacement(pPlacement); 
00209       pLadder->setHitErrorCalculator(StiSvtHitErrorCalculator::instance());
00210       pLadder->setKey(1,layer);
00211       pLadder->setKey(2,ladder);
00212       pLadder->setElossCalculator(siElossCalculator);
00213       add(layer,ladder,pLadder);
00214     } // for ladder
00215   } // for layer
00216   if (StiVMCToolKit::GetVMC()) {useVMCGeometry();}
00217   if (debug()) {
00218     cout << "StiSvtDetectorBuilder::buildDetectors list of built detectors" << endl;
00219     Int_t nlayers = _detectors.size(); 
00220     for (Int_t layer = 0; layer < nlayers; layer++) {
00221       Int_t nLadders = _detectors[layer].size();
00222       for (Int_t ladder = 0; ladder < nLadders; ladder++) {
00223         cout << "layer " << layer << "\tladder = " << ladder << "\t" << _detectors[layer][ladder]->getName() << endl;
00224       }
00225     }
00226   }
00227 }
00228 //________________________________________________________________________________
00229 void StiSvtDetectorBuilder::useVMCGeometry() {
00230   cout << "StiSvtDetectorBuilder::buildDetectors() -I- Use VMC geometry" << endl;
00231   SetCurrentDetectorBuilder(this);
00232   struct Material_t {
00233     const Char_t *name;
00234     StiMaterial    **p;
00235   };
00236   Material_t map[] = {
00237     {"AIR", &_gasMat},
00238     {"SILICON", &_siMat},
00239     {"SILICON", &_hybridMat}
00240   };
00241   Int_t M = sizeof(map)/sizeof(Material_t);
00242   for (Int_t i = 0; i < M; i++) {
00243     const TGeoMaterial *mat =  gGeoManager->GetMaterial(map[i].name); 
00244     if (! mat) continue;
00245     Double_t PotI = StiVMCToolKit::GetPotI(mat);
00246     *map[i].p = add(new StiMaterial(mat->GetName(),
00247                                     mat->GetZ(),
00248                                     mat->GetA(),
00249                                     mat->GetDensity(),
00250                                     mat->GetDensity()*mat->GetRadLen(),
00251                                     PotI));
00252   }
00253   const VolumeMap_t SvtVolumes[] = { 
00254     {"SOUM", "Outer shileding structure","HALL_1/CAVE_1/SVTT_1/SOUM_1/*","",""},
00255     {"SXRL", "Circular water feeds","HALL_1/CAVE_1/SVTT_1/SXRL_1-2/*","",""}, 
00256     {"SXR1", "Circular water feeds","HALL_1/CAVE_1/SVTT_1/SXR1_3-4/*","",""}, 
00257     {"SXR2", "Circular water feeds","HALL_1/CAVE_1/SVTT_1/SXR2_5-6/*","",""},
00258     //    {"SCBM", "Mother of All Cables","HALL_1/CAVE_1/SVTT_1/SCBM_1-2/*","",""},
00259     {"SCBL", "The bundles of cables connecting PCBs with the transition boards","HALL_1/CAVE_1/SVTT_1/SCBM_1-2/SCBL_1","",""},
00260     {"SCB1", "The bundles of cables connecting PCBs with the transition boards","HALL_1/CAVE_1/SVTT_1/SCBM_1-2/SCB1_2","",""},
00261     {"SCB2", "The bundles of cables connecting PCBs with the transition boards","HALL_1/CAVE_1/SVTT_1/SCBM_1-2/SCB2_3","",""},
00262     {"SCB3", "The bundles of cables connecting PCBs with the transition boards","HALL_1/CAVE_1/SVTT_1/SCBM_1-2/SCB3_4","",""},
00263     {"SFED", "bundles of water pipes","HALL_1/CAVE_1/SVTT_1/SCBM_1-2/SFED_1","",""},
00264     {"SFE1", "bundles of water pipes","HALL_1/CAVE_1/SVTT_1/SCBM_1-2/SFE1_2","",""},
00265     {"SFE2", "bundles of water pipes","HALL_1/CAVE_1/SVTT_1/SCBM_1-2/SFE2_3","",""},
00266     {"SPLS", "plastic of the water pipes","HALL_1/CAVE_1/SVTT_1/SCBM_1-2/SPLS_1","",""},
00267     {"SPL1", "plastic of the water pipes","HALL_1/CAVE_1/SVTT_1/SCBM_1-2/SPL1_2","",""},
00268     {"SPL2", "plastic of the water pipes","HALL_1/CAVE_1/SVTT_1/SCBM_1-2/SPL2_3","",""},
00269     {"SALM", "aluminum shield mesh","HALL_1/CAVE_1/SVTT_1/SALM_1-2","",""},
00270     {"SOSH", "SVT outer shield","HALL_1/CAVE_1/SVTT_1/SOSH_1","",""},
00271     {"SISH", "SVT inner shield","HALL_1/CAVE_1/SVTT_1/SISH_1","",""},
00272 //     {"SLYD", "layer mother","HALL_1/CAVE_1/SVTT_1/SLYD_1/*","",""},
00273 //     {"SLY1", "layer mother","HALL_1/CAVE_1/SVTT_1/SLY1_2/*","",""},
00274 //     {"SLY2", "layer mother","HALL_1/CAVE_1/SVTT_1/SLY2_3/*","",""},
00275 //     {"SLY3", "layer mother","HALL_1/CAVE_1/SVTT_1/SLY3_4/*","",""},
00276 //     {"SLY4", "layer mother","HALL_1/CAVE_1/SVTT_1/SLY4_5/*","",""},
00277 //     {"SLY5", "layer mother","HALL_1/CAVE_1/SVTT_1/SLY5_6/*","",""},
00278     {"SELE","electronics mother volume","HALL_1/CAVE_1/SVTT_1/SLYD_1/SLSD_1/SELE_1","",""},//   Weight = 0.0422223[kG]
00279 #if 1
00280     {"SEL1","electronics mother volume","HALL_1/CAVE_1/SVTT_1/SLYD_2/SLSD_1/SELE_1","",""},//   Weight = 0.0422223[kG]
00281     {"SEL2","electronics mother volume","HALL_1/CAVE_1/SVTT_1/SLYD_3/SLSD_1/SELE_1","",""},//   Weight = 0.048099[kG]
00282     {"SEL3","electronics mother volume","HALL_1/CAVE_1/SVTT_1/SLYD_4/SLSD_1/SELE_1","",""},//   Weight = 0.048099[kG]
00283     {"SEL4","electronics mother volume","HALL_1/CAVE_1/SVTT_1/SLYD_5/SLSD_1/SELE_1","",""},//   Weight = 0.0510373[kG]
00284     {"SEL5","electronics mother volume","HALL_1/CAVE_1/SVTT_1/SLYD_6/SLSD_1/SELE_1","",""},//   Weight = 0.0510373[kG]
00285 #endif
00286     //  {"SVTD", "an active wafer volume","HALL_1/CAVE_1/SVTT_1/SLY*/SLS*/SLD*/STL*/STS*/SVTD_1","svt","SVTD"}, // <+++
00287     {"SBWC", "water manifold to support cone bracket mother","HALL_1/CAVE_1/SVTT_1/SBWC_1-2/*","",""},
00288     {"SWMM", "water manifold mother","HALL_1/CAVE_1/SVTT_1/SWMM_1-2/*","",""},
00289     {"SIES", "Volume to hold inner endring screws","HALL_1/CAVE_1/SVTT_1/SIES_1-2/*","",""},
00290     {"SOES", "Volume to hold outer endring screws","HALL_1/CAVE_1/SVTT_1/SOES_1-2/*","",""},
00291     {"SBRG", "Bracket joining the end rungs","HALL_1/CAVE_1/SVTT_1/SBRG_1-2/*","",""},
00292     {"SOER", "outer end ring","HALL_1/CAVE_1/SVTT_1/SOER_1-2/*","",""},
00293     {"SIRT", "inner end ring tube piece ","HALL_1/CAVE_1/SVTT_1/SIRT_1-2","",""},
00294     {"SIRP", "inner end ring polygon piece ","HALL_1/CAVE_1/SVTT_1/SIRP_1-2","",""}
00295     //?    {"STAC", "twinax cable approximation, copper","HALL_1/CAVE_1/SVTT_1/SCON_1/STAC_1-2","",""} 
00296     // StiDetectorVolume        SVTT_1_SCON_1_STAC_1    StiDetector OBJ: TTUBE  STAC    StiCylindricalShape Rmin 21.935 Rmax 37.77 dz 86.9 
00297   };
00298   Int_t NoSvtVols = sizeof(SvtVolumes)/sizeof(VolumeMap_t);
00299   TString pathT("HALL_1/CAVE_1");
00300   TString path("");
00301   for (Int_t i = 0; i < NoSvtVols; i++) {
00302     gGeoManager->RestoreMasterVolume(); 
00303     gGeoManager->CdTop();
00304     gGeoManager->cd(pathT); path = pathT;
00305     TGeoNode *nodeT = gGeoManager->GetCurrentNode();
00306     if (! nodeT) continue;
00307     StiVMCToolKit::LoopOverNodes(nodeT, path, SvtVolumes[i].name, MakeAverageVolume);
00308   }
00309 }