X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=TPC%2FAliTPCPreprocessor.cxx;h=603bea4858567883c359cc0d50e82188fedc97a9;hb=a69e1a34d55dd7c39794d550e9efb9a4db405b4e;hp=ceba9134d170faa9c6127ab84c992437fc944b6c;hpb=342d3b42b8e44fa3d4c78b679a2dbd68249a781d;p=u%2Fmrichter%2FAliRoot.git diff --git a/TPC/AliTPCPreprocessor.cxx b/TPC/AliTPCPreprocessor.cxx index ceba9134d17..603bea48585 100644 --- a/TPC/AliTPCPreprocessor.cxx +++ b/TPC/AliTPCPreprocessor.cxx @@ -27,6 +27,9 @@ #include "AliTPCCalibPedestal.h" #include "AliTPCCalibPulser.h" #include "AliTPCCalibCE.h" +#include "AliTPCCalibRaw.h" +#include "AliTPCdataQA.h" +#include "ARVersion.h" #include "TFile.h" #include "TTree.h" #include "TGraph.h" @@ -37,16 +40,27 @@ const Int_t kValCutTemp = 100; // discard temperatures > 100 degrees const Int_t kDiffCutTemp = 5; // discard temperature differences > 5 degrees +const Double_t kHighVoltageDifference = 1e-4; // don't record High Voltage points + // differing by less than 1e-4 from + // previous point. const TString kPedestalRunType = "PEDESTAL"; // pedestal run identifier -const TString kPulserRunType = "PULSER"; // pulser run identifier +const TString kPulserRunType = "PULSER"; // pulser run identifier const TString kPhysicsRunType = "PHYSICS"; // physics run identifier -const TString kStandAloneRunType = "STANDALONE"; // standalone run identifier -const TString kCosmicRunType = "COSMIC"; // cosmic run identifier -const TString kLaserRunType = "LASER"; // laser run identifier +const TString kCosmicRunType = "COSMIC"; // cosmic run identifier +const TString kLaserRunType = "LASER"; // laser run identifier const TString kDaqRunType = "DAQ"; // DAQ run identifier const TString kAmandaTemp = "TPC_PT_%d_TEMPERATURE"; // Amanda string for temperature entries +const TString kAmandaDDL = "DDL%d"; // Amanda string for list of active DDLs +const Int_t kNumDDL = 216; // number of TPC DDLs +const Int_t kFirstDDL = 768; // identifier of first DDL //const Double_t kFitFraction = 0.7; // Fraction of DCS sensor fits required const Double_t kFitFraction = -1.0; // Don't require minimum number of fits in commissioning run +const Int_t kNumPressureSensors = 3; // number of pressure sensors +const char* kPressureSensorNames[kNumPressureSensors] = { + "CavernAtmosPressure", + "CavernAtmosPressure2", + "SurfaceAtmosPressure" }; + // // This class is the SHUTTLE preprocessor for the TPC detector. @@ -57,7 +71,8 @@ ClassImp(AliTPCPreprocessor) //______________________________________________________________________________________________ AliTPCPreprocessor::AliTPCPreprocessor(AliShuttleInterface* shuttle) : AliPreprocessor("TPC",shuttle), - fConfEnv(0), fTemp(0), fHighVoltage(0), fHighVoltageStat(0), fConfigOK(kTRUE), fROC(0) + fConfEnv(0), fTemp(0), fHighVoltage(0), fHighVoltageStat(0), fGoofie(0), + fPressure(0), fConfigOK(kTRUE), fROC(0) { // constructor fROC = AliTPCROC::Instance(); @@ -67,7 +82,6 @@ AliTPCPreprocessor::AliTPCPreprocessor(AliShuttleInterface* shuttle) : AddRunType(kPedestalRunType); AddRunType(kPulserRunType); AddRunType(kPhysicsRunType); - AddRunType(kStandAloneRunType); AddRunType(kCosmicRunType); AddRunType(kLaserRunType); AddRunType(kDaqRunType); @@ -76,7 +90,8 @@ AliTPCPreprocessor::AliTPCPreprocessor(AliShuttleInterface* shuttle) : //______________________________________________________________________________________________ AliTPCPreprocessor::AliTPCPreprocessor(const AliTPCPreprocessor& ) : AliPreprocessor("TPC",0), - fConfEnv(0), fTemp(0), fHighVoltage(0), fHighVoltageStat(0), fConfigOK(kTRUE), fROC(0) + fConfEnv(0), fTemp(0), fHighVoltage(0), fHighVoltageStat(0), fGoofie(0), + fPressure(0), fConfigOK(kTRUE), fROC(0) { Fatal("AliTPCPreprocessor", "copy constructor not implemented"); @@ -91,6 +106,9 @@ AliTPCPreprocessor::~AliTPCPreprocessor() delete fTemp; delete fHighVoltage; + delete fHighVoltageStat; + delete fGoofie; + delete fPressure; } //______________________________________________________________________________________________ AliTPCPreprocessor& AliTPCPreprocessor::operator = (const AliTPCPreprocessor& ) @@ -104,11 +122,8 @@ AliTPCPreprocessor& AliTPCPreprocessor::operator = (const AliTPCPreprocessor& ) void AliTPCPreprocessor::Initialize(Int_t run, UInt_t startTime, UInt_t endTime) { - // Creates AliTestDataDCS object -- start maps half an hour beforre actual run start - UInt_t startTimeLocal = startTime-1800; - - AliPreprocessor::Initialize(run, startTimeLocal, endTime); + AliPreprocessor::Initialize(run, startTime, endTime); AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run, TTimeStamp((time_t)startTime,0).AsString(), @@ -138,7 +153,7 @@ void AliTPCPreprocessor::Initialize(Int_t run, UInt_t startTime, fConfigOK = kFALSE; return; } - fTemp = new AliTPCSensorTempArray(startTimeLocal, fEndTime, confTree, kAmandaTemp); + fTemp = new AliTPCSensorTempArray(startTime, endTime, confTree, kAmandaTemp); fTemp->SetValCut(kValCutTemp); fTemp->SetDiffCut(kDiffCutTemp); } @@ -157,10 +172,17 @@ void AliTPCPreprocessor::Initialize(Int_t run, UInt_t startTime, fConfigOK = kFALSE; return; } - fHighVoltage = new AliDCSSensorArray(startTimeLocal, fEndTime, confTree); + time_t timeStart = (time_t)(((TString)GetRunParameter("DAQ_time_start")).Atoi()); + time_t timeEnd = (time_t)(((TString)GetRunParameter("DAQ_time_end")).Atoi()); + fHighVoltage = new AliDCSSensorArray (UInt_t(timeStart), + UInt_t(timeEnd), confTree); + } // High voltage status values + TString hvStatConf = fConfEnv->GetValue("HighVoltageStat","ON"); + hvStatConf.ToUpper(); + if (hvStatConf != "OFF" ) { confTree=0; entry=0; entry = GetFromOCDB("Config", "HighVoltageStat"); @@ -170,8 +192,43 @@ void AliTPCPreprocessor::Initialize(Int_t run, UInt_t startTime, fConfigOK = kFALSE; return; } - fHighVoltageStat = new AliDCSSensorArray(startTimeLocal, fEndTime, confTree); + fHighVoltageStat = new AliDCSSensorArray(startTime, endTime, confTree); } + + // Goofie values + + TString goofieConf = fConfEnv->GetValue("Goofie","ON"); + goofieConf.ToUpper(); + if (goofieConf != "OFF" ) { + confTree=0; + entry=0; + entry = GetFromOCDB("Config", "Goofie"); + if (entry) confTree = (TTree*) entry->GetObject(); + if ( confTree==0 ) { + Log("AliTPCPreprocsessor: Goofie Config OCDB entry missing.\n"); + fConfigOK = kFALSE; + return; + } + fGoofie = new AliDCSSensorArray(startTime, endTime, confTree); + } + + // Pressure values + + TString runType = GetRunType(); + + if( runType == kPhysicsRunType || + runType == kLaserRunType ) { + TString pressureConf = fConfEnv->GetValue("Pressure","ON"); + pressureConf.ToUpper(); + if (pressureConf != "OFF" ) { + TClonesArray * array = new TClonesArray("AliDCSSensor",kNumPressureSensors); + for(Int_t j = 0; j < kNumPressureSensors; j++) { + AliDCSSensor * sens = new ((*array)[j])AliDCSSensor; + sens->SetStringID(kPressureSensorNames[j]); + } + fPressure = new AliDCSSensorArray(startTime, endTime, array); + } + } } //______________________________________________________________________________________________ @@ -190,44 +247,65 @@ UInt_t AliTPCPreprocessor::Process(TMap* dcsAliasMap) TObject * status; UInt_t dcsResult=0; - if (errorHandling == "OFF" ) { - if (!dcsAliasMap) dcsResult=1; - if (dcsAliasMap->GetEntries() == 0 ) dcsResult=1; - status = new TParameter("dcsResult",dcsResult); - resultArray->Add(status); - } else { - if (!dcsAliasMap) return 9; - if (dcsAliasMap->GetEntries() == 0 ) return 9; - } + if (!dcsAliasMap) dcsResult=1; + if (dcsAliasMap->GetEntries() == 0 ) dcsResult=1; + status = new TParameter("dcsResult",dcsResult); + resultArray->Add(status); - - TString runType = GetRunType(); + if ( dcsResult == 0 ) { + // Temperature sensors are processed by AliTPCCalTemp - TString tempConf = fConfEnv->GetValue("Temperature","ON"); - tempConf.ToUpper(); - if (tempConf != "OFF" ) { - UInt_t tempResult = MapTemperature(dcsAliasMap); - result=tempResult; - status = new TParameter("tempResult",tempResult); - resultArray->Add(status); - } + TString tempConf = fConfEnv->GetValue("Temperature","ON"); + tempConf.ToUpper(); + if (tempConf != "OFF" ) { + UInt_t tempResult = MapTemperature(dcsAliasMap); + if ( tempConf != "TRY") result+=tempResult; + status = new TParameter("tempResult",tempResult); + resultArray->Add(status); + } // High Voltage recordings - TString hvConf = fConfEnv->GetValue("HighVoltage","ON"); - hvConf.ToUpper(); - if (hvConf != "OFF" ) { - UInt_t hvResult = MapHighVoltage(dcsAliasMap); - result+=hvResult; - status = new TParameter("hvResult",hvResult); - resultArray->Add(status); - } + TString hvConf = fConfEnv->GetValue("HighVoltage","ON"); + hvConf.ToUpper(); + if (hvConf != "OFF" ) { + UInt_t hvResult = MapHighVoltage(dcsAliasMap); + if (hvConf != "TRY") result+=hvResult; + status = new TParameter("hvResult",hvResult); + resultArray->Add(status); + } + + // Goofie values + + TString goofieConf = fConfEnv->GetValue("Goofie","ON"); + goofieConf.ToUpper(); + if (goofieConf != "OFF" ) { + UInt_t goofieResult = MapGoofie(dcsAliasMap); + if (goofieConf != "TRY") result+=goofieResult; + status = new TParameter("goofieResult",goofieResult); + resultArray->Add(status); + } + + // Pressure values + + if( runType == kPhysicsRunType || + runType == kLaserRunType ) { + + TString pressureConf = fConfEnv->GetValue("Pressure","ON"); + pressureConf.ToUpper(); + if (pressureConf != "OFF" ) { + UInt_t pressureResult = MapPressure(dcsAliasMap); + status = new TParameter("pressureResult",pressureResult); + resultArray->Add(status); + } + } + } // Other calibration information will be retrieved through FXS files // examples: // TList* fileSourcesDAQ = GetFile(AliShuttleInterface::kDAQ, "pedestals"); @@ -254,7 +332,7 @@ UInt_t AliTPCPreprocessor::Process(TMap* dcsAliasMap) if (source == "DAQHLT" ) numSources=2; UInt_t pedestalResult=0; for (Int_t i=0; i("pulserResult",pulserResult); resultArray->Add(status); } } +// raw calibration processing + + if(runType == kPhysicsRunType) { + Int_t numSources = 1; + Int_t rawSource[2] = {AliShuttleInterface::kDAQ,AliShuttleInterface::kHLT} ; + TString source = fConfEnv->GetValue("Raw","DAQ"); + source.ToUpper(); + if ( source != "OFF") { + if ( source == "HLT") rawSource[0] = AliShuttleInterface::kHLT; + if (!GetHLTStatus()) rawSource[0] = AliShuttleInterface::kDAQ; + if (source == "HLTDAQ" ) { + numSources=2; + rawSource[0] = AliShuttleInterface::kHLT; + rawSource[1] = AliShuttleInterface::kDAQ; + } + if (source == "DAQHLT" ) numSources=2; + if (source == "TRY" ) numSources=2; + UInt_t rawResult=0; + for (Int_t i=0; i("rawResult",rawResult); + resultArray->Add(status); + } + } + + + // Altro configuration + + + TString altroConf = fConfEnv->GetValue("AltroConf","ON"); + altroConf.ToUpper(); + if (altroConf != "OFF" ) { + UInt_t altroResult = ExtractAltro(AliShuttleInterface::kDCS,dcsAliasMap); + if (altroConf != "TRY" ) result+=altroResult; + status = new TParameter("altroResult",altroResult); + resultArray->Add(status); + } + // Central Electrode processing - if( runType == kPhysicsRunType || runType == kStandAloneRunType || - runType == kDaqRunType ) { + if( runType == kPhysicsRunType || + runType == kLaserRunType ) { -// if (true) { // do CE processing for all run types Int_t numSources = 1; Int_t ceSource[2] = {AliShuttleInterface::kDAQ,AliShuttleInterface::kHLT} ; TString source = fConfEnv->GetValue("CE","DAQ"); @@ -311,28 +430,61 @@ UInt_t AliTPCPreprocessor::Process(TMap* dcsAliasMap) ceSource[1] = AliShuttleInterface::kDAQ; } if (source == "DAQHLT" ) numSources=2; + if (source == "TRY" ) numSources=2; UInt_t ceResult=0; for (Int_t i=0; i("ceResult",ceResult); resultArray->Add(status); + + numSources = 1; + Int_t qaSource[2] = {AliShuttleInterface::kDAQ,AliShuttleInterface::kHLT} ; + source = fConfEnv->GetValue("QA","DAQ"); + source.ToUpper(); + if ( source != "OFF" ) { + if ( source == "HLT") qaSource[0] = AliShuttleInterface::kHLT; + if (!GetHLTStatus()) qaSource[0] = AliShuttleInterface::kDAQ; + if (source == "HLTDAQ" ) { + numSources=2; + qaSource[0] = AliShuttleInterface::kHLT; + qaSource[1] = AliShuttleInterface::kDAQ; + } + if (source == "DAQHLT" ) numSources=2; + if (source == "TRY" ) numSources=2; + UInt_t qaResult=0; + for (Int_t i=0; i("qaResult",qaResult); + resultArray->Add(status); } + } } - if (errorHandling == "OFF" ) { - AliCDBMetaData metaData; - metaData.SetBeamPeriod(0); - metaData.SetResponsible("Haavard Helstrup"); - metaData.SetComment("Preprocessor AliTPC status."); - Store("Calib", "PreprocStatus", resultArray, &metaData, 0, kFALSE); - resultArray->Delete(); - return 0; - } else { - return result; - } +// Store component status to OCDB + + AliCDBMetaData metaData; + metaData.SetBeamPeriod(0); + metaData.SetResponsible("Haavard Helstrup"); + metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH); + metaData.SetComment("Preprocessor AliTPC status."); + Store("Calib", "PreprocStatus", resultArray, &metaData, 0, kFALSE); + resultArray->Delete(); + delete resultArray; + + if (errorHandling == "OFF" ) return 0; + return result; + } //______________________________________________________________________________________________ UInt_t AliTPCPreprocessor::MapTemperature(TMap* dcsAliasMap) @@ -362,6 +514,7 @@ UInt_t AliTPCPreprocessor::MapTemperature(TMap* dcsAliasMap) AliCDBMetaData metaData; metaData.SetBeamPeriod(0); metaData.SetResponsible("Haavard Helstrup"); + metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH); metaData.SetComment("Preprocessor AliTPC data base entries."); Bool_t storeOK = Store("Calib", "Temperature", fTemp, &metaData, 0, kFALSE); @@ -372,49 +525,74 @@ UInt_t AliTPCPreprocessor::MapTemperature(TMap* dcsAliasMap) return result; } - //______________________________________________________________________________________________ -UInt_t AliTPCPreprocessor::MapHighVoltage(TMap* dcsAliasMap) +UInt_t AliTPCPreprocessor::MapPressure(TMap* dcsAliasMap) { - // extract DCS HV maps. Perform fits to save space + // extract DCS pressure maps. Perform fits to save space UInt_t result=0; - TMap *map = fHighVoltage->ExtractDCS(dcsAliasMap); + TMap *map = fPressure->ExtractDCS(dcsAliasMap); if (map) { - fHighVoltage->MakeSplineFit(map); - Double_t fitFraction = 1.0*fHighVoltage->NumFits()/fHighVoltage->NumSensors(); + fPressure->MakeSplineFit(map); + Double_t fitFraction = 1.0*fPressure->NumFits()/fPressure->NumSensors(); if (fitFraction > kFitFraction ) { - AliInfo(Form("High voltage recordings extracted, fits performed.\n")); + AliInfo(Form("Pressure values extracted, fits performed.\n")); } else { - Log ("Too few high voltage recordings fitted. \n"); + Log ("Too few pressure maps fitted. \n"); result = 9; } } else { - Log("No high voltage recordings extracted. \n"); + Log("No pressure map extracted. \n"); result=9; } delete map; + return result; +} - TMap *map2 = fHighVoltageStat->ExtractDCS(dcsAliasMap); - if (map2) { - fHighVoltageStat->StoreGraph(map2); +//______________________________________________________________________________________________ +UInt_t AliTPCPreprocessor::MapHighVoltage(TMap* dcsAliasMap) +{ + + // extract DCS HV maps. Perform fits to save space + + UInt_t result=0; + TMap *map = fHighVoltage->ExtractDCS(dcsAliasMap); + if (map) { + fHighVoltage->ClearFit(); + fHighVoltage->RemoveGraphDuplicates(kHighVoltageDifference); + // don't keep new point if too similar to previous one + fHighVoltage->SetGraph(map); } else { - Log("No high voltage status recordings extracted. \n"); + Log("No high voltage recordings extracted. \n"); result=9; } - delete map2; + delete map; - // add status maps to high voltage sensor array + TString hvStatConf = fConfEnv->GetValue("HighVoltageStat","ON"); + hvStatConf.ToUpper(); + if (hvStatConf != "OFF" ) { + TMap *map2 = fHighVoltageStat->ExtractDCS(dcsAliasMap); + if (map2) { + fHighVoltageStat->ClearFit(); + fHighVoltageStat->SetGraph(map2); + } else { + Log("No high voltage status recordings extracted. \n"); + result=9; + } + delete map2; - fHighVoltage->AddSensors(fHighVoltageStat); + // add status maps to high voltage sensor array + fHighVoltage->AddSensors(fHighVoltageStat); + } // Now store the final CDB file if ( result == 0 ) { AliCDBMetaData metaData; metaData.SetBeamPeriod(0); metaData.SetResponsible("Haavard Helstrup"); + metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH); metaData.SetComment("Preprocessor AliTPC data base entries."); Bool_t storeOK = Store("Calib", "HighVoltage", fHighVoltage, &metaData, 0, kFALSE); @@ -426,6 +604,42 @@ UInt_t AliTPCPreprocessor::MapHighVoltage(TMap* dcsAliasMap) } +//______________________________________________________________________________________________ +UInt_t AliTPCPreprocessor::MapGoofie(TMap* dcsAliasMap) +{ + + // extract DCS Goofie maps. Do not perform fits (low update rate) + + UInt_t result=0; + + TMap *map = fGoofie->ExtractDCS(dcsAliasMap); + if (map) { + fGoofie->ClearFit(); + fGoofie->SetGraph(map); + } else { + Log("No Goofie recordings extracted. \n"); + result=9; + } + delete map; + + // Now store the final CDB file + + if ( result == 0 ) { + AliCDBMetaData metaData; + metaData.SetBeamPeriod(0); + metaData.SetResponsible("Haavard Helstrup"); + metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH); + metaData.SetComment("Preprocessor AliTPC data base entries."); + + Bool_t storeOK = Store("Calib", "Goofie", fGoofie, &metaData, 0, kFALSE); + if ( !storeOK ) result=1; + + } + + return result; + +} + //______________________________________________________________________________________________ @@ -503,6 +717,7 @@ UInt_t AliTPCPreprocessor::ExtractPedestals(Int_t sourceFXS) AliCDBMetaData metaData; metaData.SetBeamPeriod(0); metaData.SetResponsible("Haavard Helstrup"); + metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH); metaData.SetComment("Preprocessor AliTPC data base entries."); Bool_t storeOK = Store("Calib", "Pedestals", calPadPed, &metaData, 0, kTRUE); @@ -511,10 +726,13 @@ UInt_t AliTPCPreprocessor::ExtractPedestals(Int_t sourceFXS) if ( !storeOK ) ++result; } } else { - Log ("Error: no entries in input file list!"); + Log ("Error: no entries in pedestal file list!"); result = 1; } + delete calPadPed; + delete calPadRMS; + return result; } @@ -608,71 +826,142 @@ UInt_t AliTPCPreprocessor::ExtractPulser(Int_t sourceFXS) AliCDBMetaData metaData; metaData.SetBeamPeriod(0); metaData.SetResponsible("Haavard Helstrup"); + metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH); metaData.SetComment("Preprocessor AliTPC data base entries."); Bool_t storeOK = Store("Calib", "Pulser", pulserObjects, &metaData, 0, kTRUE); if ( !storeOK ) ++result; } } else { - Log ("Error: no entries in input file list!"); + Log ("Error: no entries in pulser file list!"); result = 1; } + pulserObjects->Delete(); + delete pulserObjects; return result; } +//______________________________________________________________________________________________ + + +UInt_t AliTPCPreprocessor::ExtractRaw(Int_t sourceFXS) +{ + // + // Read Raw calibration file from file exchage server + // + + UInt_t result=0; + TObjArray* rawArray = new TObjArray; + + TList* list = GetFileSources(sourceFXS,"tpcCalibRaw"); + + if (list && list->GetEntries()>0) { + +// loop through all files + + UInt_t index = 0; + while (list->At(index)!=NULL) { + TObjString* fileNameEntry = (TObjString*) list->At(index); + if (fileNameEntry!=NULL) { + TString fileName = GetFile(sourceFXS, "tpcCalibRaw", + fileNameEntry->GetString().Data()); + TFile *f = TFile::Open(fileName); + if (!f) { + Log ("Error opening raw file."); + result =2; + break; + } + AliTPCCalibRaw *calRaw; + f->GetObject("tpcCalibRaw",calRaw); + if ( !calRaw ) { + Log ("No raw calibration object in file."); + result = 2; + break; + } + rawArray->Add(calRaw); + f->Close(); + } + ++index; + } // while(list) +// +// Store updated pedestal entry to OCDB +// + AliCDBMetaData metaData; + metaData.SetBeamPeriod(0); + metaData.SetResponsible("Haavard Helstrup"); + metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH); + metaData.SetComment("Preprocessor AliTPC data base entries."); + + Bool_t storeOK = Store("Calib", "Raw", rawArray, &metaData, 0, kTRUE); + if ( !storeOK ) ++result; + } else { + Log ("Error: no entries in raw file list!"); + result = 1; + } + + rawArray->Delete(); + delete rawArray; + + return result; +} +//______________________________________________________________________________________________ + UInt_t AliTPCPreprocessor::ExtractCE(Int_t sourceFXS) { // // Read Central Electrode file from file exchage server - // Keep original entry from OCDB in case no new CE calibration is available + // // - TObjArray *ceObjects=0; AliTPCCalPad *ceTmean=0; AliTPCCalPad *ceTrms=0; AliTPCCalPad *ceQmean=0; TObjArray *rocTtime=0; TObjArray *rocQtime=0; - AliCDBEntry* entry = GetFromOCDB("Calib", "CE"); - if (entry) ceObjects = (TObjArray*)entry->GetObject(); - if ( ceObjects==NULL ) { - Log("AliTPCPreprocsessor: No previous TPC central electrode entry available.\n"); - ceObjects = new TObjArray; - } + TObjArray *ceObjects= new TObjArray; + Int_t nSectors = fROC->GetNSectors(); - ceTmean = (AliTPCCalPad*)ceObjects->FindObject("CETmean"); - if ( !ceTmean ) { - ceTmean = new AliTPCCalPad("CETmean","CETmean"); - ceObjects->Add(ceTmean); - } - ceTrms = (AliTPCCalPad*)ceObjects->FindObject("CETrms"); - if ( !ceTrms ) { - ceTrms = new AliTPCCalPad("CETrms","CETrms"); - ceObjects->Add(ceTrms); - } - ceQmean = (AliTPCCalPad*)ceObjects->FindObject("CEQmean"); - if ( !ceQmean ) { - ceQmean = new AliTPCCalPad("CEQmean","CEQmean"); - ceObjects->Add(ceQmean); - } - //!new from here please have a look!!! - rocTtime = (TObjArray*)ceObjects->FindObject("rocTtime"); - if ( !rocTtime ) { - rocTtime = new TObjArray(nSectors); - rocTtime->SetName("rocTtime"); - ceObjects->Add(rocTtime); - } + ceTmean = new AliTPCCalPad("CETmean","CETmean"); + ceObjects->Add(ceTmean); + + ceTrms = new AliTPCCalPad("CETrms","CETrms"); + ceObjects->Add(ceTrms); + + ceQmean = new AliTPCCalPad("CEQmean","CEQmean"); + ceObjects->Add(ceQmean); + + rocTtime = new TObjArray(nSectors+2); // also make room for A and C side average + rocTtime->SetName("rocTtime"); + ceObjects->Add(rocTtime); - rocQtime = (TObjArray*)ceObjects->FindObject("rocQtime"); - if ( !rocQtime ) { - rocQtime = new TObjArray(nSectors); - rocQtime->SetName("rocQtime"); - ceObjects->Add(rocQtime); + rocQtime = new TObjArray(nSectors); + rocQtime->SetName("rocQtime"); + ceObjects->Add(rocQtime); + +// Temperature maps + + if (fTemp) { + AliTPCSensorTempArray *tempMap = new AliTPCSensorTempArray(*fTemp); + tempMap->SetNameTitle("TempMap","TempMap"); + ceObjects->Add(tempMap); } +// Pressure maps + + if (fPressure) { + AliDCSSensor *sensor=0, *sensorCopy=0; + for (Int_t isensor=0; isensorGetSensor(kPressureSensorNames[isensor]); + if (sensor) { + sensorCopy = new AliDCSSensor(*sensor); + sensorCopy->SetNameTitle(kPressureSensorNames[isensor],kPressureSensorNames[isensor]); + ceObjects->Add(sensorCopy); + } + } + } UInt_t result=0; @@ -697,6 +986,11 @@ UInt_t AliTPCPreprocessor::ExtractCE(Int_t sourceFXS) AliTPCCalibCE *calCE; f->GetObject("tpcCalibCE",calCE); + if (!calCE) { + Log ("No valid calibCE object."); + result=2; + break; + } // replace entries for the sectors available in the present file for (Int_t sector=0; sectorMakeGraphTimeCE(sector,0,2); // T time graph if ( grT ) rocTtime->AddAt(grT,sector); TGraph *grQ=calCE->MakeGraphTimeCE(sector,0,3); // Q time graph - if ( grQ ) rocTtime->AddAt(grQ,sector); + if ( grQ ) rocQtime->AddAt(grQ,sector); } + + TGraph *grT=calCE->MakeGraphTimeCE(-1,0,2); // A side average + if ( grT ) { + rocTtime->AddAt(grT,nSectors); + } else { + result=10; + } + grT=calCE->MakeGraphTimeCE(-2,0,2); // C side average + if ( grT ) { + rocTtime->AddAt(grT,nSectors+1); + } else { + result=10; + } delete calCE; f->Close(); } @@ -722,15 +1029,360 @@ UInt_t AliTPCPreprocessor::ExtractCE(Int_t sourceFXS) AliCDBMetaData metaData; metaData.SetBeamPeriod(0); metaData.SetResponsible("Haavard Helstrup"); + metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH); metaData.SetComment("Preprocessor AliTPC data base entries."); - Bool_t storeOK = Store("Calib", "CE", ceObjects, &metaData, 0, kTRUE); - if ( !storeOK ) ++result; - + if ( result == 0 ) { + Bool_t storeOK = Store("Calib", "CE", ceObjects, &metaData, 0, kTRUE); + if ( !storeOK ) ++result; + } else { + Log ("Warning: Average time graphs not available - no OCDB entry written"); + } } else { - Log ("Error: no entries!"); + Log ("Error: no CE entries available from FXS!"); result = 1; } + ceObjects->Delete(); + delete ceObjects; + + return result; +} +//______________________________________________________________________________________________ + +UInt_t AliTPCPreprocessor::ExtractQA(Int_t sourceFXS) +{ + // + // Read Quality Assurance file from file exchage server + // + + UInt_t result=0; + + TList* list = GetFileSources(sourceFXS,"QA"); + + if (list && list->GetEntries()>0) { + +// only one QA objetc should be available! + + AliTPCdataQA *calQA; + + UInt_t nentries = list->GetEntries(); + UInt_t index=0; + if ( nentries > 1) Log ( "More than one QA entry. First one processed"); + TObjString* fileNameEntry = (TObjString*) list->At(index); + if (fileNameEntry!=NULL) { + TString fileName = GetFile(sourceFXS, "QA", + fileNameEntry->GetString().Data()); + TFile *f = TFile::Open(fileName); + if (!f) { + Log ("Error opening QA file."); + result =2; + } else { + f->GetObject("tpcCalibQA",calQA); + if ( calQA ) { +// +// Store updated pedestal entry to OCDB +// + AliCDBMetaData metaData; + metaData.SetBeamPeriod(0); + metaData.SetResponsible("Haavard Helstrup"); + metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH); + metaData.SetComment("Preprocessor AliTPC data base entries."); + + Bool_t storeOK = Store("Calib", "QA", calQA, &metaData, 0, kFALSE); + if ( !storeOK ) ++result; + + delete calQA; + } + } + } else { + Log ("Error: no QA files on FXS!"); + result = 2; + } + } else { + Log ("Error: no QA entries in FXS list!"); + result = 1; + } + return result; +} + +//______________________________________________________________________________________________ + + +UInt_t AliTPCPreprocessor::ExtractAltro(Int_t sourceFXS, TMap* dcsMap) +{ + // + // Read Altro configuration file from file exchage server + // Keep original entry from OCDB in case no new pulser calibration is available + // + TObjArray *altroObjects=0; + AliTPCCalPad *acqStart=0; + AliTPCCalPad *zsThr=0; + AliTPCCalPad *acqStop=0; + AliTPCCalPad *FPED=0; + AliTPCCalPad *masked=0; + AliTPCCalPad *k1=0, *k2=0, *k3=0; + AliTPCCalPad *l1=0, *l2=0, *l3=0; + TMap *mapRCUconfig=0; + + AliCDBEntry* entry = GetFromOCDB("Calib", "AltroConfig"); + if (entry) altroObjects = (TObjArray*)entry->GetObject(); + if ( altroObjects==NULL ) { + Log("AliTPCPreprocsessor: No previous TPC altro calibration entry available.\n"); + altroObjects = new TObjArray; + } + + acqStart = (AliTPCCalPad*)altroObjects->FindObject("AcqStart"); + if ( !acqStart ) { + acqStart = new AliTPCCalPad("AcqStart","AcqStart"); + altroObjects->Add(acqStart); + } + zsThr = (AliTPCCalPad*)altroObjects->FindObject("ZsThr"); + if ( !zsThr ) { + zsThr = new AliTPCCalPad("ZsThr","ZsThr"); + altroObjects->Add(zsThr); + } + FPED = (AliTPCCalPad*)altroObjects->FindObject("FPED"); + if ( !FPED ) { + FPED = new AliTPCCalPad("FPED","FPED"); + altroObjects->Add(FPED); + } + acqStop = (AliTPCCalPad*)altroObjects->FindObject("AcqStop"); + if ( !acqStop ) { + acqStop = new AliTPCCalPad("AcqStop","AcqStop"); + altroObjects->Add(acqStop); + } + masked = (AliTPCCalPad*)altroObjects->FindObject("Masked"); + if ( !masked ) { + masked = new AliTPCCalPad("Masked","Masked"); + altroObjects->Add(masked); + } + k1 = (AliTPCCalPad*)altroObjects->FindObject("K1"); + if ( !k1 ) { + k1 = new AliTPCCalPad("K1","K1"); + altroObjects->Add(k1); + } + k2 = (AliTPCCalPad*)altroObjects->FindObject("K2"); + if ( !k2 ) { + k2 = new AliTPCCalPad("K2","K2"); + altroObjects->Add(k2); + } + k3 = (AliTPCCalPad*)altroObjects->FindObject("K3"); + if ( !k3 ) { + k3 = new AliTPCCalPad("K3","K3"); + altroObjects->Add(k3); + } + l1 = (AliTPCCalPad*)altroObjects->FindObject("L1"); + if ( !l1 ) { + l1 = new AliTPCCalPad("L1","L1"); + altroObjects->Add(l1); + } + l2 = (AliTPCCalPad*)altroObjects->FindObject("L2"); + if ( !l2 ) { + l2 = new AliTPCCalPad("L2","L2"); + altroObjects->Add(l2); + } + l3 = (AliTPCCalPad*)altroObjects->FindObject("L3"); + if ( !l3 ) { + l3 = new AliTPCCalPad("L3","L3"); + altroObjects->Add(l3); + } + mapRCUconfig = (TMap*)altroObjects->FindObject("RCUconfig"); + if (!mapRCUconfig) { + mapRCUconfig = new TMap(); + mapRCUconfig->SetName("RCUconfig"); + altroObjects->Add(mapRCUconfig); + } + + + UInt_t result=0; + TString idFXS[2]={"AltroConfigA","AltroConfigC"}; + + Int_t nSectors = fROC->GetNSectors(); + Bool_t changed=false; + if (altroObjects == 0 ) altroObjects = new TObjArray; + +// extract list of active DDLs + + Bool_t found; + TString arrDDL(kNumDDL); + arrDDL.Append('x',kNumDDL); + for ( Int_t iDDL = 0; iDDLFindObject(stringID.Data()); + found = false; + if ( pair ) { + TObjArray *valueSet=(TObjArray*)pair->Value(); + if ( valueSet) { + AliDCSValue *val = (AliDCSValue*)valueSet->At(0); + if (val) { + found = val->GetBool(); + if (found){ + arrDDL[iDDL] = '1'; + } else { + arrDDL[iDDL] = '0'; + } + } + } + } + } + TObjString *ddlArray = new TObjString; + ddlArray->SetString(arrDDL); + TMap *activeDDL = new TMap; + activeDDL->SetName("DDLArray"); + TObjString *key = new TObjString("DDLArray"); + activeDDL->Add(key,ddlArray); + altroObjects->Add(activeDDL); + changed=true; + + +// extract Altro configuration files + + for ( Int_t id=0; id<2; id++) { + TList* list = GetFileSources(sourceFXS,idFXS[id].Data()); + + if (list && list->GetEntries()>0) { + +// loop through all files from LDCs + + UInt_t index = 0; + while (list->At(index)!=NULL) { + TObjString* fileNameEntry = (TObjString*) list->At(index); + if (fileNameEntry!=NULL) { + TString fileName = GetFile(sourceFXS, idFXS[id].Data(), + fileNameEntry->GetString().Data()); + TFile *f = TFile::Open(fileName); + if (!f) { + char message[40]; + sprintf(message,"Error opening Altro configuration file, id = %d",id); + Log (message); + result =2; + break; + } + TObjArray *altroFXS; + f->GetObject("AltroConfig",altroFXS); + if ( !altroFXS ) { + Log ("No Altro configuration object in file."); + result = 2; + break; + } + + // replace entries for the sectors available in the present file + AliTPCCalPad *acqStartFXS=(AliTPCCalPad*)altroFXS->FindObject("AcqStart"); + AliTPCCalPad *zsThrFXS=(AliTPCCalPad*)altroFXS->FindObject("ZsThr"); + AliTPCCalPad *acqStopFXS=(AliTPCCalPad*)altroFXS->FindObject("AcqStop"); + AliTPCCalPad *FPEDFXS=(AliTPCCalPad*)altroFXS->FindObject("FPED"); + AliTPCCalPad *maskedFXS=(AliTPCCalPad*)altroFXS->FindObject("Masked"); + AliTPCCalPad *k1FXS=(AliTPCCalPad*)altroFXS->FindObject("K1"); + AliTPCCalPad *k2FXS=(AliTPCCalPad*)altroFXS->FindObject("K2"); + AliTPCCalPad *k3FXS=(AliTPCCalPad*)altroFXS->FindObject("K3"); + AliTPCCalPad *l1FXS=(AliTPCCalPad*)altroFXS->FindObject("L1"); + AliTPCCalPad *l2FXS=(AliTPCCalPad*)altroFXS->FindObject("L2"); + AliTPCCalPad *l3FXS=(AliTPCCalPad*)altroFXS->FindObject("L3"); + TMap *mapRCUconfigFXS = (TMap*)altroFXS->FindObject("RCUconfig"); + TIterator *mapFXSiter = mapRCUconfigFXS->MakeIterator(); + + changed=true; + for (Int_t sector=0; sectorGetCalROC(sector); + if ( rocAcqStart ) acqStart->SetCalROC(rocAcqStart,sector); + } + if (zsThrFXS ) { + AliTPCCalROC *rocZsThr=zsThrFXS->GetCalROC(sector); + if ( rocZsThr ) zsThr->SetCalROC(rocZsThr,sector); + } + if (acqStopFXS) { + AliTPCCalROC *rocAcqStop=acqStopFXS->GetCalROC(sector); + if ( rocAcqStop ) acqStop->SetCalROC(rocAcqStop,sector); + } + if (FPEDFXS ) { + AliTPCCalROC *rocFPED=FPEDFXS->GetCalROC(sector); + if ( rocFPED ) FPED->SetCalROC(rocFPED,sector); + } + if (maskedFXS) { + AliTPCCalROC *rocMasked=maskedFXS->GetCalROC(sector); + if ( rocMasked ) masked->SetCalROC(rocMasked,sector); + } + if (k1FXS) { + AliTPCCalROC *rocK1=k1FXS->GetCalROC(sector); + if ( rocK1 ) k1->SetCalROC(rocK1,sector); + } + if (k2FXS) { + AliTPCCalROC *rocK2=k2FXS->GetCalROC(sector); + if ( rocK2 ) k2->SetCalROC(rocK2,sector); + } + if (k3FXS) { + AliTPCCalROC *rocK3=k3FXS->GetCalROC(sector); + if ( rocK3 ) k3->SetCalROC(rocK3,sector); + } + if (l1FXS) { + AliTPCCalROC *rocL1=l1FXS->GetCalROC(sector); + if ( rocL1 ) l1->SetCalROC(rocL1,sector); + } + if (l2FXS) { + AliTPCCalROC *rocL2=l2FXS->GetCalROC(sector); + if ( rocL2 ) l2->SetCalROC(rocL2,sector); + } + if (l3FXS) { + AliTPCCalROC *rocL3=l3FXS->GetCalROC(sector); + if ( rocL3 ) l3->SetCalROC(rocL3,sector); + } + } + if (mapRCUconfigFXS) { + Int_t mapEntries = mapRCUconfigFXS->GetEntries(); + TObjString* keyFXS; + TVectorF* vecFXS; + TVectorF* vec; // nSectors = 72 (total number of inner/outer sectors) + for (Int_t i=0; iNext(); + vecFXS=(TVectorF*)mapRCUconfigFXS->GetValue(keyFXS); + vec=(TVectorF*)mapRCUconfig->GetValue(keyFXS); + if (!vec) { + vec = new TVectorF(3*nSectors); + *vec = -1; + mapRCUconfig->Add(keyFXS,vec); + } + if (vec->GetNoElements() != 3*nSectors ) { + vec->ResizeTo(3*nSectors); + } + if (id==0) { // A side + vec->SetSub(0,vecFXS->GetSub(0,nSectors/2-1)); + vec->SetSub(nSectors,vecFXS->GetSub(nSectors,2*nSectors-1)); + } else { // C side + vec->SetSub(nSectors/2,vecFXS->GetSub(nSectors/2,nSectors-1)); + vec->SetSub(2*nSectors,vecFXS->GetSub(2*nSectors,3*nSectors-1)); + } + } + } + delete altroFXS; + f->Close(); + } + ++index; + } // while(list) + } else { + Log ("Error: no entries in AltroConfig file list!"); + result = 1; + } + + } // for - id +// +// Store updated pedestal entry to OCDB +// + if (changed) { + AliCDBMetaData metaData; + metaData.SetBeamPeriod(0); + metaData.SetResponsible("Haavard Helstrup"); + metaData.SetAliRootVersion(ALIROOT_SVN_BRANCH); + metaData.SetComment("Preprocessor AliTPC data base entries."); + + Bool_t storeOK = Store("Calib", "AltroConfig", altroObjects, &metaData, 0, kFALSE); + if ( !storeOK ) ++result; + } + + altroObjects->Delete(); + delete altroObjects; + return result; }