Jpp  19.1.0
the software that should make you happy
JRandomTimesliceWriterRunByRun.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 script=${0##*/}
6 
7 # ------------------------------------------------------------------------------------------
8 #
9 # Example script to run JRandomTimesliceWriter.
10 #
11 # ------------------------------------------------------------------------------------------
12 
13 if [ -z $JPP_DIR ]; then
14  echo "Variable JPP_DIR undefined."
15  exit
16 fi
17 
18 source $JPP_DIR/setenv.sh $JPP_DIR
19 
20 set_variable: DEBUG TIMESLICE_DEBUG 1
21 set_variable: WORKDIR TIMESLICE_WORKDIR ./
22 set_variable: BACKGROUND_HZ TIMESLICE_BACKGROUND_HZ `getK40Rates`
23 set_variable: SEED TIMESLICE_SEED 0
24 set_variable: NUMBER_OF_SLICES TIMESLICE_SLICES 10
25 set_variable: RECYCLING TIMESLICE_RECYCLING "0 0"
26 set_variable: FACTOR TIMESLICE_FACTOR 1.0
27 set_variable: SAMPLER TIMESLICE_SAMPLER "1 1000"
28 
29 if do_usage $*; then
30  usage "$script <detector file> <output file> <DAQ file> <PMT parameters file>"\
31  "\nAuxiliary script to produce random time slice data."
32 fi
33 
34 if (( $# == 4 )); then
35  set_variable DETECTOR $1
36  set_variable OUTPUT_FILE $2
37  set_variable DAQ_FILE $3
38  set_variable PMT_FILE $4
39 else
40  fatal "Wrong number of arguments."
41 fi
42 
43 JPrintTree -f $DAQ_FILE -@ "type = JDAQSummaryslice" -k number_of_entries | read N
44 
45 set_variable RUNBYRUN "file=$DAQ_FILE; sampler=${SAMPLER%% *} $(($N / 2))"
46 
47 JRandomTimesliceWriter \
48  -a ${DETECTOR} \
49  -o ${OUTPUT_FILE} \
50  -n ${NUMBER_OF_SLICES} \
51  -P ${PMT_FILE} \
52  -P "QE=${FACTOR}" \
53  -S ${SEED} \
54  -B "$BACKGROUND_HZ" \
55  -N "$RECYCLING" \
56  -r "${RUNBYRUN}" \
57  -d ${DEBUG} \
58  --!