Jpp - the software that should make you happy
 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 
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 
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.0/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
do echo Generating $dir eval D
Definition: JDrawLED.sh:50
set_variable INPUT_FILE
exit
Definition: JPizza.sh:36
then for HISTOGRAM in h0 h1
Definition: JMatrixNZ.sh:69
then JPizza f
Definition: JPizza.sh:46
then echo
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
then $JPP_DIR examples JDetector JSlewing P slewing
Definition: JSlewing.sh:41
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
then awk F
* usage
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
&set_variable OUTPUT_FILE
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:45
then JCalibrateToT a
Definition: JTuneHV.sh:116
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
then echo Launching message logger in new window xterm geometry e JGetMessage H $HOST T $TAG d3
alias put_queue eval echo n
Definition: qlib.csh:19
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
then display $WORKDIR
Definition: plot-Domino.sh:127
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR tmp set_variable JDAQ_TIMESLICE JDAQTimeslice if do_usage *then usage $script< detector file >< inputfile > fi case set_variable WORKDIR
Definition: JLegolas.sh:23
set_variable DETECTOR
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 source JAcoustics sh $DETECTOR_ID CHECK_EXIT_CODE typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
Definition: JCanberra.sh:38
version
Definition: JCalibratePMT.sh:7
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:60
void scale(vector< double > &v, double c)
scale vector content