1 /**************************************************************************
2 * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
16 #include "AliITSOnlineSDDCMN.h"
22 ///////////////////////////////////////////////////////////////////
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 //
28 ///////////////////////////////////////////////////////////////////
31 ClassImp(AliITSOnlineSDDCMN)
32 //______________________________________________________________________
33 AliITSOnlineSDDCMN::AliITSOnlineSDDCMN():AliITSOnlineSDD(),fNEvents(0),fMinCorrNoise(0.),fMaxCorrNoise(0.),fNSigmaNoise(0.)
35 // default constructor
41 //______________________________________________________________________
42 AliITSOnlineSDDCMN::AliITSOnlineSDDCMN(Int_t nddl, Int_t ncarlos, Int_t sid):AliITSOnlineSDD(nddl,ncarlos,sid),fNEvents(0),fMinCorrNoise(0.),fMaxCorrNoise(0.),fNSigmaNoise(0.)
44 // default constructor
50 //______________________________________________________________________
51 AliITSOnlineSDDCMN::~AliITSOnlineSDDCMN(){
54 //______________________________________________________________________
55 void AliITSOnlineSDDCMN::Reset(){
58 for(Int_t i=0;i<fgkNAnodes;i++){
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_ddl%02dc%02d_sid%d.data",fDDL,fCarlos,fSide);
72 FILE* basf=fopen(basfilnam,"r");
74 AliWarning("Baseline file not present, Set all baselines to 50\n");
75 for(Int_t ian=0;ian<fgkNAnodes;ian++){
78 fOffsetBaseline[ian]=0;
83 Int_t n,ok,eqbase,offbase;
84 Float_t base,rms,cmn,corrnoi;
85 for(Int_t ian=0;ian<fgkNAnodes;ian++){
86 fscanf(basf,"%d %d %f %d %d %f %f %f\n",&n,&ok,&base,&eqbase,&offbase,&rms,&cmn,&corrnoi);
89 fEqBaseline[ian]=eqbase;
90 fOffsetBaseline[ian]=offbase;
96 //______________________________________________________________________
97 void AliITSOnlineSDDCMN::ValidateAnodes(){
99 for(Int_t ian=0;ian<fgkNAnodes;ian++){
100 if(!fGoodAnode[ian]) continue;
101 if(GetAnodeCorrNoise(ian)>fMaxCorrNoise || GetAnodeCorrNoise(ian)<fMinCorrNoise) fGoodAnode[ian]=0;
102 if(GetAnodeCorrNoise(ian)>fNSigmaNoise*CalcMeanNoise()) fGoodAnode[ian]=0;
106 //______________________________________________________________________
107 void AliITSOnlineSDDCMN::AddEvent(TH2F* hrawd){
110 const Int_t kTimeBins=fLastGoodTB-fFirstGoodTB+1;
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];
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];
125 for(Int_t ian=0;ian<fgkNAnodes;ian++){
126 if(!fGoodAnode[ian]) continue;
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);
135 for(Int_t ian=0;ian<fgkNAnodes;ian++){
136 if(!fGoodAnode[ian]) continue;
138 for(Int_t itb=fFirstGoodTB;itb<=fLastGoodTB;itb++){
139 sumQ+=TMath::Power(hcorrd->GetBinContent(itb+1,ian+1)-fBaseline[ian],2);
141 fSumCorrNoise[ian]+=TMath::Sqrt(sumQ/(Float_t)kTimeBins);
145 //______________________________________________________________________
146 Float_t AliITSOnlineSDDCMN::CalcMeanNoise() const{
150 for(Int_t ian=0;ian<fgkNAnodes;ian++){
151 if(!fGoodAnode[ian]) continue;
152 meanns+=GetAnodeCorrNoise(ian);
155 if(cnt>0) meanns/=(Float_t)cnt;
158 //______________________________________________________________________
159 void AliITSOnlineSDDCMN::WriteToASCII(){
161 Char_t outfilnam[100];
162 sprintf(outfilnam,"SDDbase_step2_ddl%02dc%02d_sid%d.data",fDDL,fCarlos,fSide);
163 FILE* outf=fopen(outfilnam,"w");
164 for(Int_t ian=0;ian<fgkNAnodes;ian++){
165 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));
170 //______________________________________________________________________
171 TH1F* AliITSOnlineSDDCMN::GetBaselineAnodeHisto() const {
174 sprintf(hisnam,"hbd%02dc%02ds%d",fDDL,fCarlos,fSide);
175 TH1F* h=new TH1F(hisnam,"",256,-0.5,255.5);
176 for(Int_t ian=0;ian<fgkNAnodes;ian++){
177 h->SetBinContent(ian+1,GetAnodeBaseline(ian));
181 //______________________________________________________________________
182 TH1F* AliITSOnlineSDDCMN::GetRawNoiseAnodeHisto() const {
185 sprintf(hisnam,"hnd%02dc%02ds%d",fDDL,fCarlos,fSide);
186 TH1F* h=new TH1F(hisnam,"",256,-0.5,255.5);
187 for(Int_t ian=0;ian<fgkNAnodes;ian++){
188 h->SetBinContent(ian+1,GetAnodeRawNoise(ian));
192 //______________________________________________________________________
193 TH1F* AliITSOnlineSDDCMN::GetCorrNoiseAnodeHisto() const {
196 sprintf(hisnam,"hcd%02dc%02ds%d",fDDL,fCarlos,fSide);
197 TH1F* h=new TH1F(hisnam,"",256,-0.5,255.5);
198 for(Int_t ian=0;ian<fgkNAnodes;ian++){
199 h->SetBinContent(ian+1,GetAnodeCorrNoise(ian));
203 //______________________________________________________________________
204 TH1F* AliITSOnlineSDDCMN::GetBaselineHisto() const {
207 sprintf(hisnam,"hdbd%02dc%02ds%d",fDDL,fCarlos,fSide);
208 TH1F* h=new TH1F(hisnam,"",100,0.,150.);
209 for(Int_t ian=0;ian<fgkNAnodes;ian++){
210 h->Fill(GetAnodeBaseline(ian));
214 //______________________________________________________________________
215 TH1F* AliITSOnlineSDDCMN::GetRawNoiseHisto() const {
218 sprintf(hisnam,"hdnd%02dc%02ds%d",fDDL,fCarlos,fSide);
219 TH1F* h=new TH1F(hisnam,"",100,0.,8.);
220 for(Int_t ian=0;ian<fgkNAnodes;ian++){
221 h->Fill(GetAnodeRawNoise(ian));
225 //______________________________________________________________________
226 TH1F* AliITSOnlineSDDCMN::GetCorrNoiseHisto() const {
229 sprintf(hisnam,"hdcd%02dc%02ds%d",fDDL,fCarlos,fSide);
230 TH1F* h=new TH1F(hisnam,"",100,0.,8.);
231 for(Int_t ian=0;ian<fgkNAnodes;ian++){
232 h->Fill(GetAnodeCorrNoise(ian));
236 //______________________________________________________________________
237 Bool_t AliITSOnlineSDDCMN::WriteToROOT(TFile *fil){
240 AliWarning("Invalid pointer to ROOT file");
245 sprintf(hisnam,"hgood%02dc%02ds%d",fDDL,fCarlos,fSide);
246 TH1F hgood(hisnam,"",256,-0.5,255.5);
247 sprintf(hisnam,"hbase%02dc%02ds%d",fDDL,fCarlos,fSide);
248 TH1F hbase(hisnam,"",256,-0.5,255.5);
249 sprintf(hisnam,"hnois%02dc%02ds%d",fDDL,fCarlos,fSide);
250 TH1F hnois(hisnam,"",256,-0.5,255.5);
251 sprintf(hisnam,"hcmn%02dc%02ds%d",fDDL,fCarlos,fSide);
252 TH1F hcmn(hisnam,"",256,-0.5,255.5);
253 sprintf(hisnam,"hcorn%02dc%02ds%d",fDDL,fCarlos,fSide);
254 TH1F hcorn(hisnam,"",256,-0.5,255.5);
255 for(Int_t ian=0;ian<fgkNAnodes;ian++){
256 hgood.SetBinContent(ian+1,float(IsAnodeGood(ian)));
257 hbase.SetBinContent(ian+1,GetAnodeBaseline(ian));
258 hnois.SetBinContent(ian+1,GetAnodeRawNoise(ian));
259 hcmn.SetBinContent(ian+1,GetAnodeCommonMode(ian));
260 hcorn.SetBinContent(ian+1,GetAnodeCorrNoise(ian));