// Author: Ricardo Eusebi, FNAL
// phone: (630) 840 4326 
// email: eusebi@fnal.gov
// Created: sometime in 2004
// This code is in the process of being PUBLISHED.
// DO NOT PUBLISH OR DISTRIBUTE WITHOUT THE CONSENT OF RICARDO EUSEBI

#include "ParameterDefinitions.C"

//C++ libraries
#include <iostream.h>


Parameters GetBenchmark(Int_t Bench){

  //Define the parameters for this run.

  //The resulting value
  Parameters param;

  //This are typically redefined
  param.eff.mTop    = 175;  
  param.eff.mHiggs  = 120 ;
  param.eff.mh0     = 80; // Used only when !param.eff_use_BRs && param.susy.use_TREE_LEVEL

  //This are not
  param.eff.wTop    = 1.4;  
  param.eff.wHiggs  = 1.4;  
  param.eff_use_BRs = kFALSE  ;  //kFALSE
     param.br.tHb  = 0;
     param.br.Hcs  = 0;
     param.br.HWbb = 0;
     param.br.HWh0 = 0; // Used only when param.eff_use_BRs && param.susy.use_TREE_LEVEL
     param.br.h0bb = 0; 

  param.dataX   = 1.0   ;  //1.0
  param.smear   = kTRUE ;  //kTRUE

  param.susy.use_TREE_LEVEL = kFALSE; //IRRELEVANT
  param.susy.use_CPsuperH  = kTRUE; 
    param.susy.mu  = -500 ;
    param.susy.Mg  = 1000 ;
    param.susy.M2  = 1000 ; 
    param.susy.MQ  = 1000 ; 
    param.susy.MU  = 1000 ; 
    param.susy.MD  = 1000 ; 
    param.susy.Recalculate_stop_masses=kTRUE;
    param.susy.mt1 = 1000 ; 
    param.susy.mt2 = param.susy.mt1 ; 
    param.susy.Recalculate_sbottom_masses=kTRUE;
    param.susy.mb1 = 1000 ; 
    param.susy.mb2 = param.susy.mb1 ;
    param.susy.At = 500;
    param.susy.Ab = 500;
    param.susy.thetat = 0;
    param.susy.thetab = 0;
    param.susy.ht = 0; // always recalculated 

     //The type of study
     param.StudyType = 0 ; // THIS SHOULD BE TAKEN OUT OF HERE
     // STUDY = 0  ;// For regular benchmarks, nothing is redefined 
     // STUDY = 1  ;// MINIMAL MIXING   : Calculation of At on fly (Ab=At)
     // STUDY = 2  ;// MAXIMAL MIXING   : Calculation of At on fly (Ab=At)
     // STUDY = 3  ;// LARGE mu SCENARIO: Calculation of At on fly (Ab=At)

    //Redefine these two things according to the Benchmarks
    if (Bench==1){
       param.StudyType = 0 ;
       param.susy.mu = -500;
       param.susy.At = 2000;
       param.susy.Ab = param.susy.At;

    }else if (Bench==2){
       param.StudyType = 0 ;
       param.susy.mu = -500;
       param.susy.At = -500;
       param.susy.Ab = param.susy.At;

    }else if (Bench==3){
       param.StudyType = 0 ;
       param.susy.mu = 500 ;
       param.susy.At = 500 ; 
       param.susy.Ab = param.susy.At;

    }else if (Bench==4){
       param.StudyType = 0 ;
       param.susy.mu = 500 ;
       param.susy.At = 2800;
       param.susy.Ab = param.susy.At;

    }else if (Bench==5){
      param.StudyType = 1; //Minimal Mixing
      Double_t Msusy = 1000;
      param.susy.mu  = -200 ;
      param.susy.Mg  = 0.8 * Msusy ;
      param.susy.M2  = 200  ; 
      param.susy.MQ  = Msusy ; 
      param.susy.MU  = Msusy ; 
      param.susy.MD  = Msusy ; 
      // param.susy.mt1 and param.susy.mt2  will be recalculated
      // param.susy.mb1  and param.susy.mb2 will be recalculated
      // param.susy.At and param.susy.Ab    will be recalculated
      // param.susy.thetat and thetab       will be recalculated


    }else if (Bench==6){
      param.StudyType = 2; //Maximal Mixing
      Double_t Msusy = 1000;
      param.susy.mu  = -200 ;
      param.susy.Mg  = 0.8 * Msusy ;
      param.susy.M2  = 200  ; 
      param.susy.MQ  = Msusy ; 
      param.susy.MU  = Msusy ; 
      param.susy.MD  = Msusy ; 
      // param.susy.mt1 and param.susy.mt2  will be recalculated
      // param.susy.mb1  and param.susy.mb2 will be recalculated
      // param.susy.At and param.susy.Ab    will be recalculated
      // param.susy.thetat and thetab       will be recalculated

    }else if (Bench==7){
      param.StudyType = 3; // Large mu scenario
      Double_t Msusy = 400  ;
      param.susy.mu  = 1000 ;
      param.susy.Mg  = 200  ;
      param.susy.M2  = 400  ; 
      param.susy.MQ  = Msusy ; 
      param.susy.MU  = Msusy ; 
      param.susy.MD  = Msusy ; 
      // param.susy.mt1 and param.susy.mt2  will be recalculated
      // param.susy.mb1  and param.susy.mb2 will be recalculated
      // param.susy.At and param.susy.Ab    will be recalculated
      // param.susy.thetat and thetab       will be recalculated

    }else if(Bench >=40 && Bench <50){
      //Extremely exotic benchmarks

        if (Bench==40){
          param.susy.tanbeta = 30;
          param.susy.mu      = -500;
        }else if (Bench==41){
          param.susy.tanbeta = 30;
          param.susy.mu      = +500 ;
        }else if (Bench==42){
          param.susy.tanbeta = 50;
          param.susy.mu      = -500;
        }else if (Bench==43){
          param.susy.tanbeta = 50;
          param.susy.mu      = +500;
        }
  
        param.eff.mTop = 175;    
  
        //SUSY params
        param.susy.Mg  = 1000 ;
        param.susy.M2  = 1000 ; 
        param.susy.MQ  = 1000 ; 
        param.susy.MU  = 1000 ; 
        param.susy.MD  = 1000 ; 
        param.susy.Recalculate_stop_masses=kFALSE;// Overriden by StudyType
        param.susy.mt1 = 1000 ; 
        param.susy.mt2 = param.susy.mt1 ; 
        param.susy.Recalculate_sbottom_masses=kFALSE;// Overriden by StudyType
        param.susy.mb1 = 1000 ; 
        param.susy.mb2 = param.susy.mb1 ;
        param.susy.At = 500;
        param.susy.Ab = 0;
  
        //The type of study
        param.StudyType = 0 ; 
  
    }else if (Bench >=50 && Bench<60){
      //To compare to Marcela's paper

      param.susy.Recalculate_stop_masses=kTRUE;
      param.susy.Recalculate_sbottom_masses=kTRUE;

      if (Bench==50){

        //Marcela with mu = -500 
        param.susy.mu  = -500;

      } else if (Bench==51){

        //Marcela with mu = +500 
        param.susy.mu  = +500;

      }

    } else {
      cout<<"Benchmark #"<<Bench<<" not defined!"<<endl;
    }

    return param;

}//GetBenchmark