~madteam/mg5amcnlo/series2.0 : revision 263

357

Main:timesAllowErrors = $ERR_NUM_PY8 ! Number of allowed errors

358

359

! 2) Settings related to output in init(), next(), amd stat().

360

Init:showChangedSettings = on ! Shows all non-default settings

361

360

Init:showChangedParticleData = off ! Shows all non-default particle settings

362

Next:numberCount = 100 ! print message every n events

363

Next:numberShowInfo = 2 ! print event information n times

364

Next:numberShowProcess = 1 ! print process record n times

365

361

Next:numberShowEvent = $MAXPR_PY8 ! print event record n times

366

362

Stat:showProcessLevel = on ! Process statistics

367

363

Stat:showErrors = on ! Error statistics

368

364

Check:epTolErr = 0.001 ! Momentum-conservation tolerance

369

365

370

366

! 3) Beam-parameter settings.

371

Beams:idA = $iA ! Beam identities

372

Beams:idB = $iB ! Beam identities

373

367

Beams:frameType = 4 ! LHE initialization

374

368

Beams:LHEF = events.lhe ! Input LHE file

375

369

401

395

if [ "$ISLHA5" -eq 1 ]

402

396

then

403

397

cat <<EOF >> $ifile

404

PDF:pSet = LHAPDF5:$PDFSETN

398

PDF:pSet = LHAPDF5:$PDFSETN ! PDF set used in showers and multiparton interactions

405

399

EOF

406

400

if [ $PDFGROUP = "LHAPDF" ]

407

401

then

420

414

elif [ "$ISLHA6" -eq 1 ]

421

415

then

422

416

cat <<EOF >> $ifile

423

PDF:pSet = LHAPDF6:$PDFSETN

417

PDF:pSet = LHAPDF6:$PDFSETN ! PDF set used in showers and multiparton interactions

424

418

EOF

425

419

fi

426

420

fi

432

426

EOF

433

427

fi

434

428

cat <<EOF >> $ifile

435

#PDF:pSet = 7 ! CTEQ6L

436

PDF:pSet = 13 ! NNPDF2.3 QCD+QED

429

PDF:pSet = 13 ! PDF set used in showers and multiparton interactions (13 = NNPDF2.3 QCD+QED)

437

430

EOF

438

431

fi

439

432

444

437

EOF

445

438

fi

446

439

447

cat <<EOF >> $ifile

448

ProcessLevel:all = on ! Generation

449

ProcessLevel:resonanceDecays = on ! Resonance decays

450

PartonLevel:all = on ! Paron level: if off, stops after hard process generation

451

EOF

452

440

if [ $UE_PY8 = ".FALSE." ]

453

441

then

454

442

cat <<EOF >> $ifile

455

PartonLevel:MPI = off ! Multiple interactions

443

PartonLevel:MPI = off ! Multiparton interactions

444

! ** USE on FOR REALISTIC SHOWER SIMULATION ** !!!

456

445

EOF

457

446

else

458

447

cat <<EOF >> $ifile

459

PartonLevel:MPI = on ! Multiple interactions

448

PartonLevel:MPI = on ! Multiparton interactions

460

449

EOF

461

450

fi

462

451

463

cat <<EOF >> $ifile

464

PartonLevel:ISR = on ! Initial state shower

465

PartonLevel:FSR = on ! Final state shower

466

PartonLevel:FSRinProcess = on ! Final state shower in association with the hard process

467

PartonLevel:FSRinResonances = on ! Final state shower in resonance decays

468

HadronLevel:all = on ! Hadron level: if off, stops before hadronization

469

EOF

470

471

452

if [ $HADRONIZE_PY8 = ".FALSE." ]

472

453

then

473

454

cat <<EOF >> $ifile

474

HadronLevel:Hadronize = off ! Hadronization

455

HadronLevel:Hadronize = off ! Hadronisation

456

! ** USE on FOR REALISTIC SHOWER SIMULATION ** !!!

475

457

EOF

476

458

else

477

459

cat <<EOF >> $ifile

478

HadronLevel:Hadronize = on ! Hadronization

460

HadronLevel:Hadronize = on ! Hadronisation

479

461

EOF

480

462

fi

481

463

482

464

cat <<EOF >> $ifile

483

#HadronLevel:Decay = on ! Hadron decays

484

PhaseSpace:mHatMin = 4. ! Min invariant mass

485

PhaseSpace:mHatMax = -1. ! Max invariant mass

486

PhaseSpace:pTHatMin = 0. ! Min pT in 2->2

487

PhaseSpace:pTHatMax = -1. ! Max pT in 2->2

488

PhaseSpace:pTHatMinDiverge = 1. ! If massless final state, to avoid divergences

489

PhaseSpace:useBreitWigners = on ! Masses according to Breit-Wigner

490

#PhaseSpace:pTHat3Min = 0. ! Min pT for the hardest parton in 2->3

491

PhaseSpace:pTHat3Max = -1. ! Max pT for the hardest parton in 2->3

492

PhaseSpace:pTHat5Min = 0. ! Min pT for the softest parton in 2->3

493

PhaseSpace:pTHat5Max = -1. ! Max pT for the softest parton in 2->3

494

PhaseSpace:RsepMin = 0. ! Min R separation in 2->3

495

465

496

466

! 5) Final-state shower.

497

467

TimeShower:pTmaxMatch = 1 ! Use scalup (re-check)

499

469

TimeShower:alphaSvalue = 0.118 ! Alpha_s(MZ) in final-state shower

500

470

TimeShower:alphaSorder = 1 ! Alpha_s running order in final-state shower

501

471

TimeShower:alphaEMorder = 0 ! Alpha_EM running order in final-state shower

502

TimeShower:interleave = on ! If on, FSR interleaved with ISR

503

TimeShower:allowBeamRecoil = on ! If off, no energy transfer to from ISR to FSR

504

472

TimeShower:dampenBeamRecoil = off ! Dampens the effect of beam recoil

505

473

TimeShower:globalRecoil = on ! All final-state particles recoil against the branching

506

474

TimeShower:nMaxGlobalRecoil = 1 ! Number of splittings with TimeShower:globalRecoil = on

508

476

TimeShower:nMaxGlobalBranch = 1 ! Number of FSR splittings proposed with global recoil

509

477

TimeShower:nPartonsInBorn = -1 ! Number of Born QCD final-state partons (to treat H and S differently)

510

478

TimeShower:limitPTmaxGlobal = on ! Limits pT < min(SCALUP,mDipole/2).

511

TimeShower:QCDshower = on ! QCD final-state shower

512

TimeShower:nGluonToQuark = 5 ! Number if flavors allowed in g->qqbar

479

513

480

TimeShower:QEDshowerByQ = off ! Prevent quarks from radiating photons

481

! ** USE on FOR REALISTIC SHOWER SIMULATION ** !!!

514

482

TimeShower:QEDshowerByL = off ! Prevent leptons from radiating photons

483

! ** USE on FOR REALISTIC SHOWER SIMULATION ** !!!

515

484

TimeShower:QEDshowerByGamma = off ! Prevent photons from branching

485

! ** USE on FOR REALISTIC SHOWER SIMULATION ** !!!

486

516

487

TimeShower:MEcorrections = off ! No Matrix-element corrections

517

TimeShower:MEafterFirst = off ! No Matrix-element corrections after first emission

518

TimeShower:phiPolAsym = on ! Azimuthal asymmetry induced by gluon polarization

519

488

TimeShower:alphaSuseCMW = false ! Use the CMW prescription in FSR

520

489

TimeShower:weightGluonToQuark = 1 ! Use normal Altarelli-Parisi kernels for g -> q qbar

521

490

525

494

SpaceShower:alphaSvalue = 0.118 ! Alpha_s(MZ) in final-state shower

526

495

SpaceShower:alphaSorder = 1 ! Alpha_s running order in final-state shower

527

496

SpaceShower:alphaEMorder = 0 ! Alpha_EM running order in final-state shower

528

SpaceShower:QCDshower = on ! QCD final-state shower

497

529

498

SpaceShower:QEDshowerByQ = off ! Prevent quarks from radiating photons

499

! ** USE on FOR REALISTIC SHOWER SIMULATION ** !!!

530

500

SpaceShower:QEDshowerByL = off ! Prevent leptons from radiating photons

501

! ** USE on FOR REALISTIC SHOWER SIMULATION ** !!!

502

531

503

SpaceShower:MEcorrections = off ! No Matrix-element corrections

532

SpaceShower:MEafterFirst = off ! No Matrix-element corrections after first emission

533

SpaceShower:phiPolAsym = on ! Azimuthal asymmetry induced by gluon polarization

534

SpaceShower:nQuarkIn = 5 ! Number of flavors in g->qqbar and also in incoming beams

535

504

SpaceShower:rapidityorder = off ! Do not order branchings in rapidity

536

505

SpaceShower:alphaSuseCMW = false ! Use the CMW prescription in ISR

537

506

538

507

! 7) Non-perturbative stuff.

539

508

BeamRemnants:primordialKT = off ! No primordial kT

509

! ** USE on FOR REALISTIC SHOWER SIMULATION ** !!!

540

510

541

! 8) Particle characteristics.

511

! 8) Particle features.

542

512

1:m0 = $DMASSPY ! down mass

543

513

2:m0 = $UMASSPY ! up mass

544

514

3:m0 = $SMASSPY ! strange mass

581

551

JetMatching:scheme = 1 ! MadGraph MLM-type merging

582

552

JetMatching:setMad = off ! Don't read from LHE in MG5LO format

583

553

584

585

554

EOF

586

555

fi

556

cat <<EOF >> $ifile

557

! 10) Decays.

558

EOF

587

559

588

560

if [ $PI_STABLE_PY8 = ".TRUE." ]

589

561

then

590

562

cat <<EOF >> $ifile

591

563

111:mayDecay = false ! stable pi0

592

#211:mayDecay = false ! stable pions

593

564

EOF

594

565

fi

595

566

if [ $B_STABLE_PY8 = ".TRUE." ]

659

630

EOF

660

631

fi

661

632

633

662

634

arr=()

663

635

numDM=0

664

636

iDM=0

671

643

fi

672

644

done

673

645

674

numDM=${#arr[@]}

675

for DM in "${arr[@]}"

676

do

677

DM=( $DM )

678

cat <<EOF >> ./$ifile

679

${DM[0]}:onMode = off ! First switch decays off

680

EOF

681

done

682

683

for DM in "${arr[@]}"

684

do

685

let "iDM+=1"

686

DM=( $DM )

687

num_of_part=$(echo ${#DM[@]})

688

let "num_of_part-=6"

689

690

if [ $num_of_part = 2 ]

691

then

692

cat <<EOF >> ./$ifile

693

${DM[0]}:onIfMatch = ${DM[2]} ${DM[3]} ! Decay mode $iDM of $numDM

694

EOF

695

elif [ $num_of_part = 3 ]

696

then

697

cat <<EOF >> ./$ifile

698

${DM[0]}:onifMatch = ${DM[2]} ${DM[3]} ${DM[4]} ! Decay mode $iDM of $numDM

699

EOF

700

elif [ $num_of_part = 4 ]

701

then

702

cat <<EOF >> ./$ifile

703

${DM[0]}:onifMatch = ${DM[2]} ${DM[3]} ${DM[4]} ${DM[5]} ! Decay mode $iDM of $numDM

704

EOF

705

elif [ $num_of_part = 5 ]

706

then

707

cat <<EOF >> ./$ifile

708

${DM[0]}:onifMatch = ${DM[2]} ${DM[3]} ${DM[4]} ${DM[5]} ${DM[6]} ! Decay mode $iDM of $numDM

709

EOF

710

fi

646

for DM in "${arr[@]}"

647

do

648

cat <<EOF >> ./$ifile

649

$DM ! decay mode

650

EOF

711

651

done

712

652

713

653

786

726

ExtraIncl="-I""$ExtraIncl"

787

727

fi

788

728

}

729

730

731

# Collection of Pythia8 settings that used to appear in the

732

# card, and that are supefluous because overridden or default

733

734

#Beams:idA = $iA ! Beam identities

735

#Beams:idB = $iB ! Beam identities

736

#ProcessLevel:all = on ! Generation

737

#ProcessLevel:resonanceDecays = on ! Resonance decays

738

#PartonLevel:all = on ! Paron level: if off, stops after hard

739

# ! process generation

740

#PartonLevel:ISR = on ! Initial state shower

741

#PartonLevel:FSR = on ! Final state shower

742

#PartonLevel:FSRinProcess = on ! Final state shower in association with

743

# ! the hard process

744

#PartonLevel:FSRinResonances = on ! Final state shower in resonance decays

745

#HadronLevel:all = on ! Hadron level: if off, stops before

746

# ! hadronization

747

#HadronLevel:Decay = on ! Hadron decays

748

#PhaseSpace:mHatMin = 4. ! Min invariant mass

749

#PhaseSpace:mHatMax = -1. ! Max invariant mass

750

#PhaseSpace:pTHatMin = 0. ! Min pT in 2->2

751

#PhaseSpace:pTHatMax = -1. ! Max pT in 2->2

752

#PhaseSpace:pTHatMinDiverge = 1. ! If massless final state, to avoid

753

# ! divergences

754

#PhaseSpace:useBreitWigners = on ! Masses according to Breit-Wigner

755

#PhaseSpace:pTHat3Min = 0. ! Min pT for the hardest parton in 2->3

756

#PhaseSpace:pTHat3Max = -1. ! Max pT for the hardest parton in 2->3

757

#PhaseSpace:pTHat5Min = 0. ! Min pT for the softest parton in 2->3

758

#PhaseSpace:pTHat5Max = -1. ! Max pT for the softest parton in 2->3

759

#PhaseSpace:RsepMin = 0. ! Min R separation in 2->3

760

#TimeShower:interleave = on ! If on, FSR interleaved with ISR

761

#TimeShower:allowBeamRecoil = on ! If off, no energy transfer to from

762

# ! ISR to FSR

763

#TimeShower:QCDshower = on ! QCD final-state shower

764

#TimeShower:nGluonToQuark = 5 ! Number if flavors allowed in g->qqbar

765

#TimeShower:phiPolAsym = on ! Azimuthal asymmetry induced by gluon

766

# ! polarization

767

#SpaceShower:QCDshower = on ! QCD final-state shower

768

#SpaceShower:phiPolAsym = on ! Azimuthal asymmetry induced by gluon

769

# ! polarization

770

#SpaceShower:nQuarkIn = 5 ! Number of flavors in g->qqbar and also in

771

# ! incoming beams

772

#Init:showChangedSettings = on ! Shows all non-default settings

773

#Next:numberCount = 100 ! print message every n events

774

#Next:numberShowInfo = 2 ! print event information n times

775

#Next:numberShowProcess = 1 ! print process record n times

776

#TimeShower:MEafterFirst = off ! No Matrix-element corrections after first emission

777

# ! it is already off if MEcorrections = off

778

#SpaceShower:MEafterFirst = off ! No Matrix-element corrections after first emission

779

# ! it is already off if MEcorrections = off

780