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 **************************************************************************/
17 ///////////////////////////////////////////////////////////////////////////////
19 // TPC calibration class for parameters which saved per pad //
20 // Authors: Marian Ivanov and Haavard Helstrup //
22 ///////////////////////////////////////////////////////////////////////////////
24 #include "AliTPCSensorTempArray.h"
25 #include "TLinearFitter.h"
31 ClassImp(AliTPCSensorTempArray)
34 //_____________________________________________________________________________
35 AliTPCSensorTempArray::AliTPCSensorTempArray():AliDCSSensorArray()
38 // AliTPCSensorTempArray default constructor
42 //_____________________________________________________________________________
43 AliTPCSensorTempArray::AliTPCSensorTempArray(Int_t run) : AliDCSSensorArray()
46 // Read configuration from OCDB
51 AliCDBManager::Instance()->Get("TPC/Config/Temperature",run);
52 TTree *tree = (TTree*) entry->GetObject();
53 fSensors = AliTPCSensorTemp::ReadTree(tree);
54 fSensors->BypassStreamer(kFALSE);
57 //_____________________________________________________________________________
58 AliTPCSensorTempArray::AliTPCSensorTempArray(UInt_t startTime, UInt_t endTime,
59 TTree* confTree, const TString& amandaString)
63 // AliTPCSensorTempArray constructor for Shuttle preprocessor
64 // (confTree read from OCDB)
66 fSensors = AliTPCSensorTemp::ReadTree(confTree,amandaString);
67 fSensors->BypassStreamer(kFALSE);
68 fStartTime = TTimeStamp((time_t)startTime,0);
69 fEndTime = TTimeStamp((time_t)endTime,0);
72 //_____________________________________________________________________________
73 AliTPCSensorTempArray::AliTPCSensorTempArray(const char *fname,
74 const TString& amandaString) :
78 // AliTPCSensorTempArray constructor
80 fSensors = AliTPCSensorTemp::ReadList(fname,amandaString);
81 fSensors->BypassStreamer(kFALSE);
85 //_____________________________________________________________________________
86 AliTPCSensorTempArray::AliTPCSensorTempArray(const AliTPCSensorTempArray &c):
90 // AliTPCSensorTempArray copy constructor
95 ///_____________________________________________________________________________
96 AliTPCSensorTempArray::~AliTPCSensorTempArray()
99 // AliTPCSensorTempArray destructor
103 //_____________________________________________________________________________
104 AliTPCSensorTempArray &AliTPCSensorTempArray::operator=(const AliTPCSensorTempArray &c)
107 // Assignment operator
111 new (this) AliTPCSensorTempArray(c);
112 fSensors = (TClonesArray*)c.fSensors->Clone();
118 //_____________________________________________________________________________
119 void AliTPCSensorTempArray::ReadSensors(const char *dbEntry)
122 // Read list of temperature sensors from text file
124 AliCDBEntry *entry = AliCDBManager::Instance()->Get(dbEntry);
125 TTree *tree = (TTree*) entry->GetObject();
126 fSensors = AliTPCSensorTemp::ReadTree(tree);
130 //_____________________________________________________________________________
131 AliTPCSensorTemp* AliTPCSensorTempArray::GetSensor(Int_t type, Int_t side, Int_t sector, Int_t num)
134 // Return sensor information for sensor specified by type, side, sector and num
136 Int_t nsensors = fSensors->GetEntries();
137 for (Int_t isensor=0; isensor<nsensors; isensor++) {
138 AliTPCSensorTemp *entry = (AliTPCSensorTemp*)fSensors->At(isensor);
139 if (entry->GetSide() == side &&
140 entry->GetType() == type &&
141 entry->GetSector() == sector &&
142 entry->GetNum() == num ) return entry;
146 //_____________________________________________________________________________
148 AliTPCSensorTemp* AliTPCSensorTempArray::GetSensor(Int_t IdDCS){
149 return dynamic_cast<AliTPCSensorTemp*>(AliDCSSensorArray::GetSensor(IdDCS));
151 //_____________________________________________________________________________
153 AliTPCSensorTemp* AliTPCSensorTempArray::GetSensor(Double_t x, Double_t y, Double_t z){
154 return dynamic_cast<AliTPCSensorTemp*>(AliDCSSensorArray::GetSensor(x,y,z));
156 //_____________________________________________________________________________
158 Double_t AliTPCSensorTempArray::GetTempGradientY(UInt_t timeSec, Int_t side){
160 // Extract Linear Vertical Temperature Gradient [K/cm] within the TPC on
163 // Values based on TemperatureSensors within the TPC (type: 3(TPC))
165 // FIXME: Also return residual-distribution, covariance Matrix
166 // or simply chi2 for validity check?
169 TLinearFitter fitter(3,"x0++x1++x2");
173 Int_t nsensors = fSensors->GetEntries();
174 for (Int_t isensor=0; isensor<nsensors; isensor++) { // loop over all sensors
175 AliTPCSensorTemp *entry = (AliTPCSensorTemp*)fSensors->At(isensor);
177 if (entry->GetType()==3 && entry->GetSide()==side) { // take SensorType:TPC
182 Double_t y = entry->GetValue(timeSec); // get temperature value
183 fitter.AddPoint(x,y,1); // add values to LinearFitter
189 fitter.GetParameters(param);
191 return param[2]; // return vertical (Y) tempGradient