Put 0.1 as minimal value of the gas gain (Raphaelle)
[u/mrichter/AliRoot.git] / T0 / AliT0CalibTimeEq.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, 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
16 /* $Id$ */
17
18 ///////////////////////////////////////////////////////////////////////////////
19 //                                                                           //
20 // class for T0 calibration                       TM-AC-AM_6-02-2006  
21 // equalize time shift for each time CFD channel
22 //                                                                           //
23 ///////////////////////////////////////////////////////////////////////////////
24
25 #include "AliT0CalibTimeEq.h"
26 #include "AliLog.h"
27 #include <TFile.h>
28 #include <TMath.h>
29 #include <TF1.h>
30 #include <TSpectrum.h>
31 #include <TProfile.h>
32 #include <iostream>
33
34 ClassImp(AliT0CalibTimeEq)
35
36 //________________________________________________________________
37   AliT0CalibTimeEq::AliT0CalibTimeEq():TNamed(),
38                                        fMeanVertex(0),        
39                                        fRmsVertex(0)      
40 {
41   //
42
43 }
44
45 //________________________________________________________________
46 AliT0CalibTimeEq::AliT0CalibTimeEq(const char* name):TNamed(),
47                                        fMeanVertex(0),        
48                                        fRmsVertex(0)      
49 {
50   //constructor
51
52   TString namst = "Calib_";
53   namst += name;
54   SetName(namst.Data());
55   SetTitle(namst.Data());
56 }
57
58 //________________________________________________________________
59 AliT0CalibTimeEq::AliT0CalibTimeEq(const AliT0CalibTimeEq& calibda):TNamed(calibda),            
60                                        fMeanVertex(0),        
61                                        fRmsVertex(0)      
62 {
63 // copy constructor
64   SetName(calibda.GetName());
65   SetTitle(calibda.GetName());
66
67
68 }
69
70 //________________________________________________________________
71 AliT0CalibTimeEq &AliT0CalibTimeEq::operator =(const AliT0CalibTimeEq& calibda)
72 {
73 // assignment operator
74   SetName(calibda.GetName());
75   SetTitle(calibda.GetName());
76  
77   return *this;
78 }
79
80 //________________________________________________________________
81 AliT0CalibTimeEq::~AliT0CalibTimeEq()
82 {
83   //
84   // destrictor
85 }
86 //________________________________________________________________
87 void AliT0CalibTimeEq::Reset()
88 {
89   //reset values
90
91   memset(fCFDvalue,0,120*sizeof(Float_t));
92   memset(fTimeEq,1,24*sizeof(Float_t));
93 }
94
95
96 //________________________________________________________________
97 void  AliT0CalibTimeEq::Print(Option_t*) const
98 {
99   // print time values
100
101   printf("\n    ----    PM Arrays       ----\n\n");
102   printf(" Time delay CFD \n");
103   for (Int_t i=0; i<24; i++) printf(" CFD  %f ",fTimeEq[i]);
104   printf("\n Mean Vertex %f \n", fMeanVertex);
105
106
107
108 //________________________________________________________________
109 Bool_t AliT0CalibTimeEq::ComputeOnlineParams(const char* filePhys)
110 {
111   // compute online equalized time
112   Double_t mean=0, meanver=0;
113   Double_t rms=0, rmsver=0;
114   Int_t nent=0;
115   Bool_t ok=false;
116   Int_t npeaks = 20;
117   Int_t sigma=3;
118   Bool_t down=false;
119   Int_t index[20];
120   Int_t nfound=0;
121   gFile = TFile::Open(filePhys);
122   if(!gFile) {
123     AliError("No input PHYS data found ");
124   }
125   else
126     {
127       //    gFile->ls();
128       ok=true;
129       Char_t buf1[30];
130       for (Int_t i=0; i<24; i++)
131         {
132           sprintf(buf1,"CFD1minCFD%d",i+1);
133           TH1F *cfd = (TH1F*) gFile->Get(buf1);
134           if(!cfd) AliWarning(Form("no histograms collected by PHYS DA for channel %i", i));
135           //      printf(" i = %d buf1 = %s\n", i, buf1);
136           if(cfd) {
137             TSpectrum *s = new TSpectrum(2*npeaks,1);
138             nfound = s->Search(cfd,sigma," ",0.1);
139             if(nfound!=0){
140               Float_t *xpeak = s->GetPositionX();
141               TMath::Sort(nfound, xpeak, index,down);
142               Float_t xp = xpeak[index[0]];
143               Double_t hmax = xp+3*sigma;
144               Double_t hmin = xp-3*sigma;
145               cfd->GetXaxis()->SetRangeUser(hmin-20,hmax+20);
146               mean=cfd->GetMean();
147               rms=cfd->GetRMS();
148               nent=cfd->GetEntries();
149               if(nent<500 || rms>20. ) {
150                 ok=false;
151                 AliWarning(Form("Data is not good enouph in PMT %i - mean %f rsm %f nentries %i", i,mean, rms, nent));
152                 
153               }
154             }
155             else 
156               {
157                 ok=false;
158                 AliWarning(Form("Data is not good enouph in PMT %i , no clean peak", i));
159               }
160           }
161           
162         
163            SetTimeEq(i,mean);
164            SetTimeEqRms(i,rms);
165           if (cfd) delete cfd;
166         }
167       TH1F *ver = (TH1F*) gFile->Get("hVertex");
168       if(!ver) AliWarning("no T0 histogram collected by PHYS DA ");
169       if(ver) {
170         meanver = ver->GetMean();
171         rmsver = ver->GetRMS();
172       }
173       SetMeanVertex(meanver);
174       SetRmsVertex(rmsver);
175       
176       gFile->Close();
177       delete gFile;
178
179     }
180     return ok; 
181 }
182
183