X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FPHOSLEDda.cxx;h=21fcdcbcd2f4274a5e7017e8dee363011bcc64e8;hb=c64e9808fab3161aba11f934db98417e62be2fd3;hp=bdabd7f50a95e2719b9c01c88e63eb905f4d7b5c;hpb=11999cae0d11fedd040775667ae7f00090bc15e3;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/PHOSLEDda.cxx b/PHOS/PHOSLEDda.cxx index bdabd7f50a9..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 */ @@ -27,8 +27,10 @@ extern "C" { #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 @@ -43,6 +45,9 @@ int main(int argc, char **argv) { "RIO", "TStreamerInfo()"); + AliLog::SetGlobalDebugLevel(0) ; + AliLog::SetGlobalLogLevel(AliLog::kFatal); + int status; if (argc!=2) { @@ -50,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]); @@ -62,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) { @@ -97,18 +144,35 @@ int main(int argc, char **argv) { /* init some counters */ int nevents_physics=0; int nevents_total=0; - - Int_t fMod = 2; // module 2! + AliRawReader *rawReader = NULL; - - AliPHOSRcuDA1 da1(fMod,-1,0); // DA1 for module2, no input/output file + AliPHOSRcuDA1* dAs[5]; + for(Int_t iMod=0; iMod<5; iMod++) { + dAs[iMod] = 0; + } + 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; + Int_t nFired = -1; + Int_t sigStart, sigLength; + Int_t caloFlag; + - Int_t gain = -1; - Int_t X = -1; - Int_t Z = -1; + TH1I fFiredCells("fFiredCells","Number of fired cells per event",100,0,1000); /* main loop (infinite) */ for(;;) { @@ -141,35 +205,83 @@ 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); + nFired = 0; - while(dc.NextDigit()) { - - X = dc.GetRow() - 1; - Z = dc.GetColumn() - 1; + rawReader = new AliRawReaderDate((void*)event); + AliCaloRawStreamV3 stream(rawReader,"PHOS",mapping); + AliPHOSRawFitterv0 fitter; + fitter.SubtractPedestals(kTRUE); // assume that data is non-ZS - 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(caloFlag==1 && fitter.GetEnergy()>40) + nFired++; + } - e[X][Z][gain] = dc.GetEnergy(); - t[X][Z][gain] = dc.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,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); + fFiredCells.Fill(nFired); delete rawReader; nevents_physics++; @@ -187,52 +299,59 @@ 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 localfile[128]; - sprintf(localfile,"PHOS_Module%d_LED.root",fMod); - TFile* f = new TFile(localfile,"recreate"); const TH2F* h2=0; const TH1F* h1=0; + char localfile[128]; Int_t nGood=0; // >10 entries in peak Int_t nMax=-111; // max. number of entries in peak Int_t iXmax=-1; Int_t iZmax=-1; + Int_t iModMax=-1; + + sprintf(localfile,"PHOS_LED.root"); + TFile* f = new TFile(localfile,"recreate"); - for(Int_t iX=0; iX<64; iX++) { - for(Int_t iZ=0; iZ<56; iZ++) { - - h1 = da1.GetHgLgRatioHistogram(iX,iZ); // High Gain/Low Gain ratio - if(h1) { - if(h1->GetMaximum()>10.) nGood++; - if(h1->GetMaximum()>nMax) {nMax = (Int_t)h1->GetMaximum(); iXmax=iX; iZmax=iZ;} - h1->Write(); - } + 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++) { + + h1 = dAs[iMod]->GetHgLgRatioHistogram(iX,iZ); // High Gain/Low Gain ratio + if(h1) { + if(h1->GetMaximum()>10.) nGood++; + if(h1->GetMaximum()>nMax) { + nMax = (Int_t)h1->GetMaximum(); iXmax=iX; iZmax=iZ; iModMax=iMod; + } + h1->Write(); + } - for(Int_t iGain=0; iGain<2; iGain++) { - h2 = da1.GetTimeEnergyHistogram(iX,iZ,iGain); // Time vs Energy - if(h2) h2->Write(); + for(Int_t iGain=0; iGain<2; iGain++) { + 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 */ - - for(Int_t iMod=0; iMod<5; iMod++) { - sprintf(localfile,"PHOS_Module%d_LED.root",iMod); - daqDA_FES_storeFile(localfile,"LED"); - } + 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)",fMod,iXmax,iZmax); + printf("(module,iX,iZ)=(%d,%d,%d)",iModMax,iXmax,iZmax); printf("\n"); return status;