9 # ------------------------------------------------------------------------------------------
11 # Utility script to create a ROOT TH2D PDF for E.M. Shower Position Reconstruction.
13 # ------------------------------------------------------------------------------------------
15 if [ -z $JPP_DIR ]; then
16 echo "Variable JPP_DIR undefined."
20 source $JPP_DIR/setenv.sh $JPP_DIR
22 if ( do_usage $* ); then
23 usage "$script [detector file [input file [output file]]]"
28 DetectorFile="/afs/in2p3.fr/throng/km3net/detectors/orca_115strings_av23min20mhorizontal_18OMs_alt9mvertical_v1.detx"
30 set_variable WORKDIR $JPP_DATA
31 set_variable FINAL_HISTOGRAM $WORKDIR/PDFShowerPosition.root
32 set_variable DETECTOR $DetectorFile
35 # DOWNLOAD FILES FROM IRODS AND EXTRACT THEM
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
39 MCFiles=elec-CC_3-100GeV_1-50.root.tar.gz
41 iget $MCFilesIRODS $WORKDIR
43 tar -xvzf $WORKDIR/$MCFiles -C $WORKDIR
45 FILENAME=JTE_elec-CC_3-100GeV_
47 for file in $(/bin/ls $WORKDIR/$FILENAME*)
54 set_variable INPUT_FILE $file
55 set_variable OUTPUT_HISTOGRAM $WORKDIR/shower_position_histogram_$nrun.root
58 3) OUTPUT_HISTOGRAM=$3;&
63 print_variable DETECTOR INPUT_FILE OUTPUT_HISTOGRAM
64 check_input_file $DETECTOR $INPUT_FILE
65 check_output_file $OUTPUT_HISTOGRAM
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"
76 JORCAShowerPositionPDF \
79 -o $OUTPUT_HISTOGRAM \
94 hadd $WORKDIR/$FINAL_HISTOGRAM $WORKDIR/shower_position_histogram_*.root
96 rm $WORKDIR/shower_position_histogram_*.root
98 # CONVERT HISTOGRAM INTO A PDF (NORMALISE IT)
100 set_variable OPTION "NXYW"
102 $JPP_DIR/software/JGizmo/JConvertToPDF2D \
104 -f $FINAL_HISTOGRAM":hPDF2Dist" \
107 rm $WORKDIR/$FILENAME*