X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=HMPID%2FHMPIDda.cxx;h=e5baad333c2e909185dd8af3166eb81d7ab5e9c7;hb=2bf46337a64073c10c80e20b626da602bb973a86;hp=cbc7b03be1b1ca9e18e526680063eaed17318a96;hpb=cb9b108e24a1d7363fd413dffc9baaf67a657d5e;p=u%2Fmrichter%2FAliRoot.git diff --git a/HMPID/HMPIDda.cxx b/HMPID/HMPIDda.cxx index cbc7b03be1b..e5baad333c2 100644 --- a/HMPID/HMPIDda.cxx +++ b/HMPID/HMPIDda.cxx @@ -2,14 +2,15 @@ HMPID DA for online calibration -Contact: Levente.Molnar@ba.infn.it, Giacomo.Volpe@ba.infn.it -Link: http://richpc1.ba.infn.it/~levente/Files4Public/ValidateHmpidDA/ -Run Type: PEDESTAL -- but we select on the PHYSICS_EVENTS in th HMPIDda.cxx -DA Type: LDC -Number of events needed: 1000 events -Input Files: Raw pedestal file, EXTERNAL config file: HmpidSigmaCut.txt on both HMPID LDCs -Output Files: 14 txt files including pedestal values -Trigger types used: PEDESTAL RUN (selecting on PHYSICS_EVENT) +Contact: Levente.Molnar@cern.ch, Giacomo.Volpe@ba.infn.it +Link: https://twiki.cern.ch/twiki/bin/view/ALICE/DAInstructions +Run Type: PEDESTAL -- but we select on the PHYSICS_EVENTS in the HMPIDda.cxx +DA Type: LDC +Reference Run: 78734 +Number of events needed: ~ 1000 events +Input Files: Raw pedestal file, EXTERNAL config files: HmpDaqDaConfig.txt, HmpDeadChannelMap.txt on all HMPID LDCs +Output Files: 2 x 14 txt files including pedestal and error values +Trigger types used: CALIBRATION RUN (but selecting on PHYSICS_EVENT) */ @@ -39,63 +40,57 @@ extern "C" { #include "TObject.h" #include "TPluginManager.h" +//AMORE +#include + int main(int argc, char **argv){ - int status; + int status=0; + const Char_t *hmpConfigFile = "HmpDaqDaConfig.txt"; + const Char_t *hmpDeadChannelMapFile = "HmpDeadChannelMap.txt"; + const Int_t ddlOffset = 1536; + TString hmpIn,hmpIn2; + TString feeIn; + /* log start of process */ - printf("HMPID DA program started\n"); + printf("HMP PedestalDa: HMPID DA program started\n"); - gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo", - "*", - "TStreamerInfo", - "RIO", - "TStreamerInfo()"); + gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo","*","TStreamerInfo","RIO","TStreamerInfo()"); /* check that we got some arguments = list of files */ if (argc<2) { - printf("Wrong number of arguments\n"); + printf("HMP PedestalDa: Wrong number of arguments\n"); return -1; } - /* copy locally a file from daq detector config db - status=daqDA_DB_getFile("myconfig","./myconfig.txt"); - if (status) { - printf("Failed to get config file : %d\n",status); - return -1; - } - and possibly use it */ + /* copy locally a file from daq detector config db */ + + hmpIn=Form("./%s",hmpConfigFile); + status=daqDA_DB_getFile(hmpConfigFile,hmpIn.Data()); + if (status) { printf("Failed to get HMPID config file status: %d\n",status);return -1; } + hmpIn2=Form("./%s",hmpDeadChannelMapFile); + status=daqDA_DB_getFile(hmpDeadChannelMapFile,hmpIn2.Data()); + if (status) { printf("Failed to get HMPID dead channel file status: %d\n",status);return -1; } + /* report progress */ daqDA_progressReport(10); - - - + /* define wait event timeout - 1s max */ monitorSetNowait(); monitorSetNoWaitNetworkTimeout(1000); /* set local storage of ped files for Fe2C */ - + /* init the pedestal calculation */ AliHMPIDCalib *pCal=new AliHMPIDCalib(); - /* Set the number of sigma cuts inside the file HmpidSigmaCut.txt on BOTH LDCs! */ + /* Set the number of sigma cuts inside the file HmpidSigmaCut.txt on all LDCs! */ /* If the file is NOT present then the default cut 3 will be used!*/ - - //pCal->SetSigCutFromFile("HmpidSigmaCut.txt"); - pCal->SetSigCutFromShell("HMPID_SIGMA_CUT"); //decision to make later: wether to use file or env variable... - pCal->SetDaOutFromShell("HMPID_DA_OUT"); //decision to make later: wether to use file or env variable... - pCal->SetFeeInFromShell("HMPID_FEE_IN"); //decision to make later: wether to use file or env variable... - - - /* ONLY set this option to kTRUE if you want to create the ADC dsitributions for all 161280 pads!!!!*/ - /* kTRUE is not suggested for production mode b/c of the memory consumption! */ - pCal->SetWriteHistoPads(kFALSE); //use this option for default production useage!!! - //pCal->SetWriteHistoPads(kTRUE); //only for expert debug - //pCal->SetWriteHistoPads(kTRUE,kTRUE,13); //DO NOT USE THIS OPTION! - + pCal->SetSigCutFromFile(hmpIn); + pCal->SetDeadChannelMapFromFile(hmpIn2); /* init event counter */ Int_t firstEvt=0; @@ -108,7 +103,7 @@ int main(int argc, char **argv){ status=monitorSetDataSource( argv[n] ); if (status!=0) { - printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status)); + printf("HMP PedestalDa: monitorSetDataSource() failed : %s\n",monitorDecodeError(status)); return -1; } @@ -128,13 +123,11 @@ int main(int argc, char **argv){ status=monitorGetEventDynamic((void **)&event); if (status==MON_ERR_EOF) /* end of monitoring file has been reached */ { - printf("End of monitoring file has been reached! \n"); + printf("HMP PedestalDa: End of monitoring file has been reached! \n"); break; } - - if (status!=0) { - printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status)); + printf("HMP PedestalDa: monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status)); return -1; } @@ -150,11 +143,6 @@ int main(int argc, char **argv){ if (eventT==PHYSICS_EVENT || eventT==CALIBRATION_EVENT ) { //updated: 18/02/2008 based on http://alice-ecs.web.cern.ch/alice-ecs/runtypes_3.16.html runNum=(unsigned long)event->eventRunNb; //assuming that only one run is processed at a time ldcId=(unsigned long)event->eventLdcId; - if(iEvtNcal==firstEvt && pCal->GetWritePads()==kTRUE) - { - if(pCal->GetLargePads()==kFALSE) pCal->InitFile((Int_t)ldcId); //The number for iEvtNcal should be the same as for the first value - else pCal->InitFile((Int_t)runNum); //The number for iEvtNcal should be the same as for the first value - } iEvtNcal++; AliRawReader *reader = new AliRawReaderDate((void*)event); @@ -179,7 +167,7 @@ int main(int argc, char **argv){ /* exit when last event received, no need to wait for TERM signal */ if (eventT==END_OF_RUN) { - printf("EOR event detected\n"); + printf("HMP PedestalDa: EOR event detected\n"); break; } // events loop @@ -190,41 +178,59 @@ int main(int argc, char **argv){ }//arg /* write report */ - printf("HMPID DA processed RUN #%s on LDC#%d, with %d calibration events\n",getenv("DATE_RUN_NUMBER"),ldcId,iEvtNcal); + printf("HMP PedestalDa: HMPID DA processed RUN #%s on LDC#%d, with %d calibration events\n",getenv("DATE_RUN_NUMBER"),ldcId,iEvtNcal); if (!iEvtNcal) { - printf("No calibration events have been read. Exiting\n"); + printf("HMP PedestalDa: No calibration events have been read. Exiting\n"); return -1; } /* report progress */ daqDA_progressReport(90); - TString sDaOut=pCal->GetDaOutFromShell(); - - for(Int_t nDDL=0; nDDL < AliHMPIDRawStream::kNDDL; nDDL++) { - + /* send pedestal and error files to DAQ FXS */ + for(Int_t nDDL=0; nDDL < AliHMPIDRawStream::kNDDL; nDDL++) { + feeIn=Form("thr%d.dat",ddlOffset+nDDL); /* Calculate pedestal for the given ddl, if there is no ddl go t next */ - if(!pCal->CalcPedestal(nDDL,Form("%sHmpidPedDdl%02i.txt",sDaOut.Data(),nDDL),iEvtNcal)) continue; - /* to create pedestal file as Paolo uncomment the line */ - //if(!pCal->CalcPedestalPaolo(nDDL,Form("%sHmpidPedDdl%02i.txt",sDaOut.Data(),nDDL),iEvtNcal)) continue; - if(!pCal->WriteErrors(nDDL,Form("%sHmpidErrorsDdl%02i.txt",sDaOut.Data(),nDDL),iEvtNcal)) continue; - - /* store the result file on FES */ - - status=daqDA_FES_storeFile(Form("%sHmpidPedDdl%02i.txt",sDaOut.Data(),nDDL),Form("HmpidPedDdl%02i.txt",nDDL)); - if (status) { printf("Failed to export file : %d\n",status); } - status=daqDA_FES_storeFile(Form("%sHmpidErrorsDdl%02i.txt",sDaOut.Data(),nDDL),Form("HmpidErrorsDdl%02i.txt",nDDL)); - if (status) { printf("Failed to export file : %d\n",status); } - + if(pCal->CalcPedestal(nDDL,Form("HmpidPedDdl%02i.txt",nDDL),Form("%s",feeIn.Data()),iEvtNcal)) { + status=daqDA_DB_storeFile(feeIn.Data(),feeIn.Data()); //store a single threshold file for a DDL in DAQ DB + if (status) { printf("HMP PedestalDa: Failed to store file %s in DAQ DB, status: %d\n",feeIn.Data(),status); } + status=daqDA_FES_storeFile(Form("HmpidPedDdl%02i.txt",nDDL),Form("HmpidPedDdl%02i.txt",nDDL)); //store a single pedestal file for a DDL + if (status) { printf("HMP PedestalDa: Failed to export file on FES: %d\n",status); } + }//pedestals and thresholds + if(pCal->WriteErrors(nDDL,Form("HmpidErrorsDdl%02i.txt",nDDL),iEvtNcal)) { + status=daqDA_FES_storeFile(Form("HmpidErrorsDdl%02i.txt",nDDL),Form("HmpidErrorsDdl%02i.txt",nDDL)); + if (status) { printf("HMP PedestalDa: Failed to export file : %d\n",status); } + }//errors }//nDDL - if(pCal->GetWritePads()==kTRUE) pCal->CloseFile(); + /* send files to AMORE DB */ + daqDA_progressReport(95); + Int_t statusAmoreDA=0; + amore::da::AmoreDA amoreDA(amore::da::AmoreDA::kSender); + for(Int_t iCh=0; iCh < AliHMPIDParam::kMaxCh; iCh++) { + statusAmoreDA+=amoreDA.Send(Form("fPedMeanMap%d",iCh), pCal->GetPedMeanMap(iCh)); + statusAmoreDA+=amoreDA.Send(Form("fPedSigMap%d",iCh), pCal->GetPedSigMap(iCh)); + } + for(Int_t iCh=0;iCh<=AliHMPIDParam::kMaxCh;iCh++) + { + for(Int_t iFee=0;iFee<6;iFee++) + { + statusAmoreDA+=amoreDA.Send(Form("f1DPedMean_Ch%d_FEE_%d",iCh,iFee),pCal->GetPedMean(6*iCh+iFee)); + statusAmoreDA+=amoreDA.Send(Form("f1DPedSigma_Ch%d_FEE_%d",iCh,iFee),pCal->GetPedSigma(6*iCh+iFee)); + } + } + + printf("HMP PedestalDa: num. masked pads: %d, num of currently dead pads: %d.\n",pCal->GetNumMaskedPads(),pCal->GetNumDeadPads()); + delete pCal; - if (status) return -1; + if(status) return -1; + if(statusAmoreDA) return -1; /* report progress */ daqDA_progressReport(100); + return status; } +