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

/*
PMD DA for online calibration

contact: basanta@iitb.ac.in, Satyajit.Jena@cern.ch
Link: https://twiki.cern.ch/twiki/bin/view/ALICE/DA
Reference run:/afs/cern.ch/user/s/sjena/public/run83496.raw
Run Type: PHYSICS
DA Type: MON
Number of events needed: 1 million for PB+PB, 200 milion for p+p
Input Files: PMD/Calib/Mapping/Run0_999999999_v0_s0.root, PMD_PED.root, PMD_GAIN_CONFIGFILE, pmd_gain_tempfile.dat
Output Files: PMDGAINS.root, to be exported to the DAQ FES
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 "AliPMDCalibGain.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()");


    int status = 0;


    Int_t filestatus = -1, xvar = 5;
    Int_t totevt = -1, maxevt = -1;
    Int_t hotevtsize = -1;
    Bool_t hotfilestatus = false;

    // Reads the pedestal file and keep the values in memory for subtraction

    AliPMDCalibGain *calibgain = new AliPMDCalibGain();

    // Fetch the pedestal file - PMD_PED.root 

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

    if(!status)
      {
	printf("*** Pedestal file retrieved from DB *** \n");
      }
    else
      {
	printf("*** Pedestal file NOT retrieved from DB *** \n");
	return -1;
      }
    
    Int_t pstatus = calibgain->ExtractPedestal("PMD_PED.root");

    if(pstatus == -3) return -3;

    TTree *ic    = NULL;
    TTree *meanc = NULL;

    // Retrieve the PMD_GAIN_CONFIGFILE
    status = daqDA_DB_getFile("PMD_GAIN_CONFIGFILE","PMD_GAIN_CONFIGFILE");

    FILE *fp1 = NULL;

    fp1 = fopen("PMD_GAIN_CONFIGFILE","r");

    if (fp1 == NULL)
      {
	printf("*** PMD GAIN Configfile doesn't exist,Provide one ***\n");
	return -1;
      }
    else
      {
	fscanf(fp1,"%d %d %d %d %d\n",&filestatus, &xvar, &totevt, &maxevt, &hotevtsize);
	//printf("%d %d %d %d %d\n",filestatus, xvar, totevt, maxevt, hotevtsize);
      }
    fclose(fp1);

    
    if (filestatus == 1)
      {
	// Retrieve the Temporray ascii file from DB
	status = daqDA_DB_getFile("pmd_gain_tempfile.dat","pmd_gain_tempfile.dat");
	if(!status)
	  {
	    calibgain->ReadTempFile("pmd_gain_tempfile.dat");
	  }
	else
	  {
	    printf("--- pmd_gain_tempfile.dat: not retrieved from DB ---\n");
	  }
	// Retrieve the hot cell file from DB - PMD_HOT.root
	status = daqDA_DB_getFile("PMD_HOT.root","PMD_HOT.root");
	if(!status)
	  {
	    calibgain->ExtractHotChannel("PMD_HOT.root");
	  }
	else
	  {
	    printf("--- pmd_gain_tempfile.dat: not retrieved from DB ---\n");
	  }
      }


    // decoding the events
    

    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 GAIN DA - strted generating the gain of a cell\n");  
    
    /* init some counters */
    int nevents_physics=0;
    int nevents_total=0;
    
    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");
      printf("Run Number set to Zero \n");
      //return -1;
    }
    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 */
	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;
		
	    case PHYSICS_EVENT:
		nevents_physics++;
		totevt++;
		//if(nevents_physics%100 == 0)printf("Physis Events = %d\n",nevents_physics);
		AliRawReader *rawReader = new AliRawReaderDate((void*)event);
		TObjArray *pmdddlcont = new TObjArray();
		calibgain->ProcessEvent(rawReader, pmdddlcont);

		if (totevt%hotevtsize == 0) hotfilestatus = true;
		delete pmdddlcont;
		pmdddlcont = 0x0;
		delete rawReader;
		rawReader = 0x0;
		
	}
       
	/* free resources */
	free(event);
	
    }

    /* exit when last event received, no need to wait for TERM signal */


    ic    = new TTree("ic","PMD Gain tree");
    meanc = new TTree("meanc","PMD Module mean tree");

    if (filestatus == 0)
      {
	TFile *hotRun = new TFile ("PMD_HOT.root","RECREATE");

	TTree *hot = new TTree("hot","PMD Hot cell tree");
	
	calibgain->FindHotCell(hot,xvar);
	
	hot->Write();
	hotRun->Close();

	// store the hot cell root file in the DB

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

	// store the hot cell root file in the file exchange server

	printf("root file for hot cell is created and getting exported\n");
	status = daqDA_FES_storeFile("PMD_HOT.root","PMD_HOT.root");
      }

    if (hotfilestatus)
      {
	TFile *hotRun = new TFile ("PMD_HOT.root","RECREATE");

	TTree *hot = new TTree("hot","PMD Hot cell tree");
	
	calibgain->FindHotCell(hot,xvar);
	
	hot->Write();
	hotRun->Close();

	// store the hot cell root file in the DB

	status = daqDA_DB_storeFile("PMD_HOT.root","PMD_HOT.root");
	
	// store the hot cell root file in the file exchange server

	printf("root file for hot cell is created and getting exported\n");
	status = daqDA_FES_storeFile("PMD_HOT.root","PMD_HOT.root");
      }


    fp1 = fopen("PMD_GAIN_CONFIGFILE","w+");

    if (totevt < maxevt)
      {
	printf("-----------------------------------------------\n");
	printf("***  Required Number of Events not reached  ***\n");
	printf("***  Number of Events processed = %d        ***\n",totevt);
	printf("***  Writing the intermediate ASCII file    ***\n");
	printf("-----------------------------------------------\n");

	calibgain->WriteTempFile("pmd_gain_tempfile.dat");

	// Store the Intermediate ascii file in the DB
	status = daqDA_DB_storeFile("pmd_gain_tempfile.dat","pmd_gain_tempfile.dat");

	filestatus = 1;
	fprintf(fp1,"%d %d %d %d %d\n",filestatus,xvar,totevt,maxevt,hotevtsize);
	fclose(fp1);

	// Store the configfile in the DB
	status = daqDA_DB_storeFile("PMD_GAIN_CONFIGFILE","PMD_GAIN_CONFIGFILE");

      }
    else if (totevt >= maxevt)
      {
	printf("-----------------------------------------------\n");
	printf("***  Required Number of Events reached = %d ***\n",totevt);
	printf("***  Writing the PMDGAINS.root file           ***\n");
	printf("-----------------------------------------------\n");

	calibgain->Analyse(ic, meanc);

	TFile * gainRun = new TFile ("PMDGAINS.root","RECREATE"); 
	ic->Write();
	gainRun->Close();

	TFile * meanRun = new TFile ("PMD_MEAN_SM.root","RECREATE"); 
	meanc->Write();
	meanRun->Close();


	filestatus = 0;
	totevt     = 0;
	fprintf(fp1,"%d %d %d %d %d\n",filestatus,xvar,totevt,maxevt,hotevtsize);
	fclose(fp1);

	// Store the configfile in the DB
	status = daqDA_DB_storeFile("PMD_GAIN_CONFIGFILE","PMD_GAIN_CONFIGFILE");
      }
    
    delete ic;
    ic = 0;

    delete meanc;
    meanc = 0;
    
    delete calibgain;

    /* store the result file on FES */
 
    if (filestatus == 0)
      {
	printf("root file for cell gain is created and getting exported\n");
	status = daqDA_FES_storeFile("PMDGAINS.root","PMDGAINS.root");
	printf("root file for normalised means of different modules\n");
	status = daqDA_FES_storeFile("PMD_MEAN_SM.root","PMD_MEAN_SM.root");
      }

    if (hotfilestatus)
      {
      }

    
    if (status) {
      status = -2;
    }


    return status;
}
Commit	Line	Data
e43e472f	1	/*
	2	PMD DA for online calibration
	3
8a09f5aa	4	contact: basanta@iitb.ac.in, Satyajit.Jena@cern.ch
f30b7e6c	5	Link: https://twiki.cern.ch/twiki/bin/view/ALICE/DA
53fa2405	6	Reference run:/afs/cern.ch/user/s/sjena/public/run83496.raw
e43e472f	7	Run Type: PHYSICS
	8	DA Type: MON
	9	Number of events needed: 1 million for PB+PB, 200 milion for p+p
6159ad12	10	Input Files: PMD/Calib/Mapping/Run0_999999999_v0_s0.root, PMD_PED.root, PMD_GAIN_CONFIGFILE, pmd_gain_tempfile.dat
d31b0352	11	Output Files: PMDGAINS.root, to be exported to the DAQ FES
e43e472f	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"
e43e472f	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 "AliPMDCalibGain.h"
ee562d68	30	#include "AliLog.h"
ee562d68	31	#include "AliCDBManager.h"
e43e472f	32
	33	//ROOT
	34	#include "TFile.h"
	35	#include "TH1F.h"
	36	#include "TBenchmark.h"
	37	#include "TTree.h"
	38	#include "TROOT.h"
	39	#include "TPluginManager.h"
ee562d68	40	#include "TSystem.h"
e43e472f	41
	42
	43	/* Main routine
	44	Arguments:
	45	1- monitoring data source
	46	*/
	47	int main(int argc, char **argv) {
	48
	49	/* magic line from Rene */
	50	gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
	51	"*",
	52	"TStreamerInfo",
	53	"RIO",
	54	"TStreamerInfo()");
	55
d31b0352	56
	57	int status = 0;
	58
	59
9d1486e3	60	Int_t filestatus = -1, xvar = 5;
9d1486e3	61	Int_t totevt = -1, maxevt = -1;
b2a5fcdc	62	Int_t hotevtsize = -1;
b2a5fcdc	63	Bool_t hotfilestatus = false;
1da50520	64
	65	// Reads the pedestal file and keep the values in memory for subtraction
	66
b9007181	67	AliPMDCalibGain *calibgain = new AliPMDCalibGain();
d31b0352	68
d31b0352	69	// Fetch the pedestal file - PMD_PED.root
9d1486e3	70
d31b0352	71	status = daqDA_DB_getFile("PMD_PED.root","PMD_PED.root");
	72
	73	if(!status)
	74	{
	75	printf("* Pedestal file retrieved from DB * \n");
	76	}
	77	else
	78	{
	79	printf("* Pedestal file NOT retrieved from DB * \n");
	80	return -1;
	81	}
9d1486e3	82
b9007181	83	Int_t pstatus = calibgain->ExtractPedestal("PMD_PED.root");
1da50520	84
1da50520	85	if(pstatus == -3) return -3;
e43e472f	86
9d1486e3	87	TTree *ic = NULL;
9d1486e3	88	TTree *meanc = NULL;
e43e472f	89
d31b0352	90	// Retrieve the PMD_GAIN_CONFIGFILE
	91	status = daqDA_DB_getFile("PMD_GAIN_CONFIGFILE","PMD_GAIN_CONFIGFILE");
	92
1da50520	93	FILE *fp1 = NULL;
1da50520	94
d31b0352	95	fp1 = fopen("PMD_GAIN_CONFIGFILE","r");
1da50520	96
	97	if (fp1 == NULL)
	98	{
d31b0352	99	printf("* PMD GAIN Configfile doesn't exist,Provide one *\n");
d31b0352	100	return -1;
1da50520	101	}
	102	else
	103	{
b2a5fcdc	104	fscanf(fp1,"%d %d %d %d %d\n",&filestatus, &xvar, &totevt, &maxevt, &hotevtsize);
b2a5fcdc	105	//printf("%d %d %d %d %d\n",filestatus, xvar, totevt, maxevt, hotevtsize);
1da50520	106	}
1da50520	107	fclose(fp1);
9d1486e3	108
e43e472f	109
1da50520	110	if (filestatus == 1)
1da50520	111	{
d31b0352	112	// Retrieve the Temporray ascii file from DB
	113	status = daqDA_DB_getFile("pmd_gain_tempfile.dat","pmd_gain_tempfile.dat");
	114	if(!status)
	115	{
b9007181	116	calibgain->ReadTempFile("pmd_gain_tempfile.dat");
d31b0352	117	}
	118	else
	119	{
	120	printf("--- pmd_gain_tempfile.dat: not retrieved from DB ---\n");
	121	}
9d1486e3	122	// Retrieve the hot cell file from DB - PMD_HOT.root
	123	status = daqDA_DB_getFile("PMD_HOT.root","PMD_HOT.root");
	124	if(!status)
	125	{
b9007181	126	calibgain->ExtractHotChannel("PMD_HOT.root");
9d1486e3	127	}
	128	else
	129	{
	130	printf("--- pmd_gain_tempfile.dat: not retrieved from DB ---\n");
	131	}
1da50520	132	}
1da50520	133
9d1486e3	134
e43e472f	135	// decoding the events
e43e472f	136
e43e472f	137
	138	if (argc!=2) {
	139	printf("Wrong number of arguments\n");
	140	return -1;
	141	}
	142
e43e472f	143
	144	/* define data source : this is argument 1 */
	145	status=monitorSetDataSource( argv[1] );
	146	if (status!=0) {
	147	printf("monitorSetDataSource() failed : %s\n",monitorDecodeError(status));
	148	return -1;
	149	}
	150
	151	/* declare monitoring program */
	152	status=monitorDeclareMp( __FILE__ );
	153	if (status!=0) {
	154	printf("monitorDeclareMp() failed : %s\n",monitorDecodeError(status));
	155	return -1;
	156	}
	157
	158	/* define wait event timeout - 1s max */
	159	monitorSetNowait();
	160	monitorSetNoWaitNetworkTimeout(1000);
	161
	162	/* log start of process */
1da50520	163	printf("PMD GAIN DA - strted generating the gain of a cell\n");
e43e472f	164
	165	/* init some counters */
	166	int nevents_physics=0;
	167	int nevents_total=0;
	168
	169	struct eventHeaderStruct *event;
338c285a	170	eventTypeType eventT = 0;
338c285a	171
e43e472f	172	Int_t iev=0;
ee562d68	173
ee562d68	174	// Get run number
0073c15c	175
28816992	176	int runNr = 0;
28816992	177
ee562d68	178	if (getenv("DATE_RUN_NUMBER")==0) {
ee562d68	179	printf("DATE_RUN_NUMBER not properly set.\n");
28816992	180	printf("Run Number set to Zero \n");
28816992	181	//return -1;
ee562d68	182	}
28816992	183	else
	184	{
	185	runNr = atoi(getenv("DATE_RUN_NUMBER"));
	186	}
ee562d68	187
	188	if (gSystem->AccessPathName("localOCDB/PMD/Calib/Mapping",kFileExists))
	189	{
	190	if (gSystem->mkdir("localOCDB/PMD/Calib/Mapping",kTRUE) != 0)
	191	{
	192	printf("Failed to create directory: localOCDB/PMD/Calib/Mapping");
	193	return -1;
	194	}
	195	}
6159ad12	196	status = daqDA_DB_getFile("PMD/Calib/Mapping/Run0_999999999_v0_s0.root","localOCDB/PMD/Calib/Mapping/Run0_999999999_v0_s0.root");
ee562d68	197	if (status)
ee562d68	198	{
6159ad12	199	printf("Failed to get PMD-Mapping file (PMD/Calib/Mapping/Run0_999999999_v0_s0.root) from DAQdetDB, status=%d\n", status);
ee562d68	200	return -1;
	201	}
	202
	203	// Global initializations
	204	AliLog::SetGlobalLogLevel(AliLog::kError);
	205	AliCDBManager *man = AliCDBManager::Instance();
	206	man->SetDefaultStorage("local://localOCDB");
	207	man->SetRun(runNr);
	208
e43e472f	209
	210	/* main loop (infinite) */
	211	for(;;) {
	212
	213	/* check shutdown condition */
	214	if (daqDA_checkShutdown()) {break;}
	215
	216	/* get next event (blocking call until timeout) */
	217	status=monitorGetEventDynamic((void **)&event);
	218	if (status==MON_ERR_EOF) {
	219	printf ("End of File detected\n");
	220	break; /* end of monitoring file has been reached */
	221	}
	222
	223	if (status!=0) {
	224	printf("monitorGetEventDynamic() failed : %s\n",monitorDecodeError(status));
	225	break;
	226	}
	227
	228	/* retry if got no event */
	229	if (event==NULL) {
	230	continue;
	231	}
	232
	233	iev++;
	234
	235	/* use event - here, just write event id to result file */
	236	nevents_total++;
	237	eventT=event->eventType;
	238	switch (event->eventType){
	239
	240	/* START OF RUN */
	241	case START_OF_RUN:
	242	break;
	243	/* END START OF RUN */
	244
	245	/* END OF RUN */
	246	case END_OF_RUN:
	247	break;
	248
	249	case PHYSICS_EVENT:
	250	nevents_physics++;
b2a5fcdc	251	totevt++;
1da50520	252	//if(nevents_physics%100 == 0)printf("Physis Events = %d\n",nevents_physics);
e43e472f	253	AliRawReader rawReader = new AliRawReaderDate((void)event);
338c285a	254	TObjArray *pmdddlcont = new TObjArray();
b9007181	255	calibgain->ProcessEvent(rawReader, pmdddlcont);
e43e472f	256
b2a5fcdc	257	if (totevt%hotevtsize == 0) hotfilestatus = true;
338c285a	258	delete pmdddlcont;
338c285a	259	pmdddlcont = 0x0;
e43e472f	260	delete rawReader;
	261	rawReader = 0x0;
	262
	263	}
	264
	265	/* free resources */
	266	free(event);
	267
338c285a	268	}
	269
	270	/* exit when last event received, no need to wait for TERM signal */
40a42048	271
9d1486e3	272
	273	ic = new TTree("ic","PMD Gain tree");
	274	meanc = new TTree("meanc","PMD Module mean tree");
	275
	276	if (filestatus == 0)
	277	{
	278	TFile *hotRun = new TFile ("PMD_HOT.root","RECREATE");
	279
	280	TTree *hot = new TTree("hot","PMD Hot cell tree");
	281
b9007181	282	calibgain->FindHotCell(hot,xvar);
9d1486e3	283
	284	hot->Write();
	285	hotRun->Close();
	286
	287	// store the hot cell root file in the DB
	288
	289	status = daqDA_DB_storeFile("PMD_HOT.root","PMD_HOT.root");
b2a5fcdc	290
	291	// store the hot cell root file in the file exchange server
	292
	293	printf("root file for hot cell is created and getting exported\n");
	294	status = daqDA_FES_storeFile("PMD_HOT.root","PMD_HOT.root");
9d1486e3	295	}
9d1486e3	296
b2a5fcdc	297	if (hotfilestatus)
	298	{
	299	TFile *hotRun = new TFile ("PMD_HOT.root","RECREATE");
	300
	301	TTree *hot = new TTree("hot","PMD Hot cell tree");
	302
b9007181	303	calibgain->FindHotCell(hot,xvar);
b2a5fcdc	304
	305	hot->Write();
	306	hotRun->Close();
	307
	308	// store the hot cell root file in the DB
	309
	310	status = daqDA_DB_storeFile("PMD_HOT.root","PMD_HOT.root");
	311
	312	// store the hot cell root file in the file exchange server
	313
	314	printf("root file for hot cell is created and getting exported\n");
	315	status = daqDA_FES_storeFile("PMD_HOT.root","PMD_HOT.root");
	316	}
e43e472f	317
1da50520	318
d31b0352	319	fp1 = fopen("PMD_GAIN_CONFIGFILE","w+");
1da50520	320
	321	if (totevt < maxevt)
	322	{
d31b0352	323	printf("-----------------------------------------------\n");
	324	printf("* Required Number of Events not reached *\n");
	325	printf("* Number of Events processed = %d *\n",totevt);
	326	printf("* Writing the intermediate ASCII file *\n");
	327	printf("-----------------------------------------------\n");
	328
b9007181	329	calibgain->WriteTempFile("pmd_gain_tempfile.dat");
d31b0352	330
	331	// Store the Intermediate ascii file in the DB
	332	status = daqDA_DB_storeFile("pmd_gain_tempfile.dat","pmd_gain_tempfile.dat");
1da50520	333
1da50520	334	filestatus = 1;
b2a5fcdc	335	fprintf(fp1,"%d %d %d %d %d\n",filestatus,xvar,totevt,maxevt,hotevtsize);
9d1486e3	336	fclose(fp1);
d31b0352	337
	338	// Store the configfile in the DB
	339	status = daqDA_DB_storeFile("PMD_GAIN_CONFIGFILE","PMD_GAIN_CONFIGFILE");
	340
1da50520	341	}
	342	else if (totevt >= maxevt)
	343	{
d31b0352	344	printf("-----------------------------------------------\n");
	345	printf("* Required Number of Events reached = %d *\n",totevt);
	346	printf("* Writing the PMDGAINS.root file *\n");
	347	printf("-----------------------------------------------\n");
	348
b9007181	349	calibgain->Analyse(ic, meanc);
1da50520	350
	351	TFile * gainRun = new TFile ("PMDGAINS.root","RECREATE");
	352	ic->Write();
	353	gainRun->Close();
	354
9d1486e3	355	TFile * meanRun = new TFile ("PMD_MEAN_SM.root","RECREATE");
	356	meanc->Write();
	357	meanRun->Close();
	358
	359
1da50520	360	filestatus = 0;
1da50520	361	totevt = 0;
b2a5fcdc	362	fprintf(fp1,"%d %d %d %d %d\n",filestatus,xvar,totevt,maxevt,hotevtsize);
9d1486e3	363	fclose(fp1);
d31b0352	364
	365	// Store the configfile in the DB
	366	status = daqDA_DB_storeFile("PMD_GAIN_CONFIGFILE","PMD_GAIN_CONFIGFILE");
1da50520	367	}
e43e472f	368
338c285a	369	delete ic;
338c285a	370	ic = 0;
9d1486e3	371
	372	delete meanc;
	373	meanc = 0;
1da50520	374
b9007181	375	delete calibgain;
1da50520	376
	377	/* store the result file on FES */
	378
	379	if (filestatus == 0)
	380	{
9d1486e3	381	printf("root file for cell gain is created and getting exported\n");
a5f36c63	382	status = daqDA_FES_storeFile("PMDGAINS.root","PMDGAINS.root");
9d1486e3	383	printf("root file for normalised means of different modules\n");
9d1486e3	384	status = daqDA_FES_storeFile("PMD_MEAN_SM.root","PMD_MEAN_SM.root");
1da50520	385	}
b2a5fcdc	386
	387	if (hotfilestatus)
	388	{
	389	}
	390
9d1486e3	391
1da50520	392	if (status) {
	393	status = -2;
	394	}
338c285a	395
338c285a	396
9d1486e3	397
e43e472f	398	return status;
e43e472f	399	}