1 #include <TGeoVolume.h>
3 #include <TGeoManager.h>
4 #include "AliITSURecoDet.h"
5 #include "AliITSUGeomTGeo.h"
6 #include "AliITSsegmentation.h"
7 #include "AliITSUSegmentationPix.h"
8 #include "AliITSUReconstructor.h"
11 ClassImp(AliITSURecoDet)
14 const Char_t* AliITSURecoDet::fgkBeamPipeVolName = "IP_PIPE";
17 //______________________________________________________
18 AliITSURecoDet::AliITSURecoDet(AliITSUGeomTGeo* geom, const char* name)
28 SetNameTitle(name,name);
29 fLayers.SetOwner(kTRUE); // layers belong to this array
30 fLayersActive.SetOwner(kFALSE); // this one just points on active layers in fLayers
34 //______________________________________________________
35 AliITSURecoDet::~AliITSURecoDet()
38 fLayersActive.Clear();
39 fLayers.Clear(); // owned!
42 //______________________________________________________
43 void AliITSURecoDet::Print(Option_t* opt) const
46 printf("Detector %s, %d layers, %d active layers\n",GetName(),GetNLayers(),GetNLayersActive());
47 TString opts = opt; opts.ToLower();
48 if (opts.Contains("lr")) for (int i=0;i<GetNLayers();i++) GetLayer(i)->Print(opt);
51 //______________________________________________________
52 void AliITSURecoDet::AddLayer(const AliITSURecoLayer* lr)
55 fLayers.AddLast((TObject*)lr);
58 fLayersActive.AddLast((TObject*)lr);
63 //______________________________________________________
64 Bool_t AliITSURecoDet::Build()
66 // build detector from TGeo
68 if (!fGeom) AliFatal("Geometry interface is not set");
69 int nlr = fGeom->GetNLayers();
70 if (!nlr) AliFatal("No geometry loaded");
72 // build active ITS layers
73 for (int ilr=0;ilr<nlr;ilr++) {
74 int lrTyp = fGeom->GetLayerDetTypeID(ilr);
75 // name layer according its active id, detector type and segmentation tyoe
76 AliITSURecoLayer* lra = new AliITSURecoLayer(Form("Lr%d%s%d",ilr,fGeom->GetDetTypeName(lrTyp),
77 lrTyp%AliITSUGeomTGeo::kMaxSegmPerDetType),
79 lra->SetPassive(kFALSE);
83 // build passive ITS layers
85 double rMin,rMax,zMin,zMax;
87 TGeoVolume *v = gGeoManager->GetVolume(fgkBeamPipeVolName);
88 AliITSURecoLayer* lrp = 0;
89 if (!v) AliWarning("No beam pipe found in geometry");
91 TGeoTube *t=(TGeoTube*)v->GetShape();
96 lrp = new AliITSURecoLayer("BeamPipe");
99 lrp->SetR(0.5*(rMin+rMax));
102 lrp->SetPassive(kTRUE);
108 lrp = new AliITSURecoLayer("TPC-ITSwall");
109 lrp->SetRMin(AliITSUReconstructor::GetRecoParam()->GetTPCITSWallRMin());
110 lrp->SetRMax(AliITSUReconstructor::GetRecoParam()->GetTPCITSWallRMax());
111 lrp->SetR(0.5*(lrp->GetRMin()+lrp->GetRMax()));
112 lrp->SetZMin(-AliITSUReconstructor::GetRecoParam()->GetTPCITSWallZSpanH());
113 lrp->SetZMax( AliITSUReconstructor::GetRecoParam()->GetTPCITSWallZSpanH());
114 lrp->SetMaxStep( AliITSUReconstructor::GetRecoParam()->GetTPCITSWallMaxStep());
115 lrp->SetPassive(kTRUE);
123 //______________________________________________________
124 void AliITSURecoDet::IndexLayers()
126 // sort and index layers
127 const Double_t kRMargin = 1e-2; // 100 micron margin
128 fLayersActive.Sort();
129 for (int i=0;i<fNLayersActive;i++) GetLayerActive(i)->SetActiveID(i);
131 for (int i=0;i<fNLayers;i++) GetLayer(i)->SetID(i);
133 SetRMin(GetLayer(0)->GetRMin()-kRMargin);
134 SetRMax(GetLayer(fNLayers-1)->GetRMax()+kRMargin);
138 //______________________________________________________
139 Int_t AliITSURecoDet::FindLastLayerID(Double_t r, int dir) const
141 // find the last layer which the particle moving in direction dir (1:outward,-1:inward)
142 // will traverse on its way to radius r
146 for (ilr=0;ilr<fNLayers;ilr++) {
147 AliITSURecoLayer* lr = GetLayer(ilr);
148 if ( r<lr->GetR(-dir) ) break; // this layer at least entered
150 return --ilr; // -1 will correspond to point below the smalles layer
153 for (ilr=fNLayers;ilr--;) {
154 AliITSURecoLayer* lr = GetLayer(ilr);
155 if ( r>lr->GetR(-dir) ) break; // this layer at least entered
158 return ilr<fNLayers ? ilr:-1; // -1 will correspond to point above outer layer
163 //______________________________________________________
164 Int_t AliITSURecoDet::FindFirstLayerID(Double_t r, int dir) const
166 // find the first layer which the particle moving in direction dir (1:outward,-1:inward)
167 // will traverse starting from radius r
171 for (ilr=0;ilr<fNLayers;ilr++) {
172 AliITSURecoLayer* lr = GetLayer(ilr);
173 if ( r<lr->GetR(dir) ) break; // this layer at least entered
175 return ilr<fNLayers ? ilr:-1; // -1 will correspond to point above outer leayer
178 for (ilr=fNLayers;ilr--;) {
179 AliITSURecoLayer* lr = GetLayer(ilr);
180 if ( r>lr->GetR(dir) ) break; // this layer at least entered
182 return ilr; // -1 will correspond to point below inner layer