X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FPHOSGAINda.cxx;h=803281f2c8659764d6d54a8f9605c242e4bce664;hb=70f2ec6e03381aa6cebd589cd2de7b8822279cf4;hp=bc0c994ad5c7f8294768a063a61c9d3108cacf86;hpb=62adef6ca81aaf5b4d779d23dbc808a0a8f3bb08;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/PHOSGAINda.cxx b/PHOS/PHOSGAINda.cxx index bc0c994ad5c..803281f2c86 100644 --- a/PHOS/PHOSGAINda.cxx +++ b/PHOS/PHOSGAINda.cxx @@ -1,13 +1,13 @@ /* contact: Boris.Polishchuk@cern.ch link: see comments in the $ALICE_ROOT/PHOS/AliPHOSRcuDA1.cxx -reference run: /castor/cern.ch/alice/phos/2007/10/04/18/07000008249001.1000.root -run type: PHYSICS, STANDALONE +reference run: /alice/data/2009/LHC09b_PHOS/000075883/raw/09000075883017.20.root +run type: PHYSICS DA type: MON number of events needed: 1000 -input files: RCU0.data RCU1.data RCU2.data RCU3.data -Output files: PHOS_Module2_Calib.root -trigger types used: L0 +input files: Mod0RCU0.data Mod0RCU1.data Mod0RCU2.data Mod0RCU3.data Mod1RCU0.data Mod1RCU1.data Mod1RCU2.data Mod1RCU3.data Mod2RCU0.data Mod2RCU1.data Mod2RCU2.data Mod2RCU3.data Mod3RCU0.data Mod3RCU1.data Mod3RCU2.data Mod3RCU3.data Mod4RCU0.data Mod4RCU1.data Mod4RCU2.data Mod4RCU3.data zs.txt +Output files: PHOS_Calib_Total.root contains cumulative statistics for a number of runs. +Trigger types used: PHYSICS */ @@ -21,13 +21,16 @@ extern "C" { #include #include +#include +#include #include "AliRawReader.h" #include "AliRawReaderDate.h" #include "AliPHOSRcuDA1.h" -#include "AliPHOSRawDecoder.h" +#include "AliPHOSRawFitterv0.h" #include "AliCaloAltroMapping.h" - +#include "AliCaloRawStreamV3.h" +#include "AliLog.h" /* Main routine Arguments: @@ -35,6 +38,15 @@ extern "C" { */ int main(int argc, char **argv) { + gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo", + "*", + "TStreamerInfo", + "RIO", + "TStreamerInfo()"); + + AliLog::SetGlobalDebugLevel(0) ; + AliLog::SetGlobalLogLevel(AliLog::kFatal); + int status; if (argc!=2) { @@ -42,19 +54,72 @@ int main(int argc, char **argv) { return -1; } + /* Retrieve ZS parameters from DAQ DB */ + const char* zsfile = "zs.txt"; + int failZS = daqDA_DB_getFile(zsfile, zsfile); + + Int_t offset,threshold; + + if(!failZS) { + FILE *f = fopen(zsfile,"r"); + int scan = fscanf(f,"%d %d",&offset,&threshold); + } + + /* Retrieve mapping files from DAQ DB */ + const char* mapFiles[20] = { + "Mod0RCU0.data", + "Mod0RCU1.data", + "Mod0RCU2.data", + "Mod0RCU3.data", + "Mod1RCU0.data", + "Mod1RCU1.data", + "Mod1RCU2.data", + "Mod1RCU3.data", + "Mod2RCU0.data", + "Mod2RCU1.data", + "Mod2RCU2.data", + "Mod2RCU3.data", + "Mod3RCU0.data", + "Mod3RCU1.data", + "Mod3RCU2.data", + "Mod3RCU3.data", + "Mod4RCU0.data", + "Mod4RCU1.data", + "Mod4RCU2.data", + "Mod4RCU3.data" + }; + + for(Int_t iFile=0; iFile<20; iFile++) { + int failed = daqDA_DB_getFile(mapFiles[iFile], mapFiles[iFile]); + if(failed) { + printf("Cannot retrieve file %s from DAQ DB. Exit.\n",mapFiles[iFile]); + return -1; + } + } + /* Open mapping files */ - AliAltroMapping *mapping[4]; + AliAltroMapping *mapping[20]; TString path = "./"; - path += "RCU"; + + path += "Mod"; TString path2; - for(Int_t i = 0; i < 4; i++) { + TString path3; + Int_t iMap = 0; + + for(Int_t iMod = 0; iMod < 5; iMod++) { path2 = path; - path2 += i; - path2 += ".data"; - mapping[i] = new AliCaloAltroMapping(path2.Data()); + path2 += iMod; + path2 += "RCU"; + + for(Int_t iRCU=0; iRCU<4; iRCU++) { + path3 = path2; + path3 += iRCU; + path3 += ".data"; + mapping[iMap] = new AliCaloAltroMapping(path3.Data()); + iMap++; + } } - /* define data source : this is argument 1 */ status=monitorSetDataSource( argv[1] ); if (status!=0) { @@ -80,15 +145,29 @@ int main(int argc, char **argv) { int nevents_total=0; AliRawReader *rawReader = NULL; + AliPHOSRcuDA1* dAs[5]; - AliPHOSRcuDA1 da1(2,-1); // DA1 (Calibration DA) for module2 + for(Int_t iMod=0; iMod<5; iMod++) { + dAs[iMod] = 0; + } Float_t e[64][56][2]; Float_t t[64][56][2]; - Int_t gain = -1; - Int_t X = -1; - Int_t Z = -1; + for(Int_t iX=0; iX<64; iX++) { + for(Int_t iZ=0; iZ<56; iZ++) { + for(Int_t iGain=0; iGain<2; iGain++) { + e[iX][iZ][iGain] = 0.; + t[iX][iZ][iGain] = 0.; + } + } + } + + Int_t cellX = -1; + Int_t cellZ = -1; + Int_t nBunches = 0; + Int_t sigStart, sigLength; + Int_t caloFlag; /* main loop (infinite) */ for(;;) { @@ -121,39 +200,84 @@ int main(int argc, char **argv) { if (eventT==PHYSICS_EVENT) { - for(Int_t iX=0; iX<64; iX++) { - for(Int_t iZ=0; iZ<56; iZ++) { - for(Int_t iGain=0; iGain<2; iGain++) { - e[iX][iZ][iGain] = 0.; - t[iX][iZ][iGain] = 0.; - } - } - } - rawReader = new AliRawReaderDate((void*)event); -// AliPHOSRawDecoderv1 dc(rawReader,mapping); - AliPHOSRawDecoder dc(rawReader,mapping); - dc.SubtractPedestals(kTRUE); + AliCaloRawStreamV3 stream(rawReader,"PHOS",mapping); + AliPHOSRawFitterv0 fitter; + fitter.SubtractPedestals(kTRUE); // assume that data is non-ZS - while(dc.NextDigit()) { - - X = dc.GetRow() - 1; - Z = dc.GetColumn() - 1; - - if(dc.IsLowGain()) gain = 0; - else - gain = 1; + if(!failZS) { + fitter.SubtractPedestals(kFALSE); + fitter.SetAmpOffset(offset); + fitter.SetAmpThreshold(threshold); + } + + while (stream.NextDDL()) { + while (stream.NextChannel()) { + + /* Retrieve ZS parameters from data*/ + if(failZS) { + short value = stream.GetAltroCFG1(); + bool ZeroSuppressionEnabled = (value >> 15) & 0x1; + bool AutomaticBaselineSubtraction = (value >> 14) & 0x1; + if(ZeroSuppressionEnabled) { + offset = (value >> 10) & 0xf; + threshold = value & 0x3ff; + fitter.SubtractPedestals(kFALSE); + fitter.SetAmpOffset(offset); + fitter.SetAmpThreshold(threshold); + } + } + + cellX = stream.GetCellX(); + cellZ = stream.GetCellZ(); + caloFlag = stream.GetCaloFlag(); // 0=LG, 1=HG, 2=TRU + + if(caloFlag!=0 && caloFlag!=1) continue; //TRU data! + + // In case of oscillating signals with ZS, + //a channel can have several bunches. + + nBunches = 0; + while (stream.NextBunch()) { + nBunches++; + if (nBunches > 1) continue; + sigStart = stream.GetStartTimeBin(); + sigLength = stream.GetBunchLength(); + fitter.SetChannelGeo(stream.GetModule(),cellX,cellZ,caloFlag); + fitter.Eval(stream.GetSignals(),sigStart,sigLength); + } // End of NextBunch() + + if (nBunches != 1) continue; + + e[cellX][cellZ][caloFlag] = fitter.GetEnergy(); + t[cellX][cellZ][caloFlag] = fitter.GetTime(); + } + + if(stream.GetModule()<0 || stream.GetModule()>4) continue; - e[X][Z][gain] = dc.GetEnergy(); - t[X][Z][gain] = dc.GetTime(); + if(dAs[stream.GetModule()]) + dAs[stream.GetModule()]->FillHistograms(e,t); + else { + dAs[stream.GetModule()] = new AliPHOSRcuDA1(stream.GetModule(),-1,0); + dAs[stream.GetModule()]->FillHistograms(e,t); + } + for(Int_t iX=0; iX<64; iX++) { + for(Int_t iZ=0; iZ<56; iZ++) { + for(Int_t iGain=0; iGain<2; iGain++) { + e[iX][iZ][iGain] = 0.; + t[iX][iZ][iGain] = 0.; + } + } + } + } - da1.FillHistograms(e,t); - //da1.UpdateHistoFile(); +// da1.FillHistograms(e,t); +// //da1.UpdateHistoFile(); - delete rawReader; - nevents_physics++; + delete rawReader; + nevents_physics++; } nevents_total++; @@ -168,7 +292,44 @@ int main(int argc, char **argv) { } } - for(Int_t i = 0; i < 4; i++) delete mapping[i]; + for(Int_t i = 0; i < 20; i++) delete mapping[i]; + + /* Be sure that all histograms are saved */ + + char h2name[80]; + char totfile[80]; + + //Write the Total file (accumulated statistics for number of runs) + sprintf(totfile,"PHOS_Calib_Total.root"); + TFile * ftot = new TFile(totfile,"recreate"); + +// if (!ftot->IsZombie()){ +// printf("Updating file %s.\n",ftot->GetName()); + +// for(Int_t iMod=0; iMod<5; iMod++) { +// if(!dAs[iMod]) continue; + +// printf("DA1 for module %d detected.\n",iMod); + +// for(Int_t iX=0; iX<64; iX++) { +// for(Int_t iZ=0; iZ<56; iZ++) { + +// for(Int_t iGain=0; iGain<2; iGain++) { +// sprintf(h2name,"%d_%d_%d_%d",iMod,iX,iZ,iGain); +// TH2F* h2tot = (TH2F*)ftot->Get(h2name); +// const TH2F* h2run = dAs[iMod]->GetTimeEnergyHistogram(iX,iZ,iGain); // Time vs Energy +// if(!h2tot && h2run) h2run->Write(); +// if(h2tot && h2run) { h2tot->Add(h2run); h2tot->Write(h2tot->GetName(),TObject::kWriteDelete); } +// } +// } +// } +// } +// } + + ftot->Close(); + + /* Store output files to the File Exchange Server */ + daqDA_FES_storeFile(totfile,"AMPLITUDES"); return status; }