Merge pull request swertz#20 from BrieucF/mis_template_fit

Mis template fit

analyzer/init.sh

@@ -16,8 +16,7 @@ else
 
     module load python/python27_sl6_gcc49
     module load gcc/gcc-4.9.1-sl6_amd64
-    #module load root/6.02.05-sl6_gcc49
-    source /home/fynu/swertz/soft/root-6.02.12/bin/thisroot.sh
+    module load root/6.06.02-sl6_gcc49
     export PYTHONPATH="/home/fynu/sbrochet/.local/lib/python2.7/site-packages:$PYTHONPATH"
 
 fi

analyzer/mis_cards/mis_templateFit.conf

@@ -0,0 +1,7 @@
+mode = template
+pseudoNumber = 1000
+anaCfg = examples/mis_template_fit_boxYields_allFloat.yml
+outFile = results/tests/globalFitAllFloat.root        
+
+params_0 = TT:1; DY:1; ZZ:1; WW:1; WZ:1; ZH:0
+params_1 = TT:1; DY:1; ZZ:1; WW:1; WZ:1; ZH:1

analyzer/mis_cards/mis_templateFit.yml

@@ -0,0 +1,69 @@
+# Datasets:
+
+datasets:
+  TT:
+    # 0 --> not fitted, 1 --> fitted
+    signal: 1
+    range: [0.5,1.5]
+    histname: TT_yields
+  DY:
+    signal: 1
+    range: [0.5,1.5]
+    histname: DY_yields
+  ZZ:
+    signal: 1
+    range: [-20,20]
+    histname: ZZ_yields
+  ZH:
+    signal: 1
+    range: [-50,50]
+    histname: ZH_yields
+  WW:
+    signal: 0
+    range: [-10,10]
+    histname: WW_yields
+  WZ:
+    signal: 0
+    range: [-10,10]
+    histname: WZ_yields
+
+#  OtG:
+#    # Fitted
+#    signal: 1
+#    # range: define fit range as (range[0]*nevents, range[1]*nevents)
+#    # the fit might crash if the total PDF gets negative for some fit parameters
+#    # : > has to be tuned
+#    range: [-0.2,0.2]
+#    histname: OtG_yields
+#  
+#  OG:
+#    signal: 1
+#    range: [-1.,1.]
+#    histname: OG_yields
+#  
+#  OC81qq:
+#    signal: 1
+#    range: [-20,20]
+#    histname: OC81qq_yields
+#  
+#  OC8ut:
+#    signal: 1
+#    range: [-50,50]
+#    histname: OC8ut_yields
+#
+# Analysis, i.e. fit configuration
+
+analysis:
+
+  # variable used for fitting 
+  inputvar: Box_yields
+
+  numcpu: 1
+  # options:
+  #      - get: get histogram (histname) from file (histfile)
+  #		   - verbose: print all the fit information (default off)
+  options: get
+
+  # If "get" mode: specify file containing all the histograms.
+  # If every process has a different file, put a "histfile" for each process INSTEAD of here
+  histfile: results/tests/MIS_hists.root

analyzer/python/mcstudy.py

@@ -147,7 +147,7 @@ def mcStudyMain(mcStudyFile):
         ROOT.gROOT.SetBatch(ROOT.kTRUE);
         if ROOT.gROOT.GetVersionInt() >= 60400:
             ROOT.gStyle.SetPalette(57)
-        else
+        else:
             ROOT.gStyle.SetPalette(54)
 
         cnv = ROOT.TCanvas("cnv_correlations", "Correlations", 900, 600)