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

//#define VERBOSEARGS\r
{\r
  // set job and simulation variables as : alienoutner needed for seed\r
  // root run.C  --run <x> --event <y> --aliencounter <counter> --process <proc> --minhard <min> --maxhard <max> --minpt <minpt>\r
\r
\r
  char* program = "aliroot";\r
  //  char* program = "alienaliroot";\r
  \r
  int nrun = 0;\r
  int nevent = 0;\r
  int seed = 0;\r
  int naliencounter = 0;\r
  \r
  char sseed[1024];\r
  char srun[1024];\r
  char sevent[1024];\r
  char saliencounter[1024];\r
  char sprocess[1024];\r
  char sminpthard[1024];\r
  char smaxpthard[1024];\r
  char sminptgammapi0[1024];\r
  char squench[1024];\r
  char sqhat[1024];\r
\r
  Int_t npthardbin = 0;\r
  const Int_t nPtHardBins = 3;\r
  const Int_t ptHardLo[nPtHardBins] = {15,50,100};\r
  const Int_t ptHardHi[nPtHardBins] = {50,100,-1};\r
\r
\r
\r
\r
  sprintf(srun,"");\r
  sprintf(sevent,"");\r
  sprintf(saliencounter,"");\r
  sprintf(sprocess,"");\r
  sprintf(sminpthard,"");\r
  sprintf(smaxpthard,"");\r
  sprintf(sminptgammapi0,"");\r
  sprintf(squench,"");\r
  sprintf(sqhat,"");\r
\r
  for (int i=0; i< gApplication->Argc();i++){\r
#ifdef VERBOSEARGS\r
    printf("Arg  %d:  %s\n",i,gApplication->Argv(i));\r
#endif\r
    if (!(strcmp(gApplication->Argv(i),"--run")))\r
      nrun = atoi(gApplication->Argv(i+1));\r
    sprintf(srun,"%d",nrun);\r
    if (!(strcmp(gApplication->Argv(i),"--event")))\r
      nevent = atoi(gApplication->Argv(i+1));\r
    sprintf(sevent,"%d",nevent);\r
    if (!(strcmp(gApplication->Argv(i),"--aliencounter")))\r
      naliencounter = atoi(gApplication->Argv(i+1));\r
    sprintf(saliencounter,"%d",naliencounter);\r
\r
    if (!(strcmp(gApplication->Argv(i),"--process")))\r
      sprintf(sprocess, gApplication->Argv(i+1));\r
\r
    if (!(strcmp(gApplication->Argv(i),"--pthardbin")))\r
      npthardbin = atoi(gApplication->Argv(i+1));\r
\r
    \r
    if (!(strcmp(gApplication->Argv(i),"--minhard")))\r
      sprintf(sminpthard,gApplication->Argv(i+1));\r
    \r
    if (!(strcmp(gApplication->Argv(i),"--maxhard")))\r
      sprintf(smaxpthard,gApplication->Argv(i+1));\r
\r
\r
    if (!(strcmp(gApplication->Argv(i),"--minpt")))\r
      sprintf(sminptgammapi0,gApplication->Argv(i+1));\r
\r
    if (!(strcmp(gApplication->Argv(i),"--quench")))\r
      sprintf(squench,gApplication->Argv(i+1));\r
\r
    if (!(strcmp(gApplication->Argv(i),"--qhat")))\r
      sprintf(sqhat,gApplication->Argv(i+1));\r
\r
  }\r
\r
  if(!(strcmp(sminpthard,""))){\r
    sprintf(sminpthard,"%d",ptHardLo[npthardbin%nPtHardBins]);\r
    Printf(">>> MinPtHard %s %d",sminpthard,npthardbin%nPtHardBins);\r
  }\r
  if(!(strcmp(smaxpthard,""))){\r
    sprintf(smaxpthard,"%d",ptHardHi[npthardbin%nPtHardBins]);\r
    Printf(">>> MaxPtHard %s %d",smaxpthard,npthardbin%nPtHardBins);\r
  }\r
 \r
\r
  \r
  /*\r
  char cmd[200] ; \r
  sprintf(cmd, ".! tar zxvf DBpp.tgz") ; \r
  gROOT->ProcessLine(cmd) ; \r
  sprintf(cmd, ".! tar zxvf TPCCalib_v4-16-Rev-01.tgz") ; \r
  gROOT->ProcessLine(cmd) ; \r
  */\r
  gSystem->Exec("echo $PWD ");\r
  gSystem->Exec("ls -l ");\r
\r
  seed = nrun * 100000 + naliencounter * 1000 + nevent;\r
  sprintf(sseed,"%d",seed);\r
\r
  if (seed==0) {\r
    fprintf(stderr,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");\r
    fprintf(stderr,"!!!!  WARNING! Seeding variable for MC is 0          !!!!\n");\r
    fprintf(stderr,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");\r
  } else {\r
    fprintf(stdout,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");\r
    fprintf(stdout,"!!!  MC Seed is %d \n",seed);\r
    fprintf(stdout,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");\r
  }\r
  \r
  // set the seed environment variable\r
  gSystem->Setenv("CONFIG_SEED",sseed);\r
  gSystem->Setenv("DC_RUN",srun);\r
  gSystem->Setenv("DC_EVENT",sevent);\r
  gSystem->Setenv("DC_RUN_TYPE",sprocess);//"kPyGammaJetPHOS");\r
  gSystem->Setenv("PTHARDMIN",sminpthard);//"20");\r
  gSystem->Setenv("PTHARDMAX",smaxpthard);//"30");\r
//  gSystem->Setenv("DC_RUN_TYPE","kPyJetJet");\r
//  gSystem->Setenv("PTHARDMIN","40");\r
//  gSystem->Setenv("PTHARDMAX","-1");\r
\r
  gSystem->Setenv("PTGAMMAPI0MIN",sminptgammapi0);//"1");\r
  gSystem->Setenv("QUENCHING",squench);\r
  gSystem->Setenv("QHAT",sqhat);\r
//  gSystem->Setenv("QUENCHING","0");\r
//  gSystem->Setenv("QHAT","0");\r
\r
  gSystem->Setenv("ECMS","10000");\r
//  gSystem->Setenv("ECMS","14000");\r
  gSystem->Setenv("ALIMDC_RAWDB1","./mdc1");\r
  gSystem->Setenv("ALIMDC_RAWDB2","./mdc2");\r
  gSystem->Setenv("ALIMDC_TAGDB","./mdc1/tag");\r
  gSystem->Setenv("ALIMDC_RUNDB","./mdc1/meta");\r
  cout<< "SIMRUN:: Run " << gSystem->Getenv("DC_RUN") << " Event " << gSystem->Getenv("DC_EVENT") \r
	  << " Process "   << gSystem->Getenv("DC_RUN_TYPE") \r
	  << " minpthard " << gSystem->Getenv("PTHARDMIN") \r
	  << " maxpthard " << gSystem->Getenv("PTHARDMAX") \r
	  << " minpt "     << gSystem->Getenv("PTGAMMAPI0MIN") \r
	  << endl;\r
  //gSystem->Exec("cp $ROOTSYS/etc/system.rootrc .rootrc");\r
  cout<<">>>>> SIMULATION <<<<<"<<endl;\r
\r
  gSystem->Exec("echo ALICE_ROOT $ALICE_ROOT");\r
  gSystem->Exec("echo ROOTSYS $ROOTSYS");\r
  gSystem->Exec("echo LD_LIB $LD_LIBRARY_PATH");\r
  gSystem->Exec("echo PATH $PATH");\r
  gSystem->Exec("ls -l $ALICE_ROOT/bin/");\r
  gSystem->Exec("ls -l $ALICE_ROOT/bin/*");\r
\r
  gSystem->Exec(Form("%s -b -q \"sim.C(%d)\" > sim.log 2>&1",program,nevent));\r
  \r
\r
  // residual\r
\r
\r
  gSystem->MakeDirectory("residual");\r
  gSystem->Exec("rm residual/*.root");\r
  gSystem->Exec("ln -s ${PWD}/*.root residual/");\r
\r
  // residual_trd\r
  gSystem->MakeDirectory("residual_trd");\r
  gSystem->Exec("rm residual_trd/*.root");\r
  gSystem->Exec("ln -s ${PWD}/*.root residual_trd/");\r
\r
  gSystem->ChangeDirectory("residual");\r
  cout<<">>>>> RECONSTRUCTION <<<<<"<<endl;\r
  Int_t iStatus  = 0;\r
  iStatus = gSystem->Exec(Form("%s -b -q ../rec.C > rec.log 2>&1",program));\r
  Printf("%d",iStatus);\r
  cout<<">>>>> TAG <<<<<"<<endl;\r
  //  iStatus = gSystem->Exec("aliroot -b -q ../tag.C 2>&1 tag.log");\r
  iStatus = gSystem->Exec(Form("%s -b -q ../tag.C > tag.log 2>&1",program));\r
  Printf("%d",iStatus);\r
  cout<<">>>>> CHECK ESD <<<<<"<<endl;\r
  iStatus = gSystem->Exec(Form("%s -b -q ../CheckESD.C > check.log 2>&1",program));\r
  cout<<">>>>> CREATE AOD <<<<<"<<endl;\r
  iStatus = gSystem->Exec(Form("%s -b -q ../runAODFilterJets.C > aod.log 2>&1",program));\r
  Printf("%d",iStatus);\r
\r
  gSystem->Exec(Form("mv Run*tag.root ../"));\r
  gSystem->Exec(Form("mv *RecPoints.root ../")); // for debugging\r
  TString addString("");\r
  // logs\r
  gSystem->Exec(Form("mv rec.log ../rec%s.log",addString.Data()));\r
  gSystem->Exec(Form("mv tag.log ../tag%s.log",addString.Data()));\r
  gSystem->Exec(Form("mv check.log ../check%s.log",addString.Data()));\r
  gSystem->Exec(Form("mv aod.log ../aod%s.log",addString.Data()));\r
  // roots\r
  gSystem->Exec(Form("mv AliESDs.root ../AliESDs%s.root",addString.Data()));\r
  gSystem->Exec(Form("mv AliESDfriends.root ../AliESDfriends%s.root",addString.Data()));\r
  gSystem->Exec(Form("mv AliAOD.root ../AliAOD%s.root",addString.Data()));\r
  gSystem->ChangeDirectory("../");\r
\r
\r
  \r
\r
  // residual_trd\r
\r
  gSystem->ChangeDirectory("residual_trd");\r
  /*\r
  sprintf(cmd, ".! tar zxvf ../DBpp.tgz") ; \r
  gROOT->ProcessLine(cmd) ; \r
  sprintf(cmd, ".! tar zxvf ../TPCCalib_v4-16-Rev-01.tgz") ; \r
  gROOT->ProcessLine(cmd) ; \r
  */\r
  gSystem->Exec("ls -l");\r
  cout<<">>>>> REC <<<<<"<<endl;\r
  iStatus = gSystem->Exec(Form("%s -b -q ../rec2.C > rec.log 2>&1",program));\r
  \r
  Printf("%d",iStatus);\r
  cout<<">>>>> CHECK ESD <<<<<"<<endl;\r
  iStatus = gSystem->Exec(Form("%s -b -q ../CheckESD.C > check.log 2>&1",program));\r
  Printf("%d",iStatus);\r
  cout<<">>>>> CREATE AOD <<<<<"<<endl;\r
  iStatus = gSystem->Exec(Form("%s -b -q ../runAODFilterJets.C > aod.log 2>&1",program));\r
  Printf("%d",iStatus);\r
  TString addString("TRD");\r
  // logs\r
  gSystem->Exec(Form("mv rec.log ../rec%s.log",addString.Data()));\r
  gSystem->Exec(Form("mv check.log ../check%s.log",addString.Data()));\r
  gSystem->Exec(Form("mv aod.log ../aod%s.log",addString.Data()));\r
  // roots\r
  gSystem->Exec(Form("mv AliESDs.root ../AliESDs%s.root",addString.Data()));\r
  gSystem->Exec(Form("mv AliESDfriends.root ../AliESDfriends%s.root",addString.Data()));\r
  gSystem->Exec(Form("mv AliAOD.root ../AliAOD%s.root",addString.Data()));\r
  gSystem->ChangeDirectory("../");\r
  gSystem->Exec("ls -l");\r
}\r
Commit	Line	Data
afea418a	1	//#define VERBOSEARGS\r
	2	{\r
	3	// set job and simulation variables as : alienoutner needed for seed\r
	4	// root run.C --run <x> --event <y> --aliencounter <counter> --process <proc> --minhard <min> --maxhard <max> --minpt <minpt>\r
	5	\r
	6	\r
	7	char* program = "aliroot";\r
	8	// char* program = "alienaliroot";\r
	9	\r
	10	int nrun = 0;\r
	11	int nevent = 0;\r
	12	int seed = 0;\r
	13	int naliencounter = 0;\r
	14	\r
	15	char sseed[1024];\r
	16	char srun[1024];\r
	17	char sevent[1024];\r
	18	char saliencounter[1024];\r
	19	char sprocess[1024];\r
	20	char sminpthard[1024];\r
	21	char smaxpthard[1024];\r
	22	char sminptgammapi0[1024];\r
	23	char squench[1024];\r
	24	char sqhat[1024];\r
	25	\r
	26	Int_t npthardbin = 0;\r
	27	const Int_t nPtHardBins = 3;\r
	28	const Int_t ptHardLo[nPtHardBins] = {15,50,100};\r
	29	const Int_t ptHardHi[nPtHardBins] = {50,100,-1};\r
	30	\r
	31	\r
	32	\r
	33	\r
	34	sprintf(srun,"");\r
	35	sprintf(sevent,"");\r
	36	sprintf(saliencounter,"");\r
	37	sprintf(sprocess,"");\r
	38	sprintf(sminpthard,"");\r
	39	sprintf(smaxpthard,"");\r
	40	sprintf(sminptgammapi0,"");\r
	41	sprintf(squench,"");\r
	42	sprintf(sqhat,"");\r
	43	\r
	44	for (int i=0; i< gApplication->Argc();i++){\r
	45	#ifdef VERBOSEARGS\r
	46	printf("Arg %d: %s\n",i,gApplication->Argv(i));\r
	47	#endif\r
	48	if (!(strcmp(gApplication->Argv(i),"--run")))\r
	49	nrun = atoi(gApplication->Argv(i+1));\r
	50	sprintf(srun,"%d",nrun);\r
	51	if (!(strcmp(gApplication->Argv(i),"--event")))\r
	52	nevent = atoi(gApplication->Argv(i+1));\r
	53	sprintf(sevent,"%d",nevent);\r
	54	if (!(strcmp(gApplication->Argv(i),"--aliencounter")))\r
	55	naliencounter = atoi(gApplication->Argv(i+1));\r
	56	sprintf(saliencounter,"%d",naliencounter);\r
	57	\r
	58	if (!(strcmp(gApplication->Argv(i),"--process")))\r
	59	sprintf(sprocess, gApplication->Argv(i+1));\r
	60	\r
	61	if (!(strcmp(gApplication->Argv(i),"--pthardbin")))\r
	62	npthardbin = atoi(gApplication->Argv(i+1));\r
	63	\r
	64	\r
65	if (!(strcmp(gApplication->Argv(i),"--minhard")))\r
66	sprintf(sminpthard,gApplication->Argv(i+1));\r
67	\r
68	if (!(strcmp(gApplication->Argv(i),"--maxhard")))\r
69	sprintf(smaxpthard,gApplication->Argv(i+1));\r
70	\r
71	\r
72	if (!(strcmp(gApplication->Argv(i),"--minpt")))\r
73	sprintf(sminptgammapi0,gApplication->Argv(i+1));\r
74	\r
75	if (!(strcmp(gApplication->Argv(i),"--quench")))\r
76	sprintf(squench,gApplication->Argv(i+1));\r
77	\r
78	if (!(strcmp(gApplication->Argv(i),"--qhat")))\r
79	sprintf(sqhat,gApplication->Argv(i+1));\r
80	\r
81	}\r
82	\r
83	if(!(strcmp(sminpthard,""))){\r
84	sprintf(sminpthard,"%d",ptHardLo[npthardbin%nPtHardBins]);\r
85	Printf(">>> MinPtHard %s %d",sminpthard,npthardbin%nPtHardBins);\r
86	}\r
87	if(!(strcmp(smaxpthard,""))){\r
88	sprintf(smaxpthard,"%d",ptHardHi[npthardbin%nPtHardBins]);\r
89	Printf(">>> MaxPtHard %s %d",smaxpthard,npthardbin%nPtHardBins);\r
90	}\r
91	\r
92	\r
93	\r
94	/*\r
95	char cmd[200] ; \r
96	sprintf(cmd, ".! tar zxvf DBpp.tgz") ; \r
97	gROOT->ProcessLine(cmd) ; \r
98	sprintf(cmd, ".! tar zxvf TPCCalib_v4-16-Rev-01.tgz") ; \r
99	gROOT->ProcessLine(cmd) ; \r
100	*/\r
101	gSystem->Exec("echo $PWD ");\r
102	gSystem->Exec("ls -l ");\r
103	\r
104	seed = nrun * 100000 + naliencounter * 1000 + nevent;\r
105	sprintf(sseed,"%d",seed);\r
106	\r
107	if (seed==0) {\r
108	fprintf(stderr,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");\r
109	fprintf(stderr,"!!!! WARNING! Seeding variable for MC is 0 !!!!\n");\r
110	fprintf(stderr,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");\r
111	} else {\r
112	fprintf(stdout,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");\r
113	fprintf(stdout,"!!! MC Seed is %d \n",seed);\r
114	fprintf(stdout,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");\r
115	}\r
116	\r
117	// set the seed environment variable\r
118	gSystem->Setenv("CONFIG_SEED",sseed);\r
119	gSystem->Setenv("DC_RUN",srun);\r
120	gSystem->Setenv("DC_EVENT",sevent);\r
121	gSystem->Setenv("DC_RUN_TYPE",sprocess);//"kPyGammaJetPHOS");\r
122	gSystem->Setenv("PTHARDMIN",sminpthard);//"20");\r
123	gSystem->Setenv("PTHARDMAX",smaxpthard);//"30");\r
124	// gSystem->Setenv("DC_RUN_TYPE","kPyJetJet");\r
125	// gSystem->Setenv("PTHARDMIN","40");\r
126	// gSystem->Setenv("PTHARDMAX","-1");\r
127	\r
128	gSystem->Setenv("PTGAMMAPI0MIN",sminptgammapi0);//"1");\r
129	gSystem->Setenv("QUENCHING",squench);\r
130	gSystem->Setenv("QHAT",sqhat);\r
131	// gSystem->Setenv("QUENCHING","0");\r
132	// gSystem->Setenv("QHAT","0");\r
133	\r
134	gSystem->Setenv("ECMS","10000");\r
135	// gSystem->Setenv("ECMS","14000");\r
136	gSystem->Setenv("ALIMDC_RAWDB1","./mdc1");\r
137	gSystem->Setenv("ALIMDC_RAWDB2","./mdc2");\r
138	gSystem->Setenv("ALIMDC_TAGDB","./mdc1/tag");\r
139	gSystem->Setenv("ALIMDC_RUNDB","./mdc1/meta");\r
140	cout<< "SIMRUN:: Run " << gSystem->Getenv("DC_RUN") << " Event " << gSystem->Getenv("DC_EVENT") \r
141	<< " Process " << gSystem->Getenv("DC_RUN_TYPE") \r
142	<< " minpthard " << gSystem->Getenv("PTHARDMIN") \r
143	<< " maxpthard " << gSystem->Getenv("PTHARDMAX") \r
144	<< " minpt " << gSystem->Getenv("PTGAMMAPI0MIN") \r
145	<< endl;\r
146	//gSystem->Exec("cp $ROOTSYS/etc/system.rootrc .rootrc");\r
147	cout<<">>>>> SIMULATION <<<<<"<<endl;\r
148	\r
149	gSystem->Exec("echo ALICE_ROOT $ALICE_ROOT");\r
150	gSystem->Exec("echo ROOTSYS $ROOTSYS");\r
151	gSystem->Exec("echo LD_LIB $LD_LIBRARY_PATH");\r
152	gSystem->Exec("echo PATH $PATH");\r
153	gSystem->Exec("ls -l $ALICE_ROOT/bin/");\r
154	gSystem->Exec("ls -l $ALICE_ROOT/bin/*");\r
155	\r
156	gSystem->Exec(Form("%s -b -q \"sim.C(%d)\" > sim.log 2>&1",program,nevent));\r
157	\r
158	\r
159	// residual\r
160	\r
161	\r
162	gSystem->MakeDirectory("residual");\r
163	gSystem->Exec("rm residual/*.root");\r
164	gSystem->Exec("ln -s ${PWD}/*.root residual/");\r
165	\r
166	// residual_trd\r
167	gSystem->MakeDirectory("residual_trd");\r
168	gSystem->Exec("rm residual_trd/*.root");\r
169	gSystem->Exec("ln -s ${PWD}/*.root residual_trd/");\r
170	\r
171	gSystem->ChangeDirectory("residual");\r
172	cout<<">>>>> RECONSTRUCTION <<<<<"<<endl;\r
173	Int_t iStatus = 0;\r
174	iStatus = gSystem->Exec(Form("%s -b -q ../rec.C > rec.log 2>&1",program));\r
175	Printf("%d",iStatus);\r
176	cout<<">>>>> TAG <<<<<"<<endl;\r
177	// iStatus = gSystem->Exec("aliroot -b -q ../tag.C 2>&1 tag.log");\r
178	iStatus = gSystem->Exec(Form("%s -b -q ../tag.C > tag.log 2>&1",program));\r
179	Printf("%d",iStatus);\r
180	cout<<">>>>> CHECK ESD <<<<<"<<endl;\r
181	iStatus = gSystem->Exec(Form("%s -b -q ../CheckESD.C > check.log 2>&1",program));\r
182	cout<<">>>>> CREATE AOD <<<<<"<<endl;\r
183	iStatus = gSystem->Exec(Form("%s -b -q ../runAODFilterJets.C > aod.log 2>&1",program));\r
184	Printf("%d",iStatus);\r
185	\r
186	gSystem->Exec(Form("mv Run*tag.root ../"));\r
187	gSystem->Exec(Form("mv *RecPoints.root ../")); // for debugging\r
188	TString addString("");\r
189	// logs\r
190	gSystem->Exec(Form("mv rec.log ../rec%s.log",addString.Data()));\r
191	gSystem->Exec(Form("mv tag.log ../tag%s.log",addString.Data()));\r
192	gSystem->Exec(Form("mv check.log ../check%s.log",addString.Data()));\r
193	gSystem->Exec(Form("mv aod.log ../aod%s.log",addString.Data()));\r
194	// roots\r
195	gSystem->Exec(Form("mv AliESDs.root ../AliESDs%s.root",addString.Data()));\r
196	gSystem->Exec(Form("mv AliESDfriends.root ../AliESDfriends%s.root",addString.Data()));\r
197	gSystem->Exec(Form("mv AliAOD.root ../AliAOD%s.root",addString.Data()));\r
198	gSystem->ChangeDirectory("../");\r
199	\r
200	\r
201	\r
202	\r
203	// residual_trd\r
204	\r
205	gSystem->ChangeDirectory("residual_trd");\r
206	/*\r
207	sprintf(cmd, ".! tar zxvf ../DBpp.tgz") ; \r
208	gROOT->ProcessLine(cmd) ; \r
209	sprintf(cmd, ".! tar zxvf ../TPCCalib_v4-16-Rev-01.tgz") ; \r
210	gROOT->ProcessLine(cmd) ; \r
211	*/\r
212	gSystem->Exec("ls -l");\r
213	cout<<">>>>> REC <<<<<"<<endl;\r
214	iStatus = gSystem->Exec(Form("%s -b -q ../rec2.C > rec.log 2>&1",program));\r
215	\r
216	Printf("%d",iStatus);\r
217	cout<<">>>>> CHECK ESD <<<<<"<<endl;\r
218	iStatus = gSystem->Exec(Form("%s -b -q ../CheckESD.C > check.log 2>&1",program));\r
219	Printf("%d",iStatus);\r
220	cout<<">>>>> CREATE AOD <<<<<"<<endl;\r
221	iStatus = gSystem->Exec(Form("%s -b -q ../runAODFilterJets.C > aod.log 2>&1",program));\r
222	Printf("%d",iStatus);\r
223	TString addString("TRD");\r
224	// logs\r
225	gSystem->Exec(Form("mv rec.log ../rec%s.log",addString.Data()));\r
226	gSystem->Exec(Form("mv check.log ../check%s.log",addString.Data()));\r
227	gSystem->Exec(Form("mv aod.log ../aod%s.log",addString.Data()));\r
228	// roots\r
229	gSystem->Exec(Form("mv AliESDs.root ../AliESDs%s.root",addString.Data()));\r
230	gSystem->Exec(Form("mv AliESDfriends.root ../AliESDfriends%s.root",addString.Data()));\r
231	gSystem->Exec(Form("mv AliAOD.root ../AliAOD%s.root",addString.Data()));\r
232	gSystem->ChangeDirectory("../");\r
233	gSystem->Exec("ls -l");\r
234	}\r