[u/mrichter/AliRoot.git] / TOF / TOFda.cxx

/*

TOF DA for online calibration

Contact: Chiara.Zampolli@bo.infn.it
Link: www.bo.infn.it/~zampolli
Run Type: PHYSICS
DA Type: MON
Number of events needed: depending on the run, being run-level
Input Files: TOFdaTotal.root, to be updated if existing
Output Files: TOFdaRun.root, TOFdaTotal.root, both to be exported to the DAQ FXS
Trigger types used: PHYSICS_EVENT

*/

#define FILE_TOTAL "TOFdaTotal.root"
#define FILE_RUN "TOFdaRun.root"

// DATE
#include "event.h"
#include "monitor.h"
#include "daqDA.h"

#include <stdio.h>
#include <stdlib.h>

//AliRoot
#include <AliTOFRawStream.h>
#include <AliRawReaderDate.h>
#include <AliRawReader.h>
#include <AliTOFGeometry.h>
#include <AliT0RawReader.h>
#include <AliDAQ.h>
#include <AliTOFHitData.h>
#include <AliTOFHitDataBuffer.h>
#include <AliTOFDaConfigHandler.h>

//ROOT
#include <TFile.h>
#include <TKey.h>
#include <TH2S.h>
#include <TObject.h>
#include <TMath.h>
#include <TSystem.h>
#include "TROOT.h"
#include "TPluginManager.h"
#include "TSAXParser.h"

/* Main routine
      Arguments: 
      1- monitoring data source
*/
int main(int argc, char **argv) {

  /* magic line from Rene */
  gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
                    "*",
                    "TStreamerInfo",
                    "RIO",
                    "TStreamerInfo()");

  AliTOFGeometry * geom = new AliTOFGeometry();

  static const Int_t size = AliTOFGeometry::NPadXSector()*AliTOFGeometry::NSectors();
  static const Int_t nbins = 500;
  static const Int_t binmin = -20;
  const Float_t c = 2.99792458E10; //speed of light [cm/s]
  TH1F::AddDirectory(0);
  TH2S * htofPartial = new TH2S("htof","histo with delays",
				size,-0.5,size*1.-0.5,
				nbins,binmin-0.5,nbins*1.+binmin-0.5);
  
  int status;

  /* log start of process */
  printf("TOF DA started\n");  

  /* check that we got some arguments = list of files */
  if (argc!=2) {
    printf("Wrong number of arguments\n");
    return -1;
  }

  /* retrieve config file */
  int getConfigFile = daqDA_DB_getFile("TOFPhysicsConfig.xml","TOFPhysicsConfig.xml");
  if (getConfigFile != 0){
    printf("Failed to retrieve config file from DB! returning...\n");
    return -1;
  }

  AliTOFDaConfigHandler* tofHandler = new AliTOFDaConfigHandler();
  TSAXParser *parser = new TSAXParser();
  parser->ConnectToHandler("AliTOFDaConfigHandler", tofHandler);  
  if (parser->ParseFile("./TOFPhysicsConfig.xml") != 0) {
    printf("Failed parsing config file! retunring... \n");
    return -1;
  }

  Int_t debugFlag = tofHandler->GetDebugFlag();
  printf("The debug flag is %i\n",debugFlag);
  Int_t t0Flag = tofHandler->GetT0Flag();
  printf("The T0 flag is %i\n\n",t0Flag);
  if (t0Flag) {
	  printf("The T0 time will be subtracted from the measured TOF time. So, in TDC bins: \n");
	  printf("tof = tofRaw - (rawReaderT0->GetData(51,0)+rawReaderT0->GetData(52,0))/2) \n\n");
  }
  else {
	  printf("The T0 time will not be used.\n");
	  printf("tof = tofRaw \n\n");
  }

  /* init some counters */
  int nevents_physics=0;
  int nevents_total=0;

  Int_t iev=0;

  Int_t nPDBEntriesToT = 0;
  Int_t nDBEntriesToT = 0;
  AliTOFHitData *HitData;
  Int_t dummy = -1;
  Int_t Volume[5];
  for (Int_t i=0;i<5;i++) Volume[i]=-1;
  AliTOFRawStream *rawStreamTOF = new AliTOFRawStream();
  AliTOFHitDataBuffer *DataBuffer;
  AliTOFHitDataBuffer *PackedDataBuffer;
  Int_t nDBEntries = 0;
  Int_t nPDBEntries = 0;
  
  struct eventHeaderStruct *event;
  eventTypeType eventT;

  /* define data source : this is argument 1 */  
  status=monitorSetDataSource( argv[1] );
  if (status!=0) {
    printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
    return -1;
  }

  /* declare monitoring program */
  status=monitorDeclareMp( __FILE__ );
  if (status!=0) {
    printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
    return -1;
  }

  /* define wait event timeout - 1s max */
  monitorSetNowait();
  monitorSetNoWaitNetworkTimeout(1000);
  
  /* main loop (infinite) */
  for(;;) {
    
    /* check shutdown condition */
    if (daqDA_checkShutdown()) break;
    
    /* get next event (blocking call until timeout) */
    status=monitorGetEventDynamic((void **)&event);
    if (status==MON_ERR_EOF) {
      printf ("End of File detected\n");
      break; /* end of monitoring file has been reached */
    }
    
    if (status!=0) {
      printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
      break;
    }
    
    /* retry if got no event */
    if (event==NULL) continue;

    iev++; 

   /* use event - here, just write event id to result file */
    nevents_total++;
    eventT=event->eventType;
    switch (event->eventType) {
      
      /* START OF RUN */
    case START_OF_RUN:
      break;
      /* END START OF RUN */
      
      /* END OF RUN */
    case END_OF_RUN:
      break;
      /* END END OF RUN */
      
    case PHYSICS_EVENT:
      nevents_physics++;
      AliRawReader *rawReader = new AliRawReaderDate((void*)event);
      //rawReader->RequireHeader(kFALSE);

      //T0 event
      Int_t meantime = 0;     
      AliT0RawReader *rawReaderT0 = new AliT0RawReader(rawReader,kTRUE);
      if (!rawReaderT0->Next()) {
        printf("T0: no raw data found!\n");
      } 
      else {
	/*
        Int_t allData[105][5];
        for (Int_t i=0; i<105; i++) {
	  allData[i][0]=rawReaderT0->GetData(i,0);
	}
        meantime = allData[49][0];
	*/
        //meantime = rawReaderT0->GetData(49,0); //OLD
        meantime = (Int_t)((rawReaderT0->GetData(51,0)+rawReaderT0->GetData(52,0))/2.); //Alla
	if (debugFlag > 0) {
		printf("\nT0 for the current event:\n");   // debugging purpose
		printf("time zero = %i (TDC bin)\n", meantime);   // debugging purpose
		printf("time zero = %f (ns)\n\n", (Float_t)(meantime)*24.4*1E-3);   // debugging purpose
	}
      }
      
      delete rawReaderT0;
      rawReaderT0 = 0x0;
      rawReader->Reset();
      
      //TOF event
      dummy = -1;
      for (Int_t ii=0; ii<5; ii++) Volume[ii]=-1;
      rawStreamTOF->SetRawReader(rawReader);
      //rawReader->ReadHeader();
      rawStreamTOF->ResetBuffers();
      rawStreamTOF->DecodeDDL(0, AliDAQ::NumberOfDdls("TOF") - 1,0);
      nPDBEntriesToT = 0;
      nDBEntriesToT = 0;
      for (Int_t iDDL = 0; iDDL < AliDAQ::NumberOfDdls("TOF"); iDDL++) {

	/* read decoded data */
	DataBuffer = rawStreamTOF->GetDataBuffer(iDDL);
	PackedDataBuffer = rawStreamTOF->GetPackedDataBuffer(iDDL);
	
	/* get buffer entries */
	nDBEntries = DataBuffer->GetEntries();
	nPDBEntries = PackedDataBuffer->GetEntries();
	nPDBEntriesToT+=nPDBEntries;
	nDBEntriesToT+=nDBEntries;

	//for (Int_t iHit = 0; iHit < nDBEntries; iHit++) {
	// HitData = DataBuffer->GetHit(iHit);
	  /* store volume information */
	// rawStreamTOF->EquipmentId2VolumeId(HitData, Volume);
	//}
	/* reset buffer */
	DataBuffer->Reset();

	/* read data buffer hits */
	for (Int_t iHit = 0; iHit < nPDBEntries; iHit++) {
	  HitData = PackedDataBuffer->GetHit(iHit);
	  /* add volume information */
	  HitData->SetDDLID(iDDL);
	  rawStreamTOF->EquipmentId2VolumeId(HitData, Volume);
	  if (Volume[0]==-1 ||
	      Volume[1]==-1 ||
	      Volume[2]==-1 ||
	      Volume[3]==-1 ||
	      Volume[4]==-1) continue;
	  else {
	    dummy = Volume[3];
	    Volume[3] = Volume[4];
	    Volume[4] = dummy;
	    Int_t tofRaw = (Int_t)((Double_t)HitData->GetTime()*1E3/AliTOFGeometry::TdcBinWidth());
	    Int_t tof;
	    if (!t0Flag) tof = tofRaw;
	    else tof = tofRaw - meantime;
	    Int_t index = geom->GetIndex(Volume);
	    Float_t pos[3];
	    geom->GetPosPar(Volume,pos);
	    Float_t texp = TMath::Sqrt(pos[0]*pos[0]+pos[1]*pos[1]+pos[2]*pos[2])/c*1E9; //expected time in ns
	    Float_t texpBin = texp*1E3/AliTOFGeometry::TdcBinWidth(); //expected time in number of TDC bin
	    Int_t deltabin = tof-TMath::Nint(texpBin);   //to be used with real data; rounding expected time to Int_t
	    htofPartial->Fill(index,deltabin); //channel index start from 0, bin index from 1
	    //debugging printings
	    if (debugFlag > 1) {
		    printf("tofRaw = %i, tof = %i \n",tofRaw,tof);
	    }
	    if (debugFlag > 2) {
		    printf("sector %2d, plate %1d, strip %2d, padz %1d, padx %2d \n",Volume[0],Volume[1],Volume[2],Volume[3],Volume[4]); // too verbose
		    printf("pos x = %f, pos y = %f, pos z = %f \n",pos[0],pos[1],pos[2]); // too verbose
		    printf("expected time = %f (ns)\n",texp); // too verbose
		    printf("expected time bin = %f (TDC bin)\n",texpBin); // too verbose
		    printf("measured time bin = %i (TDC bin) with %f (ns) and ACQ bit = %i \n",tof, HitData->GetTime(), HitData->GetACQ()); // too verbose
		    printf("index = %6d, deltabin = %d , filling index = %6d, and bin = %d\n",index, deltabin, index, deltabin); // too verbose
	    }

	  }
	  /* reset buffer */
	  PackedDataBuffer->Reset();
	}
      }
      //if (debugFlag) {
      //  printf(" Packed Hit Buffer Entries = %i \n",nPDBEntriesToT); // too verbose
      //  printf(" Hit Buffer Entries = %i \n",nDBEntriesToT); // too verbose
      //}

      delete rawReader;
      rawReader = 0x0;
    }
    
    /* free resources */
    free(event);
    
    /* exit when last event received, no need to wait for TERM signal */
    if (eventT==END_OF_RUN) {
      printf("EOR event detected\n");
      break;
    }

  }

  delete rawStreamTOF;
  rawStreamTOF = 0x0;

  delete geom;
  geom = 0x0;

  //write the Run level file   
  TFile * fileRun = new TFile (FILE_RUN,"RECREATE"); 
  htofPartial->Write();
  fileRun->Close();

  //write the Total file
  TH2S *htoftot = 0x0;
  TFile * filetot = 0x0;
  Bool_t isThere=kFALSE;
  const char *dirname = "./";
  TString filename = FILE_TOTAL;
  if((gSystem->FindFile(dirname,filename))!=NULL){
    isThere=kTRUE;
    printf("%s found \n",FILE_TOTAL);
  }
  if (isThere) {

    TFile * filetot1 = new TFile (FILE_TOTAL,"READ"); 
    //look for the file
    if (!filetot1->IsZombie()){
      printf("updating file %s \n",FILE_TOTAL);
      TIter next(filetot1->GetListOfKeys());
      TKey *key;
      //look for the histogram
      while ((key=(TKey*)next())){
	const char * namekey = key->GetName();
	if (strcmp(namekey,"htoftot")==0) {
	  printf(" histo found \n");
	  htoftot = (TH2S*) filetot1->Get("htoftot");
	  htoftot->AddDirectory(0);
	  htoftot->Add(htofPartial);
	  break;
	}
      }
    }
    filetot1->Close();
    delete filetot1;
    filetot1=0x0;
  }
  else {
    printf(" no %s file found \n",FILE_TOTAL);
    htoftot = new TH2S(*htofPartial);
    htoftot->SetName("htoftot");
    htoftot->AddDirectory(0);
  }

  filetot = new TFile (FILE_TOTAL,"RECREATE");
  filetot->cd();
  htoftot->Write();
  filetot->Close();
  
  delete fileRun;
  delete filetot;
  delete htofPartial;
  delete htoftot;

  fileRun = 0x0;
  filetot = 0x0;
  htofPartial = 0x0;
  htoftot = 0x0;

  /* write report */
  printf("Run #%s, received %d physics events out of %d\n",
	 getenv("DATE_RUN_NUMBER"),nevents_physics,nevents_total);

  status = 0;

  /* export file to FXS */
  if (daqDA_FES_storeFile(FILE_RUN, "RUNLevel"))
    status=-2;
  if (daqDA_FES_storeFile(FILE_TOTAL, "DELAYS"))
    status=-2;

  return status;
}
Commit	Line	Data
7b81ffe5	1	/*
7b81ffe5	2
cee5b064	3	TOF DA for online calibration
7b81ffe5	4
78c742e3	5	Contact: Chiara.Zampolli@bo.infn.it
	6	Link: www.bo.infn.it/~zampolli
	7	Run Type: PHYSICS
	8	DA Type: MON
	9	Number of events needed: depending on the run, being run-level
	10	Input Files: TOFdaTotal.root, to be updated if existing
	11	Output Files: TOFdaRun.root, TOFdaTotal.root, both to be exported to the DAQ FXS
	12	Trigger types used: PHYSICS_EVENT
dd635345	13
cee5b064	14	*/
7b81ffe5	15
cee5b064	16	#define FILE_TOTAL "TOFdaTotal.root"
cee5b064	17	#define FILE_RUN "TOFdaRun.root"
7b81ffe5	18
cee5b064	19	// DATE
67cb6c02	20	#include "event.h"
	21	#include "monitor.h"
	22	#include "daqDA.h"
	23
	24	#include <stdio.h>
	25	#include <stdlib.h>
7b81ffe5	26
7b81ffe5	27	//AliRoot
cee5b064	28	#include <AliTOFRawStream.h>
	29	#include <AliRawReaderDate.h>
	30	#include <AliRawReader.h>
	31	#include <AliTOFGeometry.h>
cee5b064	32	#include <AliT0RawReader.h>
	33	#include <AliDAQ.h>
	34	#include <AliTOFHitData.h>
	35	#include <AliTOFHitDataBuffer.h>
59f6cc12	36	#include <AliTOFDaConfigHandler.h>
7b81ffe5	37
7b81ffe5	38	//ROOT
cee5b064	39	#include <TFile.h>
	40	#include <TKey.h>
	41	#include <TH2S.h>
	42	#include <TObject.h>
	43	#include <TMath.h>
	44	#include <TSystem.h>
3bc66139	45	#include "TROOT.h"
3bc66139	46	#include "TPluginManager.h"
67cb6c02	47	#include "TSAXParser.h"
7b81ffe5	48
	49	/* Main routine
	50	Arguments:
	51	1- monitoring data source
	52	*/
	53	int main(int argc, char **argv) {
67cb6c02	54
3bc66139	55	/* magic line from Rene */
	56	gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
	57	"*",
	58	"TStreamerInfo",
	59	"RIO",
	60	"TStreamerInfo()");
	61
7b81ffe5	62	AliTOFGeometry * geom = new AliTOFGeometry();
7b81ffe5	63
7e561cb0	64	static const Int_t size = AliTOFGeometry::NPadXSector()*AliTOFGeometry::NSectors();
7b81ffe5	65	static const Int_t nbins = 500;
7b81ffe5	66	static const Int_t binmin = -20;
67cb6c02	67	const Float_t c = 2.99792458E10; //speed of light [cm/s]
7b81ffe5	68	TH1F::AddDirectory(0);
67cb6c02	69	TH2S * htofPartial = new TH2S("htof","histo with delays",
	70	size,-0.5,size*1.-0.5,
	71	nbins,binmin-0.5,nbins*1.+binmin-0.5);
7b81ffe5	72
7b81ffe5	73	int status;
67cb6c02	74
	75	/* log start of process */
	76	printf("TOF DA started\n");
	77
	78	/* check that we got some arguments = list of files */
7b81ffe5	79	if (argc!=2) {
	80	printf("Wrong number of arguments\n");
	81	return -1;
	82	}
	83
67cb6c02	84	/* retrieve config file */
59f6cc12	85	int getConfigFile = daqDA_DB_getFile("TOFPhysicsConfig.xml","TOFPhysicsConfig.xml");
67cb6c02	86	if (getConfigFile != 0){
	87	printf("Failed to retrieve config file from DB! returning...\n");
	88	return -1;
	89	}
	90
59f6cc12	91	AliTOFDaConfigHandler* tofHandler = new AliTOFDaConfigHandler();
67cb6c02	92	TSAXParser *parser = new TSAXParser();
59f6cc12	93	parser->ConnectToHandler("AliTOFDaConfigHandler", tofHandler);
59f6cc12	94	if (parser->ParseFile("./TOFPhysicsConfig.xml") != 0) {
67cb6c02	95	printf("Failed parsing config file! retunring... \n");
	96	return -1;
	97	}
	98
	99	Int_t debugFlag = tofHandler->GetDebugFlag();
4fc590b9	100	printf("The debug flag is %i\n",debugFlag);
59f6cc12	101	Int_t t0Flag = tofHandler->GetT0Flag();
4fc590b9	102	printf("The T0 flag is %i\n\n",t0Flag);
	103	if (t0Flag) {
	104	printf("The T0 time will be subtracted from the measured TOF time. So, in TDC bins: \n");
	105	printf("tof = tofRaw - (rawReaderT0->GetData(51,0)+rawReaderT0->GetData(52,0))/2) \n\n");
	106	}
	107	else {
	108	printf("The T0 time will not be used.\n");
	109	printf("tof = tofRaw \n\n");
	110	}
67cb6c02	111
	112	/* init some counters */
	113	int nevents_physics=0;
	114	int nevents_total=0;
	115
	116	Int_t iev=0;
	117
	118	Int_t nPDBEntriesToT = 0;
	119	Int_t nDBEntriesToT = 0;
	120	AliTOFHitData *HitData;
	121	Int_t dummy = -1;
	122	Int_t Volume[5];
	123	for (Int_t i=0;i<5;i++) Volume[i]=-1;
	124	AliTOFRawStream *rawStreamTOF = new AliTOFRawStream();
6ae3229a	125	AliTOFHitDataBuffer *DataBuffer;
6ae3229a	126	AliTOFHitDataBuffer *PackedDataBuffer;
67cb6c02	127	Int_t nDBEntries = 0;
	128	Int_t nPDBEntries = 0;
	129
	130	struct eventHeaderStruct *event;
	131	eventTypeType eventT;
	132
7b81ffe5	133	/* define data source : this is argument 1 */
	134	status=monitorSetDataSource( argv[1] );
	135	if (status!=0) {
	136	printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
	137	return -1;
	138	}
	139
	140	/* declare monitoring program */
	141	status=monitorDeclareMp( __FILE__ );
	142	if (status!=0) {
	143	printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
	144	return -1;
	145	}
	146
	147	/* define wait event timeout - 1s max */
	148	monitorSetNowait();
	149	monitorSetNoWaitNetworkTimeout(1000);
	150
7b81ffe5	151	/* main loop (infinite) */
	152	for(;;) {
	153
	154	/* check shutdown condition */
67cb6c02	155	if (daqDA_checkShutdown()) break;
7b81ffe5	156
	157	/* get next event (blocking call until timeout) */
	158	status=monitorGetEventDynamic((void **)&event);
	159	if (status==MON_ERR_EOF) {
	160	printf ("End of File detected\n");
	161	break; /* end of monitoring file has been reached */
	162	}
	163
	164	if (status!=0) {
	165	printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
	166	break;
	167	}
	168
	169	/* retry if got no event */
67cb6c02	170	if (event==NULL) continue;
7b81ffe5	171
	172	iev++;
	173
	174	/* use event - here, just write event id to result file */
cee5b064	175	nevents_total++;
7b81ffe5	176	eventT=event->eventType;
67cb6c02	177	switch (event->eventType) {
7b81ffe5	178
	179	/* START OF RUN */
	180	case START_OF_RUN:
	181	break;
	182	/* END START OF RUN */
	183
67cb6c02	184	/* END OF RUN */
7b81ffe5	185	case END_OF_RUN:
7b81ffe5	186	break;
67cb6c02	187	/* END END OF RUN */
7b81ffe5	188
7b81ffe5	189	case PHYSICS_EVENT:
cee5b064	190	nevents_physics++;
7b81ffe5	191	AliRawReader rawReader = new AliRawReaderDate((void)event);
cee5b064	192	//rawReader->RequireHeader(kFALSE);
7b81ffe5	193
cee5b064	194	//T0 event
cee5b064	195	Int_t meantime = 0;
719903f3	196	AliT0RawReader *rawReaderT0 = new AliT0RawReader(rawReader,kTRUE);
cee5b064	197	if (!rawReaderT0->Next()) {
cee5b064	198	printf("T0: no raw data found!\n");
dd635345	199	}
	200	else {
	201	/*
	202	Int_t allData[105][5];
	203	for (Int_t i=0; i<105; i++) {
cee5b064	204	allData[i][0]=rawReaderT0->GetData(i,0);
dd635345	205	}
cee5b064	206	meantime = allData[49][0];
dd635345	207	*/
d40e79ac	208	//meantime = rawReaderT0->GetData(49,0); //OLD
67cb6c02	209	meantime = (Int_t)((rawReaderT0->GetData(51,0)+rawReaderT0->GetData(52,0))/2.); //Alla
59f6cc12	210	if (debugFlag > 0) {
4fc590b9	211	printf("\nT0 for the current event:\n"); // debugging purpose
59f6cc12	212	printf("time zero = %i (TDC bin)\n", meantime); // debugging purpose
4fc590b9	213	printf("time zero = %f (ns)\n\n", (Float_t)(meantime)24.41E-3); // debugging purpose
59f6cc12	214	}
7b81ffe5	215	}
7b81ffe5	216
	217	delete rawReaderT0;
	218	rawReaderT0 = 0x0;
	219	rawReader->Reset();
dd635345	220
cee5b064	221	//TOF event
67cb6c02	222	dummy = -1;
	223	for (Int_t ii=0; ii<5; ii++) Volume[ii]=-1;
	224	rawStreamTOF->SetRawReader(rawReader);
	225	//rawReader->ReadHeader();
cee5b064	226	rawStreamTOF->ResetBuffers();
cee5b064	227	rawStreamTOF->DecodeDDL(0, AliDAQ::NumberOfDdls("TOF") - 1,0);
67cb6c02	228	nPDBEntriesToT = 0;
	229	nDBEntriesToT = 0;
	230	for (Int_t iDDL = 0; iDDL < AliDAQ::NumberOfDdls("TOF"); iDDL++) {
	231
cee5b064	232	/* read decoded data */
	233	DataBuffer = rawStreamTOF->GetDataBuffer(iDDL);
	234	PackedDataBuffer = rawStreamTOF->GetPackedDataBuffer(iDDL);
	235
	236	/* get buffer entries */
6ae3229a	237	nDBEntries = DataBuffer->GetEntries();
6ae3229a	238	nPDBEntries = PackedDataBuffer->GetEntries();
cee5b064	239	nPDBEntriesToT+=nPDBEntries;
cee5b064	240	nDBEntriesToT+=nDBEntries;
67cb6c02	241
67cb6c02	242	//for (Int_t iHit = 0; iHit < nDBEntries; iHit++) {
cee5b064	243	// HitData = DataBuffer->GetHit(iHit);
cee5b064	244	/* store volume information */
67cb6c02	245	// rawStreamTOF->EquipmentId2VolumeId(HitData, Volume);
cee5b064	246	//}
cee5b064	247	/* reset buffer */
6ae3229a	248	DataBuffer->Reset();
67cb6c02	249
cee5b064	250	/* read data buffer hits */
67cb6c02	251	for (Int_t iHit = 0; iHit < nPDBEntries; iHit++) {
6ae3229a	252	HitData = PackedDataBuffer->GetHit(iHit);
cee5b064	253	/* add volume information */
	254	HitData->SetDDLID(iDDL);
	255	rawStreamTOF->EquipmentId2VolumeId(HitData, Volume);
	256	if (Volume[0]==-1 \|\|
	257	Volume[1]==-1 \|\|
	258	Volume[2]==-1 \|\|
	259	Volume[3]==-1 \|\|
	260	Volume[4]==-1) continue;
	261	else {
	262	dummy = Volume[3];
	263	Volume[3] = Volume[4];
	264	Volume[4] = dummy;
59f6cc12	265	Int_t tofRaw = (Int_t)((Double_t)HitData->GetTime()*1E3/AliTOFGeometry::TdcBinWidth());
	266	Int_t tof;
	267	if (!t0Flag) tof = tofRaw;
	268	else tof = tofRaw - meantime;
96f01799	269	Int_t index = geom->GetIndex(Volume);
cee5b064	270	Float_t pos[3];
96f01799	271	geom->GetPosPar(Volume,pos);
67cb6c02	272	Float_t texp = TMath::Sqrt(pos[0]pos[0]+pos[1]pos[1]+pos[2]pos[2])/c1E9; //expected time in ns
67cb6c02	273	Float_t texpBin = texp*1E3/AliTOFGeometry::TdcBinWidth(); //expected time in number of TDC bin
cee5b064	274	Int_t deltabin = tof-TMath::Nint(texpBin); //to be used with real data; rounding expected time to Int_t
	275	htofPartial->Fill(index,deltabin); //channel index start from 0, bin index from 1
	276	//debugging printings
59f6cc12	277	if (debugFlag > 1) {
	278	printf("tofRaw = %i, tof = %i \n",tofRaw,tof);
	279	}
	280	if (debugFlag > 2) {
	281	printf("sector %2d, plate %1d, strip %2d, padz %1d, padx %2d \n",Volume[0],Volume[1],Volume[2],Volume[3],Volume[4]); // too verbose
	282	printf("pos x = %f, pos y = %f, pos z = %f \n",pos[0],pos[1],pos[2]); // too verbose
	283	printf("expected time = %f (ns)\n",texp); // too verbose
	284	printf("expected time bin = %f (TDC bin)\n",texpBin); // too verbose
	285	printf("measured time bin = %i (TDC bin) with %f (ns) and ACQ bit = %i \n",tof, HitData->GetTime(), HitData->GetACQ()); // too verbose
	286	printf("index = %6d, deltabin = %d , filling index = %6d, and bin = %d\n",index, deltabin, index, deltabin); // too verbose
	287	}
67cb6c02	288
cee5b064	289	}
cee5b064	290	/* reset buffer */
6ae3229a	291	PackedDataBuffer->Reset();
7b81ffe5	292	}
7b81ffe5	293	}
67cb6c02	294	//if (debugFlag) {
	295	// printf(" Packed Hit Buffer Entries = %i \n",nPDBEntriesToT); // too verbose
	296	// printf(" Hit Buffer Entries = %i \n",nDBEntriesToT); // too verbose
	297	//}
	298
7b81ffe5	299	delete rawReader;
7b81ffe5	300	rawReader = 0x0;
7b81ffe5	301	}
cee5b064	302
7b81ffe5	303	/* free resources */
	304	free(event);
	305
	306	/* exit when last event received, no need to wait for TERM signal */
	307	if (eventT==END_OF_RUN) {
	308	printf("EOR event detected\n");
	309	break;
	310	}
67cb6c02	311
7b81ffe5	312	}
67cb6c02	313
	314	delete rawStreamTOF;
	315	rawStreamTOF = 0x0;
	316
7b81ffe5	317	delete geom;
	318	geom = 0x0;
	319
	320	//write the Run level file
cee5b064	321	TFile * fileRun = new TFile (FILE_RUN,"RECREATE");
7b81ffe5	322	htofPartial->Write();
7b81ffe5	323	fileRun->Close();
7b81ffe5	324
cee5b064	325	//write the Total file
7b81ffe5	326	TH2S *htoftot = 0x0;
cee5b064	327	TFile * filetot = 0x0;
	328	Bool_t isThere=kFALSE;
	329	const char *dirname = "./";
	330	TString filename = FILE_TOTAL;
	331	if((gSystem->FindFile(dirname,filename))!=NULL){
	332	isThere=kTRUE;
	333	printf("%s found \n",FILE_TOTAL);
	334	}
67cb6c02	335	if (isThere) {
cee5b064	336
	337	TFile * filetot1 = new TFile (FILE_TOTAL,"READ");
	338	//look for the file
	339	if (!filetot1->IsZombie()){
	340	printf("updating file %s \n",FILE_TOTAL);
	341	TIter next(filetot1->GetListOfKeys());
	342	TKey *key;
	343	//look for the histogram
	344	while ((key=(TKey*)next())){
	345	const char * namekey = key->GetName();
67cb6c02	346	if (strcmp(namekey,"htoftot")==0) {
cee5b064	347	printf(" histo found \n");
	348	htoftot = (TH2S*) filetot1->Get("htoftot");
	349	htoftot->AddDirectory(0);
	350	htoftot->Add(htofPartial);
	351	break;
	352	}
7b81ffe5	353	}
7b81ffe5	354	}
cee5b064	355	filetot1->Close();
	356	delete filetot1;
	357	filetot1=0x0;
7b81ffe5	358	}
7b81ffe5	359	else {
cee5b064	360	printf(" no %s file found \n",FILE_TOTAL);
7b81ffe5	361	htoftot = new TH2S(*htofPartial);
	362	htoftot->SetName("htoftot");
	363	htoftot->AddDirectory(0);
	364	}
cee5b064	365
cee5b064	366	filetot = new TFile (FILE_TOTAL,"RECREATE");
7b81ffe5	367	filetot->cd();
	368	htoftot->Write();
	369	filetot->Close();
	370
	371	delete fileRun;
	372	delete filetot;
7b81ffe5	373	delete htofPartial;
	374	delete htoftot;
	375
	376	fileRun = 0x0;
	377	filetot = 0x0;
7b81ffe5	378	htofPartial = 0x0;
	379	htoftot = 0x0;
	380
	381	/* write report */
67cb6c02	382	printf("Run #%s, received %d physics events out of %d\n",
67cb6c02	383	getenv("DATE_RUN_NUMBER"),nevents_physics,nevents_total);
7b81ffe5	384
67cb6c02	385	status = 0;
7b81ffe5	386
cee5b064	387	/* export file to FXS */
67cb6c02	388	if (daqDA_FES_storeFile(FILE_RUN, "RUNLevel"))
cee5b064	389	status=-2;
67cb6c02	390	if (daqDA_FES_storeFile(FILE_TOTAL, "DELAYS"))
cee5b064	391	status=-2;
7b81ffe5	392
	393	return status;
	394	}