FairRoot/PandaRoot
PndMdtIGeometry.cxx
Go to the documentation of this file.
1 // -------------------------------------------------------------------------
2 // ----- PndMdtIGeometry source file -----
3 // author, Jifeng Hu, hu@to.infn.it, Torino University
4 // -------------------------------------------------------------------------
5 
6 #include "PndMdtIGeometry.h"
7 #include "TVector3.h"
8 #include "TRandom.h"
9 #include "TGeoManager.h"
10 #include "TGeoBBox.h"
11 #include "TMath.h"
12 #include <assert.h>
13 #include <iostream>
14 #include <math.h>
15 #include <list>
16 #include <algorithm>
17 #include "PndMdtID.h"
18 #include <fstream>
19 
20 using std::cerr;
21 using std::cout;
22 using std::endl;
23 ostream& operator<<(ostream& os, const TVector3& v3)
24 {
25  os <<'('<<v3.X()<<", "<<v3.Y()<<", "<<v3.Z()<<')';
26  return os;
27 }
28 
29 PndMdtIGeometry* PndMdtIGeometry::fInstance = 0;
30 PndMdtIGeometry* PndMdtIGeometry::Instance()
31 {
32  if(!fInstance) fInstance = new PndMdtIGeometry();
33  return fInstance;
34 }
35 // ----- Default constructor -------------------------------------------
36 PndMdtIGeometry::PndMdtIGeometry()
37  : TNamed("PndMdtIGeometry", "parameterized simulation of MDT")
38 {
39  fNumofMdt = 0;
40  fNumofStrip = 0;
41  fVerbose = 0;
42  fGoodGeometry = kFALSE;
43  fStateTip = "Success";
44  fBarrelVersion = fEndcapVersion = fMuonFilterVersion = fForwardVersion = -1;
45  fInitialized = kFALSE;
46 }
47 PndMdtIGeometry::~PndMdtIGeometry()
48 {
49  //fWirePosMap.clear();
50  //fWireLengthMap.clear();
51  //fStripPosMap.clear();
52  //fStripLengthMap.clear();
53  fGeoMap.clear();
54  fLayerInfoMap.clear();
55 }
56 // -------------------------------------------------------------------------
57 Bool_t PndMdtIGeometry::Init()
58 {
59  if(!fInitialized){
60  fInitialized = kTRUE;
61  TGeoVolume* fTopVol = gGeoManager->GetTopVolume();
62  TGeoNode* fTopNode = gGeoManager->GetTopNode();
63  fGoodGeometry = kFALSE;
64  LoadGeometry(fTopNode);
65  if(fGoodGeometry){
66  fGoodGeometry = GetStripInfo();
67  }
68  if(fVerbose>0){
69  cout<<"========================================================================="<<endl;
70  cout<<"PndMdtIGeometry::Init(), information of MDT geometry version, \n"
71  <<"\tBarrel #"<<fBarrelVersion<<endl
72  <<"\tEndcap #"<<fEndcapVersion<<endl
73  <<"\tForward #"<<fForwardVersion<<endl
74  <<"\tFilter #"<<fMuonFilterVersion<<endl;
75  cout<<"ROOT Geometry, top volume name#"<<fTopVol->GetName()<<", top node name#"<<fTopNode->GetName()<<endl;
76  cout<<"State :"<<fStateTip<<endl;
77  cout<<(fGoodGeometry?"Succeed to" : "Fail to ")<<" load in "<<fNumofMdt<<" tubes, "<<fNumofStrip<<" strips."<<endl;
78  cout<<"========================================================================="<<endl;
79  }
80  }
81  return fGoodGeometry;
82 }
83 namespace shadow {
84  struct FullPathConverter
85  {
86  FullPathConverter(TString& str) :fullPath(str){}
87  void operator()(TString aStr) { fullPath.Append(aStr); fullPath.Append("/"); }
88  private:
89  TString& fullPath;
90  };
91 }
92 void PndMdtIGeometry::LoadGeometry(TGeoNode* node)
93 {
94  static std::list<TString> stackPath;
95  TString name(node->GetName());
96  stackPath.push_back(name);
97  fGoodGeometry = kTRUE;
98  if(node->GetNdaughters() > 0)
99  {
100  Int_t nd = node->GetNdaughters();
101  for (Int_t i=0; i<nd; i++) {
102  TGeoNode* theNode = node->GetDaughter(i);
103  if(fGoodGeometry)
104  LoadGeometry(theNode);
105  }
106  }else{
107  //if(name.Contains("MDT"))
108  if(FindSensorinPath(name))
109  {
110  ++fNumofMdt;
111  TString fFullPath;
112  for_each(stackPath.begin(), stackPath.end(), shadow::FullPathConverter(fFullPath));
113  fFullPath.Remove(fFullPath.Length()-1);
114  fGoodGeometry = GetGeometryInfoV2(fFullPath, stackPath.back());
115  }
116  }
117  stackPath.pop_back();
118 }
119 //analyis the bound of one layer
120 void PndMdtIGeometry::InsertLayerInfo(Int_t layerID, const TVector3& pos, const TVector3& lwh)
121 {
122  LayerInfoMapIter it = fLayerInfoMap.find(layerID);
123  TVector3 pos1 = pos - lwh;
124  TVector3 pos2 = pos + lwh;
125  if( it == fLayerInfoMap.end())
126  {
127  LayerBoundary aLayerBound;
128  aLayerBound.minX = pos1.X();
129  aLayerBound.minY = pos1.Y();
130  aLayerBound.minZ = pos1.Z();
131  aLayerBound.maxX = pos2.X();
132  aLayerBound.maxY = pos2.Y();
133  aLayerBound.maxZ = pos2.Z();
134  fLayerInfoMap.insert(LayerInfoMapValue(layerID, aLayerBound));
135  }else{
136  LayerBoundary& aLayerBound = it->second;
137  if(pos1.X() < aLayerBound.minX) aLayerBound.minX = pos1.X();
138  if(pos1.Y() < aLayerBound.minY) aLayerBound.minY = pos1.Y();
139  if(pos1.Z() < aLayerBound.minZ) aLayerBound.minZ = pos1.Z();
140  if(pos2.X() > aLayerBound.maxX) aLayerBound.maxX = pos2.X();
141  if(pos2.Y() > aLayerBound.maxY) aLayerBound.maxY = pos2.Y();
142  if(pos2.Z() > aLayerBound.maxZ) aLayerBound.maxZ = pos2.Z();
143  }
144 }
145 Bool_t PndMdtIGeometry::GetStripInfo()
146 {
147  static const Double_t cGAPOFBOXANDSTRIP = 0.1;//cent meter
148  static const Double_t cSTRIPWIDTH = 1.; //cent meter;
149  static const Double_t cSTRIPTHICKNESS = 0.20;//cent menter
150  static const Double_t cGAPOFSTRIPS = 0.0;//cent meter;
151 
152  fStateTip = " GetStripInfo:: place strips on ";
153  fStateTip += fLayerInfoMap.size();
154  fStateTip += " layers.";
155  LayerInfoMapIter it = fLayerInfoMap.begin();
156  LayerInfoMapIter end = fLayerInfoMap.end();
157  Double_t mStripZPos;
158  Double_t mStripXPos;
159  Double_t mStripLen;
160  Double_t mDiffX;
161  Double_t mDiffY;
162  Int_t stripNO = 0;
163  Bool_t isSucc = kTRUE;
164  for(; it != end; ++it)
165  {
166  stripNO = 0;
167  Int_t layerID = it->first;
168  LayerBoundary& theBoundary = it->second;
169  Int_t mod = PndMdtID::Module(layerID);
170  Int_t sec = PndMdtID::Sector(layerID);
171  Int_t layer = PndMdtID::Layer(layerID);
172  //put strip
173  if(mod == 1)
174  {//barrel, left angle = 0
175  Double_t zPos = theBoundary.minZ;
176  while( zPos < theBoundary.maxZ)
177  {
178  mStripZPos= zPos + 0.5*cSTRIPWIDTH;
179  zPos += (cSTRIPWIDTH + cGAPOFSTRIPS);
180  mDiffY = theBoundary.maxY - theBoundary.minY;
181  mDiffX = theBoundary.maxX - theBoundary.minX;
182  Double_t cx = 0.5*(theBoundary.minX + theBoundary.maxX);
183  Double_t cy = 0.5*(theBoundary.minY + theBoundary.maxY);
184  mStripLen = TMath::Abs(mDiffX*TMath::Cos(sec*TMath::Pi()/4.)) + TMath::Abs(mDiffY*TMath::Sin(sec*TMath::Pi()/4.));
185  Double_t dx = TMath::Abs((cGAPOFBOXANDSTRIP + 0.5*cSTRIPTHICKNESS + 0.5)*TMath::Cos(sec*TMath::Pi()/4.));
186  Double_t dy = TMath::Abs((cGAPOFBOXANDSTRIP + 0.5*cSTRIPTHICKNESS + 0.5)*TMath::Sin(sec*TMath::Pi()/4.));
187  isSucc = AddStripInfo(PndMdtID::Identifier(mod, sec, layer, stripNO), TVector3(cx+dx,cy+dy,mStripZPos), mStripLen);
188  if(!isSucc) return kFALSE;
189  ++ stripNO;
190  }
191  }else if(mod == 2)//endcap
192  {
193  if(layer < 3){//layer start from 1
194  Double_t xPos = theBoundary.minY;
195  while( xPos < theBoundary.maxY)
196  {
197  mStripXPos= xPos + 0.5*cSTRIPWIDTH;
198  xPos += (cSTRIPWIDTH + cGAPOFSTRIPS);
199  mStripLen = theBoundary.maxX - theBoundary.minX;
200  Double_t cx = 0.5*(theBoundary.minX + theBoundary.maxX);
201  Double_t cz = 0.5*(theBoundary.minZ + theBoundary.maxZ);
202  Double_t dz = (cGAPOFBOXANDSTRIP + 0.5*cSTRIPTHICKNESS + 0.5);
203  isSucc = AddStripInfo(PndMdtID::Identifier(mod, sec, layer, stripNO), TVector3(cx, mStripXPos, cz+dz), mStripLen);
204  if(!isSucc) return kFALSE;
205  ++ stripNO;
206  }
207  }else{
208  Double_t xPos = theBoundary.minX;
209  while( xPos < theBoundary.maxX)
210  {
211  mStripXPos= xPos + 0.5*cSTRIPWIDTH;
212  xPos += (cSTRIPWIDTH + cGAPOFSTRIPS);
213  mStripLen = theBoundary.maxY - theBoundary.minY;
214  Double_t cy = 0.5*(theBoundary.minY + theBoundary.maxY);
215  Double_t cz = 0.5*(theBoundary.minZ + theBoundary.maxZ);
216  Double_t dz = (cGAPOFBOXANDSTRIP + 0.5*cSTRIPTHICKNESS + 0.5);
217  isSucc = AddStripInfo(PndMdtID::Identifier(mod, sec, layer, stripNO), TVector3(mStripXPos, cy, cz+dz), mStripLen);
218  if(!isSucc) return kFALSE;
219  ++ stripNO;
220  }
221  }
222  }
223  else
224  {
225  Double_t xPos = theBoundary.minX;
226  while( xPos < theBoundary.maxX)
227  {
228  mStripXPos= xPos + 0.5*cSTRIPWIDTH;
229  xPos += (cSTRIPWIDTH + cGAPOFSTRIPS);
230  mStripLen = theBoundary.maxY - theBoundary.minY;
231  Double_t cy = 0.5*(theBoundary.minY + theBoundary.maxY);
232  Double_t cz = 0.5*(theBoundary.minZ + theBoundary.maxZ);
233  Double_t dz = (cGAPOFBOXANDSTRIP + 0.5*cSTRIPTHICKNESS + 0.5);
234  isSucc = AddStripInfo(PndMdtID::Identifier(mod, sec, layer, stripNO), TVector3(mStripXPos, cy, cz+dz), mStripLen);
235  if(!isSucc) return kFALSE;
236  ++ stripNO;
237  }
238  }
239  if(fVerbose>2)
240  cout<<"layer ID #"<<layerID<<", mod#"<<mod<<", sec#"<<sec<<", layer#"<<layer
241  <<", X#"<<theBoundary.minX<<", "<<theBoundary.maxX
242  <<", Y#"<<theBoundary.minY<<", "<<theBoundary.maxY
243  <<", Z#"<<theBoundary.minZ<<", "<<theBoundary.maxZ<<", "<<stripNO<<" strips."<<endl;
244  fNumofStrip += stripNO;
245  }
246  return kTRUE;
247 }
248 Bool_t PndMdtIGeometry::GetGeometryInfoV2(const TString& fFullPath, const TString& Name)
249 {
250  Int_t iMod(-1);
251  Int_t iSec(-1);
252  Int_t iSec1(-1);//for barrel special section
253  Int_t iLayer(-1);
254  Int_t iLayer1(-1);
255  Int_t iBox(-1);
256  Int_t iWire(-1);
257  Bool_t isSucc;
258  //read this tube information
259  isSucc = gGeoManager->cd(fFullPath.Data());
260  if(isSucc){
261  Int_t iholder;
262  if(fFullPath.BeginsWith("cave_1/Mdt_1/")){
263  fBarrelVersion = fEndcapVersion = fMuonFilterVersion = fForwardVersion = 0;
264  sscanf(Name.Data(),"MDT%is%il%ib%iw%i_%i", &iMod, &iSec, &iLayer, &iBox, &iWire, &iWire);
265  }
266  else if(fFullPath.BeginsWith("cave_1/Mdt")){
267  if(fFullPath.Contains("MdtBarrel")){//
268  fBarrelVersion = 1;
269  sscanf(fFullPath.Data(),"cave_1/MdtBarrel_0/MdtBarrelOct%i_%i/MdtBarrelOct%iLayer%i_%i/BP%is%il%ib%iw%i_%i/BA%i_%i/MDT%is%il%ib%iw%i_%i"
270  , &iSec1, &iSec, &iholder
271  , &iLayer, &iLayer1
272  , &iMod, &iholder, &iholder, &iBox, &iholder, &iholder
273  , &iholder, &iholder
274  , &iholder, &iholder, &iholder, &iholder, &iholder, &iWire);
275  }else if(fFullPath.Contains("MdtEndcap")){
276  fEndcapVersion = 1;
277  sscanf(fFullPath.Data(),"cave_1/MdtEndcap_0/MdtEndcapLayer%i_%i/BP%is%il%ib%iw%i_%i/BA%i_%i/MDT%is%il%ib%iw%i_%i"
278  , &iLayer1, &iLayer //MdtEndcapLayer ...
279  , &iMod, &iSec, &iholder, &iholder, &iholder, &iBox//BP...
280  , &iholder, &iholder //BA ...
281  , &iholder, &iholder, &iholder, &iholder, &iholder, &iWire);//MDT ...
282  }else if(fFullPath.Contains("MdtMF")){
283  fMuonFilterVersion = 1;
284  sscanf(fFullPath.Data(),"cave_1/MdtMF_0/MdtMFLayer%i_%i/BP%is%il%ib%iw%i_%i/BA%i_%i/MDT%is%il%ib%iw%i_%i"
285  , &iLayer1, &iLayer
286  , &iMod, &iSec, &iholder, &iBox, &iholder, &iholder
287  , &iholder, &iholder
288  , &iholder, &iholder, &iholder, &iholder, &iholder, &iWire);
289  }else if(fFullPath.Contains("MdtForward")){
290  fForwardVersion = 1;
291  sscanf(fFullPath.Data(),"cave_1/MdtForward_0/MdtForwardLayer%i_%i/BP%is%il%ib%iw%i_%i/BA%i_%i/MDT%is%il%ib%iw%i_%i"
292  , &iLayer1, &iLayer
293  , &iMod, &iSec, &iholder, &iBox, &iholder, &iholder
294  , &iholder, &iholder
295  , &iholder, &iholder, &iholder, &iholder, &iholder, &iWire);
296  }else{
297  cout<<"PndMdtIGeometry::GetGeometryInfoV2, Unknown geometry pattern."<<endl;
298  }
299  }
300  else if(fFullPath.BeginsWith("cave_1/Mdt")){
301  fBarrelVersion = fEndcapVersion = fMuonFilterVersion = fForwardVersion = 0;
302  sscanf(Name.Data(),"MDT%is%il%ib%iw%i_%i", &iMod, &iSec, &iLayer, &iBox, &iWire, &iWire);
303  }else {
304  cout<<"PndMdtIGeometry::GetGeometryInfoV2, Unknown geometry pattern."<<endl;
305  }
306  TVector3 itsPos(0,0,0);
307  TGeoBBox* theShape = (TGeoBBox*)gGeoManager->GetCurrentVolume()->GetShape();
308  TVector3 itsOrigin(theShape->GetOrigin());
309  gGeoManager->LocalToMaster(&itsOrigin[0], &itsPos[0]);
310  TVector3 itsLWH(theShape->GetDX(), theShape->GetDY(), theShape->GetDZ());
311  Double_t wirelen ;
312  if( iMod == 1 ) wirelen = 2.*itsLWH.Z() ;
313  else if(iMod == 2 && iLayer < 2) wirelen = 2.*itsLWH.X();
314  else wirelen = 2.*itsLWH.Y();
315 
316  if(fVerbose > 2){
317  cout<<fFullPath<<endl;
318  cout<<"("<<iMod<<", "<<iSec<<", "<<iLayer<<", "<<iBox<<", "<<iWire<<")" <<itsPos<<", "<<itsLWH<<endl;
319  }
320 
321  TGeoMatrix* matrix = gGeoManager->GetCurrentMatrix()->MakeClone();
322  Int_t detID = PndMdtID::Identifier(iMod, iSec, iLayer, iBox, iWire);
323  isSucc = AddTubeInfo(detID, itsPos, wirelen, matrix);
324  //
325  Int_t layerID = PndMdtID::LayerID(iMod, iSec, iLayer);
326  InsertLayerInfo(layerID, itsPos, itsLWH);
327 
328  AddMatrix(layerID, matrix);
329  }else{
330  fStateTip = "gGeoManager:: nonexisted path";
331  cerr<< "PndMdtIGeometry::GetGeometryInfoV2, "<<fStateTip<<endl;;
332  }
333  return isSucc;
334 }
335 Bool_t PndMdtIGeometry::AddMatrix(Int_t layerID, TGeoMatrix* matrix)
336 {
337  if(fMatrixMap.find(layerID) == fMatrixMap.end()){
338  if(fVerbose>1)
339  cout <<"PndMdtIGeometry::AddMatrix(), MDT layer ID#"<<layerID<<" insert"<<endl;
340  fMatrixMap.insert(std::pair<Int_t, TGeoMatrix*>(layerID, matrix));
341  }
342  return kTRUE;
343 }
344 Bool_t PndMdtIGeometry::AddTubeInfo(Int_t detID, const TVector3& center, Double_t wirelen, TGeoMatrix* matrix)
345 {
346  if(fGeoMap.find(detID) != fGeoMap.end()){
347  fStateTip = "AddTubeInfo:: detID not unique";
348  cout << "PndMdtIGeometry::"<<fStateTip<<", "<<detID<<endl;;
349  return kFALSE;
350  }else{
351  InfoType _info;
352  _info.Position = center;
353  if(TMath::Abs(_info.Position.X()) < 1.e-8) _info.Position.SetX(0.);
354  if(TMath::Abs(_info.Position.Y()) < 1.e-8) _info.Position.SetY(0.);
355  if(TMath::Abs(_info.Position.Z()) < 1.e-8) _info.Position.SetZ(0.);
356  _info.Length = wirelen;
357  _info.Matrix = matrix;
358  if(fVerbose>1)
359  cout <<"PndMdtIGeometry::AddTubeInfo(), MDT detector ID#"<<detID<<" insert"<<endl;
360  fGeoMap.insert(std::pair<Int_t, InfoType>(detID, _info));
361  }
362  return kTRUE;
363 }
364 Bool_t PndMdtIGeometry::AddStripInfo(Int_t detID, TVector3 stripPos, Double_t stripLen, TGeoMatrix* matrix)
365 {
366  if(fGeoMap.find(detID) != fGeoMap.end()){
367  fStateTip = "AddStripInfo:: detID not unique";
368  cout << "PndMdtIGeometry::"<<fStateTip<<", "<<detID<<endl;;
369  return kFALSE;
370  }else{
371  InfoType _info;
372  _info.Position = stripPos;
373  _info.Length = stripLen;
374  _info.Matrix = matrix;
375  if(fVerbose>1)
376  cout <<"PndMdtIGeometry::AddStripInfo(), MDT detector ID#"<<detID<<" insert"<<endl;
377  fGeoMap.insert(std::pair<Int_t, InfoType>(detID, _info));
378  }
379  return kTRUE;
380 }
381 void PndMdtIGeometry::Print() const
382 {
383  std::map<Int_t, InfoType>::const_iterator cit = fGeoMap.begin();;
384  std::map<Int_t, InfoType>::const_iterator cend = fGeoMap.end();;
385 
386  std::ofstream ofs1("Tube.map", std::ios::out);
387  ofs1<<"============================================================"<<endl;
388  ofs1<<"==>Tube Placement"<<endl;
389  ofs1<<"============================================================"<<endl;
390  ofs1<<"Mod\t"<<"Sec\t"<<"Layer\t"<<"Box\t"<<"NO.\t"<<"X\t"<<"Y\t"<<"Z\t"<<endl;
391  for(; cit != cend; ++cit){
392  if(PndMdtID::isWire(cit->first))
393  {
394  ofs1<<PndMdtID::Module(cit->first)<<"\t"
395  <<PndMdtID::Sector(cit->first)<<"\t"
396  <<PndMdtID::Layer(cit->first)<<"\t"
397  <<PndMdtID::Box(cit->first)<<"\t"
398  <<PndMdtID::Wire(cit->first)<<"\t"
399  <<cit->second.Position.X()<<"\t"
400  <<cit->second.Position.Y()<<"\t"
401  <<cit->second.Position.Z()<<"\t"<<endl;
402  }
403  }
404  ofs1.close();
405  std::ofstream ofs2("Strip.map", std::ios::out);
406  ofs2<<"============================================================"<<endl;
407  ofs2<<"==>Strip Placement"<<endl;
408  ofs2<<"============================================================"<<endl;
409  ofs2<<"Mod\t"<<"Sec\t"<<"Layer\t"<<"NO.\t"<<"X\t"<<"Y\t"<<"Z\t"<<endl;
410  cit = fGeoMap.begin();;
411  for(; cit != cend; ++cit){
412  if(PndMdtID::isStrip(cit->first))
413  {
414  ofs2<<PndMdtID::Module(cit->first)<<"\t"
415  <<PndMdtID::Sector(cit->first)<<"\t"
416  <<PndMdtID::Layer(cit->first)<<"\t"
417  <<PndMdtID::Strip(cit->first)<<"\t"
418  <<cit->second.Position.X()<<"\t"
419  <<cit->second.Position.Y()<<"\t"
420  <<cit->second.Position.Z()<<"\t"<<endl;
421  }
422  }
423  ofs2.close();
424 }
425 
426 ClassImp(PndMdtIGeometry)
pos
TVector3 pos
Definition: GammaSpectraAnalysis.C:43
PndMdtID::Box
static Short_t Box(Int_t detID)
Definition: PndMdtID.h:23
dy
double dy
Definition: runGeaneTestSimBox.cxx:45
PndMdtID::isWire
static Bool_t isWire(Int_t detID)
Definition: PndMdtID.h:26
PndMdtID::Module
static Short_t Module(Int_t detID)
Definition: PndMdtID.h:20
shadow::FullPathConverter::operator()
void operator()(TString aStr)
Definition: PndMdtIGeometry.cxx:87
PndMdtID::Layer
static Short_t Layer(Int_t detID)
Definition: PndMdtID.h:22
i
Int_t i
Definition: run_full.C:25
PndMdtIGeometry::fGeoMap
std::map< Int_t, InfoType > fGeoMap
Definition: PndMdtIGeometry.h:143
CAMath::Sin
static T Sin(const T &x)
Definition: PndCAMath.h:42
PndMdtIGeometry::AddStripInfo
Bool_t AddStripInfo(Int_t detID, TVector3 stripPos, Double_t stripLen, TGeoMatrix *matrix=0)
Definition: PndMdtIGeometry.cxx:364
operator<<
std::ostream & operator<<(std::ostream &o, const PndEventInfo &a)
Definition: PndEventInfo.cxx:68
PndMdtID::LayerID
static Int_t LayerID(Int_t iMod, Int_t iOct, Int_t iLayer)
Definition: PndMdtID.h:16
dz
Double_t dz
Definition: createRootGeoFileBwEnd_FwEndCapGeo.C:24
gGeoManager
TGeoManager * gGeoManager
Definition: createRootGeoFile.C:97
CAMath::Cos
static T Cos(const T &x)
Definition: PndCAMath.h:43
PndMdtIGeometry::fInstance
static PndMdtIGeometry * fInstance
Definition: PndMdtIGeometry.h:179
PndMdtIGeometry::AddTubeInfo
Bool_t AddTubeInfo(Int_t detID, const TVector3 &center, Double_t wirelen, TGeoMatrix *matrix)
Definition: PndMdtIGeometry.cxx:344
PndMdtID::Identifier
static Int_t Identifier(Int_t iMod, Int_t iOct, Int_t iLayer, Int_t iBox, Int_t iWire)
Definition: PndMdtID.h:10
CAMath::Abs
static T Abs(const T &x)
Definition: PndCAMath.h:39
PndMdtID::Sector
static Short_t Sector(Int_t detID)
Definition: PndMdtID.h:21
std::ostream
basic_ostream< char, char_traits< char > > ostream
Definition: PndFTSCATrackParamVector.h:1258
Double_t
Double_t
Definition: SimCompleteLinkDef.h:6
PndMdtIGeometry::LoadGeometry
void LoadGeometry(TGeoNode *node)
Definition: PndMdtIGeometry.cxx:92
PndMdtIGeometry::FindSensorinPath
Bool_t FindSensorinPath(TString path)
Definition: PndMdtIGeometry.h:119
out
TFile * out
Definition: reco_muo.C:20
PndMdtID::Strip
static Short_t Strip(Int_t detID)
Definition: PndMdtID.h:25
name
TString name
Definition: createRootGeoFile.C:82
PndMdtIGeometry::InsertLayerInfo
void InsertLayerInfo(Int_t layerID, const TVector3 &pos, const TVector3 &lwh)
Definition: PndMdtIGeometry.cxx:120
dx
double dx
Definition: runGeaneTestSimBox.cxx:44
PndMdtIGeometry::LayerInfoMapValue
std::pair< Int_t, LayerBoundary > LayerInfoMapValue
Definition: PndMdtIGeometry.h:159
PndMdtID::isStrip
static Bool_t isStrip(Int_t detID)
Definition: PndMdtID.h:27
cy
int cy
Definition: dedx_bands.C:17
PndMdtIGeometry::Instance
static PndMdtIGeometry * Instance()
Definition: PndMdtIGeometry.cxx:30
PndMdtIGeometry::LayerInfoMapIter
std::map< Int_t, LayerBoundary >::iterator LayerInfoMapIter
Definition: PndMdtIGeometry.h:158
PndMdtID::Wire
static Short_t Wire(Int_t detID)
Definition: PndMdtID.h:24
ClassImp
ClassImp(PndAnaContFact)
layer
Int_t layer
Definition: reco_muo.C:36
cx
int cx
Definition: dedx_bands.C:16
PndMdtIGeometry::fMatrixMap
std::map< Int_t, TGeoMatrix * > fMatrixMap
Definition: PndMdtIGeometry.h:144
Pi
Double_t Pi
Definition: PndGeoHypGeCrystal.cxx:10
PndMdtIGeometry::AddMatrix
Bool_t AddMatrix(Int_t layerID, TGeoMatrix *matrix)
Definition: PndMdtIGeometry.cxx:335
PndMdtIGeometry::fLayerInfoMap
std::map< Int_t, LayerBoundary > fLayerInfoMap
Definition: PndMdtIGeometry.h:157
PndMdtIGeometry::Print
void Print() const
Definition: PndMdtIGeometry.cxx:381
PndMdtIGeometry::GetGeometryInfoV2
Bool_t GetGeometryInfoV2(const TString &fFullPath, const TString &Name)
Definition: PndMdtIGeometry.cxx:248
Generated on Wed Apr 3 2019 08:43:58 for FairRoot/PandaRoot by   doxygen 1.8.5