Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JMultiplicityK40.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 version=1.0
6 script=${0##*/}
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Utility script to make coincidence rate plots.
11 #
12 # ------------------------------------------------------------------------------------------
13 
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 
21 source $JPP_DIR/setenv.sh $JPP_DIR
22 
23 
24 set_variable DIR $JPP_DIR/examples/JPhysics/
27 
28 if ( do_usage $* ); then
29  usage "$script [working directory]"
30 fi
31 
32 case $# in
33  1) WORKDIR=$1;;
34 esac
35 
36 if [[ ! -f $WORKDIR/k40\[KM3NeT\].root ]]; then
37 
38  $DIR/JMultiplicityK40 \
39  -n 100000000 \
40  -F 4 \
41  -G 2 \
42  -D "0.216 5" \
43  -A 1 \
44  -U \
45  -o $WORKDIR/k40\[KM3NeT\].root
46 fi
47 
48 for (( a = 0.0; $a <= 1.0; a += 0.5 )); do
49 
50  set_variable LABEL `printf '%3.1f' $a`
51 
52  if [[ ! -f $WORKDIR/k40\[${LABEL}\].root ]]; then
53 
54  $DIR/JMultiplicityK40 \
55  -n 100000000 \
56  -F 4 \
57  -G 2 \
58  -D "0.216 5" \
59  -A 3 -a $a \
60  -U \
61  -o $WORKDIR/k40\[${LABEL}\].root
62  fi
63 done
64 
65 INPUT_FILES=(`ls $WORKDIR/k40\[*\].root`)
66 
68 
69 JPlot1D \
70  -f "${INPUT_FILES//\.root/.root:pmt}" \
71  -\> "cos(#theta)" \
72  -\^ "acceptance" \
73  -y "0 1.6" \
74  -L TR \
75  -T "PMT" \
76  -o pmt.gif
77 
78 M=2
79 
80 JPlot1D \
81  -f "${INPUT_FILES//\.root/.root:M.${M}}" \
82  -\> "D [m]" \
83  -\^ "rate [Hz]" \
84  -x "0 5" \
85  -y "1 1e2" -Y \
86  -T "M = ${M}" \
87  -o k40.gif
then usage for INPUT_FILE in $INPUT_FILES[*]
do echo Generating $dir eval D
Definition: JDrawLED.sh:50
do $JPP JMEstimator M
Definition: JMEstimator.sh:37
void set_variable(const std::string &name, const std::string &value)
Set environment variable.
exit
Definition: JPizza.sh:36
do set_array DAQHEADER JPrintDAQHeader f
Definition: JTuneHV.sh:79
do echo Generating $dir eval DIR
Definition: JPlotNPE1D.sh:52
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTORXY_TXT $WORKDIR $DETECTORXY_TXT tail read X Y CHI2 RMS printf optimum n $X $Y $CHI2 $RMS awk v Y
then echo
do set_variable OUTPUT_DIRECTORY $WORKDIR T
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
then awk F
* usage
then fatal Not enough arguments fi set_variable DETECTOR $argv[1] INPUT_FILES
then for DETECTOR in Antares KM3NeT
Definition: JMultiPMT.sh:42
then $JPP_DIR software JCalibrate JCalibrateToT a
Definition: JTuneHV.sh:108
alias put_queue eval echo n
Definition: qlib.csh:19
then display $WORKDIR
Definition: plot-Domino.sh:127
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 source JAcoustics sh $DETECTOR_ID typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
Definition: JCanberra.sh:36
version
Definition: JCalibratePMT.sh:7
then usage $script< detector file >< inputfile > fi case set_variable WORKDIR
Definition: JLegolas.sh:28
source $JPP_DIR setenv csh $JPP_DIR eval JShellParser o a A
esac done
Definition: JAddHDE.sh:21
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62