Some warning going to error and viceversa:
[u/mrichter/AliRoot.git] / STEER / STEER / AliLHCClockPhase.cxx
CommitLineData
8c1f1006 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
30d5d647 16#include <Riostream.h>
8c1f1006 17
18#include "AliLHCClockPhase.h"
19#include "AliDCSValue.h"
20#include "AliLog.h"
21
66b0310c 22using std::endl;
23using std::cout;
8c1f1006 24ClassImp(AliLHCClockPhase)
25
26//______________________________________________________________________________
27AliLHCClockPhase::AliLHCClockPhase():
28 TObject(),
29 fPhaseB1(),
30 fPhaseB2()
31{
32 // default constructor
33 // ...
34 fPhaseB1.SetOwner();
35 fPhaseB2.SetOwner();
36}
37
38//______________________________________________________________________________
39void AliLHCClockPhase::AddPhaseB1DP(UInt_t timestamp, Float_t phase)
40{
41 // Add a phase beam1 measurement
42 // to the array of data-points
43
44 fPhaseB1.AddLast(new AliDCSValue(phase,timestamp));
45}
46
47//______________________________________________________________________________
48void AliLHCClockPhase::AddPhaseB2DP(UInt_t timestamp, Float_t phase)
49{
50 // Add a phase beam2 measurement
51 // to the array of data-points
52
53 fPhaseB2.AddLast(new AliDCSValue(phase,timestamp));
54}
55
56//______________________________________________________________________________
57const AliDCSValue* AliLHCClockPhase::GetPhaseB1DP(Int_t index) const
58{
59 // Get the value of the phase
60 // The argument is the DP index
61 AliDCSValue *value = (AliDCSValue*)fPhaseB1.At(index);
62 if (!value) {
63 AliFatal(Form("Invalid index of the beam1 data point: %d (0 -> %d)",
64 index,fPhaseB1.GetEntries()));
65 return NULL;
66 }
67 return value;
68}
69
70//______________________________________________________________________________
71const AliDCSValue* AliLHCClockPhase::GetPhaseB2DP(Int_t index) const
72{
73 // Get the value of the phase
74 // The argument is the DP index
75 AliDCSValue *value = (AliDCSValue*)fPhaseB2.At(index);
76 if (!value) {
77 AliFatal(Form("Invalid index of the beam2 data point: %d (0 -> %d)",
78 index,fPhaseB2.GetEntries()));
79 return NULL;
80 }
81 return value;
82}
30d5d647 83
84//______________________________________________________________________________
85Float_t AliLHCClockPhase::GetMeanPhaseB1() const
86{
87 // Get mean beam1 phase shift.
88 // The mean is calculated using all the data-points.
89 Int_t n = 0;
90 Float_t phase = 0;
91 for(Int_t i = 0; i < fPhaseB1.GetEntries(); ++i) {
92 AliDCSValue *value = (AliDCSValue*)fPhaseB1.At(i);
93 if (value) {
94 phase += value->GetFloat();
95 ++n;
96 }
97 }
98 if (n == 0) {
99 AliError("No beam1 measurements found! Assuming 0 phase shift!");
100 return 0;
101 }
102 phase /= n;
103 return phase;
104}
105
106//______________________________________________________________________________
107Float_t AliLHCClockPhase::GetMeanPhaseB2() const
108{
109 // Get mean beam2 phase shift.
110 // The mean is calculated using all the data-points.
111 Int_t n = 0;
112 Float_t phase = 0;
113 for(Int_t i = 0; i < fPhaseB2.GetEntries(); ++i) {
114 AliDCSValue *value = (AliDCSValue*)fPhaseB2.At(i);
115 if (value) {
116 phase += value->GetFloat();
117 ++n;
118 }
119 }
120 if (n == 0) {
121 AliError("No beam2 measurements found! Assuming 0 phase shift!");
122 return 0;
123 }
124 phase /= n;
125 return phase;
126}
127
128//______________________________________________________________________________
129Float_t AliLHCClockPhase::GetMeanPhase() const
130{
131 // Get mean phase shift.
132 // Beam1 and beam2 phases are calculated using
133 // all beam1 and beam2 data-points. The two phases are
134 // then averaged in order to get the mean value.
135 return (GetMeanPhaseB1() + GetMeanPhaseB2())/2;
136}
137
138//______________________________________________________________________________
139Float_t AliLHCClockPhase::GetPhaseB1(UInt_t timestamp) const
140{
141 // Get beam1 phase shift using the
142 // event time-stamp as an argument.
143 // This methods loops over data-points and
144 // returns the value of the closest (in time)
145 // measurement found.
146 Long64_t deltat = 0xffffffff;
147 Float_t phase = 0;
148 for(Int_t i = 0; i < fPhaseB1.GetEntries(); ++i) {
149 AliDCSValue *value = (AliDCSValue*)fPhaseB1.At(i);
150 if (value) {
151 if (TMath::Abs((Long64_t)timestamp-(Long64_t)value->GetTimeStamp()) <= deltat) {
152 phase = value->GetFloat();
153 deltat = TMath::Abs((Long64_t)timestamp-(Long64_t)value->GetTimeStamp());
154 }
155 }
156 }
157 if (deltat == 0xffffffff) {
158 AliError(Form("Can't get the beam1 phase shift at time-stamp = %u",timestamp));
159 return 0;
160 }
161 return phase;
162}
163
164//______________________________________________________________________________
165Float_t AliLHCClockPhase::GetPhaseB2(UInt_t timestamp) const
166{
167 // Get beam2 phase shift using the
168 // event time-stamp as an argument.
169 // This methods loops over data-points and
170 // returns the value of the closest (in time)
171 // measurement found.
172 Long64_t deltat = 0xffffffff;
173 Float_t phase = 0;
174 for(Int_t i = 0; i < fPhaseB2.GetEntries(); ++i) {
175 AliDCSValue *value = (AliDCSValue*)fPhaseB2.At(i);
176 if (value) {
177 if (TMath::Abs((Long64_t)timestamp-(Long64_t)value->GetTimeStamp()) <= deltat) {
178 phase = value->GetFloat();
179 deltat = TMath::Abs((Long64_t)timestamp-(Long64_t)value->GetTimeStamp());
180 }
181 }
182 }
183 if (deltat == 0xffffffff) {
184 AliError(Form("Can't get the beam2 phase shift at time-stamp = %u",timestamp));
185 return 0;
186 }
187 return phase;
188}
189
190//______________________________________________________________________________
191Float_t AliLHCClockPhase::GetPhase(UInt_t timestamp) const
192{
193 // Get mean phase shift using the event time-stamp as
194 // an argument.
195 // Beam1 and beam2 phases are calculated using
196 // the closest beam1 and beam2 data-points. The two phases are
197 // then averaged in order to get the mean value.
198 return (GetPhaseB1(timestamp) + GetPhaseB2(timestamp))/2;
199}
200
201//______________________________________________________________________________
202void AliLHCClockPhase::Print( const Option_t* ) const
203{
204 // Print all the data-points for beam1 and beam2
205 // as well as the mean phases
206 cout << "AliLHCClockPhase object:" << endl;
207
208 cout << "Beam1:" << endl;
209 for(Int_t i = 0; i < fPhaseB1.GetEntries(); ++i) {
210 AliDCSValue *value = (AliDCSValue*)fPhaseB1.At(i);
211 if (value) cout << "TS=" << value->GetTimeStamp() << " " << value->GetFloat() << endl;
212 }
213 cout << "Beam1 mean phase=" << GetMeanPhaseB1() << " ns" << endl;
214
215 cout << "Beam2:" << endl;
216 for(Int_t i = 0; i < fPhaseB2.GetEntries(); ++i) {
217 AliDCSValue *value = (AliDCSValue*)fPhaseB2.At(i);
218 if (value) cout << "TS=" << value->GetTimeStamp() << " " << value->GetFloat() << endl;
219 }
220 cout << "Beam2 mean phase=" << GetMeanPhaseB2() << " ns" << endl;
221
222 cout << "Mean phase (beam1 & beam2) =" << GetMeanPhase() << " ns" << endl;
223}