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(startTime);
69 fEndTime = TTimeStamp(endTime);
72 //_____________________________________________________________________________
73 AliTPCSensorTempArray::AliTPCSensorTempArray(const char *fname) :
77 // AliTPCSensorTempArray constructor
79 fSensors = AliTPCSensorTemp::ReadList(fname);
80 fSensors->BypassStreamer(kFALSE);
84 //_____________________________________________________________________________
85 AliTPCSensorTempArray::AliTPCSensorTempArray(const AliTPCSensorTempArray &c):
89 // AliTPCSensorTempArray copy constructor
94 ///_____________________________________________________________________________
95 AliTPCSensorTempArray::~AliTPCSensorTempArray()
98 // AliTPCSensorTempArray destructor
102 //_____________________________________________________________________________
103 AliTPCSensorTempArray &AliTPCSensorTempArray::operator=(const AliTPCSensorTempArray &c)
106 // Assignment operator
110 new (this) AliTPCSensorTempArray(c);
111 fSensors = (TClonesArray*)c.fSensors->Clone();
117 //_____________________________________________________________________________
118 void AliTPCSensorTempArray::ReadSensors(const char *dbEntry)
121 // Read list of temperature sensors from text file
123 AliCDBEntry *entry = AliCDBManager::Instance()->Get(dbEntry);
124 TTree *tree = (TTree*) entry->GetObject();
125 fSensors = AliTPCSensorTemp::ReadTree(tree);
129 //_____________________________________________________________________________
130 AliTPCSensorTemp* AliTPCSensorTempArray::GetSensor(Int_t type, Int_t side, Int_t sector, Int_t num)
133 // Return sensor information for sensor specified by type, side, sector and num
135 Int_t nsensors = fSensors->GetEntries();
136 for (Int_t isensor=0; isensor<nsensors; isensor++) {
137 AliTPCSensorTemp *entry = (AliTPCSensorTemp*)fSensors->At(isensor);
138 if (entry->GetSide() == side &&
139 entry->GetType() == type &&
140 entry->GetSector() == sector &&
141 entry->GetNum() == num ) return entry;
145 //_____________________________________________________________________________
147 AliTPCSensorTemp* AliTPCSensorTempArray::GetSensor(Int_t IdDCS){
148 return dynamic_cast<AliTPCSensorTemp*>(AliDCSSensorArray::GetSensor(IdDCS));
150 //_____________________________________________________________________________
152 AliTPCSensorTemp* AliTPCSensorTempArray::GetSensor(Double_t x, Double_t y, Double_t z){
153 return dynamic_cast<AliTPCSensorTemp*>(AliDCSSensorArray::GetSensor(x,y,z));
155 //_____________________________________________________________________________
157 Double_t AliTPCSensorTempArray::GetTempGradientY(UInt_t timeSec, Int_t side){
159 // Extract Linear Vertical Temperature Gradient [K/cm] within the TPC on
162 // Values based on TemperatureSensors within the TPC (type: 3(TPC))
164 // FIXME: Also return residual-distribution, covariance Matrix
165 // or simply chi2 for validity check?
168 TLinearFitter fitter(3,"x0++x1++x2");
172 Int_t nsensors = fSensors->GetEntries();
173 for (Int_t isensor=0; isensor<nsensors; isensor++) { // loop over all sensors
174 AliTPCSensorTemp *entry = (AliTPCSensorTemp*)fSensors->At(isensor);
176 if (entry->GetType()==3 && entry->GetSide()==side) { // take SensorType:TPC
181 Double_t y = entry->GetValue(timeSec); // get temperature value
182 fitter.AddPoint(x,y,1); // add values to LinearFitter
188 fitter.GetParameters(param);
190 return param[2]; // return vertical (Y) tempGradient