/*
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
+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
+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
*/
#include "AliPHOSRawFitterv0.h"
#include "AliCaloAltroMapping.h"
#include "AliCaloRawStreamV3.h"
-
+#include "AliLog.h"
/* Main routine
Arguments:
"RIO",
"TStreamerInfo()");
+ AliLog::SetGlobalDebugLevel(0) ;
+ AliLog::SetGlobalLogLevel(AliLog::kFatal);
+
int status;
if (argc!=2) {
}
/* Retrieve mapping files from DAQ DB */
- const char* mapFiles[4] = {"RCU0.data","RCU1.data","RCU2.data","RCU3.data"};
+ 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<4; iFile++) {
+ 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]);
}
/* 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 */
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.
fitter.Eval(stream.GetSignals(),sigStart,sigLength);
} // End of NextBunch()
- if (nBunches>1) continue;
+ if (nBunches != 1) continue;
e[cellX][cellZ][caloFlag] = fitter.GetEnergy();
t[cellX][cellZ][caloFlag] = fitter.GetTime();
}
-
+
+ if(stream.GetModule()<0 || stream.GetModule()>4) continue;
+
if(dAs[stream.GetModule()])
dAs[stream.GetModule()]->FillHistograms(e,t);
- else
- dAs[stream.GetModule()] = new AliPHOSRcuDA1(stream.GetModule(),-1);
+ 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 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];
- char localfile[128];
+ //Write the Total file (accumulated statistics for number of runs)
+ sprintf(totfile,"PHOS_Calib_Total.root");
+ TFile * ftot = new TFile(totfile,"recreate");
- for(Int_t iMod=0; iMod<5; iMod++) {
- if(!dAs[iMod]) continue;
+// if (!ftot->IsZombie()){
+// printf("Updating file %s.\n",ftot->GetName());
+
+// for(Int_t iMod=0; iMod<5; iMod++) {
+// if(!dAs[iMod]) continue;
- dAs[iMod]->UpdateHistoFile();
- dAs[iMod]->SetWriteToFile(kFALSE);
+// printf("DA1 for module %d detected.\n",iMod);
- /* Store output files to the File Exchange Server */
- sprintf(localfile,"PHOS_Module%d_Calib.root",iMod);
- daqDA_FES_storeFile(localfile,"AMPLITUDES");
- }
+// 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;
}