Added OutsideThreshold function. Corrections to TTimeStamp syntax. Set startTime...
[u/mrichter/AliRoot.git] / STEER / AliDCSSensor.cxx
CommitLineData
7264822f 1/**************************************************************************
2 * Copyright(c) 2006-07, 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
17////////////////////////////////////////////////////////////////////////////////
18// //
19// Class describing TPC temperature sensors (including pointers to graphs/fits//
20// Authors: Marian Ivanov, Haavard Helstrup and Martin Siska //
21// //
22////////////////////////////////////////////////////////////////////////////////
23
24
7264822f 25#include "AliDCSSensor.h"
e27404da 26ClassImp(AliDCSSensor)
8cb8848e 27
8cb8848e 28const Double_t kSecInHour = 3600.; // seconds in one hour
29
7264822f 30
31
32AliDCSSensor::AliDCSSensor():
33 fId(),
34 fIdDCS(0),
24cd438b 35 fStringID(),
7264822f 36 fStartTime(0),
8cb8848e 37 fEndTime(0),
7264822f 38 fGraph(0),
39 fFit(0),
40 fX(0),
41 fY(0),
42 fZ(0)
43{
44 //
45 // Standard constructor
46 //
47}
48
49AliDCSSensor::AliDCSSensor(const AliDCSSensor& source) :
50 TNamed(source),
51 fId(source.fId),
52 fIdDCS(source.fIdDCS),
24cd438b 53 fStringID(source.fStringID),
7264822f 54 fStartTime(source.fStartTime),
8cb8848e 55 fEndTime(source.fEndTime),
7264822f 56 fGraph(source.fGraph),
57 fFit(source.fFit),
58 fX(source.fX),
59 fY(source.fY),
60 fZ(source.fZ)
61//
62// Copy constructor
63//
64{ }
65
66AliDCSSensor& AliDCSSensor::operator=(const AliDCSSensor& source){
67//
68// assignment operator
69//
70 if (&source == this) return *this;
71 new (this) AliDCSSensor(source);
dc0184e7 72
73 return *this;
7264822f 74}
75
76//_____________________________________________________________________________
dc0184e7 77Double_t AliDCSSensor::GetValue(UInt_t timeSec)
7264822f 78{
79 //
80 // Get temperature value for actual sensor
81 // timeSec given as offset from start-of-run measured in seconds
82 //
67a165ed 83 Bool_t inside;
bf84b403 84 return Eval(TTimeStamp((time_t)(fStartTime+timeSec)),inside);
7264822f 85}
86//_____________________________________________________________________________
87Double_t AliDCSSensor::GetValue(TTimeStamp time)
88{
89 // Get temperature value for actual sensor
90 // time given as absolute TTimeStamp
91 //
67a165ed 92 Bool_t inside;
93 return Eval(time, inside);
8cb8848e 94}
95
96//_____________________________________________________________________________
97
67a165ed 98Double_t AliDCSSensor::Eval(const TTimeStamp& time, Bool_t inside) const
8cb8848e 99{
100 //
101 // Return temperature at given time
67a165ed 102 // If time < start of map return value at start of map, inside = false
103 // If time > end of map return value at end of map, inside = false
8cb8848e 104
105 UInt_t timeSec = time.GetSec();
106 UInt_t diff = timeSec-fStartTime;
67a165ed 107 inside = true;
8cb8848e 108
67a165ed 109 if ( timeSec < fStartTime ) {
110 inside=false;
111 diff=0;
112 }
113 if ( timeSec > fEndTime ) {
114 inside=false;
115 diff = fEndTime-fStartTime;
116 }
8cb8848e 117
118 Double_t timeHour = diff/kSecInHour;
119 if ( fFit ) {
120 return fFit->Eval(timeHour);
121 } else {
67a165ed 122 return -99;
8cb8848e 123 }
124}
125
126TGraph* AliDCSSensor::MakeGraph(Int_t nPoints) const
127{
128 //
129 // Make graph from start time to end time of DCS values
130 //
131
132 UInt_t stepTime = (fEndTime-fStartTime)/nPoints;
133
134 if ( !fFit ) return 0;
135
136 Double_t *x = new Double_t[nPoints+1];
137 Double_t *y = new Double_t[nPoints+1];
138 for (Int_t ip=0; ip<nPoints; ip++) {
139 x[ip] = fStartTime+ip*stepTime;
140 y[ip] = fFit->Eval(ip*stepTime/kSecInHour);
141 }
142
143 TGraph *graph = new TGraph(nPoints,x,y);
144 delete [] x;
145 delete [] y;
146
147 graph->GetXaxis()->SetTimeDisplay(1);
148 graph->GetXaxis()->SetLabelOffset(0.02);
149 graph->GetXaxis()->SetTimeFormat("#splitline{%d/%m}{%H:%M}");
82dfb6c8 150
8cb8848e 151 return graph;
7264822f 152}
153
24cd438b 154//_____________________________________________________________________________
155
156TClonesArray * AliDCSSensor::ReadTree(TTree* tree) {
157 //
158 // read values from ascii file
159 //
82dfb6c8 160
24cd438b 161 Int_t nentries = tree->GetEntries();
82dfb6c8 162
24cd438b 163 char stringId[100];
164 Int_t num=0;
165 Int_t idDCS=0;
166 Double_t x=0;
167 Double_t y=0;
168 Double_t z=0;
169
170 tree->SetBranchAddress("StringID",&stringId);
171 tree->SetBranchAddress("IdDCS",&idDCS);
172 tree->SetBranchAddress("Num",&num);
173 tree->SetBranchAddress("X",&x);
174 tree->SetBranchAddress("Y",&y);
175 tree->SetBranchAddress("Z",&z);
176
177 // firstSensor = (Int_t)tree->GetMinimum("ECha");
178 // lastSensor = (Int_t)tree->GetMaximum("ECha");
179
180 TClonesArray * array = new TClonesArray("AliDCSSensor",nentries);
181 printf ("nentries = %d\n",nentries);
182
183 for (Int_t isensor=0; isensor<nentries; isensor++){
184 AliDCSSensor * sens = new ((*array)[isensor])AliDCSSensor;
185 tree->GetEntry(isensor);
186 sens->SetId(isensor);
187 sens->SetIdDCS(idDCS);
188 sens->SetStringID(TString(stringId));
189 sens->SetX(x);
190 sens->SetY(y);
191 sens->SetZ(z);
192
193 }
194 return array;
195}
7264822f 196