X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PHOS%2FPHOSGAINda.cxx;h=803281f2c8659764d6d54a8f9605c242e4bce664;hb=51c3974f6891429a15c5210aad2796c0b60931da;hp=06598d05a10e50742c9d207421735e90d2ed70e9;hpb=2905f08d979f40df2b4759a689f9bdda8f9d8cc9;p=u%2Fmrichter%2FAliRoot.git diff --git a/PHOS/PHOSGAINda.cxx b/PHOS/PHOSGAINda.cxx index 06598d05a10..803281f2c86 100644 --- a/PHOS/PHOSGAINda.cxx +++ b/PHOS/PHOSGAINda.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 +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 */ @@ -30,7 +30,7 @@ extern "C" { #include "AliPHOSRawFitterv0.h" #include "AliCaloAltroMapping.h" #include "AliCaloRawStreamV3.h" - +#include "AliLog.h" /* Main routine Arguments: @@ -44,6 +44,9 @@ int main(int argc, char **argv) { "RIO", "TStreamerInfo()"); + AliLog::SetGlobalDebugLevel(0) ; + AliLog::SetGlobalLogLevel(AliLog::kFatal); + int status; if (argc!=2) { @@ -63,9 +66,30 @@ int main(int argc, char **argv) { } /* 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]); @@ -74,15 +98,26 @@ 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 */ @@ -178,7 +213,21 @@ int main(int argc, char **argv) { 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 @@ -198,7 +247,7 @@ int main(int argc, char **argv) { 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(); @@ -209,7 +258,7 @@ int main(int argc, char **argv) { if(dAs[stream.GetModule()]) dAs[stream.GetModule()]->FillHistograms(e,t); else { - dAs[stream.GetModule()] = new AliPHOSRcuDA1(stream.GetModule(),-1); + dAs[stream.GetModule()] = new AliPHOSRcuDA1(stream.GetModule(),-1,0); dAs[stream.GetModule()]->FillHistograms(e,t); } @@ -243,22 +292,44 @@ 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 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; }