Jpp  15.0.2
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JMultiplicityK40.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 make coincidence rate plots.
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 DIR $JPP_DIR/examples/JPhysics/
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
36 esac
37 
38 if [[ ! -f $WORKDIR/k40\[KM3NeT\].root ]]; then
39 
40  $DIR/JMultiplicityK40 \
41  -n 100000000 \
42  -F 4 \
43  -G 2 \
44  -D "0.216 5" \
45  -A 1 \
46  -U \
47  -o $WORKDIR/k40\[KM3NeT\].root
48 fi
49 
50 for (( a = 0.0; $a <= 1.0; a += 0.5 )); do
51 
52  set_variable LABEL `printf '%3.1f' $a`
53 
54  if [[ ! -f $WORKDIR/k40\[${LABEL}\].root ]]; then
55 
56  $DIR/JMultiplicityK40 \
57  -n 100000000 \
58  -F 4 \
59  -G 2 \
60  -D "0.216 5" \
61  -A 3 -a $a \
62  -U \
63  -o $WORKDIR/k40\[${LABEL}\].root
64  fi
65 done
66 
67 INPUT_FILES=(`ls $WORKDIR/k40\[*\].root`)
68 
70 
71 JPlot1D \
72  -f "${INPUT_FILES//\.root/.root:pmt}" \
73  -\> "cos(#theta)" \
74  -\^ "acceptance" \
75  -y "0 1.6" \
76  -L TR \
77  -T "PMT" \
78  -o pmt.$FORMAT $BATCH
79 
80 M=2
81 
82 JPlot1D \
83  -f "${INPUT_FILES//\.root/.root:M.${M}}" \
84  -\> "D [m]" \
85  -\^ "rate [Hz]" \
86  -x "0 5" \
87  -y "1 1e2" -Y \
88  -L TR \
89  -T "M = ${M}" \
90  -o k40.$FORMAT $BATCH
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR set_variable DEBUG set_variable WORKDIR if do_usage *then usage for INPUT_FILE in $INPUT_FILES[*]
do $JPP JMEstimator M
Definition: JMEstimator.sh:37
then JMuonPostfit f
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
Definition: JPizza.sh:36
do echo Generating $dir eval DIR
Definition: JPlotNPE1D.sh:54
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
const int n
Definition: JPolint.hh:660
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 for DETECTOR in Antares KM3NeT
Definition: JMultiPMT.sh:45
then JCalibrateToT a
Definition: JTuneHV.sh:116
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:128
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
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/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
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
version
Definition: JCalibratePMT.sh:7
source $JPP_DIR setenv csh $JPP_DIR eval JShellParser o a A