0th version of time-dependent basic calibration
authordsilverm <dsilverm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 6 Jul 2008 21:10:25 +0000 (21:10 +0000)
committerdsilverm <dsilverm@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 6 Jul 2008 21:10:25 +0000 (21:10 +0000)
EMCAL/AliEMCALCalibTimeDep.cxx [new file with mode: 0644]
EMCAL/AliEMCALCalibTimeDep.h [new file with mode: 0644]

diff --git a/EMCAL/AliEMCALCalibTimeDep.cxx b/EMCAL/AliEMCALCalibTimeDep.cxx
new file mode 100644 (file)
index 0000000..8ecd8a8
--- /dev/null
@@ -0,0 +1,277 @@
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+/* $Id: AliEMCALCalibTimeDep.cxx $ */
+
+//_________________________________________________________________________
+///*-- Author: 
+///////////////////////////////////////////////////////////////////////////////
+//                                                                           //
+// class for EMCAL time-dep calibration                                      //
+//                                                                           //
+///////////////////////////////////////////////////////////////////////////////
+
+#include <iostream>
+#include <TGraphSmooth.h>
+#include "AliLog.h"
+#include "AliCDBEntry.h"
+#include "AliCDBManager.h"
+#include "AliEMCALSensorTempArray.h"
+#include "AliEMCALCalibTimeDep.h"
+
+/* first a bunch of constants.. */
+const double fkSecToHour = 1.0/3600.0; // conversion factor from seconds to hours
+
+// some global variables for APD handling; assume the same for all, at least to start with
+const double fkTempSlope = 0.017; // 1.7% per deg. C, seems about right for all APDs, from studies at Catania, and by Rachid
+const double fkNormTemp = 20; // let's say that 20 degrees C is normal as default
+
+using namespace std;
+
+ClassImp(AliEMCALCalibTimeDep)
+
+//________________________________________________________________
+AliEMCALCalibTimeDep::AliEMCALCalibTimeDep() :
+  fRun(0),
+  fStartTime(0),
+  fEndTime(0),
+  fMinTemp(0),
+  fMaxTemp(0),
+  fMinTime(0),
+  fMaxTime(0),
+  fRefTemp(fkNormTemp), 
+  fTempArray(NULL)
+{
+  // Constructor
+}
+
+//________________________________________________________________
+AliEMCALCalibTimeDep::AliEMCALCalibTimeDep(const AliEMCALCalibTimeDep& calibt) :
+  fRun(calibt.GetRunNumber()),
+  fStartTime(calibt.GetStartTime()),
+  fEndTime(calibt.GetEndTime()),
+  fMinTemp(calibt.GetMinTemp()),
+  fMaxTemp(calibt.GetMaxTemp()),
+  fMinTime(calibt.GetMinTime()),
+  fMaxTime(calibt.GetMaxTime()),
+  fRefTemp(calibt.GetRefTemp()),
+  fTempArray(calibt.GetTempArray())
+{
+  // copy constructor
+}
+
+
+//________________________________________________________________
+AliEMCALCalibTimeDep &AliEMCALCalibTimeDep::operator =(const AliEMCALCalibTimeDep& calibt)
+{
+  // assignment operator; use copy ctor
+  if (&calibt == this) return *this;
+
+  new (this) AliEMCALCalibTimeDep(calibt);
+  return *this;
+}
+
+//________________________________________________________________
+AliEMCALCalibTimeDep::~AliEMCALCalibTimeDep()
+{
+  // Destructor
+}
+
+//________________________________________________________________
+void  AliEMCALCalibTimeDep::Reset() 
+{
+  // clear variables to default
+  fRun = 0;
+  fStartTime = 0;
+  fEndTime = 0;
+  fMinTemp = 0;
+  fMaxTemp = 0;
+  fMinTime = 0;
+  fMaxTime = 0;
+  fRefTemp = fkNormTemp;
+  fTempArray = NULL;
+  return;
+}
+
+//________________________________________________________________
+void  AliEMCALCalibTimeDep::Print() const
+{
+  // print some info
+  cout << endl << " AliEMCALCalibTimeDep::Print() " << endl;
+  // basic variables, all 'publicly available' also
+  cout << " VARIABLE DUMP: " << endl
+       << " GetStartTime() " << GetStartTime() << endl
+       << " GetEndTime() " << GetEndTime() << endl
+       << " GetMinTemp() " << GetMinTemp() << endl
+       << " GetMaxTemp() " << GetMaxTemp() << endl
+       << " GetRefTemp() " << GetRefTemp() << endl;
+  // run ranges
+  cout << " RUN INFO: " << endl
+       << " length (in hours) " << GetLengthOfRunInHours() << endl
+       << " range of temperature measurements (in hours) " << GetRangeOfTempMeasureInHours()
+       << " (in deg. C) " << GetRangeOfTempMeasureInDegrees()
+       << endl;
+  // range in correction values
+  double corrAtMinTemp = GetCorrection( fMinTemp );
+  double corrAtMaxTemp = GetCorrection( fMaxTemp );
+  double corrMaxMinDiff = 100*(corrAtMinTemp - corrAtMaxTemp);
+  cout << " CORRECTION INFO : " << endl
+       << " corrAtMinTemp " << corrAtMinTemp << endl
+       << " corrAtMaxTemp " << corrAtMaxTemp << endl
+       << " corrMaxMinDiff (~%) [=(corrAtMin - corrAtMax)*100] " 
+       << corrMaxMinDiff << endl;
+
+  return;
+}
+//________________________________________________________________ 
+double AliEMCALCalibTimeDep::GetLengthOfRunInHours() const
+{
+  return (fEndTime - fStartTime)*fkSecToHour;
+}
+//________________________________________________________________ 
+double AliEMCALCalibTimeDep::GetRangeOfTempMeasureInHours() const
+{
+  return (fMaxTime - fMinTime)*fkSecToHour;
+}
+//________________________________________________________________ 
+double AliEMCALCalibTimeDep::GetRangeOfTempMeasureInDegrees() const
+{
+  return (fMaxTemp - fMinTemp);
+}
+
+//________________________________________________________________
+void AliEMCALCalibTimeDep::Initialize(Int_t run, 
+                                     UInt_t startTime, UInt_t endTime)
+{
+  Reset(); // start fresh
+
+  fRun = run;
+  fStartTime = startTime;
+  fEndTime = endTime;
+  
+  // collect the needed information
+  GetTemperatureInfo(); // temperature readings during the run
+
+  return;
+}
+
+//________________________________________________________________
+double AliEMCALCalibTimeDep::GetTemperature(int secSinceRunStart) const
+{// return estimate for all SuperModules and sensors, that had data 
+
+  // first convert from seconds to hours..
+  double timeHour = secSinceRunStart * fkSecToHour;
+
+  //  return fTempSpline->Eval(timeHour);
+  return timeHour; // DStmp - FIXME - just return time for now
+}
+
+//________________________________________________________________
+double AliEMCALCalibTimeDep::GetTemperatureSM(int imod, int secSinceRunStart) const
+{// return estimate for this one SuperModule, if it had data 
+
+  // first convert from seconds to hours..
+  double timeHour = secSinceRunStart * fkSecToHour;
+
+  //  return fTempSpline->Eval(timeHour);
+  return timeHour; // DStmp - FIXME - just return time for now
+}
+
+//________________________________________________________________
+double AliEMCALCalibTimeDep::GetTemperatureSMSensor(int imod, int isens, int secSinceRunStart) const
+{// return estimate for this one SuperModule and sensor, if it had data 
+  // first convert from seconds to hours..
+  double timeHour = secSinceRunStart * fkSecToHour;
+  
+  //  return fTempSpline->Eval(timeHour);
+  return timeHour; // DStmp - FIXME - just return time for now
+}
+
+//________________________________________________________________
+double AliEMCALCalibTimeDep::GetCorrection(double temperature) const
+{ // gain decreases with increasing temperature
+  double gain = (1.0 - (temperature-fRefTemp)*fkTempSlope);  
+  // i.e. multiplicative correction to ADC increases with increasing temperature
+  return 1.0/gain;
+}
+
+/* Next comes the method that does the work in picking up all the needed info..*/
+//________________________________________________________________
+void AliEMCALCalibTimeDep::GetTemperatureInfo() 
+{
+  // pick up Preprocessor output, based on fRun (most recent version)
+  AliCDBEntry* entry = AliCDBManager::Instance()->Get("EMCAL/Calib/Temperature", fRun);
+  if (entry) {
+    fTempArray = (AliEMCALSensorTempArray *) entry->GetObject();
+  }
+
+  if (fTempArray) { // DStmp - FIXME - not implemented what should be done here yet
+    AliInfo( Form("NumSensors %d - IdDCS: first %d last %d",
+                 fTempArray->NumSensors(),
+                 fTempArray->GetFirstIdDCS(), fTempArray->GetLastIdDCS() ) );
+
+    /* // below is just examples on what could be done to access the data..
+  AliEMCALSensorTemp *o = arr->GetSensor(0);
+  o->Print();
+  cout << " side " << o->GetSide()
+       << " sector " << o->GetSector()
+       << " num " << o->GetNum()
+       << " startTime " << o->GetStartTime()
+       << " endTime " << o->GetEndTime()
+       << endl;
+
+  AliSplineFit *f = o->GetFit();
+  int np = 0;
+
+  if (f) {
+    f->SplineFit(0);
+
+    np = f->GetKnots();
+    cout << " np " << np << endl;
+    Double_t *x = f->GetX();
+    Double_t *y0 = f->GetY0();
+    Double_t *y1 = f->GetY1();
+    for (int i=0; i<np; i++) {
+      cout << " i " << i
+           << " x " << x[i]
+           << " y0 " << y0[i]
+           << " y1 " << y1[i]
+           << endl;
+      //      g->SetPoint(i, x[i], y0[i]);
+    }
+    double start = x[0];
+    double stop = x[np-1];
+
+    TGraph *g = f->MakeGraph(start, stop, np);
+    f->MakeSmooth(g, ratio, "normal");
+    np = f->GetKnots();
+    cout << " np " << np << endl;
+    TGraph *gSmooth = f->MakeGraph(start, stop, np);
+
+  // compare the raw and smoothed versions and select one..
+    // plot in 2 canvases/pads - CONTINUE HERE
+
+    g->Draw("ALP");
+    gSmooth->Draw("ALP");
+    f->Draw();
+    }
+
+     */
+
+  }
+  
+  return;
+}
diff --git a/EMCAL/AliEMCALCalibTimeDep.h b/EMCAL/AliEMCALCalibTimeDep.h
new file mode 100644 (file)
index 0000000..4daae97
--- /dev/null
@@ -0,0 +1,74 @@
+#ifndef ALIEMCALCALIBTIMEDEP_H
+#define ALIEMCALCALIBTIMEDEP_H
+
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/* $Id: AliEMCALCalibTimeDep.h $ */
+
+/////////////////////////////////////////////////
+// class for EMCAL time-dependence calibration //
+/////////////////////////////////////////////////
+
+#include "TObject.h"
+
+class AliEMCALSensorTempArray;
+
+class AliEMCALCalibTimeDep : public TObject {
+
+ public:
+  AliEMCALCalibTimeDep(); //! ctor
+  AliEMCALCalibTimeDep(const AliEMCALCalibTimeDep &calibt); //! copy ctor
+  AliEMCALCalibTimeDep& operator= (const AliEMCALCalibTimeDep &calibt); //! 
+  virtual ~AliEMCALCalibTimeDep(); //! dtor
+  virtual void Print() const; 
+
+  void Initialize(Int_t run, UInt_t startTime, UInt_t endTime);//!
+
+  // access pointer to the basic temperature object: AliEMCALSensorTempArray 
+  AliEMCALSensorTempArray  * GetTempArray() const { return fTempArray; } //
+
+  // simple getters
+  int GetRunNumber() const { return fRun; } //
+  int GetStartTime() const { return fStartTime; } // 
+  int GetEndTime() const { return fEndTime; } // 
+  double GetLengthOfRunInHours() const; // 
+
+  // range of temperature readings/values during the run 
+  double GetMinTemp() const { return fMinTemp; } // 
+  double GetMaxTemp() const { return fMaxTemp; } // 
+  int GetMinTime() const { return fMinTime; } // 
+  int GetMaxTime() const { return fMaxTime; } // 
+  double GetRangeOfTempMeasureInHours() const; // 
+  double GetRangeOfTempMeasureInDegrees() const; // 
+
+  // reference temperature, that we normalize to
+  double GetRefTemp() const { return fRefTemp; } //
+  void SetRefTemp(double refTemp) { fRefTemp = refTemp; } //
+
+  // basic calibration info
+  double GetTemperature(int secSinceStart) const; // for all sensors, all SuperModules
+  double GetTemperatureSM(int imod, int secSinceStart) const; // for all sensors, in a SuperModule
+  double GetTemperatureSMSensor(int imod, int isens, int secSinceStart) const; // for a sensor, in a SuperModule
+  double GetCorrection(double temperature) const; //
+
+ private:
+
+  void Reset();
+  void GetTemperatureInfo(); // pick up Preprocessor output
+  //
+  int fRun;
+  int fStartTime;
+  int fEndTime;
+  double fMinTemp;
+  double fMaxTemp;
+  int fMinTime;
+  int fMaxTime;
+  double fRefTemp;
+  //
+  AliEMCALSensorTempArray  *fTempArray;     // CDB class for temperature sensors
+  //
+  ClassDef(AliEMCALCalibTimeDep,1)    // EMCAL Calibration data
+};
+
+#endif