Jpp  18.1.0
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Functions | Variables
JEvtReweightGSeaGenParameterScan.sh File Reference

Go to the source code of this file.

Functions

then if (is_CCNikhef)
 
skip elif ((BINFRAC< 1.0))
 
fi for ((j=1;j<=${#PARAMETER_FILES};++j))
 
then fatal Number of input parameters (=${N_PARS}) is incompatible with given formula
 
fi unset PAREQ for ((k=0;k<${#PARAMETERS};++k))
 
do if (($PRIMARY==$NUTYPE))
 
fi done done if ((${#PARINPUTS}==0))
 

Variables

 version
 
then echo Variable JPP_DIR
undefined exit fi source
$JPP_DIR setenv sh $JPP_DIR
source $JPP_DIR software
scripts qlib sh set_variable 
IFS0
 
then usage $script< detector
file >< inputfile >
< particlePDGcode
< reweightingformula
< parameterfile >
< oscillationprobabilitytable >
n Note that if more than one
input file is given all
optional parameters need to be
specified fi case set_variable 
WORKDIR
 
&set_variable OSCPARAMETERS
 
&set_variable OSCPROBTABLE
 
set_array PARAMETER_FILES
 
set_array FORMULAE
 
set_array NUTYPES
 
set_array INPUT_FILES
 
set_variable DETX_FILE
 
 while [[-f ${argv[$i]}]]
 
then echo Submitting
reweighting and histogram
comparison jobs to nikhef stbc
batch 
queue
 
then fatal Please specify a
valid directory for local
scratch files other 
than
 
then echo printf s n
 
then rm
 
do if [["$(JPrintHeader -f ${INPUT_FILE} -k simul --)"!=" 0"]]
 
then set_array INPUT_MC_FILES
 
then $JPP_DIR examples JAAnet
JSortEvt 
f
 
then echo
 
then echo Variable JPP_DIR undefined
 
 exit
 
 fi
 
then submit_job
 
do set_variable MC_MERGE_FILE
 
then rm i
 
then continue
 
do set_variable NUTYPE
 
do set_array PAREQS
 
then set_array PARINPUTS
 
fi set_variable REWEIGHTED_FILE
 

Function Documentation

then if ( is_CCNikhef  )
skip elif ( (BINFRAC< 1.0)  )
fi for ( (j=1;j<=${#PARAMETER_FILES};++j)  )
then fatal Number of input parameters ( )

Definition at line 366 of file JEvtReweightGSeaGenParameterScan.sh.

366 {N_PARS}) is incompatible with given formula: ${FORMULA}"
is
Definition: JDAQCHSM.chsm:167
&set_variable FORMULA
fi unset PAREQ for ( (k=0;k<${#PARAMETERS};++k )
do if ( ($PRIMARY==$NUTYPE )
fi done done if ( (${#PARINPUTS}==0)  )

Variable Documentation

version
Initial value:
=1.0
script=${0##*/}
if [ -z $JPP_DIR ]

Definition at line 5 of file JEvtReweightGSeaGenParameterScan.sh.

then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR source $JPP_DIR software scripts qlib sh set_variable IFS0
Initial value:
{IFS}
set_variable: DEBUG REWEIGHTING_DEBUG 2
set_variable+ BATCH_PROCESSING REWEIGHTING_BATCH_PROCESSING 1
set_variable+ BATCH_LOGGING REWEIGHTING_BATCH_LOGGING 1
set_variable: QUEUE_SIZE REWEIGHTING_QUEUE_SIZE 3
set_variable: QUEUE_NAME REWEIGHTING_QUEUE_NAME ${TMPDIR:-/tmp}/queue
set_variable: JOB_SIZE REWEIGHTING_JOB_SIZE 1
set_variable: TEST_TYPE1D REWEIGHTING_TEST_TYPE1D 0
set_variable: TEST_TYPE1D_PARAMETERS REWEIGHTING_TEST_TYPE1D_PARAMETERS 0.05
set_variable: TEST_TYPE2D REWEIGHTING_TEST_TYPE2D 2
set_variable: TEST_TYPE2D_PARAMETERS REWEIGHTING_TEST_TYPE2D_PARAMETERS 0.05
set_variable: SCRATCHDIR REWEIGHTING_SCRATCHDIR ${TMPDIR:-/tmp}
set_variable: N_TRIGGERED_HITS_RANGE REWEIGHTING_N_TRIGGERED_HITS_RANGE "20 $((2**31))"
set_variable: COSTH_RANGE REWEIGHTING_COSTH_RANGE "-1.0 1.0"
set_variable: ENERGY_RANGE REWEIGHTING_ENERGY_RANGE "0.0 $((2.0**31))"
set_variable: RECOTYPE REWEIGHTING_RECOTYPE "track"
if do_usage $*
usr share Modules init zsh export TMPDIR
&set_variable OSCPARAMETERS
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:21
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR
Definition: JLegolas.sh:20
then echo Submitting reweighting and histogram comparison jobs to nikhef stbc batch queue
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62

Definition at line 24 of file JEvtReweightGSeaGenParameterScan.sh.

then usage $script<detector file><inputfile><particlePDGcode<reweightingformula<parameterfile><oscillationprobabilitytable> n Note that if more than one input file is given all optional parameters need to be specified fi case set_variable WORKDIR

Definition at line 57 of file JEvtReweightGSeaGenParameterScan.sh.

& set_variable OSCPARAMETERS

Definition at line 58 of file JEvtReweightGSeaGenParameterScan.sh.

& set_variable OSCPROBTABLE

Definition at line 59 of file JEvtReweightGSeaGenParameterScan.sh.

set_array PARAMETER_FILES

Definition at line 60 of file JEvtReweightGSeaGenParameterScan.sh.

set_array FORMULAE

Definition at line 61 of file JEvtReweightGSeaGenParameterScan.sh.

do set_array NUTYPES
Initial value:
{NUTYPES[@]} ${argv[$i]}
set_array FORMULAE
${FORMULAE[@]} ${argv[$((i+1))]}
set_array PARAMETER_FILES
${PARAMETER_FILES[@]} ${argv[$((i+2))]}
(( i+=3 ))
set_variable WORKDIR ${argv[-1]}
esac
typeset -A NULABELS
set_array NULABELS
-16 "anutau"
-14 "anumu"
-12 "anue"
12 "nue"
14 "numu"
16 "nutau"
if [[ ${BATCH_PROCESSING} == 1 ]]
&set_variable OSCPROBTABLE
&set_variable OSCPARAMETERS
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:21
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR
Definition: JLegolas.sh:20
set_array PARAMETER_FILES
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
esac done
Definition: JAddHDE.sh:21
do alias $i

Definition at line 62 of file JEvtReweightGSeaGenParameterScan.sh.

do set_array INPUT_FILES
Initial value:
{INPUT_FILES[@]} ${argv[$i]}
(( ++i ))
while (( i < $# - 2 ))
do sput $STRING $FLOOR sget MODULE INPUT_FILES
Definition: JCanberra.sh:77
esac done
Definition: JAddHDE.sh:21
do alias $i

Definition at line 63 of file JEvtReweightGSeaGenParameterScan.sh.

*set_variable DETX_FILE
Initial value:
{argv[1]}
i=2

Definition at line 64 of file JEvtReweightGSeaGenParameterScan.sh.

while[[-f ${argv[$i]}]]

Definition at line 69 of file JEvtReweightGSeaGenParameterScan.sh.

then fatal No batch processing option available yet for CC Lyon else echo Submitting reweighting and histogram comparison jobs to local batch queue

Definition at line 111 of file JEvtReweightGSeaGenParameterScan.sh.

then fatal Please specify a valid directory for local scratch files other than
Initial value:
{TMPDIR:-/tmp}\""
elif ( is_CCLyon )
usr share Modules init zsh export TMPDIR
skip elif((BINFRAC< 1.0))

Definition at line 112 of file JEvtReweightGSeaGenParameterScan.sh.

then echo printf s n
Initial value:
echo "----------------------------"
echo "Remove the above regular files? (y/n)"
read REMOVE_MC_MERGE_FILES
if [[ ${REMOVE_MC_MERGE_FILES} == "y" || ${REMOVE_MC_MERGE_FILES} == "Y" ||
${REMOVE_MC_MERGE_FILES} == "yes" || ${REMOVE_MC_MERGE_FILES} == "Yes" || ${REMOVE_MC_MERGE_FILES} == "YES" ]]
std::istream & read(std::istream &in, JTestSummary &summary, const char delimiter= ' ')
Read test summary.
fi set_array MC_MERGE_FILES(find ${WORKDIR}-maxdepth 1-name"*.merge.root") if(($
then echo
Definition: JQAQC.sh:90

Definition at line 137 of file JEvtReweightGSeaGenParameterScan.sh.

then rm
Initial value:
echo "Dividing into DAQ- and MC-files..."
set_variable INPUT_FILE
do sput $STRING $FLOOR sget MODULE INPUT_FILES
Definition: JCanberra.sh:77
fi set_array MC_MERGE_FILES(find ${WORKDIR}-maxdepth 1-name"*.merge.root") if(($
then echo
Definition: JQAQC.sh:90
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY JAcoustics sh $DETECTOR_ID source JAcousticsToolkit sh CHECK_EXIT_CODE typeset A EMITTERS get_tripods $WORKDIR tripod txt EMITTERS get_transmitters $WORKDIR transmitter txt EMITTERS for EMITTER in
Definition: JCanberra.sh:46

Definition at line 147 of file JEvtReweightGSeaGenParameterScan.sh.

do unset PARINPUTS if[[${PARAMETERSETS:0:1}==\#]]

Definition at line 158 of file JEvtReweightGSeaGenParameterScan.sh.

then set_array INPUT_MC_FILES
Initial value:
else # DAQ-file
set_array INPUT_DAQ_FILES ${INPUT_DAQ_FILES[@]} ${INPUT_FILE}
if (( ${#INPUT_MC_FILES} <= 0 ))
then usage $script[< detector identifier >< run range >]< QA/QCfile > nExample script to produce data quality plots nWhen a detector identifier and run range are data are downloaded from the database nand subsequently stored in the given QA QC file
Definition: JDataQuality.sh:19
Q DAQ
Definition: JDataQuality.sh:59
set_variable INPUT_FILE
then set_array INPUT_MC_FILES
esac done
Definition: JAddHDE.sh:21

Definition at line 159 of file JEvtReweightGSeaGenParameterScan.sh.

then $JPP_DIR examples JAAnet JSortEvt f

Definition at line 175 of file JEvtReweightGSeaGenParameterScan.sh.

then echo
Initial value:
{NAME}\$ ${TEST_TYPE1D} ${TEST_TYPE1D_PARAMETERS}" >> ${STEERING_FILE}
set_array HISTOGRAMS1D ${HISTOGRAMS1D[@]} ${NAME}
else
echo "${NAME}\$ ${TEST_TYPE2D} ${TEST_TYPE2D_PARAMETERS}" >> ${STEERING_FILE}
set_array HISTOGRAMS2D ${HISTOGRAMS2D[@]} ${NAME}
fi
fi
done
# Process MC-files with standard weights
for (( i=1; i<=${#MC_MERGE_FILES}; ++i)); do
set_variable MC_MERGE_FILE ${MC_MERGE_FILES[$i]}
set_variable POSTFIT_FILE ${MC_MERGE_FILE:r}.postfit.root
set_variable COMPARISON_FILE ${POSTFIT_FILE:r}.comparison.root
echo "Submitting histogram comparison job for ${MC_MERGE_FILE:t}"
configure_job \
${SCRATCHDIR}/${COMPARISON_FILE:t:r}.job\%.sh ${JOB_SIZE}
cat >> ${EXEC_FILE} <<EOF
${START_JOB+#!/bin/zsh}
${START_JOB+"if [ -z \${JPP_DIR} ]
do set_variable MC_MERGE_FILE

Definition at line 253 of file JEvtReweightGSeaGenParameterScan.sh.

then echo Variable JPP_DIR undefined

Definition at line 279 of file JEvtReweightGSeaGenParameterScan.sh.

exit

Definition at line 279 of file JEvtReweightGSeaGenParameterScan.sh.

fi

Definition at line 281 of file JEvtReweightGSeaGenParameterScan.sh.

then submit_job
Initial value:
{EXEC_FILE} ${BATCH_PROCESSING:+${QUEUE_NAME} ${SCRATCHDIR} ${JOBIDS[1]}}
set_variable PARDELIM ';'
set_variable PARSET_LIST ${WORKDIR}/parametersets.txt
paste -d"$PARDELIM" ${PARAMETER_FILES[@]} > ${PARSET_LIST}
for (( i=1; i<=${#MC_MERGE_FILES}; ++i))
reset_job_config()
Definition: qlib.sh:148
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:21
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR
Definition: JLegolas.sh:20
set_array PARAMETER_FILES
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
fi set_array MC_MERGE_FILES(find ${WORKDIR}-maxdepth 1-name"*.merge.root") if(($
esac done
Definition: JAddHDE.sh:21

Definition at line 316 of file JEvtReweightGSeaGenParameterScan.sh.

do set_variable MC_MERGE_FILE
Initial value:
echo "Submitting reweighting job for ${MC_MERGE_FILE:t}"
if [[ -f ${MC_MERGE_FILE:r}.reweighting.comparison.meta.txt ]]
o $QUALITY_ROOT d $DEBUG!CHECK_EXIT_CODE JPlot1D f
Definition: JDataQuality.sh:76
data_type r[M+1]
Definition: JPolint.hh:779
fi set_array MC_MERGE_FILES(find ${WORKDIR}-maxdepth 1-name"*.merge.root") if(($
then echo
Definition: JQAQC.sh:90
do set_variable MC_MERGE_FILE
do alias $i

Definition at line 332 of file JEvtReweightGSeaGenParameterScan.sh.

then rm i
Initial value:
{MC_MERGE_FILE:r}.reweighting.comparison.meta.txt
set_array PRIMARIES
$(JPrintHeader -f ${MC_MERGE_FILE} -k primary -- | cut -d'=' -f2)
set_array PRIMARIES
${PRIMARIES[@]}
$(JPrintHeader -f ${MC_MERGE_FILE} -k flux | cut -d' ' -f3)
while IFS="$PARDELIM" read -r PARAMETERSETS
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
std::istream & read(std::istream &in, JTestSummary &summary, const char delimiter= ' ')
Read test summary.
o $QUALITY_ROOT d $DEBUG!CHECK_EXIT_CODE JPlot1D f
Definition: JDataQuality.sh:76
data_type r[M+1]
Definition: JPolint.hh:779
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
double f3(const double x, const double y, const double z)
3D function.
Definition: JPolynome3D.cc:23
do set_variable MC_MERGE_FILE

Definition at line 337 of file JEvtReweightGSeaGenParameterScan.sh.

then continue

Definition at line 351 of file JEvtReweightGSeaGenParameterScan.sh.

do set_variable NUTYPE
Initial value:
{NUTYPES[$j]}
set_array PARAMETERS $(echo $PARAMETERSETS | cut -d"$PARDELIM" -f$j)
set_variable N_FPAR $(echo -n $FORMULA | grep -oP "(^|[^x])\[[0-9]*\]" | cut -d'[' -f2- | sort -u | wc -l)
if (( ${#PARAMETERS} != N_FPAR ))
esac set_variable N_FPAR(echo-n ${FORMULA}|grep-oP"(^|[^x])\[[0-9]*\]"|cut-d'['-f2-|sort-u|wc-l) set_variable N_PARS $(($(wc-w< $
o $QUALITY_ROOT d $DEBUG!CHECK_EXIT_CODE JPlot1D f
Definition: JDataQuality.sh:76
const int n
Definition: JPolint.hh:697
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:21
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable TRIPOD $argv[2] set_variable TX $argv[3] set_variable TY $argv[4] if[[!-f $DETECTOR]]
Definition: JFootprint.sh:28
&set_variable FORMULA
do grep(chi2\[.\]\|gradient\)" $INPUT_FILE > $WORKDIR/stage.txt awk 'BEGIN
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
double u[N+1]
Definition: JPolint.hh:776
then echo
Definition: JQAQC.sh:90

Definition at line 357 of file JEvtReweightGSeaGenParameterScan.sh.

do set_array PAREQS
Initial value:
{PAREQS[@]} "p${k}=${PARAMETERS[$((k+1))]};"
for PRIMARY in ${PRIMARIES[@]}
do set_array PAREQS
&set_variable PRIMARY
Definition: JMuonPostfit.sh:46
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY JAcoustics sh $DETECTOR_ID source JAcousticsToolkit sh CHECK_EXIT_CODE typeset A EMITTERS get_tripods $WORKDIR tripod txt EMITTERS get_transmitters $WORKDIR transmitter txt EMITTERS for EMITTER in
Definition: JCanberra.sh:46
esac done
Definition: JAddHDE.sh:21

Definition at line 373 of file JEvtReweightGSeaGenParameterScan.sh.

then set_array PARINPUTS
Initial value:
{PARINPUTS[@]} "${NULABELS[${NUTYPE}]}=\"${FORMULA}; ${PAREQS[@]}\""
break
then set_array PARINPUTS

Definition at line 381 of file JEvtReweightGSeaGenParameterScan.sh.

fi set_variable REWEIGHTED_FILE

Definition at line 389 of file JEvtReweightGSeaGenParameterScan.sh.