]>
Commit | Line | Data |
---|---|---|
2c82db0e | 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 | ||
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 | ||
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 | } |