ParticleGun Problem

What happens, when, and how to reproduce the problem:

==> Run a Pythia6PtGun to generate Z0->mumu, with pT[0,20]GeV/c, in CMSSW_3_3_4

* The python file is this.

* turned off initial/final state radiation, and modified the pythiaDefault_cff.py also to include mumu decay only, like this.

* run this analyzer with this python file.

* In the analyzer, look at the generated Z0 pT distribution (left figure bellow), and then, loop over the daughter muons, pair them, and retrieve the Z0 parent (right figure bellow).


pt_genZ0.gif pt_genDimuon.gif

Basically, what it happens, see also the bellow printout, is that though the parent Z0 has a p!=0, the muons are decayed back-to-back, as if the Z0 is at rest (the boost part of the decay was not propagated properly). It happens only with the Z0 (tried with JPsi and don't get into this issue), and it's not a PYTHIA problem (I ran PYTHIA 6-422, and it never gets into this situation, so the decayer itself is fine).

pt_z0Gen = 12.602                    px = -7.59486   py = 10.0563    pz = 254.636 p = 254.94 E = 268.55 mass = 84.3823
E_mu1 = 184.443 pt1=33.5841 px1 = 6.15032  py1 = -33.0161 pz1 = 181.36   eta1 = 2.38804 phi1 = -1.38662
E_mu2 = 84.1065 pt2=45.1056 px2 = -13.6769 py2 = 42.982    pz2 = 70.9885 eta2 = 1.23503 phi2 = 1.87887
pt_z0Dimu = 12.4888                px = -7.52661   py = 9.96591    pz = 252.348 p = 252.657 E = 268.55 mass = 91.0124

pt_z0Gen = 18.302                    px = -13.6041   py = 12.2437    pz = -2.07556 p = 18.4198 E = 92.4163 mass = 90.5621
E_mu1 = 46.6812 pt1=40.4225 px1 = -11.1639 py1 = 38.8504  pz1 = -23.3484 eta1 = -0.549529 phi1 = 1.85061
E_mu2 = 46.6812 pt2=40.4225 px2 = 11.1639  py2 = -38.8504 pz2 = 23.3484 eta2 = 0.549529 phi2 = -1.29098
pt_z0Dimu = 0                           px = 0               py = 0               pz = 0              p = 0            E = 93.3625 mass = 93.3625

pt_z0Gen = 13.9235                  px = 12.453      py = -6.22797   pz = 42.1019  p = 44.3445 E = 84.7429 mass = 72.2144
E_mu1 = 46.0783 pt1=14.1707 px1 = 4.48367  py1 = -13.4426 pz1 = -43.845 eta1 = -1.8477 phi1 = -1.24886
E_mu2 = 46.0783 pt2=14.1707 px2 = -4.48367 py2 = 13.4426  pz2 = 43.845  eta2 = 1.8477 phi2 = 1.89273
pt_z0Dimu = 0                           px = 0                py = 0              pz = 0              p = 0            E = 92.1566 mass = 92.1566

Topic attachments
I Attachment History Action Size Date Who Comment
Texttxt DecayGunPtZ2MuMu_cfi.py.txt r1 manage 1.0 K 2009-11-26 - 18:33 CameliaMironov  
Unknown file formatcc Z0RecoMuonAnalyzer.cc r1 manage 7.2 K 2009-11-26 - 18:33 CameliaMironov  
Texttxt gunRunner.py.txt r1 manage 2.0 K 2009-11-26 - 18:33 CameliaMironov  
GIFgif pt_genDimuon.gif r1 manage 7.1 K 2009-11-26 - 18:32 CameliaMironov  
GIFgif pt_genZ0.gif r1 manage 10.2 K 2009-11-26 - 18:32 CameliaMironov  
Texttxt pythiaDefault_cff.py.txt r1 manage 1.8 K 2009-11-26 - 18:33 CameliaMironov  
Edit | Attach | Watch | Print version | History: r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r1 - 2009-11-26 - CameliaMironov
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback