Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
JCharge.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2#
3# \author mdejong & bjjung
4#
5version=1.0
6script=${0##*/}
7
8# ------------------------------------------------------------------------------------------
9#
10# Utility script to plot charge distribution from JDETECTOR::JPMTDefaultSignalProcessor.
11#
12# ------------------------------------------------------------------------------------------
13
14
15if [ -z $JPP_DIR ]; then
16 echo "Variable JPP_DIR undefined."
17 exit
18fi
19
20
21source $JPP_DIR/setenv.sh $JPP_DIR
22
23
24set_variable DEBUG 3
25set_variable WORKDIR ${TMPDIR:-/tmp}
26set_variable OUTPUT_FILE
27set_variable PMT_FILE "gain=0.9"
28set_variable NPE 1
29set_variable X_RANGE "0 $((3*NPE))"
30set_variable: FORMAT GRAPHICS_FORMAT gif
31set_variable+ BATCH GRAPHICS_BATCH -B
32
33
34if do_usage $*; then
35 usage "$script [output file] [pmt file]"
36fi
37
38case $# in
39 2) set_variable PMT_FILE $2;&
40 1) set_variable OUTPUT_FILE $1;;
41esac
42
43
44$JPP_DIR/examples/JDetector/JCharge \
45 -o $WORKDIR/charge.root \
46 -n 10000000 \
47 -N $NPE \
48 -P $PMT_FILE \
49 -d $DEBUG
50
51for TEST in 1 2; do
52
53 JPlot1D \
54 -f $WORKDIR/charge.root:${TEST}.1 \
55 -f $WORKDIR/charge.root:${TEST}.0 \
56 -x $X_RANGE \
57 -> "Charge [npe]" \
58 -\^ "Probability density [npe^{-1}]" \
59 -S 0.4 \
60 -o charge.$FORMAT $BATCH
61
62 JSum1D \
63 -f $WORKDIR/charge.root:${TEST}.0 \
64 -f $WORKDIR/charge.root:${TEST}.1 \
65 -o $WORKDIR/sum.root -B
66
67 JPlot1D \
68 -f $WORKDIR/sum.root:${TEST}.1 \
69 -f $WORKDIR/sum.root:${TEST}.0 \
70 -x $X_RANGE \
71 -y "1e-3 1.1" -Y \
72 -> "Charge [npe]" \
73 -\^ "Probability" \
74 -O "][" \
75 -S 0.4 \
76 -G XY \
77 -o charge.$FORMAT $BATCH
78done
79
80rm -f $WORKDIR/{charge,sum}.root