[u/mrichter/AliRoot.git] / PMD / PMDPEDESTALda.cxx

/*
PMD DA for online calibration

contact: basanta@phy.iitb.ac.in
Link: https://twiki.cern.ch/twiki/bin/view/ALICE/DA
Reference Run:/afs/cern.ch/user/s/sjena/public/run83496.raw
Run Type: PEDESTAL
DA Type: LDC
Number of events needed: 1000
Input Files: PMD/Calib/Mapping/Run0_999999999_v0_s0.root
Output Files: pmd_ped.root, to be exported to the DAQ FXS, pedestal230*.ped
Trigger types used: PHYSICS_EVENT

*/
extern "C" {
#include <daqDA.h>
}

#include "event.h"
#include "monitor.h"

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

//AliRoot
#include "AliRawReaderDate.h"
#include "AliPMDCalibPedestal.h"
#include "AliLog.h"
#include "AliCDBManager.h"

//ROOT
#include "TFile.h"
#include "TH1F.h"
#include "TBenchmark.h"
#include "TTree.h"
#include "TROOT.h"
#include "TPluginManager.h"
#include "TSystem.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()");

    
    AliPMDCalibPedestal *calibped = new AliPMDCalibPedestal();

    TTree *ped = NULL;

    //    TH1F::AddDirectory(0);
  
      
    // decoding the events
  
    int status;
    
    if (argc!=2) {
	printf("Wrong number of arguments\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("PMD PED DA - started generating mean and rms of each channel\n");  
    
    /* init some counters */


    struct eventHeaderStruct *event;
    eventTypeType eventT = 0;
    Int_t iev=0;

    // Get run number

    int runNr = 0;
    if (getenv("DATE_RUN_NUMBER")==0) {
      printf("DATE_RUN_NUMBER not properly set.\n");
      //return -1;
      printf("Run Number is set to Zero \n");
    }
    else
      {
	runNr = atoi(getenv("DATE_RUN_NUMBER"));
      }

    if (gSystem->AccessPathName("localOCDB/PMD/Calib/Mapping",kFileExists))
      {
	if (gSystem->mkdir("localOCDB/PMD/Calib/Mapping",kTRUE) != 0)
	  {
	    printf("Failed to create directory: localOCDB/PMD/Calib/Mapping");
	    return -1;
	  }
      }
    status = daqDA_DB_getFile("PMD/Calib/Mapping/Run0_999999999_v0_s0.root","localOCDB/PMD/Calib/Mapping/Run0_999999999_v0_s0.root");
    if (status)
      {
	printf("Failed to get PMD-Mapping file (PMD/Calib/Mapping/Run0_999999999_v0_s0.root) from DAQdetDB, status=%d\n", status);
	return -1;
      }

    // Global initializations
    AliLog::SetGlobalLogLevel(AliLog::kError);
    AliCDBManager *man = AliCDBManager::Instance();
    man->SetDefaultStorage("local://localOCDB");
    man->SetRun(runNr);

    
    /* 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;
	switch (event->eventType){
	    
	    /* START OF RUN */
	    case START_OF_RUN:
		break;
		/* END START OF RUN */
		
		/* END OF RUN */
	    case END_OF_RUN:
		break;
		
	    case PHYSICS_EVENT:

		AliRawReader *rawReader = new AliRawReaderDate((void*)event);
		TObjArray *pmdddlcont = new TObjArray();
		calibped->ProcessEvent(rawReader,pmdddlcont);
		
		delete pmdddlcont;
		pmdddlcont = 0;
		delete rawReader;
		rawReader = 0x0;
	}
	
	/* free resources */
	free(event);
	
	/* exit when last event received, no need to wait for TERM signal */

	// either END_OF_RUN or 1000 events

	if (iev == 1000)
	  {
	    printf("-- 1000 pedestal events processed : terminating --\n");
	    printf("-- eventT is set to END_OF_RUN --\n");
	    eventT = END_OF_RUN;
	    break;
	  }

    }

    printf(" Total number of events processed = %i \n",iev);

    ped = new TTree("ped","PMD Pedestal tree");
    
    if (eventT==END_OF_RUN) {
      printf("EOR event detected\n");
      calibped->Analyse(ped);
    }
    
    TFile * pedRun = new TFile ("PMD_PED.root","RECREATE"); 
    ped->Write();
    pedRun->Close();

    delete ped;
    ped = 0;

    delete calibped;

/* store the pedestal file in database */

    status = daqDA_DB_storeFile("PMD_PED.root","PMD_PED.root");

    if (!status)
      {
	printf("--- PEDESTAL root FILE IS IN THE DAQ DATABASE\n");
      }
    else
      {
	printf("--- Storing the root file into the database failed\n");
      }

/* store the pedestal file on FES */
 
    status = daqDA_FES_storeFile("PMD_PED.root","PMD_PED.root");
    if (status) {
	status = -2;
    }


    return status;
}
Commit	Line	Data
	1	/*
	2	PMD DA for online calibration
	3
	4	contact: basanta@phy.iitb.ac.in
	5	Link: https://twiki.cern.ch/twiki/bin/view/ALICE/DA
	6	Reference Run:/afs/cern.ch/user/s/sjena/public/run83496.raw
	7	Run Type: PEDESTAL
	8	DA Type: LDC
	9	Number of events needed: 1000
	10	Input Files: PMD/Calib/Mapping/Run0_999999999_v0_s0.root
	11	Output Files: pmd_ped.root, to be exported to the DAQ FXS, pedestal230*.ped
	12	Trigger types used: PHYSICS_EVENT
	13
	14	*/
	15	extern "C" {
	16	#include <daqDA.h>
	17	}
	18
	19	#include "event.h"
	20	#include "monitor.h"
	21
	22	#include <Riostream.h>
	23	#include <stdio.h>
	24	#include <stdlib.h>
	25
	26	//AliRoot
	27	#include "AliRawReaderDate.h"
	28	#include "AliPMDCalibPedestal.h"
	29	#include "AliLog.h"
	30	#include "AliCDBManager.h"
	31
	32	//ROOT
	33	#include "TFile.h"
	34	#include "TH1F.h"
	35	#include "TBenchmark.h"
	36	#include "TTree.h"
	37	#include "TROOT.h"
	38	#include "TPluginManager.h"
	39	#include "TSystem.h"
	40
	41	/* Main routine
	42	Arguments:
	43	1- monitoring data source
	44	*/
	45	int main(int argc, char **argv) {
	46
	47	/* magic line from Rene */
	48	gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
	49	"*",
	50	"TStreamerInfo",
	51	"RIO",
	52	"TStreamerInfo()");
	53
	54
	55	AliPMDCalibPedestal *calibped = new AliPMDCalibPedestal();
	56
	57	TTree *ped = NULL;
	58
	59	// TH1F::AddDirectory(0);
	60
	61
	62	// decoding the events
	63
	64	int status;
	65
	66	if (argc!=2) {
	67	printf("Wrong number of arguments\n");
	68	return -1;
	69	}
	70
	71
	72	/* define data source : this is argument 1 */
	73	status=monitorSetDataSource( argv[1] );
	74	if (status!=0) {
	75	printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
	76	return -1;
	77	}
	78
	79	/* declare monitoring program */
	80	status=monitorDeclareMp( __FILE__ );
	81	if (status!=0) {
	82	printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
	83	return -1;
	84	}
	85
	86	/* define wait event timeout - 1s max */
	87	monitorSetNowait();
	88	monitorSetNoWaitNetworkTimeout(1000);
	89
	90	/* log start of process */
	91	printf("PMD PED DA - started generating mean and rms of each channel\n");
	92
	93	/* init some counters */
	94
	95
	96	struct eventHeaderStruct *event;
	97	eventTypeType eventT = 0;
	98	Int_t iev=0;
	99
	100	// Get run number
	101
	102	int runNr = 0;
	103	if (getenv("DATE_RUN_NUMBER")==0) {
	104	printf("DATE_RUN_NUMBER not properly set.\n");
	105	//return -1;
	106	printf("Run Number is set to Zero \n");
	107	}
	108	else
	109	{
	110	runNr = atoi(getenv("DATE_RUN_NUMBER"));
	111	}
	112
	113	if (gSystem->AccessPathName("localOCDB/PMD/Calib/Mapping",kFileExists))
	114	{
	115	if (gSystem->mkdir("localOCDB/PMD/Calib/Mapping",kTRUE) != 0)
	116	{
	117	printf("Failed to create directory: localOCDB/PMD/Calib/Mapping");
	118	return -1;
	119	}
	120	}
	121	status = daqDA_DB_getFile("PMD/Calib/Mapping/Run0_999999999_v0_s0.root","localOCDB/PMD/Calib/Mapping/Run0_999999999_v0_s0.root");
	122	if (status)
	123	{
	124	printf("Failed to get PMD-Mapping file (PMD/Calib/Mapping/Run0_999999999_v0_s0.root) from DAQdetDB, status=%d\n", status);
	125	return -1;
	126	}
	127
	128	// Global initializations
	129	AliLog::SetGlobalLogLevel(AliLog::kError);
	130	AliCDBManager *man = AliCDBManager::Instance();
	131	man->SetDefaultStorage("local://localOCDB");
	132	man->SetRun(runNr);
	133
	134
	135	/* main loop (infinite) */
	136	for(;;) {
	137
	138	/* check shutdown condition */
	139	if (daqDA_checkShutdown()) {break;}
	140
	141	/* get next event (blocking call until timeout) */
	142	status=monitorGetEventDynamic((void **)&event);
	143	if (status==MON_ERR_EOF) {
	144	printf ("End of File detected\n");
	145	break; /* end of monitoring file has been reached */
	146	}
	147
	148	if (status!=0) {
	149	printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
	150	break;
	151	}
	152
	153	/* retry if got no event */
	154	if (event==NULL) {
	155	continue;
	156	}
	157
	158	iev++;
	159
	160	/* use event - here, just write event id to result file */
	161	eventT=event->eventType;
	162	switch (event->eventType){
	163
	164	/* START OF RUN */
	165	case START_OF_RUN:
	166	break;
	167	/* END START OF RUN */
	168
	169	/* END OF RUN */
	170	case END_OF_RUN:
	171	break;
	172
	173	case PHYSICS_EVENT:
	174
	175	AliRawReader rawReader = new AliRawReaderDate((void)event);
	176	TObjArray *pmdddlcont = new TObjArray();
	177	calibped->ProcessEvent(rawReader,pmdddlcont);
	178
	179	delete pmdddlcont;
	180	pmdddlcont = 0;
	181	delete rawReader;
	182	rawReader = 0x0;
	183	}
	184
	185	/* free resources */
	186	free(event);
	187
	188	/* exit when last event received, no need to wait for TERM signal */
	189
	190	// either END_OF_RUN or 1000 events
	191
	192	if (iev == 1000)
	193	{
	194	printf("-- 1000 pedestal events processed : terminating --\n");
	195	printf("-- eventT is set to END_OF_RUN --\n");
	196	eventT = END_OF_RUN;
	197	break;
	198	}
	199
	200	}
	201
	202	printf(" Total number of events processed = %i \n",iev);
	203
	204	ped = new TTree("ped","PMD Pedestal tree");
	205
	206	if (eventT==END_OF_RUN) {
	207	printf("EOR event detected\n");
	208	calibped->Analyse(ped);
	209	}
	210
	211	TFile * pedRun = new TFile ("PMD_PED.root","RECREATE");
	212	ped->Write();
	213	pedRun->Close();
	214
	215	delete ped;
	216	ped = 0;
	217
	218	delete calibped;
	219
	220	/* store the pedestal file in database */
	221
	222	status = daqDA_DB_storeFile("PMD_PED.root","PMD_PED.root");
	223
	224	if (!status)
	225	{
	226	printf("--- PEDESTAL root FILE IS IN THE DAQ DATABASE\n");
	227	}
	228	else
	229	{
	230	printf("--- Storing the root file into the database failed\n");
	231	}
	232
	233	/* store the pedestal file on FES */
	234
	235	status = daqDA_FES_storeFile("PMD_PED.root","PMD_PED.root");
	236	if (status) {
	237	status = -2;
	238	}
	239
	240
	241
	242
	243	return status;
	244	}