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 **************************************************************************/
18 #include "AliITSsegmentationSDD.h"
20 #include "AliITSgeom.h"
22 #include "AliITSresponse.h"
25 ClassImp(AliITSsegmentationSDD)
26 //------------------------------
27 AliITSsegmentationSDD::AliITSsegmentationSDD(AliITSgeom* geom, AliITSresponse *resp){
37 //_____________________________________________________________________________
38 AliITSsegmentationSDD::AliITSsegmentationSDD(){
39 // standard constructor
48 //_____________________________________________________________________________
49 AliITSsegmentationSDD& AliITSsegmentationSDD::operator=(AliITSsegmentationSDD &source){
51 if(this==&source) return *this;
52 this->fNsamples = source.fNsamples;
53 this->fNanodes = source.fNanodes;
54 this->fPitch = source.fPitch;
55 this->fTimeStep = source.fTimeStep;
56 this->fDx = source.fDx;
57 this->fDz = source.fDz;
58 this->fDy = source.fDy;
59 this->fCorr = new TF1(*(source.fCorr));
60 this->fGeom = source.fGeom; // Just copy the pointer
61 this->fResponse = source.fResponse; //Just copy the pointer
64 //___________________________________________________________________________
65 AliITSsegmentationSDD::AliITSsegmentationSDD(AliITSsegmentationSDD &source){
69 //------------------------------
70 void AliITSsegmentationSDD::Init(){
71 // Standard initilisation routine
75 //fGeom = ((AliITS*)gAlice->GetModule("ITS"))->GetITSgeom();
77 AliITSgeomSDD *gsdd = (AliITSgeomSDD *) (fGeom->GetShape(3,1,1));
79 const Float_t kconv=10000.;
80 fDz = 2.*kconv*gsdd->GetDz();
81 fDx = kconv*gsdd->GetDx();
82 fDy = 2.*kconv*gsdd->GetDy();
85 //------------------------------
86 void AliITSsegmentationSDD::
87 Neighbours(Int_t iX, Int_t iZ, Int_t* Nlist, Int_t Xlist[8], Int_t Zlist[8]){
88 // returns neighbours for use in Cluster Finder routines and the like
90 if(iX >= fNanodes) printf("iX > fNanodes %d %d\n",iX,fNanodes);
91 if(iZ >= fNsamples) printf("iZ > fNsamples %d %d\n",iZ,fNsamples);
94 if(iX && (iX != fNanodes/2)) Xlist[2]=iX-1;
96 if ((iX !=fNanodes/2 -1) && (iX != fNanodes)) Xlist[3]=iX+1;
100 if (iZ < fNsamples) Zlist[1]=iZ+1;
102 Zlist[2]=Zlist[3]=iZ;
104 //------------------------------
105 void AliITSsegmentationSDD::GetPadIxz(Float_t x,Float_t z,Int_t &timebin,Int_t &anode){
106 // Returns cell coordinates (time sample,anode) for given real local coordinates (x,z)
108 // expects x, z in cm
110 const Float_t kconv=10000; // cm->um
112 Float_t speed=fResponse->DriftSpeed();
113 Int_t na = fNanodes/2;
114 Float_t driftpath=fDx-TMath::Abs(kconv*x);
115 timebin=(Int_t)(driftpath/speed/fTimeStep);
116 anode=(Int_t)(kconv*z/fPitch) + na/2;
117 if (x > 0) anode += na;
124 //------------------------------
125 void AliITSsegmentationSDD::GetPadCxz(Int_t timebin,Int_t anode,Float_t &x ,Float_t &z){
126 // Transform from cell to real local coordinates
128 // returns x, z in cm
130 const Float_t kconv=10000; // um->cm
132 Float_t speed=fResponse->DriftSpeed();
133 Int_t na = fNanodes/2;
134 Float_t driftpath=(timebin+1)*fTimeStep*speed;
135 if (anode >= na) x=(fDx-driftpath)/kconv;
136 else x = -(fDx-driftpath)/kconv;
137 if (anode >= na) anode-=na;
138 z=((anode+1)*fPitch-fDz/2)/kconv;
142 //------------------------------
143 void AliITSsegmentationSDD::GetPadTxz(Float_t &x,Float_t &z){
144 // Get anode and time bucket as floats - numbering from 0
146 // expects x, z in cm
148 const Float_t kconv=10000; // cm->um
151 Float_t speed=fResponse->DriftSpeed();
152 //Int_t na = fNanodes/2;
153 Float_t driftpath=fDx-TMath::Abs(kconv*x);
154 x=driftpath/speed/fTimeStep;
156 // z=kconv*z/fPitch + (float)na/2;
157 //if (x0 > 0) z += (float)na;
160 //------------------------------
161 void AliITSsegmentationSDD::GetLocal(Int_t module,Float_t *g ,Float_t *l){
162 // returns local coordinates from global
165 //fGeom = ((AliITS*)gAlice->GetModule("ITS"))->GetITSgeom();
167 fGeom->GtoL(module,g,l);
169 //------------------------------
170 void AliITSsegmentationSDD::GetGlobal(Int_t module,Float_t *l ,Float_t *g){
171 // return global coordinates from local
174 //fGeom = ((AliITS*)gAlice->GetModule("ITS"))->GetITSgeom();
177 fGeom->LtoG(module,l,g);