CalcHEP

Results are fluctuating:

Asked by Najimuddin Khan

Respected Sir,
  I am trying to generate a process p p -> u d DM DM in Calchep-3.6.25 for our model, p = u d. I have checked that the results very much fluctuating with command ../bin subpro_cycle 1000 1000000. Also the no of events were very less(30,100,237..etc.) than expected.

I would be very grateful to you if kindly solve the issue.

Question information

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

One problem could be absence of cuts for transverse momenta of jets ( u,d).
Other problem could be a small number of points used in the second cycle of calculation by subpro_cycle.
It could be improved by

./set_vegas
Error! 4 parameters expected
1 - number Vegas iterations at firts loop
2 - number of integral calls for earch iteration of the first loop
3 - number Vegas iterations at second loop
4 - number of integral calls for earch iteration of the second loop
5 - number of sub-cubes to prepare event generator

By default number of calls at second loop is too small for 2->4 process.

May be process contains some poles. It also can lead to problems.
For which model you perform calculations?

Revision history for this message
Najimuddin Khan (khanphysics-123) said : 		#2

Respected Sir,
I am working with IDM.
On Apr 26, 2016 4:47 PM, "Alexander Pukhov" <
<email address hidden>> wrote:

> Your question #292405 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/292405
>
> Status: Open => Needs information
>
> Alexander Pukhov requested more information:
> One problem could be absence of cuts for transverse momenta of jets
> ( u,d).
> Other problem could be a small number of points used in the second cycle
> of calculation by subpro_cycle.
> It could be improved by
>
> ./set_vegas
> Error! 4 parameters expected
> 1 - number Vegas iterations at firts loop
> 2 - number of integral calls for earch iteration of the first loop
> 3 - number Vegas iterations at second loop
> 4 - number of integral calls for earch iteration of the second loop
> 5 - number of sub-cubes to prepare event generator
>
> By default number of calls at second loop is too small for 2->4
> process.
>
> May be process contains some poles. It also can lead to problems.
> For which model you perform calculations?
>
> --
> To answer this request for more information, you can either reply to
> this email or enter your reply at the following page:
> https://answers.launchpad.net/calchep/+question/292405
>
> You received this question notification because you asked the question.
>

Revision history for this message
Najimuddin Khan (khanphysics-123) said : 		#3

Respected Sir,
 In this process...I have set the cuts:
 T(u) and T(d) minimum bound 10 GeV and Rapidity(Y) -5 to +5 for the same. It is also not improving after run "./set_vegas"

Thanks

Revision history for this message
Alexander Pukhov (pukhov) said : 		#4

Are you going to calculate DM production via W/Z fusion?

Alexander Pukhov

On 04/27/2016 12:27 PM, Najimuddin Khan wrote:
> Question #292405 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/292405
>
> Najimuddin Khan gave more information on the question:
> Respected Sir,
> In this process...I have set the cuts:
> T(u) and T(d) minimum bound 10 GeV and Rapidity(Y) -5 to +5 for the same. It is also not improving after run "./set_vegas"
>
>
> Thanks
>

Revision history for this message
Najimuddin Khan (khanphysics-123) said : 		#5

Yes Sir,
 I want to calculate dark matter pair with vector bosons fusion. Is there
any loss of gauge invariance? Or is there any other problem? In my
process(pp -> DM DM u d), I have considered that the p consist only with u
and d. I got there are 2 sub processes with 30+30 diagrams. I have also
checked in the graphics mode taking different nSess or nCalls and got
different results in different sessions with error(5-20%) for each sub
processes.

Thanks

On Wed, Apr 27, 2016 at 5:57 PM, Alexander Pukhov <
<email address hidden>> wrote:

> Your question #292405 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/292405
>
> Status: Open => Answered
>
> Alexander Pukhov proposed the following answer:
> Are you going to calculate DM production via W/Z fusion?
>
> Alexander Pukhov
>
>
> On 04/27/2016 12:27 PM, Najimuddin Khan wrote:
> > Question #292405 on CalcHEP changed:
> > https://answers.launchpad.net/calchep/+question/292405
> >
> > Najimuddin Khan gave more information on the question:
> > Respected Sir,
> > In this process...I have set the cuts:
> > T(u) and T(d) minimum bound 10 GeV and Rapidity(Y) -5 to +5 for the
> same. It is also not improving after run "./set_vegas"
> >
> >
> > Thanks
> >
>
> --
> 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/292405/+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/292405
>
> You received this question notification because you asked the question.
>

--
From
Najimuddin Khan
PhD Student
IIT Indore,
Khandwa Road,
Madhya Pradesh,
India-452017

Revision history for this message
Alexander Pukhov (pukhov) said : 		#6

OK, I'll check now.

Best
     Alexander Pukhov

On 04/27/2016 03:51 PM, Najimuddin Khan wrote:
> Question #292405 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/292405
>
> Status: Answered => Open
>
> Najimuddin Khan is still having a problem:
> Yes Sir,
> I want to calculate dark matter pair with vector bosons fusion. Is there
> any loss of gauge invariance? Or is there any other problem? In my
> process(pp -> DM DM u d), I have considered that the p consist only with u
> and d. I got there are 2 sub processes with 30+30 diagrams. I have also
> checked in the graphics mode taking different nSess or nCalls and got
> different results in different sessions with error(5-20%) for each sub
> processes.
>
>
> Thanks
>
> On Wed, Apr 27, 2016 at 5:57 PM, Alexander Pukhov <
> <email address hidden>> wrote:
>
>> Your question #292405 on CalcHEP changed:
>> https://answers.launchpad.net/calchep/+question/292405
>>
>> Status: Open => Answered
>>
>> Alexander Pukhov proposed the following answer:
>> Are you going to calculate DM production via W/Z fusion?
>>
>> Alexander Pukhov
>>
>>
>> On 04/27/2016 12:27 PM, Najimuddin Khan wrote:
>>> Question #292405 on CalcHEP changed:
>>> https://answers.launchpad.net/calchep/+question/292405
>>>
>>> Najimuddin Khan gave more information on the question:
>>> Respected Sir,
>>> In this process...I have set the cuts:
>>> T(u) and T(d) minimum bound 10 GeV and Rapidity(Y) -5 to +5 for the
>> same. It is also not improving after run "./set_vegas"
>>>
>>> Thanks
>>>
>> --
>> 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/292405/+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/292405
>>
>> You received this question notification because you asked the question.
>>
>

Revision history for this message
Alexander Pukhov (pukhov) said : 		#7

Yes, I see. there is 30% error. I'll try to find a way to solve it.

Best
    Alexander Pukhov

On 04/27/2016 03:51 PM, Najimuddin Khan wrote:
> Question #292405 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/292405
>
> Status: Answered => Open
>
> Najimuddin Khan is still having a problem:
> Yes Sir,
> I want to calculate dark matter pair with vector bosons fusion. Is there
> any loss of gauge invariance? Or is there any other problem? In my
> process(pp -> DM DM u d), I have considered that the p consist only with u
> and d. I got there are 2 sub processes with 30+30 diagrams. I have also
> checked in the graphics mode taking different nSess or nCalls and got
> different results in different sessions with error(5-20%) for each sub
> processes.
>
>
> Thanks
>
> On Wed, Apr 27, 2016 at 5:57 PM, Alexander Pukhov <
> <email address hidden>> wrote:
>
>> Your question #292405 on CalcHEP changed:
>> https://answers.launchpad.net/calchep/+question/292405
>>
>> Status: Open => Answered
>>
>> Alexander Pukhov proposed the following answer:
>> Are you going to calculate DM production via W/Z fusion?
>>
>> Alexander Pukhov
>>
>>
>> On 04/27/2016 12:27 PM, Najimuddin Khan wrote:
>>> Question #292405 on CalcHEP changed:
>>> https://answers.launchpad.net/calchep/+question/292405
>>>
>>> Najimuddin Khan gave more information on the question:
>>> Respected Sir,
>>> In this process...I have set the cuts:
>>> T(u) and T(d) minimum bound 10 GeV and Rapidity(Y) -5 to +5 for the
>> same. It is also not improving after run "./set_vegas"
>>>
>>> Thanks
>>>
>> --
>> 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/292405/+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/292405
>>
>> You received this question notification because you asked the question.
>>
>

Revision history for this message
Najimuddin Khan (khanphysics-123) said : 		#8

Thank you sir
On Apr 27, 2016 9:57 PM, "Alexander Pukhov" <
<email address hidden>> wrote:

> Your question #292405 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/292405
>
> Alexander Pukhov proposed the following answer:
> Yes, I see. there is 30% error. I'll try to find a way to solve it.
>
> Best
> Alexander Pukhov
>
>
> On 04/27/2016 03:51 PM, Najimuddin Khan wrote:
> > Question #292405 on CalcHEP changed:
> > https://answers.launchpad.net/calchep/+question/292405
> >
> > Status: Answered => Open
> >
> > Najimuddin Khan is still having a problem:
> > Yes Sir,
> > I want to calculate dark matter pair with vector bosons fusion. Is
> there
> > any loss of gauge invariance? Or is there any other problem? In my
> > process(pp -> DM DM u d), I have considered that the p consist only
> with u
> > and d. I got there are 2 sub processes with 30+30 diagrams. I have also
> > checked in the graphics mode taking different nSess or nCalls and got
> > different results in different sessions with error(5-20%) for each sub
> > processes.
> >
> >
> > Thanks
> >
> > On Wed, Apr 27, 2016 at 5:57 PM, Alexander Pukhov <
> > <email address hidden>> wrote:
> >
> >> Your question #292405 on CalcHEP changed:
> >> https://answers.launchpad.net/calchep/+question/292405
> >>
> >> Status: Open => Answered
> >>
> >> Alexander Pukhov proposed the following answer:
> >> Are you going to calculate DM production via W/Z fusion?
> >>
> >> Alexander Pukhov
> >>
> >>
> >> On 04/27/2016 12:27 PM, Najimuddin Khan wrote:
> >>> Question #292405 on CalcHEP changed:
> >>> https://answers.launchpad.net/calchep/+question/292405
> >>>
> >>> Najimuddin Khan gave more information on the question:
> >>> Respected Sir,
> >>> In this process...I have set the cuts:
> >>> T(u) and T(d) minimum bound 10 GeV and Rapidity(Y) -5 to +5 for
> the
> >> same. It is also not improving after run "./set_vegas"
> >>>
> >>> Thanks
> >>>
> >> --
> >> 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/292405/+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/292405
> >>
> >> 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://answers.launchpad.net/calchep/+question/292405/+confirm?answer_id=6
>
> 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/292405
>
> You received this question notification because you asked the question.
>

Revision history for this message
Best Alexander Pukhov (pukhov) said : 		#9

Here bad convergence of integration is caused by poles W->u,D;
Z->q,q, H3-> q,q,~X, and HC->q,q,~X
To smooth them I apply regularization

  #Kinematical_scheme
12 -> 34 , 56
34 -> 3 , 4
56 -> 5 , 6

#Regularization
*** Table ***
  Regularization
  Momentum |> Mass <|> Width <| Power|
34 |MZ |wZ |2
34 |MW |wW |2
345 |MH3 |wH3 |2
346 |MH3 |wH3 |2
345 |MHC |wH3 |2
346 |MHC |wHC |2
=========================================
For integration I used 2*( 5*1000008) points.
In results I have got precision of integration about 1%.
Protocol files for all 20 subprocesses are attached. session.dat with
all settings is attached as well.

In principle such regularization should be done automatically. I'll try
to realize this option.

Other comments. I have deleted small huU and hdD vertices which needs
only for direct detection in micrOMEGAs.
Also I removed loop induced AAh and GGh. I have tested that GGh is
small. ( about 1E-6pb without 2 jets).

Best
     Alexander Pukhov

On 04/27/2016 07:37 PM, Najimuddin Khan wrote:
> Question #292405 on CalcHEP changed:
> https://answers.launchpad.net/calchep/+question/292405
>
> Status: Answered => Open
>
> Najimuddin Khan is still having a problem:
> Thank you sir
> On Apr 27, 2016 9:57 PM, "Alexander Pukhov" <
> <email address hidden>> wrote:
>
>> Your question #292405 on CalcHEP changed:
>> https://answers.launchpad.net/calchep/+question/292405
>>
>> Alexander Pukhov proposed the following answer:
>> Yes, I see. there is 30% error. I'll try to find a way to solve it.
>>
>> Best
>> Alexander Pukhov
>>
>>
>> On 04/27/2016 03:51 PM, Najimuddin Khan wrote:
>>> Question #292405 on CalcHEP changed:
>>> https://answers.launchpad.net/calchep/+question/292405
>>>
>>> Status: Answered => Open
>>>
>>> Najimuddin Khan is still having a problem:
>>> Yes Sir,
>>> I want to calculate dark matter pair with vector bosons fusion. Is
>> there
>>> any loss of gauge invariance? Or is there any other problem? In my
>>> process(pp -> DM DM u d), I have considered that the p consist only
>> with u
>>> and d. I got there are 2 sub processes with 30+30 diagrams. I have also
>>> checked in the graphics mode taking different nSess or nCalls and got
>>> different results in different sessions with error(5-20%) for each sub
>>> processes.
>>>
>>>
>>> Thanks
>>>
>>> On Wed, Apr 27, 2016 at 5:57 PM, Alexander Pukhov <
>>> <email address hidden>> wrote:
>>>
>>>> Your question #292405 on CalcHEP changed:
>>>> https://answers.launchpad.net/calchep/+question/292405
>>>>
>>>> Status: Open => Answered
>>>>
>>>> Alexander Pukhov proposed the following answer:
>>>> Are you going to calculate DM production via W/Z fusion?
>>>>
>>>> Alexander Pukhov
>>>>
>>>>
>>>> On 04/27/2016 12:27 PM, Najimuddin Khan wrote:
>>>>> Question #292405 on CalcHEP changed:
>>>>> https://answers.launchpad.net/calchep/+question/292405
>>>>>
>>>>> Najimuddin Khan gave more information on the question:
>>>>> Respected Sir,
>>>>> In this process...I have set the cuts:
>>>>> T(u) and T(d) minimum bound 10 GeV and Rapidity(Y) -5 to +5 for
>> the
>>>> same. It is also not improving after run "./set_vegas"
>>>>> Thanks
>>>>>
>>>> --
>>>> 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/292405/+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/292405
>>>>
>>>> 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://answers.launchpad.net/calchep/+question/292405/+confirm?answer_id=6
>>
>> 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/292405
>>
>> You received this question notification because you asked the question.
>>

Revision history for this message
Najimuddin Khan (khanphysics-123) said : 		#10

Thanks Alexander Pukhov, that solved my question.