Macro files for the X(3872) Lineshape Scan Study RN-QCD-2016-002
-----------------------------------------------------------------

* Macros Fig...C: These macros reproduce the plots in the note RN-QCD-2016-002 with the corresponding plot number
  Just setup root and execute 'root -l Fig...C' in a shell; they require the macro PandaSmartLabel.C
  
  - Fig06abc : m(2mu) vs. m(2pi) for signal (left), generic bkg (middle), NR J/psi 2pi bkg (right); masses after 4C fit
  
  - Fig06def : m(2mu) + m(2pi) vs. m(2pi) for signal (left), generic bkg (middle), NR J/psi 2pi bkg (right); masses after 4C fit 
  
  - Fig06ghi : m(2e) vs. m(2pi) for signal (left), generic bkg (middle), non-resonant J/psi 2pi bkg (right) ; masses after 4C fit
  
  - Fig06jkl : m(2e) + m(2pi) vs. m(2pi) for signal (left), generic bkg (middle), NR J/psi 2pi bkg (right); masses after 4C fit 
  
  - Fig07abc : m(2mu2pi) after preselection for signal(black) + generic bkg(red);
               m(2mu) after presel. for signal(black) and generic bkg(red) with gaussian fit (mu = 3097.4 MeV and sigma = 7.0 MeV) to core of signal(blue)
               m(2mu) after final selection for signal(black), generic bkg(red -> dashed line scaled according to cross section), non-resonant J/psi 2pi bkg(blue); masses after 4C fit  
               
  - Fig07def : m(2e2pi) after preselection for signal(black) + generic bkg(red);
               m(2e) after presel. for signal(black) and generic bkg(red) with gaussian fit (mu = 3094.7 MeV and sigma = 8.4 MeV) to core of signal(blue);
               m(2e) after final selection for signal(black), generic bkg(red -> dashed line scaled according to cross section), non-resonant J/psi 2pi bkg(blue); masses after 4C fit 
  
  - Fig13ab  : m(2mu) candidate distribution for signal (left) with double gaussian fit (blue) and generic bkg (right) with pol2 fit (red)
 
  - Fig14b   : example scan graph for BW study, input Gamma_0 = 130keV, 20 scan positions x 2 days data taking, HESRr mode 
  
  - Fig14c   : distribution of Gamma_fit - Gamma_0 for 300 toy experiments like shown on Fig14b
  
  - Fig15ab  : sensitivity (left) and bias (right) of BW study for various input Gamma_0, HESRr (black) and HL (red) modes
  
  - Fig15c   : measured E_f vs. input E_f,0 for molecule study in HESRr mode; magenta lines = threshold E_f,th = -8.56 MeV of bound/virtual state; blue line f(x) = x
  
  - Fig15d   : measured E_f vs. input E_f,0 for molecule study in HL mode; magenta lines = threshold E_f,th = -8.56 MeV of bound/virtual state; blue line f(x) = x
  
  - Fig15e   : mis-identification probability for nature of state as function of input E_f,0 in HL(red) and HESRr(black) modes
  
  - Fig16a   : example scan graph for BW study with Gamma_0 = 100keV, HESRr mode
  
  - Fig16b   : example scan graph for BW study with Gamma_0 = 500keV, HESRr mode
  
  - Fig16c   : example scan graph for BW study with Gamma_0 = 50keV, HL mode
  
  - Fig16d   : scan graph for BW study to be compared with Martin Galuskas analysis -> dp/p=2E-5, 20 scan points x 2 days, HESRr lumi, Gamma_0 = 100keV  
  
  - Fig17a   : example scan graph for molecule study with E_f,0 = -7MeV, HESRr mode

  - Fig17b   : example scan graph for molecule study with E_f,0 = -10MeV, HESRr mode

  - Fig17b   : example scan graph for molecule study with E_f,0 = -7.5MeV, HL mode

  - Fig17b   : scan graph for molecule study to be compared with Martin Galuskas analysis -> dp/p=2E-5, 20 scan points x 2 days, HESRr lumi, E_f,0 = -8.3MeV
  
  - Fig18ab  : sensitivity (left) and bias (right) of BW study for various input Gamma_0, HL (red) mode only (cf. Fig15ab)
  
  - Fig18cd  : sensitivity (left) and bias (right) of BW study for various input Gamma_0, HESRr (red) mode only (cf. Fig15ab)
  
  - Fig18e   : mis-identification probability for nature of state as function of input E_f,0, HL (red) mode only (cf. Fig15e)
  
  - Fig18f   : mis-identification probability for nature of state as function of input E_f,0, HESRr (red) mode only (cf. Fig15e)
  
  

* xscan_final2.C: This macro was used to simulate a scan study

  USAGE:
  void xscan_final2(<xsec>, <dEHR>, <steps>, <days>, <ntoy>, <nsave>, <fitmode>, <opt>, <pref>)
  
    <xsec>     : X peak cross section in nb;      default = 100
    <dEHR>     : HR mode dE in keV;               default = 83.9
    <steps>    : number of scan energy positions; default = 40
    <days>     : days of data taking/position;    default = 2
    <nsave>    : # scan example plots to save;    default = 1
    <fitmode>  : 0=cutncount, 1=ML J/psi;         default = 1
    <opt>      : option string;                   default = "BW Mol"
                 'BW' = Breit Wigner study, 
                 'Mol' = molecule study, 
                 'fixwin' = fixed scan window, 
                 '!HL' = no HL mode, 
                 '!HR' = no HR mode, 
                 '!HESRr' = no HESRr mode
    <pref>     : prefix string for output;        default = "TestScan_"
    

* xplots_final.C: Macro to produce the reco plots from data. Requires the analysis ROOT files, located in /u/kgoetzen/data/Xscan/