1 /**************************************************************************
2 * Copyright(c) 2007-2009, 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 <Riostream.h>
20 #include "AliITSCalibrationSDD.h"
23 //////////////////////////////////////////////////////
24 // Calibration class for set:ITS //
25 // Specific subdetector implementation //
26 // for silicon drift detectors //
29 //////////////////////////////////////////////////////
31 const Float_t AliITSCalibrationSDD::fgkTemperatureDefault = 296.;
32 const Float_t AliITSCalibrationSDD::fgkNoiseDefault = 2.38;
33 const Float_t AliITSCalibrationSDD::fgkGainDefault = 1.;
34 const Float_t AliITSCalibrationSDD::fgkBaselineDefault = 20.;
35 //______________________________________________________________________
39 ClassImp(AliITSCalibrationSDD)
41 AliITSCalibrationSDD::AliITSCalibrationSDD():
56 // default constructor
59 for(Int_t ian=0;ian<fgkWings*fgkChannels*fgkChips;ian++){
60 fBaseline[ian]=fgkBaselineDefault;
61 fNoise[ian]=fgkNoiseDefault;
63 SetNoiseAfterElectronics(ian);
65 for(Int_t iw=0;iw<fgkWings;iw++){
66 SetZSLowThreshold(iw);
67 SetZSHighThreshold(iw);
68 for(Int_t icp=0;icp<fgkChips;icp++){
69 Int_t chipindex=iw*fgkChips+icp;
70 fIsChipBad[chipindex]=kFALSE;
73 SetTemperature(fgkTemperatureDefault);
76 //______________________________________________________________________
77 AliITSCalibrationSDD::AliITSCalibrationSDD(const char *dataType):
95 for(Int_t ian=0;ian<fgkWings*fgkChannels*fgkChips;ian++){
96 fBaseline[ian]=fgkBaselineDefault;
97 fNoise[ian]=fgkNoiseDefault;
99 SetNoiseAfterElectronics(ian);
101 for(Int_t iw=0;iw<fgkWings;iw++){
102 SetZSLowThreshold(iw);
103 SetZSHighThreshold(iw);
104 for(Int_t icp=0;icp<fgkChips;icp++){
105 Int_t chipindex=iw*fgkChips+icp;
106 fIsChipBad[chipindex]=kFALSE;
110 SetTemperature(fgkTemperatureDefault);
111 SetDataType(dataType);
113 //_____________________________________________________________________
114 AliITSCalibrationSDD::~AliITSCalibrationSDD(){
117 if(fMapAW0) delete fMapAW0;
118 if(fMapAW1) delete fMapAW1;
119 if(fMapTW0) delete fMapTW0;
120 if(fMapTW1) delete fMapTW1;
121 if(fDrSpeed0) delete fDrSpeed0;
122 if(fDrSpeed1) delete fDrSpeed1;
125 //______________________________________________________________________
126 void AliITSCalibrationSDD::GiveCompressParam(Int_t cp[4]) const {
127 // give compression param
133 //_____________________________________________________________________
134 void AliITSCalibrationSDD::SetBadChannel(Int_t i,Int_t anode){
135 //Set bad anode (set gain=0 for these channels);
137 if(anode<0 || anode >fgkChannels*fgkChips*fgkWings-1){
138 AliError("Wrong anode number");
141 fBadChannels[i]=anode;
144 //______________________________________________________________________
145 void AliITSCalibrationSDD::GetCorrections(Float_t z, Float_t x, Float_t &devz, Float_t &devx, AliITSsegmentationSDD* seg){
146 //correction of coordinates using the maps stored in the DB
147 Int_t nSide=seg->GetSideFromLocalX(x);
149 // if(nSide==0) devz=fMapAW0->GetCorrection(z,x,seg);
150 // else devz=fMapAW1->GetCorrection(z,x,seg);
152 if(nSide==0) devx=fMapTW0->GetCorrection(z,x,seg);
153 else devx=fMapTW1->GetCorrection(z,x,seg);
156 //______________________________________________________________________
157 void AliITSCalibrationSDD::GetShiftsForSimulation(Float_t z, Float_t x, Float_t &devz, Float_t &devx, AliITSsegmentationSDD* seg){
158 //correction of coordinates using the maps stored in the DB
159 Int_t nSide=seg->GetSideFromLocalX(x);
161 // if(nSide==0) devz=fMapAW0->GetCorrection(z,x,seg);
162 // else devz=fMapAW1->GetCorrection(z,x,seg);
164 if(nSide==0) devx=fMapTW0->GetShiftForSimulation(z,x,seg);
165 else devx=fMapTW1->GetShiftForSimulation(z,x,seg);
168 //______________________________________________________________________
169 void AliITSCalibrationSDD::PrintGains() const{
172 if( GetDeadChips() == 0 &&
173 GetDeadChannels() == 0 )
176 // Print Electronics Gains
177 cout << "**************************************************" << endl;
178 cout << " Print Electronics Gains " << endl;
179 cout << "**************************************************" << endl;
181 // Print SDD electronic gains
182 for(Int_t ian=0; ian<fgkWings*fgkChips*fgkChannels;ian++){
183 printf("Gain for channel %d = %f\n",ian,fGain[ian]);
187 //______________________________________________________________________
188 void AliITSCalibrationSDD::Print(){
189 // Print SDD response Parameters
191 cout << "**************************************************" << endl;
192 cout << " Silicon Drift Detector Response Parameters " << endl;
193 cout << "**************************************************" << endl;
194 cout << "Hardware compression parameters: " << endl;
195 cout << "Noise before electronics (arbitrary units): " << fNoise[0] << endl;
196 cout << "Baseline (ADC units): " << fBaseline[0] << endl;
197 cout << "Noise after electronics (ADC units): " << fNoiseAfterEl[0] << endl;
198 cout << "Temperature: " << Temperature() << " K " << endl;