PHOS/AliPHOSCalibData.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 //                                                                           //
  20 // class for PHOS calibration                                                 //
  21 //                                                                           //
  22 ///////////////////////////////////////////////////////////////////////////////
  23
  24 #include "TMath.h"
  25 #include "TRandom.h"
  26 #include "AliPHOSCalibData.h"
  27 #include "AliCDBManager.h"
  28 #include "AliCDBStorage.h"
  29 #include "AliCDBId.h"
  30 #include "AliCDBEntry.h"
  31
  32 ClassImp(AliPHOSCalibData)
  33
  34 //________________________________________________________________
  35   AliPHOSCalibData::AliPHOSCalibData():
  36     TNamed(), fCalibDataEmc(0x0), fCalibDataCpv(0x0)
  37 {
  38   // Default constructor
  39
  40   fEmcDataPath="PHOS/Calib/EmcGainPedestals";
  41   fCpvDataPath="PHOS/Calib/CpvGainPedestals";
  42
  43   AliCDBEntry* entryEmc = AliCDBManager::Instance()->Get(fEmcDataPath.Data());
  44   if(entryEmc)
  45     fCalibDataEmc = (AliPHOSEmcCalibData*)entryEmc->GetObject();
  46
  47   AliCDBEntry* entryCpv = AliCDBManager::Instance()->Get(fCpvDataPath.Data());
  48   if(entryCpv)
  49     fCalibDataCpv = (AliPHOSCpvCalibData*)entryCpv->GetObject();
  50
  51 }
  52
  53 //________________________________________________________________
  54 AliPHOSCalibData::AliPHOSCalibData(Int_t runNumber) :
  55   TNamed("phosCalib","PHOS Calibration Data Manager"),
  56   fCalibDataEmc(0x0), fCalibDataCpv(0x0)
  57 {
  58   // Constructor
  59
  60   fEmcDataPath="PHOS/Calib/EmcGainPedestals";
  61   fCpvDataPath="PHOS/Calib/CpvGainPedestals";
  62
  63   AliCDBEntry* entryEmc = AliCDBManager::Instance()->Get(fEmcDataPath.Data(),runNumber);
  64   if(entryEmc)
  65     fCalibDataEmc = (AliPHOSEmcCalibData*)entryEmc->GetObject();
  66
  67   AliCDBEntry* entryCpv = AliCDBManager::Instance()->Get(fCpvDataPath.Data(),runNumber);
  68   if(entryCpv)
  69     fCalibDataCpv = (AliPHOSCpvCalibData*)entryCpv->GetObject();
  70
  71 }
  72
  73 //________________________________________________________________
  74 AliPHOSCalibData::~AliPHOSCalibData()
  75 {
  76   // Destructor
  77
  78 }
  79
  80 //________________________________________________________________
  81 void AliPHOSCalibData::Reset()
  82 {
  83   // Set all pedestals to 0 and all ADC channels to 1
  84
  85   fCalibDataEmc->Reset();
  86   fCalibDataCpv->Reset();
  87 }
  88
  89 //________________________________________________________________
  90 void  AliPHOSCalibData::Print(Option_t *option) const
  91 {
  92   if (fCalibDataEmc) fCalibDataEmc->Print(option);
  93   if (fCalibDataCpv) fCalibDataCpv->Print(option);
  94 }
  95
  96 //________________________________________________________________
  97 void AliPHOSCalibData::CreateNew()
  98 {
  99   if(fCalibDataEmc) delete fCalibDataEmc;
 100   fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
 101
 102   if(fCalibDataCpv) delete fCalibDataCpv;
 103   fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
 104
 105 }
 106
 107 //________________________________________________________________
 108 Bool_t AliPHOSCalibData::WriteEmc(Int_t firstRun, Int_t lastRun, AliCDBMetaData *md)
 109 {
 110
 111   if(!fCalibDataEmc) return kFALSE;
 112
 113   AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS");
 114   if(!storage)
 115     storage = AliCDBManager::Instance()->GetDefaultStorage();
 116
 117   if(storage) {
 118     AliCDBId id(fEmcDataPath.Data(),firstRun,lastRun);
 119     storage->Put(fCalibDataEmc,id, md);
 120     return kTRUE;
 121   }
 122   else
 123     return kFALSE;
 124
 125 }
 126
 127 //________________________________________________________________
 128 Bool_t AliPHOSCalibData::WriteCpv(Int_t firstRun, Int_t lastRun, AliCDBMetaData *md)
 129 {
 130
 131   if(!fCalibDataCpv) return kFALSE;
 132
 133   AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS");
 134   if(!storage)
 135     storage = AliCDBManager::Instance()->GetDefaultStorage();
 136
 137   if(storage) {
 138     AliCDBId id(fCpvDataPath.Data(),firstRun,lastRun);
 139     storage->Put(fCalibDataCpv,id, md);
 140     return kTRUE;
 141   }
 142   else
 143     return kFALSE;
 144
 145 }
 146
 147 //________________________________________________________________
 148 Float_t AliPHOSCalibData::GetADCchannelEmc(Int_t module, Int_t column, Int_t row) const
 149 {
 150   //module, column,raw should follow the internal PHOS convention:
 151   //module 1:5, column 1:56, row 1:64
 152
 153   if(fCalibDataEmc)
 154     return fCalibDataEmc->GetADCchannelEmc(module,column,row);
 155   else
 156     return 0.0015; // default width of one EMC ADC channel in GeV
 157 }
 158
 159 Float_t AliPHOSCalibData::GetADCpedestalEmc(Int_t module, Int_t column, Int_t row) const
 160 {
 161   if(fCalibDataEmc)
 162     return fCalibDataEmc->GetADCpedestalEmc(module,column,row);
 163   else
 164     return 0.005; // default EMC ADC pedestal
 165 }
 166
 167 void AliPHOSCalibData::SetADCchannelEmc(Int_t module, Int_t column, Int_t row, Float_t value)
 168 {
 169   if(!fCalibDataEmc)
 170     fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
 171
 172   fCalibDataEmc->SetADCchannelEmc(module,column,row,value);
 173 }
 174
 175 void AliPHOSCalibData::SetADCpedestalEmc(Int_t module, Int_t column, Int_t row, Float_t value)
 176 {
 177   if(!fCalibDataEmc)
 178     fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
 179
 180   fCalibDataEmc->SetADCpedestalEmc(module,column,row,value);
 181 }
 182
 183 //________________________________________________________________
 184 Float_t AliPHOSCalibData::GetADCchannelCpv(Int_t module, Int_t column, Int_t row) const
 185 {
 186   //module, column,raw should follow the internal CPV convention:
 187   //module 1:5, column 1:64, row 1:128
 188
 189   if(fCalibDataCpv)
 190     return fCalibDataCpv->GetADCchannelCpv(module,column,row);
 191   else
 192     return 0.0012; // default width of one ADC channel in CPV 'popugais'
 193 }
 194
 195 Float_t AliPHOSCalibData::GetADCpedestalCpv(Int_t module, Int_t column, Int_t row) const
 196 {
 197   if(fCalibDataCpv)
 198     return fCalibDataCpv->GetADCpedestalCpv(module,column,row);
 199   else
 200     return 0.012; // default CPV ADC pedestal
 201 }
 202
 203 void AliPHOSCalibData::SetADCchannelCpv(Int_t module, Int_t column, Int_t row, Float_t value)
 204 {
 205   if(!fCalibDataCpv)
 206     fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
 207
 208   fCalibDataCpv->SetADCchannelCpv(module,column,row,value);
 209 }
 210
 211 void AliPHOSCalibData::SetADCpedestalCpv(Int_t module, Int_t column, Int_t row, Float_t value)
 212 {
 213   if(!fCalibDataCpv)
 214     fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
 215
 216   fCalibDataCpv->SetADCpedestalCpv(module,column,row,value);
 217 }
 218
 219 //________________________________________________________________
 220 void AliPHOSCalibData::RandomEmc()
 221 {
 222
 223   if(fCalibDataEmc) delete fCalibDataEmc;
 224   fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
 225
 226   TRandom rn;
 227   rn.SetSeed(0); //the seed is set to the current  machine clock
 228
 229   Float_t ADCchanelEmc,ADCpedestalEmc;
 230
 231   for(Int_t module=1; module<6; module++) {
 232     for(Int_t column=1; column<57; column++) {
 233       for(Int_t row=1; row<65; row++) {
 234         ADCchanelEmc=rn.Uniform(0.00075,0.00375); // Cmax/Cmin = 5, (Cmax-Cmin)/2 = 0.0015
 235         ADCpedestalEmc=rn.Uniform(0.0045,0.0055); //+-10% spread of pedestals from 0.005
 236         fCalibDataEmc->SetADCchannelEmc(module,column,row,ADCchanelEmc);
 237         fCalibDataEmc->SetADCpedestalEmc(module,column,row,ADCpedestalEmc);
 238       }
 239     }
 240   }
 241
 242 }
 243
 244 //________________________________________________________________
 245 void AliPHOSCalibData::RandomCpv()
 246 {
 247
 248   if(fCalibDataCpv) delete fCalibDataCpv;
 249   fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
 250
 251   TRandom rn;
 252   rn.SetSeed(0); //the seed is set to the current  machine clock
 253
 254   Float_t ADCchanelCpv,ADCpedestalCpv;
 255
 256   for(Int_t module=1; module<6; module++) {
 257     for(Int_t column=1; column<65; column++) {
 258       for(Int_t row=1; row<129; row++) {
 259         ADCchanelCpv=TMath::Abs(rn.Uniform(0.0009,0.0015)); // 0.0012 +- 25%
 260         ADCpedestalCpv=rn.Uniform(0.0048,0.0192); // Ped[max]/Ped[min] = 4, <Ped> = 0.012
 261         fCalibDataCpv->SetADCchannelCpv(module,column,row,ADCchanelCpv);
 262         fCalibDataCpv->SetADCpedestalCpv(module,column,row,ADCpedestalCpv);
 263       }
 264     }
 265   }
 266 }