StandardModel/BosonJetLL.c at master · bin-wu-phys/StandardModel · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
#include <math.h>
#include "BosonJetLL.h"
#include <iostream>
#include <fstream>
#include <sstream>
#include <iomanip>
#include <gsl/gsl_sf_bessel.h>
#include <gsl/gsl_sf_dilog.h>

//#define SCREEN_ONLY

BosonJetLL::BosonJetLL():BosonJetNLL(){
  initLL();
}

BosonJetLL::BosonJetLL(double Ecm, Particle *boson):BosonJetNLL(Ecm, boson){
  initLL();
}

void BosonJetLL::initLL(){
  _ngl=false;
}

double BosonJetLL::bHqg(){
  cout << "We are here..." << endl;
  return -(_cf*(2.0*log(-_muh*_muh/_uhat))+_nc*log(_muh*_muh*_uhat/(_shat*_that)));
}

//filename

void BosonJetLL::getFilename(stringstream &filename, double ymin, double ymax){
  filename << "R" << _R << "pT" << ymin << "." << ymax << "LL";
  if(_css){
    filename << "CSS";
  }
  if(_nonpert){
    filename << "NP";
  }
  filename <<"muh" << _muhV << "mub" << _muxTV << "muJ" << _muJV << "mut" << _mutV;
  filename <<".dat";
}