make_xml_mcc1.0.sh
Go to the documentation of this file.
1 #! /bin/bash
2 #----------------------------------------------------------------------
3 #
4 # Name: make_xml_mcc1.0.sh
5 #
6 # Purpose: Make xml files for mcc 1.0. This script loops over all
7 # generator-level fcl files in the source area of the currently
8 # setup version of lbnecode (that is, under
9 # $LBNECODE_DIR/source/fcl/lbne35t/gen), and makes a corresponding xml
10 # project file in the local directory.
11 #
12 # Usage:
13 #
14 # make_xml_mcc1.0.sh [-h|--help] [-r <release>] [-u|--user <user>] [--local <dir|tar>] [--nev <n>] [--nevjob <n>] [--nevgjob <n>]
15 #
16 # Options:
17 #
18 # -h|--help - Print help.
19 # -r <release> - Use the specified larsoft/lbnecode release.
20 # -u|--user <user> - Use users/<user> as working and output directories
21 # (default is to use lbnepro).
22 # --local <dir|tar> - Specify larsoft local directory or tarball (xml
23 # tag <local>...</local>).
24 # --nev <n> - Specify number of events for all samples. Otherwise
25 # use hardwired defaults.
26 # --nevjob <n> - Specify the default number of events per job.
27 # --nevgjob <n> - Specify the maximum number of events per gen/g4 job.
28 #
29 #----------------------------------------------------------------------
30 
31 # Parse arguments.
32 
33 rel=v03_06_00_01
34 userdir=lbnepro
35 userbase=$userdir
36 nevarg=0
37 nevjob=0
38 nevgjobarg=0
39 local=''
40 
41 while [ $# -gt 0 ]; do
42  case "$1" in
43 
44  # User directory.
45 
46  -h|--help )
47  echo "Usage: make_xml_mcc5.0.sh [-h|--help] [-r <release>] [-u|--user <user>] [--local <dir|tar>] [--nev <n>] [--nevjob <n>] [--nevgjob <n>]"
48  exit
49  ;;
50 
51  # Release.
52 
53  -r )
54  if [ $# -gt 1 ]; then
55  rel=$2
56  shift
57  fi
58  ;;
59 
60  # User.
61 
62  -u|--user )
63  if [ $# -gt 1 ]; then
64  userdir=users/$2
65  userbase=$2
66  shift
67  fi
68  ;;
69 
70  # Local release.
71 
72  --local )
73  if [ $# -gt 1 ]; then
74  local=$2
75  shift
76  fi
77  ;;
78 
79  # Total number of events.
80 
81  --nev )
82  if [ $# -gt 1 ]; then
83  nevarg=$2
84  shift
85  fi
86  ;;
87 
88  # Number of events per job.
89 
90  --nevjob )
91  if [ $# -gt 1 ]; then
92  nevjob=$2
93  shift
94  fi
95  ;;
96 
97  # Number of events per gen/g4 job.
98 
99  --nevgjob )
100  if [ $# -gt 1 ]; then
101  nevgjobarg=$2
102  shift
103  fi
104  ;;
105 
106  esac
107  shift
108 done
109 
110 # Get qualifier.
111 
112 qual=e5
113 ver=`echo $rel | cut -c2-3`
114 if [ $ver -gt 2 ]; then
115  qual=e6
116 fi
117 
118 # Delete existing xml files.
119 
120 rm -f *.xml
121 
122 find $LBNECODE_DIR/source/fcl/lbne35t/gen -name \*.fcl | while read fcl
123 do
124  if ! echo $fcl | grep -q common; then
125  newprj=`basename $fcl .fcl`
126  newxml=${newprj}.xml
127  filt=1
128  generator=SingleGen
129  if echo $newprj | grep -q cosmics; then
130  generator=CRY
131  fi
132  if echo $newprj | grep -q textfile; then
133  generator=TextFileGen
134  fi
135 
136  # Make xml file.
137 
138  echo "Making ${newprj}.xml"
139 
140  # Generator
141 
142  genfcl=`basename $fcl`
143 
144  # G4
145 
146  g4fcl=standard_g4_lbne35t.fcl
147 
148  # Detsim (optical + tpc).
149 
150  detsimfcl=standard_detsim_lbne35t.fcl
151  if echo $newprj | grep -q milliblock; then
152  detsimfcl=standard_detsim_lbne35t_milliblock.fcl
153  fi
154  # Reco 2D
155 
156 # reco2dfcl=standard_reco_uboone_2D.fcl
157 
158  # Reco 3D
159 
160 # reco3dfcl=standard_reco_uboone_3D.fcl
161 
162  # Reco
163  recofcl=standard_reco_lbne35t.fcl
164  if echo $newprj | grep -q milliblock; then
165  recofcl=standard_reco_lbne35t_milliblock.fcl
166  fi
167 
168  # Merge/Analysis
169 
170  mergefcl=standard_ana_lbne35t.fcl
171 
172  # Set number of events per job.
173  if [ $nevjob -eq 0 ]; then
174  if [ $newprj = prodcosmics_lbne35t_milliblock ]; then
175  nevjob=100
176  elif [ $newprj = prodcosmics_lbne35t_onewindow ]; then
177  nevjob=100
178  elif [ $newprj = prodtextfile_lbne35t ]; then
179  nevjob=1000
180  else
181  nevjob=100
182  fi
183  fi
184  # Set number of gen/g4 events per job.
185 
186  nevgjob=$nevgjobarg
187  if [ $nevgjob -eq 0 ]; then
188  if echo $newprj | grep -q dirt; then
189  if echo $newprj | grep -q cosmic; then
190  nevgjob=200
191  else
192  nevgjob=2000
193  fi
194  else
195  nevgjob=nevjob
196  fi
197  fi
198 
199  # Set number of events.
200 
201  nev=$nevarg
202  if [ $nev -eq 0 ]; then
203  if [ $newprj = prodcosmics_lbne35t_milliblock ]; then
204  nev=10000
205  elif [ $newprj = prodcosmics_lbne35t_onewindow ]; then
206  nev=10000
207  elif [ $newprj = prodtextfile_lbne35t ]; then
208  nev=1000000
209  else
210  nev=10000
211  fi
212  fi
213  nev=$(( $nev * $filt ))
214 
215  # Calculate the number of worker jobs.
216 
217  njob1=$(( $nev / $nevgjob )) # Pre-filter (gen, g4)
218  njob2=$(( $nev / $nevjob / $filt )) # Post-filter (detsim and later)
219  if [ $njob1 -lt $njob2 ]; then
220  njob1=$njob2
221  fi
222 
223  cat <<EOF > $newxml
224 <?xml version="1.0"?>
225 
226 <!-- Production Project -->
227 
228 <!DOCTYPE project [
229 <!ENTITY release "$rel">
230 <!ENTITY file_type "mc">
231 <!ENTITY run_type "physics">
232 <!ENTITY name "$newprj">
233 <!ENTITY tag "mcc1.0">
234 ]>
235 
236 <project name="&name;">
237 
238  <!-- Group -->
239  <group>lbne</group>
240 
241  <!-- Project size -->
242  <numevents>$nev</numevents>
243 
244  <!-- Operating System -->
245  <os>SL6</os>
246 
247  <!-- Batch resources -->
248  <resource>DEDICATED,OPPORTUNISTIC</resource>
249 
250  <!-- Larsoft information -->
251  <larsoft>
252  <tag>&release;</tag>
253  <qual>${qual}:prof</qual>
254 EOF
255  echo "local=$local"
256  if [ x$local != x ]; then
257  echo " <local>${local}</local>" >> $newxml
258  fi
259  cat <<EOF >> $newxml
260  </larsoft>
261 
262  <!-- lbne35t metadata parameters -->
263 
264  <parameter name ="MCName">${newprj}</parameter>
265  <parameter name ="MCDetectorType">35t</parameter>
266  <parameter name ="MCGenerators">${generator}</parameter>
267 
268  <!-- Project stages -->
269 
270  <stage name="gen">
271  <fcl>$genfcl</fcl>
272 EOF
273  if echo $newprj | grep -q textfile; then
274  echo " <inputmode>textfile</inputmode>" >> $newxml
275  echo " <inputlist>/lbne/app/users/tjyang/larsoft_dev/job/AntiMuonCutEvents.txt</inputlist>" >> $newxml
276  fi
277  cat <<EOF >> $newxml
278  <outdir>/lbne/data2/${userdir}/&release;/gen/&name;</outdir>
279  <workdir>/lbne/app/users/${userbase}/&release;/gen/&name;</workdir>
280  <numjobs>$njob1</numjobs>
281  <datatier>generated</datatier>
282  <defname>&name;_&tag;_gen</defname>
283  </stage>
284 
285  <stage name="g4">
286  <fcl>$g4fcl</fcl>
287  <outdir>/lbne/data2/${userdir}/&release;/g4/&name;</outdir>
288  <workdir>/lbne/app/users/${userbase}/&release;/g4/&name;</workdir>
289  <numjobs>$njob1</numjobs>
290  <datatier>simulated</datatier>
291  <defname>&name;_&tag;_g4</defname>
292  </stage>
293 
294 EOF
295  if [ x$detsimfcl != x ]; then
296  cat <<EOF >> $newxml
297  <stage name="detsim">
298  <fcl>$detsimfcl</fcl>
299  <outdir>/lbne/data2/${userdir}/&release;/detsim/&name;</outdir>
300  <workdir>/lbne/app/users/${userbase}/&release;/detsim/&name;</workdir>
301  <numjobs>$njob2</numjobs>
302  <datatier>detector-simulated</datatier>
303  <defname>&name;_&tag;_detsim</defname>
304  </stage>
305 
306 EOF
307  fi
308  if [ x$optsimfcl != x ]; then
309  cat <<EOF >> $newxml
310  <stage name="optsim">
311  <fcl>$optsimfcl</fcl>
312  <outdir>/lbne/data2/${userdir}/&release;/optsim/&name;</outdir>
313  <workdir>/lbne/app/users/${userbase}/&release;/optsim/&name;</workdir>
314  <numjobs>$njob2</numjobs>
315  <datatier>optical-simulated</datatier>
316  <defname>&name;_&tag;_optsim</defname>
317  </stage>
318 
319 EOF
320  fi
321  if [ x$tpcsimfcl != x ]; then
322  cat <<EOF >> $newxml
323  <stage name="tpcsim">
324  <fcl>$tpcsimfcl</fcl>
325  <outdir>/lbne/data2/${userdir}/&release;/tpcsim/&name;</outdir>
326  <workdir>/lbne/app/users/${userbase}/&release;/tpcsim/&name;</workdir>
327  <numjobs>$njob2</numjobs>
328  <datatier>tpc-simulated</datatier>
329  <defname>&name;_&tag;_tpcsim</defname>
330  </stage>
331 
332 EOF
333  fi
334  cat <<EOF >> $newxml
335  <stage name="reco">
336  <fcl>$recofcl</fcl>
337  <outdir>/lbne/data2/${userdir}/&release;/reco/&name;</outdir>
338  <workdir>/lbne/app/users/${userbase}/&release;/reco/&name;</workdir>
339  <numjobs>$njob2</numjobs>
340  <datatier>full-reconstructed</datatier>
341  <defname>&name;_&tag;_reco</defname>
342  </stage>
343 
344  <stage name="mergeana">
345  <fcl>$mergefcl</fcl>
346  <outdir>/lbne/data2/${userdir}/&release;/mergeana/&name;</outdir>
347  <workdir>/lbne/app/users/${userbase}/&release;/mergeana/&name;</workdir>
348  <numjobs>$njob2</numjobs>
349  <targetsize>8000000000</targetsize>
350  <datatier>full-reconstructed</datatier>
351  <defname>&name;_&tag;</defname>
352  </stage>
353 
354  <!-- file type -->
355  <filetype>&file_type;</filetype>
356 
357  <!-- run type -->
358  <runtype>&run_type;</runtype>
359 
360 </project>
361 EOF
362 
363  fi
364 
365 done
Generated by   doxygen 1.8.11