/*
contact: Boris.Polishchuk@cern.ch
link: http://aliceinfo.cern.ch/static/phpBB3/viewtopic.php?f=4&t=17
-reference run: /castor/cern.ch/alice/phos/2007/10/02/13/07000008232001.10.root
-run type: STANDALONE
+reference run: /alice/data/2009/LHC09c_PHOS/000098979/raw
+run type: LED
DA type: MON
number of events needed: 1000
-input files: RCU0.data RCU1.data RCU2.data RCU3.data
-Output files: PHOS_Module2_BCM.root
-Trigger types used: CALIBRATION_EVENT or PHYSICS
+number of events needed: 1000
+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
+Output files: PHOS_BCM.root
+Trigger types used: CALIBRATION_EVENT
*/
#include <TSystem.h>
#include <TROOT.h>
#include <TPluginManager.h>
+#include <TH1.h>
+#include <TH2.h>
#include "AliRawReader.h"
#include "AliRawReaderDate.h"
#include "AliPHOSDA2.h"
-#include "AliPHOSRawFitterv1.h"
+#include "AliPHOSRawFitterv3.h"
#include "AliCaloAltroMapping.h"
#include "AliCaloRawStreamV3.h"
+#include "AliLog.h"
/* Main routine
"RIO",
"TStreamerInfo()");
+ AliLog::SetGlobalDebugLevel(0) ;
+ AliLog::SetGlobalLogLevel(AliLog::kFatal);
+
int status;
if (argc!=2) {
return -1;
}
-
+ short offset=-1;
+ short threshold=-1;
+
/* Retrieve mapping files from DAQ DB */
- const char* mapFiles[4] = {"RCU0.data","RCU1.data","RCU2.data","RCU3.data"};
-
- for(Int_t iFile=0; iFile<4; iFile++) {
+ 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) {
+ 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] );
int nevents_total=0;
AliRawReader *rawReader = NULL;
-
- AliPHOSDA2* da2 = new AliPHOSDA2(2); // DA2 ("Checking for bad channels") for module2
+ AliPHOSDA2* dAs[5];
- Float_t q[64][56][2];
+ for(Int_t iMod=0; iMod<5; iMod++) {
+ dAs[iMod] = 0;
+ }
- Int_t gain = -1;
+ Float_t q[64][56][2];
+
Int_t cellX = -1;
Int_t cellZ = -1;
Int_t nBunches = 0;
Int_t nFired = -1;
Int_t sigStart, sigLength;
+ Int_t caloFlag;
/* main loop (infinite) */
for(;;) {
rawReader = new AliRawReaderDate((void*)event);
AliCaloRawStreamV3 stream(rawReader,"PHOS",mapping);
- AliPHOSRawFitterv0 fitter();
+ AliPHOSRawFitterv3 fitter;
fitter.SubtractPedestals(kTRUE); // assume that data is non-ZS
while (stream.NextDDL()) {
while (stream.NextChannel()) {
+ /* Retrieve ZS parameters from data*/
+ 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 = fRawStream.GetStartTimeBin();
- sigLength = fRawStream.GetBunchLength();
- fitter.SetSamples(fRawStream->GetSignals(),sigStart,sigLength);
+ sigStart = stream.GetStartTimeBin();
+ sigLength = stream.GetBunchLength();
+ fitter.SetChannelGeo(stream.GetModule(),cellX,cellZ,caloFlag);
+ fitter.Eval(stream.GetSignals(),sigStart,sigLength);
+ q[cellX][cellZ][caloFlag] = fitter.GetSignalQuality();
} // End of NextBunch()
- fitter.SetNBunches(nBunches);
- fitter.SetChannelGeo(module,cellX,cellZ,caloFlag);
- fitter.Eval();
-
- if (nBunches>1 || caloFlag!=0 || caloFlag!=1 || fitter.GetSignalQuality()>1) continue;
-
- q[cellX][cellZ][caloFlag] = fitter.GetSignalQuality();
-
- if(gain && dc.GetEnergy()>40)
+ if(caloFlag==1 && fitter.GetEnergy()>40)
nFired++;
}
}
-
- da2->FillQualityHistograms(q);
- da2->FillFiredCellsHistogram(nFired);
- //da1.UpdateHistoFile();
+ if(stream.GetModule()<0 || stream.GetModule()>4) continue;
+
+ if(dAs[stream.GetModule()]) {
+ dAs[stream.GetModule()]->FillQualityHistograms(q);
+ dAs[stream.GetModule()]->FillFiredCellsHistogram(nFired);
+ }
+ else {
+ dAs[stream.GetModule()] = new AliPHOSDA2(stream.GetModule(),0);
+ dAs[stream.GetModule()]->FillQualityHistograms(q);
+ dAs[stream.GetModule()]->FillFiredCellsHistogram(nFired);
+ }
+
delete rawReader;
nevents_physics++;
}
}
}
- 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 */
- delete da2;
+
+ char lnam[128];
+ char ltitl[128];
+
+ char hnam[128];
+ char htitl[128];
+
+ const TH1F* hist1=0;
+ TH2F* maps[2];
+
+ TFile* f = new TFile("PHOS_BCM.root","recreate");
+
+ for(Int_t iMod=0; iMod<5; iMod++) {
+ if(!dAs[iMod]) continue;
- /* Store output files to the File Exchange Server */
- daqDA_FES_storeFile("PHOS_Module2_BCM.root","BAD_CHANNELS");
+ printf("DA2 for module %d detected.\n",iMod);
+ sprintf(lnam,"gmaplow%d",iMod);
+ sprintf(ltitl,"Quality map for Low gain in Module %d",iMod);
+
+ sprintf(hnam,"gmaphigh%d",iMod);
+ sprintf(htitl,"Quality map for High gain in Module %d",iMod);
+
+ maps[0] = new TH2F(lnam, ltitl, 64,0.,64.,56,0.,56.);
+ maps[1] = new TH2F(hnam, htitl, 64,0.,64.,56,0.,56.);
+
+ for(Int_t iX=0; iX<64; iX++) {
+ for(Int_t iZ=0; iZ<56; iZ++) {
+
+ for(Int_t iGain=0; iGain<2; iGain++) {
+ hist1 = dAs[iMod]->GetQualityHistogram(iX,iZ,iGain);
+ if(hist1) {
+ hist1->Write();
+ Double_t mean = hist1->GetMean();
+ maps[iGain]->SetBinContent(iX+1,iZ+1,mean);
+ }
+ }
+ }
+ }
+
+ maps[0]->Write(); delete maps[0];
+ maps[1]->Write(); delete maps[1];
+
+ }
+
+ f->Close();
+
+ if(offset>0 && threshold>0)
+ printf("ZS parameters: offset %d, threshold %d.\n",offset,threshold);
+
+ /* Store output files to the File Exchange Server */
+ daqDA_FES_storeFile("PHOS_BCM.root","BAD_CHANNELS");
+
return status;
}