1 |
ceballos |
1.1 |
#include "MitPhysics/Mods/interface/PwhgWrapper.h"
|
2 |
|
|
|
3 |
|
|
|
4 |
|
|
extern"C" {
|
5 |
|
|
void pwhg_cphto_reweight_(double *mh, double *gh, double *mt, int *BWflag, double *m, double *w);
|
6 |
|
|
}
|
7 |
|
|
|
8 |
|
|
|
9 |
|
|
double pwhegwrapper::getweight(double mh,double gh,double mt,double m,int BWflag){
|
10 |
|
|
/*
|
11 |
|
|
c INPUT
|
12 |
|
|
c mh : Higgs boson mass (used in the POWHEG BOX generation)
|
13 |
|
|
c gh : Higgs boson width (used in the POWHEG BOX generation)
|
14 |
|
|
c mt : top quark mass
|
15 |
|
|
c BWflag : 0 if the sample to reweight was produced with fixed Higgs width
|
16 |
|
|
c 1 if the sample to reweight was produced with running Higgs
|
17 |
|
|
c width (this is the default in the POWHEG BOX)
|
18 |
|
|
c m : virtuality of the produced Higgs boson resonance
|
19 |
|
|
c OUTPUT
|
20 |
|
|
c w : the reweighting factor
|
21 |
|
|
*/
|
22 |
|
|
double w;
|
23 |
|
|
|
24 |
|
|
pwhg_cphto_reweight_(&mh, &gh, &mt, &BWflag, &m, &w);
|
25 |
|
|
|
26 |
|
|
return w;
|
27 |
|
|
}
|