Jpp  15.0.1
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JFilter.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 test trigger algorithms.
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 
25 set_variable WORKDIR $JPP_DATA
26 set_variable INPUT_FILE $WORKDIR/trigger_efficiency+background.root
27 set_variable OUTPUT_FILE $WORKDIR/filter\%.root
28 set_variable DETECTOR $WORKDIR/km3net_reference.detx
29 set_variable HISTOGRAM 1
30 
31 
32 if do_usage $*; then
33  usage "$script [detector file [input file [output file]]]"
34 fi
35 
36 case $# in
40 esac
41 
42 
43 if (( 1 )); then
44 
45  check_input_file $DETECTOR $INPUT_FILE
46 
47  for OPTION in A B C D E F; do
48 
49  FILE="${OUTPUT_FILE/\%/[${OPTION}]}"
50 
51  print_variable DETECTOR INPUT_FILE FILE
52 
53  $JPP_DIR/examples/JTrigger/JFilter \
54  -a $DETECTOR \
55  -f $INPUT_FILE \
56  -o $FILE \
57  -c $OPTION \
58  -T 20 \
59  -R 200 \
60  -H $HISTOGRAM \
61  -d $DEBUG
62 
63  done
64 
65 fi
66 
67 
68 if (( 1 )); then
69 
70  X_LABEL[1]="Number of hits"
71  X_LABEL[2]="E_{#nu} [GeV]"
72 
73  X_OPTION[1]=""
74  X_OPTION[2]="-XX"
75 
76  X_RANGE[1]="-x 1 50"
77  X_RANGE[2]="-x 2 7"
78 
79  FILE_E=""
80  FILE_P=""
81 
82  for OPTION in A B C D E F; do
83  FILE_E+=" ${OUTPUT_FILE/\%/[${OPTION}]}:he"
84  FILE_P+=" ${OUTPUT_FILE/\%/[${OPTION}]}:hp"
85  done
86 
87  FILE_E=${FILE_E## }
88  FILE_P=${FILE_P## }
89 
90  JPlot1D \
91  -f $FILE_E \
92  -y "0 1.1" \
93  -\^ "efficiency" \
94  -\> $X_LABEL[$HISTOGRAM] \
95  $X_OPTION[$HISTOGRAM] \
96  $X_RANGE[$HISTOGRAM] \
97  -L TR
98 
99 
100  JPlot1D \
101  -f $FILE_P \
102  -y "0 1.1" \
103  -\^ "purity" \
104  -\> $X_LABEL[$HISTOGRAM] \
105  $X_OPTION[$HISTOGRAM] \
106  $X_RANGE[$HISTOGRAM] \
107  -L TR
108 
109 fi
then usage $script< detector >< run >< outputfile > fi case set_variable OPTION
Definition: JTuna.sh:26
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
set_variable INPUT_FILE
static const double H
Planck constant [eV s].
exit
Definition: JPizza.sh:36
then usage E
Definition: JMuonPostfit.sh:35
then echo
static const double C
Physics constants.
do FILE
Definition: JFilter.sh:49
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
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
then usage $script[distance] fi case set_variable R
Definition: JDrawLED.sh:43
&set_variable OUTPUT_FILE
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
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
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
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: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
then usage eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY kill_child_processes_at_exit attach getModule a $DETECTOR typeset Z STRING typeset Z FLOOR typeset A X_LABEL typeset A YMAX X_LABEL[twist]
Definition: plot-compass.sh:45
version
Definition: JCalibratePMT.sh:7
*fatal Wrong option $OPTION
source $JPP_DIR setenv csh $JPP_DIR eval JShellParser o a A