3157 int ZZWWdec_onevertexflag;
3159 outname = Form(
"testME_ProdDec_MCFM_JHUGen_WBFZZ_Comparison_Ping_%i_%i_%i_%s.out", motherflavor, isZZWW, vbfvhchannel, (hasInterf ?
"4l" :
"2l2l"));
3160 ZZWWdec_onevertexflag=2;
3162 else if(decZZWW==2){
3163 outname = Form(
"testME_ProdDec_MCFM_JHUGen_WBFWW_Comparison_Ping_%i_%i_%i.out", motherflavor, isZZWW, vbfvhchannel);
3164 ZZWWdec_onevertexflag=1;
3167 if (hasInterf==1 && decZZWW!=1)
return;
3174 for (
int ii=0; ii<
nmsq; ii++){
for (
int jj=0; jj<
nmsq; jj++) mearray[ii][jj]=0; }
3176 void add(
const mcfmme& other){
3177 proddecme+=other.proddecme;
3178 for (
int ii=0; ii<
nmsq; ii++){
for (
int jj=0; jj<
nmsq; jj++) mearray[ii][jj]+=(other.mearray)[ii][jj]; }
3180 void multiplyarray(
const float val){
3181 for (
int ii=0; ii<
nmsq; ii++){
for (
int jj=0; jj<
nmsq; jj++) mearray[ii][jj]*=val; }
3184 for (
int ii=0; ii<
nmsq; ii++){
for (
int jj=0; jj<
nmsq; jj++)
cout <<
'\t' << mearray[ii][jj];
cout << endl; }
3196 for (
int ii=0; ii<
nmsq; ii++){
for (
int jj=0; jj<
nmsq; jj++) mearray[ii][jj]=0; }
3198 void multiplyarray(
const float val){
3199 for (
int ii=0; ii<
nmsq; ii++){
for (
int jj=0; jj<
nmsq; jj++) mearray[ii][jj]*=val; }
3202 for (
int ii=0; ii<
nmsq; ii++){
for (
int jj=0; jj<
nmsq; jj++)
cout <<
'\t' << mearray[ii][jj];
cout << endl; }
3208 float wPOLE=4.07e-3;
3213 int idMother[2]={ 0 };
3214 if (vbfvhchannel<=1){
3216 if (motherflavor==0){
3219 else if (motherflavor==1){
3220 idMother[0]=2; idMother[1]=1;
3222 (isZZWW==1 && vbfvhchannel==0) ||
3223 (isZZWW==2 && vbfvhchannel==0);
3225 else if (motherflavor==2){
3226 idMother[0]=-2; idMother[1]=-1;
3228 (isZZWW==1 && vbfvhchannel==0) ||
3229 (isZZWW==2 && vbfvhchannel==0);
3232 else if (motherflavor==3){
3233 idMother[0]=2; idMother[1]=-1;
3235 (isZZWW==1 && vbfvhchannel==0) ||
3236 (isZZWW==2 && vbfvhchannel>=1) ||
3237 (isZZWW==0 && vbfvhchannel==0) ||
3238 (isZZWW==0 && vbfvhchannel>=1);
3240 else if (motherflavor==4){
3241 idMother[0]=-2; idMother[1]=1;
3243 (isZZWW==1 && vbfvhchannel==0) ||
3244 (isZZWW==2 && vbfvhchannel>=1) ||
3245 (isZZWW==0 && vbfvhchannel==0) ||
3246 (isZZWW==0 && vbfvhchannel==1);
3249 else if (motherflavor==5){
3250 idMother[0]=2; idMother[1]=-2;
3252 (isZZWW==1 && vbfvhchannel==0) ||
3253 (isZZWW==2 && vbfvhchannel==0) ||
3254 (isZZWW==1 && vbfvhchannel>=1);
3256 else if (motherflavor==6){
3257 idMother[0]=-2; idMother[1]=2;
3259 (isZZWW==1 && vbfvhchannel==0) ||
3260 (isZZWW==2 && vbfvhchannel==0) ||
3261 (isZZWW==1 && vbfvhchannel>=1);
3263 else if (motherflavor==7){
3264 idMother[0]=1; idMother[1]=-1;
3266 (isZZWW==1 && vbfvhchannel==0) ||
3267 (isZZWW==2 && vbfvhchannel==0) ||
3268 (isZZWW==1 && vbfvhchannel>=1);
3270 else if (motherflavor==8){
3271 idMother[0]=-1; idMother[1]=1;
3273 (isZZWW==1 && vbfvhchannel==0) ||
3274 (isZZWW==2 && vbfvhchannel==0) ||
3275 (isZZWW==1 && vbfvhchannel>=1);
3279 else if (motherflavor==9){
3280 idMother[0]=4; idMother[1]=3;
3282 (isZZWW==1 && vbfvhchannel==0) ||
3283 (isZZWW==2 && vbfvhchannel==0);
3285 else if (motherflavor==10){
3286 idMother[0]=4; idMother[1]=1;
3288 (isZZWW==1 && vbfvhchannel==0) ||
3289 (isZZWW==2 && vbfvhchannel==0);
3291 else if (motherflavor==11){
3292 idMother[0]=2; idMother[1]=3;
3294 (isZZWW==1 && vbfvhchannel==0) ||
3295 (isZZWW==2 && vbfvhchannel==0);
3297 else if (motherflavor==12){
3298 idMother[0]=-4; idMother[1]=-3;
3300 (isZZWW==1 && vbfvhchannel==0) ||
3301 (isZZWW==2 && vbfvhchannel==0);
3303 else if (motherflavor==13){
3304 idMother[0]=-4; idMother[1]=-1;
3306 (isZZWW==1 && vbfvhchannel==0) ||
3307 (isZZWW==2 && vbfvhchannel==0);
3309 else if (motherflavor==14){
3310 idMother[0]=-2; idMother[1]=-3;
3312 (isZZWW==1 && vbfvhchannel==0) ||
3313 (isZZWW==2 && vbfvhchannel==0);
3316 else if (motherflavor==15){
3317 idMother[0]=4; idMother[1]=-3;
3319 (isZZWW==1 && vbfvhchannel==0) ||
3320 (isZZWW==2 && vbfvhchannel>=1) ||
3321 (isZZWW==0 && vbfvhchannel==0) ||
3322 (isZZWW==0 && vbfvhchannel>=1);
3324 else if (motherflavor==16){
3325 idMother[0]=4; idMother[1]=-1;
3327 (isZZWW==1 && vbfvhchannel==0) ||
3328 (isZZWW==2 && vbfvhchannel>=1) ||
3329 (isZZWW==0 && vbfvhchannel==0) ||
3330 (isZZWW==0 && vbfvhchannel>=1);
3332 else if (motherflavor==17){
3333 idMother[0]=2; idMother[1]=-3;
3335 (isZZWW==1 && vbfvhchannel==0) ||
3336 (isZZWW==2 && vbfvhchannel>=1) ||
3337 (isZZWW==0 && vbfvhchannel==0) ||
3338 (isZZWW==0 && vbfvhchannel>=1);
3340 else if (motherflavor==18){
3341 idMother[0]=-4; idMother[1]=3;
3343 (isZZWW==1 && vbfvhchannel==0) ||
3344 (isZZWW==2 && vbfvhchannel>=1) ||
3345 (isZZWW==0 && vbfvhchannel==0) ||
3346 (isZZWW==0 && vbfvhchannel==1);
3348 else if (motherflavor==19){
3349 idMother[0]=-4; idMother[1]=1;
3351 (isZZWW==1 && vbfvhchannel==0) ||
3352 (isZZWW==2 && vbfvhchannel>=1) ||
3353 (isZZWW==0 && vbfvhchannel==0) ||
3354 (isZZWW==0 && vbfvhchannel==1);
3356 else if (motherflavor==20){
3357 idMother[0]=-2; idMother[1]=3;
3359 (isZZWW==1 && vbfvhchannel==0) ||
3360 (isZZWW==2 && vbfvhchannel>=1) ||
3361 (isZZWW==0 && vbfvhchannel==0) ||
3362 (isZZWW==0 && vbfvhchannel==1);
3365 else if (motherflavor==21){
3366 idMother[0]=4; idMother[1]=-4;
3368 (isZZWW==1 && vbfvhchannel==0) ||
3369 (isZZWW==2 && vbfvhchannel==0) ||
3370 (isZZWW==1 && vbfvhchannel>=1);
3372 else if (motherflavor==22){
3373 idMother[0]=-4; idMother[1]=4;
3375 (isZZWW==1 && vbfvhchannel==0) ||
3376 (isZZWW==2 && vbfvhchannel==0) ||
3377 (isZZWW==1 && vbfvhchannel>=1);
3379 else if (motherflavor==23){
3380 idMother[0]=3; idMother[1]=-3;
3382 (isZZWW==1 && vbfvhchannel==0) ||
3383 (isZZWW==2 && vbfvhchannel==0) ||
3384 (isZZWW==1 && vbfvhchannel>=1);
3386 else if (motherflavor==24){
3387 idMother[0]=-3; idMother[1]=3;
3389 (isZZWW==1 && vbfvhchannel==0) ||
3390 (isZZWW==2 && vbfvhchannel==0) ||
3391 (isZZWW==1 && vbfvhchannel>=1);
3393 else if (motherflavor==25){
3394 idMother[0]=5; idMother[1]=-5;
3396 (isZZWW==1 && vbfvhchannel==0) ||
3397 (isZZWW==2 && vbfvhchannel==0) ||
3398 (isZZWW==1 && vbfvhchannel>=1);
3400 else if (motherflavor==26){
3401 idMother[0]=-5; idMother[1]=5;
3403 (isZZWW==1 && vbfvhchannel==0) ||
3404 (isZZWW==2 && vbfvhchannel==0) ||
3405 (isZZWW==1 && vbfvhchannel>=1);
3407 else if (motherflavor<=31){
3408 idMother[0]=-(motherflavor-26); idMother[1]=0;
3410 (isZZWW==1 && vbfvhchannel==0) ||
3411 (isZZWW==2 && vbfvhchannel==0) ||
3412 (isZZWW==1 && vbfvhchannel>=1) ||
3413 (isZZWW==2 && vbfvhchannel>=1);
3415 else if (motherflavor<=36){
3416 idMother[0]=(motherflavor-31); idMother[1]=0;
3418 (isZZWW==1 && vbfvhchannel==0) ||
3419 (isZZWW==2 && vbfvhchannel==0) ||
3420 (isZZWW==1 && vbfvhchannel>=1) ||
3421 (isZZWW==2 && vbfvhchannel>=1);
3424 else if (vbfvhchannel==2){
3426 if (motherflavor==3){
3427 idMother[0]=14; idMother[1]=-13;
3429 (isZZWW==2 && vbfvhchannel>=1) ||
3430 (isZZWW==0 && vbfvhchannel>=1);
3432 else if (motherflavor==4){
3433 idMother[0]=-14; idMother[1]=13;
3435 (isZZWW==2 && vbfvhchannel>=1) ||
3436 (isZZWW==0 && vbfvhchannel>=1);
3439 else if (motherflavor==5){
3440 idMother[0]=14; idMother[1]=-14;
3442 (isZZWW==1 && vbfvhchannel>=1);
3444 else if (motherflavor==6){
3445 idMother[0]=-14; idMother[1]=14;
3447 (isZZWW==1 && vbfvhchannel>=1);
3449 else if (motherflavor==7){
3450 idMother[0]=13; idMother[1]=-13;
3452 (isZZWW==1 && vbfvhchannel>=1);
3454 else if (motherflavor==8){
3455 idMother[0]=-13; idMother[1]=13;
3457 (isZZWW==1 && vbfvhchannel>=1);
3463 else if (vbfvhchannel==1){
3465 (idMother[0]==-idMother[1] && idMother[0]!=0)
3467 (idMother[0]==-idMother[1] && idMother[0]==0 && isZZWW==1)
3469 (idMother[0]!=-idMother[1] && (idMother[0]==0 || idMother[1]==0) && isZZWW==1)
3472 (TMath::Sign(1, idMother[0])==-TMath::Sign(1, idMother[1]) && abs(idMother[0])%2!=abs(idMother[1])%2)
3474 (idMother[0]==-idMother[1] && idMother[0]==0 && isZZWW==2)
3476 ((idMother[0]==0 || idMother[1]==0) && isZZWW==2)
3480 else if (vbfvhchannel==2){
3482 else if (TMath::Sign(1, idMother[0])==-TMath::Sign(1, idMother[1]) && abs(idMother[0])%2!=abs(idMother[1])%2)
prod=
TVar::Lep_WH_S;
3488 ofstream tout(outname.Data());
3489 streambuf* coutbuf =
cout.rdbuf();
3490 cout.rdbuf(tout.rdbuf());
3493 melaptr.reset(
new Mela(erg_tev, mPOLE, verbosity));
3497 mela.setVerbosity(verbosity);
3499 jhume p_prod_0mplus_dec_0mplus_VAJHU;
3500 jhume p_prod_0hplus_dec_0hplus_VAJHU;
3501 jhume p_prod_fa2_dec_fa2_VAJHU;
3502 jhume p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU;
3503 jhume p_prod_0minus_dec_0minus_VAJHU;
3504 jhume p_prod_fa3_dec_fa3_VAJHU;
3505 jhume p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU;
3506 jhume p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU;
3507 jhume p_prod_0mplusL2_dec_0mplusL2_VAJHU;
3508 jhume p_prod_fL2_dec_fL2_VAJHU;
3509 jhume p_prod_0minusZA_dec_0minusZA_VAJHU;
3510 jhume p_prod_fa3ZA_dec_fa3ZA_VAJHU;
3511 jhume p_prod_0minusAA_dec_0minusAA_VAJHU;
3512 jhume p_prod_fa3AA_dec_fa3AA_VAJHU;
3514 mcfmme p_bkg_VAMCFM, p_bkg_VAMCFM_rssum;
3515 mcfmme p_prod_0mplus_dec_0mplus_VAMCFM;
3516 mcfmme p_prod_0hplus_dec_0hplus_VAMCFM;
3517 mcfmme p_prod_0minus_dec_0minus_VAMCFM;
3518 mcfmme p_prod_0mplusL2_dec_0mplusL2_VAMCFM;
3519 mcfmme p_prod_fL2_dec_fL2_VAMCFM;
3520 mcfmme p_prod_fa2_dec_fa2_VAMCFM;
3521 mcfmme p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAMCFM;
3522 mcfmme p_prod_fa3_dec_fa3_VAMCFM;
3523 mcfmme p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAMCFM;
3524 mcfmme p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAMCFM;
3525 mcfmme p_prod_0minusZA_dec_0minusZA_VAMCFM;
3526 mcfmme p_prod_fa3ZA_dec_fa3ZA_VAMCFM;
3527 mcfmme p_prod_0minusAA_dec_0minusAA_VAMCFM;
3528 mcfmme p_prod_fa3AA_dec_fa3AA_VAMCFM;
3535 float costhetastar=0;
3542 float mVstar=0, mV=0;
3544 float pingMom[8][4]={
3545 { 0, 0, 865.37881546721542, 865.37881546721542 },
3546 { 0, 0, -624.03396598421773, 624.03396598421773 },
3547 { 7.6145299215002638, -17.259247740062808, 9.4660586470659975, 21.106135714241464 },
3548 { 90.901719112641416, -69.683681833050798, 32.066319224729980, 118.94194752090492 },
3549 { 78.476352131782917, -35.264818847819797, -8.8615639484695272, 86.490881645951262 },
3550 { 191.68369742375290, -197.85205601463366, 100.99437243828194, 293.40746273989180 },
3551 { -131.59521398083137, 330.56000090294270, 437.01695094737875, 563.53440884737279 },
3552 { -237.08108460884614, -10.500196467375645, -329.33728782598945, 405.93194498307093 }
3555 int GenLep1Id=0, GenLep2Id=0, GenLep3Id=0, GenLep4Id=0;
3564 GenLep3Id=GenLep1Id;
3565 GenLep4Id=GenLep2Id;
3576 int idOrdered[4] ={ GenLep1Id, GenLep2Id, GenLep3Id, GenLep4Id };
3579 for (
unsigned int ip=0; ip<2; ip++){
3583 TLorentzVector(pingMom[ip][0], pingMom[ip][1], pingMom[ip][2], pingMom[ip][3])
3588 for (
unsigned int ip=2; ip<6; ip++){
3592 TLorentzVector(pingMom[ip][0], pingMom[ip][1], pingMom[ip][2], pingMom[ip][3])
3597 for (
unsigned int ip=6; ip<8; ip++){
3601 TLorentzVector(pingMom[ip][0], pingMom[ip][1], pingMom[ip][2], pingMom[ip][3])
3700 double invckm_ud=1, invckm_us=0, invckm_cd=0, invckm_cs=1, invckm_ts=0, invckm_tb=1, invckm_ub=0, invckm_cb=0, invckm_td=0;
3701 TUtil::SetCKMElements(&invckm_ud, &invckm_us, &invckm_cd, &invckm_cs, &invckm_ts, &invckm_tb, &invckm_ub, &invckm_cb, &invckm_td);
3707 mela.computeP(p_prod_0mplus_dec_0mplus_VAJHU.decme,
false);
3710 mela.computeP(p_prod_0hplus_dec_0hplus_VAJHU.decme,
false);
3714 mela.computeP(p_prod_fa2_dec_fa2_VAJHU.decme,
false);
3718 mela.computeP(p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.decme,
false);
3721 mela.computeP(p_prod_0minus_dec_0minus_VAJHU.decme,
false);
3725 mela.computeP(p_prod_fa3_dec_fa3_VAJHU.decme,
false);
3729 mela.computeP(p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.decme,
false);
3733 mela.computeP(p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.decme,
false);
3737 mela.computeP(p_prod_0minusZA_dec_0minusZA_VAJHU.decme,
false);
3740 mela.computeP(p_prod_0minusAA_dec_0minusAA_VAJHU.decme,
false);
3744 mela.computeP(p_prod_fa3ZA_dec_fa3ZA_VAJHU.decme,
false);
3748 mela.computeP(p_prod_fa3AA_dec_fa3AA_VAJHU.decme,
false);
3750 p_prod_0mplusL2_dec_0mplusL2_VAJHU.decme=p_prod_0mplus_dec_0mplus_VAJHU.decme;
3751 p_prod_fL2_dec_fL2_VAJHU.decme=p_prod_0mplus_dec_0mplus_VAJHU.decme;
3755 mela.computeP(p_prod_0mplusL2_dec_0mplusL2_VAJHU.decme,
false);
3759 mela.computeP(p_prod_fL2_dec_fL2_VAJHU.decme,
false);
3761 p_prod_0minusZA_dec_0minusZA_VAJHU.decme=p_prod_0mplus_dec_0mplus_VAJHU.decme;
3762 p_prod_0minusAA_dec_0minusAA_VAJHU.decme=p_prod_0mplus_dec_0mplus_VAJHU.decme;
3763 p_prod_fa3ZA_dec_fa3ZA_VAJHU.decme=p_prod_0mplus_dec_0mplus_VAJHU.decme;
3764 p_prod_fa3AA_dec_fa3AA_VAJHU.decme=p_prod_0mplus_dec_0mplus_VAJHU.decme;
3767 bool computeL2WWprod=(isZZWW==2 || (decZZWW==2 && vbfvhchannel==0));
3768 bool computeJHUZA=
true;
3769 if (vbfvhchannel>=1){
3770 if (isZZWW==2 && decZZWW==1){
3772 else if (isZZWW==1 && decZZWW==2){
3774 else if (isZZWW==1 && decZZWW==1){
3776 else computeJHUZA=
false;
3779 if (isZZWW==2 && decZZWW==1){
3781 else if (isZZWW==1 && decZZWW==2){
3783 else if (isZZWW==1 && decZZWW==1){
3785 else computeJHUZA=
false;
3787 if (vbfvhchannel==0){
3791 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_VV_1][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
3793 mela.computeProdP(p_prod_0mplus_dec_0mplus_VAJHU.prodme,
false);
3794 mela.getIORecord()->getUnweightedMEArray(p_prod_0mplus_dec_0mplus_VAJHU.mearray);
3796 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_VV_2][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
3798 mela.computeProdP(p_prod_0hplus_dec_0hplus_VAJHU.prodme,
false);
3799 mela.getIORecord()->getUnweightedMEArray(p_prod_0hplus_dec_0hplus_VAJHU.mearray);
3803 mela.computeProdP(p_prod_fa2_dec_fa2_VAJHU.prodme,
false);
3804 mela.getIORecord()->getUnweightedMEArray(p_prod_fa2_dec_fa2_VAJHU.mearray);
3808 mela.computeProdP(p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.prodme,
false);
3809 mela.getIORecord()->getUnweightedMEArray(p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.mearray);
3811 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_VV_4][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
3813 mela.computeProdP(p_prod_0minus_dec_0minus_VAJHU.prodme,
false);
3814 mela.getIORecord()->getUnweightedMEArray(p_prod_0minus_dec_0minus_VAJHU.mearray);
3818 mela.computeProdP(p_prod_fa3_dec_fa3_VAJHU.prodme,
false);
3819 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3_dec_fa3_VAJHU.mearray);
3823 mela.computeProdP(p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.prodme,
false);
3824 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.mearray);
3828 mela.computeProdP(p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.prodme,
false);
3829 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.mearray);
3831 if (computeL2WWprod){
3834 mela.computeProdP(p_prod_0mplusL2_dec_0mplusL2_VAJHU.prodme,
false);
3835 mela.getIORecord()->getUnweightedMEArray(p_prod_0mplusL2_dec_0mplusL2_VAJHU.mearray);
3839 mela.computeProdP(p_prod_fL2_dec_fL2_VAJHU.prodme,
false);
3840 mela.getIORecord()->getUnweightedMEArray(p_prod_fL2_dec_fL2_VAJHU.mearray);
3843 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_VV_1][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
3845 mela.computeProdP(p_prod_0mplusL2_dec_0mplusL2_VAJHU.prodme,
false);
3846 mela.getIORecord()->getUnweightedMEArray(p_prod_0mplusL2_dec_0mplusL2_VAJHU.mearray);
3848 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_VV_1][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
3850 mela.computeProdP(p_prod_fL2_dec_fL2_VAJHU.prodme,
false);
3851 mela.getIORecord()->getUnweightedMEArray(p_prod_fL2_dec_fL2_VAJHU.mearray);
3854 if (isZZWW==2 && decZZWW==1){
3857 mela.differentiate_HWW_HZZ=
true;
3859 else if (isZZWW==1 && decZZWW==2){
3862 mela.differentiate_HWW_HZZ=
true;
3864 else if (isZZWW==1 && decZZWW==1){
3866 mela.differentiate_HWW_HZZ=
true;
3868 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_ZA_4][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
3870 mela.computeProdP(p_prod_0minusZA_dec_0minusZA_VAJHU.prodme,
false);
3871 mela.getIORecord()->getUnweightedMEArray(p_prod_0minusZA_dec_0minusZA_VAJHU.mearray);
3873 if (isZZWW==2 && decZZWW==1){
3876 mela.differentiate_HWW_HZZ=
true;
3878 else if (isZZWW==1 && decZZWW==2){
3881 mela.differentiate_HWW_HZZ=
true;
3883 else if (isZZWW==1 && decZZWW==1){
3885 mela.differentiate_HWW_HZZ=
true;
3887 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_AA_4][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
3889 mela.computeProdP(p_prod_0minusAA_dec_0minusAA_VAJHU.prodme,
false);
3890 mela.getIORecord()->getUnweightedMEArray(p_prod_0minusAA_dec_0minusAA_VAJHU.mearray);
3892 if (isZZWW==2 && decZZWW==1){
3896 mela.differentiate_HWW_HZZ=
true;
3898 else if (isZZWW==1 && decZZWW==2){
3902 mela.differentiate_HWW_HZZ=
true;
3904 else if (isZZWW==1 && decZZWW==1){
3907 mela.differentiate_HWW_HZZ=
true;
3911 mela.computeProdP(p_prod_fa3ZA_dec_fa3ZA_VAJHU.prodme,
false);
3912 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3ZA_dec_fa3ZA_VAJHU.mearray);
3914 if (isZZWW==2 && decZZWW==1){
3918 mela.differentiate_HWW_HZZ=
true;
3920 else if (isZZWW==1 && decZZWW==2){
3924 mela.differentiate_HWW_HZZ=
true;
3926 else if (isZZWW==1 && decZZWW==1){
3929 mela.differentiate_HWW_HZZ=
true;
3933 mela.computeProdP(p_prod_fa3AA_dec_fa3AA_VAJHU.prodme,
false);
3934 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3AA_dec_fa3AA_VAJHU.mearray);
3937 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_VV_1][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
3939 mela.computeProdP(p_prod_0minusZA_dec_0minusZA_VAJHU.prodme,
false);
3940 mela.getIORecord()->getUnweightedMEArray(p_prod_0minusZA_dec_0minusZA_VAJHU.mearray);
3942 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_VV_1][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
3944 mela.computeProdP(p_prod_0minusAA_dec_0minusAA_VAJHU.prodme,
false);
3945 mela.getIORecord()->getUnweightedMEArray(p_prod_0minusAA_dec_0minusAA_VAJHU.mearray);
3947 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_VV_1][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
3949 mela.computeProdP(p_prod_fa3ZA_dec_fa3ZA_VAJHU.prodme,
false);
3950 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3ZA_dec_fa3ZA_VAJHU.mearray);
3952 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_VV_1][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
3954 mela.computeProdP(p_prod_fa3AA_dec_fa3AA_VAJHU.prodme,
false);
3955 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3AA_dec_fa3AA_VAJHU.mearray);
3971 mela.computeProdP(p_prod_0mplus_dec_0mplus_VAJHU.prodme,
false);
3972 mela.getIORecord()->getUnweightedMEArray(p_prod_0mplus_dec_0mplus_VAJHU.mearray);
3976 mela.computeProdP(p_prod_0hplus_dec_0hplus_VAJHU.prodme,
false);
3977 mela.getIORecord()->getUnweightedMEArray(p_prod_0hplus_dec_0hplus_VAJHU.mearray);
3981 mela.computeProdP(p_prod_fa2_dec_fa2_VAJHU.prodme,
false);
3982 mela.getIORecord()->getUnweightedMEArray(p_prod_fa2_dec_fa2_VAJHU.mearray);
3986 mela.computeProdP(p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.prodme,
false);
3987 mela.getIORecord()->getUnweightedMEArray(p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.mearray);
3991 mela.computeProdP(p_prod_0minus_dec_0minus_VAJHU.prodme,
false);
3992 mela.getIORecord()->getUnweightedMEArray(p_prod_0minus_dec_0minus_VAJHU.mearray);
3996 mela.computeProdP(p_prod_fa3_dec_fa3_VAJHU.prodme,
false);
3997 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3_dec_fa3_VAJHU.mearray);
4001 mela.computeProdP(p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.prodme,
false);
4002 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.mearray);
4006 mela.computeProdP(p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.prodme,
false);
4007 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.mearray);
4011 mela.computeProdP(p_prod_0mplusL2_dec_0mplusL2_VAJHU.prodme,
false);
4012 mela.getIORecord()->getUnweightedMEArray(p_prod_0mplusL2_dec_0mplusL2_VAJHU.mearray);
4015 mela.computeProdP(p_prod_fL2_dec_fL2_VAJHU.prodme,
false);
4016 mela.getIORecord()->getUnweightedMEArray(p_prod_fL2_dec_fL2_VAJHU.mearray);
4020 mela.computeProdP(p_prod_0minusZA_dec_0minusZA_VAJHU.prodme,
false);
4021 mela.getIORecord()->getUnweightedMEArray(p_prod_0minusZA_dec_0minusZA_VAJHU.mearray);
4024 mela.computeProdP(p_prod_0minusAA_dec_0minusAA_VAJHU.prodme,
false);
4025 mela.getIORecord()->getUnweightedMEArray(p_prod_0minusAA_dec_0minusAA_VAJHU.mearray);
4028 mela.computeProdP(p_prod_fa3ZA_dec_fa3ZA_VAJHU.prodme,
false);
4029 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3ZA_dec_fa3ZA_VAJHU.mearray);
4032 mela.computeProdP(p_prod_fa3AA_dec_fa3AA_VAJHU.prodme,
false);
4033 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3AA_dec_fa3AA_VAJHU.mearray);
4037 mela.computeProdP(p_prod_0minusZA_dec_0minusZA_VAJHU.prodme,
false);
4038 mela.getIORecord()->getUnweightedMEArray(p_prod_0minusZA_dec_0minusZA_VAJHU.mearray);
4041 mela.computeProdP(p_prod_0minusAA_dec_0minusAA_VAJHU.prodme,
false);
4042 mela.getIORecord()->getUnweightedMEArray(p_prod_0minusAA_dec_0minusAA_VAJHU.mearray);
4045 mela.computeProdP(p_prod_fa3ZA_dec_fa3ZA_VAJHU.prodme,
false);
4046 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3ZA_dec_fa3ZA_VAJHU.mearray);
4049 mela.computeProdP(p_prod_fa3AA_dec_fa3AA_VAJHU.prodme,
false);
4050 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3AA_dec_fa3AA_VAJHU.mearray);
4064 mela.computeProdP(p_prod_0mplus_dec_0mplus_VAJHU.prodme,
false);
4065 mela.getIORecord()->getUnweightedMEArray(p_prod_0mplus_dec_0mplus_VAJHU.mearray);
4068 mela.computeProdP(p_prod_0hplus_dec_0hplus_VAJHU.prodme,
false);
4069 mela.getIORecord()->getUnweightedMEArray(p_prod_0hplus_dec_0hplus_VAJHU.mearray);
4073 mela.computeProdP(p_prod_fa2_dec_fa2_VAJHU.prodme,
false);
4074 mela.getIORecord()->getUnweightedMEArray(p_prod_fa2_dec_fa2_VAJHU.mearray);
4078 mela.computeProdP(p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.prodme,
false);
4079 mela.getIORecord()->getUnweightedMEArray(p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.mearray);
4082 mela.computeProdP(p_prod_0minus_dec_0minus_VAJHU.prodme,
false);
4083 mela.getIORecord()->getUnweightedMEArray(p_prod_0minus_dec_0minus_VAJHU.mearray);
4087 mela.computeProdP(p_prod_fa3_dec_fa3_VAJHU.prodme,
false);
4088 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3_dec_fa3_VAJHU.mearray);
4092 mela.computeProdP(p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.prodme,
false);
4093 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.mearray);
4097 mela.computeProdP(p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.prodme,
false);
4098 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.mearray);
4101 mela.computeProdP(p_prod_0minusZA_dec_0minusZA_VAJHU.prodme,
false);
4102 mela.getIORecord()->getUnweightedMEArray(p_prod_0minusZA_dec_0minusZA_VAJHU.mearray);
4105 mela.computeProdP(p_prod_0minusAA_dec_0minusAA_VAJHU.prodme,
false);
4106 mela.getIORecord()->getUnweightedMEArray(p_prod_0minusAA_dec_0minusAA_VAJHU.mearray);
4109 mela.computeProdP(p_prod_fa3ZA_dec_fa3ZA_VAJHU.prodme,
false);
4110 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3ZA_dec_fa3ZA_VAJHU.mearray);
4113 mela.computeProdP(p_prod_fa3AA_dec_fa3AA_VAJHU.prodme,
false);
4114 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3AA_dec_fa3AA_VAJHU.mearray);
4118 mela.computeProdP(p_prod_0mplusL2_dec_0mplusL2_VAJHU.prodme,
false);
4119 mela.getIORecord()->getUnweightedMEArray(p_prod_0mplusL2_dec_0mplusL2_VAJHU.mearray);
4122 mela.computeProdP(p_prod_fL2_dec_fL2_VAJHU.prodme,
false);
4123 mela.getIORecord()->getUnweightedMEArray(p_prod_fL2_dec_fL2_VAJHU.mearray);
4130 double propagator = 1./(pow(pow(mzz, 2)-pow(mPOLE, 2), 2)+pow(mPOLE*wPOLE, 2));
4131 if (vbfvhchannel>=1){
4132 mela.getIORecord()->getHiggsMassWidth(mh, gah, 0);
4133 propagator /= 1./(pow(pow(mzz, 2)-pow(mh, 2), 2) + pow(mh*gah, 2));
4135 p_prod_0mplus_dec_0mplus_VAJHU.prodme *=
propagator; p_prod_0mplus_dec_0mplus_VAJHU.multiplyarray(
propagator);
4136 p_prod_0minus_dec_0minus_VAJHU.prodme *=
propagator; p_prod_0minus_dec_0minus_VAJHU.multiplyarray(
propagator);
4137 p_prod_0hplus_dec_0hplus_VAJHU.prodme *=
propagator; p_prod_0hplus_dec_0hplus_VAJHU.multiplyarray(
propagator);
4138 p_prod_0mplusL2_dec_0mplusL2_VAJHU.prodme *=
propagator; p_prod_0mplusL2_dec_0mplusL2_VAJHU.multiplyarray(
propagator);
4139 p_prod_0minusZA_dec_0minusZA_VAJHU.prodme *=
propagator; p_prod_0minusZA_dec_0minusZA_VAJHU.multiplyarray(
propagator);
4140 p_prod_0minusAA_dec_0minusAA_VAJHU.prodme *=
propagator; p_prod_0minusAA_dec_0minusAA_VAJHU.multiplyarray(
propagator);
4141 p_prod_fa2_dec_fa2_VAJHU.prodme *=
propagator; p_prod_fa2_dec_fa2_VAJHU.multiplyarray(
propagator);
4142 p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.prodme *=
propagator; p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.multiplyarray(
propagator);
4143 p_prod_fa3_dec_fa3_VAJHU.prodme *=
propagator; p_prod_fa3_dec_fa3_VAJHU.multiplyarray(
propagator);
4144 p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.prodme *=
propagator; p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.multiplyarray(
propagator);
4145 p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.prodme *=
propagator; p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.multiplyarray(
propagator);
4146 p_prod_fL2_dec_fL2_VAJHU.prodme *=
propagator; p_prod_fL2_dec_fL2_VAJHU.multiplyarray(
propagator);
4147 p_prod_fa3ZA_dec_fa3ZA_VAJHU.prodme *=
propagator; p_prod_fa3ZA_dec_fa3ZA_VAJHU.multiplyarray(
propagator);
4148 p_prod_fa3AA_dec_fa3AA_VAJHU.prodme *=
propagator; p_prod_fa3AA_dec_fa3AA_VAJHU.multiplyarray(
propagator);
4150 p_prod_0mplus_dec_0mplus_VAJHU.proddecme = p_prod_0mplus_dec_0mplus_VAJHU.prodme*p_prod_0mplus_dec_0mplus_VAJHU.decme; p_prod_0mplus_dec_0mplus_VAJHU.multiplyarray(p_prod_0mplus_dec_0mplus_VAJHU.decme);
4151 if (isZZWW==ZZWWdec_onevertexflag){
4152 p_prod_0hplus_dec_0hplus_VAJHU.proddecme = p_prod_0hplus_dec_0hplus_VAJHU.prodme*p_prod_0mplus_dec_0mplus_VAJHU.decme; p_prod_0hplus_dec_0hplus_VAJHU.multiplyarray(p_prod_0mplus_dec_0mplus_VAJHU.decme);
4153 p_prod_fa2_dec_fa2_VAJHU.proddecme = p_prod_fa2_dec_fa2_VAJHU.prodme*p_prod_0mplus_dec_0mplus_VAJHU.decme; p_prod_fa2_dec_fa2_VAJHU.multiplyarray(p_prod_0mplus_dec_0mplus_VAJHU.decme);
4154 p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.proddecme = p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.prodme*p_prod_0mplus_dec_0mplus_VAJHU.decme; p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.multiplyarray(p_prod_0mplus_dec_0mplus_VAJHU.decme);
4155 p_prod_0minus_dec_0minus_VAJHU.proddecme = p_prod_0minus_dec_0minus_VAJHU.prodme*p_prod_0mplus_dec_0mplus_VAJHU.decme; p_prod_0minus_dec_0minus_VAJHU.multiplyarray(p_prod_0mplus_dec_0mplus_VAJHU.decme);
4156 p_prod_fa3_dec_fa3_VAJHU.proddecme = p_prod_fa3_dec_fa3_VAJHU.prodme*p_prod_0mplus_dec_0mplus_VAJHU.decme; p_prod_fa3_dec_fa3_VAJHU.multiplyarray(p_prod_0mplus_dec_0mplus_VAJHU.decme);
4157 p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.proddecme = p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.prodme*p_prod_0mplus_dec_0mplus_VAJHU.decme; p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.multiplyarray(p_prod_0mplus_dec_0mplus_VAJHU.decme);
4158 p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.proddecme = p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.prodme*p_prod_0mplus_dec_0mplus_VAJHU.decme; p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.multiplyarray(p_prod_0mplus_dec_0mplus_VAJHU.decme);
4161 p_prod_0hplus_dec_0hplus_VAJHU.proddecme = p_prod_0hplus_dec_0hplus_VAJHU.prodme*p_prod_0hplus_dec_0hplus_VAJHU.decme; p_prod_0hplus_dec_0hplus_VAJHU.multiplyarray(p_prod_0hplus_dec_0hplus_VAJHU.decme);
4162 p_prod_fa2_dec_fa2_VAJHU.proddecme = p_prod_fa2_dec_fa2_VAJHU.prodme*p_prod_fa2_dec_fa2_VAJHU.decme; p_prod_fa2_dec_fa2_VAJHU.multiplyarray(p_prod_fa2_dec_fa2_VAJHU.decme);
4163 p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.proddecme = p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.prodme*p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.decme; p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.multiplyarray(p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.decme);
4164 p_prod_0minus_dec_0minus_VAJHU.proddecme = p_prod_0minus_dec_0minus_VAJHU.prodme*p_prod_0minus_dec_0minus_VAJHU.decme; p_prod_0minus_dec_0minus_VAJHU.multiplyarray(p_prod_0minus_dec_0minus_VAJHU.decme);
4165 p_prod_fa3_dec_fa3_VAJHU.proddecme = p_prod_fa3_dec_fa3_VAJHU.prodme*p_prod_fa3_dec_fa3_VAJHU.decme; p_prod_fa3_dec_fa3_VAJHU.multiplyarray(p_prod_fa3_dec_fa3_VAJHU.decme);
4166 p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.proddecme = p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.prodme*p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.decme; p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.multiplyarray(p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.decme);
4167 p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.proddecme = p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.prodme*p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.decme; p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.multiplyarray(p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.decme);
4170 p_prod_0mplusL2_dec_0mplusL2_VAJHU.proddecme = p_prod_0mplusL2_dec_0mplusL2_VAJHU.prodme*p_prod_0mplusL2_dec_0mplusL2_VAJHU.decme; p_prod_0mplusL2_dec_0mplusL2_VAJHU.multiplyarray(p_prod_0mplusL2_dec_0mplusL2_VAJHU.decme);
4171 p_prod_fL2_dec_fL2_VAJHU.proddecme = p_prod_fL2_dec_fL2_VAJHU.prodme*p_prod_fL2_dec_fL2_VAJHU.decme; p_prod_fL2_dec_fL2_VAJHU.multiplyarray(p_prod_fL2_dec_fL2_VAJHU.decme);
4173 p_prod_0minusZA_dec_0minusZA_VAJHU.proddecme = p_prod_0minusZA_dec_0minusZA_VAJHU.prodme*p_prod_0minusZA_dec_0minusZA_VAJHU.decme; p_prod_0minusZA_dec_0minusZA_VAJHU.multiplyarray(p_prod_0minusZA_dec_0minusZA_VAJHU.decme);
4174 p_prod_0minusAA_dec_0minusAA_VAJHU.proddecme = p_prod_0minusAA_dec_0minusAA_VAJHU.prodme*p_prod_0minusAA_dec_0minusAA_VAJHU.decme; p_prod_0minusAA_dec_0minusAA_VAJHU.multiplyarray(p_prod_0minusAA_dec_0minusAA_VAJHU.decme);
4175 p_prod_fa3ZA_dec_fa3ZA_VAJHU.proddecme = p_prod_fa3ZA_dec_fa3ZA_VAJHU.prodme*p_prod_fa3ZA_dec_fa3ZA_VAJHU.decme; p_prod_fa3ZA_dec_fa3ZA_VAJHU.multiplyarray(p_prod_fa3ZA_dec_fa3ZA_VAJHU.decme);
4176 p_prod_fa3AA_dec_fa3AA_VAJHU.proddecme = p_prod_fa3AA_dec_fa3AA_VAJHU.prodme*p_prod_fa3AA_dec_fa3AA_VAJHU.decme; p_prod_fa3AA_dec_fa3AA_VAJHU.multiplyarray(p_prod_fa3AA_dec_fa3AA_VAJHU.decme);
4190 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_VV_1][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
4192 mela.computeProdDecP(p_prod_0mplus_dec_0mplus_VAMCFM.proddecme,
false);
4193 mela.getIORecord()->getUnweightedMEArray(p_prod_0mplus_dec_0mplus_VAMCFM.mearray);
4195 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_VV_2][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
4197 mela.computeProdDecP(p_prod_0hplus_dec_0hplus_VAMCFM.proddecme,
false);
4198 mela.getIORecord()->getUnweightedMEArray(p_prod_0hplus_dec_0hplus_VAMCFM.mearray);
4202 mela.computeProdDecP(p_prod_fa2_dec_fa2_VAMCFM.proddecme,
false);
4203 mela.getIORecord()->getUnweightedMEArray(p_prod_fa2_dec_fa2_VAMCFM.mearray);
4207 mela.computeProdDecP(p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAMCFM.proddecme,
false);
4208 mela.getIORecord()->getUnweightedMEArray(p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAMCFM.mearray);
4210 if (isZZWW!=2){
mela.selfDHzzcoupl[0][
gHIGGS_VV_4][0]=1;
if (isZZWW==1)
mela.differentiate_HWW_HZZ=
true; }
4212 mela.computeProdDecP(p_prod_0minus_dec_0minus_VAMCFM.proddecme,
false);
4213 mela.getIORecord()->getUnweightedMEArray(p_prod_0minus_dec_0minus_VAMCFM.mearray);
4217 mela.computeProdDecP(p_prod_fa3_dec_fa3_VAMCFM.proddecme,
false);
4218 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3_dec_fa3_VAMCFM.mearray);
4222 mela.computeProdDecP(p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAMCFM.proddecme,
false);
4223 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAMCFM.mearray);
4227 mela.computeProdDecP(p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAMCFM.proddecme,
false);
4228 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAMCFM.mearray);
4231 bool computeL2WWproddec=(isZZWW==2 || decZZWW==2);
4232 if (computeL2WWproddec){
4233 mela.differentiate_HWW_HZZ=
true;
4234 if (isZZWW==2 && decZZWW!=2){
4235 if (vbfvhchannel==0){
4246 else if (isZZWW!=2 && decZZWW==2){
4247 if (vbfvhchannel==0){
4272 mela.computeProdDecP(p_prod_0mplusL2_dec_0mplusL2_VAMCFM.proddecme,
false);
4273 mela.getIORecord()->getUnweightedMEArray(p_prod_0mplusL2_dec_0mplusL2_VAMCFM.mearray);
4275 mela.differentiate_HWW_HZZ=
true;
4276 if (isZZWW==2 && decZZWW!=2){
4277 if (vbfvhchannel==0){
4288 else if (isZZWW!=2 && decZZWW==2){
4289 if (vbfvhchannel==0){
4306 mela.computeProdDecP(p_prod_fL2_dec_fL2_VAMCFM.proddecme,
false);
4307 mela.getIORecord()->getUnweightedMEArray(p_prod_fL2_dec_fL2_VAMCFM.mearray);
4314 for (
int ii=0; ii<
nmsq; ii++){
4315 for (
int jj=0; jj<
nmsq; jj++){
4316 p_prod_0mplusL2_dec_0mplusL2_VAMCFM.mearray[ii][jj]=p_prod_0mplus_dec_0mplus_VAMCFM.mearray[ii][jj]; p_prod_0mplusL2_dec_0mplusL2_VAMCFM.proddecme=p_prod_0mplus_dec_0mplus_VAMCFM.proddecme;
4317 p_prod_fL2_dec_fL2_VAMCFM.mearray[ii][jj]=p_prod_0mplus_dec_0mplus_VAMCFM.mearray[ii][jj]; p_prod_fL2_dec_fL2_VAMCFM.proddecme=p_prod_0mplus_dec_0mplus_VAMCFM.proddecme;
4323 bool testZAcoupl=
true;
4324 if (vbfvhchannel>=1){
4325 if (isZZWW==2 && decZZWW==1){
4329 else if (isZZWW==1 && decZZWW==2){
4333 else if (isZZWW==1 && decZZWW==1){
4337 else testZAcoupl=
false;
4340 if (isZZWW==2 && decZZWW==1){
4344 else if (isZZWW==1 && decZZWW==2){
4348 else if (isZZWW==1 && decZZWW==1){
4352 else testZAcoupl=
false;
4356 if (vbfvhchannel>=1){
4357 if (isZZWW==2 && decZZWW==1){
4360 mela.differentiate_HWW_HZZ=
true;
4362 else if (isZZWW==1 && decZZWW==2){
4364 mela.differentiate_HWW_HZZ=
true;
4366 else if (isZZWW==1 && decZZWW==1){
4368 mela.differentiate_HWW_HZZ=
true;
4372 if (isZZWW==2 && decZZWW==1){
4375 mela.differentiate_HWW_HZZ=
true;
4377 else if (isZZWW==1 && decZZWW==2){
4380 mela.differentiate_HWW_HZZ=
true;
4382 else if (isZZWW==1 && decZZWW==1){
4384 mela.differentiate_HWW_HZZ=
true;
4387 mela.computeProdDecP(p_prod_0minusZA_dec_0minusZA_VAMCFM.proddecme,
false);
4388 mela.getIORecord()->getUnweightedMEArray(p_prod_0minusZA_dec_0minusZA_VAMCFM.mearray);
4390 if (vbfvhchannel>=1){
4391 if (isZZWW==2 && decZZWW==1){
4394 mela.differentiate_HWW_HZZ=
true;
4396 else if (isZZWW==1 && decZZWW==2){
4398 mela.differentiate_HWW_HZZ=
true;
4400 else if (isZZWW==1 && decZZWW==1){
4402 mela.differentiate_HWW_HZZ=
true;
4406 if (isZZWW==2 && decZZWW==1){
4409 mela.differentiate_HWW_HZZ=
true;
4411 else if (isZZWW==1 && decZZWW==2){
4414 mela.differentiate_HWW_HZZ=
true;
4416 else if (isZZWW==1 && decZZWW==1){
4418 mela.differentiate_HWW_HZZ=
true;
4421 mela.computeProdDecP(p_prod_0minusAA_dec_0minusAA_VAMCFM.proddecme,
false);
4422 mela.getIORecord()->getUnweightedMEArray(p_prod_0minusAA_dec_0minusAA_VAMCFM.mearray);
4424 if (vbfvhchannel>=1){
4425 if (isZZWW==2 && decZZWW==1){
4429 mela.differentiate_HWW_HZZ=
true;
4431 else if (isZZWW==1 && decZZWW==2){
4434 mela.differentiate_HWW_HZZ=
true;
4436 else if (isZZWW==1 && decZZWW==1){
4439 mela.differentiate_HWW_HZZ=
true;
4443 if (isZZWW==2 && decZZWW==1){
4447 mela.differentiate_HWW_HZZ=
true;
4449 else if (isZZWW==1 && decZZWW==2){
4453 mela.differentiate_HWW_HZZ=
true;
4455 else if (isZZWW==1 && decZZWW==1){
4458 mela.differentiate_HWW_HZZ=
true;
4461 mela.computeProdDecP(p_prod_fa3ZA_dec_fa3ZA_VAMCFM.proddecme,
false);
4462 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3ZA_dec_fa3ZA_VAMCFM.mearray);
4464 if (vbfvhchannel>=1){
4465 if (isZZWW==2 && decZZWW==1){
4469 mela.differentiate_HWW_HZZ=
true;
4471 else if (isZZWW==1 && decZZWW==2){
4474 mela.differentiate_HWW_HZZ=
true;
4476 else if (isZZWW==1 && decZZWW==1){
4479 mela.differentiate_HWW_HZZ=
true;
4483 if (isZZWW==2 && decZZWW==1){
4487 mela.differentiate_HWW_HZZ=
true;
4489 else if (isZZWW==1 && decZZWW==2){
4493 mela.differentiate_HWW_HZZ=
true;
4495 else if (isZZWW==1 && decZZWW==1){
4498 mela.differentiate_HWW_HZZ=
true;
4501 mela.computeProdDecP(p_prod_fa3AA_dec_fa3AA_VAMCFM.proddecme,
false);
4502 mela.getIORecord()->getUnweightedMEArray(p_prod_fa3AA_dec_fa3AA_VAMCFM.mearray);
4505 for (
int ii=0; ii<
nmsq; ii++){
4506 for (
int jj=0; jj<
nmsq; jj++){
4507 p_prod_0minusZA_dec_0minusZA_VAMCFM.mearray[ii][jj]=p_prod_0mplus_dec_0mplus_VAMCFM.mearray[ii][jj]; p_prod_0minusZA_dec_0minusZA_VAMCFM.proddecme=p_prod_0mplus_dec_0mplus_VAMCFM.proddecme;
4508 p_prod_0minusAA_dec_0minusAA_VAMCFM.mearray[ii][jj]=p_prod_0mplus_dec_0mplus_VAMCFM.mearray[ii][jj]; p_prod_0minusAA_dec_0minusAA_VAMCFM.proddecme=p_prod_0mplus_dec_0mplus_VAMCFM.proddecme;
4509 p_prod_fa3ZA_dec_fa3ZA_VAMCFM.mearray[ii][jj]=p_prod_0mplus_dec_0mplus_VAMCFM.mearray[ii][jj]; p_prod_fa3ZA_dec_fa3ZA_VAMCFM.proddecme=p_prod_0mplus_dec_0mplus_VAMCFM.proddecme;
4510 p_prod_fa3AA_dec_fa3AA_VAMCFM.mearray[ii][jj]=p_prod_0mplus_dec_0mplus_VAMCFM.mearray[ii][jj]; p_prod_fa3AA_dec_fa3AA_VAMCFM.proddecme=p_prod_0mplus_dec_0mplus_VAMCFM.proddecme;
4517 mela.computeProdDecP(p_bkg_VAMCFM.proddecme,
false);
4518 mela.getIORecord()->getUnweightedMEArray(p_bkg_VAMCFM.mearray);
4519 if (motherflavor==0 && vbfvhchannel<2){
4520 for (
int r=-5;
r<=5;
r++){
4521 for (
int s=-5; s<=5; s++){
4530 mcfmme p_bkg_VAMCFM_rsindiv;
4531 mela.computeProdDecP(p_bkg_VAMCFM_rsindiv.proddecme,
false);
4532 mela.getIORecord()->getUnweightedMEArray(p_bkg_VAMCFM_rsindiv.mearray);
4533 p_bkg_VAMCFM_rssum.add(p_bkg_VAMCFM_rsindiv);
4544 cout <<
"Production variables:\n";
4545 cout <<
"\tmJJ = " << mjj << endl;
4546 cout <<
"\tPhi = " << Phi << endl;
4547 cout <<
"\tJHUGen (mass, width): (" << mh <<
", " << gah <<
")" << endl;
4549 cout <<
"Bkg" << endl;
4550 cout <<
"\tMCFM ME: " << p_bkg_VAMCFM.proddecme << endl;
4551 cout <<
"0mplus" << endl;
4552 cout <<
"\tJHUGen decay-alone: " << p_prod_0mplus_dec_0mplus_VAJHU.decme << endl;
4553 cout <<
"\tJHUGen prod.-alone: " << p_prod_0mplus_dec_0mplus_VAJHU.prodme << endl;
4554 cout <<
"\tJHUGen ME: " << p_prod_0mplus_dec_0mplus_VAJHU.proddecme << endl;
4555 cout <<
"\tMCFM ME: " << p_prod_0mplus_dec_0mplus_VAMCFM.proddecme << endl;
4556 cout <<
"0hplus" << endl;
4557 cout <<
"\tJHUGen decay-alone: " << p_prod_0hplus_dec_0hplus_VAJHU.decme << endl;
4558 cout <<
"\tJHUGen prod.-alone: " << p_prod_0hplus_dec_0hplus_VAJHU.prodme << endl;
4559 cout <<
"\tJHUGen ME: " << p_prod_0hplus_dec_0hplus_VAJHU.proddecme << endl;
4560 cout <<
"\tMCFM ME: " << p_prod_0hplus_dec_0hplus_VAMCFM.proddecme << endl;
4561 cout <<
"fa2" << endl;
4562 cout <<
"\tJHUGen decay-alone: " << p_prod_fa2_dec_fa2_VAJHU.decme << endl;
4563 cout <<
"\tJHUGen prod.-alone: " << p_prod_fa2_dec_fa2_VAJHU.prodme << endl;
4564 cout <<
"\tJHUGen ME: " << p_prod_fa2_dec_fa2_VAJHU.proddecme << endl;
4565 cout <<
"\tMCFM ME: " << p_prod_fa2_dec_fa2_VAMCFM.proddecme << endl;
4566 cout <<
"fa2, phia2=90" << endl;
4567 cout <<
"\tJHUGen decay-alone: " << p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.decme << endl;
4568 cout <<
"\tJHUGen prod.-alone: " << p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.prodme << endl;
4569 cout <<
"\tJHUGen ME: " << p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.proddecme << endl;
4570 cout <<
"\tMCFM ME: " << p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAMCFM.proddecme << endl;
4571 cout <<
"0minus" << endl;
4572 cout <<
"\tJHUGen decay-alone: " << p_prod_0minus_dec_0minus_VAJHU.decme << endl;
4573 cout <<
"\tJHUGen prod.-alone: " << p_prod_0minus_dec_0minus_VAJHU.prodme << endl;
4574 cout <<
"\tJHUGen ME: " << p_prod_0minus_dec_0minus_VAJHU.proddecme << endl;
4575 cout <<
"\tMCFM ME: " << p_prod_0minus_dec_0minus_VAMCFM.proddecme << endl;
4576 cout <<
"fa3" << endl;
4577 cout <<
"\tJHUGen decay-alone: " << p_prod_fa3_dec_fa3_VAJHU.decme << endl;
4578 cout <<
"\tJHUGen prod.-alone: " << p_prod_fa3_dec_fa3_VAJHU.prodme << endl;
4579 cout <<
"\tJHUGen ME: " << p_prod_fa3_dec_fa3_VAJHU.proddecme << endl;
4580 cout <<
"\tMCFM ME: " << p_prod_fa3_dec_fa3_VAMCFM.proddecme << endl;
4581 cout <<
"fa3, phia3=90" << endl;
4582 cout <<
"\tJHUGen decay-alone: " << p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.decme << endl;
4583 cout <<
"\tJHUGen prod.-alone: " << p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.prodme << endl;
4584 cout <<
"\tJHUGen ME: " << p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.proddecme << endl;
4585 cout <<
"\tMCFM ME: " << p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAMCFM.proddecme << endl;
4586 cout <<
"fa3, phia3=270" << endl;
4587 cout <<
"\tJHUGen decay-alone: " << p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.decme << endl;
4588 cout <<
"\tJHUGen prod.-alone: " << p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.prodme << endl;
4589 cout <<
"\tJHUGen ME: " << p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.proddecme << endl;
4590 cout <<
"\tMCFM ME: " << p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAMCFM.proddecme << endl;
4591 cout <<
"0mplusL2" << endl;
4592 cout <<
"\tJHUGen decay-alone: " << p_prod_0mplusL2_dec_0mplusL2_VAJHU.decme << endl;
4593 cout <<
"\tJHUGen prod.-alone: " << p_prod_0mplusL2_dec_0mplusL2_VAJHU.prodme << endl;
4594 cout <<
"\tJHUGen ME: " << p_prod_0mplusL2_dec_0mplusL2_VAJHU.proddecme << endl;
4595 cout <<
"\tMCFM ME: " << p_prod_0mplusL2_dec_0mplusL2_VAMCFM.proddecme << endl;
4596 cout <<
"fL2" << endl;
4597 cout <<
"\tJHUGen decay-alone: " << p_prod_fL2_dec_fL2_VAJHU.decme << endl;
4598 cout <<
"\tJHUGen prod.-alone: " << p_prod_fL2_dec_fL2_VAJHU.prodme << endl;
4599 cout <<
"\tJHUGen ME: " << p_prod_fL2_dec_fL2_VAJHU.proddecme << endl;
4600 cout <<
"\tMCFM ME: " << p_prod_fL2_dec_fL2_VAMCFM.proddecme << endl;
4601 cout <<
"0minusZA" << endl;
4602 cout <<
"\tJHUGen decay-alone: " << p_prod_0minusZA_dec_0minusZA_VAJHU.decme << endl;
4603 cout <<
"\tJHUGen prod.-alone: " << p_prod_0minusZA_dec_0minusZA_VAJHU.prodme << endl;
4604 cout <<
"\tJHUGen ME: " << p_prod_0minusZA_dec_0minusZA_VAJHU.proddecme << endl;
4605 cout <<
"\tMCFM ME: " << p_prod_0minusZA_dec_0minusZA_VAMCFM.proddecme << endl;
4606 cout <<
"fa3ZA" << endl;
4607 cout <<
"\tJHUGen decay-alone: " << p_prod_fa3ZA_dec_fa3ZA_VAJHU.decme << endl;
4608 cout <<
"\tJHUGen prod.-alone: " << p_prod_fa3ZA_dec_fa3ZA_VAJHU.prodme << endl;
4609 cout <<
"\tJHUGen ME: " << p_prod_fa3ZA_dec_fa3ZA_VAJHU.proddecme << endl;
4610 cout <<
"\tMCFM ME: " << p_prod_fa3ZA_dec_fa3ZA_VAMCFM.proddecme << endl;
4611 cout <<
"0minusAA" << endl;
4612 cout <<
"\tJHUGen decay-alone: " << p_prod_0minusAA_dec_0minusAA_VAJHU.decme << endl;
4613 cout <<
"\tJHUGen prod.-alone: " << p_prod_0minusAA_dec_0minusAA_VAJHU.prodme << endl;
4614 cout <<
"\tJHUGen ME: " << p_prod_0minusAA_dec_0minusAA_VAJHU.proddecme << endl;
4615 cout <<
"\tMCFM ME: " << p_prod_0minusAA_dec_0minusAA_VAMCFM.proddecme << endl;
4616 cout <<
"fa3AA" << endl;
4617 cout <<
"\tJHUGen decay-alone: " << p_prod_fa3AA_dec_fa3AA_VAJHU.decme << endl;
4618 cout <<
"\tJHUGen prod.-alone: " << p_prod_fa3AA_dec_fa3AA_VAJHU.prodme << endl;
4619 cout <<
"\tJHUGen ME: " << p_prod_fa3AA_dec_fa3AA_VAJHU.proddecme << endl;
4620 cout <<
"\tMCFM ME: " << p_prod_fa3AA_dec_fa3AA_VAMCFM.proddecme << endl;
4622 cout <<
"Arrays:" << endl;
4623 cout <<
"0mplus" << endl;
4624 cout <<
"\tJHUGen" << endl;
4625 p_prod_0mplus_dec_0mplus_VAJHU.printarray();
4626 cout <<
"\tMCFM" << endl;
4627 p_prod_0mplus_dec_0mplus_VAMCFM.printarray();
4628 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4629 for (
int ii=0; ii<
nmsq; ii++){
4630 for (
int jj=0; jj<
nmsq; jj++){
4632 if (p_prod_0mplus_dec_0mplus_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_0mplus_dec_0mplus_VAJHU.mearray[ii][jj]/p_prod_0mplus_dec_0mplus_VAMCFM.mearray[ii][jj];
4637 cout <<
"0hplus" << endl;
4638 cout <<
"\tJHUGen" << endl;
4639 p_prod_0hplus_dec_0hplus_VAJHU.printarray();
4640 cout <<
"\tMCFM" << endl;
4641 p_prod_0hplus_dec_0hplus_VAMCFM.printarray();
4642 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4643 for (
int ii=0; ii<
nmsq; ii++){
4644 for (
int jj=0; jj<
nmsq; jj++){
4646 if (p_prod_0hplus_dec_0hplus_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_0hplus_dec_0hplus_VAJHU.mearray[ii][jj]/p_prod_0hplus_dec_0hplus_VAMCFM.mearray[ii][jj];
4651 cout <<
"fa2" << endl;
4652 cout <<
"\tJHUGen" << endl;
4653 p_prod_fa2_dec_fa2_VAJHU.printarray();
4654 cout <<
"\tMCFM" << endl;
4655 p_prod_fa2_dec_fa2_VAMCFM.printarray();
4656 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4657 for (
int ii=0; ii<
nmsq; ii++){
4658 for (
int jj=0; jj<
nmsq; jj++){
4660 if (p_prod_fa2_dec_fa2_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_fa2_dec_fa2_VAJHU.mearray[ii][jj]/p_prod_fa2_dec_fa2_VAMCFM.mearray[ii][jj];
4665 cout <<
"fa2, phia2=90" << endl;
4666 cout <<
"\tJHUGen" << endl;
4667 p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.printarray();
4668 cout <<
"\tMCFM" << endl;
4669 p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAMCFM.printarray();
4670 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4671 for (
int ii=0; ii<
nmsq; ii++){
4672 for (
int jj=0; jj<
nmsq; jj++){
4674 if (p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAJHU.mearray[ii][jj]/p_prod_fa2_phia2_90_dec_fa2_phia2_90_VAMCFM.mearray[ii][jj];
4679 cout <<
"0minus" << endl;
4680 cout <<
"\tJHUGen" << endl;
4681 p_prod_0minus_dec_0minus_VAJHU.printarray();
4682 cout <<
"\tMCFM" << endl;
4683 p_prod_0minus_dec_0minus_VAMCFM.printarray();
4684 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4685 for (
int ii=0; ii<
nmsq; ii++){
4686 for (
int jj=0; jj<
nmsq; jj++){
4688 if (p_prod_0minus_dec_0minus_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_0minus_dec_0minus_VAJHU.mearray[ii][jj]/p_prod_0minus_dec_0minus_VAMCFM.mearray[ii][jj];
4693 cout <<
"fa3" << endl;
4694 cout <<
"\tJHUGen" << endl;
4695 p_prod_fa3_dec_fa3_VAJHU.printarray();
4696 cout <<
"\tMCFM" << endl;
4697 p_prod_fa3_dec_fa3_VAMCFM.printarray();
4698 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4699 for (
int ii=0; ii<
nmsq; ii++){
4700 for (
int jj=0; jj<
nmsq; jj++){
4702 if (p_prod_fa3_dec_fa3_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_fa3_dec_fa3_VAJHU.mearray[ii][jj]/p_prod_fa3_dec_fa3_VAMCFM.mearray[ii][jj];
4707 cout <<
"fa3, phia3=90" << endl;
4708 cout <<
"\tJHUGen" << endl;
4709 p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.printarray();
4710 cout <<
"\tMCFM" << endl;
4711 p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAMCFM.printarray();
4712 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4713 for (
int ii=0; ii<
nmsq; ii++){
4714 for (
int jj=0; jj<
nmsq; jj++){
4716 if (p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAJHU.mearray[ii][jj]/p_prod_fa3_phia3_90_dec_fa3_phia3_90_VAMCFM.mearray[ii][jj];
4721 cout <<
"fa3, phia3=270" << endl;
4722 cout <<
"\tJHUGen" << endl;
4723 p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.printarray();
4724 cout <<
"\tMCFM" << endl;
4725 p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAMCFM.printarray();
4726 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4727 for (
int ii=0; ii<
nmsq; ii++){
4728 for (
int jj=0; jj<
nmsq; jj++){
4730 if (p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAJHU.mearray[ii][jj]/p_prod_fa3_phia3_270_dec_fa3_phia3_270_VAMCFM.mearray[ii][jj];
4735 cout <<
"0mplusL2" << endl;
4736 cout <<
"\tJHUGen" << endl;
4737 p_prod_0mplusL2_dec_0mplusL2_VAJHU.printarray();
4738 cout <<
"\tMCFM" << endl;
4739 p_prod_0mplusL2_dec_0mplusL2_VAMCFM.printarray();
4740 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4741 for (
int ii=0; ii<
nmsq; ii++){
4742 for (
int jj=0; jj<
nmsq; jj++){
4744 if (p_prod_0mplusL2_dec_0mplusL2_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_0mplusL2_dec_0mplusL2_VAJHU.mearray[ii][jj]/p_prod_0mplusL2_dec_0mplusL2_VAMCFM.mearray[ii][jj];
4749 cout <<
"fL2" << endl;
4750 cout <<
"\tJHUGen" << endl;
4751 p_prod_fL2_dec_fL2_VAJHU.printarray();
4752 cout <<
"\tMCFM" << endl;
4753 p_prod_fL2_dec_fL2_VAMCFM.printarray();
4754 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4755 for (
int ii=0; ii<
nmsq; ii++){
4756 for (
int jj=0; jj<
nmsq; jj++){
4758 if (p_prod_fL2_dec_fL2_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_fL2_dec_fL2_VAJHU.mearray[ii][jj]/p_prod_fL2_dec_fL2_VAMCFM.mearray[ii][jj];
4763 cout <<
"0minusZA" << endl;
4764 cout <<
"\tJHUGen" << endl;
4765 p_prod_0minusZA_dec_0minusZA_VAJHU.printarray();
4766 cout <<
"\tMCFM" << endl;
4767 p_prod_0minusZA_dec_0minusZA_VAMCFM.printarray();
4768 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4769 for (
int ii=0; ii<
nmsq; ii++){
4770 for (
int jj=0; jj<
nmsq; jj++){
4772 if (p_prod_0minusZA_dec_0minusZA_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_0minusZA_dec_0minusZA_VAJHU.mearray[ii][jj]/p_prod_0minusZA_dec_0minusZA_VAMCFM.mearray[ii][jj];
4777 cout <<
"fa3ZA" << endl;
4778 cout <<
"\tJHUGen" << endl;
4779 p_prod_fa3ZA_dec_fa3ZA_VAJHU.printarray();
4780 cout <<
"\tMCFM" << endl;
4781 p_prod_fa3ZA_dec_fa3ZA_VAMCFM.printarray();
4782 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4783 for (
int ii=0; ii<
nmsq; ii++){
4784 for (
int jj=0; jj<
nmsq; jj++){
4786 if (p_prod_fa3ZA_dec_fa3ZA_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_fa3ZA_dec_fa3ZA_VAJHU.mearray[ii][jj]/p_prod_fa3ZA_dec_fa3ZA_VAMCFM.mearray[ii][jj];
4791 cout <<
"0minusAA" << endl;
4792 cout <<
"\tJHUGen" << endl;
4793 p_prod_0minusAA_dec_0minusAA_VAJHU.printarray();
4794 cout <<
"\tMCFM" << endl;
4795 p_prod_0minusAA_dec_0minusAA_VAMCFM.printarray();
4796 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4797 for (
int ii=0; ii<
nmsq; ii++){
4798 for (
int jj=0; jj<
nmsq; jj++){
4800 if (p_prod_0minusAA_dec_0minusAA_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_0minusAA_dec_0minusAA_VAJHU.mearray[ii][jj]/p_prod_0minusAA_dec_0minusAA_VAMCFM.mearray[ii][jj];
4805 cout <<
"fa3AA" << endl;
4806 cout <<
"\tJHUGen" << endl;
4807 p_prod_fa3AA_dec_fa3AA_VAJHU.printarray();
4808 cout <<
"\tMCFM" << endl;
4809 p_prod_fa3AA_dec_fa3AA_VAMCFM.printarray();
4810 cout <<
"\tJHUGen/MCFM Ratio" << endl;
4811 for (
int ii=0; ii<
nmsq; ii++){
4812 for (
int jj=0; jj<
nmsq; jj++){
4814 if (p_prod_fa3AA_dec_fa3AA_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_fa3AA_dec_fa3AA_VAJHU.mearray[ii][jj]/p_prod_fa3AA_dec_fa3AA_VAMCFM.mearray[ii][jj];
4819 cout <<
"Bkg" << endl;
4820 cout <<
"\tMCFM" << endl;
4821 p_bkg_VAMCFM.printarray();
4822 cout <<
"\tMCFM 0mplus/Bkg Ratio" << endl;
4823 for (
int ii=0; ii<
nmsq; ii++){
4824 for (
int jj=0; jj<
nmsq; jj++){
4826 if (p_bkg_VAMCFM.mearray[ii][jj]!=0.)
cout << p_prod_0mplus_dec_0mplus_VAJHU.mearray[ii][jj]/p_bkg_VAMCFM.mearray[ii][jj];
4831 if (motherflavor==0 && vbfvhchannel<2){
4832 cout <<
"Bkg manual sum" << endl;
4833 cout <<
"\tMCFM" << endl;
4834 p_bkg_VAMCFM_rssum.printarray();
4835 cout <<
"\tMCFM Bkg (re-sum)/Bkg Ratio" << endl;
4836 bool wrongRatio=
false;
4837 for (
int ii=0; ii<
nmsq; ii++){
4838 for (
int jj=0; jj<
nmsq; jj++){
4841 if (p_bkg_VAMCFM.mearray[ii][jj]!=0.) rr = p_bkg_VAMCFM_rssum.mearray[ii][jj]/p_bkg_VAMCFM.mearray[ii][jj];
4842 else if (p_bkg_VAMCFM_rssum.mearray[ii][jj]!=0.) rr = -9999;
4844 if (rr!=0. && rr!=4.) wrongRatio=
true;
4850 for (
int r=-5;
r<=5;
r++){
4851 for (
int s=-5; s<=5; s++){
4860 mcfmme p_bkg_VAMCFM_rsindiv;
4861 mela.computeProdDecP(p_bkg_VAMCFM_rsindiv.proddecme,
false);
4862 mela.getIORecord()->getUnweightedMEArray(p_bkg_VAMCFM_rsindiv.mearray);
4863 if (p_bkg_VAMCFM_rsindiv.proddecme>0.){
4865 mela.computeProdDecP(p_bkg_VAMCFM_rsindiv.proddecme,
false);
4867 cout <<
"Outgoing id1, id2 = " <<
r <<
" , " << s << endl;
4869 for (
int ii=0; ii<
nmsq; ii++){
4870 for (
int jj=0; jj<
nmsq; jj++){
4872 cout << p_bkg_VAMCFM_rsindiv.mearray[ii][jj];
4878 mela.setVerbosity(verbosity);
4890 mela.resetInputEvent();
4893 TUtil::SetCKMElements(&bkpvckm_ud, &bkpvckm_us, &bkpvckm_cd, &bkpvckm_cs, &bkpvckm_ts, &bkpvckm_tb, &bkpvckm_ub, &bkpvckm_cb, &bkpvckm_td);
4896 cout.rdbuf(coutbuf);
4898 mela.setVerbosity(bkpverbosity);