Jpp
JORCAShowerPositionPDF.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author adomi
4 #
5 
6 version=1.0
7 script=${0##*/}
8 
9 # ------------------------------------------------------------------------------------------
10 #
11 # Utility script to create a ROOT TH2D PDF for E.M. Shower Position Reconstruction.
12 #
13 # ------------------------------------------------------------------------------------------
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 source $JPP_DIR/setenv.sh $JPP_DIR
21 
22 if ( do_usage $* ); then
23  usage "$script [detector file [input file [output file]]]"
24 fi
25 
26 SUB=${${RUN}:0:1}
27 
28 DetectorFile="/afs/in2p3.fr/throng/km3net/detectors/orca_115strings_av23min20mhorizontal_18OMs_alt9mvertical_v1.detx"
29 
30 set_variable WORKDIR $JPP_DATA
31 set_variable FINAL_HISTOGRAM $WORKDIR/PDFShowerPosition.root
32 set_variable DETECTOR $DetectorFile
33 set_variable DEBUG 1
34 
35 # DOWNLOAD FILES FROM IRODS AND EXTRACT THEM
36 
37 MCFilesIRODS=/in2p3/km3net/mc/atm_neutrino/KM3NeT_ORCA_115_23m_9m/v1.1.1/JTE/elec-CC_3-100GeV_1-50.root.tar.gz
38 
39 MCFiles=elec-CC_3-100GeV_1-50.root.tar.gz
40 
41 iget $MCFilesIRODS $WORKDIR
42 
43 tar -xvzf $WORKDIR/$MCFiles -C $WORKDIR
44 
45 FILENAME=JTE_elec-CC_3-100GeV_
46 
47 for file in $(/bin/ls $WORKDIR/$FILENAME*)
48 do
49 
50 nrun=$((nrun + 1))
51 
52 echo $nrun
53 
54 set_variable INPUT_FILE $file
55 set_variable OUTPUT_HISTOGRAM $WORKDIR/shower_position_histogram_$nrun.root
56 
57 case $# in
58  3) OUTPUT_HISTOGRAM=$3;&
59  2) INPUT_FILE=$2;&
60  1) DETECTOR=$1;;
61 esac
62 
63 print_variable DETECTOR INPUT_FILE OUTPUT_HISTOGRAM
64 check_input_file $DETECTOR $INPUT_FILE
65 check_output_file $OUTPUT_HISTOGRAM
66 
67 if (( 1 )); then
68 
69 # default input values
70  set_variable T_RES_NS "-100 100"
71  set_variable D_M "0 80"
72  set_variable COS_T "-1 1"
73 
74  timer_start
75 
76  JORCAShowerPositionPDF \
77  -a $DETECTOR \
78  -f $INPUT_FILE \
79  -o $OUTPUT_HISTOGRAM \
80  -T ${T_RES_NS} \
81  -D ${D_M} \
82  -C ${COS_T} \
83  -d ${DEBUG} --!
84 
85  timer_stop
86  timer_print
87 
88 fi
89 
90 done
91 
92 rm $FINAL_HISTOGRAM
93 
94 hadd $WORKDIR/$FINAL_HISTOGRAM $WORKDIR/shower_position_histogram_*.root
95 
96 rm $WORKDIR/shower_position_histogram_*.root
97 
98 # CONVERT HISTOGRAM INTO A PDF (NORMALISE IT)
99 
100 set_variable OPTION "NXYW"
101 
102 $JPP_DIR/software/JGizmo/JConvertToPDF2D \
103  -O $OPTION \
104  -f $FINAL_HISTOGRAM":hPDF2Dist" \
105  -o $OUTPUT_PDF
106 
107 rm $WORKDIR/$FILENAME*
108