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),fLowThreshold(0),fHighThreshold(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 nddl, Int_t ncarlos, Int_t sid):AliITSOnlineSDD(nddl,ncarlos,sid),fNEvents(0),fLowThreshold(0),fHighThreshold(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   TString basfilnam;
  71   basfilnam.Form("SDDbase_step1_ddl%02dc%02d_sid%d.data",fDDL,fCarlos,fSide);
  72   FILE* basf=fopen(basfilnam.Data(),"r");
  73   if(basf==0){
  74     AliWarning(Form("Baseline file not present (ddl %d  carlos %d side %d, Set all baselines to 50\n",fDDL,fCarlos,fSide));
  75     for(Int_t ian=0;ian<fgkNAnodes;ian++){
  76       fBaseline[ian]=50.;
  77       fEqBaseline[ian]=50;
  78       fOffsetBaseline[ian]=0;
  79       fGoodAnode[ian]=1;
  80     }
  81     return;
  82   }
  83   fscanf(basf,"%d\n",&fHighThreshold);
  84   fscanf(basf,"%d\n",&fLowThreshold);
  85   Int_t n,ok,eqbase,offbase;
  86   Float_t base,rms,cmn,corrnoi;
  87   for(Int_t ian=0;ian<fgkNAnodes;ian++){
  88     fscanf(basf,"%d %d %f %d %d %f %f %f\n",&n,&ok,&base,&eqbase,&offbase,&rms,&cmn,&corrnoi);
  89     fGoodAnode[ian]=ok;
  90     fBaseline[ian]=base;
  91     fEqBaseline[ian]=eqbase;
  92     fOffsetBaseline[ian]=offbase;
  93     fRawNoise[ian]=rms;
  94     fCMN[ian]=cmn;
  95   }
  96   fclose(basf);
  97 }
  98 //______________________________________________________________________
  99 void  AliITSOnlineSDDCMN::ValidateAnodes(){
 100   //
 101   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 102     if(!fGoodAnode[ian]) continue;
 103     if(GetAnodeCorrNoise(ian)>fMaxCorrNoise || GetAnodeCorrNoise(ian)<fMinCorrNoise) fGoodAnode[ian]=0;
 104     if(GetAnodeCorrNoise(ian)>fNSigmaNoise*CalcMeanNoise()) fGoodAnode[ian]=0;
 105   }
 106 }
 107
 108 //______________________________________________________________________
 109 TH2F* AliITSOnlineSDDCMN::GetCleanEvent(TH2F* hrawd) const {
 110   //
 111   TH2F* hcorrd=new TH2F("hcorrd","",hrawd->GetNbinsX(),hrawd->GetXaxis()->GetXmin(),hrawd->GetXaxis()->GetXmax(),hrawd->GetNbinsY(),hrawd->GetYaxis()->GetXmin(),hrawd->GetYaxis()->GetXmax());
 112   for(Int_t itb=fFirstGoodTB;itb<=fLastGoodTB;itb++){
 113     Float_t sumEven=0., sumOdd=0.;
 114     Int_t countEven=0, countOdd=0;
 115     for(Int_t ian=0;ian<fgkNAnodes;ian+=2){
 116       if(!fGoodAnode[ian]) continue;
 117       sumEven+=hrawd->GetBinContent(itb+1,ian+1)-fBaseline[ian];
 118       countEven++;
 119     }
 120     for(Int_t ian=1;ian<fgkNAnodes;ian+=2){
 121       if(!fGoodAnode[ian]) continue;
 122       sumOdd+=hrawd->GetBinContent(itb+1,ian+1)-fBaseline[ian];
 123       countOdd++;
 124     }
 125     for(Int_t ian=0;ian<fgkNAnodes;ian++){
 126       if(!fGoodAnode[ian]) continue;
 127       Float_t meanN;
 128       if(ian%2==0) meanN=sumEven/(Float_t)countEven;
 129       else meanN=sumOdd/(Float_t)countOdd;
 130       Float_t cntCorr=hrawd->GetBinContent(itb+1,ian+1)-fCMN[ian]*meanN;
 131       hcorrd->SetBinContent(itb+1,ian+1,cntCorr);
 132     }
 133   }
 134   return hcorrd;
 135 }
 136 //______________________________________________________________________
 137 void AliITSOnlineSDDCMN::AddEvent(TH2F* hrawd){
 138   //
 139   fNEvents++;
 140   TH2F* hcorrd=GetCleanEvent(hrawd);
 141
 142   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 143     if(!fGoodAnode[ian]) continue;
 144     Float_t sumQ=0.;
 145     Int_t cnt=0;
 146     for(Int_t itb=fFirstGoodTB;itb<=fLastGoodTB;itb++){
 147       Float_t cntdiff=hcorrd->GetBinContent(itb+1,ian+1)-fBaseline[ian];
 148       sumQ+=cntdiff*cntdiff;
 149       cnt++;
 150     }
 151     fSumCorrNoise[ian]+=TMath::Sqrt(sumQ/(Float_t)cnt);
 152   }
 153   delete hcorrd;
 154 }
 155 //______________________________________________________________________
 156 Float_t AliITSOnlineSDDCMN::CalcMeanNoise() const{
 157   //
 158   Float_t meanns=0.;
 159   Int_t cnt=0;
 160   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 161     if(!fGoodAnode[ian]) continue;
 162     meanns+=GetAnodeCorrNoise(ian);
 163     cnt++;
 164   }
 165   if(cnt>0) meanns/=(Float_t)cnt;
 166   return meanns;
 167 }
 168 //______________________________________________________________________
 169 void AliITSOnlineSDDCMN::WriteToASCII(){
 170   //
 171   TString outfilnam;
 172   outfilnam.Form("SDDbase_step2_ddl%02dc%02d_sid%d.data",fDDL,fCarlos,fSide);
 173   FILE* outf=fopen(outfilnam.Data(),"w");
 174   fprintf(outf,"%d\n",fHighThreshold);
 175   fprintf(outf,"%d\n",fLowThreshold);
 176   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 177     fprintf(outf,"%d %d %8.3f %d %d %8.3f %8.3f %8.3f\n",ian,IsAnodeGood(ian),GetAnodeBaseline(ian),GetAnodeEqualizedBaseline(ian),GetAnodeBaselineOffset(ian),GetAnodeRawNoise(ian),GetAnodeCommonMode(ian),GetAnodeCorrNoise(ian));
 178   }
 179   fclose(outf);
 180 }
 181
 182 //______________________________________________________________________
 183 TH1F* AliITSOnlineSDDCMN::GetBaselineAnodeHisto() const {
 184   //
 185   TString hisnam;
 186   hisnam.Form("hbase%02dc%02ds%d",fDDL,fCarlos,fSide);
 187   TH1F* h=new TH1F(hisnam.Data(),"",256,-0.5,255.5);
 188   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 189     h->SetBinContent(ian+1,GetAnodeBaseline(ian));
 190   }
 191   return h;
 192 }
 193 //______________________________________________________________________
 194 TH1F* AliITSOnlineSDDCMN::GetRawNoiseAnodeHisto() const {
 195   //
 196   TString hisnam;
 197   hisnam.Form("hnois%02dc%02ds%d",fDDL,fCarlos,fSide);
 198   TH1F* h=new TH1F(hisnam.Data(),"",256,-0.5,255.5);
 199   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 200     h->SetBinContent(ian+1,GetAnodeRawNoise(ian));
 201   }
 202   return h;
 203 }
 204 //______________________________________________________________________
 205 TH1F* AliITSOnlineSDDCMN::GetCorrNoiseAnodeHisto() const {
 206   //
 207   TString hisnam;
 208   hisnam.Form("hcorn%02dc%02ds%d",fDDL,fCarlos,fSide);
 209   TH1F* h=new TH1F(hisnam.Data(),"",256,-0.5,255.5);
 210   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 211     h->SetBinContent(ian+1,GetAnodeCorrNoise(ian));
 212   }
 213   return h;
 214 }
 215 //______________________________________________________________________
 216 TH1F* AliITSOnlineSDDCMN::GetCMNCoefAnodeHisto() const {
 217 //
 218   TString hisnam;
 219   hisnam.Form("hcmn%02dc%02ds%d",fDDL,fCarlos,fSide);
 220   TH1F* h=new TH1F(hisnam.Data(),"",256,-0.5,255.5);
 221   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 222     h->SetBinContent(ian+1,GetAnodeCommonMode(ian));
 223   }
 224   return h;
 225 }
 226 //______________________________________________________________________
 227 TH1F* AliITSOnlineSDDCMN::GetStatusAnodeHisto() const {
 228 //
 229   TString hisnam;
 230   hisnam.Form("hgood%02dc%02ds%d",fDDL,fCarlos,fSide);
 231   TH1F* h=new TH1F(hisnam.Data(),"",256,-0.5,255.5);
 232   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 233     h->SetBinContent(ian+1,float(IsAnodeGood(ian)));
 234   }
 235   return h;
 236 }
 237 //______________________________________________________________________
 238 TH1F* AliITSOnlineSDDCMN::GetBaselineHisto() const {
 239   //
 240   TString hisnam;
 241   hisnam.Form("hdbd%02dc%02ds%d",fDDL,fCarlos,fSide);
 242   TH1F* h=new TH1F(hisnam.Data(),"",100,0.,150.);
 243   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 244     h->Fill(GetAnodeBaseline(ian));
 245   }
 246   return h;
 247 }
 248 //______________________________________________________________________
 249 TH1F* AliITSOnlineSDDCMN::GetRawNoiseHisto() const {
 250   //
 251   TString hisnam;
 252   hisnam.Form("hdnd%02dc%02ds%d",fDDL,fCarlos,fSide);
 253   TH1F* h=new TH1F(hisnam.Data(),"",100,0.,8.);
 254   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 255     h->Fill(GetAnodeRawNoise(ian));
 256   }
 257   return h;
 258 }
 259 //______________________________________________________________________
 260 TH1F* AliITSOnlineSDDCMN::GetCorrNoiseHisto() const {
 261   //
 262   TString hisnam;
 263   hisnam.Form("hdcd%02dc%02ds%d",fDDL,fCarlos,fSide);
 264   TH1F* h=new TH1F(hisnam.Data(),"",100,0.,8.);
 265   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 266     h->Fill(GetAnodeCorrNoise(ian));
 267   }
 268   return h;
 269 }
 270 //______________________________________________________________________
 271 Bool_t AliITSOnlineSDDCMN::WriteToROOT(TFile *fil){
 272   //
 273   if(fil==0){
 274     AliWarning("Invalid pointer to ROOT file");
 275     return kFALSE;
 276   }
 277   TString hisnam;
 278   fil->cd();
 279   hisnam.Form("hgood%02dc%02ds%d",fDDL,fCarlos,fSide);
 280   TH1F hgood(hisnam.Data(),"",256,-0.5,255.5);
 281   hisnam.Form("hbase%02dc%02ds%d",fDDL,fCarlos,fSide);
 282   TH1F hbase(hisnam.Data(),"",256,-0.5,255.5);
 283   hisnam.Form("hnois%02dc%02ds%d",fDDL,fCarlos,fSide);
 284   TH1F hnois(hisnam.Data(),"",256,-0.5,255.5);
 285   hisnam.Form("hcmn%02dc%02ds%d",fDDL,fCarlos,fSide);
 286   TH1F hcmn(hisnam.Data(),"",256,-0.5,255.5);
 287   hisnam.Form("hcorn%02dc%02ds%d",fDDL,fCarlos,fSide);
 288   TH1F hcorn(hisnam.Data(),"",256,-0.5,255.5);
 289   for(Int_t ian=0;ian<fgkNAnodes;ian++){
 290     hgood.SetBinContent(ian+1,float(IsAnodeGood(ian)));
 291     hbase.SetBinContent(ian+1,GetAnodeBaseline(ian));
 292     hnois.SetBinContent(ian+1,GetAnodeRawNoise(ian));
 293     hcmn.SetBinContent(ian+1,GetAnodeCommonMode(ian));
 294     hcorn.SetBinContent(ian+1,GetAnodeCorrNoise(ian));
 295   }
 296   hgood.Write();
 297   hbase.Write();
 298   hnois.Write();
 299   hcmn.Write();
 300   hcorn.Write();
 301   return kTRUE;
 302 }