(Haavard)

[u/mrichter/AliRoot.git] / STEER / AliDCSSensor.cxx

diff --git a/STEER/AliDCSSensor.cxx b/STEER/AliDCSSensor.cxx
index 7f51c247173a279e305e13921288b9fccb8bd4be..bf194802968f656b12aa68aa25e7ad11ac0033b1 100644 (file)
--- a/STEER/AliDCSSensor.cxx
+++ b/STEER/AliDCSSensor.cxx
@@ -36,13 +36,19 @@
 
 
 #include "AliDCSSensor.h"
-ClassImp(AliDCSSensor)
+ClassImp(AliDCSSensor);
+
+const Double_t kSmall = -9e99;     // invalid small value
+const Double_t kLarge = 9e99;      // invalid large value
+const Double_t kSecInHour = 3600.; // seconds in one hour
+
 
 
 AliDCSSensor::AliDCSSensor():
   fId(),
   fIdDCS(0),
   fStartTime(0),
+  fEndTime(0),
   fGraph(0),
   fFit(0),
   fX(0),
@@ -59,6 +65,7 @@ AliDCSSensor::AliDCSSensor(const AliDCSSensor& source) :
    fId(source.fId),
    fIdDCS(source.fIdDCS),
    fStartTime(source.fStartTime),
+   fEndTime(source.fEndTime),
    fGraph(source.fGraph),
    fFit(source.fFit),
    fX(source.fX),
@@ -86,8 +93,12 @@ Double_t AliDCSSensor::GetValue(UInt_t timeSec)
  // Get temperature value for actual sensor
  //  timeSec given as offset from start-of-run measured in seconds
  //
- Double_t timeHrs = timeSec/3600.0;
- return  fFit->Eval(timeHrs,0);
+ Double_t timeHrs = timeSec/kSecInHour;
+ if (fFit) {
+  return  fFit->Eval(timeHrs,0);
+ } else { 
+  return kSmall;
+ }
 }
 //_____________________________________________________________________________
 Double_t AliDCSSensor::GetValue(TTimeStamp time) 
@@ -95,8 +106,63 @@ Double_t AliDCSSensor::GetValue(TTimeStamp time)
  // Get temperature value for actual sensor
  //  time given as absolute TTimeStamp
  //
- Double_t timeHrs = (time.GetSec() - fStartTime)/3600.0;
- return fFit->Eval(timeHrs,0);
+ Double_t timeHrs = (time.GetSec() - fStartTime)/kSecInHour;
+ if (fFit) {
+  return  fFit->Eval(timeHrs,0);
+ } else { 
+  return kSmall;
+ }
+}
+
+//_____________________________________________________________________________
+
+Double_t AliDCSSensor::Eval(const TTimeStamp& time) const
+{
+  // 
+  // Return temperature at given time
+  //  If time < start of map  return kSmall
+  //  If time > end of map    return kLarge
+  
+  UInt_t timeSec = time.GetSec();
+  UInt_t diff = timeSec-fStartTime;
+  
+  if ( timeSec < fStartTime ) return kSmall;
+  if ( timeSec > fEndTime ) return kLarge;
+ 
+  Double_t timeHour = diff/kSecInHour;
+  if ( fFit ) {
+     return fFit->Eval(timeHour); 
+  } else {
+     return kSmall;
+  }
+}
+
+TGraph* AliDCSSensor::MakeGraph(Int_t nPoints) const
+{
+  //
+  // Make graph from start time to end time of DCS values 
+  //
+
+  UInt_t stepTime = (fEndTime-fStartTime)/nPoints;
+  
+  if ( !fFit ) return 0;
+
+  Double_t *x = new Double_t[nPoints+1];
+  Double_t *y = new Double_t[nPoints+1];
+  for (Int_t ip=0; ip<nPoints; ip++) {
+    x[ip] = fStartTime+ip*stepTime;
+    y[ip] = fFit->Eval(ip*stepTime/kSecInHour);
+  }
+  
+  TGraph *graph = new TGraph(nPoints,x,y);
+  delete [] x;
+  delete [] y;
+  
+  graph->GetXaxis()->SetTimeDisplay(1);
+  graph->GetXaxis()->SetLabelOffset(0.02);
+  graph->GetXaxis()->SetTimeFormat("#splitline{%d/%m}{%H:%M}");
+  
+  return graph;
 }