ITS/AliITSOnlineSDDCMN.cxx

   1 /**************************************************************************
   2  * Copyright(c) 2007-2009, 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 #include "AliITSOnlineSDDCMN.h"
  16 #include "AliLog.h"
  17 #include <TH2F.h>
  18 #include <TMath.h>
  19
  20
  21 ///////////////////////////////////////////////////////////////////
  22 //                                                               //
  23 // Implementation of the class used for analysis of SDD noise    //
  24 // corrected for common mode                                     //
  25 // Origin: F.Prino, Torino, prino@to.infn.it                     //
  26 //                                                               //
  27 ///////////////////////////////////////////////////////////////////
  28
  29
  30 ClassImp(AliITSOnlineSDDCMN)
  31 //______________________________________________________________________
  32   AliITSOnlineSDDCMN::AliITSOnlineSDDCMN():AliITSOnlineSDD(),fNEvents(0),fMinCorrNoise(0.),fMaxCorrNoise(0.),fNSigmaNoise(0.)
  33 {
  34   // default constructor
  35   Reset();
  36   SetMinNoise();
  37   SetMaxNoise();
  38   SetNSigmaNoise();
  39 }
  40 //______________________________________________________________________
  41   AliITSOnlineSDDCMN::AliITSOnlineSDDCMN(Int_t mod, Int_t sid):AliITSOnlineSDD(mod,sid),fNEvents(0),fMinCorrNoise(0.),fMaxCorrNoise(0.),fNSigmaNoise(0.)
  42 {
  43   // default constructor
  44   Reset();
  45   SetMinNoise();
  46   SetMaxNoise();
  47   SetNSigmaNoise();
  48 }
  49 //______________________________________________________________________
  50 AliITSOnlineSDDCMN::~AliITSOnlineSDDCMN(){
  51   // Destructor
  52 }
  53 //______________________________________________________________________
  54 void AliITSOnlineSDDCMN::Reset(){
  55   fNEvents=0;
  56   for(Int_t i=0;i<fgkNAnodes;i++){
  57     fGoodAnode[i]=1;
  58     fBaseline[i]=0.;
  59     fCMN[i]=0.;
  60     fSumCorrNoise[i]=0.;
  61   }
  62   ReadBaselines();
  63 }
  64 //______________________________________________________________________
  65 void AliITSOnlineSDDCMN::ReadBaselines(){
  66   // assume baselines and good anodes are taken from previous run
  67   Char_t basfilnam[100];
  68   sprintf(basfilnam,"SDDbase_step1_mod%03d_sid%d.data",fModuleId,fSide);
  69   FILE* basf=fopen(basfilnam,"r");
  70   if(basf==0){
  71     AliWarning("Baselinefile not present, Set all baselines to 50\n");
  72     for(Int_t ian=0;ian<fgkNAnodes;ian++){
  73       fBaseline[ian]=50.;
  74       fGoodAnode[ian]=1;
  75     }
  76     return;
  77   }
  78   Int_t n,ok;
  79   Float_t base,rms,cmn,corrnoi;
  80   for(Int_t ian=0;ian<fgkNAnodes;ian++){
  81     fscanf(basf,"%d %d %f %f %f %f\n",&n,&ok,&base,&rms,&cmn,&corrnoi);
  82     fGoodAnode[ian]=ok;
  83     fBaseline[ian]=base;
  84     fRawNoise[ian]=rms;
  85     fCMN[ian]=cmn;
  86   }
  87   fclose(basf);
  88 }
  89 //______________________________________________________________________
  90 void  AliITSOnlineSDDCMN::ValidateAnodes(){
  91   for(Int_t ian=0;ian<fgkNAnodes;ian++){
  92     if(!fGoodAnode[ian]) continue;
  93     if(GetAnodeCorrNoise(ian)>fMaxCorrNoise || GetAnodeCorrNoise(ian)<fMinCorrNoise) fGoodAnode[ian]=0;
  94     if(GetAnodeCorrNoise(ian)>fNSigmaNoise*CalcMeanNoise()) fGoodAnode[ian]=0;
  95   }
  96 }
  97
  98 //______________________________________________________________________
  99 void AliITSOnlineSDDCMN::AddEvent(TH2F* hrawd){
 100   //
 101   fNEvents++;
 102   Float_t tbmax=(Float_t)hrawd->GetNbinsX();
 103   TH2F* hcorrd=new TH2F("hcorrd","",hrawd->GetNbinsX(),hrawd->GetXaxis()->GetXmin(),hrawd->GetXaxis()->GetXmax(),hrawd->GetNbinsY(),hrawd->GetYaxis()->GetXmin(),hrawd->GetYaxis()->GetXmax());
 104   for(Int_t itb=0;itb<tbmax;itb++){
 105     Float_t sumEven=0., sumOdd=0.;
 106     Int_t countEven=0, countOdd=0;
 107     for(Int_t ian=0;ian<fgkNAnodes;ian+=2){
 108       if(!fGoodAnode[ian]) continue;
 109       sumEven+=hrawd->GetBinContent(itb+1,ian+1)-fBaseline[ian];
 110       countEven++;
 111     }
 112     for(Int_t ian=1;ian<fgkNAnodes;ian+=2){
 113       if(!fGoodAnode[ian]) continue;
 114       sumOdd+=hrawd->GetBinContent(itb+1,ian+1)-fBaseline[ian];
 115       countOdd++;
 116     }
 117     for(Int_t ian=0;ian<fgkNAnodes;ian++){
 118       if(!fGoodAnode[ian]) continue;
 119       Float_t meanN;
 120       if(ian%2==0) meanN=sumEven/(Float_t)countEven;
 121       else meanN=sumOdd/(Float_t)countOdd;
 122       Float_t cntCorr=hrawd->GetBinContent(itb+1,ian+1)-fCMN[ian]*meanN;
 123       hcorrd->SetBinContent(itb+1,ian+1,cntCorr);
 124     }
 125   }
 126
 127   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 128     if(!fGoodAnode[ian]) continue;
 129     Float_t sumQ=0.;
 130      for(Int_t itb=0;itb<tbmax;itb++){
 131       sumQ+=TMath::Power(hcorrd->GetBinContent(itb+1,ian+1)-fBaseline[ian],2);
 132     }
 133     fSumCorrNoise[ian]+=TMath::Sqrt(sumQ/tbmax);
 134   }
 135   delete hcorrd;
 136 }
 137 //______________________________________________________________________
 138 Float_t AliITSOnlineSDDCMN::CalcMeanNoise(){
 139   //
 140   Float_t meanns=0.;
 141   Int_t cnt=0;
 142   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 143     if(!fGoodAnode[ian]) continue;
 144     meanns+=GetAnodeCorrNoise(ian);
 145     cnt++;
 146   }
 147   if(cnt>0) meanns/=(Float_t)cnt;
 148   return meanns;
 149 }
 150 //______________________________________________________________________
 151 void AliITSOnlineSDDCMN::WriteToFXS(){
 152   //
 153   Char_t outfilnam[100];
 154   sprintf(outfilnam,"SDDbase_step2_mod%03d_sid%d.data",fModuleId,fSide);
 155   FILE* outf=fopen(outfilnam,"w");
 156   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 157     fprintf(outf,"%d %d %8.3f %8.3f %8.3f %8.3f\n",ian,IsAnodeGood(ian),GetAnodeBaseline(ian),GetAnodeRawNoise(ian),GetAnodeCommonMode(ian),GetAnodeCorrNoise(ian));
 158   }
 159   fclose(outf);
 160 }