8 # ------------------------------------------------------------------------------------------
10 # Utility script to determine muon energy correction base on Monte Carlo simulations.
12 # ------------------------------------------------------------------------------------------
14 if [ -z $JPP_DIR ];
then
15 echo "Variable JPP_DIR undefined."
19 source $JPP_DIR/setenv.sh $JPP_DIR
28 set_variable: FORMAT GRAPHICS_FORMAT gif
32 usage "$script <working directory> <option>"\
33 "\nPossible options: run fit interpolate"
37 fatal
"Wrong number of arguments."
44 fatal
"Invalid option $OPTION."
57 for ((
i = 1;
$i <= $NUMBER_OF_RUNS; ++
i ));
do
64 if [[ $OPTION ==
"fit" || $OPTION ==
"interpolate" ]];
then
66 if ( ! reuse_file
$WORKDIR/energy_correction.root );
then
71 fatal
"No input file in directory $WORKDIR"
79 if [[ $OPTION ==
"fit" ]];
then
82 -
f "$INPUT_FILES[*]" \
84 -
F "[0] + [1]*x + [2]*x*x + [3]*x*x*x" \
93 elif [[ $OPTION ==
"interpolate" ]];
then
96 -
f "$INPUT_FILES[*]" \
105 -\>
"E_{fit} [GeV]" \
106 -\<
"E_{true} [GeV]" \
112 -
o ExE.$FORMAT $BATCH
120 -\>
"log(E_{fit}/GeV)" \
121 -\^
"log(E_{true}/GeV)" \
then usage $script< detector >< run >< outputfile > fi case set_variable OPTION
then usage $script< detector identifier >< startdate\"YYYY-MM-DDHH:MM:SS\"><finaldate\"YYYY-MM-DDHH:MM:SS\"><QA/QCfile> fi case set_variable QAQC_TXT $argv[4]
then usage $script< input file >[option[primary[working directory]]] nWhere option can be E
o $QUALITY_ROOT d $DEBUG!CHECK_EXIT_CODE JPlot1D f
usr share Modules init zsh export TMPDIR
fatal Invalid number of arguments esac let YMIN
then fatal Wrong number of arguments fi JConvertDetectorFormat a o
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
do set_variable OUTPUT_DIRECTORY $WORKDIR T
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR
do JSirene sh $DETECTOR $JPP_DATA genhen km3net_wpd_V2_0 evt gz $WORKDIR sirene_
skip elif((BINFRAC< 1.0))
do sput $STRING $FLOOR sget MODULE INPUT_FILES
then usage $script< input file >[option[primary[working directory]]] nWhere option can be N
fi set_variable FILENAME $WORKDIR
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
then set_variable MODULE getModule a $DETECTOR L $STRING $FLOOR JEditDetector a $DETECTOR M $MODULE add $X o $DETECTOR else echo No update of detector $DETECTOR
do set_variable MODULE getModule a $WORKDIR detector_a datx L $STRING JEditDetector a $WORKDIR detector_a datx M $MODULE setz o $WORKDIR detector_a datx JEditDetector a $WORKDIR detector_b datx M $MODULE setz o $WORKDIR detector_b datx done echo Output stored at $WORKDIR detector_a datx and $WORKDIR tripod_a txt JDrawDetector2D a $WORKDIR detector_a datx a $WORKDIR detector_b datx L BL o detector $FORMAT $BATCH JDrawDetector2D T $WORKDIR tripod_a txt T $WORKDIR tripod_b txt L BL o tripod $FORMAT $BATCH JCompareDetector a $WORKDIR detector_a datx b $WORKDIR detector_b datx o $WORKDIR abc root &dev null for KEY in X Y Z
then eval ls $DIR $ID $TYPE _
do set_variable DETECTOR_TXT $WORKDIR detector
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
*fatal Wrong option $OPTION
#define DEBUG(A)
Message macros.