memory leak fixed
[u/mrichter/AliRoot.git] / T0 / AliT0Parameters.cxx
index a752429..4d82cd6 100644 (file)
 #include "AliT0CalibData.h"   
 #include "AliT0CalibWalk.h"   
 #include "AliT0CalibTimeEq.h"   
-#include "AliT0LookUpValue.h"
+#include "AliT0CalibLatency.h"   
 #include "AliT0LookUpKey.h"
+#include "AliT0LookUpValue.h"
 #include <AliCDBManager.h>        
 #include <AliCDBEntry.h>          
+#include <AliCDBStorage.h>  
 #include <TMath.h>
+#include <TSystem.h>
+//#include <Riostream.h>
 #include <TGeoManager.h>
 #include <TGeoPhysicalNode.h>
+#include <TGeoMatrix.h>
 #include <AliGeomManager.h>
 
 AliT0CalibTimeEq* AliT0Parameters::fgCalibData = 0;
 AliT0CalibData* AliT0Parameters::fgLookUp = 0;
 AliT0CalibWalk* AliT0Parameters::fgSlewCorr =0;
+AliT0CalibLatency *AliT0Parameters::fgLatency=0;
 //====================================================================
 ClassImp(AliT0Parameters)
 #if 0
@@ -69,14 +75,23 @@ AliT0Parameters::AliT0Parameters()
    fAmpLEDRec(0), 
    fPMTeff(),
    fWalk(0),
+   fQTC(0),
+   fAmpLED(0),
    fTimeDelayCFD(0), 
  //  fTimeV0(0), 
    fTimeDelayTVD(0),
    fMeanT0(512),
+   fMeanVertex(0),
+   fLatencyHPTDC(0),
+   fLatencyL1(0),
+   fLatencyL1A(0),
+   fLatencyL1C(0),
    fLookUp(0),
    fNumberOfTRMs(2),
-   fCalibentry(), fLookUpentry(),fSlewCorr()
-
+   fCalibentry(), 
+   fLookUpentry(),
+   fSlewCorr(),
+   fLatency()
   
 {
   // Default constructor 
@@ -133,7 +148,18 @@ AliT0Parameters::Init()
     fIsInit = kFALSE;
     return;
   }
-  fIsInit = kTRUE;
+  //latency
+  
+ fLatency  = stor->Get("T0/Calib/Latency");
+  if (fLatency){
+    fgLatency  = (AliT0CalibLatency*)fLatency->GetObject();
+  }
+  else {
+     AliWarning(" !!! no latency  in CDB "); 
+    return;
+  }
+  
+fIsInit = kTRUE;
 }
 
 
@@ -145,33 +171,46 @@ void AliT0Parameters::InitIfOnline()
 // for switching to this one should write
   // AliT0RawReader myrawreader(rawReader);
 //     myrawreader.SetOnlineMode(kTRUE);
- if (fIsInit) return;
-   //standart configuration (used for simulation)
-   //Int_t trm=0; Int_t tdc=0; Int_t chain=0; Int_t channel=0;
+  
+  if (fIsInit) return;
+  //standart configuration (used for simulation)
+  //Int_t trm=0; Int_t tdc=0; Int_t chain=0; Int_t channel=0;
   // configuration for test Jun07.
-   fgLookUp = new AliT0CalibData("T0");
-
+  fgLookUp = new AliT0CalibData("T0");
+  
   fNumberOfTRMs = 1;
- fgLookUp-> SetNumberOfTRMs(fNumberOfTRMs);
+  fgLookUp-> SetNumberOfTRMs(fNumberOfTRMs);
   Int_t trm=7; Int_t tdc=0; Int_t chain=0; Int_t channel=0;
   for (Int_t ik=0; ik<105; ik++)
-        {
-         AliT0LookUpKey * lookkey= new AliT0LookUpKey();
-         AliT0LookUpValue * lookvalue= new AliT0LookUpValue();
-
-
-          lookvalue->SetTRM(trm);
-          lookvalue->SetTDC(tdc);
-          lookvalue->SetChain(chain);
-          lookvalue->SetChannel(channel);
-          lookkey->SetKey(ik);
-         if (channel<6) channel +=2;
-         else {channel = 0; tdc++;}
-         if(ik==57) { tdc=0; channel=0; chain = 1;}
-
-         fgLookUp->GetMapLookup()->Add((TObject*)lookvalue,(TObject*)lookkey); 
-       }
+    {
+      AliT0LookUpKey * lookkey= new AliT0LookUpKey();
+      AliT0LookUpValue * lookvalue= new AliT0LookUpValue();
+      lookvalue->SetTRM(trm);
+      lookvalue->SetTDC(tdc);
+      lookvalue->SetChain(chain);
+      lookvalue->SetChannel(channel);
+      lookkey->SetKey(ik);
+      fgLookUp->GetMapLookup()->Add((TObject*)lookvalue,(TObject*)lookkey);    
+      if (channel<6) channel +=2;
+      else {channel = 0; tdc++;}
+      if(ik==56) { tdc=0; channel=0; chain = 1;}
+    }
+  //2 recent open channels
+  trm=7; tdc=14; channel=2; chain=0; 
+  for (Int_t ik=105; ik<107; ik++)
+    {
+      AliT0LookUpKey * lookkey= new AliT0LookUpKey();
+      AliT0LookUpValue * lookvalue= new AliT0LookUpValue();
+      
+      lookvalue->SetTRM(trm);
+      lookvalue->SetTDC(tdc);
+      lookvalue->SetChain(chain);
+      lookvalue->SetChannel(channel);
+      lookkey->SetKey(ik);
+      fgLookUp->GetMapLookup()->Add((TObject*)lookvalue,(TObject*)lookkey);    
+      if (channel<6) channel +=2;
+      else {channel = 0; tdc++;}
+    }
   
   fIsInit=kTRUE;
 }
@@ -190,10 +229,70 @@ AliT0Parameters::GetTimeDelayCFD(Int_t ipmt)
   return fgCalibData->GetTimeEq(ipmt);
 }
 //__________________________________________________________________
+Float_t
+AliT0Parameters::GetCFD(Int_t ipmt) 
+  {
+  // return  CFD channel
+   
+    return fgCalibData->GetCFDvalue(ipmt,0);
+}
+
+//__________________________________________________________________
+Float_t
+AliT0Parameters::GetLatencyHPTDC() 
+  {
+  // return LatencyHPTDC for CFD channel
+  if (!fLatency) 
+    {
+      fLatencyHPTDC=9000.;
+      return fLatencyHPTDC;
+    }
+   
+  return fgLatency->GetLatencyHPTDC();
+}
+//__________________________________________________________________
+Float_t
+AliT0Parameters::GetLatencyL1() 
+  {
+  // return time delay for CFD channel
+   
+  return fgLatency->GetLatencyL1();
+}
+
+//__________________________________________________________________
+Float_t
+AliT0Parameters::GetLatencyL1A() 
+  {
+  // return time delay for CFD channel
+   
+  return fgLatency->GetLatencyL1A();
+}
+
+//__________________________________________________________________
+Float_t
+AliT0Parameters::GetLatencyL1C() 
+  {
+  // return time delay for CFD channel
+   
+  return fgLatency->GetLatencyL1C();
+}
+//__________________________________________________________________
+
+Float_t
+AliT0Parameters:: GetMeanVertex()
+{ 
+  if (!fCalibentry) 
+    {
+      fMeanVertex=0;
+      return fMeanVertex;
+    }
+   
+  return fgCalibData->GetMeanVertex();
+}
+//__________________________________________________________________
 
 TGraph *AliT0Parameters::GetAmpLEDRec(Int_t ipmt) const
 {
- // get walk correction by LED-CFD vs CFD graph
    if (!fSlewCorr) {
      AliError("No slewing correction is available!");
      return  (TGraph*)fAmpLEDRec.At(ipmt); 
@@ -205,7 +304,6 @@ TGraph *AliT0Parameters::GetAmpLEDRec(Int_t ipmt) const
 
 TGraph *AliT0Parameters::GetWalk(Int_t ipmt) const
 {
-  // get walk correction by QTC vs CFD
   if (!fSlewCorr) {
     AliError("No walk correction is available!");
     return  (TGraph*)fWalk.At(ipmt); 
@@ -215,21 +313,31 @@ TGraph *AliT0Parameters::GetWalk(Int_t ipmt) const
 
 //__________________________________________________________________
 
-Float_t AliT0Parameters::GetWalkVal(Int_t ipmt, Float_t mv) const
+TGraph *AliT0Parameters::GetQTC(Int_t ipmt) const
 {
- // get walk correction by QTC vs CFD
-   if (!fSlewCorr) {
-    return ((TGraph*)fWalk.At(ipmt))->Eval(mv); 
+  if (!fSlewCorr) {
+    AliError("No walk correction is available!");
+    //    return  (TGraph*)fQTC.At(ipmt); 
+   return  0; 
   } 
-  return fgSlewCorr -> GetWalkVal(ipmt, mv) ;
+  return fgSlewCorr -> GetQTC(ipmt) ;
 }
 
+//__________________________________________________________________
+TGraph *AliT0Parameters::GetAmpLED(Int_t ipmt) const
+{
+  if (!fSlewCorr) {
+    AliError("No walk correction is available!");
+    //    return  (TGraph*)fQTC.At(ipmt); 
+   return  0; 
+  } 
+  return fgSlewCorr -> GetAmpLED(ipmt) ;
+}
 
 //__________________________________________________________________
 void 
 AliT0Parameters::SetPMTeff(Int_t ipmt)
 {
-  // PMT quantum efficiency
   Float_t lambda[50];
   Float_t eff[50 ] = {0,        0,       0.23619,  0.202909, 0.177913, 
                    0.175667, 0.17856, 0.190769, 0.206667, 0.230286,
@@ -251,7 +359,6 @@ AliT0Parameters::SetPMTeff(Int_t ipmt)
 Int_t 
 AliT0Parameters::GetChannel(Int_t trm,  Int_t tdc, Int_t chain, Int_t channel)
 {
-    // get logical number of channel according physical address
 
   if (fgLookUp) {
     AliT0LookUpValue key(trm,tdc,chain,channel);
@@ -273,8 +380,6 @@ AliT0Parameters::GetChannel(Int_t trm,  Int_t tdc, Int_t chain, Int_t channel)
 //__________________________________________________________________
 TMap *AliT0Parameters::GetMapLookup()
 {
-  // get LookUp Table
-
   if (!fgLookUp){
     cout<<" No look up table in OCDB";
     return 0;