Jpp - the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JCalibrateK40-PMT.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 JCalibrateK40/JFitK40.
12 #
13 # ------------------------------------------------------------------------------------------
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 source $JPP_DIR/setenv.sh $JPP_DIR
21 
22 set_variable DIR $JPP_DIR/examples/JCalibrate
24 set_variable WORKDIR /tmp/
25 
26 if do_usage $*; then
27  usage "$script [working directory]"
28 fi
29 
30 case $# in
31  1) set_variable WORKDIR $1;;
32 esac
33 
35 set_variable BACKGROUND_HZ 5e3 500
36 set_variable numberOfSlices 6000
37 set_variable DAQ_TIMESLICE JDAQTimesliceL1
38 
40 set_variable numberOfLines 1
41 set_variable numberOfFloors 1
42 
43 JDetector \
44  -D 1 \
45  -@ "detectorType = $detectorType" \
46  -@ "numberOfLines = $numberOfLines" \
47  -@ "numberOfFloors = $numberOfFloors" \
48  -o $DETECTOR \
49  -d $DEBUG
50 
51 set_variable MODULE `getModule -a $DETECTOR -L "1 1"`
52 
53 if ( ! reuse_file $WORKDIR/timeslice.root ); then
54 
55  JRandomTimesliceWriter \
56  -a $DETECTOR \
57  -o $WORKDIR/timeslice.root \
58  -B "$BACKGROUND_HZ" \
59  -n $numberOfSlices \
60  -d $DEBUG --!
61 fi
62 
63 if ( ! reuse_file $WORKDIR/trigger_processor.root ); then
64 
65  JTriggerProcessor \
66  -a $DETECTOR \
67  -f $WORKDIR/timeslice.root \
68  -o $WORKDIR/trigger_processor.root \
69  -C JDAQTimesliceL0 \
70  -@"ctMin=-1.0" \
71  -@"TMaxLocal_ns=25" \
72  -@"writeL1=1" \
73  -d $DEBUG
74 fi
75 
76 if (( 1 )); then
77 
78  JCalibrateK40 \
79  -f $WORKDIR/trigger_processor.root \
80  -C $DAQ_TIMESLICE \
81  -a $DETECTOR \
82  -o $WORKDIR/monitor.root \
83  -d $DEBUG
84 
85  JMergeCalibrateK40 \
86  -f $WORKDIR/monitor.root \
87  -o $WORKDIR/merge.root \
88  -d $DEBUG
89 
90  JFitK40 \
91  -a $DETECTOR \
92  -f $WORKDIR/merge.root \
93  -o $WORKDIR/fit.root \
94  -w \
95  -d $DEBUG
96 
97  for (( PMT = 0; PMT < 31; PMT += 1 )); do
98 
99  rm -f $WORKDIR/trigger_processor-\[${PMT}\].root
100 
101  $DIR/JRemovePMT \
102  -f $WORKDIR/trigger_processor.root \
103  -o $WORKDIR/trigger_processor-\[${PMT}\].root \
104  -P "-1 $PMT" \
105  -d $DEBUG
106 
107  JCalibrateK40 \
108  -f $WORKDIR/trigger_processor-\[${PMT}\].root \
109  -C $DAQ_TIMESLICE \
110  -a $DETECTOR \
111  -o $WORKDIR/monitor-\[${PMT}\].root \
112  -d $DEBUG
113 
114  JMergeCalibrateK40 \
115  -f $WORKDIR/monitor-\[${PMT}\].root \
116  -o $WORKDIR/merge-\[${PMT}\].root \
117  -d $DEBUG
118 
119  JFitK40 \
120  -a $DETECTOR \
121  -f $WORKDIR/merge-\[${PMT}\].root \
122  -o $WORKDIR/fit-\[${PMT}\].root \
123  -w \
124  -d $DEBUG
125  done
126 fi
127 
128 INPUT_FILES=(`ls -U $WORKDIR/fit.root $WORKDIR/fit-\[*\].root`)
129 
130 JVariance1D \
131  -f "${INPUT_FILES//\.root/.root:${MODULE}.1QE}" \
132  -x "101 -0.1 +0.1" \
133  -F first \
134  -o $WORKDIR/variance.root
135 
136 JScale1D \
137  -f $WORKDIR/variance.root:h0 \
138  -o $WORKDIR/scale.root \
139  -F "1.0/GetEntries" \
140  -O "nosw2"
141 
142 JPlot1D \
143  -f $WORKDIR/scale.root:h0 \
144  -y "1.0e-3 1.0" -Y \
145  -\> "#DeltaQE" \
146  -\^ "number of PMTs" \
147  -s 1010 \
148  -T "" -o QE-summary.gif
149 
150 
151 JPlot1D \
152  -w 1000x500 \
153  -f $WORKDIR/fit.root:${MODULE}.1QE \
154  -f $WORKDIR/fit-\[0\].root:${MODULE}.1QE \
155  -f $WORKDIR/fit-\[10\].root:${MODULE}.1QE \
156  -f $WORKDIR/fit-\[20\].root:${MODULE}.1QE \
157  -f $WORKDIR/fit-\[30\].root:${MODULE}.1QE \
158  -y "0.7 1.5" \
159  -\> "PMT" \
160  -\^ "QE" \
161  -L TR \
162  -T "" \
163  -o QE.gif
164 
data_type w[N+1][M+1]
Definition: JPolint.hh:741
do echo Generating $dir eval D
Definition: JDrawLED.sh:50
JFitK40_t< TF2 > JFitK40
Type definition for backward compatibility.
Definition: JFitK40.hh:806
static const JPBS_t PMT(3, 4, 2, 3)
PBS of photo-multiplier tube (PMT)
exit
Definition: JPizza.sh:36
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 The file $DIR KM3NeT_00000001_00000000 root already please rename or remove it first
then JPizza f
Definition: JPizza.sh:46
then echo
static const double C
Physics constants.
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
do set_variable OUTPUT_DIRECTORY $WORKDIR T
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
then awk F
* usage
set_array INPUT_FILES
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
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
static void monitor()
Definition: DataQueue.cc:32
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
const JModule & getModule(const JDetector &detector, const JModuleLocation &location)
find module with a given string and floor number
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:60
esac done
Definition: JAddHDE.sh:21
void scale(vector< double > &v, double c)
scale vector content