ITS/AliITSPreprocessorSDD.cxx

   1 /////////////////////////////////////////
   2 // Class for SDD digits preprocessing  //
   3 //                                     //
   4 //                                     //
   5 ////////////////////////////////////////
   6
   7 #include "AliITSPreprocessorSDD.h"
   8 #include "AliITSCalibrationSDD.h"
   9 #include "AliShuttleInterface.h"
  10 #include "AliCDBMetaData.h"
  11 #include "TObjArray.h"
  12 #include "TH1F.h"
  13 #include "AliLog.h"
  14 #include "TFile.h"
  15
  16 const Int_t AliITSPreprocessorSDD::fgkNumberOfSDD = 260;
  17 const Int_t AliITSPreprocessorSDD::fgkNumberOfChannels = 512;
  18 const TString AliITSPreprocessorSDD::fgkNameHistoPedestals = "hpedestal";
  19 const TString AliITSPreprocessorSDD::fgkNameHistoNoise = "hnoise";
  20 ClassImp(AliITSPreprocessorSDD)
  21
  22
  23 UInt_t AliITSPreprocessorSDD::Process(TMap*/* dcsAliasMap*/){
  24
  25   //preprocessing.
  26
  27   UInt_t result = 0;
  28   const char* filename = GetFile(kDAQ,"PEDESTALS","GDC");
  29   const char* filenamen= GetFile(kDAQ,"NOISE","GDC");
  30   const char* filenamed= GetFile(kDAQ,"DEADCHANNELS","GDC");
  31   TFile* f1 = TFile::Open(filename,"r");
  32   TFile* f2 = TFile::Open(filenamen,"r");
  33   Char_t namehisto[20];
  34   Char_t namehisto2[20];
  35
  36   FILE* filed = fopen(filenamed,"read");
  37   Int_t numOfBadChannels[fgkNumberOfSDD];
  38   Int_t** badCh = new Int_t*[fgkNumberOfSDD];
  39
  40   Char_t row[50];
  41   Int_t nSDD=0;
  42   Char_t str[20];
  43   char dims[1];
  44   Int_t dim;
  45   sprintf(str,"MODULE=%d",0);
  46   while(!feof(filed)){
  47     fscanf(filed,"%s\n",row);
  48     if(strcmp(row,str)==0){
  49       fscanf(filed,"%s %d\n",dims,&dim);
  50       badCh[nSDD] = new Int_t[dim];
  51       numOfBadChannels[nSDD]=dim;
  52       for(Int_t ibad=0;ibad<dim;ibad++){
  53         fscanf(filed,"%d\n",&badCh[nSDD][ibad]);
  54       }
  55     }
  56     nSDD++;
  57     sprintf(str,"MODULE=%d",nSDD);
  58   }
  59
  60
  61
  62   AliCDBMetaData *md1= new AliCDBMetaData(); // metaData describing the object
  63   md1->SetObjectClassName("AliITSCalibration");
  64   md1->SetResponsible("Elisabetta Crescio");
  65   md1->SetBeamPeriod(0);
  66   md1->SetAliRootVersion("head September 2005"); //root version
  67   md1->SetComment("This is a test");
  68
  69   TObjArray respSDD(260);
  70   respSDD.SetOwner(kFALSE);
  71
  72   for(Int_t imod=0;imod<fgkNumberOfSDD;imod++){
  73     AliITSCalibrationSDD *cal = new AliITSCalibrationSDD("simulated");
  74     cal->SetDeadChannels(numOfBadChannels[imod]);
  75     for(Int_t ich=0;ich<numOfBadChannels[imod];ich++){
  76       cal->SetBadChannel(ich,badCh[imod][ich]);
  77     }
  78     sprintf(namehisto,"%s_%d",fgkNameHistoPedestals.Data(),imod);
  79     sprintf(namehisto2,"%s_%d",fgkNameHistoNoise.Data(),imod);
  80     TH1F* hbas = (TH1F*)f1->Get(namehisto);
  81     TH1F* hnoi = (TH1F*)f2->Get(namehisto2);
  82     for(Int_t ien=0;ien<fgkNumberOfChannels;ien++){
  83       cal->SetBaseline(ien,hbas->GetBinContent(ien+1));
  84       cal->SetNoiseAfterElectronics(ien,hnoi->GetBinContent(ien+1));
  85     }
  86     respSDD.Add(cal);
  87   }
  88
  89   result = Store("Calib","Data",&respSDD,md1);
  90
  91   for(Int_t i=0;i<fgkNumberOfSDD;i++){
  92     delete badCh[i];
  93   }
  94   delete [] badCh;
  95   f1->Close();
  96   f2->Close();
  97   fclose(filed);
  98   return result;
  99
 100 }