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"
25 #include "AliT0LookUpValue.h"
29 #include <TObjArray.h>
36 #include <Riostream.h>
37 #include <TSpectrum.h>
38 #include <TVirtualFitter.h>
43 ClassImp(AliT0CalibWalk)
45 //________________________________________________________________
46 AliT0CalibWalk::AliT0CalibWalk(): TNamed(),
53 //________________________________________________________________
54 AliT0CalibWalk::AliT0CalibWalk(const char* name):TNamed(),
58 TString namst = "Calib_";
60 SetName(namst.Data());
61 SetTitle(namst.Data());
65 //________________________________________________________________
66 AliT0CalibWalk::AliT0CalibWalk(const AliT0CalibWalk& calibda) :
73 SetName(calibda.GetName());
74 SetTitle(calibda.GetName());
79 //________________________________________________________________
80 AliT0CalibWalk &AliT0CalibWalk::operator =(const AliT0CalibWalk& calibda)
82 // assignment operator
83 SetName(calibda.GetName());
84 SetTitle(calibda.GetName());
89 //________________________________________________________________
90 AliT0CalibWalk::~AliT0CalibWalk()
94 //________________________________________________________________
95 void AliT0CalibWalk::SetWalk(Int_t ipmt)
99 Int_t x[70000], y[70000], index[70000];
100 Float_t time[10000],amplitude[10000];
104 const char * filename = gSystem->ExpandPathName("$ALICE_ROOT/T0/data/CFD-Amp.txt");
105 ifstream inFile(filename);
106 if(!inFile) AliError(Form("Cannot open file %s !",filename));
109 while(getline(inFile,buffer)){
116 cout<<" number of data "<<i<<endl;
118 TMath::Sort(i, y, index,down);
119 Int_t amp=0, iin=0, isum=0, sum=0;
121 for (Int_t ii=0; ii<i; ii++)
132 time[isum] = Float_t (sum/(iin));
134 time[isum] =Float_t (x[ind]);
135 amplitude[isum] = Float_t (amp);
145 TGraph* gr = new TGraph(isum, amplitude, time);
146 fWalk.AddAtAndExpand(gr,ipmt);
149 //________________________________________________________________
151 void AliT0CalibWalk::SetAmpLEDRec(Int_t ipmt)
154 Float_t x[100], y[100];
157 const char * filename = gSystem->ExpandPathName("$ALICE_ROOT/T0/data/CFD-LED.txt");
158 ifstream inFile(filename);
159 if(!inFile) {AliError(Form("Cannot open file %s !",filename));}
164 while(getline(inFile,buffer)){
170 Float_t y1[100], x1[100];
171 for (Int_t ir=0; ir<i; ir++){
172 y1[ir]=y[i-ir]; x1[ir]=x[i-ir];}
173 TGraph* gr = new TGraph(i,y1,x1);
174 fAmpLEDRec.AddAtAndExpand(gr,ipmt);
178 //________________________________________________________________
180 void AliT0CalibWalk::MakeWalkCorrGraph(const char *laserFile)
182 TFile *gFile = TFile::Open(laserFile);
183 gSystem->Load("libSpectrum");
190 Char_t buf1[10], buf2[10],title[10];
192 for (Int_t i=0; i<24; i++)
194 sprintf(buf1,"T0_C_%i_CFD",i+1);
195 sprintf(buf2,"CFD_QTC%i",i+1);
196 // cout<<buf1<<" "<<buf2<<endl;
197 TH2F *qtc_cfd = (TH2F*) gFile->Get(buf2);
198 TH1F *cfd = (TH1F*) gFile->Get(buf1);
199 TSpectrum *s = new TSpectrum(2*npeaks,1.);
200 Int_t nfound = s->Search(cfd,sigma," ",0.05);
201 // cout<<"Found "<<nfound<<" peaks sigma "<<sigma<<endl;
204 Float_t *xpeak = s->GetPositionX();
205 TMath::Sort(nfound, xpeak, index,down);
206 Float_t xp = xpeak[index[0]];
207 Float_t hmax = xp+10*sigma;
208 Float_t hmin = xp-10*sigma;
209 Int_t nbins= qtc_cfd->GetXaxis()->GetNbins();
210 TProfile *pr_y = qtc_cfd->ProfileX();
211 pr_y->SetMaximum(hmax);
212 pr_y->SetMinimum(hmin);
214 Double_t *xx = new Double_t[np];
215 Double_t *yy = new Double_t[np];
218 for (Int_t ip=1; ip<nbins; ip++)
222 if (pr_y->GetBinContent(ip) !=0)
224 yg +=pr_y->GetBinContent(ip);
230 xx[ip/20] = Float_t (pr_y->GetBinCenter(ip));
236 sprintf(title,"Walk %i",i+1);
237 gr[i] = new TGraph(np,xx,yy);
238 gr[i]->SetTitle(title);
239 gr[i]->SetMinimum(hmin);
240 gr[i]->SetMaximum(hmax);
241 gr[i]->SetMarkerStyle(7);
242 fWalk.AddAtAndExpand(gr[i],i);