STEER/AliDCSSensor.cxx

   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
  25 #include "AliDCSSensor.h"
  26 #include "TDatime.h"
  27 ClassImp(AliDCSSensor)
  28
  29 const Double_t kSecInHour = 3600.; // seconds in one hour
  30
  31
  32
  33 AliDCSSensor::AliDCSSensor():
  34   fId(),
  35   fIdDCS(0),
  36   fStringID(),
  37   fStartTime(0),
  38   fEndTime(0),
  39   fGraph(0),
  40   fFit(0),
  41   fX(0),
  42   fY(0),
  43   fZ(0)
  44 {
  45   //
  46   //  Standard constructor
  47   //
  48 }
  49
  50 AliDCSSensor::AliDCSSensor(const AliDCSSensor& source) :
  51    TNamed(source),
  52    fId(source.fId),
  53    fIdDCS(source.fIdDCS),
  54    fStringID(source.fStringID),
  55    fStartTime(source.fStartTime),
  56    fEndTime(source.fEndTime),
  57    fGraph(source.fGraph),
  58    fFit(source.fFit),
  59    fX(source.fX),
  60    fY(source.fY),
  61    fZ(source.fZ)
  62 //
  63 //  Copy constructor
  64 //
  65 { }
  66
  67 AliDCSSensor& AliDCSSensor::operator=(const AliDCSSensor& source){
  68 //
  69 // assignment operator
  70 //
  71   if (&source == this) return *this;
  72   new (this) AliDCSSensor(source);
  73
  74   return *this;
  75 }
  76
  77 //_____________________________________________________________________________
  78 Double_t AliDCSSensor::GetValue(UInt_t timeSec)
  79 {
  80  //
  81  // Get temperature value for actual sensor
  82  //  timeSec given as offset from start-of-map measured in seconds
  83  //  *NOTE* In the current TPC setup, start-of-map is defined as the
  84  //         first measured point for each sensor. This will be different
  85  //         for each sensor in the array. If you want to get a value at the
  86  //         same absolute time, use AliDCSSensor::GetValue(TTimeStamp time)
  87  //         or AliDCSSensorArray::GetValue (UInt_t timeSec, Int_t sensor)
  88  //         which measure offsets with respect to the (global) start-of-run
  89  //
  90  Bool_t inside=kTRUE;
  91  return Eval(TTimeStamp((time_t)(fStartTime+timeSec),0),inside);
  92 }
  93 //_____________________________________________________________________________
  94 Double_t AliDCSSensor::GetValue(TTimeStamp time)
  95 {
  96  // Get temperature value for actual sensor
  97  //  time given as absolute TTimeStamp
  98  //
  99  Bool_t inside=kTRUE;
 100  return Eval(time, inside);
 101 }
 102
 103 //_____________________________________________________________________________
 104
 105 Double_t AliDCSSensor::Eval(const TTimeStamp& time, Bool_t inside) const
 106 {
 107   //
 108   // Return temperature at given time
 109   //  If time < start of map  return value at start of map, inside = false
 110   //  If time > end of map    return value at end of map, inside = false
 111
 112   UInt_t timeSec = time.GetSec();
 113   UInt_t diff = timeSec-fStartTime;
 114   inside = true;
 115
 116   if ( timeSec < fStartTime ) {
 117      inside=false;
 118      diff=0;
 119   }
 120   if ( timeSec > fEndTime ) {
 121      inside=false;
 122      diff = fEndTime-fStartTime;
 123   }
 124
 125   Double_t timeHour = diff/kSecInHour;
 126   if ( fFit ) {
 127      return fFit->Eval(timeHour);
 128   } else {
 129      if ( fGraph ) {
 130        return EvalGraph(timeHour);
 131      } else {
 132        return -99;
 133      }
 134   }
 135 }
 136 //_____________________________________________________________________________
 137 Double_t AliDCSSensor::EvalGraph(const Double_t& timeHour) const
 138 {
 139   //
 140   // Extract last value in graph observed before time given by timeHour
 141   //
 142
 143   // return -99 if point specified is before beginning of graph
 144   Double_t x=0; Double_t y=0;
 145   fGraph->GetPoint(0,x,y);
 146   if ( timeHour < x ) return -99;
 147
 148   // return previous point when first time > timeHour is observed
 149
 150   Int_t npoints = fGraph->GetN();
 151   for (Int_t i=1; i<npoints; i++) {
 152      fGraph->GetPoint(i,x,y);
 153      if ( timeHour < x ) {
 154        fGraph->GetPoint(i-1,x,y);
 155        return y;
 156      }
 157   }
 158
 159   // return last point if all times are < timeHour
 160   return y;
 161 }
 162
 163
 164 //_____________________________________________________________________________
 165 TGraph* AliDCSSensor::MakeGraph(Int_t nPoints, Bool_t debug) const
 166 {
 167   //
 168   // Make graph from start time to end time of DCS values
 169   //
 170
 171
 172
 173   UInt_t stepTime = (fEndTime-fStartTime)/nPoints;
 174
 175   if (debug==kTRUE) {
 176      printf ("Start time %d, End time %d, step time %d\n",
 177      fStartTime,fEndTime,stepTime);
 178      TTimeStamp t((time_t)fStartTime,0); t.Print();
 179      TTimeStamp t2((time_t)fEndTime,0); t2.Print();
 180   }
 181
 182   if ( !fFit ) return 0;
 183
 184   Double_t *x = new Double_t[nPoints+1];
 185   Double_t *y = new Double_t[nPoints+1];
 186   for (Int_t ip=0; ip<nPoints; ip++) {
 187     x[ip] = (time_t)(fStartTime+ip*stepTime);
 188     y[ip] = fFit->Eval(ip*stepTime/kSecInHour);
 189     if (debug==kTRUE) {
 190      TTimeStamp t3((time_t)x[ip],0);
 191      printf ("x=%f, y=%f  ",x[ip],y[ip]);
 192      t3.Print();
 193     }
 194   }
 195
 196   TGraph *graph = new TGraph(nPoints,x,y);
 197   delete [] x;
 198   delete [] y;
 199
 200   graph->GetXaxis()->SetTimeDisplay(1);
 201   graph->GetXaxis()->SetLabelOffset(0.02);
 202   graph->GetXaxis()->SetTimeFormat("#splitline{%d/%m}{%H:%M}");
 203
 204   return graph;
 205 }
 206
 207 //_____________________________________________________________________________
 208
 209 TClonesArray * AliDCSSensor::ReadTree(TTree* tree) {
 210   //
 211   // read values from ascii file
 212   //
 213
 214   Int_t nentries = tree->GetEntries();
 215
 216   char stringId[100];
 217   Int_t num=0;
 218   Int_t idDCS=0;
 219   Double_t x=0;
 220   Double_t y=0;
 221   Double_t z=0;
 222
 223   tree->SetBranchAddress("StringID",&stringId);
 224   tree->SetBranchAddress("IdDCS",&idDCS);
 225   tree->SetBranchAddress("Num",&num);
 226   tree->SetBranchAddress("X",&x);
 227   tree->SetBranchAddress("Y",&y);
 228   tree->SetBranchAddress("Z",&z);
 229
 230   // firstSensor = (Int_t)tree->GetMinimum("ECha");
 231   // lastSensor = (Int_t)tree->GetMaximum("ECha");
 232
 233   TClonesArray * array = new TClonesArray("AliDCSSensor",nentries);
 234    printf ("nentries = %d\n",nentries);
 235
 236   for (Int_t isensor=0; isensor<nentries; isensor++){
 237     AliDCSSensor * sens = new ((*array)[isensor])AliDCSSensor;
 238     tree->GetEntry(isensor);
 239     sens->SetId(isensor);
 240     sens->SetIdDCS(idDCS);
 241     sens->SetStringID(TString(stringId));
 242     sens->SetX(x);
 243     sens->SetY(y);
 244     sens->SetZ(z);
 245
 246   }
 247   return array;
 248 }
 249