**************************************************************************/
/*
-$Log$
+$Log: AliT0Preprocessor.cxx,v $
+Revision 1.8 2007/12/07 15:22:51 alla
+bug fixed by Alberto
+
+Revision 1.7 2007/12/06 16:35:24 alla
+new bugs fixed by Tomek
+
Revision 1.5 2007/11/23 19:28:52 alla
bug fixed
Version 1.1 2006/10
Preliminary test version (T.Malkiewicz)
*/
-
// T0 preprocessor:
-// 1) takes data from DCS and passes it to the class AliTOFDataDCS for processing and writes the result to the Reference DB.
-// 2) takes data form DAQ (both from Laser Calibration and Physics runs), processes it, and stores either to OCDB or to Reference DB.
+// 1) takes data from DCS and passes it to the class AliTOFDataDCS
+// for processing and writes the result to the Reference DB.
+// 2) takes data form DAQ (both from Laser Calibration and Physics runs),
+// processes it, and stores either to OCDB or to Reference DB.
+
#include "AliT0Preprocessor.h"
#include "AliT0DataDCS.h"
#include "AliCDBMetaData.h"
#include "AliDCSValue.h"
+#include "AliCDBEntry.h"
#include "AliLog.h"
#include <TTimeStamp.h>
#include <TObjString.h>
#include <TNamed.h>
#include "AliT0Dqclass.h"
+#include "TClass.h"
ClassImp(AliT0Preprocessor)
//____________________________________________________
-AliT0Preprocessor::AliT0Preprocessor(AliShuttleInterface* shuttle) : AliPreprocessor("T00", shuttle), fData(0)
+AliT0Preprocessor::AliT0Preprocessor(AliShuttleInterface* shuttle) :
+ AliPreprocessor("T00", shuttle),
+ fData(0)
{
//constructor
+ AddRunType("PHYSICS");
+ AddRunType("STANDALONE");
+ // AddRunType("LASER");
}
//____________________________________________________
AliT0Preprocessor::~AliT0Preprocessor()
{
+ //destructor
delete fData;
fData = 0;
}
void AliT0Preprocessor::Initialize(Int_t run, UInt_t startTime, UInt_t endTime)
{
+ // Creates AliT0DataDCS object
AliPreprocessor::Initialize(run, startTime, endTime);
AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run, TTimeStamp(startTime).AsString(), TTimeStamp(endTime).AsString()));
- fData = new AliT0DataDCS(fRun, fStartTime, fEndTime);
+ fData = new AliT0DataDCS(fRun, fStartTime, fEndTime, GetStartTimeDCSQuery(), GetEndTimeDCSQuery());
}
//____________________________________________________
-UInt_t AliT0Preprocessor::Process(TMap* dcsAliasMap )
-{
- // T0 preprocessor return codes:
- // return=0 : all ok
- // return=1 : no DCS input data
- // return=2 : failed to store DCS data
- // return=3 : no Laser data (Walk correction)
- // return=4 : failed to store OCDB time equalized data
- // return=5 : no DAQ input for OCDB
- // return=6 : failed to retrieve DAQ data from OCDB
- // return=7 : failed to store T0 OCDB data
+Bool_t AliT0Preprocessor::ProcessDCS(){
+ // Check whether DCS should be processed or not...
+ TString runType = GetRunType();
+ Log(Form("ProcessDCS - RunType: %s",runType.Data()));
+
+ if((runType == "STANDALONE")||
+ (runType == "PHYSICS") ) {
+ // || (runType == "LASER")){
+ // return kFALSE;
+ return kTRUE;
+ }else{
+ return kFALSE;
+ }
+}
+//____________________________________________________
+
+UInt_t AliT0Preprocessor::ProcessDCSDataPoints(TMap* dcsAliasMap){
+ // Fills data into AliT0DataDCS object
+ Log("Processing DCS DP");
Bool_t resultDCSMap=kFALSE;
Bool_t resultDCSStore=kFALSE;
- Bool_t resultLaser=kFALSE;
- Bool_t resultOnline=kFALSE;
-
+
if(!dcsAliasMap)
{
Log("No DCS input data");
}
else
{
- /*
resultDCSMap=fData->ProcessData(*dcsAliasMap);
if(!resultDCSMap)
{
}
else
{
- Float_t *meanScaler[24] = fData->GetScalerMean();
-
AliCDBMetaData metaDataDCS;
metaDataDCS.SetBeamPeriod(0);
metaDataDCS.SetResponsible("Tomasz Malkiewicz");
metaDataDCS.SetComment("This preprocessor fills an AliTODataDCS object.");
AliInfo("Storing DCS Data");
- resultDCSStore = Store("Calib","DCSData",meanScaler, &metaDataDCS);
+ resultDCSStore = StoreReferenceData("Calib","DCSData",fData, &metaDataDCS);
if (!resultDCSStore)
{
Log("Some problems occurred while storing DCS data results in ReferenceDB");
return 2;// return error Code for processed DCS data not stored
}
-
}
- */
}
+ return 0;
+}
+//____________________________________________________
- // processing DAQ
+UInt_t AliT0Preprocessor::ProcessLaser(){
+ // Processing data from DAQ Standalone run
+ Log("Processing Laser calibration - Walk Correction");
+
+ //Retrieve the last T0 calibration object
- TString runType = GetRunType();
+ Float_t parqtcold[24][2], parledold[24][2],parqtcnew[24][2], parlednew[24][2] , goodled[24][2], goodqtc[24][2];
- if(runType == "T0_STANDALONE_LASER")
- {
- TList* list = GetFileSources(kDAQ, "LASER");
- if (list)
- {
- TIter iter(list);
- TObjString *source;
- while ((source = dynamic_cast<TObjString *> (iter.Next())))
- {
- const char *laserFile = GetFile(kDAQ, "LASER", source->GetName());
- if (laserFile)
+ // std::cout<<"sizeof "<<sizeof(parqtcold)<<std::endl;
+ memset(parqtcold, 0, sizeof(parqtcold));
+ memset(parqtcnew, 0, sizeof(parqtcnew));
+ memset(parledold, 0, sizeof(parledold));
+ memset(parlednew, 0, sizeof(parlednew));
+ Int_t iStore=0;
+
+
+ AliT0CalibWalk* clb=0;
+ AliCDBEntry* entryCalib = GetFromOCDB("Calib", "Slewing_Walk");
+ if(!entryCalib)
+ Log(Form("Cannot find any AliCDBEntry for [Calib, SlewingWalk]!"));
+ else {
+ clb =dynamic_cast<AliT0CalibWalk*>(entryCalib->GetObject());
+ for(Int_t i=0; i<24; i++)
+ {
+ for(Int_t ipar=0; ipar<2; ipar++)
+ {
+ // std::cout<<"parqtcold "<<parqtcold[i][ipar]<<std::endl;
+ parqtcold[i][ipar] = clb->GetQTCpar(i,ipar);
+ parledold[i][ipar] = clb->GetLEDpar(i, ipar);
+ goodqtc[i][ipar] = 999;
+ goodled[i][ipar] = 999;
+ // std:: cout<<" old "<<i<<" "<<ipar<<" qtc "<< parqtcold[i][ipar]<<" led "<<parledold[i][ipar]<< std::endl;
+ }
+ }
+ }
+
+
+ Bool_t resultLaser=kFALSE;
+ //processing DAQ
+ TList* list = GetFileSources(kDAQ, "LASER");
+ if (list)
+ {
+ TIter iter(list);
+ TObjString *source;
+ while ((source = dynamic_cast<TObjString *> (iter.Next())))
+ {
+ const char *laserFile = GetFile(kDAQ, "LASER", source->GetName());
+ if (laserFile)
{
Log(Form("File with Id LASER found in source %s!", source->GetName()));
AliT0CalibWalk *laser = new AliT0CalibWalk();
laser->MakeWalkCorrGraph(laserFile);
- AliCDBMetaData metaData;
- metaData.SetBeamPeriod(0);
+ //check difference with what was before
+ if(laser && clb ){
+ iStore = 1;
+ for(Int_t i=0; i<24; i++)
+ {
+ for(Int_t ifit=0; ifit<2; ifit++)
+ {
+ parqtcnew[i][ifit] = laser->GetQTCpar(i,ifit);
+ if( parqtcold[i][ifit] != 0 && parqtcnew[i][ifit] !=0)
+ {
+ goodqtc[i][ifit] =
+ (parqtcnew[i][ifit] - parqtcold[i][ifit])/parqtcold[i][ifit];
+ // std::cout<<"qtc "<<i<<" "<<ifit<<" "<< goodqtc[i][ifit]<< std::endl;
+ }
+ parlednew[i][ifit] = laser->GetLEDpar(i,ifit);
+ if(parledold[i][ifit] != 0 && parlednew[i][ifit]!= 0 )
+ {
+ goodled[i][ifit]=
+ (parlednew[i][ifit] - parledold[i][ifit])/parledold[i][ifit];
+ // std::cout<<"led "<<i<<" "<<ifit<<" "<< goodled[i][ifit]<< std::endl;
+ }
+ // if(TMath::Abs(goodqtc[i][ifit])>0.1 ||
+ // TMath::Abs(goodled[i][ifit])>0.1)
+ // iStore = 0;
+ }
+ }
+ }
+
+ AliCDBMetaData metaData;
+ metaData.SetBeamPeriod(0);
metaData.SetResponsible("Tomek&Michal");
metaData.SetComment("Walk correction from laser runs.");
- //TObjArray* arrLaser = laser->GetfWalk();
- resultLaser=Store("Calib","Slewing_Walk", laser, &metaData, 0, 1);
+ if( iStore>0)
+ resultLaser=Store("Calib","Slewing_Walk", laser, &metaData, 0, 1);
delete laser;
+ Log(Form("resultLaser = %d",resultLaser));
}
else
{
Log("No Laser Data stored");
return 3;//return error code for failure in storing Laser Data
}
- }
- }
- else if(runType == "PHYSICS")
- {
- TList* listPhys = GetFileSources(kDAQ, "PHYSICS");
- if (listPhys)
+ } else {
+ Log("No sources found for id LASER!");
+ return 1;
+ }
+ return 0;
+}
+
+//____________________________________________________
+
+UInt_t AliT0Preprocessor::ProcessPhysics(){
+ //Processing data from DAQ Physics run
+ Log("Processing Physics");
+
+ Bool_t resultOnline=kFALSE;
+ //processing DAQ
+ TList* listPhys = GetFileSources(kDAQ, "PHYSICS");
+ if (listPhys)
{
TIter iter(listPhys);
TObjString *sourcePhys;
{
AliT0CalibTimeEq *online = new AliT0CalibTimeEq();
online->Reset();
- online->ComputeOnlineParams("CFD", 20, 4., filePhys);
+ Bool_t writeok = online->ComputeOnlineParams(filePhys);
AliCDBMetaData metaData;
metaData.SetBeamPeriod(0);
- metaData.SetResponsible("Tomek&Michal");
+ metaData.SetResponsible("Alla Maevskaya");
metaData.SetComment("Time equalizing result.");
- resultOnline = Store("Calib","TimeDelay", online, &metaData, 0, 1);
+
+ if (writeok) resultOnline = Store("Calib","TimeDelay", online, &metaData, 0, 1);
+ else {
+
+ Log(Form("writeok = %d not enough data for equalizing",resultOnline));
+ return 0;
+ }
+ Log(Form("resultOnline = %d",resultOnline));
delete online;
}
- else
+ else
{
Log(Form("Could not find file with Id PHYSICS in source %s!", sourcePhys->GetName()));
return 1;
}
+
}
if (!resultOnline)
{
- Log("No Laser Data stored");
+ Log("No Data stored");
return 4;//return error code for failure in storing OCDB Data
}
- }
- }
-
- return 0;
+ } else {
+ Log("No sources found for id PHYSICS!");
+ return 1;
+ }
+ return 0;
}
+//____________________________________________________
+UInt_t AliT0Preprocessor::Process(TMap* dcsAliasMap )
+{
+ // T0 preprocessor return codes:
+ // return=0 : all ok
+ // return=1 : no DCS input data
+ // return=2 : failed to store DCS data
+ // return=3 : no Laser data (Walk correction)
+ // return=4 : failed to store OCDB time equalized data
+ // return=5 : no DAQ input for OCDB
+ // return=6 : failed to retrieve DAQ data from OCDB
+ // return=7 : failed to store T0 OCDB data
+ // return=8 : not enough data for equalizing
+ Bool_t dcsDP = ProcessDCS();
+ Log(Form("dcsDP = %d",dcsDP));
+ TString runType = GetRunType();
+ Log(Form("RunType: %s",runType.Data()));
+ //processing
+ if(runType == "STANDALONE"){
+ if(dcsDP==1){
+ Int_t iresultDCS = ProcessDCSDataPoints(dcsAliasMap);
+ return iresultDCS;
+ }
+ }
+ /*
+ if(runType == "LASER"){
+ Int_t iresultLaser = ProcessLaser();
+ if(dcsDP==1){
+ Int_t iresultDCS = ProcessDCSDataPoints(dcsAliasMap);
+ return iresultDCS;
+ }
+
+ Log(Form("iresultLaser = %d",iresultLaser));
+ return iresultLaser;
+ }
+ */
+ else if(runType == "PHYSICS"){
+ Int_t iresultPhysics = ProcessPhysics();
+ if(dcsDP==1){
+ Int_t iresultDCS = ProcessDCSDataPoints(dcsAliasMap);
+ return iresultDCS;
+ }
+ Log(Form("iresultPhysics = %d",iresultPhysics));
+ return iresultPhysics;
+ }
+
+
+
+ return 0;
+}