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 **************************************************************************/
16 /* $Id: AliT0CalibSeasonTimeShift.cxx 42881 2010-08-16 10:59:14Z alla $ */
18 ///////////////////////////////////////////////////////////////////////////////
20 // class for T0 calibration TM-AC-AM_6-02-2006
21 // equalize time shift for each time CFD channel
23 ///////////////////////////////////////////////////////////////////////////////
25 #include "AliT0CalibSeasonTimeShift.h"
33 ClassImp(AliT0CalibSeasonTimeShift)
35 //________________________________________________________________
36 AliT0CalibSeasonTimeShift::AliT0CalibSeasonTimeShift():TNamed()
39 for (Int_t i=0; i<4; i++)
40 fMeanPar[i] = fSigmaPar[i] = 0;
43 //________________________________________________________________
44 AliT0CalibSeasonTimeShift::AliT0CalibSeasonTimeShift(const char* name):TNamed()
48 TString namst = "Calib_";
50 SetName(namst.Data());
51 SetTitle(namst.Data());
53 for (Int_t i=0; i<4; i++)
54 fMeanPar[i] = fSigmaPar[i] = 0;
58 //________________________________________________________________
59 AliT0CalibSeasonTimeShift::AliT0CalibSeasonTimeShift(const AliT0CalibSeasonTimeShift& calibda):TNamed(calibda)
63 SetName(calibda.GetName());
64 SetTitle(calibda.GetName());
65 ((AliT0CalibSeasonTimeShift &) calibda).Copy(*this);
70 //________________________________________________________________
71 AliT0CalibSeasonTimeShift &AliT0CalibSeasonTimeShift::operator =(const AliT0CalibSeasonTimeShift& calibda)
73 // assignment operator
74 SetName(calibda.GetName());
75 SetTitle(calibda.GetName());
76 if (this != &calibda) ((AliT0CalibSeasonTimeShift &) calibda).Copy(*this);
81 //________________________________________________________________
82 AliT0CalibSeasonTimeShift::~AliT0CalibSeasonTimeShift()
89 //________________________________________________________________
90 void AliT0CalibSeasonTimeShift::Print(Option_t*) const
94 printf("\n ---- T0 results ----\n\n");
95 printf(" (T0A+T0C)/2 = %f; T0A = %f; T0C = %f; resolution = %f \n", fMeanPar[0], fMeanPar[1],fMeanPar[2],fMeanPar[3]);
96 printf(" sigma(T0A+T0C)/2 = %f; sigma(T0 = %f; sigma(T0C) = %f; sigma(resolution) = %f \n" , fSigmaPar[0], fSigmaPar[1], fSigmaPar[2],fSigmaPar[3]);
100 //________________________________________________________________
101 Bool_t AliT0CalibSeasonTimeShift::SetT0Par(Float_t par[4],Float_t spar[4])
104 for (Int_t i=0; i<4; i++)
106 fMeanPar[i] = par[i];
107 fSigmaPar[i] = spar[i];
108 if ( fSigmaPar[i] == 0 || fSigmaPar[i] > 500) ok = false;
113 //________________________________________________________________
114 Int_t AliT0CalibSeasonTimeShift::SetT0Par(const char* filePhys, Float_t *cdbtime)
116 // compute online equalized time
120 TObjArray * tzeroObj = NULL;
122 gFile = TFile::Open(filePhys);
124 AliError("No input PHYS data found ");
129 // TDirectory *dr = (TDirectory*) gFile->Get("T0Calib");
130 tzeroObj = dynamic_cast<TObjArray*>(gFile->Get("T0Calib"));
131 TString histname[4]={"fTzeroORAplusORC", "fTzeroORA", "fTzeroORC", "fResolution"};
132 for (Int_t i=0; i<4; i++)
134 if(cfd) cfd->Reset();
136 cfd = (TH1F*)tzeroObj->FindObject( histname[i].Data());
138 cfd = (TH1F*)gFile ->Get(histname[i].Data());
142 AliError(Form("no histograms collected for %s", histname[i].Data()));
146 if( cfd->GetEntries() == 0) {
148 AliError(Form("%s histogram is empty", histname[i].Data()));
151 GetMeanAndSigma(cfd, mean, sigma);
152 if (sigma == 0 || sigma > 500 || cfd->GetEntries()<200 ){
154 fMeanPar[i] = cdbtime[i];
157 if ( sigma > 0 && sigma < 500 && cfd->GetEntries()>200)
160 fSigmaPar[i] = sigma;
169 //________________________________________________________________________
170 void AliT0CalibSeasonTimeShift::GetMeanAndSigma(TH1F* hist, Float_t &mean, Float_t &sigma) {
172 const double window =2.; //fit window
174 double meanEstimate, sigmaEstimate;
176 maxBin = hist->GetMaximumBin(); //position of maximum
177 meanEstimate = hist->GetBinCenter( maxBin); // mean of gaussian sitting in maximum
178 sigmaEstimate = hist->GetRMS();
179 TF1* fit= new TF1("fit","gaus", meanEstimate - window*sigmaEstimate, meanEstimate + window*sigmaEstimate);
180 fit->SetParameters(hist->GetBinContent(maxBin), meanEstimate, sigmaEstimate);
181 hist->Fit("fit","RQ","Q");
183 mean = (Float_t) fit->GetParameter(1);
184 sigma = (Float_t) fit->GetParameter(2);