1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
19 #include "AliITSsegmentationSDD.h"
21 #include "AliITSgeom.h"
23 #include "AliITSresponse.h"
25 ClassImp(AliITSsegmentationSDD)
26 //----------------------------------------------------------------------
27 AliITSsegmentationSDD::AliITSsegmentationSDD(AliITSgeom* geom,
28 AliITSresponse *resp){
38 //______________________________________________________________________
39 AliITSsegmentationSDD::AliITSsegmentationSDD(){
40 // standard constructor
49 //______________________________________________________________________
50 AliITSsegmentationSDD& AliITSsegmentationSDD::operator=(AliITSsegmentationSDD
53 if(this==&source) return *this;
54 this->fNsamples = source.fNsamples;
55 this->fNanodes = source.fNanodes;
56 this->fPitch = source.fPitch;
57 this->fTimeStep = source.fTimeStep;
58 this->fDx = source.fDx;
59 this->fDz = source.fDz;
60 this->fDy = source.fDy;
61 this->fCorr = new TF1(*(source.fCorr));
62 this->fGeom = source.fGeom; // Just copy the pointer
63 this->fResponse = source.fResponse; //Just copy the pointer
66 //______________________________________________________________________
67 AliITSsegmentationSDD::AliITSsegmentationSDD(AliITSsegmentationSDD &source){
71 //----------------------------------------------------------------------
72 void AliITSsegmentationSDD::Init(){
73 // Standard initilisation routine
77 //fGeom = ((AliITS*)gAlice->GetModule("ITS"))->GetITSgeom();
79 AliITSgeomSDD *gsdd = (AliITSgeomSDD *) (fGeom->GetShape(3,1,1));
81 const Float_t kconv=10000.;
82 fDz = 2.*kconv*gsdd->GetDz();
83 fDx = kconv*gsdd->GetDx();
84 fDy = 2.*kconv*gsdd->GetDy();
87 //----------------------------------------------------------------------
88 void AliITSsegmentationSDD::
89 Neighbours(Int_t iX, Int_t iZ, Int_t* Nlist, Int_t Xlist[8], Int_t Zlist[8]){
90 // returns neighbours for use in Cluster Finder routines and the like
92 if(iX >= fNanodes) printf("iX > fNanodes %d %d\n",iX,fNanodes);
93 if(iZ >= fNsamples) printf("iZ > fNsamples %d %d\n",iZ,fNsamples);
96 if(iX && (iX != fNanodes/2)) Xlist[2]=iX-1;
98 if ((iX !=fNanodes/2 -1) && (iX != fNanodes)) Xlist[3]=iX+1;
100 if(iZ) Zlist[0]=iZ-1;
102 if (iZ < fNsamples) Zlist[1]=iZ+1;
104 Zlist[2]=Zlist[3]=iZ;
106 //----------------------------------------------------------------------
107 void AliITSsegmentationSDD::GetPadIxz(Float_t x,Float_t z,
108 Int_t &timebin,Int_t &anode){
109 // Returns cell coordinates (time sample,anode) for given real local
112 // expects x, z in cm
114 const Float_t kconv=10000; // cm->um
116 Float_t speed=fResponse->DriftSpeed();
117 Int_t na = fNanodes/2;
118 Float_t driftpath=fDx-TMath::Abs(kconv*x);
119 timebin=(Int_t)(driftpath/speed/fTimeStep);
120 anode=(Int_t)(kconv*z/fPitch) + na/2;
121 if (x > 0) anode += na;
127 //----------------------------------------------------------------------
128 void AliITSsegmentationSDD::GetPadCxz(Int_t timebin,Int_t anode,
129 Float_t &x ,Float_t &z){
130 // Transform from cell to real local coordinates
131 // returns x, z in cm
133 const Float_t kconv=10000; // um->cm
135 Float_t speed=fResponse->DriftSpeed();
136 Int_t na = fNanodes/2;
137 Float_t driftpath=(timebin+1)*fTimeStep*speed;
138 if (anode >= na) x=(fDx-driftpath)/kconv;
139 else x = -(fDx-driftpath)/kconv;
140 if (anode >= na) anode-=na;
141 z=((anode+1)*fPitch-fDz/2)/kconv;
144 //----------------------------------------------------------------------
145 void AliITSsegmentationSDD::GetPadTxz(Float_t &x,Float_t &z){
146 // Get anode and time bucket as floats - numbering from 0
148 // expects x, z in cm
150 const Float_t kconv=10000; // cm->um
153 Float_t speed=fResponse->DriftSpeed();
154 Int_t na = fNanodes/2;
155 Float_t driftpath=fDx-TMath::Abs(kconv*x);
156 x=driftpath/speed/fTimeStep;
157 z=kconv*z/fPitch + (float)na/2;
161 //----------------------------------------------------------------------
162 void AliITSsegmentationSDD::GetLocal(Int_t module,Float_t *g ,Float_t *l){
163 // returns local coordinates from global
166 //fGeom = ((AliITS*)gAlice->GetModule("ITS"))->GetITSgeom();
168 fGeom->GtoL(module,g,l);
170 //----------------------------------------------------------------------
171 void AliITSsegmentationSDD::GetGlobal(Int_t module,Float_t *l ,Float_t *g){
172 // return global coordinates from local
175 //fGeom = ((AliITS*)gAlice->GetModule("ITS"))->GetITSgeom();
178 fGeom->LtoG(module,l,g);