Fix script bugs from EventDisplay_v4 issue #352

scripts/ANALYSIS.anasum.sh

@@ -87,8 +87,8 @@ sed -e "s|FILELIST|$FLIST|" \
 chmod u+x "$FSCRIPT.sh"
 echo "$FSCRIPT.sh"
 
-SUBC="$(dirname "$0")/helper_scripts/UTILITY.readSubmissionCommand.sh"
-# run locally or on cluster
+SUBC=`$( dirname "$0" )/helper_scripts/UTILITY.readSubmissionCommand.sh`
+# run locally or on cluster (expand shell variables in submission string)
 SUBC=`eval "echo \"$SUBC\""`
 if [[ $SUBC == *"ERROR"* ]]; then
     echo "$SUBC"
@@ -113,5 +113,3 @@ elif [[ $SUBC == *parallel* ]]; then
 elif [[ $SUBC == *simple* ]]; then
 	"$FSCRIPT.sh" |& tee $FSCRIPT.log
 fi
-
-exit

scripts/ANALYSIS.anasum_allcuts.sh

@@ -63,15 +63,15 @@ do
         fi
         echo $CDIR
         mkdir -p "$TMPDIR/${CDIR}"
-        ./ANALYSIS.anasum_parallel_from_runlist.sh ${RUNL} \
+        $(dirname "$0")/ANALYSIS.anasum_parallel_from_runlist.sh ${RUNL} \
             "$TMPDIR/${CDIR}" \
             ${C} \
             ${IGNORETYPE} \
             $VERITAS_EVNDISP_AUX_DIR/ParameterFiles/ANASUM.runparameter \
             $PREDIR $SKIPIFPROCESSED | tee -a ${TMPLOG}
     elif [[ $RUNTYPE == "V2DL3" ]]; then
         mkdir -p "$TMPDIR/v2dl3_${C}"
-         ./ANALYSIS.v2dl3.sh ${RUNL} \
+         $(dirname "$0")/ANALYSIS.v2dl3.sh ${RUNL} \
              "$TMPDIR/v2dl3_${C}" \
              ${C} | tee -a ${TMPLOG}
     else

scripts/ANALYSIS.anasum_combine.sh

@@ -76,6 +76,15 @@ SUBSCRIPT=$( dirname "$0" )"/helper_scripts/ANALYSIS.anasum_combine_sub"
 FSCRIPT="$LOGDIR/anasum_combine-$DATE-RUN$RUN-$(date +%s)"
 echo "Run script written to $FSCRIPT"
 
+# Check that anasum output files exist before combining
+NANASUMFILES=$(find "$DDIR" -maxdepth 1 -name "*.anasum.root" | wc -l)
+if [[ $NANASUMFILES -eq 0 ]]; then
+    echo "Error: no .anasum.root files found in $DDIR"
+    echo "Run ANALYSIS.anasum_parallel_from_runlist.sh and wait for completion before combining."
+    exit 1
+fi
+echo "Found $NANASUMFILES anasum output files in $DDIR"
+
 sed -e "s|RRUNLIST|$RUNLIST|" \
     -e "s|DDDIR|$DDIR|" \
     -e "s|RRUNP|$RUNP|" \

scripts/ANALYSIS.mscw_energy.sh

@@ -12,7 +12,7 @@ if [ "$#" -lt 2 ]; then
 echo "
 MSCW_ENERGY data analysis: submit jobs from a simple run list
 
-ANALYSIS.mscw_energy.sh <runlist> [output directory] [evndisp directory] [preprocessing skip] [Rec ID] [ATM]
+ANALYSIS.mscw_energy.sh <runlist> [output directory] [evndisp directory] [preprocessing skip] [Rec ID]
 
 required parameters:
 
@@ -32,9 +32,6 @@ optional parameters:
     [Rec ID]                reconstruction ID. Default 0
                             (see EVNDISP.reconstruction.runparameter)
 
-    [ATM]                   set atmosphere ID (overwrite the value from the evndisp stage)
-
-
 The analysis type (cleaning method; direction reconstruction) is read from the \$VERITAS_ANALYSIS_TYPE environmental
 variable (e.g., AP_DISP, NN_DISP; here set to: \"$VERITAS_ANALYSIS_TYPE\").
 
@@ -60,7 +57,6 @@ RUNLIST=$1
 [[ "$3" ]] && INPUTDIR=$3 || INPUTDIR="$VERITAS_PREPROCESSED_DATA_DIR/${VERITAS_ANALYSIS_TYPE:0:2}/evndisp"
 [[ "$4" ]] && SKIP=$4 || SKIP=1
 [[ "$5" ]] && ID=$5 || ID=0
-[[ "$6" ]] && FORCEDATMO=$6
 DISPBDT="1"
 
 # Read runlist

scripts/ANALYSIS.v2dl3.sh

@@ -54,7 +54,7 @@ DATE=`date +"%y%m%d"`
 LOGDIR="$VERITAS_USER_LOG_DIR/V2DL3-${DATE}-$(uuidgen)/"
 mkdir -p "$LOGDIR"
 echo -e "Log files will be written to:\n $LOGDIR"
-rm -f ${LOGIDR}/x* 2>/dev/null
+rm -f ${LOGDIR}/x* 2>/dev/null
 
 # split run list into smaller run lists
 sort -u "${RUNLIST}" -o "${LOGDIR}/$(basename "${RUNLIST}")"
@@ -98,11 +98,11 @@ do
             JOBID=$( echo "$JOBID" | grep -oP "Your job [0-9.-:]+" | awk '{ print $3 }' )
         fi
 
-        echo "RUN $AFILE JOBID $JOBID"
-        echo "RUN $AFILE SCRIPT $FSCRIPT.sh"
+        echo "RUN $J JOBID $JOBID"
+        echo "RUN $J SCRIPT $FSCRIPT.sh"
         if [[ $SUBC != */dev/null* ]] ; then
-            echo "RUN $AFILE OLOG $FSCRIPT.sh.o$JOBID"
-            echo "RUN $AFILE ELOG $FSCRIPT.sh.e$JOBID"
+            echo "RUN $J OLOG $FSCRIPT.sh.o$JOBID"
+            echo "RUN $J ELOG $FSCRIPT.sh.e$JOBID"
         fi
     elif [[ $SUBC == *condor* ]]; then
         $(dirname "$0")/helper_scripts/UTILITY.condorSubmission.sh $FSCRIPT.sh $h_vmem $tmpdir_size
@@ -116,7 +116,7 @@ do
         $SUBC $FSCRIPT.sh
     elif [[ $SUBC == *parallel* ]]; then
         echo "$FSCRIPT.sh &> $FSCRIPT.log" >> "$LOGDIR/runscripts.$TIMETAG.dat"
-        echo "RUN $AFILE OLOG $FSCRIPT.log"
+        echo "RUN $J OLOG $FSCRIPT.log"
     elif [[ "$SUBC" == *simple* ]] ; then
 	    "$FSCRIPT.sh" |& tee "$FSCRIPT.log"
 	fi

scripts/IRF.combine_lookup_table_parts.sh

@@ -74,7 +74,7 @@ echo -e "Output files will be written to:\n $ODIR"
 mkdir -p "$ODIR"
 chmod g+w "$ODIR"
 
-if [[ -f $OFILE ]]; then
+if [[ -f "$ODIR/$OFILE" ]]; then
     echo "ERROR: table file $ODIR/$OFILE exists; move it or delete it"
     exit 1
 fi

scripts/IRF.generate_lookup_table_parts.sh

@@ -112,18 +112,18 @@ if [[ $SUBC == *"ERROR"* ]]; then
     exit 1
 fi
 if [[ $SUBC == *qsub* ]]; then
-    JOBID=`$SUBC $FSCRIPT.sh`
+    JOBID=`$SUBC $FSCRIPT`
     echo "JOBID: $JOBID"
 elif [[ $SUBC == *condor* ]]; then
-    $(dirname "$0")/helper_scripts/UTILITY.condorSubmission.sh $FSCRIPT.sh $h_vmem $tmpdir_size
+    $(dirname "$0")/helper_scripts/UTILITY.condorSubmission.sh $FSCRIPT $h_vmem $tmpdir_size
     echo "-------------------------------------------------------------------------------"
     echo "Job submission using HTCondor - run the following script to submit jobs:"
     echo "$EVNDISPSCRIPTS/helper_scripts/submit_scripts_to_htcondor.sh ${LOGDIR} submit"
     echo "-------------------------------------------------------------------------------"
 elif [[ $SUBC == *sbatch* ]]; then
-    $SUBC $FSCRIPT.sh
+    $SUBC $FSCRIPT
 elif [[ $SUBC == *parallel* ]]; then
-    echo "$FSCRIPT.sh &> $FSCRIPT.log" >> "$LOGDIR/runscripts.dat"
+    echo "$FSCRIPT &> $FSCRIPT.log" >> "$LOGDIR/runscripts.dat"
 elif [[ "$SUBC" == *simple* ]]; then
-    "$FSCRIPT.sh" | tee "$FSCRIPT.log"
+    "$FSCRIPT" | tee "$FSCRIPT.log"
 fi

scripts/IRF.mscw_energy_MC.sh

@@ -137,18 +137,18 @@ if [[ $SUBC == *"ERROR"* ]]; then
     exit 1
 fi
 if [[ $SUBC == *qsub* ]]; then
-    JOBID=`$SUBC $FSCRIPT.sh`
+    JOBID=`$SUBC $FSCRIPT`
     echo "JOBID: $JOBID"
 elif [[ $SUBC == *condor* ]]; then
-    $(dirname "$0")/helper_scripts/UTILITY.condorSubmission.sh $FSCRIPT.sh $h_vmem $tmpdir_size
+    $(dirname "$0")/helper_scripts/UTILITY.condorSubmission.sh $FSCRIPT $h_vmem $tmpdir_size
     echo "-------------------------------------------------------------------------------"
     echo "Job submission using HTCondor - run the following script to submit jobs:"
     echo "$EVNDISPSCRIPTS/helper_scripts/submit_scripts_to_htcondor.sh ${LOGDIR} submit"
     echo "-------------------------------------------------------------------------------"
 elif [[ $SUBC == *sbatch* ]]; then
-    $SUBC $FSCRIPT.sh
+    $SUBC $FSCRIPT
 elif [[ $SUBC == *parallel* ]]; then
-    echo "$FSCRIPT.sh &> $FSCRIPT.log" >> "$LOGDIR/runscripts.dat"
+    echo "$FSCRIPT &> $FSCRIPT.log" >> "$LOGDIR/runscripts.dat"
 elif [[ "$SUBC" == *simple* ]]; then
-    "$FSCRIPT.sh" | tee "$FSCRIPT.log"
+    "$FSCRIPT" | tee "$FSCRIPT.log"
 fi

scripts/IRF.optimizeTMVAforGammaHadronSeparation.sh

@@ -80,6 +80,7 @@ fi
 
 if [[ ! -e $VERITAS_EVNDISP_AUX_DIR/EffectiveAreas/${EFFFILE} ]]; then
     echo "Error - effective area file not found ${EFFFILE}"
+    exit 1
 fi
 if [[ ${EPOCH:0:2} == "V4" ]] || [[ ${EPOCH:0:2} == "V5" ]]; then
     RUNPAR="$VERITAS_EVNDISP_AUX_DIR/ParameterFiles/TMVA.BDT.V4.runparameter"
@@ -88,8 +89,8 @@ else
 fi
 #####################################
 # energy bins
-if grep -q "^* ENERGYBINS" "$RUNPAR"; then
-    ENBINS=$( cat "$RUNPAR" | grep "^* ENERGYBINS 1" | sed -e 's/* ENERGYBINS 1//' | sed -e 's/ /\n/g')
+if grep -q "^\* ENERGYBINS" "$RUNPAR"; then
+    ENBINS=$( cat "$RUNPAR" | grep "^\* ENERGYBINS" | sed -e 's/\* ENERGYBINS//' | sed -e 's/ /\n/g')
     declare -a EBINARRAY=( $ENBINS ) #convert to array
     count1=1
     NENE=$((${#EBINARRAY[@]}-$count1)) #get number of bins

scripts/IRF.production.sh

@@ -308,15 +308,15 @@ for VX in $EPOCH; do
        if [[ $IRFTYPE == "ANAXGBANGRES" ]]; then
             MSCWDIR="$VERITAS_IRFPRODUCTION_DIR/$EDVERSION/${ANATYPE}/${SIMTYPE}/${VX}_ATM${ATM}_gamma/MSCW_RECID${RECID}_DISP"
             echo "XGB reconstruction reading from $MSCWDIR"
-            ./IRF.dispXGB.sh "stereo_analysis" "${MSCWDIR}" "${MSCWDIR}"
+            $(dirname "$0")/IRF.dispXGB.sh "stereo_analysis" "${MSCWDIR}" "${MSCWDIR}"
             continue
        fi
        #############################################
        # Classification XGB based on MSCW files
        if [[ $IRFTYPE == "ANAXGBGH" ]]; then
             MSCWDIR="$VERITAS_IRFPRODUCTION_DIR/$EDVERSION/${ANATYPE}/${SIMTYPE}/${VX}_ATM${ATM}_gamma/MSCW_RECID${RECID}_DISP"
             echo "XGB classification reading from $MSCWDIR"
-            ./IRF.dispXGB.sh "classification" "${MSCWDIR}" "${MSCWDIR}"
+            $(dirname "$0")/IRF.dispXGB.sh "classification" "${MSCWDIR}" "${MSCWDIR}"
             continue
        fi
        #############################################
@@ -327,62 +327,58 @@ for VX in $EPOCH; do
            ODIR="$VERITAS_IRFPRODUCTION_DIR/$EDVERSION/${ANATYPE}/${SIMTYPE}/${VX}_ATM${ATM}_gamma/TrainXGBGammaHadron"
            echo "XGB Classification Training"
            echo "${BCKDIR}" "${RUNPAR}" "${ODIR}" "${SIMTYPE}" "${VX}" "${ATM}"
-           ./IRF.trainXGBforGammaHadronSeparationTraining.sh "${BCKDIR}" "${RUNPAR}" "${ODIR}" "${SIMTYPE}" "${VX}" "${ATM}"
+           $(dirname "$0")/IRF.trainXGBforGammaHadronSeparationTraining.sh "${BCKDIR}" "${RUNPAR}" "${ODIR}" "${SIMTYPE}" "${VX}" "${ATM}"
            continue
        fi
        #############################################
        # MVA training
        # train per epoch and atmosphere and for each cut
        # (cut as sizesecondmax cut is applied)
        if [[ $IRFTYPE == "TRAINTMVA" ]] || [[ $IRFTYPE == "OPTIMIZETMVA" ]]; then
-            for VX in $EPOCH; do
-                for ATM in $ATMOS; do
-                    for C in ${CUTTYPES[@]}; do
-                        echo "Training/optimising TMVA for $C cuts, ${VX} ATM${ATM}"
-                        BDTDIR="$VERITAS_IRFPRODUCTION_DIR/$EDVERSION/${ANATYPE}/BDTtraining"
-                        MVADIR="${BDTDIR}/GammaHadronBDTs_${VX:0:2}/${VX}_ATM${ATM}/${C/PointSource-/}/"
-                        # list of background files
-                        TRAINDIR="${BDTDIR}/mscw_${VX:0:2}/"
-                        if [[ $DISPBDT == "1" ]]; then
-                            TRAINDIR="${BDTDIR}/mscw_${VX:0:2}_DISP/"
-                            MVADIR="${BDTDIR}/GammaHadronBDTs_${VX:0:2}_DISP/${VX}_ATM${ATM}/${C/PointSource-/}/"
-                        fi
-                        mkdir -p -v "${MVADIR}"
-                        if [[ $IRFTYPE == "TRAINTMVA" ]]; then
-                            # retrieve size cut
-                            CUTFIL="$VERITAS_EVNDISP_AUX_DIR"/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-${C}-TMVA-Preselection.dat
-                            echo "CUTFILE: $CUTFIL"
-                            SIZECUT=`grep "* sizesecondmax" $CUTFIL | grep ${EPOCH:0:2} | awk '{print $3}' | sort -u`
-                            if [ -z "$SIZECUT" ]
-                            then
-                                echo "No size cut found; skipping cut $C"
-                                continue
-                            fi
-                            echo "Size cut applied: $SIZECUT"
-                            RUNPAR="TMVA.BDT.runparameter"
-                            if [[ ${EPOCH:0:2} == "V4" ]] || [[ ${EPOCH:0:2} == "V5" ]]; then
-                                cp -f "$VERITAS_EVNDISP_AUX_DIR"/ParameterFiles/TMVA.BDT.V4.runparameter "$MVADIR"/"$RUNPAR"
-                            else
-                                cp -f "$VERITAS_EVNDISP_AUX_DIR"/ParameterFiles/"$RUNPAR" "$MVADIR"/"$RUNPAR"
-                            fi
-                            sed -i "s/TMVASIZECUT/${SIZECUT}/" "$MVADIR"/"$RUNPAR"
-                            if [[ $CUTFIL = *"NTel3"* ]]; then
-                                sed -i "s/NImages>1/NImages>2/" "$MVADIR"/"$RUNPAR"
-                            fi
-                            ./IRF.trainTMVAforGammaHadronSeparation.sh \
-                                         "${TRAINDIR}" \
-                                         "$MVADIR"/"$RUNPAR" \
-                                         "${MVADIR}" BDT ${SIMTYPE} ${VX} "${ATM}"
-                         # Cut optimization
-                         elif [[ $IRFTYPE == "OPTIMIZETMVA" ]]; then
-                             echo "OPTIMIZE TMVA $C ${BDTDIR}/BackgroundRates/${VX:0:2}"
-                             ./IRF.optimizeTMVAforGammaHadronSeparation.sh \
-                                 "${BDTDIR}/BackgroundRates/${VX:0:2}" \
-                                 "${C/PointSource-/}" \
-                                 ${SIMTYPE} ${VX} "${ATM}"
-                         fi
-                    done
-                done
+            for C in ${CUTTYPES[@]}; do
+                echo "Training/optimising TMVA for $C cuts, ${VX} ATM${ATM}"
+                BDTDIR="$VERITAS_IRFPRODUCTION_DIR/$EDVERSION/${ANATYPE}/BDTtraining"
+                MVADIR="${BDTDIR}/GammaHadronBDTs_${VX:0:2}/${VX}_ATM${ATM}/${C/PointSource-/}/"
+                # list of background files
+                TRAINDIR="${BDTDIR}/mscw_${VX:0:2}/"
+                if [[ $DISPBDT == "1" ]]; then
+                    TRAINDIR="${BDTDIR}/mscw_${VX:0:2}_DISP/"
+                    MVADIR="${BDTDIR}/GammaHadronBDTs_${VX:0:2}_DISP/${VX}_ATM${ATM}/${C/PointSource-/}/"
+                fi
+                mkdir -p -v "${MVADIR}"
+                if [[ $IRFTYPE == "TRAINTMVA" ]]; then
+                    # retrieve size cut
+                    CUTFIL="$VERITAS_EVNDISP_AUX_DIR"/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-${C}-TMVA-Preselection.dat
+                    echo "CUTFILE: $CUTFIL"
+                    SIZECUT=`grep "* sizesecondmax" $CUTFIL | grep ${VX:0:2} | awk '{print $3}' | sort -u`
+                    if [ -z "$SIZECUT" ]
+                    then
+                        echo "No size cut found; skipping cut $C"
+                        continue
+                    fi
+                    echo "Size cut applied: $SIZECUT"
+                    RUNPAR="TMVA.BDT.runparameter"
+                    if [[ ${VX:0:2} == "V4" ]] || [[ ${VX:0:2} == "V5" ]]; then
+                        cp -f "$VERITAS_EVNDISP_AUX_DIR"/ParameterFiles/TMVA.BDT.V4.runparameter "$MVADIR"/"$RUNPAR"
+                    else
+                        cp -f "$VERITAS_EVNDISP_AUX_DIR"/ParameterFiles/"$RUNPAR" "$MVADIR"/"$RUNPAR"
+                    fi
+                    sed -i "s/TMVASIZECUT/${SIZECUT}/" "$MVADIR"/"$RUNPAR"
+                    if [[ $CUTFIL = *"NTel3"* ]]; then
+                        sed -i "s/NImages>1/NImages>2/" "$MVADIR"/"$RUNPAR"
+                    fi
+                    $(dirname "$0")/IRF.trainTMVAforGammaHadronSeparation.sh \
+                                 "${TRAINDIR}" \
+                                 "$MVADIR"/"$RUNPAR" \
+                                 "${MVADIR}" BDT ${SIMTYPE} ${VX} "${ATM}"
+                 # Cut optimization
+                 elif [[ $IRFTYPE == "OPTIMIZETMVA" ]]; then
+                     echo "OPTIMIZE TMVA $C ${BDTDIR}/BackgroundRates/${VX:0:2}"
+                     $(dirname "$0")/IRF.optimizeTMVAforGammaHadronSeparation.sh \
+                         "${BDTDIR}/BackgroundRates/${VX:0:2}" \
+                         "${C/PointSource-/}" \
+                         ${SIMTYPE} ${VX} "${ATM}"
+                 fi
             done
             continue
        fi

scripts/IRF.selectRunsForGammaHadronSeparationTraining.sh

@@ -57,7 +57,7 @@ if [[ "${RUNPAR##*.}" == "json" ]]; then
         }
     }')
 else
-    ZEBINS=$( cat "$RUNPAR" | grep "^* ZENBINS " | sed -e 's/* ZENBINS//' | sed -e 's/ /\n/g')
+    ZEBINS=$(grep "^\* ZENBINS " "$RUNPAR" | sed -e 's/\* ZENBINS//' | sed -e 's/ /\n/g')
 fi
 echo "Zenith angle definition: $ZEBINS"
 declare -a ZEBINARRAY=( $ZEBINS ) #convert to array

scripts/IRF.trainXGBforAngularReconstructionBinned.sh

@@ -73,7 +73,7 @@ MSCWLIST="$ODIR/xgbFiles.list"
 rm -f ${MSCWLIST}
 touch ${MSCWLIST}
 
-INDIR="$VERITAS_IRFPRODUCTION_DIR/$EDVERSION/${ANALYSIS_TYPE}/$SIMTYPE/${EPOCH}_ATM${ATM}_gamma/MSCW_RECID0_DISP"
+INDIR="$VERITAS_IRFPRODUCTION_DIR/$EDVERSION/${ANALYSIS_TYPE}/$SIMTYPE/${EPOCH}_ATM${ATM}_gamma/MSCW_RECID${RECID}_DISP"
 
 STEREO_PAR="$VERITAS_EVNDISP_AUX_DIR/ParameterFiles/XGB-stereo-parameter.json"
 TRAIN_ANGLES=$(jq -r ".zenith[] | select(.id==\"$ZA\") | .train | join(\" \")" $STEREO_PAR)

scripts/db_scripts/db_combine_dqm_fits.py

@@ -49,7 +49,7 @@ def combine_tables(file_list, output_file):
 
         combined_tables.extend(batch_tables)
 
-    final_combined_table = vstack(combined_tables, join_type="exact")
+    final_combined_table = vstack(combined_tables, join_type="outer")
     logging.info(f"Writing combined table to {output_file}")
     final_combined_table.write(output_file, overwrite=True)
 

scripts/db_scripts/db_run.sh

@@ -74,6 +74,7 @@ get_start_time()
     else
         field_name="data_start_time"
     fi
+    start_time_index=0
     while IFS="|" read -ra a; do
         if [[ ${a[0]} == "run_id" ]]; then
             for (( j=0; j<${#a[@]}; j++ ));
@@ -97,6 +98,7 @@ get_end_time()
     else
         field_name="data_end_time"
     fi
+    end_time_index=0
     while IFS="|" read -ra a; do
         if [[ ${a[0]} == "run_id" ]]; then
             for (( j=0; j<${#a[@]}; j++ ));
@@ -170,6 +172,7 @@ get_source_id()
     else
         OFIL=$(cat $OFIL)
     fi
+    source_index=0
     while IFS="|" read -ra a; do
         if [[ ${a[0]} == "run_id" ]]; then
             for (( j=0; j<${#a[@]}; j++ ));
@@ -289,9 +292,11 @@ read_pointing()
 
 read_run_from_DB runinfo
 read_run_from_DB rundqm
-# don't test and read if tar file exists
-# (implementation of testing missing)
-if [[ ! -e $(getDBTextFileDirectory ${RUN}).tar.gz ]] || [[ ${OVERWRITE} == "1" ]]; then
+# Re-read laser data if tarball is missing OR if sidecar .laserrun file is absent
+RUND="$(getDBTextFileDirectory ${RUN})"
+if [[ ! -e "${RUND}.tar.gz" ]] || \
+   [[ $(get_file_status ${RUN} "${RUND}/${RUN}.laserrun") == "0" ]] || \
+   [[ ${OVERWRITE} == "1" ]]; then
     read_laser_run_and_dqm
     read_laser_calibration
 fi