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

/*

DAcase2.c

This program connects to the DAQ data source passed as argument
and populates local "./result.txt" file with the ids of events received
during the run.

The program exits when being asked to shut down (daqDA_checkshutdown)
or End of Run event.

Messages on stdout are exported to DAQ log system.

DA for ZDC standalon pedestal runs
contact: Chiara.Oppedisano@cern.ch

*/

#include <stdio.h>
#include <Riostream.h>

// DATE
#include <daqDA.h>
#include <event.h>
#include <monitor.h>

//ROOT
#include <TRandom.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TProfile.h>
#include <TF1.h>
#include <TFile.h>

//AliRoot
#include <AliRawReaderDate.h>
#include <AliZDCRawStream.h>


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

  TH1F::AddDirectory(0);
  // --- Histograms for ADC pedestals 
  //     [22 signal channels x 2 gain chains + 2 reference PTMs]
  //
  TH1F *hPed[44], *hPedOutOfTime[44];
  TH2F *hPedCorr[44];
  //
  char namhist1[50], namhist2[50], namhist3[50];
  for(Int_t j=0; j<44; j++){
     if(j<10){
       sprintf(namhist1,"PedZN1_%d",j);
       sprintf(namhist2,"PedZN1OutOfTime_%d",j);
       sprintf(namhist3,"PedCorrZN1_%d",j);
     }
     else if(j>=10 && j<20){
       sprintf(namhist1,"PedZP1_%d",j-10);
       sprintf(namhist2,"PedZP1OutOfTime_%d",j-10);
       sprintf(namhist3,"PedCorrZP1_%d",j-10);
     }
     else if(j>=20 && j<24){
       sprintf(namhist1,"PedZEM_%d",j-20);
       sprintf(namhist2,"PedZEMOutOfTime_%d",j-20);
       sprintf(namhist3,"PedCorrZEM_%d",j-20);
     }
     else if(j>=24 && j<33){
       sprintf(namhist1,"PedZN2_%d",j-24);
       sprintf(namhist2,"PedZN2OutOfTime_%d",j-24);
       sprintf(namhist3,"PedCorrZN2_%d",j-24);
     }
     else if(j>=33 && j<43){
       sprintf(namhist1,"PedZP2_%d",j-33);
       sprintf(namhist2,"PedZP2OutOfTime_%d",j-33);
       sprintf(namhist3,"PedCorrZP2_%d",j-33);
     }
     hPed[j] = new TH1F(namhist1, namhist1, 100,0., 200.);
     hPedOutOfTime[j] = new TH1F(namhist2, namhist2, 100,0., 200.);
     hPedCorr[j] = new TH2F(namhist3,namhist3,100,0.,200.,100,0.,200.);
  }

  int status;
  
  if (argc!=2) {
    printf("Wrong number of arguments\n");
    return -1;
  }


  /* open result file */
  FILE *fp=NULL;
  fp=fopen("./result.txt","a");
  if (fp==NULL) {
    printf("Failed to open file\n");
    return -1;
  }
  

  /* 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);
  

  /* log start of process */
  printf("ZDC PEDESTAL monitoring program started\n");  
  
  /* init some counters */
  int nevents_physics=0;
  int nevents_total=0;

  struct equipmentStruct *equipment;
  int *eventEnd;
  int *eventData;
  int *equipmentEnd;
  int *equipmentData;
  int *equipmentID;

  struct eventHeaderStruct *event;
  eventTypeType eventT;
  Int_t iev=0;
  
  /* 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 */
    eventT=event->eventType;
    
    if(eventT==PHYSICS_EVENT){
      //fprintf(fp,"Run #%lu, event size: %lu, BC:%u, Orbit:%u, Period:%u\n",
      //
      // Initalize raw-data reading and decoding
      AliRawReader *reader = new AliRawReaderDate((void*)event);
      const AliRawDataHeader* header = reader->GetDataHeader();
      if(header) {
         UChar_t message = header->GetL1TriggerMessage();
	 if(message & 0x40000){ // DEDICATED PEDESTAL RUN
	    printf("\t L1 message -> PEDESTAL raw data\n");
	    continue;
	 }
	 else{
	    printf("\t L1 message -> NO PEDESTAL raw data found\n");
	    return -1;
	 }
      }
      //Commented until we won't have Raw Data Header...
      /*else{
         //printf("\t ERROR! No Raw Data Header found!!!\n");
	 //return -1;
      }*/
      //
    
      AliZDCRawStream *rawStreamZDC = new AliZDCRawStream(reader);
      //
      if (!rawStreamZDC->Next()) printf(" \t No raw data found!! \n");
      Int_t counter=0;
      Int_t RawADC[44], RawADCoot[44];
      for(Int_t j=0; j<44; j++){
         RawADC[j]=0;
         RawADCoot[j]=0;
      }
      while(rawStreamZDC->Next()){
        Int_t index=-1;
        if(rawStreamZDC->IsADCDataWord()){
    	  if(rawStreamZDC->GetSector(0)==1 || rawStreamZDC->GetSector(0)==2){ // *** ZN1, ZP1
	    index = 10*(rawStreamZDC->GetSector(0)-1)+rawStreamZDC->GetSector(1)+5*rawStreamZDC->GetADCGain();
	  }
	  else if(rawStreamZDC->GetSector(0)==3){ // *** ZEM 
	    index = 10*(rawStreamZDC->GetSector(0)-1)+(rawStreamZDC->GetSector(1)-1)+2*rawStreamZDC->GetADCGain();
	  }
	  else if(rawStreamZDC->GetSector(0)==4 || rawStreamZDC->GetSector(0)==5){ // *** ZN2, ZP2
	    index = 10*(rawStreamZDC->GetSector(0)-2)+rawStreamZDC->GetSector(1)+5*rawStreamZDC->GetADCGain()+4;
	  }
	  if(counter<44){
	    hPed[index]->Fill(rawStreamZDC->GetADCValue()); 
	    RawADC[counter] = rawStreamZDC->GetADCValue();
	  }
	  else{ 
	    hPedOutOfTime[index]->Fill(rawStreamZDC->GetADCValue());
	    RawADCoot[counter-44] = rawStreamZDC->GetADCValue();
	  }
	  counter++;
         }//IsADCDataWord()
	 //
	 if(counter == 88){ // Last ADC channel
           for(Int_t k=0; k<44; k++){
              hPedCorr[k]->Fill(RawADCoot[k], RawADC[k]);
           }
	 }
       }
       //
       nevents_physics++;

    }
    
    nevents_total++;


    /* 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;
    }
  }
  
  /* Analysis of the histograms */
  //
  FILE *fileShuttle;
  fileShuttle = fopen("ZDCPedestal.dat","w");
  //
  Float_t MeanPed[44], MeanPedWidth[44], 
   	MeanPedOOT[44], MeanPedWidthOOT[44],
   	CorrCoeff0[44], CorrCoeff1[44];
  // --- Out-of-time pedestals
  TF1 *ADCfunc[44];
  for(Int_t i=0; i<44; i++){
     hPed[i]->Fit("gaus","Q");
     ADCfunc[i] = hPed[i]->GetFunction("gaus");
     MeanPed[i] = ADCfunc[i]->GetParameter(1);
     MeanPedWidth[i] = ADCfunc[i]->GetParameter(2);
     fprintf(fileShuttle,"\t%f\t%f\n",MeanPed[i],MeanPedWidth[i]);
     //printf("\t MeanPed[%d] = %f\n",i, MeanPed[i]);
  }
  // --- Out-of-time pedestals
  TF1 *ADCootfunc[44];
  for(Int_t i=0; i<44; i++){
     hPedOutOfTime[i]->Fit("gaus","Q");
     ADCootfunc[i] = hPedOutOfTime[i]->GetFunction("gaus");
     MeanPedOOT[i] = ADCootfunc[i]->GetParameter(1);
     MeanPedWidthOOT[i] = ADCootfunc[i]->GetParameter(2);
     fprintf(fileShuttle,"\t%f\t%f\n",MeanPedOOT[i],MeanPedWidthOOT[i]);
     //printf("\t MeanPedOOT[%d] = %f\n",i, MeanPedOOT[i]);
  }
  //
  // --- Fit of correlations
  TProfile* hPedCorrProf[44];
  TF1 *ffunc[44];
  char namhist4[50];
  for(int i=0;i<44;i++) {
     sprintf(namhist4,"ADCvsOOT%d_Prof",i);
     hPedCorrProf[i] = hPedCorr[i]->ProfileX(namhist4,-1,-1,"S");
     hPedCorrProf[i]->SetName(namhist4);
     hPedCorrProf[i]->Fit("pol1","Q");
     ffunc[i] = hPedCorrProf[i]->GetFunction("pol1");
     CorrCoeff0[i] = ffunc[i]->GetParameter(0);
     CorrCoeff1[i] = ffunc[i]->GetParameter(1);
     fprintf(fileShuttle,"\t%f\t%f\n",CorrCoeff0[i],CorrCoeff1[i]);
     //printf("\t CorrCoeff0[%d] = %f, CorrCoeff1[%d] = %f\n",i, CorrCoeff0[i], i, CorrCoeff1[i]);
  }    
  //						       
  fclose(fileShuttle);
  

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

  /* close result file */
  fclose(fp);


  return status;
}
Commit	Line	Data
f3ca8792	1	/*
	2
	3	DAcase2.c
	4
	5	This program connects to the DAQ data source passed as argument
	6	and populates local "./result.txt" file with the ids of events received
	7	during the run.
	8
	9	The program exits when being asked to shut down (daqDA_checkshutdown)
	10	or End of Run event.
	11
	12	Messages on stdout are exported to DAQ log system.
	13
	14	DA for ZDC standalon pedestal runs
	15	contact: Chiara.Oppedisano@cern.ch
	16
	17	*/
	18
	19	#include <stdio.h>
	20	#include <Riostream.h>
	21
	22	// DATE
f3ca8792	23	#include <daqDA.h>
f3ca8792	24	#include <event.h>
	25	#include <monitor.h>
	26
	27	//ROOT
	28	#include <TRandom.h>
	29	#include <TH1F.h>
	30	#include <TH2F.h>
	31	#include <TProfile.h>
	32	#include <TF1.h>
	33	#include <TFile.h>
	34
	35	//AliRoot
	36	#include <AliRawReaderDate.h>
	37	#include <AliZDCRawStream.h>
	38
	39
	40	/* Main routine
	41	Arguments:
	42	1- monitoring data source
	43	*/
	44	int main(int argc, char **argv) {
	45
	46	TH1F::AddDirectory(0);
	47	// --- Histograms for ADC pedestals
	48	// [22 signal channels x 2 gain chains + 2 reference PTMs]
	49	//
	50	TH1F hPed[44], hPedOutOfTime[44];
	51	TH2F *hPedCorr[44];
	52	//
	53	char namhist1[50], namhist2[50], namhist3[50];
	54	for(Int_t j=0; j<44; j++){
	55	if(j<10){
	56	sprintf(namhist1,"PedZN1_%d",j);
	57	sprintf(namhist2,"PedZN1OutOfTime_%d",j);
	58	sprintf(namhist3,"PedCorrZN1_%d",j);
	59	}
	60	else if(j>=10 && j<20){
	61	sprintf(namhist1,"PedZP1_%d",j-10);
	62	sprintf(namhist2,"PedZP1OutOfTime_%d",j-10);
	63	sprintf(namhist3,"PedCorrZP1_%d",j-10);
	64	}
	65	else if(j>=20 && j<24){
	66	sprintf(namhist1,"PedZEM_%d",j-20);
	67	sprintf(namhist2,"PedZEMOutOfTime_%d",j-20);
	68	sprintf(namhist3,"PedCorrZEM_%d",j-20);
	69	}
	70	else if(j>=24 && j<33){
	71	sprintf(namhist1,"PedZN2_%d",j-24);
	72	sprintf(namhist2,"PedZN2OutOfTime_%d",j-24);
	73	sprintf(namhist3,"PedCorrZN2_%d",j-24);
	74	}
	75	else if(j>=33 && j<43){
	76	sprintf(namhist1,"PedZP2_%d",j-33);
	77	sprintf(namhist2,"PedZP2OutOfTime_%d",j-33);
	78	sprintf(namhist3,"PedCorrZP2_%d",j-33);
	79	}
	80	hPed[j] = new TH1F(namhist1, namhist1, 100,0., 200.);
	81	hPedOutOfTime[j] = new TH1F(namhist2, namhist2, 100,0., 200.);
	82	hPedCorr[j] = new TH2F(namhist3,namhist3,100,0.,200.,100,0.,200.);
	83	}
	84
	85	int status;
	86
	87	if (argc!=2) {
88	printf("Wrong number of arguments\n");
89	return -1;
90	}
91
92
93	/* open result file */
94	FILE *fp=NULL;
95	fp=fopen("./result.txt","a");
96	if (fp==NULL) {
97	printf("Failed to open file\n");
98	return -1;
99	}
100
101
102	/* define data source : this is argument 1 */
103	status = monitorSetDataSource( argv[1] );
104	if(status!=0) {
105	printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
106	return -1;
107	}
108
109
110	/* declare monitoring program */
111	status = monitorDeclareMp( __FILE__ );
112	if (status!=0) {
113	printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
114	return -1;
115	}
116
117
118	/* define wait event timeout - 1s max */
119	monitorSetNowait();
120	monitorSetNoWaitNetworkTimeout(1000);
121
122
123	/* log start of process */
124	printf("ZDC PEDESTAL monitoring program started\n");
125
126	/* init some counters */
127	int nevents_physics=0;
128	int nevents_total=0;
129
130	struct equipmentStruct *equipment;
131	int *eventEnd;
132	int *eventData;
133	int *equipmentEnd;
134	int *equipmentData;
135	int *equipmentID;
1a43342e	136
f3ca8792	137	struct eventHeaderStruct *event;
	138	eventTypeType eventT;
	139	Int_t iev=0;
	140
	141	/* main loop (infinite) */
	142	for(;;) {
	143
	144	/* check shutdown condition */
	145	if (daqDA_checkShutdown()) {break;}
	146
	147	/* get next event (blocking call until timeout) */
	148	status=monitorGetEventDynamic((void **)&event);
	149	if (status==MON_ERR_EOF) {
	150	printf ("End of File detected\n");
	151	break; /* end of monitoring file has been reached */
	152	}
	153
	154	if (status!=0) {
	155	printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
	156	break;
	157	}
	158
	159	/* retry if got no event */
	160	if (event==NULL) {
	161	continue;
	162	}
	163
	164	iev++;
	165
	166	/* use event - here, just write event id to result file */
	167	eventT=event->eventType;
	168
	169	if(eventT==PHYSICS_EVENT){
	170	//fprintf(fp,"Run #%lu, event size: %lu, BC:%u, Orbit:%u, Period:%u\n",
	171	//
	172	// Initalize raw-data reading and decoding
	173	AliRawReader reader = new AliRawReaderDate((void)event);
1a43342e	174	const AliRawDataHeader* header = reader->GetDataHeader();
	175	if(header) {
	176	UChar_t message = header->GetL1TriggerMessage();
	177	if(message & 0x40000){ // DEDICATED PEDESTAL RUN
	178	printf("\t L1 message -> PEDESTAL raw data\n");
	179	continue;
	180	}
	181	else{
	182	printf("\t L1 message -> NO PEDESTAL raw data found\n");
	183	return -1;
	184	}
	185	}
	186	//Commented until we won't have Raw Data Header...
	187	/*else{
	188	//printf("\t ERROR! No Raw Data Header found!!!\n");
	189	//return -1;
	190	}*/
f3ca8792	191	//
1a43342e	192
1a43342e	193	AliZDCRawStream *rawStreamZDC = new AliZDCRawStream(reader);
f3ca8792	194	//
1a43342e	195	if (!rawStreamZDC->Next()) printf(" \t No raw data found!! \n");
f3ca8792	196	Int_t counter=0;
	197	Int_t RawADC[44], RawADCoot[44];
	198	for(Int_t j=0; j<44; j++){
	199	RawADC[j]=0;
	200	RawADCoot[j]=0;
	201	}
	202	while(rawStreamZDC->Next()){
	203	Int_t index=-1;
	204	if(rawStreamZDC->IsADCDataWord()){
	205	if(rawStreamZDC->GetSector(0)==1 \|\| rawStreamZDC->GetSector(0)==2){ // *** ZN1, ZP1
	206	index = 10(rawStreamZDC->GetSector(0)-1)+rawStreamZDC->GetSector(1)+5rawStreamZDC->GetADCGain();
	207	}
	208	else if(rawStreamZDC->GetSector(0)==3){ // *** ZEM
	209	index = 10(rawStreamZDC->GetSector(0)-1)+(rawStreamZDC->GetSector(1)-1)+2rawStreamZDC->GetADCGain();
	210	}
	211	else if(rawStreamZDC->GetSector(0)==4 \|\| rawStreamZDC->GetSector(0)==5){ // *** ZN2, ZP2
	212	index = 10(rawStreamZDC->GetSector(0)-2)+rawStreamZDC->GetSector(1)+5rawStreamZDC->GetADCGain()+4;
	213	}
	214	if(counter<44){
	215	hPed[index]->Fill(rawStreamZDC->GetADCValue());
	216	RawADC[counter] = rawStreamZDC->GetADCValue();
	217	}
	218	else{
	219	hPedOutOfTime[index]->Fill(rawStreamZDC->GetADCValue());
	220	RawADCoot[counter-44] = rawStreamZDC->GetADCValue();
	221	}
	222	counter++;
	223	}//IsADCDataWord()
	224	//
	225	if(counter == 88){ // Last ADC channel
	226	for(Int_t k=0; k<44; k++){
	227	hPedCorr[k]->Fill(RawADCoot[k], RawADC[k]);
	228	}
	229	}
	230	}
	231	//
	232	nevents_physics++;
1a43342e	233
f3ca8792	234	}
	235
	236	nevents_total++;
	237
	238
	239	/* free resources */
	240	free(event);
	241
	242	/* exit when last event received, no need to wait for TERM signal */
	243	if (eventT==END_OF_RUN) {
	244	printf("EOR event detected\n");
	245	break;
	246	}
	247	}
	248
	249	/* Analysis of the histograms */
	250	//
	251	FILE *fileShuttle;
	252	fileShuttle = fopen("ZDCPedestal.dat","w");
	253	//
	254	Float_t MeanPed[44], MeanPedWidth[44],
	255	MeanPedOOT[44], MeanPedWidthOOT[44],
	256	CorrCoeff0[44], CorrCoeff1[44];
	257	// --- Out-of-time pedestals
	258	TF1 *ADCfunc[44];
	259	for(Int_t i=0; i<44; i++){
	260	hPed[i]->Fit("gaus","Q");
	261	ADCfunc[i] = hPed[i]->GetFunction("gaus");
	262	MeanPed[i] = ADCfunc[i]->GetParameter(1);
	263	MeanPedWidth[i] = ADCfunc[i]->GetParameter(2);
	264	fprintf(fileShuttle,"\t%f\t%f\n",MeanPed[i],MeanPedWidth[i]);
	265	//printf("\t MeanPed[%d] = %f\n",i, MeanPed[i]);
	266	}
	267	// --- Out-of-time pedestals
	268	TF1 *ADCootfunc[44];
	269	for(Int_t i=0; i<44; i++){
	270	hPedOutOfTime[i]->Fit("gaus","Q");
	271	ADCootfunc[i] = hPedOutOfTime[i]->GetFunction("gaus");
	272	MeanPedOOT[i] = ADCootfunc[i]->GetParameter(1);
	273	MeanPedWidthOOT[i] = ADCootfunc[i]->GetParameter(2);
	274	fprintf(fileShuttle,"\t%f\t%f\n",MeanPedOOT[i],MeanPedWidthOOT[i]);
	275	//printf("\t MeanPedOOT[%d] = %f\n",i, MeanPedOOT[i]);
	276	}
	277	//
	278	// --- Fit of correlations
	279	TProfile* hPedCorrProf[44];
	280	TF1 *ffunc[44];
	281	char namhist4[50];
	282	for(int i=0;i<44;i++) {
	283	sprintf(namhist4,"ADCvsOOT%d_Prof",i);
	284	hPedCorrProf[i] = hPedCorr[i]->ProfileX(namhist4,-1,-1,"S");
	285	hPedCorrProf[i]->SetName(namhist4);
	286	hPedCorrProf[i]->Fit("pol1","Q");
	287	ffunc[i] = hPedCorrProf[i]->GetFunction("pol1");
	288	CorrCoeff0[i] = ffunc[i]->GetParameter(0);
	289	CorrCoeff1[i] = ffunc[i]->GetParameter(1);
	290	fprintf(fileShuttle,"\t%f\t%f\n",CorrCoeff0[i],CorrCoeff1[i]);
	291	//printf("\t CorrCoeff0[%d] = %f, CorrCoeff1[%d] = %f\n",i, CorrCoeff0[i], i, CorrCoeff1[i]);
	292	}
	293	//
	294	fclose(fileShuttle);
	295
	296
	297	/* write report */
298	fprintf(fp,"Run #%s, received %d physics events out of %d\n",getenv("DATE_RUN_NUMBER"),nevents_physics,nevents_total);
299
300	/* close result file */
301	fclose(fp);
302
303
304	return status;
305	}