Jpp  16.0.2
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JTransitTime.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 version=1.0
6 script=${0##*/}
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Utility script to plot transition time spread from JDETECTOR::JPMTAnalogueSignaProcessor.
11 #
12 # ------------------------------------------------------------------------------------------
13 
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 
21 source $JPP_DIR/setenv.sh $JPP_DIR
22 
23 
24 set_variable DEBUG 3
25 set_variable WORKDIR ${TMPDIR:-/tmp}/
26 set_variable DIR $JPP_DIR/examples/JDetector
27 set_variable OUTPUT_FILE $WORKDIR/tts%.root
28 set_variable NUMBER_OF_ENTRIES 10000000
29 set_variable: FORMAT GRAPHICS_FORMAT gif
30 set_variable+ BATCH GRAPHICS_BATCH -B
31 
32 
33 if do_usage $*; then
34  usage "$script [output file]"
35 fi
36 
37 case $# in
38  1) set_variable OUTPUT_FILE $1;;
39 esac
40 
41 typeset -A TTS
42 
43 for OPTION in 1 38 41 67; do
44 
45  $DIR/JTransitTime \
46  -o ${OUTPUT_FILE/\%/[${OPTION}]} \
47  -n 0 \
48  -O $OPTION \
49  -n $NUMBER_OF_ENTRIES \
50  -d $DEBUG
51 
52  JFit \
53  -f ${OUTPUT_FILE/\%/[${OPTION}]}:pmt \
54  -F "[0]*exp(-0.5*(x-[1])*(x-[1])/([2]*[2]))" \
55  -@ "p0 = 0.8*GetMaximum" \
56  -@ "p1 = 0.0" \
57  -@ "p2 = 2.0" \
58  -o $WORKDIR/fit.root >& /dev/null
59 
60  set_array P2 `JPrintFit -f $WORKDIR/fit.root:pmt -@p2`
61 
62  TTS[${OPTION}]=$P2[1]
63 
64  JPlot1D \
65  -f ${OUTPUT_FILE/\%/[${OPTION}]}:tts \
66  -f ${OUTPUT_FILE/\%/[${OPTION}]}:pmt \
67  -y "1e-5 1.0" -Y \
68  -\> "time [ns]" \
69  -\^ "probability" \
70  -L TR \
71  -T "" \
72  -o TTS_${OPTION}.$FORMAT $BATCH
73 done
74 
75 JPlot1D \
76  -f ${OUTPUT_FILE/\%/[1]}:tts \
77  -f ${OUTPUT_FILE/\%/[38]}:tts \
78  -f ${OUTPUT_FILE/\%/[41]}:tts \
79  -f ${OUTPUT_FILE/\%/[67]}:tts \
80  -y "1e-5 1.0" -Y \
81  -\> "time [ns]" \
82  -\^ "probability" \
83  -L TR \
84  -T "" \
85  -o TTS.$FORMAT $BATCH
86 
87 for key value in ${(@kv)TTS}; do
88  printf "option %3d sigma %6.3f [ns]\n" $key $value
89 done
OPTION
then usage $script< detector >< run >< outputfile > fi case set_variable OPTION
Definition: JTuna.sh:26
f
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
exit
exit
Definition: JPizza.sh:36
DIR
do echo Generating $dir eval DIR
Definition: JPlotNPE1D.sh:54
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
JTOOLS::n
const int n
Definition: JPolint.hh:676
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
TMPDIR
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR if do_usage *then usage $script[(input file)+] fi set_variable DEBUG set_variable WORKDIR TMPDIR
Definition: JEnergyCorrection.sh:21
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
done
set_variable NUMBER_OF_ITERATIONS set_variable EPSILON cat acoustics_fit_parameters txt<< EOF $CONFIGURATION[*]Nmin=3;sigma_s=100.0e-6;stdev=10.0;mestimator=0;fixStrings=0;EOF for STRING in $STRINGS[*];do#fit stretching and(z) position of given string set_variable DETECTOR_TMP ${TMPDIR:-/tmp}/detector_A.datx JEditDetector-a $DETECTOR-o $DETECTOR_TMP-r $STRING JEditDetector-a $DETECTOR-o $DETECTOR-k $STRING for MUL in 0.005 0.001;do DX_M=0.2 for((N=0;$N< $NUMBER_OF_ITERATIONS;++N));do CHI2[3]=$CHI2[1] fitPositionOfString $STRING Z $DX_M fitStretchingOfString $STRING $MUL if(($CHI2[3]-$CHI2[1]< $EPSILON));then break fi done if(($N >=$NUMBER_OF_ITERATIONS));then printf"warning: reached maximum number of iterations %d - converenge %7.3f\n"$N $(($CHI2[3]-$CHI2[1])) fi done JMergeDetector-a $DETECTOR-a $DETECTOR_TMP-o $DETECTOR rm-f $DETECTOR_TMP JConvertDetectorFormat-a $DETECTOR-o $DETECTOR-r-d 0 > &dev null done
Definition: acoustics-fit.sh:563
WORKDIR
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR
Definition: JLegolas.sh:20
F
then awk F
Definition: getMIIseacurrent.sh:86
usage
* usage
Definition: JLigier-local.sh:34
o
then JConvertDetectorFormat a $DETECTOR[1] o
Definition: JCompareDetector.sh:45
p2
p2
Definition: module-Z:fit.sh:74
OUTPUT_FILE
&set_variable OUTPUT_FILE
Definition: JARCACalibrateMuon.sh:29
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
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
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:42
A
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
Definition: JShellParser.csh:15
version
version
Definition: JCalibratePMT.sh:7
$OPTION
*fatal Wrong option $OPTION
Definition: getMIIseacurrent.sh:66
null
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
Definition: examples/JLogger/JLogger.sh:33
Generated by   doxygen 1.8.5