Jpp  15.0.4
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JCalibrateToT.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 JCalibrateToT/JFitToT.
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 /tmp/
27 set_variable: FORMAT GRAPHICS_FORMAT gif
28 set_variable+ BATCH GRAPHICS_BATCH -B
29 
30 if do_usage $*; then
31  usage "$script [working directory]"
32 fi
33 
34 case $# in
35  1) set_variable WORKDIR $1;;
36 esac
37 
38 set_variable DETECTOR $WORKDIR/detector.detx
39 set_variable BACKGROUND_HZ 5e3 500 50 5 0.5
40 set_variable numberOfSlices 1000
41 
42 if (( 1 )); then
43 
44  set_variable detectorType 1
45  set_variable numberOfLines 1
46  set_variable numberOfFloors 1
47 
48  JDetector \
49  -D 1 \
50  -@ "detectorType = $detectorType" \
51  -@ "numberOfLines = $numberOfLines" \
52  -@ "numberOfFloors = $numberOfFloors" \
53  -o $DETECTOR -V V2 \
54  -d $DEBUG
55 
56 fi
57 
58 if (( 1 )); then
59 
60  debug "Generating random background."
61 
62  if ( ! reuse_file $WORKDIR/timeslice.root ); then
63 
64  JRandomTimesliceWriter \
65  -a $DETECTOR \
66  -o $WORKDIR/timeslice.root \
67  -B "$BACKGROUND_HZ" \
68  -n $numberOfSlices \
69  -d $DEBUG --!
70 
71  fi
72 fi
73 
74 if (( 1 )); then
75 
76  debug "Processing data."
77 
78  if ( ! reuse_file $WORKDIR/trigger_processor.root ); then
79 
80  JTriggerProcessor \
81  -a $DETECTOR \
82  -f $WORKDIR/timeslice.root \
83  -o $WORKDIR/trigger_processor.root \
84  -@"ctMin=-1.0" \
85  -@"TMaxLocal_ns=50" \
86  -@"writeL1=1" \
87  -d $DEBUG
88 
89  fi
90 fi
91 
92 if (( 0 )); then
93 
94  debug "Monitoring data."
95 
96  JCalibrateToT \
97  -f $WORKDIR/trigger_processor.root \
98  -o $WORKDIR/monitor.root \
99  -C JDAQTimesliceL1 \
100  -a $DETECTOR \
101  -d $DEBUG
102 
103 fi
104 
105 if (( 0 )); then
106 
107  debug "Fitting data."
108  debug "Detector: $DETECTOR"
109 
110  JFitToT \
111  -f ${WORKDIR}/monitor.root \
112  -o ${WORKDIR}/fit.root \
113  -a $DETECTOR \
114  -w \
115  -d $DEBUG
116 
117 fi
118 
119 if (( 0 )); then
120 
121  set_variable MODULE 101
122  set_variable PMT 0
123 
124  JPlot1D \
125  -f ${WORKDIR}/fit.root:${MODULE}.${PMT}.1ToT \
126  -x "0 40" \
127  -y "1e-3 0.2" -Y \
128  -\> "time over threshold [ns]" \
129  -\^ "probability" \
130  -T "" \
131  -o tot.${MODULE}.${PMT}.$FORMAT $BATCH
132 
133 fi
134 
135 if (( 1 )); then
136 
137  set_variable MODULE 101
138  set_variable PMT 0
139 
140  JPlot1D \
141  -f ${WORKDIR}/fit.root:${MODULE}.1gain\$ \
142  -y "0.0 2.0" \
143  -\> "PMT" \
144  -\^ "gain" \
145  -T "" \
146  -o tot.${MODULE}.gain.$FORMAT $BATCH
147 
148 fi
149 
150 if (( 1 )); then
151 
152  set_variable MODULE 101
153  set_variable PMT 0
154 
155  JPlot1D \
156  -f ${WORKDIR}/fit.root:${MODULE}.1gainspread \
157  -y "0.0 1.0" \
158  -\> "PMT" \
159  -\^ "gain spread" \
160  -T "" \
161  -o tot.${MODULE}.gainspread.$FORMAT $BATCH
162 
163 fi
detectorType
remove to overwrite exit fi JDetector D $ID detectorType
Definition: software/JDataFilter/JDataFilter.sh:115
JTOOLS::w
data_type w[N+1][M+1]
Definition: JPolint.hh:741
f
then JMuonPostfit f
Definition: JMuonPostfit.sh:118
JDATABASE::PBS::PMT
static const JPBS_t PMT(3, 4, 2, 3)
PBS of photo-multiplier tube (PMT)
WORKDIR
then usage $script< detector file >< inputfile > nUtility script to create PDF and CDF of transition time distribution fi case set_variable WORKDIR
Definition: JLegolas.sh:31
exit
exit
Definition: JPizza.sh:36
Y
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
Definition: detector-XY:fit1d.sh:33
echo
then echo
Definition: JDOMDAQDriver.sh:56
JPHYSICS::C
static const double C
Physics constants.
Definition: JPhysics/JConstants.hh:25
JTOOLS::n
const int n
Definition: JPolint.hh:660
set_variable
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
o
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
debug
int debug
debug level
Definition: JSirene.cc:63
usage
* usage
Definition: JLigier-local.sh:34
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
a
then JCalibrateToT a
Definition: JTuneHV.sh:116
$DETECTOR
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
Definition: module-Z:fit.sh:118
DEBUG
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
Definition: JCalibratePMT.sh:21
$WORKDIR
then display $WORKDIR
Definition: plot-Domino.sh:128
monitor
static void monitor()
Definition: DataQueue.cc:32
DETECTOR
set_variable DETECTOR
Definition: JARCACalibrateMuon.sh:31
detector
do set_variable DETECTOR_TXT $WORKDIR detector
Definition: detector-XY:fit1d.sh:38
in
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:41
D
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
version
version
Definition: JCalibratePMT.sh:7
Generated by   doxygen 1.8.5