Reweight VBS for aQGC

Hi,

This is a problem which only occurs with decay chain syntax and when some particles in the production are also present in the decay part since in that case
You have an ambiguity on how to evaluate the matrix-element if you only read the final state particles (which is what I do for the re-weighting).

If your Z decay is always Standard Model, then my suggestion would be to make such decay inside MadSpin, since Madspin is called after the reweighing.

Cheers,

Olivier

> On 14 Sep 2017, at 20:11, Tatsumi NITTA <email address hidden> wrote:
>
> Question #658027 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/658027
>
> Description changed to:
> Hello MG5 experts,
>
> We are trying to run MG5 internal reweight for pp>ZZjj>llqqjj process
> with anomalous gauge coupling implemented by EFT
> (http://feynrules.irmp.ucl.ac.be/attachment/wiki/AnomalousGaugeCoupling/all.tgz)
>
> It seems success for pp>ZZjj (Z is not decay), but not success for pp>ZZjj>llqqjj.
> Reweighting scheme was started, but seems skip "REWEIGHT: Event nb ???" process,
> finally does not get any reweighting results.
> Whole MG5 event generation process itself is finish as usual.
>
>
> - Successfully run without reweight generation.
> - We just add reweight_card.dat in Cards directory.
>
> Using MG5 v2.4.3 via atlas framework (we also checked standalone MG5
> v2.5.5 output is exactly same)
>
> I would very appreciate if you could any suggestion or advice to success
> running reweight for pp>ZZjj>llqqjj final state.
>
> Best regards,
> Tatsumi
>
> [proc_card_mg5.dat]
>
> import model sm-no_b_mass
> define p = g u c d s b u~ c~ d~ s~ b~
> define j = g u c d s b u~ c~ d~ s~ b~
> define l+ = e+ mu+ ta+
> define l- = e- mu- ta-
> define vl = ve vm vt
> define vl~ = ve~ vm~ vt~
> define w = w+ w-
> import model SM_LS0_LS1_LS2_UFO
> generate p p > j j z z QCD=0 QED=4 NP=1, z > l+ l-, z > j j
> output -f
>
>
> [Part of reweight_card.dat]
> launch --rwgt_name=FS0_1em12
> set anoinputs 1 1e-12
> set anoinputs 2 0e-12
> set anoinputs 3 0e-12
>
>
> [stdout]
>
> ===== for pp>ZZjj>llqqjj (not success) ===========
> 15:22:48 reweight -from_cards
> 16:09:45 REWEIGHT: load configuration from /tmp/local/scratch/tnitta/output_20170903_01/999998/0/PROC_SM_LS0_LS1_LS2_UFO_0/Cards/me5_configuration.txt
> 16:09:45 REWEIGHT: load configuration from /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.4.3.atlas/x86_64-slc6-gcc47-opt/input/mg5_configuration.txt
> 16:09:45 REWEIGHT: load configuration from /tmp/local/scratch/tnitta/output_20170903_01/999998/0/PROC_SM_LS0_LS1_LS2_UFO_0/Cards/me5_configuration.txt
> 16:09:45 REWEIGHT: Extracting the banner ...
> 16:09:45 REWEIGHT: process: p p > j j z z QCD=0 QED=4 NP=1, z > l+ l-, z > j j
> 16:09:45 REWEIGHT: options:
> 16:09:45 REWEIGHT: Running Reweighting
> 16:09:45 REWEIGHT: detected model: SM_LS0_LS1_LS2_UFO. Loading...
> 16:09:45 REWEIGHT: generating the square matrix element for reweighting
> 16:09:45 REWEIGHT: generate p p > j j z z QCD=0 QED=4 NP=1, z > l+ l-, z > j j ;
> 16:09:45 REWEIGHT: Done 2765
> 16:09:45 ^[[1;31mCRITICAL: tag, Pdir, data['id2path'][tag][1] = ^[[0m ((-2, 2), (-11, -4, -2, 2, 4, 11)) /tmp/local/scratch/tnitta/output_20170903_01/999998/0/PROC_SM\
> _LS0_LS1_LS2_UFO_0/rw_me/SubProcesses/P1_uux_ccxzz_z_epem_z_uux /tmp/local/scratch/tnitta/output_20170903_01/999998/0/PROC_SM_LS0_LS1_LS2_UFO_0/rw_me/SubProcesses/P1_uux_uuxzz_z_epem_z_ccx ^[[1;30m[reweight_interface.py at line 1405]^[[0m ^[[0m
> 16:09:45 REWEIGHT: Original cross-section: 0.0006401657 +- 3.602229e-06 pb
> 16:09:45 REWEIGHT: Computed cross-section:
> 16:09:45 REWEIGHT: storring files of previous run
> 16:09:45 REWEIGHT: Done
> 16:09:45 REWEIGHT:
> 16:09:45 REWEIGHT:
> 16:09:45 decay_events -from_cards
> 16:09:45 INFO: storring files of previous run
> 16:09:45 INFO: Done
>
> ===== for pp>ZZjj (success) ===========
> 08:29:49 reweight -from_cards
> 08:41:46 REWEIGHT: load configuration from /tmp/local/scratch/tnitta/output_20170903_01/999997/0/PROC_SM_LS0_LS1_LS2_UFO_0/Cards/me5_configuration.txt
> 08:41:46 REWEIGHT: load configuration from /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.4.3\
> .atlas/x86_64-slc6-gcc47-opt/input/mg5_configuration.txt
> 08:41:46 REWEIGHT: load configuration from /tmp/local/scratch/tnitta/output_20170903_01/999997/0/PROC_SM_LS0_LS1_LS2_UFO_0/Cards/me5_configuration.txt
> 08:41:46 REWEIGHT: Extracting the banner ...
> 08:41:46 REWEIGHT: process: p p > j j z z QCD=0 QED=4
> 08:41:46 REWEIGHT: options: NP=1
> 08:41:46 REWEIGHT: Running Reweighting
> 08:41:46 REWEIGHT: detected model: SM_LS0_LS1_LS2_UFO. Loading...
> 08:41:46 REWEIGHT: generating the square matrix element for reweighting
> 08:41:46 REWEIGHT: generate p p > j j z z QCD=0 QED=4 NP=1 ;
> 08:41:46 REWEIGHT: Done 82.39
> 08:41:46 >REWEIGHT: starts to compute weight for events with the following modification to the param_card:
> 08:41:46 REWEIGHT: set param_card anoinputs 1 1e-12 # orig: 1e-11
> 08:41:46
> 08:41:46 REWEIGHT: Event nb 0 current time: 08h31
> 08:41:46 REWEIGHT: Event nb 10 30.3s
> 08:41:46 REWEIGHT: Event nb 20 58.3s
> .... (just increasing nb)
> 08:41:46 REWEIGHT: Event nb 20000 9m 44s
> 08:41:46 REWEIGHT: All event done (nb_event: 22000) 9m 52s
> 08:41:46 REWEIGHT: Event /tmp/local/scratch/tnitta/output_20170903_01/999997/0/PROC_SM_LS0_LS1_LS2_UFO_0/Events/aQGCZZllqqjj/unweighted_events.lhe have now the additional weight
> ..... (other reweighting was running)
> 08:47:00 REWEIGHT: new cross-section is : 0.474127 pb (indicative error: 3.88879 pb)
> 08:47:00 REWEIGHT: Original cross-section: 0.1278308 +- 0.000190918 pb (cross-section from sum of weights: 0.1278307668)
> 08:47:00 REWEIGHT: Computed cross-section:
> 08:47:00 REWEIGHT: FS0_100em12 : 0.474127374395 +- 3.8887855492 pb
> 08:47:00 REWEIGHT: FS0_10em12 : 0.1278307668 +- 0.000190918 pb
> 08:47:00 REWEIGHT: FS0_1em12 : 0.129502029186 +- 0.000629392776534 pb
> 08:47:00 REWEIGHT: FS0_50em12 : 0.200256816538 +- 0.153826482808 pb
> 08:47:00 REWEIGHT: FS1_100em12 : 0.829225220773 +- 35.2995620702 pb
> 08:47:00 REWEIGHT: FS1_10em12 : 0.129504123001 +- 0.000881874133513 pb
> 08:47:00 REWEIGHT: FS1_1em12 : 0.129331000522 +- 0.000609002771555 pb
> 08:47:00 REWEIGHT: FS1_50em12 : 0.283815660737 +- 1.85810115963 pb
> 08:47:00 REWEIGHT: FS2_100em12 : 0.474127374395 +- 3.8887855492 pb
> 08:47:00 REWEIGHT: FS2_10em12 : 0.1278307668 +- 0.000190918 pb
> 08:47:00 REWEIGHT: FS2_1em12 : 0.129502029186 +- 0.000629392776534 pb
> 08:47:00 REWEIGHT: FS2_50em12 : 0.200256816538 +- 0.153826482808 pb
> 08:47:00 REWEIGHT: quit rwgt
> 08:47:00 REWEIGHT: storring files of previous run
> 08:47:00 REWEIGHT: Done
> 08:47:00 REWEIGHT:
> 08:47:00 REWEIGHT:
> 08:47:00 decay_events -from_cards
> 08:47:00 INFO: storring files of previous run
> 08:47:00 INFO: Done
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.