~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_uux_ttx%matrix_3.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_uux_ttx%matrix_3.f

9

9

C      orders.inc (NLO_ORDERS)

10

10

C     

11

11

C     

12

 

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

13

 

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

 

12

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

 

13

C     Process: c c~ > t t~ a [ 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=3)

 

22

      PARAMETER (NSQAMPSO=1)

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: u u~ > t t~ g [ real = QED QCD ] QCD^2<=4 QED^2<=2

113

 

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

 

112

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

 

113

C     Process: c c~ > t t~ a [ 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=3)

 

123

      PARAMETER (NSQAMPSO=1)

124

124

C     

125

125

C     ARGUMENTS 

126

126

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: u u~ > t t~ g [ real = QED QCD ] QCD^2<=4 QED^2<=2

238

 

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

 

237

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

 

238

C     Process: c c~ > t t~ a [ 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=13)

 

245

      PARAMETER (NGRAPHS=4)

246

246

      INTEGER    NWAVEFUNCS, NCOLOR

247

 

      PARAMETER (NWAVEFUNCS=10, NCOLOR=4)

 

247

      PARAMETER (NWAVEFUNCS=7, NCOLOR=2)

248

248

      INTEGER NAMPSO, NSQAMPSO

249

 

      PARAMETER (NAMPSO=2, NSQAMPSO=3)

 

249

      PARAMETER (NAMPSO=1, NSQAMPSO=1)

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(0)

 

271

      COMPLEX*16 TMP_JAMP(3)

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,  4) /1.200000000000000D+01

280

 

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

281

 

     $ +00/

282

 

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

283

 

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

284

 

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

285

 

     $ +00/

286

 

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

287

 

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

288

 

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

289

 

     $ +00/

290

 

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

291

 

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

292

 

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

293

 

     $ +01/

294

 

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

 

279

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

 

280

     $ ,3.000000000000000D+00/

 

281

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

 

282

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

 

283

     $ ,9.000000000000000D+00/

 

284

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

295

285

C     ----------

296

286

C     BEGIN CODE

297

287

C     ----------

301

291

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

302

292

      CALL VXXXXX(P(0,5),ZERO,NHEL(5),+1*IC(5),W(1,5))

303

293

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

304

 

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

 

294

      CALL FFV1_1(W(1,3),W(1,5),GC_2,MDL_MT,MDL_WT,W(1,7))

305

295

C     Amplitude(s) for diagram number 1

306

 

      CALL VVV1_0(W(1,6),W(1,7),W(1,5),GC_10,AMP(1))

307

 

      CALL FFV1P0_3(W(1,1),W(1,2),GC_2,ZERO,ZERO,W(1,8))

308

 

      CALL FFV1_1(W(1,3),W(1,5),GC_11,MDL_MT,MDL_WT,W(1,9))

 

296

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

 

297

      CALL FFV1_2(W(1,4),W(1,5),GC_2,MDL_MT,MDL_WT,W(1,7))

309

298

C     Amplitude(s) for diagram number 2

310

 

      CALL FFV1_0(W(1,4),W(1,9),W(1,8),GC_2,AMP(2))

 

299

      CALL FFV1_0(W(1,7),W(1,3),W(1,6),GC_11,AMP(2))

 

300

      CALL FFV1_2(W(1,1),W(1,5),GC_2,ZERO,ZERO,W(1,7))

 

301

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

311

302

C     Amplitude(s) for diagram number 3

312

 

      CALL FFV1_0(W(1,4),W(1,9),W(1,6),GC_11,AMP(3))

313

 

      CALL FFV2_5_3(W(1,1),W(1,2),-GC_50,GC_58,MDL_MZ,MDL_WZ,W(1,10))

 

303

      CALL FFV1_0(W(1,7),W(1,2),W(1,6),GC_11,AMP(3))

 

304

      CALL FFV1_1(W(1,2),W(1,5),GC_2,ZERO,ZERO,W(1,7))

314

305

C     Amplitude(s) for diagram number 4

315

 

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

316

 

      CALL FFV1_2(W(1,4),W(1,5),GC_11,MDL_MT,MDL_WT,W(1,9))

317

 

C     Amplitude(s) for diagram number 5

318

 

      CALL FFV1_0(W(1,9),W(1,3),W(1,8),GC_2,AMP(5))

319

 

C     Amplitude(s) for diagram number 6

320

 

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

321

 

C     Amplitude(s) for diagram number 7

322

 

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

323

 

      CALL FFV1_2(W(1,1),W(1,5),GC_11,ZERO,ZERO,W(1,9))

324

 

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

325

 

C     Amplitude(s) for diagram number 8

326

 

      CALL FFV1_0(W(1,9),W(1,2),W(1,10),GC_2,AMP(8))

327

 

C     Amplitude(s) for diagram number 9

328

 

      CALL FFV1_0(W(1,9),W(1,2),W(1,7),GC_11,AMP(9))

329

 

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

330

 

C     Amplitude(s) for diagram number 10

331

 

      CALL FFV2_5_0(W(1,9),W(1,2),W(1,6),-GC_50,GC_58,AMP(10))

332

 

      CALL FFV1_1(W(1,2),W(1,5),GC_11,ZERO,ZERO,W(1,9))

333

 

C     Amplitude(s) for diagram number 11

334

 

      CALL FFV1_0(W(1,1),W(1,9),W(1,10),GC_2,AMP(11))

335

 

C     Amplitude(s) for diagram number 12

336

 

      CALL FFV1_0(W(1,1),W(1,9),W(1,7),GC_11,AMP(12))

337

 

C     Amplitude(s) for diagram number 13

338

 

      CALL FFV2_5_0(W(1,1),W(1,9),W(1,6),-GC_50,GC_58,AMP(13))

339

 

C     JAMPs contributing to orders QCD=3 QED=0

340

 

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

341

 

     $ -01)*AMP(6)

342

 

      JAMP(2,1) = ((0.000000000000000D+00,5.000000000000000D-01))

343

 

     $ *AMP(1)+(-5.000000000000000D-01)*AMP(3)+(-5.000000000000000D-01)

344

 

     $ *AMP(9)

345

 

      JAMP(3,1) = ((0.000000000000000D+00,-5.000000000000000D-01))

346

 

     $ *AMP(1)+(-5.000000000000000D-01)*AMP(6)+(-5.000000000000000D-01)

347

 

     $ *AMP(12)

348

 

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

349

 

     $ -01)*AMP(12)

350

 

C     JAMPs contributing to orders QCD=1 QED=2

351

 

      JAMP(1,2) = (-1.000000000000000D+00)*AMP(2)+(-1.000000000000000D

352

 

     $ +00)*AMP(4)+(-1.000000000000000D+00)*AMP(5)+(

353

 

     $ -1.000000000000000D+00)*AMP(7)

354

 

      JAMP(2,2) = 0D0

355

 

      JAMP(3,2) = 0D0

356

 

      JAMP(4,2) = (-1.000000000000000D+00)*AMP(8)+(-1.000000000000000D

357

 

     $ +00)*AMP(10)+(-1.000000000000000D+00)*AMP(11)+(

358

 

     $ -1.000000000000000D+00)*AMP(13)

 

306

      CALL FFV1_0(W(1,1),W(1,7),W(1,6),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) = (1.666666666666667D-01)*TMP_JAMP(3)

 

312

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

359

313

 

360

314

      DO I=0,NSQAMPSO

361

315

        RES(I)=0D0

397

351

C     

398

352

      IMPLICIT NONE

399

353

      INTEGER NAMPSO, NSQAMPSO

400

 

      PARAMETER (NAMPSO=2, NSQAMPSO=3)

 

354

      PARAMETER (NAMPSO=1, NSQAMPSO=1)

401

355

      INTEGER NSPLITORDERS

402

356

      PARAMETER (NSPLITORDERS=2)

403

357

C     

409

363

C     

410

364

      INTEGER I, SQORDERS(NSPLITORDERS)

411

365

      INTEGER AMPSPLITORDERS(NAMPSO,NSPLITORDERS)

412

 

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

413

 

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

 

366

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

414

367

C     

415

368

C     FUNCTION

416

369

C     

435

388

C     

436

389

      IMPLICIT NONE

437

390

      INTEGER NSQAMPSO

438

 

      PARAMETER (NSQAMPSO=3)

 

391

      PARAMETER (NSQAMPSO=1)

439

392

      INTEGER NSPLITORDERS

440

393

      PARAMETER (NSPLITORDERS=2)

441

394

C     

448

401

      INTEGER I,J

449

402

      INTEGER SQSPLITORDERS(NSQAMPSO,NSPLITORDERS)

450

403

C     the values listed below are for QCD, QED

451

 

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

452

 

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

453

 

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

 

404

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

454

405

C     

455

406

C     BEGIN CODE

456

407

C     

480

431

C     

481

432

      IMPLICIT NONE

482

433

      INTEGER NSQAMPSO

483

 

      PARAMETER (NSQAMPSO=3)

 

434

      PARAMETER (NSQAMPSO=1)

484

435

      INTEGER NSPLITORDERS

485

436

      PARAMETER (NSPLITORDERS=2)

486

437

C     

493

444

      INTEGER I

494

445

      INTEGER SQSPLITORDERS(NSQAMPSO,NSPLITORDERS)

495

446

C     the values listed below are for QCD, QED

496

 

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

497

 

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

498

 

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

 

447

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

499

448

C     

500

449

C     BEGIN CODE

501

450

C     

523

472

C     

524

473

      IMPLICIT NONE

525

474

      INTEGER NSQAMPSO

526

 

      PARAMETER (NSQAMPSO=3)

 

475

      PARAMETER (NSQAMPSO=1)

527

476

      INTEGER NSQSO

528

477

 

529

478

      NSQSO=NSQAMPSO

• Older »

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