X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;ds=sidebyside;f=PMD%2FAliPMDPreprocessor.cxx;h=c8fa2ab84ea092b581546eba9937952bcd3e0ef1;hb=7ed629b1db0203ea00322199231b2cecb889ef6b;hp=6eac322a9a386c0a85da622016311721d2a656d6;hpb=7fb72bea484d6071ae14b21555dbf8ab55d36e40;p=u%2Fmrichter%2FAliRoot.git diff --git a/PMD/AliPMDPreprocessor.cxx b/PMD/AliPMDPreprocessor.cxx index 6eac322a9a3..c8fa2ab84ea 100644 --- a/PMD/AliPMDPreprocessor.cxx +++ b/PMD/AliPMDPreprocessor.cxx @@ -1,125 +1,546 @@ -// --- ROOT system +/*************************************************************************** + * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * + * * + * Author: The ALICE Off-line Project. * + * Contributors are mentioned in the code where appropriate. * + * * + * Permission to use, copy, modify and distribute this software and its * + * documentation strictly for non-commercial purposes is hereby granted * + * without fee, provided that the above copyright notice appears in all * + * copies and that both the copyright notice and this permission notice * + * appear in the supporting documentation. The authors make no claims * + * about the suitability of this software for any purpose. It is * + * provided "as is" without express or implied warranty. * + **************************************************************************/ + +/********************************************* + + * PMD Preproccessor Source Code * + + + 0 --> Run Successesfully + 1 --> No pmd Alias is available + +**********************************************/ + + #include #include +#include +#include -#include "AliPMDPreprocessor.h" -#include "AliPMDCalibData.h" #include "AliLog.h" #include "AliShuttleInterface.h" #include "AliCDBMetaData.h" -#include -#include -#include -#include - +#include "AliPMDCalibData.h" +#include "AliPMDHotData.h" +#include "AliPMDMeanSm.h" +#include "AliPMDPedestal.h" +#include "AliPMDPreprocessor.h" ClassImp(AliPMDPreprocessor) -//______________________________________________________________________________________________ AliPMDPreprocessor::AliPMDPreprocessor(AliShuttleInterface* shuttle) : AliPreprocessor("PMD", shuttle) { - // constructor + AddRunType("PHYSICS"); + AddRunType("PEDESTAL"); } -//______________________________________________________________________________________________ +//________________________________________________________ AliPMDPreprocessor::~AliPMDPreprocessor() { - // destructor + } -//______________________________________________________________________________________________ -void AliPMDPreprocessor::Initialize(Int_t run, UInt_t startTime, - UInt_t endTime) +//________________________________________________________ +void AliPMDPreprocessor::Initialize(Int_t run, UInt_t startTime, UInt_t endTime) { - // Creates AliPMDDataDAQ object - AliPreprocessor::Initialize(run, startTime, endTime); + + AliInfo(Form("\nRun : %d \nStart Time: %s \nEnd Time : %s", run, + TTimeStamp(startTime).AsString(), + TTimeStamp(endTime).AsString())); + + fRun = run; + fStartTime = startTime; + fEndTime = endTime; + +} - AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run, - TTimeStamp(startTime).AsString(), - TTimeStamp(endTime).AsString())); +//________________________________________________________ +Bool_t AliPMDPreprocessor::ProcessDAQ() +{ + TString sRunType = GetRunType(); + Log(Form("RunType %s",sRunType.Data())); + if (sRunType !="PHYSICS" || sRunType != "PEDESTAL") + { + return kFALSE; + } + return kTRUE; +} - fRun = run; - fStartTime = startTime; - fEndTime = endTime; +//________________________________________________________ +Bool_t AliPMDPreprocessor::StorePmdPED() +{ + + AliPMDPedestal *pedestal = new AliPMDPedestal(); + TList* gfsPmdPed = GetFileSources(kDAQ, "PMD_PED.root"); + + if(!gfsPmdPed) + { + Log(Form("PMDPED: No Shuttle List for PMD PED ")); + return kFALSE; + } + else + { + AliInfo("PMDPED: Here's the list of sources for PMD_PED.root"); + gfsPmdPed->Print(); + + TIter iter(gfsPmdPed); + TObjString* srcPed; + TString nameOfFile; + UInt_t resultPmdPed = 0; + + while((srcPed = dynamic_cast (iter.Next()))) + { + nameOfFile = GetFile(kDAQ, "PMD_PED.root", srcPed->GetName()); + if(nameOfFile.Length() == 0) + { + Log(Form("PMDPED: Error retrieving file from source %s failed!", srcPed->GetName())); + delete gfsPmdPed; + return kFALSE; + } + + Log(Form("PMDPED: File with id PMD_PED.root got from %s", srcPed->GetName())); + + Int_t det, sm, row, col; + Float_t mean, rms; + + TFile *opnFile= new TFile(nameOfFile.Data()); + if(!opnFile || !opnFile->IsOpen()) + { + Log(Form("PMDPED: Error opening file with Id PMD_PED.root from source %s!", srcPed->GetName())); + return kFALSE; + } + + TTree *tree = dynamic_cast (opnFile->Get("ped")); + if (!tree) + { + Log("PMDPED: Could not find object \"ped\" in PED file!"); + return kFALSE; + } + + tree->SetBranchAddress("det", &det); + tree->SetBranchAddress("sm", &sm); + tree->SetBranchAddress("row", &row); + tree->SetBranchAddress("col", &col); + tree->SetBranchAddress("mean", &mean); + tree->SetBranchAddress("rms", &rms); + + Int_t nEntries = (Int_t) tree->GetEntries(); + for(Int_t i = 0; i < nEntries; i++) + { + tree->GetEntry(i); + pedestal->SetPedMeanRms(det,sm,row,col,mean,rms); + } + opnFile->Close(); + delete opnFile; + } + + AliCDBMetaData mdPED; + mdPED.SetBeamPeriod(0); + mdPED.SetResponsible("Satyajit Jena"); + mdPED.SetComment("PMDPED: PMD preprocessor"); + + resultPmdPed = Store("Calib","Ped", pedestal, &mdPED,0,kTRUE); + delete pedestal; + if(resultPmdPed==0) + { + Log("PMDPED: Error storing"); + return kFALSE; + } + else + { + return kTRUE; + } + } + } -//______________________________________________________________________________________________ -UInt_t AliPMDPreprocessor::Process(TMap* pdaqAliasMap) +//________________________________________________________ +Bool_t AliPMDPreprocessor::StorePmdGAIN() +{ +TList* gfsPmdGain = GetFileSources(kDAQ, "PMDGAINS.root"); + +if(!gfsPmdGain) + { + Log(Form("PMDGAIN: No sources found for PMDGAINS.root!")); + return kFALSE; + } + else + { + AliInfo("PMDGAIN: Here's the list of sources for PMDGAINS.root"); + gfsPmdGain->Print(); + + TIter iter(gfsPmdGain); + TObjString* source; + TString nameOfFile; + UInt_t result = 0; + + while((source = dynamic_cast (iter.Next()))) + { + nameOfFile = GetFile(kDAQ, "PMDGAINS.root", source->GetName()); + if(nameOfFile.Length() == 0) + { + Log(Form("PMDGAIN: Error retrieving file from source %s failed!", source->GetName())); + delete gfsPmdGain; + return kFALSE; + } + + Log(Form("PMDGAIN: File with id PMDGAINS.root got from %s", source->GetName())); + + Int_t det, sm, row, col; + Float_t gain; + + TFile *opnFile = new TFile(nameOfFile.Data()); + if (!opnFile || !opnFile->IsOpen()) + { + Log(Form("PMDGAIN: Error opening file with Id PMDGAINS.root from source %s!", source->GetName())); + return kFALSE; + } + + TTree *tree = dynamic_cast (opnFile->Get("ic")); + if (!tree) + { + Log("PMDGAIN: Could not find object \"ic\" in DAQ file!"); + return kFALSE; + } + + tree->SetBranchAddress("det", &det); + tree->SetBranchAddress("sm", &sm); + tree->SetBranchAddress("row", &row); + tree->SetBranchAddress("col", &col); + tree->SetBranchAddress("gain", &gain); + + Int_t nEntries = (Int_t) tree->GetEntries(); + AliPMDCalibData *calibda = new AliPMDCalibData(); + + for(Int_t i = 0; i < nEntries; i++) + { + tree->GetEntry(i); + calibda->SetGainFact(det,sm,row,col,gain); + } + opnFile->Close(); + delete opnFile; + + AliCDBMetaData mdGAIN; + mdGAIN.SetBeamPeriod(0); + mdGAIN.SetResponsible("Satyajit Jena"); + mdGAIN.SetComment("PMDGAIN: PMD GAIN Data"); + result = Store("Calib","Gain", calibda, &mdGAIN); + delete calibda; + } + + if (result==0) + { + Log("PMDGAIN: Error storing"); + return kFALSE; + } + else + { + return kTRUE; + } + + } +} + + +//___________________________________________________ +Bool_t AliPMDPreprocessor::StorePmdHOT() { + + AliPMDHotData *hotda = new AliPMDHotData(); + TList* fsPmdHot = GetFileSources(kDAQ, "PMD_HOT.root"); + + if(!fsPmdHot) + { + Log(Form("PMDHOT: No sources found for PMD_HOT.root!")); + return kFALSE; + } + else + { + AliInfo("PMDHOT: Here's the list of sources for PMD_HOT.root"); + fsPmdHot->Print(); + + TIter iter(fsPmdHot); + TObjString* source; + UInt_t hotresult = 0; + TString nameOfFile; - if(!pdaqAliasMap) return 1; + while((source = dynamic_cast (iter.Next()))) + { + nameOfFile = GetFile(kDAQ, "PMD_HOT.root", source->GetName()); + if(nameOfFile.Length() == 0) + { + Log(Form("PMDHOT: Error retrieving file from source %s failed!", source->GetName())); + delete fsPmdHot; + return kFALSE; + } + + Log(Form("PMDHOT: File with id PMD_HOT.root got from %s", source->GetName())); + + Int_t det, sm, row, col; + Float_t flag; + + TFile *opnFile = new TFile(nameOfFile.Data()); + if(!opnFile || !opnFile->IsOpen()) + { + Log(Form("PMDHOT: Error opening file with Id PMD_HOT.root from source %s!", source->GetName())); + return kFALSE; + } + + TTree *tree = dynamic_cast (opnFile->Get("hot")); + if (!tree) + { + Log("PMDHOT: Could not find object \"hot\" in DAQ file!"); + return kFALSE; + } + + tree->SetBranchAddress("det", &det); + tree->SetBranchAddress("sm", &sm); + tree->SetBranchAddress("row", &row); + tree->SetBranchAddress("col", &col); + tree->SetBranchAddress("flag", &flag); + + Int_t nEntries = (Int_t) tree->GetEntries(); + for(Int_t j = 0; j < nEntries; j++) + { + tree->GetEntry(j); + hotda->SetHotChannel(det,sm,row,col,flag); + } + opnFile->Close(); + delete opnFile; + } + + AliCDBMetaData metaData; + metaData.SetBeamPeriod(0); + metaData.SetResponsible("Satyajit Jena"); + metaData.SetComment("PMDHOT: PMD preprocessor"); + hotresult = Store("Calib","Hot", hotda, &metaData); + delete hotda; + if(hotresult==0) + { + Log("PMDHOT: Error storing"); + return kFALSE; + } + else + { + return kTRUE; + } + + } +} + +//________________________________________________________ +Bool_t AliPMDPreprocessor::StorePmdMEAN() +{ + AliPMDMeanSm *smmeanda = new AliPMDMeanSm(); + TList* gfsPmdMean = GetFileSources(kDAQ, "PMD_MEAN_SM.root"); + + if(!gfsPmdMean) + { + Log(Form("PMDMEAN: No sources found for PMD_MEAN_SM.root!")); + return kFALSE; + } + else + { + AliInfo("PMDMEAN: Here's the list of sources for PMD_MEAN_SM.root"); + gfsPmdMean->Print(); + + TIter iter(gfsPmdMean); + TObjString* sourc; + UInt_t storeMeanData = 0; + TString filenam; + + while((sourc=dynamic_cast (iter.Next()))) + { + filenam = GetFile(kDAQ, "PMD_MEAN_SM.root", sourc->GetName()); + if(filenam.Length() == 0) + { + Log(Form("PMDMEAN: Error retrieving file from source %s failed!", sourc->GetName())); + delete gfsPmdMean; + return kFALSE; + } + + Log(Form("PMDMEAN: File with id PMD_MEAN_SM.root got from %s", sourc->GetName())); + + Int_t det = 0, sm = 0; + Float_t smmean = 0.; + + TFile *opnFile = new TFile(filenam.Data()); + if(!opnFile || !opnFile->IsOpen()) + { + Log(Form("PMDMEAN: Error opening file with Id PMD_MEAN_SM.root from source %s!", sourc->GetName())); + return kFALSE; + } + + TTree *tree = dynamic_cast (opnFile->Get("mean")); + if (!tree) + { + Log("PMDMEAN: Could not find object \"hot\" in DAQ file!"); + return kFALSE; + } + + tree->SetBranchAddress("det", &det); + tree->SetBranchAddress("sm", &sm); + tree->SetBranchAddress("smmean", &smmean); + + Int_t nEntries = (Int_t) tree->GetEntries(); + for(Int_t j = 0; j < nEntries; j++) + { + tree->GetEntry(j); + smmeanda->SetMeanSm(det,sm,smmean); + } + opnFile->Close(); + delete opnFile; + } + + AliCDBMetaData mdMEAN; + mdMEAN.SetBeamPeriod(0); + mdMEAN.SetResponsible("Satyajit Jena"); + mdMEAN.SetComment("PMDMEAN: PMD preprocessor"); + + storeMeanData = Store("Calib","SMMEAN", smmeanda, &mdMEAN); + delete smmeanda; + + if(storeMeanData==0) + { + Log("PMDMEAN: Error storing"); + return kFALSE; + } + else + { + return kTRUE; + } + } + +} + +//_____________________________________________________________ +Bool_t AliPMDPreprocessor::StorePmdDCS(TMap *sDaqAM) +{ + + AliCDBMetaData mdDCS; + mdDCS.SetResponsible("Satyajit Jena"); + mdDCS.SetComment("DCS data for PMD"); + + Bool_t resStore = kFALSE; + resStore = StoreReferenceData("DCS","Data",sDaqAM,&mdDCS); + + if(resStore == 0) + { + Log("PMDDP: Error storing"); + return kFALSE; + } + else + { + return kTRUE; + } + +} + +//_____________________________________________________________________ + +UInt_t AliPMDPreprocessor::Process(TMap* pmdDaqAliasMap) +{ + + if(!pmdDaqAliasMap) + { + return 1; + } + + TString runType = GetRunType(); + + if(runType == "PEDESTAL") + { + + Log(Form("------------------ PMD Pedestal --------------")); + Bool_t pmdPed = StorePmdPED(); + if(!pmdPed) + { + Log(Form("ERROR: Couldn't write PMD pedestal file")); + return 1; + } + else + { + Log(Form("Storing of PMD Pedestal File is Successful")); + return 0; + } + } + + else if (runType == "PHYSICS") + { + Log(Form("------------------- PMD GAIN----------------")); + Bool_t pmdGAIN = StorePmdGAIN(); + if (!pmdGAIN) + { + Log(Form("ERROR: Couldn't write PMD GAIN file")); + + } + else + { + Log(Form("Storing of PMD GAIN File is Successful")); + + } - AliPMDCalibData *calibda = new AliPMDCalibData(); - - TList* filesources = GetFileSources(kDAQ, "PMDGAINS"); - - if(!filesources) { - Log(Form("No sources found for PMDGAINS!")); - return 1; - } - - AliInfo("Here's the list of sources for PMDGAINS"); - filesources->Print(); - - TIter iter(filesources); - TObjString* source; - int i=0; - UInt_t result = 0; - TString filename; - while((source=dynamic_cast (iter.Next()))){ - filename = GetFile(kDAQ, "PMDGAINS", source->GetName()); - if(filename.Length() == 0) { - Log(Form("Error retrieving file from source %s failed!", source->GetName())); - delete filesources; - return 1; - } - - Log(Form("File with id PMDGAINS got from %s", source->GetName())); - Int_t DET,SM,ROW,COL; - Float_t GAIN; - TFile *f= new TFile(filename.Data()); - if(!f || !f->IsOpen()) - { - Log(Form("Error opening file with Id PMDGAINS from source %s!", source->GetName())); - return 1; - } - TTree *tree = dynamic_cast (f->Get("ic")); - if (!tree) - { - Log("Could not find object \"ic\" in DAQ file!"); - return 1; - } - - tree->SetBranchAddress("DET", &DET); - tree->SetBranchAddress("SM", &SM); - tree->SetBranchAddress("ROW", &ROW); - tree->SetBranchAddress("COL", &COL); - tree->SetBranchAddress("GAIN", &GAIN); - Int_t nEntries = (Int_t) tree->GetEntries(); - for(Int_t i = 0; i < nEntries; i++) - { - tree->GetEntry(i); -// if(DET>1 || SM>23 || ROW>95 || COL>95) { -// printf("Error! gain[%d,%d,%d,%d] = %f\n",DET,SM,ROW,COL,GAIN); -// continue; - // } - calibda->SetGainFact(DET,SM,ROW,COL,GAIN); - } - f->Close(); - delete f; - } - - //Now we have to store the final CDB file - AliCDBMetaData metaData; - metaData.SetBeamPeriod(0); - metaData.SetComment("test PMD preprocessor"); - - result = Store("Calib","Data", calibda, &metaData); - - delete calibda; - return result; + Log(Form("------------------- PMD HOT ----------------")); + Bool_t pmdHOT = StorePmdHOT(); + if (!pmdHOT) + { + Log(Form("ERROR: Couldn't write PMD HOT file")); + + } + else + { + Log(Form("Storing of PMD HOT File is Successful")); + + } + + Log(Form("------------------- SM MEAN ----------------")); + Bool_t pmdMEAN = StorePmdMEAN(); + if (!pmdMEAN) + { + Log(Form("ERROR: Couldn't write PMD SM MEAN file")); + + } + else + { + Log(Form("Storing of PMD SM MEAN File is Successful")); + + } + + Log(Form("------------------- DCS DP -----------------")); + Bool_t pmdDCS = StorePmdDCS(pmdDaqAliasMap); + if (!pmdDCS) + { + Log(Form("ERROR: Couldn't write PMD DCS DP file")); + + } + else + { + Log(Form("Storing of PMD DCS dp is Successul")); + + } + + if (pmdGAIN || pmdHOT || pmdMEAN || pmdDCS) + return 0; + else + return 1; + + } + return 2; } +