X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FPHOSLEDda.cxx;h=21fcdcbcd2f4274a5e7017e8dee363011bcc64e8;hb=a7141e6cbf981d328596cd2ebeccd26988c335ed;hp=1a6bd30851c1913c6614c1a3181ff429e6867ab9;hpb=d62e649ae2c6fc1e9bab27ba012195d8154bef65;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/PHOSLEDda.cxx b/PHOS/PHOSLEDda.cxx index 1a6bd30851c..21fcdcbcd2f 100644 --- a/PHOS/PHOSLEDda.cxx +++ b/PHOS/PHOSLEDda.cxx @@ -1,12 +1,12 @@ /* 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: STANDALONE +reference run: /alice/data/2009/LHC09b_PHOS/000075883/raw/09000075883017.20.root +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_LED.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_LED.root. Trigger types used: CALIBRATION_EVENT */ @@ -30,6 +30,7 @@ extern "C" { #include "AliPHOSRawFitterv0.h" #include "AliCaloAltroMapping.h" #include "AliCaloRawStreamV3.h" +#include "AliLog.h" /* Main routine @@ -44,6 +45,9 @@ int main(int argc, char **argv) { "RIO", "TStreamerInfo()"); + AliLog::SetGlobalDebugLevel(0) ; + AliLog::SetGlobalLogLevel(AliLog::kFatal); + int status; if (argc!=2) { @@ -51,10 +55,42 @@ 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); + + short 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[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]); @@ -63,18 +99,28 @@ int main(int argc, char **argv) { } /* 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) { @@ -108,7 +154,16 @@ int main(int argc, char **argv) { Float_t e[64][56][2]; Float_t t[64][56][2]; - + + 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; @@ -156,14 +211,36 @@ int main(int argc, char **argv) { AliCaloRawStreamV3 stream(rawReader,"PHOS",mapping); AliPHOSRawFitterv0 fitter; fitter.SubtractPedestals(kTRUE); // assume that data is non-ZS + + 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()) { @@ -171,15 +248,11 @@ int main(int argc, char **argv) { if (nBunches > 1) continue; sigStart = stream.GetStartTimeBin(); sigLength = stream.GetBunchLength(); - fitter.SetSamples(stream.GetSignals(),sigStart,sigLength); + fitter.SetChannelGeo(stream.GetModule(),cellX,cellZ,caloFlag); + fitter.Eval(stream.GetSignals(),sigStart,sigLength); } // End of NextBunch() - fitter.SetNBunches(nBunches); - fitter.SetChannelGeo(stream.GetModule(),cellX,cellZ,caloFlag); - fitter.Eval(); - - if(nBunches>1) continue; -// if (nBunches>1 || caloFlag!=0 || caloFlag!=1 || fitter.GetSignalQuality()>1) continue; + if(nBunches != 1) continue; e[cellX][cellZ][caloFlag] = fitter.GetEnergy(); t[cellX][cellZ][caloFlag] = fitter.GetTime(); @@ -188,10 +261,14 @@ int main(int argc, char **argv) { nFired++; } + if(stream.GetModule()<0 || stream.GetModule()>4) continue; + if(dAs[stream.GetModule()]) dAs[stream.GetModule()]->FillHistograms(e,t); - else + 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++) { @@ -222,7 +299,7 @@ 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 */ @@ -236,12 +313,13 @@ int main(int argc, char **argv) { Int_t iZmax=-1; Int_t iModMax=-1; + sprintf(localfile,"PHOS_LED.root"); + TFile* f = new TFile(localfile,"recreate"); + for(Int_t iMod=0; iMod<5; iMod++) { if(!dAs[iMod]) continue; printf("DA1 for module %d detected.\n",iMod); - sprintf(localfile,"PHOS_Module%d_LED.root",iMod); - TFile* f = new TFile(localfile,"recreate"); for(Int_t iX=0; iX<64; iX++) { for(Int_t iZ=0; iZ<56; iZ++) { @@ -259,17 +337,18 @@ int main(int argc, char **argv) { h2 = dAs[iMod]->GetTimeEnergyHistogram(iX,iZ,iGain); // Time vs Energy if(h2) h2->Write(); } - + } } - fFiredCells.Write(); - f->Close(); - - /* Store output files to the File Exchange Server */ - daqDA_FES_storeFile(localfile,"LED"); + fFiredCells.Write(); } + f->Close(); + + /* Store output files to the File Exchange Server */ + daqDA_FES_storeFile(localfile,"LED"); + printf("%d physics events of %d total processed.\n",nevents_physics,nevents_total); printf("%d histograms has >10 entries in maximum, max. is %d entries ",nGood,nMax); printf("(module,iX,iZ)=(%d,%d,%d)",iModMax,iXmax,iZmax);