X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FPHOSLEDda.cxx;h=21fcdcbcd2f4274a5e7017e8dee363011bcc64e8;hb=a7141e6cbf981d328596cd2ebeccd26988c335ed;hp=5e2c3c36b4bb81886d6aa2bd260455e328403365;hpb=8083e8192f1a02daf3c26eeea7d521c48e2ce548;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/PHOSLEDda.cxx b/PHOS/PHOSLEDda.cxx index 5e2c3c36b4b..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"}; - - 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) { 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) { @@ -98,21 +144,33 @@ 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]; - - Int_t gain = -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 nFired = -1; Int_t sigStart, sigLength; + Int_t caloFlag; + TH1I fFiredCells("fFiredCells","Number of fired cells per event",100,0,1000); @@ -147,57 +205,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.; - } - } - } - nFired = 0; - + rawReader = new AliRawReaderDate((void*)event); AliCaloRawStreamV3 stream(rawReader,"PHOS",mapping); - AliPHOSRawFitterv0 fitter(); + 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()) { 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); } // 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; + if(nBunches != 1) continue; e[cellX][cellZ][caloFlag] = fitter.GetEnergy(); t[cellX][cellZ][caloFlag] = fitter.GetTime(); - + if(caloFlag==1 && fitter.GetEnergy()>40) nFired++; } + + 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++; } @@ -214,49 +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(); } - fFiredCells.Write(); f->Close(); /* Store output files to the File Exchange Server */ - daqDA_FES_storeFile("PHOS_Module2_LED.root","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;