~maddevelopers/mg5amcnlo/3.3.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksrealew/%SubProcesses%P0_dxd_ttx%matrix_5.f

• browse files at revision 980
• compare with another revision
• download diff
• download tarball
• view history from revision 980

Show diffs side-by-side

added

removed

tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksrealew/%SubProcesses%P0_dxd_ttx%matrix_5.f

9

9

C      orders.inc (NLO_ORDERS)

10

10

C     

11

11

C     

12

 

C     Process: d~ a > t t~ d~ [ real = QED QCD ] QCD^2<=4 QED^2<=2

13

 

C     Process: s~ a > t t~ s~ [ real = QED QCD ] QCD^2<=4 QED^2<=2

 

12

C     Process: d~ g > t t~ d~ [ real = QCD QED ] QCD^2<=4 QED^2<=2

 

13

C     Process: s~ g > t t~ s~ [ real = QCD QED ] QCD^2<=4 QED^2<=2

14

14

C     

15

15

C     

16

16

C     CONSTANTS

19

19

      INTEGER    NEXTERNAL

20

20

      PARAMETER (NEXTERNAL=5)

21

21

      INTEGER NSQAMPSO

22

 

      PARAMETER (NSQAMPSO=1)

 

22

      PARAMETER (NSQAMPSO=3)

23

23

C     

24

24

C     ARGUMENTS 

25

25

C     

109

109

C     and helicities

110

110

C     for the point in phase space P(0:3,NEXTERNAL)

111

111

C     

112

 

C     Process: d~ a > t t~ d~ [ real = QED QCD ] QCD^2<=4 QED^2<=2

113

 

C     Process: s~ a > t t~ s~ [ real = QED QCD ] QCD^2<=4 QED^2<=2

 

112

C     Process: d~ g > t t~ d~ [ real = QCD QED ] QCD^2<=4 QED^2<=2

 

113

C     Process: s~ g > t t~ s~ [ real = QCD QED ] QCD^2<=4 QED^2<=2

114

114

C     

115

115

      IMPLICIT NONE

116

116

C     

120

120

      INTEGER     NCOMB

121

121

      PARAMETER ( NCOMB=32)

122

122

      INTEGER NSQAMPSO

123

 

      PARAMETER (NSQAMPSO=1)

 

123

      PARAMETER (NSQAMPSO=3)

124

124

C     

125

125

C     ARGUMENTS 

126

126

C     

168

168

      DATA GOODHEL/NCOMB*.FALSE./

169

169

      INTEGER NTRY

170

170

      DATA NTRY/0/

171

 

      DATA IDEN/12/

 

171

      DATA IDEN/96/

172

172

C     ----------

173

173

C     BEGIN CODE

174

174

C     ----------

234

234

C     Returns amplitude squared summed/avg over colors

235

235

C     for the point with external lines W(0:6,NEXTERNAL)

236

236

C     

237

 

C     Process: d~ a > t t~ d~ [ real = QED QCD ] QCD^2<=4 QED^2<=2

238

 

C     Process: s~ a > t t~ s~ [ real = QED QCD ] QCD^2<=4 QED^2<=2

 

237

C     Process: d~ g > t t~ d~ [ real = QCD QED ] QCD^2<=4 QED^2<=2

 

238

C     Process: s~ g > t t~ s~ [ real = QCD QED ] QCD^2<=4 QED^2<=2

239

239

C     

240

240

      IMPLICIT NONE

241

241

C     

242

242

C     CONSTANTS

243

243

C     

244

244

      INTEGER    NGRAPHS

245

 

      PARAMETER (NGRAPHS=4)

 

245

      PARAMETER (NGRAPHS=13)

246

246

      INTEGER    NWAVEFUNCS, NCOLOR

247

 

      PARAMETER (NWAVEFUNCS=8, NCOLOR=2)

 

247

      PARAMETER (NWAVEFUNCS=12, NCOLOR=4)

248

248

      INTEGER NAMPSO, NSQAMPSO

249

 

      PARAMETER (NAMPSO=1, NSQAMPSO=1)

 

249

      PARAMETER (NAMPSO=2, NSQAMPSO=3)

250

250

      REAL*8     ZERO

251

251

      PARAMETER (ZERO=0D0)

252

252

      COMPLEX*16 IMAG1

268

268

      REAL*8  CF(NCOLOR,NCOLOR)

269

269

      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP(NCOLOR,NAMPSO), W(8

270

270

     $ ,NWAVEFUNCS)

271

 

      COMPLEX*16 TMP_JAMP(3)

 

271

      COMPLEX*16 TMP_JAMP(0)

272

272

C     

273

273

C     FUNCTION

274

274

C     

276

276

C     

277

277

C     COLOR DATA

278

278

C     

279

 

      DATA (CF(I,  1),I=  1,  2) /9.000000000000000D+00

280

 

     $ ,3.000000000000000D+00/

281

 

C     1 T(1,4) T(3,5)

282

 

      DATA (CF(I,  2),I=  1,  2) /3.000000000000000D+00

283

 

     $ ,9.000000000000000D+00/

284

 

C     1 T(1,5) T(3,4)

 

279

      DATA (CF(I,  1),I=  1,  4) /1.200000000000000D+01

 

280

     $ ,4.000000000000000D+00,0.000000000000000D+00,4.000000000000000D

 

281

     $ +00/

 

282

C     1 T(1,4) T(2,3,5)

 

283

      DATA (CF(I,  2),I=  1,  4) /4.000000000000000D+00

 

284

     $ ,1.200000000000000D+01,4.000000000000000D+00,0.000000000000000D

 

285

     $ +00/

 

286

C     1 T(1,5) T(2,3,4)

 

287

      DATA (CF(I,  3),I=  1,  4) /0.000000000000000D+00

 

288

     $ ,4.000000000000000D+00,1.200000000000000D+01,4.000000000000000D

 

289

     $ +00/

 

290

C     1 T(2,1,4) T(3,5)

 

291

      DATA (CF(I,  4),I=  1,  4) /4.000000000000000D+00

 

292

     $ ,0.000000000000000D+00,4.000000000000000D+00,1.200000000000000D

 

293

     $ +01/

 

294

C     1 T(2,1,5) T(3,4)

285

295

C     ----------

286

296

C     BEGIN CODE

287

297

C     ----------

290

300

      CALL OXXXXX(P(0,3),MDL_MT,NHEL(3),+1*IC(3),W(1,3))

291

301

      CALL IXXXXX(P(0,4),MDL_MT,NHEL(4),-1*IC(4),W(1,4))

292

302

      CALL IXXXXX(P(0,5),ZERO,NHEL(5),-1*IC(5),W(1,5))

293

 

      CALL FFV1_1(W(1,1),W(1,2),GC_1,ZERO,ZERO,W(1,6))

294

 

      CALL FFV1P0_3(W(1,4),W(1,3),GC_11,ZERO,ZERO,W(1,7))

 

303

      CALL FFV1_1(W(1,1),W(1,2),GC_11,ZERO,ZERO,W(1,6))

 

304

      CALL FFV1P0_3(W(1,4),W(1,3),GC_2,ZERO,ZERO,W(1,7))

295

305

C     Amplitude(s) for diagram number 1

296

 

      CALL FFV1_0(W(1,5),W(1,6),W(1,7),GC_11,AMP(1))

297

 

      CALL FFV1P0_3(W(1,5),W(1,1),GC_11,ZERO,ZERO,W(1,6))

298

 

      CALL FFV1_1(W(1,3),W(1,2),GC_2,MDL_MT,MDL_WT,W(1,8))

 

306

      CALL FFV1_0(W(1,5),W(1,6),W(1,7),GC_1,AMP(1))

 

307

      CALL FFV1P0_3(W(1,4),W(1,3),GC_11,ZERO,ZERO,W(1,8))

299

308

C     Amplitude(s) for diagram number 2

300

 

      CALL FFV1_0(W(1,4),W(1,8),W(1,6),GC_11,AMP(2))

301

 

      CALL FFV1_2(W(1,4),W(1,2),GC_2,MDL_MT,MDL_WT,W(1,8))

 

309

      CALL FFV1_0(W(1,5),W(1,6),W(1,8),GC_11,AMP(2))

 

310

      CALL FFV2_5_3(W(1,4),W(1,3),-GC_50,GC_58,MDL_MZ,MDL_WZ,W(1,9))

302

311

C     Amplitude(s) for diagram number 3

303

 

      CALL FFV1_0(W(1,8),W(1,3),W(1,6),GC_11,AMP(3))

304

 

      CALL FFV1_2(W(1,5),W(1,2),GC_1,ZERO,ZERO,W(1,8))

 

312

      CALL FFV2_3_0(W(1,5),W(1,6),W(1,9),GC_50,GC_58,AMP(3))

 

313

      CALL FFV1P0_3(W(1,5),W(1,1),GC_1,ZERO,ZERO,W(1,6))

 

314

      CALL FFV1_1(W(1,3),W(1,2),GC_11,MDL_MT,MDL_WT,W(1,10))

305

315

C     Amplitude(s) for diagram number 4

306

 

      CALL FFV1_0(W(1,8),W(1,1),W(1,7),GC_11,AMP(4))

307

 

C     JAMPs contributing to orders QCD=2 QED=1

308

 

      TMP_JAMP(2) = AMP(3) +  AMP(4)  ! used 2 times

309

 

      TMP_JAMP(1) = AMP(1) +  AMP(2)  ! used 2 times

310

 

      TMP_JAMP(3) = TMP_JAMP(2) +  TMP_JAMP(1)  ! used 2 times

311

 

      JAMP(1,1) = (5.000000000000000D-01)*TMP_JAMP(3)

312

 

      JAMP(2,1) = (-1.666666666666667D-01)*TMP_JAMP(3)

 

316

      CALL FFV1_0(W(1,4),W(1,10),W(1,6),GC_2,AMP(4))

 

317

      CALL FFV1P0_3(W(1,5),W(1,1),GC_11,ZERO,ZERO,W(1,11))

 

318

C     Amplitude(s) for diagram number 5

 

319

      CALL FFV1_0(W(1,4),W(1,10),W(1,11),GC_11,AMP(5))

 

320

      CALL FFV2_3_3(W(1,5),W(1,1),GC_50,GC_58,MDL_MZ,MDL_WZ,W(1,12))

 

321

C     Amplitude(s) for diagram number 6

 

322

      CALL FFV2_5_0(W(1,4),W(1,10),W(1,12),-GC_50,GC_58,AMP(6))

 

323

      CALL FFV1_2(W(1,4),W(1,2),GC_11,MDL_MT,MDL_WT,W(1,10))

 

324

C     Amplitude(s) for diagram number 7

 

325

      CALL FFV1_0(W(1,10),W(1,3),W(1,6),GC_2,AMP(7))

 

326

C     Amplitude(s) for diagram number 8

 

327

      CALL FFV1_0(W(1,10),W(1,3),W(1,11),GC_11,AMP(8))

 

328

C     Amplitude(s) for diagram number 9

 

329

      CALL FFV2_5_0(W(1,10),W(1,3),W(1,12),-GC_50,GC_58,AMP(9))

 

330

C     Amplitude(s) for diagram number 10

 

331

      CALL VVV1_0(W(1,11),W(1,2),W(1,8),GC_10,AMP(10))

 

332

      CALL FFV1_2(W(1,5),W(1,2),GC_11,ZERO,ZERO,W(1,11))

 

333

C     Amplitude(s) for diagram number 11

 

334

      CALL FFV1_0(W(1,11),W(1,1),W(1,7),GC_1,AMP(11))

 

335

C     Amplitude(s) for diagram number 12

 

336

      CALL FFV1_0(W(1,11),W(1,1),W(1,8),GC_11,AMP(12))

 

337

C     Amplitude(s) for diagram number 13

 

338

      CALL FFV2_3_0(W(1,11),W(1,1),W(1,9),GC_50,GC_58,AMP(13))

 

339

C     JAMPs contributing to orders QCD=3 QED=0

 

340

      JAMP(1,1) = (5.000000000000000D-01)*AMP(5)+((0.000000000000000D

 

341

     $ +00,5.000000000000000D-01))*AMP(10)+(5.000000000000000D-01)

 

342

     $ *AMP(12)

 

343

      JAMP(2,1) = (-1.666666666666667D-01)*AMP(5)+(-1.666666666666667D

 

344

     $ -01)*AMP(8)

 

345

      JAMP(3,1) = (5.000000000000000D-01)*AMP(2)+(5.000000000000000D

 

346

     $ -01)*AMP(8)+((0.000000000000000D+00,-5.000000000000000D-01))

 

347

     $ *AMP(10)

 

348

      JAMP(4,1) = (-1.666666666666667D-01)*AMP(2)+(-1.666666666666667D

 

349

     $ -01)*AMP(12)

 

350

C     JAMPs contributing to orders QCD=1 QED=2

 

351

      JAMP(1,2) = 0D0

 

352

      JAMP(2,2) = AMP(4)+AMP(6)+AMP(7)+AMP(9)

 

353

      JAMP(3,2) = 0D0

 

354

      JAMP(4,2) = AMP(1)+AMP(3)+AMP(11)+AMP(13)

313

355

 

314

356

      DO I=0,NSQAMPSO

315

357

        RES(I)=0D0

351

393

C     

352

394

      IMPLICIT NONE

353

395

      INTEGER NAMPSO, NSQAMPSO

354

 

      PARAMETER (NAMPSO=1, NSQAMPSO=1)

 

396

      PARAMETER (NAMPSO=2, NSQAMPSO=3)

355

397

      INTEGER NSPLITORDERS

356

398

      PARAMETER (NSPLITORDERS=2)

357

399

C     

363

405

C     

364

406

      INTEGER I, SQORDERS(NSPLITORDERS)

365

407

      INTEGER AMPSPLITORDERS(NAMPSO,NSPLITORDERS)

366

 

      DATA (AMPSPLITORDERS(  1,I),I=  1,  2) /    2,    1/

 

408

      DATA (AMPSPLITORDERS(  1,I),I=  1,  2) /    3,    0/

 

409

      DATA (AMPSPLITORDERS(  2,I),I=  1,  2) /    1,    2/

367

410

C     

368

411

C     FUNCTION

369

412

C     

388

431

C     

389

432

      IMPLICIT NONE

390

433

      INTEGER NSQAMPSO

391

 

      PARAMETER (NSQAMPSO=1)

 

434

      PARAMETER (NSQAMPSO=3)

392

435

      INTEGER NSPLITORDERS

393

436

      PARAMETER (NSPLITORDERS=2)

394

437

C     

401

444

      INTEGER I,J

402

445

      INTEGER SQSPLITORDERS(NSQAMPSO,NSPLITORDERS)

403

446

C     the values listed below are for QCD, QED

404

 

      DATA (SQSPLITORDERS(  1,I),I=  1,  2) /    4,    2/

 

447

      DATA (SQSPLITORDERS(  1,I),I=  1,  2) /    6,    0/

 

448

      DATA (SQSPLITORDERS(  2,I),I=  1,  2) /    4,    2/

 

449

      DATA (SQSPLITORDERS(  3,I),I=  1,  2) /    2,    4/

405

450

C     

406

451

C     BEGIN CODE

407

452

C     

431

476

C     

432

477

      IMPLICIT NONE

433

478

      INTEGER NSQAMPSO

434

 

      PARAMETER (NSQAMPSO=1)

 

479

      PARAMETER (NSQAMPSO=3)

435

480

      INTEGER NSPLITORDERS

436

481

      PARAMETER (NSPLITORDERS=2)

437

482

C     

444

489

      INTEGER I

445

490

      INTEGER SQSPLITORDERS(NSQAMPSO,NSPLITORDERS)

446

491

C     the values listed below are for QCD, QED

447

 

      DATA (SQSPLITORDERS(  1,I),I=  1,  2) /    4,    2/

 

492

      DATA (SQSPLITORDERS(  1,I),I=  1,  2) /    6,    0/

 

493

      DATA (SQSPLITORDERS(  2,I),I=  1,  2) /    4,    2/

 

494

      DATA (SQSPLITORDERS(  3,I),I=  1,  2) /    2,    4/

448

495

C     

449

496

C     BEGIN CODE

450

497

C     

472

519

C     

473

520

      IMPLICIT NONE

474

521

      INTEGER NSQAMPSO

475

 

      PARAMETER (NSQAMPSO=1)

 

522

      PARAMETER (NSQAMPSO=3)

476

523

      INTEGER NSQSO

477

524

 

478

525

      NSQSO=NSQAMPSO

• Older »

Loggerhead is a web-based interface for Breezy
Version: 2.0.1