ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/cbrown/Development/DistributedModelCalculations/ShapeLimits/CreateModel.sh
(Generate patch)

Comparing UserCode/cbrown/Development/DistributedModelCalculations/ShapeLimits/CreateModel.sh (file contents):
Revision 1.3 by buchmann, Fri Apr 27 07:09:25 2012 UTC vs.
Revision 1.4 by buchmann, Mon May 7 16:43:30 2012 UTC

# Line 6 | Line 6 | if [ $# -lt 2 ]; then
6   fi
7  
8   echo -e "Have been supplied with the following arguments \n $1 \n $2 \n $3"
9 < ORIGDIR=`pwd`
9 > ORIGDIR="$( cd "$( dirname "${BASH_SOURCE[0]}" )" && pwd )"
10 >
11   cd ~/final_production_2011/CMSSW_4_2_8/src/HiggsAnalysis/
12   origscramarch=$SCRAM_ARCH
13   origbase=$CMSSW_BASE
# Line 53 | Line 54 | echo "                         ~~~~~~~~~
54   #$CMSSW_BASE/bin/$SCRAM_ARCH/combine -M HybridNew --frequentist --testStat LHC $1/model.root -H ProfileLikelihood
55  
56   echo "Ready to run limit!"
57 +
58 + modelsource="model.root"
59 + if [[ $5 -gt 0 ]]; then
60 +  modelsource="susydatacard.txt"
61 + fi
62 +  
63   if [[ $3 -gt 1 ]]; then
64   echo "Asymptotic is $3 (0) therefore doing FULL LIMITS"
65   # command="$CMSSW_BASE/bin/$SCRAM_ARCH/combine -M HybridNew -t1 --frequentist --testStat LHC $1/model.root -H ProfileLikelihood"
66 < command="$CMSSW_BASE/bin/$SCRAM_ARCH/combine -M HybridNew --rule CLs --testStat LHC --rMin $5 --rMax $6 --freq -T 1000 -i 8 --saveHybridResult --saveToys -v 1 -s 0 $1/model.root "
66 > command="$CMSSW_BASE/bin/$SCRAM_ARCH/combine -M HybridNew --rule CLs --testStat LHC -H ProfileLikelihood --freq -T 1000 -i 8 --saveHybridResult --saveToys -v 1 -s 0 $1/${modelsource} "
67   else
68   if [[ $3 -lt 1 ]]; then
69      echo "Asymptotic is $3 (0) therefore doing ASYMPTOTIC LIMITS"
70 <    command="$CMSSW_BASE/bin/$SCRAM_ARCH/combine -M Asymptotic  -t1 $1/model.root -H ProfileLikelihood"
70 >    command="$CMSSW_BASE/bin/$SCRAM_ARCH/combine -M Asymptotic  -t1 $1/${modelsource} -H ProfileLikelihood"
71   else
72      echo "Asymptotic is $3 (0) therefore doing ASYMPTOTIC LIMITS but WITH initial guess"
73 <    command="$CMSSW_BASE/bin/$SCRAM_ARCH/combine -M Asymptotic  --rMax $6 -t1 $1/model.root -H ProfileLikelihood"
73 >    command="$CMSSW_BASE/bin/$SCRAM_ARCH/combine -M Asymptotic  -t1 $1/${modelsource} -H ProfileLikelihood"
74   fi
75   fi
76  
# Line 94 | Line 101 | g++ $ORIGDIR/ReadAndSave.C -o ReadAndSav
101   if [[ $3 -gt 1 ]]; then
102    startfile=$outfile
103    outfile="allmerged.root"
104 <  $CMSSW_BASE/bin/$SCRAM_ARCH/combine $1/model.root -M HybridNew --freq --grid=$startfile --expectedFromGrid 0.16
105 <  $CMSSW_BASE/bin/$SCRAM_ARCH/combine $1/model.root -M HybridNew --freq --grid=$startfile --expectedFromGrid 0.025
106 <  $CMSSW_BASE/bin/$SCRAM_ARCH/combine $1/model.root -M HybridNew --freq --grid=$startfile --expectedFromGrid 0.5
107 <  $CMSSW_BASE/bin/$SCRAM_ARCH/combine $1/model.root -M HybridNew --freq --grid=$startfile --expectedFromGrid 0.84
108 <  $CMSSW_BASE/bin/$SCRAM_ARCH/combine $1/model.root -M HybridNew --freq --grid=$startfile --expectedFromGrid 0.975
109 <  hadd $outfile *.root
104 >  $CMSSW_BASE/bin/$SCRAM_ARCH/combine $1/${modelsource} -M HybridNew --freq --grid=$startfile --expectedFromGrid 0.16
105 >  $CMSSW_BASE/bin/$SCRAM_ARCH/combine $1/${modelsource} -M HybridNew --freq --grid=$startfile --expectedFromGrid 0.025
106 >  $CMSSW_BASE/bin/$SCRAM_ARCH/combine $1/${modelsource} -M HybridNew --freq --grid=$startfile --expectedFromGrid 0.5
107 >  $CMSSW_BASE/bin/$SCRAM_ARCH/combine $1/${modelsource} -M HybridNew --freq --grid=$startfile --expectedFromGrid 0.84
108 >  $CMSSW_BASE/bin/$SCRAM_ARCH/combine $1/${modelsource} -M HybridNew --freq --grid=$startfile --expectedFromGrid 0.975
109 >  hadd $outfile higgs*.root
110   fi
111  
112   ./ReadAndSave.exec $1/$outfile $1/ShapeDropletResult.txt

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines