wiki:WMassMeasurmentExample

Version 1 (modified by trac, 6 years ago) (diff)

--

Example of application/running of the code

Prepare the directory

  1. Duplication of Template directory $>cp -r Template First_test
  2. Change diretory: $> cd First_test
  3. Switch to MadWeight mode: $> ./bin/PassToMadWeight
  4. Parametrize the proc_card.dat. For this example the decay chain is the following one: pp>(W+>e+ve) and the transfer functions is dbl_gauss_pt_jet (the one foreseen for PGS in the CMS-like configuration). You can download this proc_card here.
  5. Generate the process: $> ./bin/newprocess
  6. Put the events sample in the ./Events directory under the name input.lhco. Our events sample is generated with MadGraph/Pythia/PGS: and is available here.

Prepare the run

  1. Configure the MadWeight_card.dat. In the given card, we will run in single machine for 20 events and 5 different mass hyppothesis. As they are no divergencies in the square matrix element, we don't want to put any cut in the evaluation of the integrals. So that we have desactivated the cut.
  2. Check the param_card.dat Check that you parameter not modified via the MadWeight_card are coherent. In this case we use default file.
  3. Check the run_card.dat Choose the name of the run (here we use "fermi") , the energy of the beam,... Note that the cuts defined here will not be used (as we have asked in the MadWeight_card.dat). Default files is used
  4. Define/Check transfer_card.dat. Only the jet part are often consider as process dependant, as far as we don't have jet, the defaulf parametrization (delta function for the electron) should be sufficient.

Launch MadWeight

In this example we will launch MW by step in order to explain in detail each step.

Note the log are done with version: 2.5.4

  1. TracMath macro processor has detected an error. Please fix the problem before continuing.


    The command:

    '/usr/bin/pdflatex -interaction=nonstopmode 86d7aebdf94d6ea05c0fb41c96e25ea3b0f9ce59.tex'
    failed with the following output:
    "This is pdfTeXk, Version 3.141592-1.40.3 (Web2C 7.5.6)\n %&-line parsing enabled.\nentering extended mode\n(./86d7aebdf94d6ea05c0fb41c96e25ea3b0f9ce59.tex\nLaTeX2e <2005/12/01>\nBabel <v3.8h> and hyphenation patterns for english, usenglishmax, dumylang, noh\nyphenation, arabic, basque, bulgarian, coptic, welsh, czech, slovak, german, ng\nerman, danish, esperanto, spanish, catalan, galician, estonian, farsi, finnish,\n french, greek, monogreek, ancientgreek, croatian, hungarian, interlingua, ibyc\nus, indonesian, icelandic, italian, latin, mongolian, dutch, norsk, polish, por\ntuguese, pinyin, romanian, russian, slovenian, uppersorbian, serbian, swedish, \nturkish, ukenglish, ukrainian, loaded.\n(/usr/share/texmf/tex/latex/base/article.cls\nDocument Class: article 2005/09/16 v1.4f Standard LaTeX document class\n(/usr/share/texmf/tex/latex/base/size10.clo))\n(/usr/share/texmf/tex/latex/base/inputenc.sty\n(/usr/share/texmf/tex/latex/base/utf8.def\n(/usr/share/texmf/tex/latex/base/t1enc.dfu)\n(/usr/share/texmf/tex/latex/base/ot1enc.dfu)\n(/usr/share/texmf/tex/latex/base/omsenc.dfu)))\n(/usr/share/texmf/tex/latex/cmap/cmap.sty) (./type1ec.sty\n(/usr/share/texmf/tex/latex/base/t1cmr.fd))\n(/usr/share/texmf/tex/latex/base/fontenc.sty\n(/usr/share/texmf/tex/latex/base/t1enc.def)<<t1.cmap>>)\n(/usr/share/texmf/tex/latex/amsmath/amsmath.sty\nFor additional information on amsmath, use the `?' option.\n(/usr/share/texmf/tex/latex/amsmath/amstext.sty\n(/usr/share/texmf/tex/latex/amsmath/amsgen.sty))\n(/usr/share/texmf/tex/latex/amsmath/amsbsy.sty)\n(/usr/share/texmf/tex/latex/amsmath/amsopn.sty))\n(/usr/share/texmf/tex/latex/amscls/amsthm.sty)\n(/usr/share/texmf/tex/latex/amsfonts/amssymb.sty\n(/usr/share/texmf/tex/latex/amsfonts/amsfonts.sty))\n(/usr/share/texmf/tex/latex/tools/bm.sty)\n(/usr/share/texmf/tex/latex/preview/preview.sty\n(/usr/share/texmf/tex/latex/preview/prtightpage.def))\n(./86d7aebdf94d6ea05c0fb41c96e25ea3b0f9ce59.aux)\nPreview: Fontsize 10pt\nPreview: PDFoutput 1\n(/usr/share/texmf/tex/latex/amsfonts/umsa.fd)\n(/usr/share/texmf/tex/latex/amsfonts/umsb.fd)\n! Misplaced alignment tab character &.\n<recently read> &\n                 \nl.14 $&\n       gt;./bin/madweight.py -1 or $\nPreview: Tightpage -32891 -32891 32891 32891\n[1{/usr/share/texmf/fonts/map/pdftex/updmap/pdftex.map}]\n(./86d7aebdf94d6ea05c0fb41c96e25ea3b0f9ce59.aux) )\n(see the transcript file for additional information)</usr/share/texmf/fonts/typ\ne1/bluesky/cm/cmmi10.pfb></usr/share/texmf/fonts/type1/bluesky/cm/cmr10.pfb></u\nsr/share/texmf/fonts/type1/bluesky/cm/cmsy10.pfb>\nOutput written on 86d7aebdf94d6ea05c0fb41c96e25ea3b0f9ce59.pdf (1 page, 11866 b\nytes).\nTranscript written on 86d7aebdf94d6ea05c0fb41c96e25ea3b0f9ce59.log.\n"

    &gt;./bin/madweight.py param:<br /> Create the different param_card and files related to the different cards.<br /> output: {{{

no mapping card found starting the creation of the param_card define new mapping file define new mapping file we have created 5 param_card's create file ./Source/madweight_card.inc create file ./Source/MadWeight_File/Transfer_Fct/transfer_card.inc

}}}

  1. TracMath macro processor has detected an error. Please fix the problem before continuing.


    The command:

    '/usr/bin/pdflatex -interaction=nonstopmode e56c557c1cfadfd70e1768f95297ebb97049c084.tex'
    failed with the following output:
    "This is pdfTeXk, Version 3.141592-1.40.3 (Web2C 7.5.6)\n %&-line parsing enabled.\nentering extended mode\n(./e56c557c1cfadfd70e1768f95297ebb97049c084.tex\nLaTeX2e <2005/12/01>\nBabel <v3.8h> and hyphenation patterns for english, usenglishmax, dumylang, noh\nyphenation, arabic, basque, bulgarian, coptic, welsh, czech, slovak, german, ng\nerman, danish, esperanto, spanish, catalan, galician, estonian, farsi, finnish,\n french, greek, monogreek, ancientgreek, croatian, hungarian, interlingua, ibyc\nus, indonesian, icelandic, italian, latin, mongolian, dutch, norsk, polish, por\ntuguese, pinyin, romanian, russian, slovenian, uppersorbian, serbian, swedish, \nturkish, ukenglish, ukrainian, loaded.\n(/usr/share/texmf/tex/latex/base/article.cls\nDocument Class: article 2005/09/16 v1.4f Standard LaTeX document class\n(/usr/share/texmf/tex/latex/base/size10.clo))\n(/usr/share/texmf/tex/latex/base/inputenc.sty\n(/usr/share/texmf/tex/latex/base/utf8.def\n(/usr/share/texmf/tex/latex/base/t1enc.dfu)\n(/usr/share/texmf/tex/latex/base/ot1enc.dfu)\n(/usr/share/texmf/tex/latex/base/omsenc.dfu)))\n(/usr/share/texmf/tex/latex/cmap/cmap.sty) (./type1ec.sty\n(/usr/share/texmf/tex/latex/base/t1cmr.fd))\n(/usr/share/texmf/tex/latex/base/fontenc.sty\n(/usr/share/texmf/tex/latex/base/t1enc.def)<<t1.cmap>>)\n(/usr/share/texmf/tex/latex/amsmath/amsmath.sty\nFor additional information on amsmath, use the `?' option.\n(/usr/share/texmf/tex/latex/amsmath/amstext.sty\n(/usr/share/texmf/tex/latex/amsmath/amsgen.sty))\n(/usr/share/texmf/tex/latex/amsmath/amsbsy.sty)\n(/usr/share/texmf/tex/latex/amsmath/amsopn.sty))\n(/usr/share/texmf/tex/latex/amscls/amsthm.sty)\n(/usr/share/texmf/tex/latex/amsfonts/amssymb.sty\n(/usr/share/texmf/tex/latex/amsfonts/amsfonts.sty))\n(/usr/share/texmf/tex/latex/tools/bm.sty)\n(/usr/share/texmf/tex/latex/preview/preview.sty\n(/usr/share/texmf/tex/latex/preview/prtightpage.def))\n(./e56c557c1cfadfd70e1768f95297ebb97049c084.aux)\nPreview: Fontsize 10pt\nPreview: PDFoutput 1\n(/usr/share/texmf/tex/latex/amsfonts/umsa.fd)\n(/usr/share/texmf/tex/latex/amsfonts/umsb.fd)\n! Misplaced alignment tab character &.\n<recently read> &\n                 \nl.14 $&\n       gt; ./bin/madweight.py -2 or $\nPreview: Tightpage -32891 -32891 32891 32891\n[1{/usr/share/texmf/fonts/map/pdftex/updmap/pdftex.map}]\n(./e56c557c1cfadfd70e1768f95297ebb97049c084.aux) )\n(see the transcript file for additional information)</usr/share/texmf/fonts/typ\ne1/bluesky/cm/cmmi10.pfb></usr/share/texmf/fonts/type1/bluesky/cm/cmr10.pfb></u\nsr/share/texmf/fonts/type1/bluesky/cm/cmsy10.pfb>\nOutput written on e56c557c1cfadfd70e1768f95297ebb97049c084.pdf (1 page, 11909 b\nytes).\nTranscript written on e56c557c1cfadfd70e1768f95297ebb97049c084.log.\n"

    &gt; ./bin/madweight.py analyzer:<br /> Analyzes the Feynman diagrams, the transfer function and creates the fortran code for the integration.<br /> output, indicates how MW analyzes your decay chains. and which change of variable it will apply on it. {{{ treating MW_P1_uxd_e-vex directory

structure of the configuration 1: particle: 3 pid : 11 level: 1 mother: -1 twin: 4 particle: 4 pid : -12 level: 1 mother: -1 twin: 3 particle: -1 pid : -24 level: 1 channel: S des: 4 3 mass/width: 40.0/2.04759951 1 ECS('s) 1 propagator(s) 1 missing particles(s) detail : b(1) 1 blob(s) associated Blob details: main 3

blob generation: [0
3 :]

treating MW_P1_dux_e-vex directory structure of the configuration 1: particle: 3 pid : 11 level: 1 mother: -1 twin: 4 particle: 4 pid : -12 level: 1 mother: -1 twin: 3 particle: -1 pid : -24 level: 1 channel: S des: 4 3 mass/width: 40.0/2.04759951 1 ECS('s) 1 propagator(s) 1 missing particles(s) detail : b(1) 1 blob(s) associated Blob details: main 3

blob generation: [0
3 :]

}}}

  1. TracMath macro processor has detected an error. Please fix the problem before continuing.


    The command:

    '/usr/bin/pdflatex -interaction=nonstopmode 9fd3a85ae4d350d82d80a4a6b9a76b6419e0fcf4.tex'
    failed with the following output:
    "This is pdfTeXk, Version 3.141592-1.40.3 (Web2C 7.5.6)\n %&-line parsing enabled.\nentering extended mode\n(./9fd3a85ae4d350d82d80a4a6b9a76b6419e0fcf4.tex\nLaTeX2e <2005/12/01>\nBabel <v3.8h> and hyphenation patterns for english, usenglishmax, dumylang, noh\nyphenation, arabic, basque, bulgarian, coptic, welsh, czech, slovak, german, ng\nerman, danish, esperanto, spanish, catalan, galician, estonian, farsi, finnish,\n french, greek, monogreek, ancientgreek, croatian, hungarian, interlingua, ibyc\nus, indonesian, icelandic, italian, latin, mongolian, dutch, norsk, polish, por\ntuguese, pinyin, romanian, russian, slovenian, uppersorbian, serbian, swedish, \nturkish, ukenglish, ukrainian, loaded.\n(/usr/share/texmf/tex/latex/base/article.cls\nDocument Class: article 2005/09/16 v1.4f Standard LaTeX document class\n(/usr/share/texmf/tex/latex/base/size10.clo))\n(/usr/share/texmf/tex/latex/base/inputenc.sty\n(/usr/share/texmf/tex/latex/base/utf8.def\n(/usr/share/texmf/tex/latex/base/t1enc.dfu)\n(/usr/share/texmf/tex/latex/base/ot1enc.dfu)\n(/usr/share/texmf/tex/latex/base/omsenc.dfu)))\n(/usr/share/texmf/tex/latex/cmap/cmap.sty) (./type1ec.sty\n(/usr/share/texmf/tex/latex/base/t1cmr.fd))\n(/usr/share/texmf/tex/latex/base/fontenc.sty\n(/usr/share/texmf/tex/latex/base/t1enc.def)<<t1.cmap>>)\n(/usr/share/texmf/tex/latex/amsmath/amsmath.sty\nFor additional information on amsmath, use the `?' option.\n(/usr/share/texmf/tex/latex/amsmath/amstext.sty\n(/usr/share/texmf/tex/latex/amsmath/amsgen.sty))\n(/usr/share/texmf/tex/latex/amsmath/amsbsy.sty)\n(/usr/share/texmf/tex/latex/amsmath/amsopn.sty))\n(/usr/share/texmf/tex/latex/amscls/amsthm.sty)\n(/usr/share/texmf/tex/latex/amsfonts/amssymb.sty\n(/usr/share/texmf/tex/latex/amsfonts/amsfonts.sty))\n(/usr/share/texmf/tex/latex/tools/bm.sty)\n(/usr/share/texmf/tex/latex/preview/preview.sty\n(/usr/share/texmf/tex/latex/preview/prtightpage.def))\n(./9fd3a85ae4d350d82d80a4a6b9a76b6419e0fcf4.aux)\nPreview: Fontsize 10pt\nPreview: PDFoutput 1\n(/usr/share/texmf/tex/latex/amsfonts/umsa.fd)\n(/usr/share/texmf/tex/latex/amsfonts/umsb.fd)\n! Misplaced alignment tab character &.\n<recently read> &\n                 \nl.14 $&\n       gt; ./bin/madweight.py -3 or $\nPreview: Tightpage -32891 -32891 32891 32891\n[1{/usr/share/texmf/fonts/map/pdftex/updmap/pdftex.map}]\n(./9fd3a85ae4d350d82d80a4a6b9a76b6419e0fcf4.aux) )\n(see the transcript file for additional information)</usr/share/texmf/fonts/typ\ne1/bluesky/cm/cmmi10.pfb></usr/share/texmf/fonts/type1/bluesky/cm/cmr10.pfb></u\nsr/share/texmf/fonts/type1/bluesky/cm/cmsy10.pfb>\nOutput written on 9fd3a85ae4d350d82d80a4a6b9a76b6419e0fcf4.pdf (1 page, 12216 b\nytes).\nTranscript written on 9fd3a85ae4d350d82d80a4a6b9a76b6419e0fcf4.log.\n"

    &gt; ./bin/madweight.py compilation: Compilation of the code. (Long and pointless output)
  2. $&gt; ./bin/madweight.py -4 or ./bin/madweight.py event: Check if the events in the input.lhco files are coherent for each subprocesses. {{{ name : fermi

time init Lhco_filter 0.0028760433197 time begin verif event Lhco_filter 0.0147750377655 time verif event Lhco_filter 2.26735496521 4617 selected events for ./SubProcesses/MW_P1_uxd_e-vex subprocess time Lhco_filter 2.27038884163 time init Lhco_filter 0.00105285644531 time begin verif event Lhco_filter 0.0399279594421 time verif event Lhco_filter 2.30534601212 4617 selected events for ./SubProcesses/MW_P1_dux_e-vex subprocess time Lhco_filter 2.30747294426 schedullar done

}}}

  1. $&gt; ./bin/madweight.py -5 or ./bin/madweight.py dir: Creates the directory for each parralel run. The name of the directory is for cross-section computation: ./SubProcesses/P1_uxd_e-vex/fermi/card_1/ for weight computation:./SubProcesses/MW_P1_uxd_e-vex/fermi/card_1/event_0/ {{{ name : fermi schedullar creating all directories P1_uxd_e-vex : deleting old run directory /data/0/O/omatt/MG_ME/TEST_W
create_dir: 100%
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Time: 00:00:00

created 5 directories P1_dux_e-vex : deleting old run directory /data/0/O/omatt/MG_ME/TEST_W

create_dir: 100%
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Time: 00:00:00

created 5 directories MW_P1_dux_e-vex : deleting old run directory /data/0/O/omatt/MG_ME/TEST_W mkdir: cannot create directory `./SubProcesses/MW_P1_dux_e-vex/fermi/': File exists

create_dir: 100%
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Time: 00:00:00

created 100 directories MW_P1_uxd_e-vex : deleting old run directory /data/0/O/omatt/MG_ME/TEST_W mkdir: cannot create directory `./SubProcesses/MW_P1_uxd_e-vex/fermi/': File exists

create_dir: 100%
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Time: 00:00:00

created 100 directories done

}}}

  1. $&gt; ./bin/madweight.py -6 or ./bin/madweight.py launch: Launches the computation of the weights on the cluster. (here by default on a single machine. {{{

mkdir: cannot create directory `./Events/fermi': File exists

compiling done compiling done schedullar launch: create scheduling file

donereation/Submission: 0%
ETA: --:--:--

make: `f77_madweight' is up to date. make: `f77_madweight' is up to date. make: `madevent' is up to date. make: `madevent' is up to date. launch ../../../f77_madweight ...(lot of launch ../../../comp_madweight)... launch ../../madevent.py

}}}

  1. $&gt; ./bin/madweight.py -7 or ./bin/madweight.py control: control the cluster status. This is cluster related, so ntohing happen in this case. {{{

name : fermi schedullar done

}}}

  1. $&gt; ./bin/madweight.py -8 or ./bin/madweight.py collect: collect all the different value {{{

name : fermi schedullar done collecting data collecting cross section results collect suceeded in 5 / 5 directories collect suceeded in 5 / 5 directories collecting output weight results collect in process MW_P1_dux_e-vex suceeded in 100 / 100 directories collect in process MW_P1_uxd_e-vex suceeded in 100 / 100 directories /data/0/O/omatt/MG_ME/TEST_MW2 /data/0/O/omatt/MG_ME/TEST_MW2 }}}

  1. $&gt; ./bin/madweight.py -9 or ./bin/madweight.py plot: compute the likelihood and try to make plot (here I don't have gnuplot on my cluster, so it fails) {{{

charge file ./Events/fermi/fermi_weights.out load cross section term ./Events/fermi/fermi_cross_weights.out no acceptance term loaded

}}} At this stage all the output are in ./Events/fermi/ For example the computation of the likelihood (whitout the AcceptanceTerm) gives: {{{

  1. 40 228.955053585 0.000715345946119
  2. 60 188.48197688 0.000750658557815
  3. 80 152.390896814 0.000965491294664
  4. 100 159.014840519 0.00113245384668
  5. 120 168.413249969 0.00123941979677

}}} The first column is the card_number, the second the W mass, the third the value of the likelihood, and the last one the error associate. In this simple example, we can claim that the correct mass is around 80 GeV. (In this case, the error is dominated by the binning)

-- Main.OlivierMattelaer - 14 Sep 2009

Attachments