Jpp  18.5.0
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JPlotNPE-PDG.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author lquinn
4 #
5 version=1.0
6 script=${0##*/}
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Utility script to plot integrals of tabulated PDFs.
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
25 
26 set_variable NORTH "0.00 0.00"
27 set_variable EAST "1.57 0.00"
28 set_variable SOUTH "3.14 0.00"
29 set_variable WEST "1.57 3.14"
30 
34 set_variable E 1.0
35 set_variable YMAX 1e+2
36 set_variable YMIN 1e-5
37 
38 set_variable PDF $JPP_DATA/J%p.dat
39 set_variable CDF $JPP_DATA/I%p.dat
40 
41 set_variable: FORMAT GRAPHICS_FORMAT gif
42 set_variable+ BATCH GRAPHICS_BATCH -B
43 
44 if do_usage $*; then
45  usage "$script [energy]"
46 fi
47 
48 case $# in
49  1) set_variable E $1;;
50 esac
51 
52 
53 for dir in NORTH EAST SOUTH WEST; do
54 
55  echo "Generating $dir."
56 
57  eval D=\$$dir
58 
59  if (( 1 )); then
60 
61  $DIR/JPlotNPE-PDG \
62  -P ${PDF/\%/${TYPE}} \
63  -C ${CDF/\%/${TYPE}} \
64  -D "${D}" \
65  -E "${E}" \
66  -o $WORKDIR/npe\[${dir}\].root \
67  -d $DEBUG --!
68 
69  fi
70 
71  if (( 0 )); then
72 
73  for FUNCTION in pdf npe cdf; do
74 
75  JPlot2D \
76  -f $WORKDIR/npe\[${dir}\].root:${FUNCTION} \
77  -z "$YMIN $YMAX" -Z \
78  -T "${dir} E = ${E} GeV" \
79  -O "colz" \
80  -\> "D [m]" \
81  -\< "cos #theta_{0}" \
82  -o ${FUNCTION}-${dir}-${E}GeV.$FORMAT $BATCH
83  done
84  fi
85 
86  if (( 1 )); then
87 
88  set_variable YMIN 0.90
89  set_variable YMAX 1.10
90 
91  for FUNCTION in npe cdf; do
92 
93  JOpera2D \
94  -f $WORKDIR/npe\[${dir}\].root:${FUNCTION} \
95  -f $WORKDIR/npe\[${dir}\].root:pdf \
96  -u Divide \
97  -o $WORKDIR/ratio-${FUNCTION}-${dir}.root
98 
99  JPlot2D \
100  -f $WORKDIR/ratio-${FUNCTION}-${dir}.root:Divide \
101  -z "$YMIN $YMAX" \
102  -T "${FUNCTION}/pdf ${dir}" \
103  -\> "D [m]" \
104  -\< "cos #theta_{0}" \
105  -O COLZ \
106  -o ratio-${FUNCTION}-${dir}.$FORMAT $BATCH
107  done
108  fi
109 done
&set_variable TYPE
Definition: getLs.sh:22
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR set_variable DEBUG set_variable PDF $JPP_DATA J p dat set_variable HISTOGRAM if do_usage *then usage $script< detector file >< inputfile >< outputfile > fi case set_variable PDF
then usage $script< input file >[option[primary[working directory]]] nWhere option can be E
Definition: JMuonPostfit.sh:40
version
Definition: JEditTuneHV.sh:5
o $QUALITY_ROOT d $DEBUG!CHECK_EXIT_CODE JPlot1D f
Definition: JDataQuality.sh:76
usr share Modules init zsh export TMPDIR
exit
Definition: JPizza.sh:36
then for dir in NORTH EAST SOUTH WEST
Definition: JDrawLED.sh:49
then
Definition: datalogs.sh:31
fatal Invalid number of arguments esac let YMIN
Definition: testToT.sh:39
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DETECTOR_ID set_array RUNS if do_usage *then usage $script< source directory > nThe source directory corresponds to a pre calibrated D0ARCA009 detector fi case set_variable DIR
then fatal Wrong number of arguments fi JConvertDetectorFormat a o
static const double C
Physics constants.
do set_variable OUTPUT_DIRECTORY $WORKDIR T
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
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 for FUNCTION in pdf npe cdf
Definition: JPlotNPE-PDG.sh:73
fi set_variable FILENAME $WORKDIR
Definition: datalogs.sh:34
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
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 JAcoustics sh $DETECTOR_ID source JAcousticsToolkit sh CHECK_EXIT_CODE typeset A EMITTERS get_tripods $WORKDIR tripod txt EMITTERS get_transmitters $WORKDIR transmitter txt EMITTERS for EMITTER in
Definition: JCanberra.sh:48
do set_variable MODULE getModule a $WORKDIR detector_a datx L $STRING JEditDetector a $WORKDIR detector_a datx M $MODULE setz o $WORKDIR detector_a datx JEditDetector a $WORKDIR detector_b datx M $MODULE setz o $WORKDIR detector_b datx done echo Output stored at $WORKDIR detector_a datx and $WORKDIR tripod_a txt JDrawDetector2D a $WORKDIR detector_a datx a $WORKDIR detector_b datx L BL o detector $FORMAT $BATCH JDrawDetector2D T $WORKDIR tripod_a txt T $WORKDIR tripod_b txt L BL o tripod $FORMAT $BATCH JCompareDetector a $WORKDIR detector_a datx b $WORKDIR detector_b datx o $WORKDIR abc root &dev null for KEY in X Y Z
double u[N+1]
Definition: JPolint.hh:865
then echo
Definition: JQAQC.sh:96
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:62
esac done
Definition: JAddHDE.sh:21
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62