../bin/subproc_cycle returns all zeros for cross sections of p,p->t,T

Asked by Suzanne Rosenzweig

I am working on a practice assignment for a course and we are asked to find the total cross section of p, p -> t, T (we take p to be composite of G, u, U, d, D for simplicity). However, even after setting the parton distribution function, the results of running the script ../bin/subproc_cycle yield all zeros for the individual and total cross sections. Any advice? Thank you!

Question information

Language:
English Edit question
Status:
Answered
For:
CalcHEP Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Alexander Belyaev (alexander.belyaev) said :
#1

Hello,

(sorry I do not see your name)

what is the energy of the collision did you put in?

Please cut and paste the content of one of your prt files which are
located in results folder.

Regards,

Alexander

On 14/10/2019 12:04, Suzanne Rosenzweig wrote:
> New question #685111 on CalcHEP:
> https://answers.launchpad.net/calchep/+question/685111
>
> I am working on a practice assignment for a course and we are asked to find the total cross section of p, p -> t, T (we take p to be composite of G, u, U, d, D for simplicity). However, even after setting the parton distribution function, the results of running the script ../bin/subproc_cycle yield all zeros for the individual and total cross sections. Any advice? Thank you!
>
--
______________________________________________________________________
Prof. Alexander S Belyaev (<email address hidden>)
https://www.hep.phys.soton.ac.uk/content/alexander-belyaev

School of Physics & Astronomy, University of Southampton
Office 5047, SO17 1BJ, TEL: +44 23805 98509, FAX: +44 23805 93910
.....................................................................
Particle Physics Department, Rutherford Appleton Laboratory
Didcot, OX11 0QX, TEL: +44 12354 45562, FAX: +44 12354 46733
.....................................................................
CERN, CH-1211 Geneva 23, Switzerland
Office 40/1-B20, Mailbox: E27910, TEL: +41 2276 71642
______________________________________________________________________

Revision history for this message
Suzanne Rosenzweig (srosenzweig) said :
#2

Thank you for answering! My name is Suzanne. The CM energy is 13 TeV (so p=6500 GeV). Here is one of the prt files:

> cat prt_12

#Subprocess 2 ( u, U -> t, T )
#Session_number 12
#Initial_state inP1=6.500000E+03 inP2=6.500000E+03
 Polarizations= { 0.000000E+00 0.000000E+00 }
  StrFun1="PDT:cteq6l1(proton)" 2212
  StrFun2="PDT:cteq6l1(proton)" 2212

#Physical_Parameters
        EE = 3.133300000000000E-01
        SW = 4.740000000000000E-01
         Q = 1.000000000000000E+02
        MW = 8.038500000000001E+01
       Mtp = 1.725000000000000E+02
      McMc = 1.230000000000000E+00
      MbMb = 4.250000000000000E+00
  alphaSMZ = 1.184000000000000E-01
        Ml = 1.777000000000000E+00
        Mh = 1.250000000000000E+02
----

#Breit-Wigner BW range 2.700000
t-channel widths 0
GI trick in s- 0
GI trick in t- 1
#VVdecays 1

#alphaQCD alphaPDF=1 alpha(MZ)=1.184000E-01 NF=5 Order=2 MbMb=4.200000E+00 Mtp=1.730000E+02
#QCDscales
 Renormalization M1
 Factorization1 Qren
 Factorization2 Qpdf1
 Showering Qren

#Composites
*** Table ***
 Composites
  Name |> Comma separated list of particles <|
p* |G,d,D,u,U
=========================================================================

#Cuts
*** Table ***
 Cuts
!| Parameter |> Min bound <|> Max bound <|
%|T(p*) | 50 |
%|J(p*,p*) | 0.5 |
============================================

#END

#Subprocess 2 ( u, U -> t, T )
#Session_number 12
#Initial_state inP1=6.500000E+03 inP2=6.500000E+03
 Polarizations= { 0.000000E+00 0.000000E+00 }
  StrFun1="PDT:cteq6l1(proton)" 2212
  StrFun2="PDT:cteq6l1(proton)" 2212

#Physical_Parameters
        EE = 3.133300000000000E-01
        SW = 4.740000000000000E-01
         Q = 1.000000000000000E+02
        MW = 8.038500000000001E+01
       Mtp = 1.725000000000000E+02
      McMc = 1.230000000000000E+00
      MbMb = 4.250000000000000E+00
  alphaSMZ = 1.184000000000000E-01
        Ml = 1.777000000000000E+00
        Mh = 1.250000000000000E+02
----

#Breit-Wigner BW range 2.700000
t-channel widths 0
GI trick in s- 0
GI trick in t- 1
#VVdecays 1

#alphaQCD alphaPDF=1 alpha(MZ)=1.184000E-01 NF=5 Order=2 MbMb=4.200000E+00 Mtp=1.730000E+02
#QCDscales
 Renormalization M1
 Factorization1 Qren
 Factorization2 Qpdf1
 Showering Qren

#Composites
*** Table ***
 Composites
  Name |> Comma separated list of particles <|
p* |G,d,D,u,U
=========================================================================

#Cuts
*** Table ***
 Cuts
!| Parameter |> Min bound <|> Max bound <|
%|T(p*) | 50 |
%|J(p*,p*) | 0.5 |
============================================

#END
   1 0.0000E+00 0.00E+00 98568
   2 0.0000E+00 0.00E+00 98568
   3 0.0000E+00 0.00E+00 98568
   4 0.0000E+00 0.00E+00 98568
   5 0.0000E+00 0.00E+00 98568
 < > 0.0000E+00 0.00E+00 492840 0
   1 0.0000E+00 0.00E+00 98568
   2 0.0000E+00 0.00E+00 98568
   3 0.0000E+00 0.00E+00 98568
   4 0.0000E+00 0.00E+00 98568
   5 0.0000E+00 0.00E+00 98568
 < > 0.0000E+00 0.00E+00 492840 0

Can you help with this problem?

Provide an answer of your own, or ask Suzanne Rosenzweig for more information if necessary.

To post a message you must log in.