[u/mrichter/AliRoot.git] / ZDC / ZDCMAPPINGda.cxx

/*\r
\r
This program reads the DAQ data files passed as argument using the monitoring library.\r
\r
The program reports about its processing progress.\r
\r
Messages on stdout are exported to DAQ log system.\r
\r
DA to write mapping for ADC modules and VME scaler\r
\r
Contact: Chiara.Oppedisano@to.infn.it\r
Link: \r
Run Type: PHYSICS\r
DA Type: MON\r
Number of events needed: at least 50% \r
Input Files:  none\r
Output Files: ZDCChMapping.dat\r
Trigger Types Used: different trigger types are used\r
\r
*/\r
\r
#define MAPDATA_FILE  "ZDCChMapping.dat"\r
#define TDCDATA_FILE  "ZDCTDCCalib.dat"\r
#define TDCHISTO_FILE "ZDCTDCHisto.root"\r
\r
#include <stdio.h>\r
#include <stdlib.h>\r
#include <Riostream.h>\r
#include <signal.h>\r
\r
// DATE\r
#include <event.h>\r
#include <monitor.h>\r
#include <daqDA.h>\r
\r
//ROOT\r
#include <TROOT.h>\r
#include <TPluginManager.h>\r
#include <TH1F.h>\r
#include <TH2F.h>\r
#include <TProfile.h>\r
#include <TF1.h>\r
#include <TFile.h>\r
#include <TFitter.h>\r
#include "TMinuitMinimizer.h"\r
\r
//AliRoot\r
#include <AliRawReaderDate.h>\r
#include <AliRawEventHeaderBase.h>\r
#include <AliZDCRawStream.h>\r
\r
int main(int argc, char **argv) {\r
\r
  // needed for streamer application\r
  gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",\r
                                        "*",\r
                                        "TStreamerInfo",\r
                                        "RIO",\r
                                        "TStreamerInfo()"); \r
  // needed for Minuit plugin\r
  gROOT->GetPluginManager()->AddHandler("ROOT::Math::Minimizer",\r
                                        "Minuit",\r
                                        "TMinuitMinimizer",\r
                                        "Minuit",\r
                                        "TMinuitMinimizer(const char*)");\r
\r
  //TMinuitMinimizer m; \r
  TVirtualFitter::SetDefaultFitter("Minuit");\r
\r
  char *monitor_table[] = { "ALL", "no", "PHY", "yes", "SOD", "all", NULL };\r
  int err = monitorDeclareTable(monitor_table);\r
  if(err){\r
    printf("monitorDeclareTable() failed: %s\n", monitorDecodeError(err));\r
    return -1;\r
  } \r
  \r
  int status=0, nphys=0;\r
  int const kNModules = 10;\r
  int const kNChannels = 24;\r
  int const kNScChannels = 32;  \r
  int const kZDCTDCGeo = 4;\r
  \r
  int itdc=0, iprevtdc=-1, ihittdc=0;\r
  float tdcData[6], tdcL0=-999.;	\r
  for(int ij=0; ij<6; ij++) tdcData[ij]=-999.;\r
  \r
  Int_t iMod=-1;\r
  Int_t modGeo[kNModules], modType[kNModules],modNCh[kNModules];\r
  for(Int_t kl=0; kl<kNModules; kl++){\r
     modGeo[kl]=modType[kl]=modNCh[kl]=0;\r
  }\r
  \r
  Int_t ich=0;\r
  Int_t adcMod[2*kNChannels], adcCh[2*kNChannels], sigCode[2*kNChannels];\r
  Int_t det[2*kNChannels], sec[2*kNChannels];\r
  for(Int_t y=0; y<2*kNChannels; y++){\r
    adcMod[y]=adcCh[y]=sigCode[y]=det[y]=sec[y]=0;\r
  }\r
  \r
  Int_t iScCh=0;\r
  Int_t scMod[kNScChannels], scCh[kNScChannels], scSigCode[kNScChannels];\r
  Int_t scDet[kNScChannels], scSec[kNScChannels];\r
  for(Int_t y=0; y<kNScChannels; y++){\r
    scMod[y]=scCh[y]=scSigCode[y]=scDet[y]=scSec[y]=0;\r
  }\r
      \r
  Int_t itdcCh=0;\r
  Int_t tdcMod[kNScChannels], tdcCh[kNScChannels], tdcSigCode[kNScChannels];\r
  Int_t tdcDet[kNScChannels], tdcSec[kNScChannels];\r
  for(Int_t y=0; y<kNScChannels; y++){\r
    tdcMod[y]=tdcCh[y]=tdcSigCode[y]=tdcDet[y]=tdcSec[y]=-1;\r
  }\r
 \r
  TH1F * hTDC[6]={0x0,0x0,0x0,0x0,0x0,0x0};\r
  for(int it=0; it<6; it++){\r
    if(it==0) hTDC[it] = new TH1F("TDCZEM1","TDC ZEM1",200, -200., 200.);\r
    else if(it==1) hTDC[it] = new TH1F("TDCZEM2","TDC ZEM2",200, -200., 200.);\r
    else if(it==2) hTDC[it] = new TH1F("TDCZNC", "TDC ZNC", 200, -200., 200.);\r
    else if(it==3) hTDC[it] = new TH1F("TDCZPC", "TDC ZPC", 200, -200., 200.);\r
    else if(it==4) hTDC[it] = new TH1F("TDCZNA", "TDC ZNA", 200, -200., 200.);\r
    else if(it==5) hTDC[it] = new TH1F("TDCZPA", "TDC ZPA", 200, -200., 200.);\r
  }\r
  \r
  /* log start of process */\r
  printf("\n ZDC MAPPING program started\n");  \r
  signal(SIGSEGV, SIG_DFL);\r
\r
  /* check that we got some arguments = list of files */\r
  if(argc<2) {\r
    printf("Wrong number of arguments\n");\r
    return -1;\r
  }\r
  \r
  FILE *mapFile4Shuttle;\r
\r
  /* read the data files */\r
  for(int n=1;n<argc;n++){\r
   \r
    status=monitorSetDataSource( argv[n] );\r
    if (status!=0) {\r
      printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));\r
      return -1;\r
    }\r
    \r
    /* declare monitoring program */\r
    status=monitorDeclareMp( __FILE__ );\r
    if (status!=0) {\r
      printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));\r
      return -1;\r
    }\r
    monitorSetNowait();\r
    monitorSetNoWaitNetworkTimeout(1000);\r
\r
\r
    int iev = 0;\r
    //Bool_t sodRead = kFALSE;\r
    //while(!sodRead){\r
    \r
    /* loop on events (infinite) */\r
    for(;;) {\r
\r
      struct eventHeaderStruct *event;\r
      eventTypeType eventT;\r
 \r
      /* check shutdown condition */\r
      if (daqDA_checkShutdown()) {break;}\r
\r
      /* get next event */\r
      status=monitorGetEventDynamic((void **)&event);\r
      if(status==MON_ERR_EOF){\r
        printf ("End of File detected\n");\r
        break; /* end of monitoring file has been reached */\r
      }\r
      if(status!=0) {\r
        printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));\r
        return -1;\r
      }\r
\r
      /* retry if got no event */\r
      if(event==NULL) continue;\r
      \r
      // Initalize raw-data reading and decoding\r
      AliRawReader *reader = new AliRawReaderDate((void*)event);\r
      reader->Reset();\r
      reader->Select("ZDC");\r
      // --- Reading event header\r
      //UInt_t evtype = reader->GetType();\r
      //printf("\t ZDCMAPPINGda -> run # %d\n",reader->GetRunNumber());\r
      //\r
      AliZDCRawStream *rawStreamZDC = new AliZDCRawStream(reader);\r
        \r
	\r
      /* use event - here, just write event id to result file */\r
      eventT=event->eventType;\r
      \r
      if(eventT==START_OF_DATA){\r
	  	\r
	rawStreamZDC->SetSODReading(kTRUE);\r
	\r
	// --------------------------------------------------------\r
	// --- Writing ascii data file for the Shuttle preprocessor\r
        mapFile4Shuttle = fopen(MAPDATA_FILE,"w");\r
	//if(rawStreamZDC->Next()){\r
	  while((rawStreamZDC->Next())){\r
            if(rawStreamZDC->IsHeaderMapping()){ // mapping header\r
	       iMod++;\r
	       modGeo[iMod]  = rawStreamZDC->GetADCModule();\r
	       modType[iMod] = rawStreamZDC->GetModType();\r
	       modNCh[iMod]  = rawStreamZDC->GetADCNChannels();\r
	    }\r
            if(rawStreamZDC->IsChMapping()){ \r
	      if(modType[iMod]==1){ // ADC mapping ----------------------\r
	        adcMod[ich]  = rawStreamZDC->GetADCModFromMap(ich);\r
	        adcCh[ich]   = rawStreamZDC->GetADCChFromMap(ich);\r
	        sigCode[ich] = rawStreamZDC->GetADCSignFromMap(ich);\r
	        det[ich]     = rawStreamZDC->GetDetectorFromMap(ich);\r
	        sec[ich]     = rawStreamZDC->GetTowerFromMap(ich);\r
	        ich++;\r
	      }\r
	      else if(modType[iMod]==2){ // VME scaler mapping -------------\r
	        scMod[iScCh]     = rawStreamZDC->GetScalerModFromMap(iScCh);\r
	        scCh[iScCh]      = rawStreamZDC->GetScalerChFromMap(iScCh);\r
	        scSigCode[iScCh] = rawStreamZDC->GetScalerSignFromMap(iScCh);\r
	        scDet[iScCh]     = rawStreamZDC->GetScDetectorFromMap(iScCh);\r
	        scSec[iScCh]     = rawStreamZDC->GetScTowerFromMap(iScCh);\r
	        iScCh++;\r
	      }\r
	      else if(modType[iMod]==6 && modGeo[iMod]==4){ // ZDC TDC mapping --------------------\r
	        tdcMod[itdcCh]     = rawStreamZDC->GetTDCModFromMap(itdcCh);\r
	        tdcCh[itdcCh]      = rawStreamZDC->GetTDCChFromMap(itdcCh);\r
	        tdcSigCode[itdcCh] = rawStreamZDC->GetTDCSignFromMap(itdcCh);\r
	        itdcCh++;\r
	      }\r
	    }\r
 	  }\r
	  // Writing data on output FXS file\r
	  for(Int_t is=0; is<2*kNChannels; is++){\r
	     fprintf(mapFile4Shuttle,"\t%d\t%d\t%d\t%d\t%d\t%d\n",\r
	       is,adcMod[is],adcCh[is],sigCode[is],det[is],sec[is]);\r
    	     //printf("  Mapping DA -> %d ADC: mod %d ch %d, code %d det %d, sec %d\n",\r
	     //  is,adcMod[is],adcCh[is],sigCode[is],det[is],sec[is]);\r
	  }\r
	  for(Int_t is=0; is<kNScChannels; is++){\r
	     fprintf(mapFile4Shuttle,"\t%d\t%d\t%d\t%d\t%d\t%d\n",\r
	       is,scMod[is],scCh[is],scSigCode[is],scDet[is],scSec[is]);\r
 	     //if(scMod[is]!=-1) printf("  Mapping DA -> %d Scaler: mod %d ch %d, code %d det %d, sec %d\n",\r
	     //  is,scMod[is],scCh[is],scSigCode[is],scDet[is],scSec[is]);\r
	  }\r
	  for(Int_t is=0; is<kNScChannels; is++){\r
	     fprintf(mapFile4Shuttle,"\t%d\t%d\t%d\t%d\n",\r
	       is,tdcMod[is],tdcCh[is],tdcSigCode[is]);\r
 	     //if(tdcMod[is]!=-1) printf("  Mapping DA -> %d TDC: mod %d ch %d, code %d\n",\r
	     //  is,tdcMod[is],tdcCh[is],tdcSigCode[is]);\r
	  }\r
	  for(Int_t is=0; is<kNModules; is++){\r
	     fprintf(mapFile4Shuttle,"\t%d\t%d\t%d\n",\r
	     modGeo[is],modType[is],modNCh[is]);\r
	     //printf("  Mapping DA -> Module mapping: geo %d type %d #ch %d\n",\r
	     //  modGeo[is],modType[is],modNCh[is]);\r
	  }\r
	  \r
	//} //if (rawstream)\r
        fclose(mapFile4Shuttle);\r
      }// SOD event\r
      else if(eventT==PHYSICS_EVENT){ \r
\r
        for(int ij=0; ij<6; ij++) tdcData[ij]=-999.;\r
	rawStreamZDC->SetSODReading(kTRUE);\r
\r
  	// ----- Setting ch. mapping -----\r
	for(int jk=0; jk<2*kNChannels; jk++){\r
	  rawStreamZDC->SetMapADCMod(jk, adcMod[jk]);\r
	  rawStreamZDC->SetMapADCCh(jk, adcCh[jk]);\r
	  rawStreamZDC->SetMapADCSig(jk, sigCode[jk]);\r
	  rawStreamZDC->SetMapDet(jk, det[jk]);\r
	  rawStreamZDC->SetMapTow(jk, sec[jk]);\r
	}\r
	\r
  	while(rawStreamZDC->Next()){\r
         if(rawStreamZDC->GetADCModule()!=kZDCTDCGeo) continue; //skipping ADCs, scalers and trigger cards\r
\r
          if(rawStreamZDC->GetADCModule()==kZDCTDCGeo && rawStreamZDC->IsZDCTDCDatum()==kTRUE){\r
             //\r
	     itdc = rawStreamZDC->GetChannel(); \r
             if(itdc==iprevtdc) ihittdc++;\r
             else ihittdc=0;\r
             iprevtdc=itdc;\r
	     //\r
	     if(((itdc>=8 && itdc<=13) || itdc==15) && ihittdc==0){\r
               if(itdc!=15){\r
	         tdcData[itdc-8] = 0.025*rawStreamZDC->GetZDCTDCDatum();\r
	         //\r
	         //printf("   ev.%d **** TDC%d %1.0f ns  \n",nphys,itdc, tdcData[itdc-8]);\r
	       }\r
	       //\r
	       else if(itdc==15){\r
	          tdcL0 = 0.025*rawStreamZDC->GetZDCTDCDatum();\r
	          //\r
		  //printf("   ev.%d ++++ TDCL0 %1.0f ns  \n",nphys,tdcL0);\r
		  //\r
		  for(int ic=0; ic<6; ic++){\r
		    if(tdcData[ic]!=-999. && tdcL0!=-999.){\r
		      hTDC[ic]->Fill(tdcData[ic]-tdcL0);\r
		      //printf(" ev.%d -> Filling histo%d: %f ns\n",nphys,ic, tdcData[ic]-tdcL0);\r
		    }\r
		  }\r
               }\r
	     }\r
	  }\r
        }\r
	\r
 	nphys++;\r
      \r
        delete rawStreamZDC;\r
        rawStreamZDC = 0x0;	\r
\r
      }//(if PHYSICS_EVENT) \r
      else if(eventT==END_OF_RUN){\r
        printf("End Of Run detected\n");\r
        break;\r
      }\r
      \r
      iev++; \r
\r
      /* free resources */\r
      free(event);\r
 \r
    } // event loop    \r
      \r
  }\r
  printf(" \n # of physics events analyzed = %d\n\n", nphys);\r
\r
  /* Analysis of the histograms */\r
  //\r
  FILE *fileShuttle;\r
  fileShuttle = fopen(TDCDATA_FILE,"w");\r
  //\r
  Float_t xUp=0., xLow=0., deltaX=0;\r
  Int_t binMax=0, nBinsx=0;\r
  Float_t mean[6]={0.,0.,0.,0.,0.,0.}, sigma[6]={0.,0.,0.,0.,0.,0.};\r
  TF1 *fitfun[6]={0x0,0x0,0x0,0x0,0x0,0x0};\r
  for(int k=0; k<6; k++){\r
    if(hTDC[k]->GetEntries()!=0){\r
       binMax = hTDC[k]->GetMaximumBin();\r
       if(binMax<=1){\r
         printf("\n WARNING! Something wrong with det %d histo \n\n", k);\r
         continue;\r
       }\r
       // \r
       xUp = (hTDC[k]->GetXaxis())->GetXmax();\r
       xLow = (hTDC[k]->GetXaxis())->GetXmin();\r
       deltaX = xUp-xLow;\r
       nBinsx = (hTDC[k]->GetXaxis())->GetNbins();\r
       //printf(" xMax = %f\n", xLow+binMax*deltaX/nBinsx);\r
       hTDC[k]->Fit("gaus","Q","",xLow+binMax*deltaX/nBinsx*0.75,xLow+binMax*deltaX/nBinsx*1.25);\r
       fitfun[k] = hTDC[k]->GetFunction("gaus");\r
       mean[k] = (Float_t) (fitfun[k]->GetParameter(1));\r
       sigma[k] = (Float_t) (fitfun[k]->GetParameter(2));\r
       //printf("\t Mean value from fit = %1.2f\n", mean[k]);\r
       //\r
       fprintf(fileShuttle,"\t%f\t%f\n",mean[k], sigma[k]);\r
     }\r
     else printf("  WARNING! TDC histo %d has no entries and can't be processed!\n",k);\r
  }\r
  //						       \r
  fclose(fileShuttle);\r
  \r
  TFile *histofile = new TFile(TDCHISTO_FILE,"RECREATE");\r
  histofile->cd();\r
  for(int k=0; k<6; k++) hTDC[k]->Write();						       \r
  histofile->Close();\r
  //\r
  for(Int_t j=0; j<6; j++) delete hTDC[j];\r
  \r
  /* store the result files on FES */\r
  // [1] File with mapping\r
  status = daqDA_FES_storeFile(MAPDATA_FILE, "MAPPING");\r
  if(status){\r
    printf("Failed to export file : %d\n",status);\r
    return -1;\r
  }\r
  // [2] File with TDC data\r
  status = daqDA_FES_storeFile(TDCDATA_FILE, "TDCDATA");\r
  if(status){\r
    printf("Failed to export TDC data file to DAQ FES\n");\r
    return -1;\r
  }\r
  // [3] File with TDC histos\r
  status = daqDA_FES_storeFile(TDCHISTO_FILE, "TDCHISTOS");\r
  if(status){\r
    printf("Failed to export TDC histos file to DAQ FES\n");\r
    return -1;\r
  }\r
\r
  return status;\r
}\r
Commit	Line	Data
76493b40	1	/*\r
	2	\r
	3	This program reads the DAQ data files passed as argument using the monitoring library.\r
	4	\r
	5	The program reports about its processing progress.\r
	6	\r
	7	Messages on stdout are exported to DAQ log system.\r
	8	\r
	9	DA to write mapping for ADC modules and VME scaler\r
	10	\r
	11	Contact: Chiara.Oppedisano@to.infn.it\r
	12	Link: \r
	13	Run Type: PHYSICS\r
	14	DA Type: MON\r
	15	Number of events needed: at least 50% \r
	16	Input Files: none\r
	17	Output Files: ZDCChMapping.dat\r
	18	Trigger Types Used: different trigger types are used\r
	19	\r
	20	*/\r
	21	\r
	22	#define MAPDATA_FILE "ZDCChMapping.dat"\r
	23	#define TDCDATA_FILE "ZDCTDCCalib.dat"\r
	24	#define TDCHISTO_FILE "ZDCTDCHisto.root"\r
	25	\r
	26	#include <stdio.h>\r
	27	#include <stdlib.h>\r
	28	#include <Riostream.h>\r
	29	#include <signal.h>\r
	30	\r
	31	// DATE\r
	32	#include <event.h>\r
	33	#include <monitor.h>\r
	34	#include <daqDA.h>\r
	35	\r
	36	//ROOT\r
	37	#include <TROOT.h>\r
	38	#include <TPluginManager.h>\r
	39	#include <TH1F.h>\r
	40	#include <TH2F.h>\r
	41	#include <TProfile.h>\r
	42	#include <TF1.h>\r
	43	#include <TFile.h>\r
	44	#include <TFitter.h>\r
	45	#include "TMinuitMinimizer.h"\r
	46	\r
	47	//AliRoot\r
	48	#include <AliRawReaderDate.h>\r
	49	#include <AliRawEventHeaderBase.h>\r
	50	#include <AliZDCRawStream.h>\r
	51	\r
	52	int main(int argc, char **argv) {\r
	53	\r
45b889f8	54	// needed for streamer application\r
76493b40	55	gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",\r
45b889f8	56	"*",\r
	57	"TStreamerInfo",\r
	58	"RIO",\r
	59	"TStreamerInfo()"); \r
	60	// needed for Minuit plugin\r
	61	gROOT->GetPluginManager()->AddHandler("ROOT::Math::Minimizer",\r
	62	"Minuit",\r
	63	"TMinuitMinimizer",\r
	64	"Minuit",\r
	65	"TMinuitMinimizer(const char*)");\r
76493b40	66	\r
45b889f8	67	//TMinuitMinimizer m; \r
76493b40	68	TVirtualFitter::SetDefaultFitter("Minuit");\r
76493b40	69	\r
76493b40	70	char *monitor_table[] = { "ALL", "no", "PHY", "yes", "SOD", "all", NULL };\r
	71	int err = monitorDeclareTable(monitor_table);\r
	72	if(err){\r
	73	printf("monitorDeclareTable() failed: %s\n", monitorDecodeError(err));\r
	74	return -1;\r
	75	} \r
	76	\r
	77	int status=0, nphys=0;\r
	78	int const kNModules = 10;\r
	79	int const kNChannels = 24;\r
	80	int const kNScChannels = 32; \r
616b80bc	81	int const kZDCTDCGeo = 4;\r
76493b40	82	\r
	83	int itdc=0, iprevtdc=-1, ihittdc=0;\r
	84	float tdcData[6], tdcL0=-999.; \r
	85	for(int ij=0; ij<6; ij++) tdcData[ij]=-999.;\r
	86	\r
	87	Int_t iMod=-1;\r
	88	Int_t modGeo[kNModules], modType[kNModules],modNCh[kNModules];\r
	89	for(Int_t kl=0; kl<kNModules; kl++){\r
	90	modGeo[kl]=modType[kl]=modNCh[kl]=0;\r
	91	}\r
	92	\r
	93	Int_t ich=0;\r
	94	Int_t adcMod[2kNChannels], adcCh[2kNChannels], sigCode[2*kNChannels];\r
	95	Int_t det[2kNChannels], sec[2kNChannels];\r
	96	for(Int_t y=0; y<2*kNChannels; y++){\r
	97	adcMod[y]=adcCh[y]=sigCode[y]=det[y]=sec[y]=0;\r
	98	}\r
	99	\r
	100	Int_t iScCh=0;\r
	101	Int_t scMod[kNScChannels], scCh[kNScChannels], scSigCode[kNScChannels];\r
	102	Int_t scDet[kNScChannels], scSec[kNScChannels];\r
	103	for(Int_t y=0; y<kNScChannels; y++){\r
	104	scMod[y]=scCh[y]=scSigCode[y]=scDet[y]=scSec[y]=0;\r
	105	}\r
	106	\r
	107	Int_t itdcCh=0;\r
	108	Int_t tdcMod[kNScChannels], tdcCh[kNScChannels], tdcSigCode[kNScChannels];\r
	109	Int_t tdcDet[kNScChannels], tdcSec[kNScChannels];\r
	110	for(Int_t y=0; y<kNScChannels; y++){\r
	111	tdcMod[y]=tdcCh[y]=tdcSigCode[y]=tdcDet[y]=tdcSec[y]=-1;\r
	112	}\r
	113	\r
616b80bc	114	TH1F * hTDC[6]={0x0,0x0,0x0,0x0,0x0,0x0};\r
45b889f8	115	for(int it=0; it<6; it++){\r
	116	if(it==0) hTDC[it] = new TH1F("TDCZEM1","TDC ZEM1",200, -200., 200.);\r
	117	else if(it==1) hTDC[it] = new TH1F("TDCZEM2","TDC ZEM2",200, -200., 200.);\r
	118	else if(it==2) hTDC[it] = new TH1F("TDCZNC", "TDC ZNC", 200, -200., 200.);\r
	119	else if(it==3) hTDC[it] = new TH1F("TDCZPC", "TDC ZPC", 200, -200., 200.);\r
	120	else if(it==4) hTDC[it] = new TH1F("TDCZNA", "TDC ZNA", 200, -200., 200.);\r
	121	else if(it==5) hTDC[it] = new TH1F("TDCZPA", "TDC ZPA", 200, -200., 200.);\r
76493b40	122	}\r
	123	\r
	124	/* log start of process */\r
	125	printf("\n ZDC MAPPING program started\n"); \r
	126	signal(SIGSEGV, SIG_DFL);\r
	127	\r
	128	/* check that we got some arguments = list of files */\r
45b889f8	129	if(argc<2) {\r
76493b40	130	printf("Wrong number of arguments\n");\r
	131	return -1;\r
	132	}\r
	133	\r
	134	FILE *mapFile4Shuttle;\r
	135	\r
	136	/* read the data files */\r
45b889f8	137	for(int n=1;n<argc;n++){\r
76493b40	138	\r
	139	status=monitorSetDataSource( argv[n] );\r
	140	if (status!=0) {\r
	141	printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));\r
	142	return -1;\r
	143	}\r
	144	\r
	145	/* declare monitoring program */\r
	146	status=monitorDeclareMp( __FILE__ );\r
	147	if (status!=0) {\r
	148	printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));\r
	149	return -1;\r
	150	}\r
	151	monitorSetNowait();\r
	152	monitorSetNoWaitNetworkTimeout(1000);\r
	153	\r
	154	\r
	155	int iev = 0;\r
	156	//Bool_t sodRead = kFALSE;\r
	157	//while(!sodRead){\r
	158	\r
	159	/* loop on events (infinite) */\r
	160	for(;;) {\r
	161	\r
	162	struct eventHeaderStruct *event;\r
	163	eventTypeType eventT;\r
	164	\r
	165	/* check shutdown condition */\r
	166	if (daqDA_checkShutdown()) {break;}\r
	167	\r
	168	/* get next event */\r
	169	status=monitorGetEventDynamic((void **)&event);\r
	170	if(status==MON_ERR_EOF){\r
	171	printf ("End of File detected\n");\r
	172	break; /* end of monitoring file has been reached */\r
	173	}\r
	174	if(status!=0) {\r
	175	printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));\r
	176	return -1;\r
	177	}\r
	178	\r
	179	/* retry if got no event */\r
	180	if(event==NULL) continue;\r
	181	\r
	182	// Initalize raw-data reading and decoding\r
	183	AliRawReader reader = new AliRawReaderDate((void)event);\r
	184	reader->Reset();\r
	185	reader->Select("ZDC");\r
	186	// --- Reading event header\r
	187	//UInt_t evtype = reader->GetType();\r
	188	//printf("\t ZDCMAPPINGda -> run # %d\n",reader->GetRunNumber());\r
	189	//\r
	190	AliZDCRawStream *rawStreamZDC = new AliZDCRawStream(reader);\r
	191	\r
	192	\r
	193	/* use event - here, just write event id to result file */\r
	194	eventT=event->eventType;\r
	195	\r
	196	if(eventT==START_OF_DATA){\r
	197	\r
	198	rawStreamZDC->SetSODReading(kTRUE);\r
	199	\r
	200	// --------------------------------------------------------\r
	201	// --- Writing ascii data file for the Shuttle preprocessor\r
202	mapFile4Shuttle = fopen(MAPDATA_FILE,"w");\r
203	//if(rawStreamZDC->Next()){\r
204	while((rawStreamZDC->Next())){\r
205	if(rawStreamZDC->IsHeaderMapping()){ // mapping header\r
206	iMod++;\r
207	modGeo[iMod] = rawStreamZDC->GetADCModule();\r
208	modType[iMod] = rawStreamZDC->GetModType();\r
209	modNCh[iMod] = rawStreamZDC->GetADCNChannels();\r
210	}\r
211	if(rawStreamZDC->IsChMapping()){ \r
212	if(modType[iMod]==1){ // ADC mapping ----------------------\r
213	adcMod[ich] = rawStreamZDC->GetADCModFromMap(ich);\r
214	adcCh[ich] = rawStreamZDC->GetADCChFromMap(ich);\r
215	sigCode[ich] = rawStreamZDC->GetADCSignFromMap(ich);\r
216	det[ich] = rawStreamZDC->GetDetectorFromMap(ich);\r
217	sec[ich] = rawStreamZDC->GetTowerFromMap(ich);\r
218	ich++;\r
219	}\r
220	else if(modType[iMod]==2){ // VME scaler mapping -------------\r
221	scMod[iScCh] = rawStreamZDC->GetScalerModFromMap(iScCh);\r
222	scCh[iScCh] = rawStreamZDC->GetScalerChFromMap(iScCh);\r
223	scSigCode[iScCh] = rawStreamZDC->GetScalerSignFromMap(iScCh);\r
224	scDet[iScCh] = rawStreamZDC->GetScDetectorFromMap(iScCh);\r
225	scSec[iScCh] = rawStreamZDC->GetScTowerFromMap(iScCh);\r
226	iScCh++;\r
227	}\r
228	else if(modType[iMod]==6 && modGeo[iMod]==4){ // ZDC TDC mapping --------------------\r
229	tdcMod[itdcCh] = rawStreamZDC->GetTDCModFromMap(itdcCh);\r
230	tdcCh[itdcCh] = rawStreamZDC->GetTDCChFromMap(itdcCh);\r
231	tdcSigCode[itdcCh] = rawStreamZDC->GetTDCSignFromMap(itdcCh);\r
232	itdcCh++;\r
233	}\r
234	}\r
235	}\r
236	// Writing data on output FXS file\r
237	for(Int_t is=0; is<2*kNChannels; is++){\r
238	fprintf(mapFile4Shuttle,"\t%d\t%d\t%d\t%d\t%d\t%d\n",\r
239	is,adcMod[is],adcCh[is],sigCode[is],det[is],sec[is]);\r
240	//printf(" Mapping DA -> %d ADC: mod %d ch %d, code %d det %d, sec %d\n",\r
241	// is,adcMod[is],adcCh[is],sigCode[is],det[is],sec[is]);\r
242	}\r
243	for(Int_t is=0; is<kNScChannels; is++){\r
244	fprintf(mapFile4Shuttle,"\t%d\t%d\t%d\t%d\t%d\t%d\n",\r
245	is,scMod[is],scCh[is],scSigCode[is],scDet[is],scSec[is]);\r
246	//if(scMod[is]!=-1) printf(" Mapping DA -> %d Scaler: mod %d ch %d, code %d det %d, sec %d\n",\r
247	// is,scMod[is],scCh[is],scSigCode[is],scDet[is],scSec[is]);\r
248	}\r
249	for(Int_t is=0; is<kNScChannels; is++){\r
250	fprintf(mapFile4Shuttle,"\t%d\t%d\t%d\t%d\n",\r
251	is,tdcMod[is],tdcCh[is],tdcSigCode[is]);\r
252	//if(tdcMod[is]!=-1) printf(" Mapping DA -> %d TDC: mod %d ch %d, code %d\n",\r
253	// is,tdcMod[is],tdcCh[is],tdcSigCode[is]);\r
254	}\r
255	for(Int_t is=0; is<kNModules; is++){\r
256	fprintf(mapFile4Shuttle,"\t%d\t%d\t%d\n",\r
257	modGeo[is],modType[is],modNCh[is]);\r
258	//printf(" Mapping DA -> Module mapping: geo %d type %d #ch %d\n",\r
259	// modGeo[is],modType[is],modNCh[is]);\r
260	}\r
261	\r
262	//} //if (rawstream)\r
263	fclose(mapFile4Shuttle);\r
264	}// SOD event\r
265	else if(eventT==PHYSICS_EVENT){ \r
266	\r
267	for(int ij=0; ij<6; ij++) tdcData[ij]=-999.;\r
268	rawStreamZDC->SetSODReading(kTRUE);\r
269	\r
270	// ----- Setting ch. mapping -----\r
45b889f8	271	for(int jk=0; jk<2*kNChannels; jk++){\r
76493b40	272	rawStreamZDC->SetMapADCMod(jk, adcMod[jk]);\r
	273	rawStreamZDC->SetMapADCCh(jk, adcCh[jk]);\r
	274	rawStreamZDC->SetMapADCSig(jk, sigCode[jk]);\r
	275	rawStreamZDC->SetMapDet(jk, det[jk]);\r
	276	rawStreamZDC->SetMapTow(jk, sec[jk]);\r
	277	}\r
	278	\r
	279	while(rawStreamZDC->Next()){\r
616b80bc	280	if(rawStreamZDC->GetADCModule()!=kZDCTDCGeo) continue; //skipping ADCs, scalers and trigger cards\r
	281	\r
	282	if(rawStreamZDC->GetADCModule()==kZDCTDCGeo && rawStreamZDC->IsZDCTDCDatum()==kTRUE){\r
76493b40	283	//\r
76493b40	284	itdc = rawStreamZDC->GetChannel(); \r
616b80bc	285	if(itdc==iprevtdc) ihittdc++;\r
	286	else ihittdc=0;\r
	287	iprevtdc=itdc;\r
	288	//\r
45b889f8	289	if(((itdc>=8 && itdc<=13) \|\| itdc==15) && ihittdc==0){\r
616b80bc	290	if(itdc!=15){\r
	291	tdcData[itdc-8] = 0.025*rawStreamZDC->GetZDCTDCDatum();\r
	292	//\r
45b889f8	293	//printf(" ev.%d **** TDC%d %1.0f ns \n",nphys,itdc, tdcData[itdc-8]);\r
616b80bc	294	}\r
	295	//\r
	296	else if(itdc==15){\r
	297	tdcL0 = 0.025*rawStreamZDC->GetZDCTDCDatum();\r
	298	//\r
45b889f8	299	//printf(" ev.%d ++++ TDCL0 %1.0f ns \n",nphys,tdcL0);\r
616b80bc	300	//\r
	301	for(int ic=0; ic<6; ic++){\r
	302	if(tdcData[ic]!=-999. && tdcL0!=-999.){\r
	303	hTDC[ic]->Fill(tdcData[ic]-tdcL0);\r
45b889f8	304	//printf(" ev.%d -> Filling histo%d: %f ns\n",nphys,ic, tdcData[ic]-tdcL0);\r
616b80bc	305	}\r
	306	}\r
	307	}\r
76493b40	308	}\r
	309	}\r
	310	}\r
	311	\r
45b889f8	312	nphys++;\r
	313	\r
	314	delete rawStreamZDC;\r
	315	rawStreamZDC = 0x0; \r
	316	\r
76493b40	317	}//(if PHYSICS_EVENT) \r
	318	else if(eventT==END_OF_RUN){\r
	319	printf("End Of Run detected\n");\r
	320	break;\r
	321	}\r
	322	\r
76493b40	323	iev++; \r
	324	\r
	325	/* free resources */\r
	326	free(event);\r
	327	\r
	328	} // event loop \r
	329	\r
	330	}\r
	331	printf(" \n # of physics events analyzed = %d\n\n", nphys);\r
	332	\r
	333	/* Analysis of the histograms */\r
	334	//\r
	335	FILE *fileShuttle;\r
	336	fileShuttle = fopen(TDCDATA_FILE,"w");\r
	337	//\r
	338	Float_t xUp=0., xLow=0., deltaX=0;\r
	339	Int_t binMax=0, nBinsx=0;\r
616b80bc	340	Float_t mean[6]={0.,0.,0.,0.,0.,0.}, sigma[6]={0.,0.,0.,0.,0.,0.};\r
616b80bc	341	TF1 *fitfun[6]={0x0,0x0,0x0,0x0,0x0,0x0};\r
45b889f8	342	for(int k=0; k<6; k++){\r
76493b40	343	if(hTDC[k]->GetEntries()!=0){\r
	344	binMax = hTDC[k]->GetMaximumBin();\r
	345	if(binMax<=1){\r
	346	printf("\n WARNING! Something wrong with det %d histo \n\n", k);\r
	347	continue;\r
	348	}\r
	349	// \r
	350	xUp = (hTDC[k]->GetXaxis())->GetXmax();\r
	351	xLow = (hTDC[k]->GetXaxis())->GetXmin();\r
	352	deltaX = xUp-xLow;\r
	353	nBinsx = (hTDC[k]->GetXaxis())->GetNbins();\r
	354	//printf(" xMax = %f\n", xLow+binMax*deltaX/nBinsx);\r
	355	hTDC[k]->Fit("gaus","Q","",xLow+binMaxdeltaX/nBinsx0.75,xLow+binMaxdeltaX/nBinsx1.25);\r
	356	fitfun[k] = hTDC[k]->GetFunction("gaus");\r
	357	mean[k] = (Float_t) (fitfun[k]->GetParameter(1));\r
	358	sigma[k] = (Float_t) (fitfun[k]->GetParameter(2));\r
	359	//printf("\t Mean value from fit = %1.2f\n", mean[k]);\r
	360	//\r
	361	fprintf(fileShuttle,"\t%f\t%f\n",mean[k], sigma[k]);\r
	362	}\r
	363	else printf(" WARNING! TDC histo %d has no entries and can't be processed!\n",k);\r
	364	}\r
	365	// \r
	366	fclose(fileShuttle);\r
	367	\r
	368	TFile *histofile = new TFile(TDCHISTO_FILE,"RECREATE");\r
	369	histofile->cd();\r
	370	for(int k=0; k<6; k++) hTDC[k]->Write(); \r
	371	histofile->Close();\r
	372	//\r
	373	for(Int_t j=0; j<6; j++) delete hTDC[j];\r
	374	\r
	375	/* store the result files on FES */\r
	376	// [1] File with mapping\r
	377	status = daqDA_FES_storeFile(MAPDATA_FILE, "MAPPING");\r
	378	if(status){\r
	379	printf("Failed to export file : %d\n",status);\r
	380	return -1;\r
	381	}\r
	382	// [2] File with TDC data\r
	383	status = daqDA_FES_storeFile(TDCDATA_FILE, "TDCDATA");\r
	384	if(status){\r
	385	printf("Failed to export TDC data file to DAQ FES\n");\r
	386	return -1;\r
	387	}\r
	388	// [3] File with TDC histos\r
	389	status = daqDA_FES_storeFile(TDCHISTO_FILE, "TDCHISTOS");\r
	390	if(status){\r
	391	printf("Failed to export TDC histos file to DAQ FES\n");\r
	392	return -1;\r
	393	}\r
	394	\r
	395	return status;\r
	396	}\r