Jpp  16.0.0-rc.1
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JMultiPMT.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 # Example script to study L0 and L0 hit probabilities.
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 WORKDIR $JPP_DATA
26 set_variable E_GEV 1e3
27 
28 set_variable NORTH "0.00 0.00"
29 set_variable EAST "4.71 0.00"
30 set_variable SOUTH "3.14 0.00"
31 set_variable WEST "1.57 0.00"
32 
33 set_variable: FORMAT GRAPHICS_FORMAT gif
34 set_variable+ BATCH GRAPHICS_BATCH -B
35 
36 if do_usage $*; then
37  usage "$script"
38 fi
39 
40 
41 # production of PDF tables
42 
43 if (( 0 )); then
44 
46 
47  make -C $JPP_DIR/software/JPhysics clean
48  make -C $JPP_DIR/software/JPhysics $DETECTOR
49 
50  make -p $WORKDIR/$DETECTOR/
51 
52  JSirene.sh -W $WORKDIR/$DETECTOR/ -P
53 
54  done
55 
56 # wait until all PDF tables have been produced
57 
58  exit
59 
60 fi
61 
62 if (( 0 )); then
63 
64  for DETECTOR in Antares KM3NeT; do
65 
66  JSirene.sh -W $WORKDIR/$DETECTOR/ -M
67 
68  done
69 fi
70 
71 
72 
73 if (( 1 )); then
74 
75  for DETECTOR in Antares KM3NeT; do
76 
77  for dir in NORTH EAST SOUTH WEST; do
78 
79  echo "Generating histograms for ${DETECTOR} direction ${dir}."
80 
81  eval D=\$$dir
82 
83  OUTPUT_FILE=\[${DETECTOR}\]_${dir}.root
84 
85  $JPP_DIR/examples/JPhysics/JMultiPMT \
86  -f $WORKDIR/${DETECTOR}/J%p.dat \
87  -o $OUTPUT_FILE \
88  -O $DETECTOR \
89  -E $E_GEV \
90  -D "$D" \
91  -d $DEBUG
92 
93  done
94  done
95 
96 fi
97 
98 
99 if (( 1 )); then
100 
101  for type in L0 L1; do
102 
103  for dir in NORTH EAST SOUTH WEST; do
104 
105  JPlot1D \
106  -f \[Antares\]_${dir}.root:${type} \
107  -f \[KM3NeT\]_${dir}.root:${type} \
108  -T "${dir}" -L TR \
109  -o ${type}-${dir}.$FORMAT $BATCH
110 
111  done
112  done
113 fi
do $JPP JMEstimator M
Definition: JMEstimator.sh:37
then usage E
Definition: JMuonPostfit.sh:35
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
exit
Definition: JPizza.sh:36
then for dir in NORTH EAST SOUTH WEST
Definition: JDrawLED.sh:49
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
case $OPTION in clean clean
do set_variable OUTPUT_DIRECTORY $WORKDIR T
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
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
* usage
then JConvertDetectorFormat a $DETECTOR[1] o
then for DETECTOR in Antares KM3NeT
Definition: JMultiPMT.sh:45
&set_variable OUTPUT_FILE
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
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 usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
static const JGeographicalLocation Antares(42, 48, 06, 10)
then display $WORKDIR
Definition: plot-Domino.sh:128
set_variable 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:42
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
version
Definition: JCalibratePMT.sh:7
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:62