Skip to content

Commit

Permalink
update
Browse files Browse the repository at this point in the history
  • Loading branch information
@wenqiang-gu
wenqiang-gu committed Oct 19, 2020
1 parent fc3fd38 commit ce7810a
Show file tree
Hide file tree
Showing 13 changed files with 37 additions and 5 deletions.
5 changes: 5 additions & 0 deletions ...t/grid/endscript_DataWCP_reco1.5-reco2.sh → ...t/grid/endscript_DataWCP_reco1.5-reco2.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -140,6 +140,10 @@ echo "++++++++++++++++++++++++++++"
echo "start reco 2"
input_celltree="nuselDATA_WCP.root" # not a celltree any more

# update the eventcount when it exists dummy files from reco 1.5
N_nuseldummy=$(ls nulseldummy*root | wc -l)
eventcount=$((eventcount-N_nuseldummy))

date | tee -a ../wirecell.log
touch WCP_STM.log
touch WCP_analysis.log
Expand Down Expand Up @@ -186,6 +190,7 @@ do

if [ ! -e nue_${input2}.root ]; then
echo "++> $input2 no output." | tee -a WCP_analysis.log
exit 205
else
echo "$n event analysis starts." | tee -a WCP_analysis.log
IFS=$'\n' #dviding symbol
Expand Down
11 changes: 8 additions & 3 deletions ...rid/endscript_OverlayWCP_reco1.5-reco2.sh → ...rid/endscript_OverlayWCP_reco1.5-reco2.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -104,7 +104,7 @@ do
echo $input2 | tee -a ../wirecell.log
mv $imagingfile imaging_$input2.root
#prod-nusel-eval ./input_data_files/ChannelWireGeometry_v2.txt imaging_$input2.root -d1 -p1 2>&1 | tee -a ../wirecell.log
prod-wire-cell-matching-nusel ./input_data_files/ChannelWireGeometry_v2.txt imaging_$input2.root -d1 -z1 2>&1 | tee -a ../wirecell.log
prod-wire-cell-matching-nusel ./input_data_files/ChannelWireGeometry_v2.txt imaging_$input2.root -d1 2>&1 | tee -a ../wirecell.log
### Exception: "No points! Quit!" -- this is a "warning" instead of "error"!
warning=`tail -n 1 ../wirecell.log`
echo "$warning"
Expand Down Expand Up @@ -140,13 +140,17 @@ echo "++++++++++++++++++++++++++++"
echo "start reco 2"
input_celltree="nuselOVERLAY_WCP.root" # not a celltree any more

# update the eventcount when it exists dummy files from reco 1.5
N_nuseldummy=$(ls nulseldummy*root | wc -l)
eventcount=$((eventcount-N_nuseldummy))

date | tee -a ../wirecell.log
touch WCP_STM.log
touch WCP_analysis.log
for ((n=0; n<$(($eventcount)); n++))
do
echo "$n event processing starts." | tee -a ../wirecell.log
wire-cell-prod-stm ./input_data_files/ChannelWireGeometry_v2.txt $input_celltree $n -d0 -o0 -g2 -z1 2>&1 | tee -a ../wirecell.log
wire-cell-prod-stm ./input_data_files/ChannelWireGeometry_v2.txt $input_celltree $n -d0 -o0 -g2 2>&1 | tee -a ../wirecell.log
echo "$n event processing finished." | tee -a ../wirecell.log
echo "$n event cosmic rejection starts." | tee -a ../wirecell.log
for stmfile in stm*.root
Expand Down Expand Up @@ -181,11 +185,12 @@ do
#event=`echo $runinfo | awk -F "_" '{print $3}'`

echo "$n event nue starts." | tee -a WCP_analysis.log
wire-cell-prod-nue ./input_data_files/ChannelWireGeometry_v2.txt $input_celltree $n -d0 -z1 2>&1 | tee -a WCP_analysis.log
wire-cell-prod-nue ./input_data_files/ChannelWireGeometry_v2.txt $input_celltree $n -d0 2>&1 | tee -a WCP_analysis.log
echo "$n event nue finished." | tee -a WCP_analysis.log

if [ ! -e nue_${input2}.root ]; then
echo "++> $input2 no output." | tee -a WCP_analysis.log
exit 205
else
echo "$n event analysis starts." | tee -a WCP_analysis.log
IFS=$'\n' #dviding symbol
Expand Down
0 ...hicl/run_celltreeub_overlay_port_prod.fcl → ...hicl/run_celltreeub_overlay_port_prod.fcl
View file
File renamed without changes.
0 ..._celltreeub_overlay_wiremod_port_prod.fcl → ..._celltreeub_overlay_wiremod_port_prod.fcl
View file
File renamed without changes.
0 ...t/grid/fhicl/run_celltreeub_port_prod.fcl → ...t/grid/fhicl/run_celltreeub_port_prod.fcl
View file
File renamed without changes.
0 ...ript/grid/fhicl/run_slimmed_port_data.fcl → ...ript/grid/fhicl/run_slimmed_port_data.fcl
View file
File renamed without changes.
0 ...t/grid/fhicl/run_slimmed_port_overlay.fcl → ...t/grid/fhicl/run_slimmed_port_overlay.fcl
View file
File renamed without changes.
0 ...script/grid/reco1.5-reco2_DataWCPport.xml → ...script/grid/reco1.5-reco2_DataWCPport.xml
View file
File renamed without changes.
0 ...ipt/grid/reco1.5-reco2_OverlayWCPport.xml → ...ipt/grid/reco1.5-reco2_OverlayWCPport.xml
View file
File renamed without changes.
5 changes: 5 additions & 0 deletions ...teractive/runOne_DataWCP_reco1.5-reco2.sh → ...teractive/runOne_DataWCP_reco1.5-reco2.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -179,6 +179,10 @@ echo "++++++++++++++++++++++++++++"
echo "start reco 2"
input_celltree="nuselDATA_WCP.root" # not a celltree any more

# update the eventcount when it exists dummy files from reco 1.5
N_nuseldummy=$(ls nulseldummy*root | wc -l)
eventcount=$((eventcount-N_nuseldummy))

date | tee -a ../wirecell.log
touch WCP_STM.log
touch WCP_analysis.log
Expand Down Expand Up @@ -225,6 +229,7 @@ do

if [ ! -e nue_${input2}.root ]; then
echo "++> $input2 no output." | tee -a WCP_analysis.log
exit 205
else
echo "$n event analysis starts." | tee -a WCP_analysis.log
IFS=$'\n' #dviding symbol
Expand Down
5 changes: 5 additions & 0 deletions ...active/runOne_OverlayWCP_reco1.5-reco2.sh → ...active/runOne_OverlayWCP_reco1.5-reco2.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -179,6 +179,10 @@ echo "++++++++++++++++++++++++++++"
echo "start reco 2"
input_celltree="nuselOVERLAY_WCP.root" # not a celltree any more

# update the eventcount when it exists dummy files from reco 1.5
N_nuseldummy=$(ls nulseldummy*root | wc -l)
eventcount=$((eventcount-N_nuseldummy))

date | tee -a ../wirecell.log
touch WCP_STM.log
touch WCP_analysis.log
Expand Down Expand Up @@ -225,6 +229,7 @@ do

if [ ! -e nue_${input2}.root ]; then
echo "++> $input2 no output." | tee -a WCP_analysis.log
exit 205
else
echo "$n event analysis starts." | tee -a WCP_analysis.log
IFS=$'\n' #dviding symbol
Expand Down
8 changes: 8 additions & 0 deletions project/reco1.5-reco2/endscript_OverlayWCP_reco1.5-reco2-LEE.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -64,6 +64,9 @@ fi
eventcount=`grep event_count celltreeOVERLAY*.json | awk '{print substr($2, 1, length($2)-1)}'`
# eventcount=1
echo "eventcount:" $eventcount | tee -a wirecell.log
if [ $eventcount -eq 0 ]; then
exit 0
fi

echo "Make WCP work directory" | tee -a wirecell.log
echo `pwd`
Expand Down Expand Up @@ -137,6 +140,10 @@ echo "++++++++++++++++++++++++++++"
echo "start reco 2"
input_celltree="nuselOVERLAY_WCP.root" # not a celltree any more

# update the eventcount when it exists dummy files from reco 1.5
N_nuseldummy=$(ls nulseldummy*root | wc -l)
eventcount=$((eventcount-N_nuseldummy))

date | tee -a ../wirecell.log
touch WCP_STM.log
touch WCP_analysis.log
Expand Down Expand Up @@ -183,6 +190,7 @@ do

if [ ! -e nue_${input2}.root ]; then
echo "++> $input2 no output." | tee -a WCP_analysis.log
exit 205
else
echo "$n event analysis starts." | tee -a WCP_analysis.log
IFS=$'\n' #dviding symbol
Expand Down
8 changes: 6 additions & 2 deletions project/reco1.5-reco2/endscript_OverlayWCP_reco1.5-reco2.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -64,6 +64,9 @@ fi
eventcount=`grep event_count celltreeOVERLAY*.json | awk '{print substr($2, 1, length($2)-1)}'`
# eventcount=1
echo "eventcount:" $eventcount | tee -a wirecell.log
if [ $eventcount -eq 0 ]; then
exit 0
fi

echo "Make WCP work directory" | tee -a wirecell.log
echo `pwd`
Expand Down Expand Up @@ -137,8 +140,9 @@ echo "++++++++++++++++++++++++++++"
echo "start reco 2"
input_celltree="nuselOVERLAY_WCP.root" # not a celltree any more

xx <= nulseldummy*root
eventcount=$eventcount - xx
# update the eventcount when it exists dummy files from reco 1.5
N_nuseldummy=$(ls nulseldummy*root | wc -l)
eventcount=$((eventcount-N_nuseldummy))

date | tee -a ../wirecell.log
touch WCP_STM.log
Expand Down

0 comments on commit ce7810a

Please sign in to comment.