pp → tH− + X

Asked by Basit Ali

I have to calculate cross section for process pp → tH− + X? what will be the batch file ?
also plz guide me for the following process.
 pp → t¯tb(¯b)+X?
I have implemented 2HDM-II-BSM and created batch file like this.
Model: 2HDM-II-BSM
Model changed: False
Gauge: Feynman
Process: p,p->t,H-,j

Composite: j=u,U,d,D,s,S,c,C,b,B,G
Composite: p=u,U,d,D,s,S,c,C,b,B,G
pdf1: PDT:cteq6l1(proton)
pdf2: PDT:cteq6l1(proton)
p1: 6500
p2: 6500

Run parameter: MH-
Run begin: 50
Run step size: 50
Run n steps: 3

Run parameter: ptcut
Run begin: 100
Run step size: 50
Run n steps: 3

#### running MH3 with the mass close to MHX ###
Parameter: MH3=MHX+1

Cut parameter: T(j)
Cut invert: False
Cut min: ptcut
Cut max:

Dist parameter: M(H-)
Dist min: 200
Dist max: 1000
Dist n bins: 100
Dist title: p,p->t,H-,j
Dist x-title: M(H-)(GeV)

Number of events (per run step): 1000
Filename: pp_HXH3j

nSess_1: 5
nCalls_1: 100000
nSess_2: 5
nCalls_2: 100000

Question information

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

Dear Basit,

your batch is absolutelyu correct.

It will work.

The process

pp → t¯tb+X

etc should be implemented in exct analigy.

Best regards,
Alexander

On 08/05/2020 21:04, Basit Ali wrote:
> New question #690561 on CalcHEP:
> https://answers.launchpad.net/calchep/+question/690561
>
> I have to calculate cross section for process pp → tH− + X? what will be the batch file ?
> also plz guide me for the following process.
> pp → t¯tb(¯b)+X?
> I have implemented 2HDM-II-BSM and created batch file like this.
> Model: 2HDM-II-BSM
> Model changed: False
> Gauge: Feynman
> Process: p,p->t,H-,j
>
> Composite: j=u,U,d,D,s,S,c,C,b,B,G
> Composite: p=u,U,d,D,s,S,c,C,b,B,G
> pdf1: PDT:cteq6l1(proton)
> pdf2: PDT:cteq6l1(proton)
> p1: 6500
> p2: 6500
>
>
> Run parameter: MH-
> Run begin: 50
> Run step size: 50
> Run n steps: 3
>
> Run parameter: ptcut
> Run begin: 100
> Run step size: 50
> Run n steps: 3
>
> #### running MH3 with the mass close to MHX ###
> Parameter: MH3=MHX+1
>
> Cut parameter: T(j)
> Cut invert: False
> Cut min: ptcut
> Cut max:
>
> Dist parameter: M(H-)
> Dist min: 200
> Dist max: 1000
> Dist n bins: 100
> Dist title: p,p->t,H-,j
> Dist x-title: M(H-)(GeV)
>
> Number of events (per run step): 1000
> Filename: pp_HXH3j
>
> nSess_1: 5
> nCalls_1: 100000
> nSess_2: 5
> nCalls_2: 100000
>
--
______________________________________________________________________
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
Basit Ali (ali-0514) said :
#2

But ,Sir I am facing a problem . It shows distributions along y-axis are empty.
How to set distributions along y-axis? otherwise graph goes parallel to x-axis.

Revision history for this message
Basit Ali (ali-0514) said :
#3

But ,Sir I am facing a problem . It shows distributions along y-axis are
empty.
How to set distributions along y-axis? otherwise graph goes parallel to
x-axis.

On Sat, May 9, 2020 at 2:09 AM Alexander Belyaev <
<email address hidden>> wrote:

> Your question #690561 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/690561
>
> Status: Open => Answered
>
> Alexander Belyaev proposed the following answer:
> Dear Basit,
>
> your batch is absolutelyu correct.
>
> It will work.
>
> The process
>
> pp → t¯tb+X
>
> etc should be implemented in exct analigy.
>
> Best regards,
> Alexander
>
>
> On 08/05/2020 21:04, Basit Ali wrote:
> > New question #690561 on CalcHEP:
> > https://answers.launchpad.net/calchep/+question/690561
> >
> > I have to calculate cross section for process pp → tH− + X? what will be
> the batch file ?
> > also plz guide me for the following process.
> > pp → t¯tb(¯b)+X?
> > I have implemented 2HDM-II-BSM and created batch file like this.
> > Model: 2HDM-II-BSM
> > Model changed: False
> > Gauge: Feynman
> > Process: p,p->t,H-,j
> >
> > Composite: j=u,U,d,D,s,S,c,C,b,B,G
> > Composite: p=u,U,d,D,s,S,c,C,b,B,G
> > pdf1: PDT:cteq6l1(proton)
> > pdf2: PDT:cteq6l1(proton)
> > p1: 6500
> > p2: 6500
> >
> >
> > Run parameter: MH-
> > Run begin: 50
> > Run step size: 50
> > Run n steps: 3
> >
> > Run parameter: ptcut
> > Run begin: 100
> > Run step size: 50
> > Run n steps: 3
> >
> > #### running MH3 with the mass close to MHX ###
> > Parameter: MH3=MHX+1
> >
> > Cut parameter: T(j)
> > Cut invert: False
> > Cut min: ptcut
> > Cut max:
> >
> > Dist parameter: M(H-)
> > Dist min: 200
> > Dist max: 1000
> > Dist n bins: 100
> > Dist title: p,p->t,H-,j
> > Dist x-title: M(H-)(GeV)
> >
> > Number of events (per run step): 1000
> > Filename: pp_HXH3j
> >
> > nSess_1: 5
> > nCalls_1: 100000
> > nSess_2: 5
> > nCalls_2: 100000
> >
> --
> ______________________________________________________________________
> 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
> ______________________________________________________________________
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
> https://answers.launchpad.net/calchep/+question/690561/+confirm?answer_id=0
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/calchep/+question/690561
>
> You received this question notification because you asked the question.
>

Revision history for this message
Alexander Belyaev (alexander.belyaev) said :
#4

It gives you what you are asking.

You have asked to plot  M(H-) which is a fixed number,

say 100,

and you want to put this number into histogram with the range

200--1000

so you have will have an empty histogram.

I.e. you have got what you have asked

Best

Alexander

On 08/05/2020 22:17, Basit Ali wrote:
> Question #690561 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/690561
>
> Basit Ali posted a new comment:
> But ,Sir I am facing a problem . It shows distributions along y-axis are empty.
> How to set distributions along y-axis? otherwise graph goes parallel to x-axis.
>
--
______________________________________________________________________
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
Basit Ali (ali-0514) said :
#5

Sir,
It shows empty y-range . adjusting to {-1:1}
I want to fix this problem .so that i will able to get proper plot

On Sat, May 9, 2020 at 3:33 AM Alexander Belyaev <
<email address hidden>> wrote:

> Your question #690561 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/690561
>
> Status: Open => Answered
>
> Alexander Belyaev proposed the following answer:
> It gives you what you are asking.
>
> You have asked to plot M(H-) which is a fixed number,
>
> say 100,
>
> and you want to put this number into histogram with the range
>
> 200--1000
>
> so you have will have an empty histogram.
>
> I.e. you have got what you have asked
>
> Best
>
> Alexander
>
>
>
> On 08/05/2020 22:17, Basit Ali wrote:
> > Question #690561 on CalcHEP changed:
> > https://answers.launchpad.net/calchep/+question/690561
> >
> > Basit Ali posted a new comment:
> > But ,Sir I am facing a problem . It shows distributions along y-axis are
> empty.
> > How to set distributions along y-axis? otherwise graph goes parallel to
> x-axis.
> >
> --
> ______________________________________________________________________
> 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
> ______________________________________________________________________
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
> https://answers.launchpad.net/calchep/+question/690561/+confirm?answer_id=3
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/calchep/+question/690561
>
> You received this question notification because you asked the question.
>

Revision history for this message
Best Alexander Belyaev (alexander.belyaev) said :
#6

Warning: empty y range [0:0], adjusting to [-1:1]
is a WARNING.

This is not and ERROR.

For 2D or higher dim scans  it does not plot cross sections

at

../calchep/work/html/numerical.html

and this si reflection of this

It does make a plot (distributions) for individual points of the scan

as soon as they defined correctly (which was your problem).

For example if you would define

Dist parameter:    M(H-)
Dist min:       0
Dist max:       1000
Dist n bins:       100
Dist title:       p,p->t,H-,j
Dist x-title:       M(H-)(GeV)

i.e.

Dist min:       200 --> Dist min: 0

then you would get non empty distribution

in

.../calchep/work/html/runs/MHc50ptcut100.html

etc for example

as I have explained you in my previous message

Is this clear now?

Best

Alexander

On 09/05/2020 07:34, Basit Ali wrote:
> Question #690561 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/690561
>
> Status: Answered => Open
>
> Basit Ali is still having a problem:
> Sir,
> It shows empty y-range . adjusting to {-1:1}
> I want to fix this problem .so that i will able to get proper plot
>
>
> On Sat, May 9, 2020 at 3:33 AM Alexander Belyaev <
> <email address hidden>> wrote:
>
>> Your question #690561 on CalcHEP changed:
>> https://answers.launchpad.net/calchep/+question/690561
>>
>> Status: Open => Answered
>>
>> Alexander Belyaev proposed the following answer:
>> It gives you what you are asking.
>>
>> You have asked to plot M(H-) which is a fixed number,
>>
>> say 100,
>>
>> and you want to put this number into histogram with the range
>>
>> 200--1000
>>
>> so you have will have an empty histogram.
>>
>> I.e. you have got what you have asked
>>
>> Best
>>
>> Alexander
>>
>>
>>
>> On 08/05/2020 22:17, Basit Ali wrote:
>>> Question #690561 on CalcHEP changed:
>>> https://answers.launchpad.net/calchep/+question/690561
>>>
>>> Basit Ali posted a new comment:
>>> But ,Sir I am facing a problem . It shows distributions along y-axis are
>> empty.
>>> How to set distributions along y-axis? otherwise graph goes parallel to
>> x-axis.
>> --
>> ______________________________________________________________________
>> 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
>> ______________________________________________________________________
>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>> https://answers.launchpad.net/calchep/+question/690561/+confirm?answer_id=3
>>
>> If you still need help, you can reply to this email or go to the
>> following page to enter your feedback:
>> https://answers.launchpad.net/calchep/+question/690561
>>
>> You received this question notification because you asked the question.
>>
--
______________________________________________________________________
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
Basit Ali (ali-0514) said :
#7

Thanks Alot.