[u/mrichter/AliRoot.git] / prod / LHC08d15 / simrun.C

//#define VERBOSEARGS

{
	// set job and simulation variables as :
	// root run.C  --run <x> --event <y> --process <proc> --minhard <min> --maxhard <max> --minpt <minpt> ...
cout<<"Begin SIMRUN"<<endl;	
	int nrun = 0;
	int nevent = 0;
	int seed = 0;
	
	char sseed[1024];
	char srun[1024];
	char sevent[1024];
	char secms[1024];
	char sprocess[1024];
	char sminpthard[1024];
	char smaxpthard[1024];
	char sminptgammapi0[1024];
	char spi0gammafrag[1024];
	char squench[1024];
	char sqhat[1024];
	char smedlength[1024];
	char strig[1024];
	char syear[1024];
	
	sprintf(srun,"");
	sprintf(sevent,"");
	sprintf(secms,"");
	sprintf(sprocess,"");
	sprintf(sminpthard,"");
	sprintf(smaxpthard,"");
	sprintf(sminptgammapi0,"");
	sprintf(spi0gammafrag,"");
	sprintf(squench,"");
	sprintf(sqhat,"");
	//sprintf(smedlength,"");
	sprintf(strig,"");
	sprintf(syear,"");
	
	for (int i=0; i< gApplication->Argc();i++){
#ifdef VERBOSEARGS
		printf("Arg  %d:  %s\n",i,gApplication->Argv(i));
#endif
		if (!(strcmp(gApplication->Argv(i),"--run")))
			nrun = atoi(gApplication->Argv(i+1));
		sprintf(srun,"%d",nrun);
	
		if (!(strcmp(gApplication->Argv(i),"--event")))
			nevent = atoi(gApplication->Argv(i+1));
		sprintf(sevent,"%d",nevent);
				
		if (!(strcmp(gApplication->Argv(i),"--process")))
			sprintf(sprocess, gApplication->Argv(i+1));
				
		if (!(strcmp(gApplication->Argv(i),"--year")))
			sprintf(syear, gApplication->Argv(i+1));

		if (!(strcmp(gApplication->Argv(i),"--ecms")))
			sprintf(secms,gApplication->Argv(i+1));
						
		if (!(strcmp(gApplication->Argv(i),"--minhard")))
			sprintf(sminpthard,gApplication->Argv(i+1));
							
		if (!(strcmp(gApplication->Argv(i),"--maxhard")))
			sprintf(smaxpthard,gApplication->Argv(i+1));
	
		if (!(strcmp(gApplication->Argv(i),"--pi0gammafrag")))
			sprintf(spi0gammafrag,gApplication->Argv(i+1));
																						
		if (!(strcmp(gApplication->Argv(i),"--minpt")))
			sprintf(sminptgammapi0,gApplication->Argv(i+1));
																	
		if (!(strcmp(gApplication->Argv(i),"--quench")))
			sprintf(squench,gApplication->Argv(i+1));
										
		if (!(strcmp(gApplication->Argv(i),"--qhat")))
			sprintf(sqhat,gApplication->Argv(i+1));
											
		//if (!(strcmp(gApplication->Argv(i),"--medlength")))
		//	sprintf(smedlength,gApplication->Argv(i+1));
											
		if (!(strcmp(gApplication->Argv(i),"--trigger")))
			sprintf(strig,gApplication->Argv(i+1));
													
	}
	cout<<"GOT VARIABLES SIMRUN"<<endl;      
	//rec params for reconstruction
	//only for release < 16
	//char cmd[200] ; 
	//sprintf(cmd, ".! tar zxvf DBpp.tgz") ; 
	//gROOT->ProcessLine(cmd) ; 
  	
	seed = nrun * 100000 + nevent;
	sprintf(sseed,"%d",seed);
	
	if (seed==0) {
		fprintf(stderr,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");
		fprintf(stderr,"!!!!  WARNING! Seeding variable for MC is 0          !!!!\n");
		fprintf(stderr,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");
	} else {
		fprintf(stdout,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");
		fprintf(stdout,"!!!  MC Seed is %d \n",seed);
		fprintf(stdout,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");
	}
	
	// set the seed environment variable
	gSystem->Setenv("CONFIG_SEED",sseed);
	gSystem->Setenv("RUN",srun);
	gSystem->Setenv("EVENT",sevent);
	gSystem->Setenv("DC_RUN_TYPE",sprocess);//"kPyGammaJetPHOS");
	gSystem->Setenv("YEAR",syear);//"14000");	
	gSystem->Setenv("ECMS",secms);//"14000");
	gSystem->Setenv("PTHARDMIN",sminpthard);//"20");
	gSystem->Setenv("PTHARDMAX",smaxpthard);//"30");
	gSystem->Setenv("PI0GAMMAINDET",spi0gammafrag);//"1");
	gSystem->Setenv("PTGAMMAPI0MIN",sminptgammapi0);//"1");
	gSystem->Setenv("QUENCHING",squench);
	gSystem->Setenv("QHAT",sqhat);
	//gSystem->Setenv("MEDLENGTH",smedlength);
	gSystem->Setenv("TRIGGER",strig);
	
	cout<< "SIMRUN:: Run = " << gSystem->Getenv("RUN") << "; Event = " << gSystem->Getenv("EVENT") 
	<< "; Process = "    << gSystem->Getenv("DC_RUN_TYPE")
	<< "; E cms = " << gSystem->Getenv("ECMS") 
	<< "; Trigger configuration = "<< gSystem->Getenv("TRIGGER")
	<< "; Year = "<<gSystem->Getenv("YEAR")	<< endl;
	
	cout<< "         minpthard = " << gSystem->Getenv("PTHARDMIN") 
	<< "; maxpthard =" << gSystem->Getenv("PTHARDMAX") 
	<< "; force pi0/gamma frag in detector = "     << gSystem->Getenv("PI0GAMMAINDET") 
	<< "; minpt = "     << gSystem->Getenv("PTGAMMAPI0MIN") 
	<< "; quenching option ="<<gSystem->Getenv("QUENCHING") 
	<< "; qhat = "<<gSystem->Getenv("QHAT") 
	//<< "; medium length = "<<gSystem->Getenv("MEDIUMLENGTH")
	<< endl;
	
	cout<<">>>>> SIMULATION <<<<<"<<endl;
	gSystem->Exec("aliroot -b -q sim.C > sim.log 2>&1");
	cout<<">>>>> RECONSTRUCTION <<<<<"<<endl;
	gSystem->Exec("aliroot -b -q rec.C > rec.log 2>&1");
	cout<<">>>>> TAG <<<<<"<<endl;
	gSystem->Exec("aliroot -b -q tag.C > tag.log 2>&1");
	cout<<">>>>> CHECK ESD <<<<<"<<endl;
	gSystem->Exec("aliroot -b -q CheckESD.C > check.log 2>&1");
	
}
Commit	Line	Data
333b9f18	1	//#define VERBOSEARGS
	2
	3	{
	4	// set job and simulation variables as :
	5	// root run.C --run <x> --event <y> --process <proc> --minhard <min> --maxhard <max> --minpt <minpt> ...
	6	cout<<"Begin SIMRUN"<<endl;
	7	int nrun = 0;
	8	int nevent = 0;
	9	int seed = 0;
	10
	11	char sseed[1024];
	12	char srun[1024];
	13	char sevent[1024];
	14	char secms[1024];
	15	char sprocess[1024];
	16	char sminpthard[1024];
	17	char smaxpthard[1024];
	18	char sminptgammapi0[1024];
	19	char spi0gammafrag[1024];
	20	char squench[1024];
	21	char sqhat[1024];
	22	char smedlength[1024];
	23	char strig[1024];
	24	char syear[1024];
	25
	26	sprintf(srun,"");
	27	sprintf(sevent,"");
	28	sprintf(secms,"");
	29	sprintf(sprocess,"");
	30	sprintf(sminpthard,"");
	31	sprintf(smaxpthard,"");
	32	sprintf(sminptgammapi0,"");
	33	sprintf(spi0gammafrag,"");
	34	sprintf(squench,"");
	35	sprintf(sqhat,"");
	36	//sprintf(smedlength,"");
	37	sprintf(strig,"");
	38	sprintf(syear,"");
	39
	40	for (int i=0; i< gApplication->Argc();i++){
	41	#ifdef VERBOSEARGS
	42	printf("Arg %d: %s\n",i,gApplication->Argv(i));
	43	#endif
	44	if (!(strcmp(gApplication->Argv(i),"--run")))
	45	nrun = atoi(gApplication->Argv(i+1));
	46	sprintf(srun,"%d",nrun);
	47
	48	if (!(strcmp(gApplication->Argv(i),"--event")))
	49	nevent = atoi(gApplication->Argv(i+1));
	50	sprintf(sevent,"%d",nevent);
	51
	52	if (!(strcmp(gApplication->Argv(i),"--process")))
	53	sprintf(sprocess, gApplication->Argv(i+1));
	54
	55	if (!(strcmp(gApplication->Argv(i),"--year")))
	56	sprintf(syear, gApplication->Argv(i+1));
	57
	58	if (!(strcmp(gApplication->Argv(i),"--ecms")))
	59	sprintf(secms,gApplication->Argv(i+1));
	60
	61	if (!(strcmp(gApplication->Argv(i),"--minhard")))
	62	sprintf(sminpthard,gApplication->Argv(i+1));
	63
	64	if (!(strcmp(gApplication->Argv(i),"--maxhard")))
65	sprintf(smaxpthard,gApplication->Argv(i+1));
66
67	if (!(strcmp(gApplication->Argv(i),"--pi0gammafrag")))
68	sprintf(spi0gammafrag,gApplication->Argv(i+1));
69
70	if (!(strcmp(gApplication->Argv(i),"--minpt")))
71	sprintf(sminptgammapi0,gApplication->Argv(i+1));
72
73	if (!(strcmp(gApplication->Argv(i),"--quench")))
74	sprintf(squench,gApplication->Argv(i+1));
75
76	if (!(strcmp(gApplication->Argv(i),"--qhat")))
77	sprintf(sqhat,gApplication->Argv(i+1));
78
79	//if (!(strcmp(gApplication->Argv(i),"--medlength")))
80	// sprintf(smedlength,gApplication->Argv(i+1));
81
82	if (!(strcmp(gApplication->Argv(i),"--trigger")))
83	sprintf(strig,gApplication->Argv(i+1));
84
85	}
86	cout<<"GOT VARIABLES SIMRUN"<<endl;
87	//rec params for reconstruction
88	//only for release < 16
89	//char cmd[200] ;
90	//sprintf(cmd, ".! tar zxvf DBpp.tgz") ;
91	//gROOT->ProcessLine(cmd) ;
92
93	seed = nrun * 100000 + nevent;
94	sprintf(sseed,"%d",seed);
95
96	if (seed==0) {
97	fprintf(stderr,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");
98	fprintf(stderr,"!!!! WARNING! Seeding variable for MC is 0 !!!!\n");
99	fprintf(stderr,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");
100	} else {
101	fprintf(stdout,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");
102	fprintf(stdout,"!!! MC Seed is %d \n",seed);
103	fprintf(stdout,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");
104	}
105
106	// set the seed environment variable
107	gSystem->Setenv("CONFIG_SEED",sseed);
108	gSystem->Setenv("RUN",srun);
109	gSystem->Setenv("EVENT",sevent);
110	gSystem->Setenv("DC_RUN_TYPE",sprocess);//"kPyGammaJetPHOS");
111	gSystem->Setenv("YEAR",syear);//"14000");
112	gSystem->Setenv("ECMS",secms);//"14000");
113	gSystem->Setenv("PTHARDMIN",sminpthard);//"20");
114	gSystem->Setenv("PTHARDMAX",smaxpthard);//"30");
115	gSystem->Setenv("PI0GAMMAINDET",spi0gammafrag);//"1");
116	gSystem->Setenv("PTGAMMAPI0MIN",sminptgammapi0);//"1");
117	gSystem->Setenv("QUENCHING",squench);
118	gSystem->Setenv("QHAT",sqhat);
119	//gSystem->Setenv("MEDLENGTH",smedlength);
120	gSystem->Setenv("TRIGGER",strig);
121
122	cout<< "SIMRUN:: Run = " << gSystem->Getenv("RUN") << "; Event = " << gSystem->Getenv("EVENT")
123	<< "; Process = " << gSystem->Getenv("DC_RUN_TYPE")
124	<< "; E cms = " << gSystem->Getenv("ECMS")
125	<< "; Trigger configuration = "<< gSystem->Getenv("TRIGGER")
126	<< "; Year = "<<gSystem->Getenv("YEAR") << endl;
127
128	cout<< " minpthard = " << gSystem->Getenv("PTHARDMIN")
129	<< "; maxpthard =" << gSystem->Getenv("PTHARDMAX")
130	<< "; force pi0/gamma frag in detector = " << gSystem->Getenv("PI0GAMMAINDET")
131	<< "; minpt = " << gSystem->Getenv("PTGAMMAPI0MIN")
132	<< "; quenching option ="<<gSystem->Getenv("QUENCHING")
133	<< "; qhat = "<<gSystem->Getenv("QHAT")
134	//<< "; medium length = "<<gSystem->Getenv("MEDIUMLENGTH")
135	<< endl;
136
137	cout<<">>>>> SIMULATION <<<<<"<<endl;
138	gSystem->Exec("aliroot -b -q sim.C > sim.log 2>&1");
139	cout<<">>>>> RECONSTRUCTION <<<<<"<<endl;
140	gSystem->Exec("aliroot -b -q rec.C > rec.log 2>&1");
141	cout<<">>>>> TAG <<<<<"<<endl;
142	gSystem->Exec("aliroot -b -q tag.C > tag.log 2>&1");
143	cout<<">>>>> CHECK ESD <<<<<"<<endl;
144	gSystem->Exec("aliroot -b -q CheckESD.C > check.log 2>&1");
145
146	}