container for season time shift
[u/mrichter/AliRoot.git] / T0 / AliT0CalibSeasonTimeShift.cxx
CommitLineData
455957bc 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: AliT0CalibSeasonTimeShift.cxx 42881 2010-08-16 10:59:14Z alla $ */
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 "AliT0CalibSeasonTimeShift.h"
26#include "AliLog.h"
27#include <TFile.h>
28#include <TMath.h>
29#include <TF1.h>
30#include <TProfile.h>
31#include <iostream>
32
33ClassImp(AliT0CalibSeasonTimeShift)
34
35//________________________________________________________________
36 AliT0CalibSeasonTimeShift::AliT0CalibSeasonTimeShift():TNamed()
37{
38 //
39
40}
41
42//________________________________________________________________
43AliT0CalibSeasonTimeShift::AliT0CalibSeasonTimeShift(const char* name):TNamed()
44{
45 //constructor
46
47 TString namst = "Calib_";
48 namst += name;
49 SetName(namst.Data());
50 SetTitle(namst.Data());
51}
52
53//________________________________________________________________
54AliT0CalibSeasonTimeShift::AliT0CalibSeasonTimeShift(const AliT0CalibSeasonTimeShift& calibda):TNamed(calibda)
55
56{
57// copy constructor
58 SetName(calibda.GetName());
59 SetTitle(calibda.GetName());
60
61
62}
63
64//________________________________________________________________
65AliT0CalibSeasonTimeShift &AliT0CalibSeasonTimeShift::operator =(const AliT0CalibSeasonTimeShift& calibda)
66{
67// assignment operator
68 SetName(calibda.GetName());
69 SetTitle(calibda.GetName());
70
71 return *this;
72}
73
74//________________________________________________________________
75AliT0CalibSeasonTimeShift::~AliT0CalibSeasonTimeShift()
76{
77 //
78 // destrictor
79}
80
81
82//________________________________________________________________
83void AliT0CalibSeasonTimeShift::Print(Option_t*) const
84{
85 // print time values
86
87 printf("\n ---- T0 results ----\n\n");
88 printf(" (T0A+T0C)/2 = %f; T0A = %f; T0C = %f; resolution = %f \n", fMeanPar[0], fMeanPar[1],fMeanPar[2],fMeanPar[3]);
89 printf(" sigma(T0A+T0C)/2 = %f; sigma(T0 = %f; sigma(T0C) = %f; sigma(resolution) = %f \n" , fSigmaPar[0], fSigmaPar[1], fSigmaPar[2],fSigmaPar[3]);
90
91}
92
93//________________________________________________________________
94void AliT0CalibSeasonTimeShift::SetT0Par(Float_t par[4],Float_t spar[4])
95{
96 for (Int_t i=0; i<4; i++)
97 {
98 fMeanPar[i] = par[i];
99 fSigmaPar[i] = spar[i];
100 }
101
102}
103
104//________________________________________________________________
105void AliT0CalibSeasonTimeShift::SetT0Par(const char* filePhys)
106{
107 // compute online equalized time
108 Float_t mean[4], sigma[4];
109
110 gFile = TFile::Open(filePhys);
111 if(!gFile) {
112 AliError("No input PHYS data found ");
113 }
114 else
115 {
116 gFile->ls();
117
118 TString histname[4]={"meanAC", "meanA", "meanC", "resolution"};
119 for (Int_t i=0; i<4; i++)
120 {
121 TH1F *cfd = (TH1F*) gFile->Get(histname[i].Data());
122 if(!cfd) AliWarning(Form("no histograms collected for %s", histname[i].Data()));
123 if(cfd) {
124 TF1 *g = new TF1("g", "gaus",-2,2);
125 cfd->Fit("g"," ","Q",-2,2);
126 Double_t par[3];
127 g->GetParameters(&par[0]);
128 fMeanPar[i] = par[1];
129 fSigmaPar[i]=par[2];
130
131 }
132 }
133
134 gFile->Close();
135 delete gFile;
136
137 }
138
139}