~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_6.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_6.f

9

9

C      orders.inc (NLO_ORDERS)

10

10

C     

11

11

C     

12

 

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

13

 

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

 

12

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

 

13

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

113

 

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

 

112

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

 

113

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

168

168

      DATA GOODHEL/NCOMB*.FALSE./

169

169

      INTEGER NTRY

170

170

      DATA NTRY/0/

171

 

      DATA IDEN/96/

 

171

      DATA IDEN/12/

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

238

 

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

 

237

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

 

238

C     Process: c a > t t~ c [ 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=12, NCOLOR=4)

 

247

      PARAMETER (NWAVEFUNCS=8, 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,3,1) T(5,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,3,4) T(5,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(2,5,1) T(3,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(2,5,4) T(3,1)

 

279

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

 

280

     $ ,3.000000000000000D+00/

 

281

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

 

282

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

 

283

     $ ,9.000000000000000D+00/

 

284

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

295

285

C     ----------

296

286

C     BEGIN CODE

297

287

C     ----------

300

290

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

301

291

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

302

292

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

303

 

      CALL FFV1_2(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))

 

293

      CALL FFV1_2(W(1,1),W(1,2),GC_2,ZERO,ZERO,W(1,6))

 

294

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

305

295

C     Amplitude(s) for diagram number 1

306

 

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

307

 

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

 

296

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

 

297

      CALL FFV1P0_3(W(1,1),W(1,5),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))

308

299

C     Amplitude(s) for diagram number 2

309

 

      CALL FFV1_0(W(1,6),W(1,5),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))

 

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

311

302

C     Amplitude(s) for diagram number 3

312

 

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

313

 

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

314

 

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

 

303

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

 

304

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

315

305

C     Amplitude(s) for diagram number 4

316

 

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

317

 

      CALL FFV1P0_3(W(1,1),W(1,5),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_5_3(W(1,1),W(1,5),-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_1(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,1),W(1,11),W(1,7),GC_2,AMP(11))

335

 

C     Amplitude(s) for diagram number 12

336

 

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

337

 

C     Amplitude(s) for diagram number 13

338

 

      CALL FFV2_5_0(W(1,1),W(1,11),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(2)+(-5.000000000000000D

341

 

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

342

 

     $ *AMP(10)

343

 

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

344

 

     $ -01)*AMP(8)

345

 

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

346

 

     $ -01)*AMP(12)

347

 

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

348

 

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

349

 

     $ *AMP(12)

350

 

C     JAMPs contributing to orders QCD=1 QED=2

351

 

      JAMP(1,2) = 0D0

352

 

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

353

 

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

354

 

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

355

 

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

356

 

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

357

 

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

 

306

      CALL FFV1_0(W(1,1),W(1,8),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)

358

313

 

359

314

      DO I=0,NSQAMPSO

360

315

        RES(I)=0D0

396

351

C     

397

352

      IMPLICIT NONE

398

353

      INTEGER NAMPSO, NSQAMPSO

399

 

      PARAMETER (NAMPSO=2, NSQAMPSO=3)

 

354

      PARAMETER (NAMPSO=1, NSQAMPSO=1)

400

355

      INTEGER NSPLITORDERS

401

356

      PARAMETER (NSPLITORDERS=2)

402

357

C     

408

363

C     

409

364

      INTEGER I, SQORDERS(NSPLITORDERS)

410

365

      INTEGER AMPSPLITORDERS(NAMPSO,NSPLITORDERS)

411

 

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

412

 

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

 

366

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

413

367

C     

414

368

C     FUNCTION

415

369

C     

434

388

C     

435

389

      IMPLICIT NONE

436

390

      INTEGER NSQAMPSO

437

 

      PARAMETER (NSQAMPSO=3)

 

391

      PARAMETER (NSQAMPSO=1)

438

392

      INTEGER NSPLITORDERS

439

393

      PARAMETER (NSPLITORDERS=2)

440

394

C     

447

401

      INTEGER I,J

448

402

      INTEGER SQSPLITORDERS(NSQAMPSO,NSPLITORDERS)

449

403

C     the values listed below are for QCD, QED

450

 

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

451

 

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

452

 

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

 

404

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

453

405

C     

454

406

C     BEGIN CODE

455

407

C     

479

431

C     

480

432

      IMPLICIT NONE

481

433

      INTEGER NSQAMPSO

482

 

      PARAMETER (NSQAMPSO=3)

 

434

      PARAMETER (NSQAMPSO=1)

483

435

      INTEGER NSPLITORDERS

484

436

      PARAMETER (NSPLITORDERS=2)

485

437

C     

492

444

      INTEGER I

493

445

      INTEGER SQSPLITORDERS(NSQAMPSO,NSPLITORDERS)

494

446

C     the values listed below are for QCD, QED

495

 

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

496

 

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

497

 

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

 

447

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

498

448

C     

499

449

C     BEGIN CODE

500

450

C     

522

472

C     

523

473

      IMPLICIT NONE

524

474

      INTEGER NSQAMPSO

525

 

      PARAMETER (NSQAMPSO=3)

 

475

      PARAMETER (NSQAMPSO=1)

526

476

      INTEGER NSQSO

527

477

 

528

478

      NSQSO=NSQAMPSO

• Older »

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