#include "monitor.h"
#include <stdio.h>
#include <stdlib.h>
-#include <pthread.h>
#include <vector>
//
//Root includes
#include "TString.h"
#include "TObjString.h"
#include "TDatime.h"
+#include "TStopwatch.h"
#include "TMap.h"
#include "TGraph.h"
#include "TMath.h"
#include "AliRawReader.h"
#include "AliRawReaderDate.h"
#include "AliTPCmapper.h"
-#include "AliTPCRawStream.h"
#include "AliTPCROC.h"
#include "AliTPCCalROC.h"
#include "AliTPCCalPad.h"
//functios, implementation below
-void SendToAmoreDB(AliTPCCalibCE &calibCE, unsigned long32 runNb);
+void SendToAmoreDB(AliTPCCalibCE *calibCE, unsigned long32 runNb);
//for threaded processing
-void *processEventBuffer(void *arg);
-
-//common event processing variables for threaded processing
-std::vector<eventHeaderStruct*> eventBuffer;
-volatile int bStop = false;
-struct timespec duree_nanosleep;
-Int_t forceNevents=-1;
-Bool_t forceBufferEnds=kFALSE;
/* Main routine
*/
int main(int argc, char **argv) {
/* log start of process */
- printf("TPC CE DA started - %s\n",__FILE__);
+ printf("TPCCEda: DA started - %s\n",__FILE__);
if (argc<2) {
- printf("Wrong number of arguments\n");
+ printf("TPCCEda: Wrong number of arguments\n");
return -1;
}
- AliLog::SetClassDebugLevel("AliTPCRawStream",-5);
AliLog::SetClassDebugLevel("AliRawReaderDate",-5);
AliLog::SetClassDebugLevel("AliTPCAltroMapping",-5);
AliLog::SetModuleDebugLevel("RAW",-5);
int i,status;
status=monitorDeclareMp( __FILE__ );
if (status!=0) {
- printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
+ printf("TPCCEda: monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
return -1;
}
monitorSetNowait();
char localfile[255];
unsigned long32 runNb=0; //run number
- //
- // thread wait time
- //
- duree_nanosleep.tv_sec=0;
- duree_nanosleep.tv_nsec=1000000; //1 ms
-
+
//
// DA configuration from configuration file
//
sprintf(localfile,"./%s",CONFIG_FILE);
status = daqDA_DB_getFile(CONFIG_FILE,localfile);
if (status) {
- printf("Failed to get configuration file (%s) from DAQdetDB, status=%d\n", CONFIG_FILE, status);
+ printf("TPCCEda: Failed to get configuration file (%s) from DAQdetDB, status=%d\n", CONFIG_FILE, status);
return -1;
}
AliTPCConfigDA config(CONFIG_FILE);
// check configuration options
TString laserTriggerName("C0LSR-ABCE-NOPF-CENT");
- size_t bufferSize=30;
+ TString monitoringType("YES");
Int_t forceTriggerId=-1;
- Int_t forceNeventsStandalone=-1;
- Bool_t forceBufferEndsGlobal=kTRUE;
-// Bool_t forceBufferEndsGlobalDummy=kFALSE;
+ Int_t saveOption=2; // how to store the object. See AliTPCCalibCE::DumpToFile
+ Bool_t skipAmore=kFALSE;
+
if ( config.GetConfigurationMap()->GetValue("LaserTriggerName") ) {
laserTriggerName=config.GetConfigurationMap()->GetValue("LaserTriggerName")->GetName();
- printf("Laser trigger class name set to: %s.\n",laserTriggerName.Data());
+ printf("TPCCEda: Laser trigger class name set to: %s.\n",laserTriggerName.Data());
}
-
- if ( config.GetConfigurationMap()->GetValue("BufferSize") ) {
- bufferSize=(size_t)config.GetValue("BufferSize");
- printf("Setting event buffer size to: %d.\n",bufferSize);
+
+ if ( config.GetConfigurationMap()->GetValue("MonitoringType") ) {
+ monitoringType=config.GetConfigurationMap()->GetValue("MonitoringType")->GetName();
+ printf("TPCCEda: Monitoring type set to: %s.\n",monitoringType.Data());
}
-
- if ( config.GetConfigurationMap()->GetValue("ForceTriggerId") ) {
- forceTriggerId=TMath::Nint(config.GetValue("ForceTriggerId"));
- printf("Only processing triggers with Id: %d.\n",forceTriggerId);
+
+ if ( config.GetConfigurationMap()->GetValue("ForceLaserTriggerId") ) {
+ forceTriggerId=TMath::Nint(config.GetValue("ForceLaserTriggerId"));
+ printf("TPCCEda: Only processing triggers with Id: %d.\n",forceTriggerId);
}
- if ( config.GetConfigurationMap()->GetValue("ForceBufferEndsGlobal") ) {
- forceBufferEndsGlobal=config.GetValue("ForceBufferEndsGlobal")!=0.;
- printf("Process all buffered events in global partition: %s.\n",forceBufferEndsGlobal?"yes":"no");
+ if ( config.GetConfigurationMap()->GetValue("SaveOption") ) {
+ saveOption=TMath::Nint(config.GetValue("SaveOption"));
+ printf("TPCCEda: Saving option set to: %d.\n",saveOption);
}
-
- if ( config.GetConfigurationMap()->GetValue("ForceNMaxEvents") ) {
- forceNevents=TMath::Nint(config.GetValue("ForceNMaxEvents"));
- printf("Forcing maximum number of %d events.\n",forceNeventsStandalone);
+
+ if ( config.GetConfigurationMap()->GetValue("SkipAmore") ) {
+ skipAmore=((TObjString*)config.GetConfigurationMap()->GetValue("SkipAmore"))->GetString().Atoi();
+ printf("TPCCEda: Skip Amore set in config\n");
}
//subsribe to laser triggers only in physics partition
//if the trigger class is not available the return value is -1
//in this case we are most probably running as a standalone
// laser run and should request all events
- unsigned char *classIdptr=0;
- int retClassId=daqDA_getClassIdFromName(laserTriggerName.Data(),classIdptr);
+ unsigned char classIdptr=0;
+ int retClassId=daqDA_getClassIdFromName(laserTriggerName.Data(),&classIdptr);
if (retClassId==0){
//interleaved laser in physics runs
+ //select proper trigger class id
char c[5];
- sprintf(c,"%u",*classIdptr);
- char *table[5] = {"PHY","Y","*",c,NULL};
+ snprintf(c,sizeof(c),"%u",(unsigned int)classIdptr);
+ char *table[5] = {"PHY",const_cast<char*>(monitoringType.Data()),"*",c,NULL};
monitorDeclareTableExtended(table);
- printf("Using trigger class Id: %s\n",c);
-// forceBufferEndsGlobal=forceBufferEndsGlobalDummy;
+ printf("TPCCEda: Using monitoring table: (PHY, %s, *, %s)\n",monitoringType.Data(),c);
} else if (retClassId==-1){
- //defaul case, accept all physics events
- char *table[3] = {"PHY","Y",NULL};
+ //global partition without laser triggered events
+ //the DA should exit properly without processing
+ printf("TPCCEda: Laser trigger class '%s' was not found among trigger class names. Will stop processing.\n",laserTriggerName.Data());
+ return 0;
+ } else if (retClassId==-2){
+ //standalone case, accept all physics events
+ char *table[5] = {"PHY","Y","*","*",NULL};
monitorDeclareTableExtended(table);
- printf("Using all trigger class Ids\n");
-// forceNevents=forceNeventsStandalone;
+ printf("TPCCEda: Using all trigger class Ids\n");
} else {
- printf("Unknown return value of 'daqDA_getClassIdFromName': %d\n",retClassId);
+ printf("TPCCEda: Unknown return value of 'daqDA_getClassIdFromName': %d\n",retClassId);
return -2;
}
sprintf(c,"%d",forceTriggerId);
char *table[5] = {"PHY","Y","*",c,NULL};
monitorDeclareTableExtended(table);
-// forceBufferEndsGlobal=forceBufferEndsGlobalDummy;
}
sprintf(localfile,"./%s",MAPPING_FILE);
status = daqDA_DB_getFile(MAPPING_FILE,localfile);
if (status) {
- printf("Failed to get mapping file (%s) from DAQdetDB, status=%d\n", MAPPING_FILE, status);
+ printf("TPCCEda: Failed to get mapping file (%s) from DAQdetDB, status=%d\n", MAPPING_FILE, status);
return -1;
}
}
if (mapping == 0) {
- printf("Failed to get mapping object from %s. ...\n", MAPPING_FILE);
+ printf("TPCCEda: Failed to get mapping object from %s. ...\n", MAPPING_FILE);
return -1;
} else {
- printf("Got mapping object from %s\n", MAPPING_FILE);
+ printf("TPCCEda: Got mapping object from %s\n", MAPPING_FILE);
}
//create calibration object
- AliTPCCalibCE calibCE(config.GetConfigurationMap()); // central electrode calibration
- calibCE.SetAltroMapping(mapping->GetAltroMapping()); // Use altro mapping we got from daqDetDb
-
- //
- // start thread
- //
-// sleep(5);
- pthread_t threadId=0;
- int threadStatus=0;
- threadStatus = pthread_create( &threadId, NULL, processEventBuffer, (void*)(&calibCE));
- eventBuffer.resize(bufferSize);
- struct timespec duree_out;
+ AliTPCCalibCE *calibCE=new AliTPCCalibCE(config.GetConfigurationMap()); // central electrode calibration
+ calibCE->SetAltroMapping(mapping->GetAltroMapping()); // Use altro mapping we got from daqDetDb
+
+ //amore update interval
+ Double_t updateInterval=300; //seconds
+ Double_t valConf=config.GetValue("AmoreUpdateInterval");
+ if ( valConf>0 ) updateInterval=valConf;
+ //timer
+ TStopwatch stopWatch;
+
//===========================//
// loop over RAW data files //
//==========================//
int nevents=0;
+ int neventsOld=0;
size_t counter=0;
for ( i=1; i<argc; i++) {
/* define data source : this is argument i */
- printf("Processing file %s\n", argv[i]);
+ printf("TPCCEda: Processing file %s\n", argv[i]);
status=monitorSetDataSource( argv[i] );
if (status!=0) {
- printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
+ printf("TPCCEda: monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
return -1;
}
-
+ Bool_t hasNewData=kFALSE;
/* read until EOF */
while (true) {
struct eventHeaderStruct *event;
/* check shutdown condition */
if (daqDA_checkShutdown()) {break;}
-
- //check for predefined number of events
- if (forceNevents>0 && calibCE.GetNeventsProcessed()>=forceNevents) {
- printf("Requested number of events reached (%d).\n",forceNevents);
+
+ /* get next event (blocking call until timeout) */
+ status=monitorGetEventDynamic((void **)&event);
+ if (status==MON_ERR_EOF) {
+ printf ("TPCCEda: End of File %d detected\n",i);
+ break; /* end of monitoring file has been reached */
+ }
+
+ if (status!=0) {
+ printf("TPCCEda: monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
break;
}
- //buffer events, only read them if a buffer position is free
- if (eventBuffer[counter]==0) {
-
- /* get next event (blocking call until timeout) */
- status=monitorGetEventDynamic((void **)&event);
- if (status==MON_ERR_EOF) {
- printf ("End of File %d detected\n",i);
- break; /* end of monitoring file has been reached */
+ /* retry if got no event */
+ if (event==NULL){
+ //use time in between bursts to
+ // send the data to AMOREdb
+ if (stopWatch.RealTime()>updateInterval && hasNewData){
+ calibCE->Analyse();
+ if (!skipAmore) SendToAmoreDB(calibCE,runNb);
+ stopWatch.Start();
+ } else {
+ stopWatch.Continue();
}
-
- if (status!=0) {
- printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
- break;
+ //debug output
+ if (nevents>neventsOld){
+ printf ("TPCCEda: %d events processed, %d used\n",nevents,calibCE->GetNeventsProcessed());
+ neventsOld=nevents;
}
-
- /* retry if got no event */
- if (event==NULL)
- continue;
+ hasNewData=kFALSE;
+ continue;
+ }
/* skip start/end of run events */
- if ( (event->eventType != physicsEvent) && (event->eventType != calibrationEvent) ){
- free(event);
- continue;
- }
-
-
- // get the run number
- runNb = event->eventRunNb;
-
-// printf(" trigger (%05d-%03d) = %8.8x %8.8x - %02u\n",nevents, calibCE.GetNeventsProcessed(),
-// event->eventTriggerPattern[1], event->eventTriggerPattern[0],event->eventType);
-
-// printf("filling buffer %d\n",counter);
- eventBuffer[counter]=event;
-
- ++nevents;
- ++counter;
- if (counter >= eventBuffer.size()) {counter=0;}
- }else{
-// printf("buffer already used: %d\n",counter);
- nanosleep(&duree_nanosleep,&duree_out);
+ if ( (event->eventType != physicsEvent) && (event->eventType != calibrationEvent) ){
+ free(event);
+ continue;
}
+
+
+ // get the run number
+ runNb = event->eventRunNb;
+
+ // CE calibration
+ calibCE->ProcessEvent(event);
+ hasNewData=kTRUE;
+
+ /* free resources */
+ free(event);
+ ++nevents;
}
}
-
- //
- // wait for thread to end
- //
- if (!forceBufferEndsGlobal) bStop = true;
- else forceBufferEnds=forceBufferEndsGlobal;
- pthread_join( threadId, NULL);
-// printf("Event Processing Thread ended with: %d\n",threadStatus);
-
- //
- // free unprocessed events
- //
- for (size_t i=0;i<eventBuffer.size();++i){
- if (eventBuffer[i]) {
- free(eventBuffer[i]);
- eventBuffer[i]=0;
-// printf("freeing buffer %d\n",i);
- }
- }
-
//
// Analyse CE data and write them to rootfile
//
- calibCE.Analyse();
- printf ("%d events processed, %d used\n",nevents,calibCE.GetNeventsProcessed());
-
- TFile * fileTPC = new TFile (RESULT_FILE,"recreate");
- calibCE.Write("tpcCalibCE");
- delete fileTPC;
- printf("Wrote %s\n",RESULT_FILE);
+ printf ("TPCCEda: %d events processed, %d used\n",nevents,calibCE->GetNeventsProcessed());
+
+ //save data to file
+ calibCE->DumpToFile(RESULT_FILE,Form("name=tpcCalibCE,type=%d",saveOption));
+ printf("TPCCEda: Wrote %s\n",RESULT_FILE);
/* store the result file on FES */
-
status=daqDA_FES_storeFile(RESULT_FILE,FILE_ID);
if (status) {
status = -2;
}
+
+ if (!skipAmore){
+ printf("TPCCEda: Amore part\n");
+ calibCE->Analyse();
+ SendToAmoreDB(calibCE,runNb);
+ }
- SendToAmoreDB(calibCE,runNb);
-
+ delete calibCE;
return status;
}
-void *processEventBuffer(void *arg)
-{
- //
- // event procssing thread functio
- //
-
- //cast argument
- AliTPCCalibCE *ce=(AliTPCCalibCE*)arg;
- AliTPCCalibCE &calibCE=*ce;
-
- size_t counter=0;
- unsigned long32 runNb=0;
- Bool_t published=kTRUE;
- struct timespec duree_out;
-
- struct eventHeaderStruct *event;
-
- //wait for the first buffer to be filled
- while (!eventBuffer[0]) nanosleep(&duree_nanosleep,&duree_out);
- //loop over buffer
- while (!bStop){
-// printf("testing buffer: %d\n",counter);
- if (eventBuffer[counter]) {
- event=eventBuffer[counter];
- runNb = event->eventRunNb;
-// printf("processing buffer: %d\n",counter);
- eventBuffer[counter]=0;
- calibCE.ProcessEvent(event);
- free(event);
- published=kFALSE;
- } else {
- //in case of empty buffer publish the results it this was not done
- if (!published) {
- SendToAmoreDB(calibCE,runNb);
- published=kTRUE;
- }
- nanosleep(&duree_nanosleep,&duree_out);
- }
- ++counter;
- if (counter >= eventBuffer.size()) {
- counter=0;
- if (forceBufferEnds) break;
- }
- }
-}
-void SendToAmoreDB(AliTPCCalibCE &calibCE, unsigned long32 runNb)
+void SendToAmoreDB(AliTPCCalibCE *calibCE, unsigned long32 runNb)
{
//AMORE
- printf ("AMORE part\n");
+// printf ("AMORE part\n");
const char *amoreDANameorig=gSystem->Getenv("AMORE_DA_NAME");
//cheet a little -- temporary solution (hopefully)
//
gSystem->Setenv("AMORE_DA_NAME",Form("TPC-%s",FILE_ID));
//
// end cheet
- TGraph *grA=calibCE.MakeGraphTimeCE(-1,0,2);
- TGraph *grC=calibCE.MakeGraphTimeCE(-2,0,2);
+ TGraph *grA=calibCE->MakeGraphTimeCE(-1,0,2);
+ TGraph *grC=calibCE->MakeGraphTimeCE(-2,0,2);
TDatime time;
TObjString info(Form("Run: %u; Date: %s",runNb,time.AsSQLString()));
amore::da::AmoreDA amoreDA(amore::da::AmoreDA::kSender);
Int_t statusDA=0;
- statusDA+=amoreDA.Send("CET0",calibCE.GetCalPadT0());
- statusDA+=amoreDA.Send("CEQ",calibCE.GetCalPadQ());
- statusDA+=amoreDA.Send("CERMS",calibCE.GetCalPadRMS());
+ statusDA+=amoreDA.Send("CET0",calibCE->GetCalPadT0());
+ statusDA+=amoreDA.Send("CEQ",calibCE->GetCalPadQ());
+ statusDA+=amoreDA.Send("CERMS",calibCE->GetCalPadRMS());
statusDA+=amoreDA.Send("DriftA",grA);
statusDA+=amoreDA.Send("DriftC",grC);
statusDA+=amoreDA.Send("Info",&info);
if ( statusDA!=0 )
- printf("Waring: Failed to write one of the calib objects to the AMORE database\n");
+ printf("TPCCEda: Waring: Failed to write one of the calib objects to the AMORE database\n");
// reset env var
if (amoreDANameorig) gSystem->Setenv("AMORE_DA_NAME",amoreDANameorig);
if (grA) delete grA;