We have the following NLOPS predictions.


============================
======     VBF       =======
============================

1. Powheg Box plus Pythia (default shower)
   A. filename = pythia/pythia-pThard2.yoda
      the nominal weights correspond to the central choice.
   B. [VAR{1..6}] corresponds to the other muF,muR variations in the ME.
      An envelope must be built out of these 7 curves.
   C. Starting shower scale variation:  
      central value in pythia/pythia-pThard0.yoda
  (D.) Renormalisation shower scale variations in pythia/pythia-pThard2-psvar.yoda
      weights [PS_MUR0.5] and [PS_MUR2.0]   --> DO NOT INCLUDE

   ==> Show the central prediction with a solid line (A),
       PS starting scale variation with a dashed line (C),
       envelope for ME muF, muR variations (B)
       envelope obtained summing in quadrature (B) and |(A)-(C)|


2. Powheg Box plus Vincia
   As 1 with pythia --> vincia in the filenames.
   Show only the central prediction.

3. Powheg Box plus Herwig (angular-ordered)
   A. filename = powheg-herwig/pwg-hw7-vbf.yoda
      the nominal weights correspond to the central choice.
   B. [W{1..8}] corresponds to the other muF,muR variations in the ME.
      An envelope must be built out of these curves. (Central weights is repetead more the once, so there are only 7 vars)
   Show only the central prediction

4. Madgraph plus Herwig (angular-ordered)
   A. filename = MG5_aMCplusH7/run_mushdefault/histos.yoda
      the nominal weights correspond to the central choice.
   B. [tag=             AAA dyn=  10 muR=XXX muF=YYY]
      with AAA = 600 (LO), 602 (NLO correction), 0 (NLO)
      build the 7 point variations envelope from tag = 0 
   C. PS scale variations from central value in 
      MG5_aMCplusH7/{run_mush0.5,run_mush2}/histos.yoda
  Show only the central prediction

5. Herwig angular ordered standalone [analysis name MY_ANALYSIS -> MY_ANALYSIS_proVBF]
   A. filename = hw7/yh-runs/VBF-AO-MCatNLO.yoda  
   B. hw7/yh-runs/VBF-AO-MCatNLO-{Qdown,Qup}.yoda corresponds to the correlated muF,muR variations
   C. hw7/yh-runs/VBF-AO-MCatNLO-{down,up}.yoda corresponds to the shower starting scale variation
  Show only the central value

6. Herwig dipole standalone [analysis name MY_ANALYSIS -> MY_ANALYSIS_proVBF]
As 5 with AO --> Dipole in the file name

7. Sherpa default
   A. filename = Sherpa/Sherpa-VBF.yoda.gz
      The nominal weights correspond to the central choice.
   B. [ME.MUR=*__ME.MUF=*__ME.LHAPDF=93300] correspond to muR/muF in the ME
   C. Sherpa/Sherpa-VBF-SclVar{Up,Down}.yoda.gz corresponds to the shower starting scale variation
   (D.) [MUR=*__ME.MUF=*__LHAPDF=93300] correspond to (correlated) muR/muF variations in the ME and in the shower. In
       (D.1) Sherpa-VBF.yoda.gz                 <-- default variation (too aggressive)
       (D.2) Sherpa-VBF-largeShowerVar.yoda.gz  <-- LL scale variation (too big)
       (D.3) Sherpa-VBF-largeShowerVarV{2,3}.yoda.gz <-- NLL scale variation (they are equivalent, take e.g. 2 )
  Show only the central value

8. Sherpa Dire
   As 7, with "Sherpa/Sherpa" --> "Sherpa/Sherpa-Dire".
   (D.2 -- D.3) are missing


Comparison between NLOPS and fixed order ---> NLOPS-1 plots

Maintain the same colouring and stiles for the main panel and the first ratio plot, with powheg+pythia the central prediction (1.A)


The "matching&PS scheme" envelope showed in the main and first ratio plot is given by summing in quadrature these 3 envelopes:
   - muF, muR variation in the ME for the central prediction (1.B)
   - shower starting scale variation for the central prediction (1.C)
   - envelope of all the central predictions in 1 -- 8

In the second ratio plot show the breakdown:
- A solid coloured line for each NLOPS curve (as already done)
- A dashed coloured line with the powheg pythia colour for the different pTmax match (1.C) [label it pTmaxMatch=0]
- A solid filled band for the ME scale variation [add a label = "ME $\mu_{F,R}$ vars"]

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

Introduce a new set of plots, NLOPS-1-2, for VBF.


All the curves are normalised wrt Sherpa Default (7.A)

First panel:

- Show Sherpa default curve (7A) in red (as in the other NLOPS), solid
- Show the other 2 curves corresponding to shower starting scale variations as a dashed and a dotted line (7.C), always in red
- Show the pure muF,muR variation in the ME (7.B) as a filled band, in red

Second panel:
- Show a filled band for the correlated muF,muR variations in the shower and in the ME (7.D.1 [deep pink], 7.D.2 [violet], 7.D.3 [blue]) 

Third panel:
- Show the hashed green band corresponding to the total NLOPS error built for NLOPS-1 plots
- Show a solid red band summing in quadrature all the errors of the firt panel (muF,muR vars in ME + shower starting scale var)
- Show a solid blue band summing in quadrature the blue envelope in the previous panel (muF,muR vars in ME+PS) plus the envelope obtained varying the PS scale vars