MadAnalysis 5

ctau in madanalysis5 expert mode

Asked by Sahasransu

I was trying to plot the ctau of a BSM particle in madanalysis5 using the function ctau() of MCParticleFormat in the expert mode. However, it seems to give the value zero every time. I need help in resolving the issue. The analysis snipet is given below.

for(MAuint32 i=0; i<event.mc()->particles().size(); i++) {
    if(abs(event.mc()->particles()[i].pdgid())!=9000007) continue;
    std::cout<<event.mc()->particles()[i].pdgid()<<"\t"<<event.mc()->particles()[i].ctau()<<std::endl;
}

One of the input events looks as given below.

<event>
13 0 +3.8396570e-09 6.13616900e+02 7.81860800e-03 9.87755100e-02
        3 -1 0 0 501 0 +0.0000000000e+00 +0.0000000000e+00 +4.0189104069e+02 4.0189105338e+02 1.0100000000e-01 0.0000e+00 -1.0000e+00
       -3 -1 0 0 0 501 -0.0000000000e+00 -0.0000000000e+00 -2.3422122208e+02 2.3422124385e+02 1.0100000000e-01 0.0000e+00 1.0000e+00
       23 2 1 2 0 0 +2.8421709430e-14 -2.8421709430e-14 +1.6766981861e+02 6.3611229723e+02 6.1361688912e+02 4.7141e-14 0.0000e+00
  9000007 2 3 3 0 0 -1.7107677329e+02 +1.2626355258e+02 +1.0585629305e+02 3.2386066054e+02 2.2016000000e+02 6.2026e+00 0.0000e+00
       24 2 4 4 0 0 -3.1426600489e+01 +1.0915897758e+01 -9.2309061579e-01 3.3672636996e+01 5.1192132036e+00 8.1136e-14 0.0000e+00
 -9000007 2 3 3 0 0 +1.7107677329e+02 -1.2626355258e+02 +6.1813525561e+01 3.1225163669e+02 2.2016000000e+02 5.8829e+00 0.0000e+00
      -24 2 6 6 0 0 +2.3085737333e+01 -2.7910681004e+01 +2.3444790133e+00 3.8187407777e+01 1.1866935096e+01 2.4407e-13 0.0000e+00
  9000005 1 4 4 0 0 -1.3965017280e+02 +1.1534765482e+02 +1.0677938366e+02 2.9018802355e+02 2.0000000000e+02 0.0000e+00 1.0000e+00
      -11 1 5 5 0 0 -8.8130700251e+00 +5.4630162720e+00 +3.6465129679e-01 1.0375341965e+01 5.1100000000e-04 0.0000e+00 1.0000e+00
       12 1 5 5 0 0 -2.2613530464e+01 +5.4528814858e+00 -1.2877419126e+00 2.3297295032e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
  9000005 1 6 6 0 0 +1.4799103596e+02 -9.8352871578e+01 +5.9469046548e+01 2.7406422891e+02 2.0000000000e+02 0.0000e+00 -1.0000e+00
       11 1 7 7 0 0 +2.0557214877e+01 -2.0920030374e+01 -2.3031439023e+00 2.9420251978e+01 5.1100000000e-04 0.0000e+00 -1.0000e+00
      -12 1 7 7 0 0 +2.5285224559e+00 -6.9906506299e+00 +4.6476229156e+00 8.7671557991e+00 0.0000000000e+00 0.0000e+00 1.0000e+00
<mgrwt>
<rscale> 0 0.61361689E+03</rscale>
<asrwt>0</asrwt>
<pdfrwt beam="1"> 1 3 0.61829393E-01 0.61361689E+03</pdfrwt>
<pdfrwt beam="2"> 1 -3 0.36034037E-01 0.61361689E+03</pdfrwt>
<totfact> 0.51104293E+01</totfact>
</mgrwt>
<rwgt>
...
</event>

Question information

Language:
English Edit question
Status:
Solved
For:
MadAnalysis 5 Edit question
Assignee:
No assignee Edit question
Solved by:
Jack Y. Araz
Solved:
Last query:
Last reply:
Revision history for this message
Jack Y. Araz (jackaraz) said : 		#1

Hi

Could you send us 10 events so that we try to reproduce your results. Also please make sure that you are using the latest version of Ma5.

Cheers

—
Jack Y. Araz

> On Oct 27, 2020, at 11:50 AM, Sahasransu <email address hidden> wrote:
>
> New question #693668 on MadAnalysis 5:
> https://answers.launchpad.net/madanalysis5/+question/693668
>
> I was trying to plot the ctau of a BSM particle in madanalysis5 using the function ctau() of MCParticleFormat in the expert mode. However, it seems to give the value zero every time. I need help in resolving the issue. The analysis snipet is given below.
>
> for(MAuint32 i=0; i<event.mc()->particles().size(); i++) {
> if(abs(event.mc()->particles()[i].pdgid())!=9000007) continue;
> std::cout<<event.mc()->particles()[i].pdgid()<<"\t"<<event.mc()->particles()[i].ctau()<<std::endl;
> }
>
> One of the input events looks as given below.
>
> <event>
> 13 0 +3.8396570e-09 6.13616900e+02 7.81860800e-03 9.87755100e-02
> 3 -1 0 0 501 0 +0.0000000000e+00 +0.0000000000e+00 +4.0189104069e+02 4.0189105338e+02 1.0100000000e-01 0.0000e+00 -1.0000e+00
> -3 -1 0 0 0 501 -0.0000000000e+00 -0.0000000000e+00 -2.3422122208e+02 2.3422124385e+02 1.0100000000e-01 0.0000e+00 1.0000e+00
> 23 2 1 2 0 0 +2.8421709430e-14 -2.8421709430e-14 +1.6766981861e+02 6.3611229723e+02 6.1361688912e+02 4.7141e-14 0.0000e+00
> 9000007 2 3 3 0 0 -1.7107677329e+02 +1.2626355258e+02 +1.0585629305e+02 3.2386066054e+02 2.2016000000e+02 6.2026e+00 0.0000e+00
> 24 2 4 4 0 0 -3.1426600489e+01 +1.0915897758e+01 -9.2309061579e-01 3.3672636996e+01 5.1192132036e+00 8.1136e-14 0.0000e+00
> -9000007 2 3 3 0 0 +1.7107677329e+02 -1.2626355258e+02 +6.1813525561e+01 3.1225163669e+02 2.2016000000e+02 5.8829e+00 0.0000e+00
> -24 2 6 6 0 0 +2.3085737333e+01 -2.7910681004e+01 +2.3444790133e+00 3.8187407777e+01 1.1866935096e+01 2.4407e-13 0.0000e+00
> 9000005 1 4 4 0 0 -1.3965017280e+02 +1.1534765482e+02 +1.0677938366e+02 2.9018802355e+02 2.0000000000e+02 0.0000e+00 1.0000e+00
> -11 1 5 5 0 0 -8.8130700251e+00 +5.4630162720e+00 +3.6465129679e-01 1.0375341965e+01 5.1100000000e-04 0.0000e+00 1.0000e+00
> 12 1 5 5 0 0 -2.2613530464e+01 +5.4528814858e+00 -1.2877419126e+00 2.3297295032e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
> 9000005 1 6 6 0 0 +1.4799103596e+02 -9.8352871578e+01 +5.9469046548e+01 2.7406422891e+02 2.0000000000e+02 0.0000e+00 -1.0000e+00
> 11 1 7 7 0 0 +2.0557214877e+01 -2.0920030374e+01 -2.3031439023e+00 2.9420251978e+01 5.1100000000e-04 0.0000e+00 -1.0000e+00
> -12 1 7 7 0 0 +2.5285224559e+00 -6.9906506299e+00 +4.6476229156e+00 8.7671557991e+00 0.0000000000e+00 0.0000e+00 1.0000e+00
> <mgrwt>
> <rscale> 0 0.61361689E+03</rscale>
> <asrwt>0</asrwt>
> <pdfrwt beam="1"> 1 3 0.61829393E-01 0.61361689E+03</pdfrwt>
> <pdfrwt beam="2"> 1 -3 0.36034037E-01 0.61361689E+03</pdfrwt>
> <totfact> 0.51104293E+01</totfact>
> </mgrwt>
> <rwgt>
> ...
> </event>
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Sahasransu (arsahasransu) said : 		#2

Dear Jack,

Please find a sample with 100 events attached below. The version of Ma5
that I am using is MA5 release : 1.8.45.

Yours sincerely,
A.R.Sahasransu

On Tue, Oct 27, 2020 at 1:05 PM Jack Y. Araz <
<email address hidden>> wrote:

> Your question #693668 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/693668
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> Could you send us 10 events so that we try to reproduce your results.
> Also please make sure that you are using the latest version of Ma5.
>
> Cheers
>
> —
> Jack Y. Araz
>
> > On Oct 27, 2020, at 11:50 AM, Sahasransu <
> <email address hidden>> wrote:
> >
> > New question #693668 on MadAnalysis 5:
> > https://answers.launchpad.net/madanalysis5/+question/693668
> >
> > I was trying to plot the ctau of a BSM particle in madanalysis5 using
> the function ctau() of MCParticleFormat in the expert mode. However, it
> seems to give the value zero every time. I need help in resolving the
> issue. The analysis snipet is given below.
> >
> > for(MAuint32 i=0; i<event.mc()->particles().size(); i++) {
> > if(abs(event.mc()->particles()[i].pdgid())!=9000007) continue;
> > std::cout<<event.mc()->particles()[i].pdgid()<<"\t"<<event.mc
> ()->particles()[i].ctau()<<std::endl;
> > }
> >
> > One of the input events looks as given below.
> >
> > <event>
> > 13 0 +3.8396570e-09 6.13616900e+02 7.81860800e-03 9.87755100e-02
> > 3 -1 0 0 501 0 +0.0000000000e+00 +0.0000000000e+00
> +4.0189104069e+02 4.0189105338e+02 1.0100000000e-01 0.0000e+00 -1.0000e+00
> > -3 -1 0 0 0 501 -0.0000000000e+00 -0.0000000000e+00
> -2.3422122208e+02 2.3422124385e+02 1.0100000000e-01 0.0000e+00 1.0000e+00
> > 23 2 1 2 0 0 +2.8421709430e-14 -2.8421709430e-14
> +1.6766981861e+02 6.3611229723e+02 6.1361688912e+02 4.7141e-14 0.0000e+00
> > 9000007 2 3 3 0 0 -1.7107677329e+02 +1.2626355258e+02
> +1.0585629305e+02 3.2386066054e+02 2.2016000000e+02 6.2026e+00 0.0000e+00
> > 24 2 4 4 0 0 -3.1426600489e+01 +1.0915897758e+01
> -9.2309061579e-01 3.3672636996e+01 5.1192132036e+00 8.1136e-14 0.0000e+00
> > -9000007 2 3 3 0 0 +1.7107677329e+02 -1.2626355258e+02
> +6.1813525561e+01 3.1225163669e+02 2.2016000000e+02 5.8829e+00 0.0000e+00
> > -24 2 6 6 0 0 +2.3085737333e+01 -2.7910681004e+01
> +2.3444790133e+00 3.8187407777e+01 1.1866935096e+01 2.4407e-13 0.0000e+00
> > 9000005 1 4 4 0 0 -1.3965017280e+02 +1.1534765482e+02
> +1.0677938366e+02 2.9018802355e+02 2.0000000000e+02 0.0000e+00 1.0000e+00
> > -11 1 5 5 0 0 -8.8130700251e+00 +5.4630162720e+00
> +3.6465129679e-01 1.0375341965e+01 5.1100000000e-04 0.0000e+00 1.0000e+00
> > 12 1 5 5 0 0 -2.2613530464e+01 +5.4528814858e+00
> -1.2877419126e+00 2.3297295032e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
> > 9000005 1 6 6 0 0 +1.4799103596e+02 -9.8352871578e+01
> +5.9469046548e+01 2.7406422891e+02 2.0000000000e+02 0.0000e+00 -1.0000e+00
> > 11 1 7 7 0 0 +2.0557214877e+01 -2.0920030374e+01
> -2.3031439023e+00 2.9420251978e+01 5.1100000000e-04 0.0000e+00 -1.0000e+00
> > -12 1 7 7 0 0 +2.5285224559e+00 -6.9906506299e+00
> +4.6476229156e+00 8.7671557991e+00 0.0000000000e+00 0.0000e+00 1.0000e+00
> > <mgrwt>
> > <rscale> 0 0.61361689E+03</rscale>
> > <asrwt>0</asrwt>
> > <pdfrwt beam="1"> 1 3 0.61829393E-01 0.61361689E+03</pdfrwt>
> > <pdfrwt beam="2"> 1 -3 0.36034037E-01 0.61361689E+03</pdfrwt>
> > <totfact> 0.51104293E+01</totfact>
> > </mgrwt>
> > <rwgt>
> > ...
> > </event>
> >
> > --
> > You received this question notification because you are an answer
> > contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/693668/+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/madanalysis5/+question/693668
>
> You received this question notification because you asked the question.
>

Revision history for this message
Jack Y. Araz (jackaraz) said : 		#3

Hi

There is no attachment in your message, could you please send it as dropbox/drive link so that we can access?

cheers

—
Jack Y. Araz

> On Oct 27, 2020, at 12:15 PM, Sahasransu <email address hidden> wrote:
>
> Question #693668 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/693668
>
> Status: Answered => Open
>
> Sahasransu is still having a problem:
> Dear Jack,
>
> Please find a sample with 100 events attached below. The version of Ma5
> that I am using is MA5 release : 1.8.45.
>
> Yours sincerely,
> A.R.Sahasransu
>
> On Tue, Oct 27, 2020 at 1:05 PM Jack Y. Araz <
> <email address hidden>> wrote:
>
>> Your question #693668 on MadAnalysis 5 changed:
>> https://answers.launchpad.net/madanalysis5/+question/693668
>>
>> Status: Open => Answered
>>
>> Jack Y. Araz proposed the following answer:
>> Hi
>>
>> Could you send us 10 events so that we try to reproduce your results.
>> Also please make sure that you are using the latest version of Ma5.
>>
>> Cheers
>>
>> —
>> Jack Y. Araz
>>
>>> On Oct 27, 2020, at 11:50 AM, Sahasransu <
>> <email address hidden>> wrote:
>>>
>>> New question #693668 on MadAnalysis 5:
>>> https://answers.launchpad.net/madanalysis5/+question/693668
>>>
>>> I was trying to plot the ctau of a BSM particle in madanalysis5 using
>> the function ctau() of MCParticleFormat in the expert mode. However, it
>> seems to give the value zero every time. I need help in resolving the
>> issue. The analysis snipet is given below.
>>>
>>> for(MAuint32 i=0; i<event.mc()->particles().size(); i++) {
>>> if(abs(event.mc()->particles()[i].pdgid())!=9000007) continue;
>>> std::cout<<event.mc()->particles()[i].pdgid()<<"\t"<<event.mc
>> ()->particles()[i].ctau()<<std::endl;
>>> }
>>>
>>> One of the input events looks as given below.
>>>
>>> <event>
>>> 13 0 +3.8396570e-09 6.13616900e+02 7.81860800e-03 9.87755100e-02
>>> 3 -1 0 0 501 0 +0.0000000000e+00 +0.0000000000e+00
>> +4.0189104069e+02 4.0189105338e+02 1.0100000000e-01 0.0000e+00 -1.0000e+00
>>> -3 -1 0 0 0 501 -0.0000000000e+00 -0.0000000000e+00
>> -2.3422122208e+02 2.3422124385e+02 1.0100000000e-01 0.0000e+00 1.0000e+00
>>> 23 2 1 2 0 0 +2.8421709430e-14 -2.8421709430e-14
>> +1.6766981861e+02 6.3611229723e+02 6.1361688912e+02 4.7141e-14 0.0000e+00
>>> 9000007 2 3 3 0 0 -1.7107677329e+02 +1.2626355258e+02
>> +1.0585629305e+02 3.2386066054e+02 2.2016000000e+02 6.2026e+00 0.0000e+00
>>> 24 2 4 4 0 0 -3.1426600489e+01 +1.0915897758e+01
>> -9.2309061579e-01 3.3672636996e+01 5.1192132036e+00 8.1136e-14 0.0000e+00
>>> -9000007 2 3 3 0 0 +1.7107677329e+02 -1.2626355258e+02
>> +6.1813525561e+01 3.1225163669e+02 2.2016000000e+02 5.8829e+00 0.0000e+00
>>> -24 2 6 6 0 0 +2.3085737333e+01 -2.7910681004e+01
>> +2.3444790133e+00 3.8187407777e+01 1.1866935096e+01 2.4407e-13 0.0000e+00
>>> 9000005 1 4 4 0 0 -1.3965017280e+02 +1.1534765482e+02
>> +1.0677938366e+02 2.9018802355e+02 2.0000000000e+02 0.0000e+00 1.0000e+00
>>> -11 1 5 5 0 0 -8.8130700251e+00 +5.4630162720e+00
>> +3.6465129679e-01 1.0375341965e+01 5.1100000000e-04 0.0000e+00 1.0000e+00
>>> 12 1 5 5 0 0 -2.2613530464e+01 +5.4528814858e+00
>> -1.2877419126e+00 2.3297295032e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
>>> 9000005 1 6 6 0 0 +1.4799103596e+02 -9.8352871578e+01
>> +5.9469046548e+01 2.7406422891e+02 2.0000000000e+02 0.0000e+00 -1.0000e+00
>>> 11 1 7 7 0 0 +2.0557214877e+01 -2.0920030374e+01
>> -2.3031439023e+00 2.9420251978e+01 5.1100000000e-04 0.0000e+00 -1.0000e+00
>>> -12 1 7 7 0 0 +2.5285224559e+00 -6.9906506299e+00
>> +4.6476229156e+00 8.7671557991e+00 0.0000000000e+00 0.0000e+00 1.0000e+00
>>> <mgrwt>
>>> <rscale> 0 0.61361689E+03</rscale>
>>> <asrwt>0</asrwt>
>>> <pdfrwt beam="1"> 1 3 0.61829393E-01 0.61361689E+03</pdfrwt>
>>> <pdfrwt beam="2"> 1 -3 0.36034037E-01 0.61361689E+03</pdfrwt>
>>> <totfact> 0.51104293E+01</totfact>
>>> </mgrwt>
>>> <rwgt>
>>> ...
>>> </event>
>>>
>>> --
>>> You received this question notification because you are an answer
>>> contact for MadAnalysis 5.
>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>>
>> https://answers.launchpad.net/madanalysis5/+question/693668/+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/madanalysis5/+question/693668
>>
>> You received this question notification because you asked the question.
>>
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Sahasransu (arsahasransu) said : 		#4

Oh! I am sorry. Here is the link to cernbox.

https://cernbox.cern.ch/index.php/s/J339VobD7CBWfQJ

On Tue, Oct 27, 2020 at 1:30 PM Jack Y. Araz <
<email address hidden>> wrote:

> Your question #693668 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/693668
>
> Status: Open => Answered
>
> Jack Y. Araz proposed the following answer:
> Hi
>
> There is no attachment in your message, could you please send it as
> dropbox/drive link so that we can access?
>
> cheers
>
> —
> Jack Y. Araz
>
> > On Oct 27, 2020, at 12:15 PM, Sahasransu <
> <email address hidden>> wrote:
> >
> > Question #693668 on MadAnalysis 5 changed:
> > https://answers.launchpad.net/madanalysis5/+question/693668
> >
> > Status: Answered => Open
> >
> > Sahasransu is still having a problem:
> > Dear Jack,
> >
> > Please find a sample with 100 events attached below. The version of Ma5
> > that I am using is MA5 release : 1.8.45.
> >
> > Yours sincerely,
> > A.R.Sahasransu
> >
> > On Tue, Oct 27, 2020 at 1:05 PM Jack Y. Araz <
> > <email address hidden>> wrote:
> >
> >> Your question #693668 on MadAnalysis 5 changed:
> >> https://answers.launchpad.net/madanalysis5/+question/693668
> >>
> >> Status: Open => Answered
> >>
> >> Jack Y. Araz proposed the following answer:
> >> Hi
> >>
> >> Could you send us 10 events so that we try to reproduce your results.
> >> Also please make sure that you are using the latest version of Ma5.
> >>
> >> Cheers
> >>
> >> —
> >> Jack Y. Araz
> >>
> >>> On Oct 27, 2020, at 11:50 AM, Sahasransu <
> >> <email address hidden>> wrote:
> >>>
> >>> New question #693668 on MadAnalysis 5:
> >>> https://answers.launchpad.net/madanalysis5/+question/693668
> >>>
> >>> I was trying to plot the ctau of a BSM particle in madanalysis5 using
> >> the function ctau() of MCParticleFormat in the expert mode. However, it
> >> seems to give the value zero every time. I need help in resolving the
> >> issue. The analysis snipet is given below.
> >>>
> >>> for(MAuint32 i=0; i<event.mc()->particles().size(); i++) {
> >>> if(abs(event.mc()->particles()[i].pdgid())!=9000007) continue;
> >>> std::cout<<event.mc()->particles()[i].pdgid()<<"\t"<<event.mc
> >> ()->particles()[i].ctau()<<std::endl;
> >>> }
> >>>
> >>> One of the input events looks as given below.
> >>>
> >>> <event>
> >>> 13 0 +3.8396570e-09 6.13616900e+02 7.81860800e-03 9.87755100e-02
> >>> 3 -1 0 0 501 0 +0.0000000000e+00 +0.0000000000e+00
> >> +4.0189104069e+02 4.0189105338e+02 1.0100000000e-01 0.0000e+00
> -1.0000e+00
> >>> -3 -1 0 0 0 501 -0.0000000000e+00 -0.0000000000e+00
> >> -2.3422122208e+02 2.3422124385e+02 1.0100000000e-01 0.0000e+00
> 1.0000e+00
> >>> 23 2 1 2 0 0 +2.8421709430e-14 -2.8421709430e-14
> >> +1.6766981861e+02 6.3611229723e+02 6.1361688912e+02 4.7141e-14
> 0.0000e+00
> >>> 9000007 2 3 3 0 0 -1.7107677329e+02 +1.2626355258e+02
> >> +1.0585629305e+02 3.2386066054e+02 2.2016000000e+02 6.2026e+00
> 0.0000e+00
> >>> 24 2 4 4 0 0 -3.1426600489e+01 +1.0915897758e+01
> >> -9.2309061579e-01 3.3672636996e+01 5.1192132036e+00 8.1136e-14
> 0.0000e+00
> >>> -9000007 2 3 3 0 0 +1.7107677329e+02 -1.2626355258e+02
> >> +6.1813525561e+01 3.1225163669e+02 2.2016000000e+02 5.8829e+00
> 0.0000e+00
> >>> -24 2 6 6 0 0 +2.3085737333e+01 -2.7910681004e+01
> >> +2.3444790133e+00 3.8187407777e+01 1.1866935096e+01 2.4407e-13
> 0.0000e+00
> >>> 9000005 1 4 4 0 0 -1.3965017280e+02 +1.1534765482e+02
> >> +1.0677938366e+02 2.9018802355e+02 2.0000000000e+02 0.0000e+00
> 1.0000e+00
> >>> -11 1 5 5 0 0 -8.8130700251e+00 +5.4630162720e+00
> >> +3.6465129679e-01 1.0375341965e+01 5.1100000000e-04 0.0000e+00
> 1.0000e+00
> >>> 12 1 5 5 0 0 -2.2613530464e+01 +5.4528814858e+00
> >> -1.2877419126e+00 2.3297295032e+01 0.0000000000e+00 0.0000e+00
> -1.0000e+00
> >>> 9000005 1 6 6 0 0 +1.4799103596e+02 -9.8352871578e+01
> >> +5.9469046548e+01 2.7406422891e+02 2.0000000000e+02 0.0000e+00
> -1.0000e+00
> >>> 11 1 7 7 0 0 +2.0557214877e+01 -2.0920030374e+01
> >> -2.3031439023e+00 2.9420251978e+01 5.1100000000e-04 0.0000e+00
> -1.0000e+00
> >>> -12 1 7 7 0 0 +2.5285224559e+00 -6.9906506299e+00
> >> +4.6476229156e+00 8.7671557991e+00 0.0000000000e+00 0.0000e+00
> 1.0000e+00
> >>> <mgrwt>
> >>> <rscale> 0 0.61361689E+03</rscale>
> >>> <asrwt>0</asrwt>
> >>> <pdfrwt beam="1"> 1 3 0.61829393E-01 0.61361689E+03</pdfrwt>
> >>> <pdfrwt beam="2"> 1 -3 0.36034037E-01 0.61361689E+03</pdfrwt>
> >>> <totfact> 0.51104293E+01</totfact>
> >>> </mgrwt>
> >>> <rwgt>
> >>> ...
> >>> </event>
> >>>
> >>> --
> >>> You received this question notification because you are an answer
> >>> contact for MadAnalysis 5.
> >>
> >> --
> >> If this answers your question, please go to the following page to let us
> >> know that it is solved:
> >>
> >>
> https://answers.launchpad.net/madanalysis5/+question/693668/+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/madanalysis5/+question/693668
> >>
> >> You received this question notification because you asked the question.
> >>
> >
> > --
> > You received this question notification because you are an answer
> > contact for MadAnalysis 5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/693668/+confirm?answer_id=2
>
> 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/madanalysis5/+question/693668
>
> You received this question notification because you asked the question.
>

Revision history for this message
Best Jack Y. Araz (jackaraz) said : 		#5

Hi

I have no problem getting ctau information with freshly downloaded ma5. Here is the piece of code that I used and the output;

script:

  INFO << "===========" << endmsg;
  for (const MCParticleFormat &part : event.mc()->particles())
    {
      if (part.ctau() == 0.) continue;
      INFO << part.pdgid() << " -> " << part.ctau() << endmsg;
    }

last two output:

===========
23 -> 6.0101e-14
9000007 -> 1.8977
24 -> 6.6995e-14
-9000007 -> 0.33707
-24 -> 1.0103e-13
===========
23 -> 5.6666e-14
9000007 -> 5.2558
24 -> 1.0191e-13
-9000007 -> 6.4553
-24 -> 1.0069e-13

MA5 version : 1.8.58

Cheers

—
Jack Y. Araz

> On Oct 27, 2020, at 12:41 PM, Sahasransu <email address hidden> wrote:
>
> Question #693668 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/693668
>
> Status: Answered => Open
>
> Sahasransu is still having a problem:
> Oh! I am sorry. Here is the link to cernbox.
>
> https://cernbox.cern.ch/index.php/s/J339VobD7CBWfQJ
>
> On Tue, Oct 27, 2020 at 1:30 PM Jack Y. Araz <
> <email address hidden>> wrote:
>
>> Your question #693668 on MadAnalysis 5 changed:
>> https://answers.launchpad.net/madanalysis5/+question/693668
>>
>> Status: Open => Answered
>>
>> Jack Y. Araz proposed the following answer:
>> Hi
>>
>> There is no attachment in your message, could you please send it as
>> dropbox/drive link so that we can access?
>>
>> cheers
>>
>> —
>> Jack Y. Araz
>>
>>> On Oct 27, 2020, at 12:15 PM, Sahasransu <
>> <email address hidden>> wrote:
>>>
>>> Question #693668 on MadAnalysis 5 changed:
>>> https://answers.launchpad.net/madanalysis5/+question/693668
>>>
>>> Status: Answered => Open
>>>
>>> Sahasransu is still having a problem:
>>> Dear Jack,
>>>
>>> Please find a sample with 100 events attached below. The version of Ma5
>>> that I am using is MA5 release : 1.8.45.
>>>
>>> Yours sincerely,
>>> A.R.Sahasransu
>>>
>>> On Tue, Oct 27, 2020 at 1:05 PM Jack Y. Araz <
>>> <email address hidden>> wrote:
>>>
>>>> Your question #693668 on MadAnalysis 5 changed:
>>>> https://answers.launchpad.net/madanalysis5/+question/693668
>>>>
>>>> Status: Open => Answered
>>>>
>>>> Jack Y. Araz proposed the following answer:
>>>> Hi
>>>>
>>>> Could you send us 10 events so that we try to reproduce your results.
>>>> Also please make sure that you are using the latest version of Ma5.
>>>>
>>>> Cheers
>>>>
>>>> —
>>>> Jack Y. Araz
>>>>
>>>>> On Oct 27, 2020, at 11:50 AM, Sahasransu <
>>>> <email address hidden>> wrote:
>>>>>
>>>>> New question #693668 on MadAnalysis 5:
>>>>> https://answers.launchpad.net/madanalysis5/+question/693668
>>>>>
>>>>> I was trying to plot the ctau of a BSM particle in madanalysis5 using
>>>> the function ctau() of MCParticleFormat in the expert mode. However, it
>>>> seems to give the value zero every time. I need help in resolving the
>>>> issue. The analysis snipet is given below.
>>>>>
>>>>> for(MAuint32 i=0; i<event.mc()->particles().size(); i++) {
>>>>> if(abs(event.mc()->particles()[i].pdgid())!=9000007) continue;
>>>>> std::cout<<event.mc()->particles()[i].pdgid()<<"\t"<<event.mc
>>>> ()->particles()[i].ctau()<<std::endl;
>>>>> }
>>>>>
>>>>> One of the input events looks as given below.
>>>>>
>>>>> <event>
>>>>> 13 0 +3.8396570e-09 6.13616900e+02 7.81860800e-03 9.87755100e-02
>>>>> 3 -1 0 0 501 0 +0.0000000000e+00 +0.0000000000e+00
>>>> +4.0189104069e+02 4.0189105338e+02 1.0100000000e-01 0.0000e+00
>> -1.0000e+00
>>>>> -3 -1 0 0 0 501 -0.0000000000e+00 -0.0000000000e+00
>>>> -2.3422122208e+02 2.3422124385e+02 1.0100000000e-01 0.0000e+00
>> 1.0000e+00
>>>>> 23 2 1 2 0 0 +2.8421709430e-14 -2.8421709430e-14
>>>> +1.6766981861e+02 6.3611229723e+02 6.1361688912e+02 4.7141e-14
>> 0.0000e+00
>>>>> 9000007 2 3 3 0 0 -1.7107677329e+02 +1.2626355258e+02
>>>> +1.0585629305e+02 3.2386066054e+02 2.2016000000e+02 6.2026e+00
>> 0.0000e+00
>>>>> 24 2 4 4 0 0 -3.1426600489e+01 +1.0915897758e+01
>>>> -9.2309061579e-01 3.3672636996e+01 5.1192132036e+00 8.1136e-14
>> 0.0000e+00
>>>>> -9000007 2 3 3 0 0 +1.7107677329e+02 -1.2626355258e+02
>>>> +6.1813525561e+01 3.1225163669e+02 2.2016000000e+02 5.8829e+00
>> 0.0000e+00
>>>>> -24 2 6 6 0 0 +2.3085737333e+01 -2.7910681004e+01
>>>> +2.3444790133e+00 3.8187407777e+01 1.1866935096e+01 2.4407e-13
>> 0.0000e+00
>>>>> 9000005 1 4 4 0 0 -1.3965017280e+02 +1.1534765482e+02
>>>> +1.0677938366e+02 2.9018802355e+02 2.0000000000e+02 0.0000e+00
>> 1.0000e+00
>>>>> -11 1 5 5 0 0 -8.8130700251e+00 +5.4630162720e+00
>>>> +3.6465129679e-01 1.0375341965e+01 5.1100000000e-04 0.0000e+00
>> 1.0000e+00
>>>>> 12 1 5 5 0 0 -2.2613530464e+01 +5.4528814858e+00
>>>> -1.2877419126e+00 2.3297295032e+01 0.0000000000e+00 0.0000e+00
>> -1.0000e+00
>>>>> 9000005 1 6 6 0 0 +1.4799103596e+02 -9.8352871578e+01
>>>> +5.9469046548e+01 2.7406422891e+02 2.0000000000e+02 0.0000e+00
>> -1.0000e+00
>>>>> 11 1 7 7 0 0 +2.0557214877e+01 -2.0920030374e+01
>>>> -2.3031439023e+00 2.9420251978e+01 5.1100000000e-04 0.0000e+00
>> -1.0000e+00
>>>>> -12 1 7 7 0 0 +2.5285224559e+00 -6.9906506299e+00
>>>> +4.6476229156e+00 8.7671557991e+00 0.0000000000e+00 0.0000e+00
>> 1.0000e+00
>>>>> <mgrwt>
>>>>> <rscale> 0 0.61361689E+03</rscale>
>>>>> <asrwt>0</asrwt>
>>>>> <pdfrwt beam="1"> 1 3 0.61829393E-01 0.61361689E+03</pdfrwt>
>>>>> <pdfrwt beam="2"> 1 -3 0.36034037E-01 0.61361689E+03</pdfrwt>
>>>>> <totfact> 0.51104293E+01</totfact>
>>>>> </mgrwt>
>>>>> <rwgt>
>>>>> ...
>>>>> </event>
>>>>>
>>>>> --
>>>>> You received this question notification because you are an answer
>>>>> contact for MadAnalysis 5.
>>>>
>>>> --
>>>> If this answers your question, please go to the following page to let us
>>>> know that it is solved:
>>>>
>>>>
>> https://answers.launchpad.net/madanalysis5/+question/693668/+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/madanalysis5/+question/693668
>>>>
>>>> You received this question notification because you asked the question.
>>>>
>>>
>>> --
>>> You received this question notification because you are an answer
>>> contact for MadAnalysis 5.
>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>>
>> https://answers.launchpad.net/madanalysis5/+question/693668/+confirm?answer_id=2
>>
>> 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/madanalysis5/+question/693668
>>
>> You received this question notification because you asked the question.
>>
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Sahasransu (arsahasransu) said : 		#6

Thanks Jack Y. Araz, that solved my question.