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