**************************************************************************/
/* $Id$ */
-/* History of cvs commits:
- *
- * $Log$
- *
-*/
+
///////////////////////////////////////////////////////////////////////////////
// EMCAL Preprocessor class. It runs by Shuttle at the end of the run,
-// calculates calibration coefficients and dead/bad channels
-// to be posted in OCDB
+// calculates stuff to be posted in OCDB
//
// Author: Boris Polichtchouk, 4 October 2006
// Adapted for EMCAL by Gustavo Conesa Balbastre, October 2006
+// Updated by David Silvermyr May 2008, based on TPC code
///////////////////////////////////////////////////////////////////////////////
+//Root
+#include "TFile.h"
+#include "TTree.h"
+#include "TEnv.h"
+#include "TParameter.h"
+
+#include <TTimeStamp.h>
+
+//AliRoot
+#include "AliShuttleInterface.h"
#include "AliEMCALPreprocessor.h"
#include "AliLog.h"
+#include "AliDCSValue.h"
#include "AliCDBMetaData.h"
-#include "AliEMCALCalibData.h"
-#include "TFile.h"
-#include "TH1.h"
-#include "TMap.h"
-#include "TRandom.h"
+#include "AliEMCALTriggerTRUDCSConfig.h"
+#include "AliEMCALTriggerSTUDCSConfig.h"
+#include "AliEMCALTriggerDCSConfig.h"
+#include "AliCaloCalibPedestal.h"
+#include "AliCaloCalibSignal.h"
+#include "AliEMCALSensorTempArray.h"
-ClassImp(AliEMCALPreprocessor)
+const Int_t kValCutTemp = 100; // discard temperatures > 100 degrees
+const Int_t kDiffCutTemp = 5; // discard temperature differences > 5 degrees
+const TString kPedestalRunType = "PEDESTAL"; // pedestal run identifier
+const TString kPhysicsRunType = "PHYSICS"; // physics run identifier
+const TString kStandAloneRunType = "STANDALONE_BC"; // standalone run identifier
+const TString kAmandaTemp = "PT_%02d.Temperature"; // Amanda string for temperature entries
+//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 during commissioning
+
+const TString kMetaResponsible = "David Silvermyr";
+//legacy comments and return codes from TPC
+const TString kMetaComment = "Preprocessor AliEMCAL data base entries.";
+const int kReturnCodeNoInfo = 9;
+const int kReturnCodeNoObject = 2;
+const int kReturnCodeNoEntries = 1;
+const int kNTRU = 62; // EMCAL (10 SM + 2/3) + (DCAL 10 SM) * 3 TRU per SM
+
+ClassImp(AliEMCALPreprocessor)
+
//_______________________________________________________________________________________
AliEMCALPreprocessor::AliEMCALPreprocessor() :
-AliPreprocessor("EMCAL",0)
+ AliPreprocessor("EMC",0),
+ fConfEnv(0),
+ fTemp(0),
+ fConfigOK(kTRUE)
{
//default constructor
}
//_______________________________________________________________________________________
-AliEMCALPreprocessor::AliEMCALPreprocessor(const char* detector, AliShuttleInterface* shuttle):
-AliPreprocessor(detector,shuttle)
+AliEMCALPreprocessor::AliEMCALPreprocessor(AliShuttleInterface* shuttle):
+ AliPreprocessor("EMC",shuttle),
+ fConfEnv(0),
+ fTemp(0),
+ fConfigOK(kTRUE)
{
- // Constructor
+ // Constructor AddRunType(kPedestalRunType);
+
+ // define run types to be processed
+ AddRunType(kPedestalRunType);
+ AddRunType(kPhysicsRunType);
}
-//_______________________________________________________________________________________
-UInt_t AliEMCALPreprocessor::Process(TMap* /*valueSet*/)
+//______________________________________________________________________________________________
+AliEMCALPreprocessor::AliEMCALPreprocessor(const AliEMCALPreprocessor& ) :
+ AliPreprocessor("EMCAL",0),
+ fConfEnv(0), fTemp(0), fConfigOK(kTRUE)
{
- // process data retrieved by the Shuttle
+ Fatal("AliEMCALPreprocessor", "copy constructor not implemented");
+}
- // The fileName with the histograms which have been produced by
- // AliEMCALCalibHistoProducer.
- // It is a responsibility of the SHUTTLE framework to form the fileName
+// assignment operator; use copy ctor to make life easy.
+//______________________________________________________________________________________________
+AliEMCALPreprocessor& AliEMCALPreprocessor::operator = (const AliEMCALPreprocessor &source )
+{
+ // assignment operator; use copy ctor
+ if (&source == this) return *this;
+
+ new (this) AliEMCALPreprocessor(source);
+ return *this;
+}
- const char* fileName = GetFile(kDAQ, "AMPLITUDES", "GDC");
- AliInfo(Form("Got filename: %s",fileName));
+//____________________________________________________________________________
+AliEMCALPreprocessor::~AliEMCALPreprocessor()
+{
+ // destructor
+ if (fTemp) delete fTemp;
+}
- TFile f(fileName);
+//______________________________________________________________________________________________
+void AliEMCALPreprocessor::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);
+
+ AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run,
+ TTimeStamp((time_t)startTime,0).AsString(),
+ TTimeStamp((time_t)endTime,0).AsString()));
+
+ // Preprocessor configuration
+ AliCDBEntry* entry = GetFromOCDB("Config", "Preprocessor");
+ if (entry) fConfEnv = (TEnv*) entry->GetObject();
+ if ( fConfEnv==0 ) {
+ Log("AliEMCALPreprocessor: Preprocessor Config OCDB entry missing.\n");
+ fConfigOK = kFALSE;
+ return;
+ }
+
+ // Temperature sensors
+ TTree *confTree = 0;
+
+ TString tempConf = fConfEnv->GetValue("Temperature","ON");
+ tempConf.ToUpper();
+ if (tempConf != "OFF" ) {
+ entry = GetFromOCDB("Config", "Temperature");
+ if (entry) confTree = (TTree*) entry->GetObject();
+ if ( confTree==0 ) {
+ Log("AliEMCALPreprocessor: Temperature Config OCDB entry missing.\n");
+ fConfigOK = kFALSE;
+ return;
+ }
+ fTemp = new AliEMCALSensorTempArray(startTimeLocal, fEndTime, confTree, kAmandaTemp);
+ fTemp->SetValCut(kValCutTemp);
+ fTemp->SetDiffCut(kDiffCutTemp);
+ }
+
+ return;
+}
- if(!f.IsOpen()) {
- AliInfo(Form("File %s is not opened, something goes wrong!",fileName));
- return 0;
+//______________________________________________________________________________________________
+UInt_t AliEMCALPreprocessor::Process(TMap* dcsAliasMap)
+{
+ // Fills data into EMCAL calibrations objects
+ // Amanda servers provide information directly through dcsAliasMap
+
+ if (!fConfigOK) return kReturnCodeNoInfo;
+ UInt_t result = 0;
+ TObjArray *resultArray = new TObjArray();
+ TString errorHandling = fConfEnv->GetValue("ErrorHandling","ON");
+ errorHandling.ToUpper();
+ TObject * status;
+
+ UInt_t dcsResult=0;
+ if (errorHandling == "OFF" ) {
+ if (!dcsAliasMap) dcsResult = kReturnCodeNoEntries;
+ else if (dcsAliasMap->GetEntries() == 0 ) dcsResult = kReturnCodeNoEntries;
+ status = new TParameter<int>("dcsResult",dcsResult);
+ resultArray->Add(status);
+ }
+ else {
+ if (!dcsAliasMap) return kReturnCodeNoInfo;
+ else if (dcsAliasMap->GetEntries() == 0 ) return kReturnCodeNoInfo;
+ }
+
+
+ TString runType = GetRunType();
+
+ // Temperature sensors are processed by AliEMCALCalTemp
+ TString tempConf = fConfEnv->GetValue("Temperature","ON");
+ tempConf.ToUpper();
+ if (tempConf != "OFF" && dcsAliasMap ) {
+ UInt_t tempResult = MapTemperature(dcsAliasMap);
+ result=tempResult;
+ status = new TParameter<int>("tempResult",tempResult);
+ resultArray->Add(status);
}
+ // Trigger configuration processing: only for Physics runs
+ TString triggerConf = fConfEnv->GetValue("Trigger","ON");
+ triggerConf.ToUpper();
+ if( runType == kPhysicsRunType ) {
+ if (triggerConf != "OFF" && dcsAliasMap ) {
+ UInt_t triggerResult = MapTriggerConfig(dcsAliasMap);
+ result+=triggerResult;
+ status = new TParameter<int>("triggerResult",triggerResult);
+ resultArray->Add(status);
+ }
+ }
+
+ // Other calibration information will be retrieved through FXS files
+ // examples:
+ // TList* fileSourcesDAQ = GetFile(AliShuttleInterface::kDAQ, "pedestals");
+ // const char* fileNamePed = GetFile(AliShuttleInterface::kDAQ, "pedestals", "LDC1");
+ //
+ // TList* fileSourcesHLT = GetFile(AliShuttleInterface::kHLT, "calib");
+ // const char* fileNameHLT = GetFile(AliShuttleInterface::kHLT, "calib", "LDC1");
+
+ // PEDESTAL ENTRIES:
+
+ if ( runType == kPedestalRunType ) {
+ Int_t numSources = 1;
+ Int_t pedestalSource[2] = {AliShuttleInterface::kDAQ, AliShuttleInterface::kHLT} ;
+ TString source = fConfEnv->GetValue("Pedestal","DAQ");
+ source.ToUpper();
+ if (source != "OFF" ) {
+ if ( source == "HLT") pedestalSource[0] = AliShuttleInterface::kHLT;
+ if (!GetHLTStatus()) pedestalSource[0] = AliShuttleInterface::kDAQ;
+ if (source == "HLTDAQ" ) {
+ numSources=2;
+ pedestalSource[0] = AliShuttleInterface::kHLT;
+ pedestalSource[1] = AliShuttleInterface::kDAQ;
+ }
+ if (source == "DAQHLT" ) numSources=2;
+ UInt_t pedestalResult=0;
+ for (Int_t i=0; i<numSources; i++ ) {
+ pedestalResult = ExtractPedestals(pedestalSource[i]);
+ if ( pedestalResult == 0 ) break;
+ }
+ result += pedestalResult;
+ status = new TParameter<int>("pedestalResult",pedestalResult);
+ resultArray->Add(status);
+ }
+ }
+
+ // SIGNAL/LED ENTRIES:
+ if( runType == kPhysicsRunType ) {
+ Int_t numSources = 1;
+ Int_t signalSource[2] = {AliShuttleInterface::kDAQ,AliShuttleInterface::kHLT} ;
+ TString source = fConfEnv->GetValue("Signal","DAQ");
+ source.ToUpper();
+ if ( source != "OFF") {
+ if ( source == "HLT") signalSource[0] = AliShuttleInterface::kHLT;
+ if (!GetHLTStatus()) signalSource[0] = AliShuttleInterface::kDAQ;
+ if (source == "HLTDAQ" ) {
+ numSources=2;
+ signalSource[0] = AliShuttleInterface::kHLT;
+ signalSource[1] = AliShuttleInterface::kDAQ;
+ }
+ if (source == "DAQHLT" ) numSources=2;
+ UInt_t signalResult=0;
+ for (Int_t i=0; i<numSources; i++ ) {
+ signalResult = ExtractSignal(signalSource[i]);
+ if ( signalResult == 0 ) break;
+ }
+ result += signalResult;
+ status = new TParameter<int>("signalResult",signalResult);
+ resultArray->Add(status);
+ }
+ }
+
+
+ // overall status at the end
+ if (errorHandling == "OFF" ) {
+ AliCDBMetaData metaData;
+ metaData.SetBeamPeriod(0);
+ metaData.SetResponsible(kMetaResponsible);
+ metaData.SetComment("Preprocessor AliEMCAL status.");
+ Bool_t storeOK = Store("Calib", "PreprocStatus", resultArray, &metaData, 0, kFALSE);
+ resultArray->Delete();
+ result = 0;
+ if ( !storeOK ) result=1;
+ return result;
+ }
+ else {
+ return result;
+ }
+
+}
+//______________________________________________________________________________________________
+UInt_t AliEMCALPreprocessor::MapTemperature(TMap* dcsAliasMap)
+{ // extract DCS temperature maps. Perform fits to save space
+ UInt_t result=0;
+
+ TMap *map = fTemp->ExtractDCS(dcsAliasMap);
+ if (map) {
+ fTemp->MakeSplineFit(map);
+ Double_t fitFraction = 1.0*fTemp->NumFits()/fTemp->NumSensors();
+ if (fitFraction > kFitFraction ) {
+ AliInfo(Form("Temperature values extracted, fits performed.\n"));
+ }
+ else {
+ Log ("Too few temperature maps fitted. \n");
+ result = kReturnCodeNoInfo;
+ }
+ }
+ else {
+ Log("No temperature map extracted. \n");
+ result = kReturnCodeNoInfo;
+ }
+ delete map;
+ // Now store the final CDB file
+
+ if ( result == 0 ) { // some info was found
+ AliCDBMetaData metaData;
+ metaData.SetBeamPeriod(0);
+ metaData.SetResponsible(kMetaResponsible);
+ metaData.SetComment(kMetaComment);
+
+ Bool_t storeOK = Store("Calib", "Temperature", fTemp, &metaData, 0, kFALSE);
+ if ( !storeOK ) result=1;
+ AliInfo(Form("Temperature info stored. result %d\n", result));
+ }
+
+ return result;
+}
- const Int_t nMod=12; // 1:5 modules
- const Int_t nCol=48; //1:56 columns in each module
- Int_t nRow=24; //1:64 rows in each module
- const Int_t nRowHalfSM = 12; //Supermodules 11 and 12 are half supermodules
-
- Double_t coeff;
- char hnam[80];
- TH1F* histo=0;
-
- // Generate name of the reference histogram
- TString refHistoName = "mod";
- refHistoName += gRandom->Integer(nMod)+1;
- refHistoName += "col";
- refHistoName += gRandom->Integer(nCol)+1;
- refHistoName += "row";
- refHistoName += gRandom->Integer(nRow)+1;
- TH1F* hRef = (TH1F*)f.Get(refHistoName);
-
- // If the reference histogram does not exist or has too little statistics,
- // it is better to give up preprocessing
- if(!hRef || hRef->GetEntries()<2) {
- AliInfo(Form("Cannot get reference histogram %s",refHistoName.Data()));
- return 0;
+//______________________________________________________________________________________________
+UInt_t AliEMCALPreprocessor::MapTriggerConfig(TMap* dcsAliasMap)
+{ // extract DCS trigger info
+
+
+ AliInfo("Print DCS alias map content");
+ dcsAliasMap->Print();
+
+ AliInfo(Form("Get TRU info from DCS DPs.\n"));
+ Int_t i, iTRU;
+ const Int_t bufsize = 1000;
+ char buf[bufsize];
+
+ AliDCSValue *dcsVal;
+ TObjArray *arrL0ALGSEL, *arrPEAKFINDER, *arrGLOBALTHRESH, *arrCOSMTHRESH;
+ TObjArray *arrMASK[6];
+
+ TObjArray *arrSTUG[3][2], *arrSTUJ[3][2];
+ TObjArray *arrSTUD, *arrSTUR, *arrSTUF;
+
+ // overall object to hold STU and DCS config info
+ // DS comment: for now only holds TRU info, i.e. only partially filled
+ // (STU info only in raw data header; unfortunately not also picked up via DCS DPs)
+ AliEMCALTriggerDCSConfig *trigConfig = new AliEMCALTriggerDCSConfig();
+ // allocate space for TRU objects
+ TClonesArray *truArr = new TClonesArray("AliEMCALTriggerTRUDCSConfig", kNTRU);
+ for( iTRU = 0; iTRU < kNTRU; iTRU++){
+ new((*truArr)[iTRU]) AliEMCALTriggerTRUDCSConfig();
}
+ trigConfig->SetTRUArr(truArr);
+
+ AliEMCALTriggerSTUDCSConfig *stuConfig = new AliEMCALTriggerSTUDCSConfig();
+
+ // loop through all TRUs
+ bool debug = true; // debug flag for AliInfo printouts for each TRU
+ for( iTRU = 0; iTRU < kNTRU; iTRU++){
+ if (debug) AliInfo( Form("iTRU %d \n", iTRU) );
+ // get the shuttled values
+ snprintf( buf, bufsize, "EMC_TRU%02d_L0ALGSEL", iTRU );
+ arrL0ALGSEL = (TObjArray*) dcsAliasMap->GetValue( buf );
+ snprintf( buf, bufsize, "EMC_TRU%02d_PEAKFINDER", iTRU );
+ arrPEAKFINDER = (TObjArray*) dcsAliasMap->GetValue( buf );
+ snprintf( buf, bufsize, "EMC_TRU%02d_GLOBALTHRESH", iTRU );
+ arrGLOBALTHRESH = (TObjArray*) dcsAliasMap->GetValue( buf );
+ snprintf( buf, bufsize, "EMC_TRU%02d_COSMTHRESH", iTRU );
+ arrCOSMTHRESH = (TObjArray*) dcsAliasMap->GetValue( buf );
+
+ for( i = 0; i < 6; i++ ){
+ snprintf( buf, bufsize, "EMC_TRU%02d_MASK%d", iTRU, i );
+ arrMASK[i] = (TObjArray*) dcsAliasMap->GetValue( buf );
+ }
+
+ // fill the objects
+ AliEMCALTriggerTRUDCSConfig* truConfig = trigConfig->GetTRUDCSConfig(iTRU);
+ if( ! truConfig ){
+ AliWarning( Form("EMC TRU%02d config not retrieved!\n", iTRU ));
+ continue;
+ }
- AliEMCALCalibData calibData;
- Double_t refMean=hRef->GetMean();
- // Calculates relative calibration coefficients for all non-zero channels
-
- for(Int_t mod=0; mod<nMod; mod++) {
- if(mod > 10) nRow = nRowHalfSM ;
- for(Int_t col=0; col<nCol; col++) {
- for(Int_t row=0; row<nRow; row++) {
- sprintf(hnam,"mod%dcol%drow%d",mod,col,row);
- histo = (TH1F*)f.Get(hnam);
- //TODO: dead channels exclusion!
- if(histo) {
- coeff = histo->GetMean()/refMean;
- calibData.SetADCchannel(mod+1,col+1,row+1,1./coeff);
- AliInfo(Form("mod %d col %d row %d coeff %f\n",mod,col,row,coeff));
- }
- else
- calibData.SetADCchannel(mod+1,col+1,row+1,-111);
+ // get last entries. fill the TRU object
+ if( ! arrL0ALGSEL ){
+ AliWarning( Form("EMC_TRU%02d_L0ALGSEL alias not found!\n", iTRU ));
+ }
+ else{
+ if (debug) AliInfo( Form("arrL0ALGSEL has %d entries \n", arrL0ALGSEL->GetEntries()) );
+ if ( arrL0ALGSEL->GetEntries() > 0 ) {
+ dcsVal = (AliDCSValue *) arrL0ALGSEL->At( arrL0ALGSEL->GetEntries() - 1 );
+ if (dcsVal) {
+ truConfig->SetL0SEL( dcsVal->GetUInt() );
+ if (debug) AliInfo( Form("saving value: %u\n", dcsVal->GetUInt()) );
+ }
}
+ else
+ AliWarning( Form("EMC_TRU%02d_L0ALGSEL has no entries!\n", iTRU ));
}
- }
+ if( ! arrPEAKFINDER ){
+ AliWarning( Form("EMC_TRU%02d_PEAKFINDER alias not found!\n", iTRU ));
+ }
+ else{
+ if (debug) AliInfo( Form("arrPEAKFINDER has %d entries \n", arrPEAKFINDER->GetEntries()) );
+ if ( arrPEAKFINDER->GetEntries() > 0 ) {
+ dcsVal = (AliDCSValue *) arrPEAKFINDER->At( arrPEAKFINDER->GetEntries() - 1 );
+ if (dcsVal){
+ truConfig->SetSELPF( dcsVal->GetUInt() );
+ if (debug) AliInfo( Form("saving value: %u\n", dcsVal->GetUInt()) );
+ }
+ }
+ else
+ AliWarning( Form("EMC_TRU%02d_PEAKFINDER has no entries!\n", iTRU ));
+ }
+
+ if( ! arrGLOBALTHRESH ){
+ AliWarning( Form("EMC_TRU%02d_GLOBALTHRESH alias not found!\n", iTRU ));
+ }
+ else{
+ if (debug) AliInfo( Form("arrGLOBALTHRESH has %d entries \n", arrGLOBALTHRESH->GetEntries()) );
+ if ( arrGLOBALTHRESH->GetEntries() > 0 ) {
+ dcsVal = (AliDCSValue *) arrGLOBALTHRESH->At( arrGLOBALTHRESH->GetEntries() - 1 );
+ if (dcsVal){
+ truConfig->SetGTHRL0( dcsVal->GetUInt() );
+ if (debug) AliInfo( Form("saving value: %u\n", dcsVal->GetUInt()) );
+ }
+ }
+ else
+ AliWarning( Form("EMC_TRU%02d_GLOBALTHRESH has no entries!\n", iTRU ));
+ }
+
+ if( ! arrCOSMTHRESH ){
+ AliWarning( Form("EMC_TRU%02d_COSMTHRESH alias not found!\n", iTRU ));
+ }
+ else{
+ if (debug) AliInfo( Form("arrCOSMTHRESH has %d entries \n", arrCOSMTHRESH->GetEntries()) );
+ if ( arrCOSMTHRESH->GetEntries() > 0 ) {
+ dcsVal = (AliDCSValue *) arrCOSMTHRESH->At( arrCOSMTHRESH->GetEntries() - 1 );
+ if (dcsVal){
+ truConfig->SetL0COSM( dcsVal->GetUInt() );
+ if (debug) AliInfo( Form("saving value: %u\n", dcsVal->GetUInt()) );
+ }
+ }
+ else
+ AliWarning( Form("EMC_TRU%02d_COSMTHRESH has no entries!\n", iTRU ));
+ }
+
+ for( i = 0; i < 6; i++ ){
+ if( ! arrMASK[i] ){
+ AliWarning( Form("EMC_TRU%02d_MASK%d alias not found!\n", iTRU, i ));
+ }
+ else{
+ if (debug) AliInfo( Form("arrMASK[%d] has %d entries \n", i, arrMASK[i]->GetEntries()) );
+ if ( arrMASK[i]->GetEntries() > 0 ) {
+ dcsVal = (AliDCSValue *) arrMASK[i]->At( arrMASK[i]->GetEntries() - 1 );
+ if (dcsVal){
+ truConfig->SetMaskReg( dcsVal->GetUInt(), i );
+ if (debug) AliInfo( Form("saving value: %u\n", dcsVal->GetUInt()) );
+ }
+ }
+ else
+ AliWarning( Form("EMC_TRU%02d_MASK%d has no entries!\n", iTRU, i ));
+ }
+ }
+
+ } // TRUs
+ AliInfo(Form("TRU info retrieved.\n"));
+
+ // STU
+ for (i = 0; i < 3; i++) {
+ for (int j = 0; j < 2; j++) {
+ arrSTUG[i][j] = (TObjArray*)dcsAliasMap->GetValue(Form("EMC_STU_G%c%d", i + 65, j));
+ arrSTUJ[i][j] = (TObjArray*)dcsAliasMap->GetValue(Form("EMC_STU_J%c%d", i + 65, j));
+
+ if (!arrSTUG[i][j]) {
+ AliWarning(Form("EMC_STU_G%c%d alias not found!", i + 65, j));
+ } else {
+ if (debug) AliInfo( Form("EMC_STU_G%c%d has %d entries", i + 65, j, arrSTUG[i][j]->GetEntries()));
+ if (arrSTUG[i][j]->GetEntries() > 0) {
+ dcsVal = (AliDCSValue*)arrSTUG[i][j]->At(arrSTUG[i][j]->GetEntries() - 1);
+ if (dcsVal) {
+ stuConfig->SetG(i, j, dcsVal->GetInt());
+ if (debug) AliInfo(Form("saving value: %u", dcsVal->GetInt()));
+ }
+ } else
+ AliWarning(Form("EMC_STU_G%c%d has no entry!", i + 65, j));
+ }
+
+ if (!arrSTUJ[i][j]) {
+ AliWarning(Form("EMC_STU_J%c%d alias not found!", i + 65, j));
+ } else {
+ if (debug) AliInfo( Form("EMC_STU_J%c%d has %d entries", i + 65, j, arrSTUJ[i][j]->GetEntries()));
+ if (arrSTUJ[i][j]->GetEntries() > 0) {
+ dcsVal = (AliDCSValue*)arrSTUJ[i][j]->At(arrSTUJ[i][j]->GetEntries() - 1);
+ if (dcsVal) {
+ stuConfig->SetJ(i, j, dcsVal->GetInt());
+ if (debug) AliInfo(Form("saving value: %u", dcsVal->GetInt()));
+ }
+ } else
+ AliWarning(Form("EMC_STU_J%c%d has no entry!", i + 65, j));
+ }
+ }
+ }
+
+ arrSTUD = (TObjArray*)dcsAliasMap->GetValue("EMC_STU_GETRAW");
+ arrSTUR = (TObjArray*)dcsAliasMap->GetValue("EMC_STU_REGION");
+ arrSTUF = (TObjArray*)dcsAliasMap->GetValue("EMC_STU_FWVERS");
+
+ if (!arrSTUD) {
+ AliWarning("EMC_STU_GETRAW alias not found!");
+ } else {
+ if (debug) AliInfo(Form("EMC_STU_GETRAW has %d entries", arrSTUD->GetEntries()));
+ if (arrSTUD->GetEntries() > 0) {
+ dcsVal = (AliDCSValue*)arrSTUD->At(arrSTUD->GetEntries() - 1);
+ if (dcsVal) {
+ stuConfig->SetRawData(dcsVal->GetInt());
+ if (debug) AliInfo(Form("saving value: %u", dcsVal->GetInt()));
+ }
+ } else
+ AliWarning("EMC_STU_GETRAW has no entry!");
+ }
+
+ if (!arrSTUR) {
+ AliWarning("EMC_STU_REGION");
+ } else {
+ if (debug) AliInfo( Form("EMC_STU_REGION has %d entries", arrSTUR->GetEntries()));
+ if (arrSTUR->GetEntries() > 0) {
+ dcsVal = (AliDCSValue*)arrSTUR->At(arrSTUR->GetEntries() - 1);
+ if (dcsVal) {
+ stuConfig->SetRegion(dcsVal->GetInt());
+ if (debug) AliInfo(Form("saving value: %u", dcsVal->GetInt()));
+ }
+ } else
+ AliWarning("EMC_STU_REGION has no entry!");
+ }
+
+ if (!arrSTUF) {
+ AliWarning("EMC_STU_FWVERS");
+ } else {
+ if (debug) AliInfo(Form("EMC_STU_FWVERS has %d entries", arrSTUF->GetEntries()));
+ if (arrSTUF->GetEntries() > 0) {
+ dcsVal = (AliDCSValue*)arrSTUF->At(arrSTUF->GetEntries() - 1);
+ if (dcsVal) {
+ stuConfig->SetFw(dcsVal->GetInt());
+ if (debug) AliInfo(Form("saving value: %u", dcsVal->GetInt()));
+ }
+ } else
+ AliWarning("EMC_STU_FWVERS has no entry!");
+ }
+
+ trigConfig->SetSTUObj(stuConfig);
+
+ AliInfo(Form("STU info retrieved."));
+
+
+ // save the objects
AliCDBMetaData metaData;
- Int_t result = Store("Calib", "Data", &calibData, &metaData);
+ metaData.SetBeamPeriod(0);
+ metaData.SetResponsible(kMetaResponsible);
+ metaData.SetComment(kMetaComment);
+
+ UInt_t result=0;
+ Bool_t storeOK = Store("Calib", "Trigger", trigConfig, &metaData, 0, kFALSE);
+ if ( !storeOK ) result=1;
+ AliInfo(Form("TRU info stored. result %d\n", result));
+
+ return result;
+}
+
+//______________________________________________________________________________________________
+UInt_t AliEMCALPreprocessor::ExtractPedestals(Int_t sourceFXS)
+{
+ // Read pedestal file from file exchange server
+ // Only store if new pedestal info is available
+ //
+ UInt_t result=0;
+
+ AliCaloCalibPedestal *calibPed = new AliCaloCalibPedestal(AliCaloCalibPedestal::kEmCal);
+ calibPed->Init();
+
+ TList* list = GetFileSources(sourceFXS,"pedestals");
+ if (list && list->GetEntries()>0) {
+
+ // loop through all files from LDCs
- f.Close();
+ int changes = 0;
+ UInt_t index = 0;
+ while (list->At(index)!=NULL) {
+ TObjString* fileNameEntry = (TObjString*) list->At(index);
+ if (fileNameEntry!=NULL) {
+ TString fileName = GetFile(sourceFXS, "pedestals",
+ fileNameEntry->GetString().Data());
+ TFile *f = TFile::Open(fileName);
+ if (!f) {
+ Log ("Error opening pedestal file.");
+ result = kReturnCodeNoObject;
+ break;
+ }
+ AliCaloCalibPedestal *calPed;
+ f->GetObject("emcCalibPedestal",calPed);
+ if ( !calPed ) {
+ Log ("No pedestal calibration object in file.");
+ result = kReturnCodeNoObject;
+ break;
+ }
+ if ( calPed->GetNEvents()>0 && calPed->GetNChanFills()>0 ) {
+ // add info for the modules available in the present file
+ Bool_t status = calibPed->AddInfo(calPed);
+ if (status) { changes++; }
+ }
+
+ delete calPed;
+ f->Close();
+ }
+ index++;
+ } // while(list)
+
+ //
+ // Store updated pedestal entry to OCDB
+ //
+ if (changes>0) {
+ AliCDBMetaData metaData;
+ metaData.SetBeamPeriod(0);
+ metaData.SetResponsible(kMetaResponsible);
+ metaData.SetComment(kMetaComment);
+
+ Bool_t storeOK = StoreReferenceData("Calib", "Pedestals", calibPed, &metaData);
+ if ( !storeOK ) result++;
+ }
+ }
+ else {
+ Log ("Error: no entries in input file list!");
+ result = kReturnCodeNoEntries;
+ }
+
return result;
+}
+
+//______________________________________________________________________________________________
+UInt_t AliEMCALPreprocessor::ExtractSignal(Int_t sourceFXS)
+{ // Read signal file from file exchange server
+ // Only store if new signal info is available
+ //
+ UInt_t result=0;
+ AliCaloCalibSignal *calibSig = new AliCaloCalibSignal(AliCaloCalibSignal::kEmCal);
+
+ TList* list = GetFileSources(sourceFXS,"signal");
+ if (list && list->GetEntries()>0) {
+
+ // loop through all files from LDCs
+
+ int changes = 0;
+ UInt_t index = 0;
+ while (list->At(index)!=NULL) {
+ TObjString* fileNameEntry = (TObjString*) list->At(index);
+ if (fileNameEntry!=NULL) {
+ TString fileName = GetFile(sourceFXS, "signal",
+ fileNameEntry->GetString().Data());
+ TFile *f = TFile::Open(fileName);
+ if (!f) {
+ Log ("Error opening signal file.");
+ result = kReturnCodeNoObject;
+ break;
+ }
+ AliCaloCalibSignal *calSig;
+ f->GetObject("emcCalibSignal",calSig);
+ if ( !calSig ) {
+ Log ("No signal calibration object in file.");
+ result = kReturnCodeNoObject;
+ break;
+ }
+ if ( calSig->GetNEvents()>0 ) {
+ // add info for the modules available in the present file
+ Bool_t status = calibSig->AddInfo(calSig);
+ if (status) { changes++; }
+ }
+
+ delete calSig;
+ f->Close();
+ }
+ index++;
+ } // while(list)
+
+ //
+ // Store updated signal entry to OCDB
+ //
+ if (changes>0) {
+ AliCDBMetaData metaData;
+ metaData.SetBeamPeriod(0);
+ metaData.SetResponsible(kMetaResponsible);
+ metaData.SetComment(kMetaComment);
+
+ Bool_t storeOK = Store("Calib", "LED", calibSig, &metaData, 0, kFALSE);
+ if ( !storeOK ) result++;
+ }
+ }
+ else {
+ Log ("Error: no entries in input file list!");
+ result = kReturnCodeNoEntries;
+ }
+ return result;
}
+
+