creating lhe file for decay process
How to generate lhe file for a decay process, e.g. h-> b b in calchep?
I am able to generate lhe files for 2->2, 2->1 process, but not for 1->2 process.
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- CalcHEP Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- PARTHA KUMAR PAUL
- Solved:
- Last query:
- Last reply:
Revision history for this message
|
#1 |
Dear Partha Kumar,
I think events (.lhe folder) just can be generated for collision process just like you mentioned (2->2 and 2->1). However, it does not work for decays (1-> 2).
Sincerely,
Marco Arroyo.
Revision history for this message
|
#2 |
Sorry. I forget to answer.
Try the following way. First, create events in CalcHEP format using
CalcHEP menu system and then use
CalcHEP/
Let me now in case of problem.
Best
Alexander Pukhov
On 10/5/23 11:50, PARTHA KUMAR PAUL wrote:
> New question #708114 on CalcHEP:
> https:/
>
> How to generate lhe file for a decay process, e.g. h-> b b in calchep?
> I am able to generate lhe files for 2->2, 2->1 process, but not for 1->2 process.
>
>
Revision history for this message
|
#3 |
Thank you Prof. Pukhov, this is working. In the subproc_cycle I am giving
100 events, but I am getting only 1. I want to use it as a input in
pythia. Is it possible to generate 100 events for h->b B?
On Tue, Oct 17, 2023 at 2:15 AM Alexander Pukhov <
<email address hidden>> wrote:
> Your question #708114 on CalcHEP changed:
> https:/
>
> Alexander Pukhov proposed the following answer:
> Sorry. I forget to answer.
>
> Try the following way. First, create events in CalcHEP format using
> CalcHEP menu system and then use
>
> CalcHEP/
>
> Let me now in case of problem.
>
> Best
>
> Alexander Pukhov
>
>
> On 10/5/23 11:50, PARTHA KUMAR PAUL wrote:
> > New question #708114 on CalcHEP:
> > https:/
> >
> > How to generate lhe file for a decay process, e.g. h-> b b in calchep?
> > I am able to generate lhe files for 2->2, 2->1 process, but not for 1->2
> process.
> >
> >
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
> https:/
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https:/
>
> You received this question notification because you asked the question.
>
--
*Partha Kumar Paul*
Revision history for this message
|
#4 |
Thanks a lot for your mail.
It is indeed a bug caused by disagreement of changes. Namely here
n_calchep waits for natural number for requested number of events, but
subproc_cycle substitutes a floating point one. To fix the bug, please,
open file
CalcHEP/
The line 570 contains the corresponding instruction
case 10: correctInt(
It has to be replaced on
case 10:
{ double nEventsR=nEvents;
}
Best
Alexander Pukhov
Try the following. Set
using CalcHEP menus, and after that call
../..
Unfortunately subproc_cycle finishes with
"Freeze grid ON "
and for next call you will need to improve it via menus.
Best
Alexander Pukhov
On 10/20/23 20:35, PARTHA KUMAR PAUL wrote:
> Question #708114 on CalcHEP changed:
> https:/
>
> Status: Answered => Open
>
> PARTHA KUMAR PAUL is still having a problem:
> Thank you Prof. Pukhov, this is working. In the subproc_cycle I am giving
> 100 events, but I am getting only 1. I want to use it as a input in
> pythia. Is it possible to generate 100 events for h->b B?
>
> On Tue, Oct 17, 2023 at 2:15 AM Alexander Pukhov <
> <email address hidden>> wrote:
>
>> Your question #708114 on CalcHEP changed:
>> https:/
>>
>> Alexander Pukhov proposed the following answer:
>> Sorry. I forget to answer.
>>
>> Try the following way. First, create events in CalcHEP format using
>> CalcHEP menu system and then use
>>
>> CalcHEP/
>>
>> Let me now in case of problem.
>>
>> Best
>>
>> Alexander Pukhov
>>
>>
>> On 10/5/23 11:50, PARTHA KUMAR PAUL wrote:
>>> New question #708114 on CalcHEP:
>>> https:/
>>>
>>> How to generate lhe file for a decay process, e.g. h-> b b in calchep?
>>> I am able to generate lhe files for 2->2, 2->1 process, but not for 1->2
>> process.
>>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>> https:/
>>
>> If you still need help, you can reply to this email or go to the
>> following page to enter your feedback:
>> https:/
>>
>> You received this question notification because you asked the question.
>>
>
Revision history for this message
|
#5 |
Thank you for the reply. I am able to generate events in calcHEP format
now. Using event2lhe, I am able to convert that into lhe format (in
terminal). But how to save it in lhe format?
When I am generating events for h->2*x, I am getting multiple event files
corresponding to each decay mode. If I try to combine these .txt files into
one lhe using the event_mixer, the following error comes:
Total width of 25 is not found
There are no events
On Sun, Oct 22, 2023 at 3:41 PM Alexander Pukhov <
<email address hidden>> wrote:
> Your question #708114 on CalcHEP changed:
> https:/
>
> Status: Open => Answered
>
> Alexander Pukhov proposed the following answer:
> Thanks a lot for your mail.
>
> It is indeed a bug caused by disagreement of changes. Namely here
> n_calchep waits for natural number for requested number of events, but
> subproc_cycle substitutes a floating point one. To fix the bug, please,
> open file
>
> CalcHEP/
>
> The line 570 contains the corresponding instruction
>
> case 10: correctInt(
>
> It has to be replaced on
>
> case 10:
> { double nEventsR=nEvents;
> correctDouble(
> nEvents=nEventsR;
> break;
> }
>
> Best
>
> Alexander Pukhov
>
>
>
> Try the following. Set
>
> "Freeze grid OFF"
>
> using CalcHEP menus, and after that call
>
> ../../bin/
>
> Unfortunately subproc_cycle finishes with
>
> "Freeze grid ON "
>
> and for next call you will need to improve it via menus.
>
>
> Best
>
> Alexander Pukhov
>
>
> On 10/20/23 20:35, PARTHA KUMAR PAUL wrote:
> > Question #708114 on CalcHEP changed:
> > https:/
> >
> > Status: Answered => Open
> >
> > PARTHA KUMAR PAUL is still having a problem:
> > Thank you Prof. Pukhov, this is working. In the subproc_cycle I am giving
> > 100 events, but I am getting only 1. I want to use it as a input in
> > pythia. Is it possible to generate 100 events for h->b B?
> >
> > On Tue, Oct 17, 2023 at 2:15 AM Alexander Pukhov <
> > <email address hidden>> wrote:
> >
> >> Your question #708114 on CalcHEP changed:
> >> https:/
> >>
> >> Alexander Pukhov proposed the following answer:
> >> Sorry. I forget to answer.
> >>
> >> Try the following way. First, create events in CalcHEP format using
> >> CalcHEP menu system and then use
> >>
> >> CalcHEP/
> >>
> >> Let me now in case of problem.
> >>
> >> Best
> >>
> >> Alexander Pukhov
> >>
> >>
> >> On 10/5/23 11:50, PARTHA KUMAR PAUL wrote:
> >>> New question #708114 on CalcHEP:
> >>> https:/
> >>>
> >>> How to generate lhe file for a decay process, e.g. h-> b b in calchep?
> >>> I am able to generate lhe files for 2->2, 2->1 process, but not for
> 1->2
> >> process.
> >>>
> >> --
> >> If this answers your question, please go to the following page to let us
> >> know that it is solved:
> >>
> https:/
> >>
> >> If you still need help, you can reply to this email or go to the
> >> following page to enter your feedback:
> >> https:/
> >>
> >> You received this question notification because you asked the question.
> >>
> >
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
> https:/
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https:/
>
> You received this question notification because you asked the question.
>
--
*Partha Kumar Paul*
Revision history for this message
|
#6 |
mix_events generates events for collisions. If Higgs particle is
produced in collision and decaySLHA.txt file exists, then event_mixer
add realizes h-decay in events.
But you asks about intermediate part of this job - lhe file for events
of h-decay.
Please, confirm your request.
Best
Alexander Pukhov
On 10/23/23 09:50, PARTHA KUMAR PAUL wrote:
> Question #708114 on CalcHEP changed:
> https:/
>
> Status: Answered => Open
>
> PARTHA KUMAR PAUL is still having a problem:
> Thank you for the reply. I am able to generate events in calcHEP format
> now. Using event2lhe, I am able to convert that into lhe format (in
> terminal). But how to save it in lhe format?
> When I am generating events for h->2*x, I am getting multiple event files
> corresponding to each decay mode. If I try to combine these .txt files into
> one lhe using the event_mixer, the following error comes:
>
> Total width of 25 is not found
> There are no events
>
> On Sun, Oct 22, 2023 at 3:41 PM Alexander Pukhov <
> <email address hidden>> wrote:
>
>> Your question #708114 on CalcHEP changed:
>> https:/
>>
>> Status: Open => Answered
>>
>> Alexander Pukhov proposed the following answer:
>> Thanks a lot for your mail.
>>
>> It is indeed a bug caused by disagreement of changes. Namely here
>> n_calchep waits for natural number for requested number of events, but
>> subproc_cycle substitutes a floating point one. To fix the bug, please,
>> open file
>>
>> CalcHEP/
>>
>> The line 570 contains the corresponding instruction
>>
>> case 10: correctInt(
>>
>> It has to be replaced on
>>
>> case 10:
>> { double nEventsR=nEvents;
>> correctDouble(
>> nEvents=nEventsR;
>> break;
>> }
>>
>> Best
>>
>> Alexander Pukhov
>>
>>
>>
>> Try the following. Set
>>
>> "Freeze grid OFF"
>>
>> using CalcHEP menus, and after that call
>>
>> ../../bin/
>>
>> Unfortunately subproc_cycle finishes with
>>
>> "Freeze grid ON"
>>
>> and for next call you will need to improve it via menus.
>>
>>
>> Best
>>
>> Alexander Pukhov
>>
>>
>> On 10/20/23 20:35, PARTHA KUMAR PAUL wrote:
>>> Question #708114 on CalcHEP changed:
>>> https:/
>>>
>>> Status: Answered => Open
>>>
>>> PARTHA KUMAR PAUL is still having a problem:
>>> Thank you Prof. Pukhov, this is working. In the subproc_cycle I am giving
>>> 100 events, but I am getting only 1. I want to use it as a input in
>>> pythia. Is it possible to generate 100 events for h->b B?
>>>
>>> On Tue, Oct 17, 2023 at 2:15 AM Alexander Pukhov <
>>> <email address hidden>> wrote:
>>>
>>>> Your question #708114 on CalcHEP changed:
>>>> https:/
>>>>
>>>> Alexander Pukhov proposed the following answer:
>>>> Sorry. I forget to answer.
>>>>
>>>> Try the following way. First, create events in CalcHEP format using
>>>> CalcHEP menu system and then use
>>>>
>>>> CalcHEP/
>>>>
>>>> Let me now in case of problem.
>>>>
>>>> Best
>>>>
>>>> Alexander Pukhov
>>>>
>>>>
>>>> On 10/5/23 11:50, PARTHA KUMAR PAUL wrote:
>>>>> New question #708114 on CalcHEP:
>>>>> https:/
>>>>>
>>>>> How to generate lhe file for a decay process, e.g. h-> b b in calchep?
>>>>> I am able to generate lhe files for 2->2, 2->1 process, but not for
>> 1->2
>>>> process.
>>>> --
>>>> If this answers your question, please go to the following page to let us
>>>> know that it is solved:
>>>>
>> https:/
>>>> If you still need help, you can reply to this email or go to the
>>>> following page to enter your feedback:
>>>> https:/
>>>>
>>>> You received this question notification because you asked the question.
>>>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>> https:/
>>
>> If you still need help, you can reply to this email or go to the
>> following page to enter your feedback:
>> https:/
>>
>> You received this question notification because you asked the question.
>>
>
Revision history for this message
|
#7 |
h is not getting produced from a collision. h is there, and it is decaying
to b bbar. I want to generate 1000 events for this decay process in lhe
file format. I am able to generate in the calchep format. But not able to
convert it into lhe format.
On Sat, Oct 28, 2023 at 9:40 PM Alexander Pukhov <
<email address hidden>> wrote:
> Your question #708114 on CalcHEP changed:
> https:/
>
> Status: Open => Answered
>
> Alexander Pukhov proposed the following answer:
>
> mix_events generates events for collisions. If Higgs particle is
> produced in collision and decaySLHA.txt file exists, then event_mixer
> add realizes h-decay in events.
>
> But you asks about intermediate part of this job - lhe file for events
> of h-decay.
>
> Please, confirm your request.
>
> Best
>
> Alexander Pukhov
>
>
> On 10/23/23 09:50, PARTHA KUMAR PAUL wrote:
> > Question #708114 on CalcHEP changed:
> > https:/
> >
> > Status: Answered => Open
> >
> > PARTHA KUMAR PAUL is still having a problem:
> > Thank you for the reply. I am able to generate events in calcHEP format
> > now. Using event2lhe, I am able to convert that into lhe format (in
> > terminal). But how to save it in lhe format?
> > When I am generating events for h->2*x, I am getting multiple event files
> > corresponding to each decay mode. If I try to combine these .txt files
> into
> > one lhe using the event_mixer, the following error comes:
> >
> > Total width of 25 is not found
> > There are no events
> >
> > On Sun, Oct 22, 2023 at 3:41 PM Alexander Pukhov <
> > <email address hidden>> wrote:
> >
> >> Your question #708114 on CalcHEP changed:
> >> https:/
> >>
> >> Status: Open => Answered
> >>
> >> Alexander Pukhov proposed the following answer:
> >> Thanks a lot for your mail.
> >>
> >> It is indeed a bug caused by disagreement of changes. Namely here
> >> n_calchep waits for natural number for requested number of events, but
> >> subproc_cycle substitutes a floating point one. To fix the bug, please,
> >> open file
> >>
> >> CalcHEP/
> >>
> >> The line 570 contains the corresponding instruction
> >>
> >> case 10: correctInt(
> >>
> >> It has to be replaced on
> >>
> >> case 10:
> >> { double nEventsR=nEvents;
> >> correctDouble(
> >> nEvents=nEventsR;
> >> break;
> >> }
> >>
> >> Best
> >>
> >> Alexander Pukhov
> >>
> >>
> >>
> >> Try the following. Set
> >>
> >> "Freeze grid OFF"
> >>
> >> using CalcHEP menus, and after that call
> >>
> >> ../../bin/
> >>
> >> Unfortunately subproc_cycle finishes with
> >>
> >> "Freeze grid ON"
> >>
> >> and for next call you will need to improve it via menus.
> >>
> >>
> >> Best
> >>
> >> Alexander Pukhov
> >>
> >>
> >> On 10/20/23 20:35, PARTHA KUMAR PAUL wrote:
> >>> Question #708114 on CalcHEP changed:
> >>> https:/
> >>>
> >>> Status: Answered => Open
> >>>
> >>> PARTHA KUMAR PAUL is still having a problem:
> >>> Thank you Prof. Pukhov, this is working. In the subproc_cycle I am
> giving
> >>> 100 events, but I am getting only 1. I want to use it as a input in
> >>> pythia. Is it possible to generate 100 events for h->b B?
> >>>
> >>> On Tue, Oct 17, 2023 at 2:15 AM Alexander Pukhov <
> >>> <email address hidden>> wrote:
> >>>
> >>>> Your question #708114 on CalcHEP changed:
> >>>> https:/
> >>>>
> >>>> Alexander Pukhov proposed the following answer:
> >>>> Sorry. I forget to answer.
> >>>>
> >>>> Try the following way. First, create events in CalcHEP format using
> >>>> CalcHEP menu system and then use
> >>>>
> >>>> CalcHEP/
> >>>>
> >>>> Let me now in case of problem.
> >>>>
> >>>> Best
> >>>>
> >>>> Alexander Pukhov
> >>>>
> >>>>
> >>>> On 10/5/23 11:50, PARTHA KUMAR PAUL wrote:
> >>>>> New question #708114 on CalcHEP:
> >>>>> https:/
> >>>>>
> >>>>> How to generate lhe file for a decay process, e.g. h-> b b in
> calchep?
> >>>>> I am able to generate lhe files for 2->2, 2->1 process, but not for
> >> 1->2
> >>>> process.
> >>>> --
> >>>> If this answers your question, please go to the following page to let
> us
> >>>> know that it is solved:
> >>>>
> >>
> https:/
> >>>> If you still need help, you can reply to this email or go to the
> >>>> following page to enter your feedback:
> >>>> https:/
> >>>>
> >>>> You received this question notification because you asked the
> question.
> >>>>
> >> --
> >> If this answers your question, please go to the following page to let us
> >> know that it is solved:
> >>
> https:/
> >>
> >> If you still need help, you can reply to this email or go to the
> >> following page to enter your feedback:
> >> https:/
> >>
> >> You received this question notification because you asked the question.
> >>
> >
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
> https:/
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https:/
>
> You received this question notification because you asked the question.
>
--
*Partha Kumar Paul*
Revision history for this message
|
#8 |
So, you create file with h->b,B events, say , events_1.txt and apply
../../bin/event2lhe events_1.txt > events_1.lhe
Best
Alexander Pukhov
On 10/29/23 09:55, PARTHA KUMAR PAUL wrote:
> Question #708114 on CalcHEP changed:
> https:/
>
> Status: Answered => Open
>
> PARTHA KUMAR PAUL is still having a problem:
> h is not getting produced from a collision. h is there, and it is decaying
> to b bbar. I want to generate 1000 events for this decay process in lhe
> file format. I am able to generate in the calchep format. But not able to
> convert it into lhe format.
>
> On Sat, Oct 28, 2023 at 9:40 PM Alexander Pukhov <
> <email address hidden>> wrote:
>
>> Your question #708114 on CalcHEP changed:
>> https:/
>>
>> Status: Open => Answered
>>
>> Alexander Pukhov proposed the following answer:
>>
>> mix_events generates events for collisions. If Higgs particle is
>> produced in collision and decaySLHA.txt file exists, then event_mixer
>> add realizes h-decay in events.
>>
>> But you asks about intermediate part of this job - lhe file for events
>> of h-decay.
>>
>> Please, confirm your request.
>>
>> Best
>>
>> Alexander Pukhov
>>
>>
>> On 10/23/23 09:50, PARTHA KUMAR PAUL wrote:
>>> Question #708114 on CalcHEP changed:
>>> https:/
>>>
>>> Status: Answered => Open
>>>
>>> PARTHA KUMAR PAUL is still having a problem:
>>> Thank you for the reply. I am able to generate events in calcHEP format
>>> now. Using event2lhe, I am able to convert that into lhe format (in
>>> terminal). But how to save it in lhe format?
>>> When I am generating events for h->2*x, I am getting multiple event files
>>> corresponding to each decay mode. If I try to combine these .txt files
>> into
>>> one lhe using the event_mixer, the following error comes:
>>>
>>> Total width of 25 is not found
>>> There are no events
>>>
>>> On Sun, Oct 22, 2023 at 3:41 PM Alexander Pukhov <
>>> <email address hidden>> wrote:
>>>
>>>> Your question #708114 on CalcHEP changed:
>>>> https:/
>>>>
>>>> Status: Open => Answered
>>>>
>>>> Alexander Pukhov proposed the following answer:
>>>> Thanks a lot for your mail.
>>>>
>>>> It is indeed a bug caused by disagreement of changes. Namely here
>>>> n_calchep waits for natural number for requested number of events, but
>>>> subproc_cycle substitutes a floating point one. To fix the bug, please,
>>>> open file
>>>>
>>>> CalcHEP/
>>>>
>>>> The line 570 contains the corresponding instruction
>>>>
>>>> case 10: correctInt(
>>>>
>>>> It has to be replaced on
>>>>
>>>> case 10:
>>>> { double nEventsR=nEvents;
>>>> correctDouble(
>>>> nEvents=nEventsR;
>>>> break;
>>>> }
>>>>
>>>> Best
>>>>
>>>> Alexander Pukhov
>>>>
>>>>
>>>>
>>>> Try the following. Set
>>>>
>>>> "Freeze grid OFF"
>>>>
>>>> using CalcHEP menus, and after that call
>>>>
>>>> ../../bin/
>>>>
>>>> Unfortunately subproc_cycle finishes with
>>>>
>>>> "Freeze grid ON"
>>>>
>>>> and for next call you will need to improve it via menus.
>>>>
>>>>
>>>> Best
>>>>
>>>> Alexander Pukhov
>>>>
>>>>
>>>> On 10/20/23 20:35, PARTHA KUMAR PAUL wrote:
>>>>> Question #708114 on CalcHEP changed:
>>>>> https:/
>>>>>
>>>>> Status: Answered => Open
>>>>>
>>>>> PARTHA KUMAR PAUL is still having a problem:
>>>>> Thank you Prof. Pukhov, this is working. In the subproc_cycle I am
>> giving
>>>>> 100 events, but I am getting only 1. I want to use it as a input in
>>>>> pythia. Is it possible to generate 100 events for h->b B?
>>>>>
>>>>> On Tue, Oct 17, 2023 at 2:15 AM Alexander Pukhov <
>>>>> <email address hidden>> wrote:
>>>>>
>>>>>> Your question #708114 on CalcHEP changed:
>>>>>> https:/
>>>>>>
>>>>>> Alexander Pukhov proposed the following answer:
>>>>>> Sorry. I forget to answer.
>>>>>>
>>>>>> Try the following way. First, create events in CalcHEP format using
>>>>>> CalcHEP menu system and then use
>>>>>>
>>>>>> CalcHEP/
>>>>>>
>>>>>> Let me now in case of problem.
>>>>>>
>>>>>> Best
>>>>>>
>>>>>> Alexander Pukhov
>>>>>>
>>>>>>
>>>>>> On 10/5/23 11:50, PARTHA KUMAR PAUL wrote:
>>>>>>> New question #708114 on CalcHEP:
>>>>>>> https:/
>>>>>>>
>>>>>>> How to generate lhe file for a decay process, e.g. h-> b b in
>> calchep?
>>>>>>> I am able to generate lhe files for 2->2, 2->1 process, but not for
>>>> 1->2
>>>>>> process.
>>>>>> --
>>>>>> If this answers your question, please go to the following page to let
>> us
>>>>>> know that it is solved:
>>>>>>
>> https:/
>>>>>> If you still need help, you can reply to this email or go to the
>>>>>> following page to enter your feedback:
>>>>>> https:/
>>>>>>
>>>>>> You received this question notification because you asked the
>> question.
>>>> --
>>>> If this answers your question, please go to the following page to let us
>>>> know that it is solved:
>>>>
>> https:/
>>>> If you still need help, you can reply to this email or go to the
>>>> following page to enter your feedback:
>>>> https:/
>>>>
>>>> You received this question notification because you asked the question.
>>>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>> https:/
>>
>> If you still need help, you can reply to this email or go to the
>> following page to enter your feedback:
>> https:/
>>
>> You received this question notification because you asked the question.
>>
>
Revision history for this message
|
#9 |
Thanks a lot. It works.