Jpp  test_elongated_shower_pde
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JGedanken.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 version=1.0
5 script=${0##*/}
6 
7 # ------------------------------------------------------------------------------------------
8 #
9 # Auxiliary script for Gedanken experiment.
10 #
11 # ------------------------------------------------------------------------------------------
12 
13 if [ -z $JPP_DIR ]; then
14  echo "Variable JPP_DIR undefined."
15  exit
16 fi
17 
18 source $JPP_DIR/setenv.sh $JPP_DIR
19 
20 set_variable: DEBUG SIRENE_DEBUG 2
21 set_variable: WORKDIR SIRENE_WORKDIR ./
23 
24 # detector
25 
26 set_variable MODULE 25.0 0.0 50.0 # position
27 
28 # particle
29 
30 set_variable E_GEV 20 # energy
31 set_variable POS 0.0 0.0 0.0 # position
32 set_variable DIR 0.0 0.0 1.0 # direction
33 set_variable TYPE 13 # PDG
34 
35 # graphics
36 
37 set_variable: FORMAT GRAPHICS_FORMAT gif
38 set_variable+ BATCH GRAPHICS_BATCH -B
39 
40 if do_usage $*; then
41  usage "$script <output file> [number of events]"
42 fi
43 
44 case $# in
47 esac
48 
50 
51 JModule \
52  -D 1001 \
53  -M 1 \
54  -P "$MODULE" \
55  -o $DETECTOR \
56  -d $DEBUG
57 
58 $JPP_DIR/examples/JSirene/JGedanken \
59  -o $WORKDIR/gedanken.root \
60  -P "$POS" \
61  -D "$DIR" \
62  -E $E_GEV \
63  -T $TYPE \
64  -n $NUMBER_OF_EVENTS \
65  -d $DEBUG
66 
67 JSirene \
68  -a $DETECTOR \
69  -f $WORKDIR/gedanken.root \
70  -o $WORKDIR/sirene.root \
71  -F $JPP_DATA/I%p.dat \
72  -k \
73  -N 0 \
74  -d $DEBUG
75 
76 $JPP_DIR/examples/JSirene/JDomino \
77  -a $DETECTOR \
78  -f $WORKDIR/sirene.root \
79  -o $WORKDIR/domino.root \
80  -d $DEBUG
81 
82 JPlot1D \
83  -f $WORKDIR/"domino.root:h\[%\]" \
84  -y "1.0e-6 1.0e0" -Y \
85  -\> "#Deltat [ns]" \
86  -\^ "dP/dt npe/ns" \
87  -T "" \
88  -O HIST \
89  -o pdf.$FORMAT $BATCH
90 
91 for i in {1..31}; do
92 
93  JPlot1D \
94  -f $WORKDIR/"domino.root:h\[${i}\]" \
95  -y "1.0e-6 1.0e0" -Y \
96  -\> "#Deltat [ns]" \
97  -\^ "dP/dt npe/ns" \
98  -T "PMT[${i}]" \
99  -O HIST \
100  -o ${TMPDIR:-/tmp}/pdf_${i}.$FORMAT $BATCH
101 done
102 
103 montage \
104  -tile 8x4 \
105  -geometry +0+0 \
106  ${TMPDIR:-/tmp}/pdf_{1..31}.$FORMAT \ PDF.$FORMAT >& /dev/null
107 
JDetector detector
Definition: JRunAnalyzer.hh:23
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
then usage $script< input_file >< detector_file >< output_file(.root)> fi set_variable NUMBER_OF_TIMESLICES set_variable NUMBER_OF_SUMMARYSLICES set_variable NUMBER_OF_EVENTS set_variable ALL_PLOTS case set_variable NUMBER_OF_EVENTS
Definition: JRunAnalyzer.sh:18
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
then JShowerPostfit f $INPUT_FILE o $OUTPUT_FILE N
exit
Definition: JPizza.sh:36
do echo Generating $dir eval DIR
Definition: JPlotNPE1D.sh:54
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTORXY_TXT $WORKDIR $DETECTORXY_TXT tail read X Y CHI2 RMS printf optimum n $X $Y $CHI2 $RMS awk v Y
then echo
const int n
Definition: JPolint.hh:676
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR if do_usage *then usage $script[(input file)+] fi set_variable DEBUG set_variable WORKDIR TMPDIR
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
then awk F
* usage
then JConvertDetectorFormat a $DETECTOR[1] o
then for APP in event gandalf start energy
Definition: JMuonMCEvt.sh:44
&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:42
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
version
Definition: JCalibratePMT.sh:7
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:62