Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JSlewingK40.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 #
4 # \author mdejong
5 #
6 version=1.0
7 script=${0##*/}
8 
9 # ------------------------------------------------------------------------------------------
10 #
11 # Utility script to test JSlewingK40.
12 #
13 # ------------------------------------------------------------------------------------------
14 
15 
16 if [ -z $JPP_DIR ]; then
17  echo "Variable JPP_DIR undefined."
18  exit
19 fi
20 
21 
22 source $JPP_DIR/setenv.sh $JPP_DIR
23 
24 
25 set_variable DEBUG 2
26 set_variable WORKDIR .
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 
37 set_variable TMAX 3.0
38 set_variable DETECTOR $WORKDIR/detector.detx
39 set_variable BACKGROUND_HZ 5e3 500 50 5 0.5
40 set_variable numberOfSlices 10000
41 set_variable INPUT_FILE $WORKDIR/timeslice.root
42 set_variable OUTPUT_FILE $WORKDIR/slewing\[Monte\ Carlo\].root
43 
44 
45 if (( 1 )); then
46 
47  detectorType=1
48  numberOfLines=1
49  numberOfFloors=1
50 
51  JDetector \
52  -D 1 \
53  -@detectorType=$detectorType \
54  -@numberOfLines=$numberOfLines \
55  -@numberOfFloors=$numberOfFloors \
56  -o $DETECTOR -V V2 \
57  -d $DEBUG
58 
59 fi
60 
61 if (( 1 )); then
62 
63  if [ -f "$WORKDIR/timeslice.root" ]; then
64  rm -i $WORKDIR/timeslice.root
65  fi
66 
67  if [ ! -f "$WORKDIR/timeslice.root" ]; then
68 
69  JRandomTimesliceWriter \
70  -a $DETECTOR \
71  -o $INPUT_FILE \
72  -B "$BACKGROUND_HZ" \
73  -n $numberOfSlices \
74  -P %.slewing=1 \
75  -d $DEBUG --!
76 
77  fi
78 fi
79 
80 
81 if (( 1 )); then
82 
83  echo "Monitoring data."
84  echo "Detector: $DETECTOR"
85 
86  JSlewingK40 \
87  -f $WORKDIR/timeslice.root \
88  -a $DETECTOR \
89  -o $OUTPUT_FILE \
90  -t "20 35" \
91  -d $DEBUG
92 
93 fi
94 
95 if (( 1 )); then
96 
97  JScale1D \
98  -f ${OUTPUT_FILE}:h0 \
99  -f ${OUTPUT_FILE}:h1 \
100  -F "1/GetEntries" \
101  -O "nosw2" \
102  -o $WORKDIR/scale.root --!
103 
104 fi
105 
106 if (( 1 )); then
107 
108  JFit \
109  -f ${WORKDIR}/scale.root:h0 \
110  -o ${WORKDIR}/fit.root \
111  -F "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2])) + [3]" \
112  -@"p0 = GetMaximum" \
113  -@"p1 = 0.0" \
114  -@"p2 = 0.5*GetRMS" \
115  -@"p3 = GetMinimum" \
116  -d3
117 
118  JPlot1D \
119  -f ${WORKDIR}/fit.root:h0 \
120  -\> "#Deltat [ns]" \
121  -\^ "number of hits [au]" \
122  -o resolution.gif
123 
124 fi
125 
126 if (( 1 )); then
127 
128  JPlot1D \
129  -f $WORKDIR/scale.root:h1 \
130  -\> "time over threshold [ns]" \
131  -\^ "number of hits [au]" \
132  -o tot.gif
133 
134 fi
135 
136 if (( 1 )); then
137 
138  JFit \
139  -f ${OUTPUT_FILE}:h2 \
140  -o ${WORKDIR}/fit.root \
141  -F "[0]*exp([1]*sqrt(x) + [2]*x) + [3]" \
142  -@"p0 = GetMaximum" \
143  -@"p1 = -0.01" \
144  -@"p2 = -0.05" \
145  -@"p3 = GetMinimum" \
146  -d3
147 
148  JPlot1D \
149  -f ${WORKDIR}/fit.root:h2 \
150  -y "-10 +10" \
151  -\> "TOT [ns]" \
152  -\^ "#Deltat [ns]" \
153  -o slewing.gif
154 
155 fi
D
do echo Generating $dir eval D
Definition: JDrawLED.sh:50
JSYSTEM::set_variable
void set_variable(const std::string &name, const std::string &value)
Set environment variable.
Definition: JSystemToolkit.hh:404
DETECTOR
set_variable DETECTOR
Definition: JLegolas.sh:31
F
then check_input_file $DETECTOR $INPUT_FILE for OPTION in A B C D E F
Definition: JFilter.sh:47
in
esac print_variable DETECTOR INPUT_FILE OUTPUT_FILE CDF for TYPE in
Definition: JSirene.sh:45
exit
exit
Definition: JPizza.sh:36
h1
then for HISTOGRAM in h0 h1
Definition: JMatrixNZ.sh:69
INPUT_FILE
&set_variable INPUT_FILE
Definition: JLegolas.sh:30
a
fi JEventTimesliceWriter a
Definition: JEventProcessor.sh:71
f
then JPizza f
Definition: JPizza.sh:46
echo
then echo
Definition: JDOMDAQDriver.sh:56
d
then print_variable DETECTOR INPUT_FILE INTERMEDIATE_FILE check_input_file $DETECTOR $INPUT_FILE check_output_file $INTERMEDIATE_FILE $OUTPUT_FILE JMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JPath.sh:52
slewing
then $JPP_DIR examples JDetector JSlewing P slewing
Definition: JSlewing.sh:41
usage
* usage
Definition: JLigier-local.sh:33
$INPUT_FILE
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
d3
then echo Launching message logger in new window xterm geometry e JGetMessage H $HOST T $TAG d3
Definition: software/JLogger/JLogger.sh:39
JRECONSTRUCTION::JFit
JFIT::JFit JFit
Definition: JEvt.hh:386
n
alias put_queue eval echo n
Definition: qlib.csh:19
$WORKDIR
then display $WORKDIR
Definition: plot-Domino.sh:142
BACKGROUND_HZ
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR set_variable DEBUG set_variable WORKDIR set_variable PIPE $WORKDIR pipe dat set_variable OUTPUT_FILE $WORKDIR timeslice root set_variable DETECTOR $JPP_DATA km3net_reference detx set_variable TRIGGER $JPP_DATA trigger_parameters_arca txt set_variable PMT_FILE $JPP_DATA PMT_parameters txt set_variable BACKGROUND_HZ
Definition: JTimesliceReprocessor.sh:28
version
version
Definition: JCalibratePMT.sh:7
WORKDIR
then usage $script< detector file >< inputfile > fi case set_variable WORKDIR
Definition: JLegolas.sh:28
P
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:60
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
scale
void scale(vector< double > &v, double c)
scale vector content
Definition: JBeaconSimulator.cc:73
Generated by   doxygen 1.8.5