for pass0
[u/mrichter/AliRoot.git] / T0 / AliT0PreprocessorOffline.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, 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 // T0 preprocessor:
20 // 2) takes data after  pass0 , 
21 // processes it, and stores either to OCDB .
22
23
24 #include "AliT0PreprocessorOffline.h"
25 #include "AliT0RecoParam.h"
26 #include "AliT0CalibTimeEq.h"
27 #include "AliCDBStorage.h"
28 #include "AliCDBMetaData.h"
29 #include "AliCDBManager.h"
30 #include "AliCTPTimeParams.h"
31 #include "AliLHCClockPhase.h"
32 #include "AliT0CalibSeasonTimeShift.h"
33 #include "AliT0CalibLatency.h"
34 #include "AliCDBEntry.h"
35 #include "AliLog.h"
36
37 #include <TTimeStamp.h>
38 #include <TFile.h>
39 #include <TObjString.h>
40 #include <TNamed.h>
41 #include "TClass.h"
42
43
44 ClassImp(AliT0PreprocessorOffline)
45
46 //____________________________________________________
47 AliT0PreprocessorOffline::AliT0PreprocessorOffline():
48 TNamed("AliT0PreprocessorOffline","AliT0PreprocessorOffline"),
49   startRun(0),      
50   endRun(0),  
51   startTime(0),    
52   endTime(0),     
53   ocdbStorage("")  
54 {
55   //constructor
56 }
57 //____________________________________________________
58
59 AliT0PreprocessorOffline::~AliT0PreprocessorOffline()
60 {
61   //destructor
62
63 }
64 //____________________________________________________
65 void AliT0PreprocessorOffline::Process(TString filePhysName, Int_t ustartRun, Int_t uendRun, TString pocdbStorage)
66 {
67   CalibOffsetChannels(filePhysName, ustartRun, uendRun, pocdbStorage);
68   CalibT0sPosition(filePhysName, ustartRun, uendRun, pocdbStorage);
69 }
70 //____________________________________________________
71
72 void AliT0PreprocessorOffline::CalibOffsetChannels(TString filePhysName, Int_t ustartRun, Int_t uendRun, TString pocdbStorage)
73 {
74
75   Float_t zero_timecdb[24]={0};
76   Float_t *timecdb = zero_timecdb;
77   Float_t *cfdvalue = zero_timecdb;
78   Int_t badpmt=0;
79   //Processing data from DAQ Physics run
80   AliInfo("Processing Time Offset between channels");
81   if (pocdbStorage.Length()>0) ocdbStorage=pocdbStorage;
82   else
83     ocdbStorage="local://"+gSystem->GetFromPipe("pwd")+"/OCDB";
84   AliCDBEntry *entry = AliCDBManager::Instance()->Get("GRP/CTP/CTPtiming");
85   if (!entry) AliFatal("CTP timing parameters are not found in OCDB !");
86   AliCTPTimeParams *ctpParams = (AliCTPTimeParams*)entry->GetObject();
87   Float_t l1Delay = (Float_t)ctpParams->GetDelayL1L0()*25.0;
88   
89   AliCDBEntry *entry1 = AliCDBManager::Instance()->Get("GRP/CTP/TimeAlign");
90   if (!entry1) AliFatal("CTP time-alignment is not found in OCDB !");
91   AliCTPTimeParams *ctpTimeAlign = (AliCTPTimeParams*)entry1->GetObject();
92   l1Delay += ((Float_t)ctpTimeAlign->GetDelayL1L0()*25.0);
93   
94   AliCDBEntry *entry4 = AliCDBManager::Instance()->Get("GRP/Calib/LHCClockPhase");
95   if (!entry4) AliFatal("LHC clock-phase shift is not found in OCDB !");
96   AliLHCClockPhase *phase = (AliLHCClockPhase*)entry4->GetObject();
97   Float_t fGRPdelays = l1Delay - phase->GetMeanPhase();
98
99   AliCDBEntry *entryL = AliCDBManager::Instance()->Get("T0/Calib/Latency");
100   if (!entryL) AliFatal(" T0 Latency is not found in OCDB !");
101   AliT0CalibLatency *lat = (AliT0CalibLatency*)entryL->GetObject();
102   Float_t fLatencyHPTDC = lat->GetLatencyHPTDC();
103   Float_t fLatencyL1 = lat->GetLatencyL1();
104
105   AliCDBEntry *entryCalib = AliCDBManager::Instance()->Get("T0/Calib/TimeDelay");
106   if(!entryCalib) {
107     AliWarning(Form("Cannot find any AliCDBEntry for [Calib, TimeDelay]!"));
108   }
109   else
110     {
111       AliT0CalibTimeEq *clb = (AliT0CalibTimeEq*)entryCalib->GetObject();
112        timecdb = clb->GetTimeEq();
113        cfdvalue = clb->GetCFDvalue();
114        for (Int_t i=0; i<24; i++) {
115          if( cfdvalue[i] < 500 || cfdvalue[i] > 30000) cfdvalue[i] =( 1000.*fLatencyHPTDC - 1000.*fLatencyL1 + 1000.*fGRPdelays)/24.4;
116          printf("Calc  mean CFD time %i %f \n",i,cfdvalue[i]);
117        }
118     }
119   //AliCDBEntry *entryCalibreco = AliCDBManager::Instance()->Get("T0/Calib/RecoParam");
120   //  if(entryCalibreco) {
121   //  AliT0RecoParam *rpr = (AliT0RecoParam*) entryCalibreco->GetObject();
122     //    badpmt = rpr->GetRefPoint();
123     badpmt = -1;
124     // }
125   AliT0CalibTimeEq *offline = new AliT0CalibTimeEq();
126   Bool_t writeok = offline->ComputeOfflineParams(filePhysName.Data(), timecdb, cfdvalue, badpmt);
127   AliCDBMetaData metaData;
128   metaData.SetBeamPeriod(1);
129   metaData.SetResponsible("Alla Maevskaya");
130   metaData.SetComment("Time equalizing result with slew");
131       
132   if (writeok)  {
133     AliCDBId* id1=NULL;
134     id1=new AliCDBId("T0/Calib/TimeDelay", ustartRun, uendRun );
135     AliCDBStorage* gStorage = AliCDBManager::Instance()->GetStorage(ocdbStorage);
136     gStorage->Put(offline, (*id1), &metaData);
137   }
138   else {
139     
140     AliWarning(Form("writeok = %d data is not OK to be in OCDB",writeok));
141       }           
142   
143   delete offline;
144  
145 }
146 //-------------------------------------------------------------------------------------
147 void AliT0PreprocessorOffline::CalibT0sPosition(TString filePhysName, Int_t ustartRun, Int_t uendRun, TString pocdbStorage)
148 {
149   Float_t zero_timecdb[4]={0};
150   Float_t *timecdb = zero_timecdb;
151   if (pocdbStorage.Length()>0) ocdbStorage=pocdbStorage;
152   else
153     ocdbStorage="local://"+gSystem->GetFromPipe("pwd")+"/OCDB";
154   AliCDBEntry *entryCalib = AliCDBManager::Instance()->Get("T0/Calib/TimeAdjust");
155   if(!entryCalib) {
156     AliWarning(Form("Cannot find any AliCDBEntry for [Calib, TimeAdjust]!"));
157   }
158   else
159     {
160       AliT0CalibSeasonTimeShift *clb = (AliT0CalibSeasonTimeShift*)entryCalib->GetObject();
161       timecdb = clb->GetT0Means();
162     }
163   
164   AliT0CalibSeasonTimeShift *offline = new AliT0CalibSeasonTimeShift();
165   Bool_t writeok = offline->SetT0Par(filePhysName.Data(), timecdb);
166   AliCDBMetaData metaData;
167   metaData.SetBeamPeriod(1);
168   metaData.SetResponsible("Alla Maevskaya");
169   metaData.SetComment("Time equalizing result with slew");
170   
171   if (writeok)  {
172     AliCDBId* id1=NULL;
173     id1=new AliCDBId("T0/Calib/TimeAdjust", ustartRun, uendRun);
174     AliCDBStorage* gStorage = AliCDBManager::Instance()->GetStorage(ocdbStorage);
175     gStorage->Put(offline, (*id1), &metaData);
176   }
177
178 }