ITS/AliITSresponseSDD.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   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  **************************************************************************/
  15
  16 /* $Id$ */
  17
  18 //////////////////////////////////////////////////////
  19 //  Base response class forITS                      //
  20 //  It is used to set static data members           //
  21 //  connected to parameters equal for all           //
  22 //  the modules                                     //
  23 //                                                  //
  24 //                                                  //
  25 //////////////////////////////////////////////////////
  26
  27 #include <TMath.h>
  28
  29 #include "AliITSresponseSDD.h"
  30 #include <AliITSgeomTGeo.h>
  31
  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)
  37
  38 //_________________________________________________________________________
  39 AliITSresponseSDD::AliITSresponseSDD():
  40 TObject(),
  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;
  49   }
  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)
  52 }
  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)
  56   Int_t minMod,maxMod;
  57   if(lay==3){
  58     minMod=1;
  59     maxMod=3;
  60     if(lad>kNLaddersLay3){
  61       AliError(Form("Ladder number %d out of range",lad));
  62       return;
  63     }
  64   }else if(lay==4){
  65     minMod=1;
  66     maxMod=4;
  67     if(lad>kNLaddersLay4){
  68       AliError(Form("Ladder number %d out of range",lad));
  69       return;
  70     }
  71   }else{
  72     AliError(Form("Layer number %d out of range",lay));
  73     return;
  74   }
  75   for(Int_t iMod=minMod; iMod<=maxMod; iMod++){
  76     Int_t modIndex=AliITSgeomTGeo::GetModuleIndex(lay,lad,iMod);
  77     SetModuleTimeZero(modIndex,tzero);
  78   }
  79 }
  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)
  83   Int_t minMod,maxMod;
  84   if(lay==3){
  85     minMod=4;
  86     maxMod=6;
  87     if(lad>kNLaddersLay3){
  88       AliError(Form("Ladder number %d out of range",lad));
  89       return;
  90     }
  91   }else if(lay==4){
  92     minMod=5;
  93     maxMod=8;
  94     if(lad>kNLaddersLay4){
  95       AliError(Form("Ladder number %d out of range",lad));
  96       return;
  97     }
  98   }else{
  99     AliError(Form("Layer number %d out of range",lay));
 100     return;
 101   }
 102   for(Int_t iMod=minMod; iMod<=maxMod; iMod++){
 103     Int_t modIndex=AliITSgeomTGeo::GetModuleIndex(lay,lad,iMod);
 104     SetModuleTimeZero(modIndex,tzero);
 105   }
 106 }
 107 //_________________________________________________________________________
 108 void AliITSresponseSDD::PrintChargeCalibrationParams() const{
 109   //
 110   printf("ADC vs. drift time corr=%f\n",GetChargevsTime());
 111   printf("-------------------------------------\n");
 112   printf("Layer 3\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);
 117       printf("%7.2f   ",tz);
 118     }
 119     printf("\n");
 120   }
 121   printf("\n");
 122   printf("Layer 4\n");
 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);
 127       printf("%7.2f   ",tz);
 128     }
 129     printf("\n");
 130   }
 131 }
 132 //_________________________________________________________________________
 133 void AliITSresponseSDD::PrintTimeZeroes() const{
 134   //
 135   printf("Layer 3\n");
 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);
 140       printf("%7.2f   ",tz);
 141     }
 142     printf("\n");
 143   }
 144   printf("\n");
 145   printf("Layer 4\n");
 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);
 150       printf("%7.2f   ",tz);
 151     }
 152     printf("\n");
 153   }
 154
 155 }
 156 //_________________________________________________________________________
 157 void AliITSresponseSDD::PrintVdriftCorerctions() const{
 158   //
 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));
 161   }
 162 }