#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 <AliGeomManager.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
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
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;
}
// 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;
}
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);
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);
//__________________________________________________________________
-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,
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);
//__________________________________________________________________
TMap *AliT0Parameters::GetMapLookup()
{
- // get LookUp Table
-
if (!fgLookUp){
cout<<" No look up table in OCDB";
return 0;