X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ITS%2FAliITSsegmentationSDD.cxx;h=2821eccd0232115d08a957c67f0a0b8d113b2cf0;hb=a24be56b08ef731a32c600e49534f4c4beb9370c;hp=7dd4e2d5cf9f125fab62f0fde18c8cc2fdc56c4f;hpb=1ca7869b62ab56172aa60913e439a5ee237757ff;p=u%2Fmrichter%2FAliRoot.git diff --git a/ITS/AliITSsegmentationSDD.cxx b/ITS/AliITSsegmentationSDD.cxx index 7dd4e2d5cf9..2821eccd023 100644 --- a/ITS/AliITSsegmentationSDD.cxx +++ b/ITS/AliITSsegmentationSDD.cxx @@ -12,68 +12,102 @@ * about the suitability of this software for any purpose. It is * * provided "as is" without express or implied warranty. * **************************************************************************/ - -#include +#include #include #include "AliITSsegmentationSDD.h" -#include "AliITS.h" +// #include "AliITS.h" #include "AliITSgeom.h" -#include "AliRun.h" -#include "AliITSresponse.h" +#include "AliITSgeomSDD.h" +#include "AliITSCalibration.h" +#include "AliITSCalibrationSDD.h" +#include "AliITSresponseSDD.h" +////////////////////////////////////////////////////// +// Segmentation class for // +// drift detectors // +// // +////////////////////////////////////////////////////// +const Float_t AliITSsegmentationSDD::fgkDxDefault = 35085.; +const Float_t AliITSsegmentationSDD::fgkDzDefault = 75264.; +const Float_t AliITSsegmentationSDD::fgkDyDefault = 300.; +const Float_t AliITSsegmentationSDD::fgkPitchDefault = 294.; +const Float_t AliITSsegmentationSDD::fgkClockDefault = 40.; +const Int_t AliITSsegmentationSDD::fgkHalfNanodesDefault = 256; +const Int_t AliITSsegmentationSDD::fgkNsamplesDefault = 256; ClassImp(AliITSsegmentationSDD) -//------------------------------ -AliITSsegmentationSDD::AliITSsegmentationSDD(AliITSgeom* geom, AliITSresponse *resp){ +//---------------------------------------------------------------------- +AliITSsegmentationSDD::AliITSsegmentationSDD(AliITSgeom* geom, + AliITSCalibration *resp): +fNsamples(0), +fNanodes(0), +fPitch(0), +fTimeStep(0), +fDriftSpeed(0){ // constructor fGeom=geom; - fResponse=resp; + AliITSCalibrationSDD* sp = (AliITSCalibrationSDD*)resp; + fDriftSpeed=sp->GetDriftSpeed(); fCorr=0; - SetDetSize(); - SetPadSize(); - SetNPads(); + SetDetSize(fgkDxDefault,fgkDzDefault,fgkDyDefault); + SetPadSize(fgkPitchDefault,fgkClockDefault); + SetNPads(fgkHalfNanodesDefault,fgkNsamplesDefault); } -//_____________________________________________________________________________ -AliITSsegmentationSDD::AliITSsegmentationSDD(){ - // standard constructor - fGeom=0; - fResponse=0; - fCorr=0; - SetDetSize(); - SetPadSize(); - SetNPads(); - -} -//_____________________________________________________________________________ -AliITSsegmentationSDD& AliITSsegmentationSDD::operator=(AliITSsegmentationSDD &source){ - // Operator = - if(this==&source) return *this; - this->fNsamples = source.fNsamples; - this->fNanodes = source.fNanodes; - this->fPitch = source.fPitch; - this->fTimeStep = source.fTimeStep; - this->fDx = source.fDx; - this->fDz = source.fDz; - this->fDy = source.fDy; - this->fCorr = new TF1(*(source.fCorr)); - this->fGeom = source.fGeom; // Just copy the pointer - this->fResponse = source.fResponse; //Just copy the pointer - return *this; -} -//___________________________________________________________________________ -AliITSsegmentationSDD::AliITSsegmentationSDD(AliITSsegmentationSDD &source){ - // Copy constructor - *this = source; -} -//------------------------------ +//______________________________________________________________________ +AliITSsegmentationSDD::AliITSsegmentationSDD() : AliITSsegmentation(), +fNsamples(0), +fNanodes(0), +fPitch(0), +fTimeStep(0), +fDriftSpeed(0){ + // Default constructor + fDriftSpeed=0; + SetDetSize(fgkDxDefault,fgkDzDefault,fgkDyDefault); + SetPadSize(fgkPitchDefault,fgkClockDefault); + SetNPads(fgkHalfNanodesDefault,fgkNsamplesDefault); + +} + +//______________________________________________________________________ +void AliITSsegmentationSDD::Copy(TObject &obj) const { + // protected method. copy this to obj + AliITSsegmentation::Copy(obj); + ((AliITSsegmentationSDD& ) obj).fNsamples = fNsamples; + ((AliITSsegmentationSDD& ) obj).fNanodes = fNanodes; + ((AliITSsegmentationSDD& ) obj).fPitch = fPitch; + ((AliITSsegmentationSDD& ) obj).fTimeStep = fTimeStep; + ((AliITSsegmentationSDD& ) obj).fDriftSpeed = fDriftSpeed; +} + +//______________________________________________________________________ +AliITSsegmentationSDD& AliITSsegmentationSDD::operator=(const AliITSsegmentationSDD &source){ + // = operator + if(this==&source) return *this; + source.Copy(*this); + return *this; +} + +//____________________________________________________________________________ +AliITSsegmentationSDD::AliITSsegmentationSDD(const AliITSsegmentationSDD &source) : + AliITSsegmentation(source), +fNsamples(0), +fNanodes(0), +fPitch(0), +fTimeStep(0), +fDriftSpeed(0){ + // copy constructor + source.Copy(*this); +} + +//---------------------------------------------------------------------- void AliITSsegmentationSDD::Init(){ // Standard initilisation routine if(!fGeom) { + Fatal("Init","the pointer to the ITS geometry class (AliITSgeom) is null\n"); return; - //fGeom = ((AliITS*)gAlice->GetModule("ITS"))->GetITSgeom(); } AliITSgeomSDD *gsdd = (AliITSgeomSDD *) (fGeom->GetShape(3,1,1)); @@ -83,9 +117,9 @@ void AliITSsegmentationSDD::Init(){ fDy = 2.*kconv*gsdd->GetDy(); } -//------------------------------ +//---------------------------------------------------------------------- void AliITSsegmentationSDD:: -Neighbours(Int_t iX, Int_t iZ, Int_t* Nlist, Int_t Xlist[8], Int_t Zlist[8]){ +Neighbours(Int_t iX, Int_t iZ, Int_t* Nlist, Int_t Xlist[8], Int_t Zlist[8]) const { // returns neighbours for use in Cluster Finder routines and the like if(iX >= fNanodes) printf("iX > fNanodes %d %d\n",iX,fNanodes); @@ -102,79 +136,231 @@ Neighbours(Int_t iX, Int_t iZ, Int_t* Nlist, Int_t Xlist[8], Int_t Zlist[8]){ else Zlist[1]=iZ; Zlist[2]=Zlist[3]=iZ; } -//------------------------------ -void AliITSsegmentationSDD::GetPadIxz(Float_t x,Float_t z,Int_t &timebin,Int_t &anode){ -// Returns cell coordinates (time sample,anode) for given real local coordinates (x,z) +//---------------------------------------------------------------------- +void AliITSsegmentationSDD::GetPadIxz(Float_t x,Float_t z, + Int_t &timebin,Int_t &anode) const { +// Returns cell coordinates (time sample,anode) incremented by 1 !!!!! +// for given real local coordinates (x,z) // expects x, z in cm const Float_t kconv=10000; // cm->um - Float_t speed=fResponse->DriftSpeed(); + x *= kconv; // Convert to microns + z *= kconv; // Convert to microns Int_t na = fNanodes/2; - Float_t driftpath=fDx-TMath::Abs(kconv*x); - timebin=(Int_t)(driftpath/speed/fTimeStep); - anode=(Int_t)(kconv*z/fPitch) + na/2; + Float_t driftpath=fDx-TMath::Abs(x); + timebin=(Int_t)(driftpath/fDriftSpeed/fTimeStep); + anode=(Int_t)(z/fPitch + na/2); if (x > 0) anode += na; timebin+=1; anode+=1; } - -//------------------------------ -void AliITSsegmentationSDD::GetPadCxz(Int_t timebin,Int_t anode,Float_t &x ,Float_t &z){ +//---------------------------------------------------------------------- +void AliITSsegmentationSDD::GetPadCxz(Int_t timebin,Int_t anode, + Float_t &x ,Float_t &z) const{ // Transform from cell to real local coordinates - // returns x, z in cm + // the +0.5 means that an # and time bin # should start from 0 !!! const Float_t kconv=10000; // um->cm + // the +0.5 means that an # and time bin # should start from 0 !!! - Float_t speed=fResponse->DriftSpeed(); Int_t na = fNanodes/2; - Float_t driftpath=(timebin+1)*fTimeStep*speed; + Float_t driftpath=(timebin+0.5)*fTimeStep*fDriftSpeed; if (anode >= na) x=(fDx-driftpath)/kconv; else x = -(fDx-driftpath)/kconv; if (anode >= na) anode-=na; - z=((anode+1)*fPitch-fDz/2)/kconv; + z=((anode+0.5)*fPitch-fDz/2)/kconv; } - -//------------------------------ -void AliITSsegmentationSDD::GetPadTxz(Float_t &x,Float_t &z){ +//---------------------------------------------------------------------- +void AliITSsegmentationSDD::GetPadTxz(Float_t &x,Float_t &z) const{ // Get anode and time bucket as floats - numbering from 0 // expects x, z in cm const Float_t kconv=10000; // cm->um - //Float_t x0=x; - Float_t speed=fResponse->DriftSpeed(); - //Int_t na = fNanodes/2; + Float_t x0=x; + Int_t na = fNanodes/2; Float_t driftpath=fDx-TMath::Abs(kconv*x); - x=driftpath/speed/fTimeStep; - z=kconv*z/fPitch; - // z=kconv*z/fPitch + (float)na/2; - //if (x0 > 0) z += (float)na; + x=driftpath/fDriftSpeed/fTimeStep; + z=kconv*z/fPitch + (float)na/2; + if (x0 < 0) x = -x; } -//------------------------------ -void AliITSsegmentationSDD::GetLocal(Int_t module,Float_t *g ,Float_t *l){ +//---------------------------------------------------------------------- +void AliITSsegmentationSDD::GetLocal(Int_t module,Float_t *g ,Float_t *l) const { // returns local coordinates from global if(!fGeom) { + Fatal("GetLocal","the pointer to the ITS geometry class (AliITSgeom) is null\n"); return; - //fGeom = ((AliITS*)gAlice->GetModule("ITS"))->GetITSgeom(); } fGeom->GtoL(module,g,l); } -//------------------------------ -void AliITSsegmentationSDD::GetGlobal(Int_t module,Float_t *l ,Float_t *g){ +//---------------------------------------------------------------------- +void AliITSsegmentationSDD::GetGlobal(Int_t module,Float_t *l ,Float_t *g) const { // return global coordinates from local if(!fGeom) { - return; - //fGeom = ((AliITS*)gAlice->GetModule("ITS"))->GetITSgeom(); + Fatal("GetGlobal","the pointer to the ITS geometry class (AliITSgeom) is null\n"); } fGeom->LtoG(module,l,g); } +//---------------------------------------------------------------------- +void AliITSsegmentationSDD::Print(Option_t *opt) const { + // Print SDD segmentation Parameters + + cout << "**************************************************" << endl; + cout << " Silicon Drift Detector Segmentation Parameters " << opt << endl; + cout << "**************************************************" << endl; + cout << "Number of Time Samples: " << fNsamples << endl; + cout << "Number of Anodes: " << fNanodes << endl; + cout << "Time Step (ns): " << fTimeStep << endl; + cout << "Anode Pitch (um): " << fPitch << endl; + cout << "Full Detector Width (x): " << fDx; + cout<<" (Default is "<DeriftSpeed()=7.3mic/ns, Dpz()=512. For other values a only the +// specific numbers will change not their layout.} +// +// 0 191 0 +// 0 |----------------------|---------------------| 256 +// | a time-bins | time-bins a | +// | n | n | +// | o |___________________o_|__> X +// | d | d | +// | e | e | +// | s | s | +// 255 |----------------------|---------------------| 511 +// | +// V +// Z + Float_t dx,dz,tb; + const Float_t kconv = 1.0E-04; // converts microns to cm. + + ix = -1; // default values + iz = -1; // default values + dx = -kconv*Dx(); // lower left edge in cm. + dz = -0.5*kconv*Dz(); // lower left edge in cm. + if(x-dx) { + Warning("LocalToDet","input argument %f out of range (%f, %f)",x,dx,-dx); + return kFALSE; // outside of defined volume. + } + if(z-dz) { + Warning("LocalToDet","input argument %f out of range (%f, %f)",z,dz,-dz); + return kFALSE; // outside of defined volume. + } + tb = fDriftSpeed*fTimeStep*kconv; // compute size of time bin. + if(x>0) dx = -(dx + x)/tb; // distance from + side in time bin units + else dx = (x - dx)/tb; // distance from - side in time bin units + dz = (z - dz)/(kconv*fPitch); // distance in z in anode pitch units + ix = (Int_t) dx; // time bin + iz = (Int_t) dz; // anode + if(x>0) iz += Npz()/2; // if x>0 then + side anodes values. + return kTRUE; // Found ix and iz, return. +} +//______________________________________________________________________ +void AliITSsegmentationSDD::DetToLocal(Int_t ix,Int_t iz,Float_t &x,Float_t &z) const +{ +// Transformation from Detector time bucket and anode coordiantes to Geant +// detector centerd local coordinates (cm). +// Input: +// Int_t ix detector x time coordinate. Has the range 0<=ixDeriftSpeed()=7.3mic/ns, Dpz()=512. For other values a only the +// specific numbers will change not their layout.} +// +// 0 191 0 +// 0 |----------------------|---------------------| 256 +// | a time-bins | time-bins a | +// | n | n | +// | o |___________________o_|__> X +// | d | d | +// | e | e | +// | s | s | +// 255 |----------------------|---------------------| 511 +// | +// V +// Z + Int_t i,j; + Float_t tb; + const Float_t kconv = 1.0E-04; // converts microns to cm. + + if(iz>=Npz()/2) x = kconv*Dx(); // default value for +x side. + else x = -kconv*Dx(); // default value for -x side. + z = -0.5*kconv*Dz(); // default value. + if(ix<0 || ix>=Npx()) { + Warning("DetToLocal","input argument %d out of range (0, %d)",ix,Npx()); + return; // outside of detector + } + if(iz<0 || iz>=Npz()) { + Warning("DetToLocal","input argument %d out of range (0, %d)",iz,Npz()); + return; // outside of detctor + } + tb = fDriftSpeed*fTimeStep*kconv; // compute size of time bin. + if(iz>=Npz()/2) tb *= -1.0; // for +x side decrement frmo Dx(). + for(i=0;i=Npz()/2) iz -=Npz()/2;// If +x side don't count anodes from -x side. + for(j=0;j