1 /**************************************************************************
2 * Copyright(c) 1998-1999, 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 **************************************************************************/
18 ///////////////////////////////////////////////////////////////////////////////
20 // class T0 walk correction Alla Maevskaya alla@inr.ru 20.11.2007 //
22 ///////////////////////////////////////////////////////////////////////////////
24 #include "AliT0CalibWalk.h"
27 #include <TObjArray.h>
33 #include <Riostream.h>
34 #include <TSpectrum.h>
38 ClassImp(AliT0CalibWalk)
40 //________________________________________________________________
41 AliT0CalibWalk::AliT0CalibWalk(): TNamed(),
48 //________________________________________________________________
49 AliT0CalibWalk::AliT0CalibWalk(const char* name):TNamed(),
53 TString namst = "Calib_";
55 SetName(namst.Data());
56 SetTitle(namst.Data());
60 //________________________________________________________________
61 AliT0CalibWalk::AliT0CalibWalk(const AliT0CalibWalk& calibda) :
68 SetName(calibda.GetName());
69 SetTitle(calibda.GetName());
74 //________________________________________________________________
75 AliT0CalibWalk &AliT0CalibWalk::operator =(const AliT0CalibWalk& calibda)
77 // assignment operator
78 SetName(calibda.GetName());
79 SetTitle(calibda.GetName());
84 //________________________________________________________________
85 AliT0CalibWalk::~AliT0CalibWalk()
89 //________________________________________________________________
90 void AliT0CalibWalk::SetWalk(Int_t ipmt)
92 //read QTC walk graph from external file
95 Int_t x[70000], y[70000], index[70000];
96 Float_t time[10000],amplitude[10000];
100 const char * filename = gSystem->ExpandPathName("$ALICE_ROOT/T0/data/CFD-Amp.txt");
101 ifstream inFile(filename);
102 if(!inFile) AliError(Form("Cannot open file %s !",filename));
105 while(getline(inFile,buffer)){
112 cout<<" number of data "<<i<<endl;
114 TMath::Sort(i, y, index,down);
115 Int_t amp=0, iin=0, isum=0, sum=0;
117 for (Int_t ii=0; ii<i; ii++)
128 time[isum] = Float_t (sum/(iin));
130 time[isum] =Float_t (x[ind]);
131 amplitude[isum] = Float_t (amp);
141 TGraph* gr = new TGraph(isum, amplitude, time);
142 fWalk.AddAtAndExpand(gr,ipmt);
145 //________________________________________________________________
147 void AliT0CalibWalk::SetAmpLEDRec(Int_t ipmt)
149 // read LED walk from external file
151 Float_t x[100], y[100];
154 const char * filename = gSystem->ExpandPathName("$ALICE_ROOT/T0/data/CFD-LED.txt");
155 ifstream inFile(filename);
156 if(!inFile) {AliError(Form("Cannot open file %s !",filename));}
161 while(getline(inFile,buffer)){
167 Float_t y1[100], x1[100];
168 for (Int_t ir=0; ir<i; ir++){
169 y1[ir]=y[i-ir]; x1[ir]=x[i-ir];}
170 TGraph* gr = new TGraph(i,y1,x1);
171 fAmpLEDRec.AddAtAndExpand(gr,ipmt);
175 //________________________________________________________________
177 void AliT0CalibWalk::MakeWalkCorrGraph(const char *laserFile)
179 //make walk corerction for preprocessor
181 TFile *gFile = TFile::Open(laserFile);
182 //gSystem->Load("libSpectrum");
190 Char_t buf1[10], buf2[10], buf3[10], title[10], title2[10], titleLED[10], title2LED[10];
191 TSpectrum *s = new TSpectrum(2*npeaks,1.);
193 for (Int_t d=0; d<2; d++)
195 for (Int_t i=0; i<12; i++)
197 sprintf(buf1,"T0_C_%i_CFD",i+1);
200 sprintf(buf2,"CFD_QTC%i",i+1);
204 sprintf(buf3,"CFD_LED%i",i+1);
206 TH2F *qtccfd = (TH2F*) gFile->Get(buf2);
207 TH2F *ledcfd = (TH2F*) gFile->Get(buf3);
208 // cout<<buf1<<" "<<buf2<<endl;
209 TH1F *cfd = (TH1F*) gFile->Get(buf1);
210 Int_t nfound = s->Search(cfd,sigma,"goff",0.05);
211 // cout<<"Found "<<nfound<<" peaks sigma "<<sigma<<endl;
214 Float_t *xpeak = s->GetPositionX();
215 TMath::Sort(nfound, xpeak, index,down);
216 Float_t xp = xpeak[index[0]];
217 Float_t hmax = xp+10*sigma;
218 Float_t hmin = xp-10*sigma;
221 Int_t nbins= qtccfd->GetXaxis()->GetNbins();
222 TProfile *prY = qtccfd->ProfileX();
223 prY->SetMaximum(hmax);
224 prY->SetMinimum(hmin);
226 Double_t *xx = new Double_t[np];
227 Double_t *yy = new Double_t[np];
230 for (Int_t ip=1; ip<nbins; ip++)
234 if (prY->GetBinContent(ip) !=0)
236 yg +=prY->GetBinContent(ip);
242 xx[ip/20] = Float_t (prY->GetBinCenter(ip));
248 sprintf(title,"Walk %i",i+1);
249 gr[i] = new TGraph(np,xx,yy);
250 gr[i]->SetTitle(title);
251 gr[i]->SetMinimum(hmin);
252 gr[i]->SetMaximum(hmax);
253 gr[i]->SetMarkerStyle(7);
254 sprintf(title2,"Walk %i",i+13);
255 gr[i+12] = new TGraph(np,xx,yy);
256 gr[i+12]->SetTitle(title2);
257 gr[i+12]->SetMinimum(hmin);
258 gr[i+12]->SetMaximum(hmax);
259 gr[i+12]->SetMarkerStyle(7);
261 fWalk.AddAtAndExpand(gr[i],i);
262 fWalk.AddAtAndExpand(gr[i+12],i+12);
269 Int_t nbinsLED= ledcfd->GetXaxis()->GetNbins();
270 TProfile *prYLED = ledcfd->ProfileX();
271 prYLED->SetMaximum(hmax);
272 prYLED->SetMinimum(hmin);
273 Int_t npLED=nbinsLED/20;
274 Double_t *xxLED = new Double_t[npLED];
275 Double_t *yyLED = new Double_t[npLED];
278 for (Int_t ip=1; ip<nbinsLED; ip++)
282 if (prYLED->GetBinContent(ip) !=0)
284 ygLED +=prYLED->GetBinContent(ip);
290 xxLED[ip/20] = Float_t (prYLED->GetBinCenter(ip));
291 yyLED[ip/20] = ygLED/ngLED;
296 sprintf(titleLED,"Walk LED %i",i+1);
297 grLED[i] = new TGraph(npLED,xxLED,yyLED);
298 grLED[i]->SetTitle(titleLED);
299 grLED[i]->SetMinimum(hmin);
300 grLED[i]->SetMaximum(hmax);
301 grLED[i]->SetMarkerStyle(7);
302 sprintf(title2LED,"Walk LED%i",i+13);
303 grLED[i+12] = new TGraph(npLED,xxLED,yyLED);
304 grLED[i+12]->SetTitle(title2LED);
305 grLED[i+12]->SetMinimum(hmin);
306 grLED[i+12]->SetMaximum(hmax);
307 grLED[i+12]->SetMarkerStyle(7);
309 fAmpLEDRec.AddAtAndExpand(grLED[i],i);
310 fAmpLEDRec.AddAtAndExpand(grLED[i+12],i+12);