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 //////////////////////////////////////////////////////
19 // Base response class forITS //
20 // It is used to set static data members //
21 // connected to parameters equal for all //
25 //////////////////////////////////////////////////////
29 #include "AliITSresponseSDD.h"
30 #include <AliITSgeomTGeo.h>
32 const Float_t AliITSresponseSDD::fgkTimeOffsetDefault = 54.30;
33 const Float_t AliITSresponseSDD::fgkADC2keVDefault = 3.34;
34 const Float_t AliITSresponseSDD::fgkChargevsTimeDefault = 0.00355;
35 const Float_t AliITSresponseSDD::fgkCarlosRXClockPeriod = 25.;
36 ClassImp(AliITSresponseSDD)
38 //_________________________________________________________________________
39 AliITSresponseSDD::AliITSresponseSDD():
41 fTimeOffset(fgkTimeOffsetDefault),
42 fADC2keV(fgkADC2keVDefault),
43 fChargevsTime(fgkChargevsTimeDefault){
44 // default constructor
45 for(Int_t i=0; i<kNSDDmods;i++){
46 fTimeZero[i]=fgkTimeOffsetDefault;
47 fDeltaVDrift[i] = fDeltaVDrift[i+kNSDDmods] = 0.;
48 fADCtokeV[i]=fgkADC2keVDefault;
50 SetVDCorr2Side(kTRUE); // default for new objects will be separate corrections for 2 sides (bwd compatible)
51 // SetVDCorrMult(kTRUE); // default for new objects will have multiplicative correction v'=(1+corr)*v (bwd compatible)
53 //_________________________________________________________________________
54 void AliITSresponseSDD::SetHalfLadderATimeZero(Int_t lay, Int_t lad, Float_t tzero){
55 // Sets time Zero for all modules of a ladder on side A (Z>0)
60 if(lad>kNLaddersLay3){
61 AliError(Form("Ladder number %d out of range",lad));
67 if(lad>kNLaddersLay4){
68 AliError(Form("Ladder number %d out of range",lad));
72 AliError(Form("Layer number %d out of range",lay));
75 for(Int_t iMod=minMod; iMod<=maxMod; iMod++){
76 Int_t modIndex=AliITSgeomTGeo::GetModuleIndex(lay,lad,iMod);
77 SetModuleTimeZero(modIndex,tzero);
80 //_________________________________________________________________________
81 void AliITSresponseSDD::SetHalfLadderCTimeZero(Int_t lay, Int_t lad, Float_t tzero){
82 // Sets time Zero for all modules of a ladder on side C (Z<0)
87 if(lad>kNLaddersLay3){
88 AliError(Form("Ladder number %d out of range",lad));
94 if(lad>kNLaddersLay4){
95 AliError(Form("Ladder number %d out of range",lad));
99 AliError(Form("Layer number %d out of range",lay));
102 for(Int_t iMod=minMod; iMod<=maxMod; iMod++){
103 Int_t modIndex=AliITSgeomTGeo::GetModuleIndex(lay,lad,iMod);
104 SetModuleTimeZero(modIndex,tzero);
107 //_________________________________________________________________________
108 void AliITSresponseSDD::PrintChargeCalibrationParams() const{
110 printf("ADC vs. drift time corr=%f\n",GetChargevsTime());
111 printf("-------------------------------------\n");
113 for(Int_t ilad=1; ilad<=14; ilad++){
114 for(Int_t idet=1; idet<=6;idet++){
115 Int_t modIndex=AliITSgeomTGeo::GetModuleIndex(3,ilad,idet);
116 Float_t tz=GetADCtokeV(modIndex);
123 for(Int_t ilad=1; ilad<=22; ilad++){
124 for(Int_t idet=1; idet<=8;idet++){
125 Int_t modIndex=AliITSgeomTGeo::GetModuleIndex(4,ilad,idet);
126 Float_t tz=GetADCtokeV(modIndex);
132 //_________________________________________________________________________
133 void AliITSresponseSDD::PrintTimeZeroes() const{
136 for(Int_t ilad=1; ilad<=14; ilad++){
137 for(Int_t idet=1; idet<=6;idet++){
138 Int_t modIndex=AliITSgeomTGeo::GetModuleIndex(3,ilad,idet);
139 Float_t tz=GetTimeZero(modIndex);
146 for(Int_t ilad=1; ilad<=22; ilad++){
147 for(Int_t idet=1; idet<=8;idet++){
148 Int_t modIndex=AliITSgeomTGeo::GetModuleIndex(4,ilad,idet);
149 Float_t tz=GetTimeZero(modIndex);
156 //_________________________________________________________________________
157 void AliITSresponseSDD::PrintVdriftCorerctions() const{
159 for(Int_t iMod=240; iMod<500; iMod++){
160 printf("Module %d dVleft=%f dVright=%f\n",iMod,GetDeltaVDrift(iMod,0),GetDeltaVDrift(iMod,1));