Some commands in ma5

Dear Safinaz,

> - First how to select events where the ratio of the transverse momentum of the leading photon to the invariant mass of γγ pair, pT/ mγγ > 0.35 ?
>
> I tried : select select PT(a[1])/M(a a) > 0.35, or select PT(a)/M(a a)
>> 0.35, but it gives me an error
You cannot do that. You need to use the expert mode of the program to implement such a cut.

> -Second I try to make jet clustering for b-jets by FastJet, but when I
> type :
>
> set main.clustering.algorithm = antikt
You should type
  set main.fastsim.algorithm = antikt
instead. Then it will work.

>
> - Third It's written in this paper :
> G. Aad et al. [ATLAS Collaboration], Phys. Rev. Lett. 114, 081802 (2015).
>
> that the selected Photons are required to be isolated: the energy in
> the calorimeter [3,42] within a cone of size =0.4
>
> How can this justified in MA ?
I am not sure about what you want to know, but you can either implement photon isolation with simple deltaR cuts in the normal mode of the program, or implement more sophisticated cuts in the expert mode of the program.

Regards

Benjamin

Thanks Benjamin Fuks, that solved my question.

Hi sir,

I think the question still opened .. Regarding the first point for defining new observable in MA, I read in this question

https://answers.launchpad.net/madanalysis5/+question/252806

that I can include this ratio as a new observable in "observable_list.py " in observable folder instead of going to the expert mode ..
actually i tried that , I defined in that list :

PT(a[1])/PT(a[2]) = ObservableBase( name = 'ratio',
                     args = [ArgumentType.COMBINATION],
                     combination = CombinationType.DEFAULT,
                     plot_auto = False,
                     plot_nbins = 100,
                     plot_xmin = 0.,
                     plot_xmax = 1000.,
                     plot_unitX = '1',
                     code_parton = 'pt()',
                     code_hadron = 'pt()',
                     code_reco = 'pt()',
                     cut_event = True,
                     cut_candidate = True,
                     tex = 'ratio'
                  )

So are these right arguments ..

Also for the second part of the answer " and including it in the physics services (in SampleAnalyzer). " can you please declare how is it ? there are in SampleAnalyzer/commons/service , many .h files, which one should I add ..

Many thanks,
Safinaz

Hi sir,

I think the question still opened .. Regarding the first point for defining new observable in MA, I read in this question

https://answers.launchpad.net/madanalysis5/+question/252806

that I can include this ratio as a new observable in "observable_list.py " in observable folder instead of going to the expert mode ..
actually i tried that , I defined in that list :

PT(a[1])/PT(a[2]) = ObservableBase( name = 'ratio',
                     args = [ArgumentType.COMBINATION],
                     combination = CombinationType.DEFAULT,
                     plot_auto = False,
                     plot_nbins = 100,
                     plot_xmin = 0.,
                     plot_xmax = 1000.,
                     plot_unitX = '1',
                     code_parton = 'pt()',
                     code_hadron = 'pt()',
                     code_reco = 'pt()',
                     cut_event = True,
                     cut_candidate = True,
                     tex = 'ratio'
                  )

So are these right arguments ..

Also for the second part of the answer " and including it in the physics services (in SampleAnalyzer). " can you please declare how is it ? there are in SampleAnalyzer/commons/service , many .h files, which one should I add ..

Many thanks,
Safinaz

Hi Safinaz,

In Physics.cpp.

Cheers,

Benjamin

PS: btw: 'PT(a[1])/PT(a[2])' will not work. This must be a string (a keyword).

On 06 Jul 2016, at 03:36 , Safinaz <email address hidden> wrote:

> Question #295994 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/295994
>
> Safinaz posted a new comment:
> Hi sir,
>
> I think the question still opened .. Regarding the first point for
> defining new observable in MA, I read in this question
>
> https://answers.launchpad.net/madanalysis5/+question/252806
>
> that I can include this ratio as a new observable in "observable_list.py " in observable folder instead of going to the expert mode ..
> actually i tried that , I defined in that list :
>
>
> PT(a[1])/PT(a[2]) = ObservableBase( name = 'ratio',
> args = [ArgumentType.COMBINATION],
> combination = CombinationType.DEFAULT,
> plot_auto = False,
> plot_nbins = 100,
> plot_xmin = 0.,
> plot_xmax = 1000.,
> plot_unitX = '1',
> code_parton = 'pt()',
> code_hadron = 'pt()',
> code_reco = 'pt()',
> cut_event = True,
> cut_candidate = True,
> tex = 'ratio'
> )
>
> So are these right arguments ..
>
>
> Also for the second part of the answer " and including it in the physics services (in SampleAnalyzer). " can you please declare how is it ? there are in SampleAnalyzer/commons/service , many .h files, which one should I add ..
>
> Many thanks,
> Safinaz
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

# PS: btw: 'PT(a[1])/PT(a[2])' will not work. This must be a string (a keyword).

So this way won't work for this variable ..

Can I ensure from other thing please ..

for photon isolation in p p > 2 a 2b , as I read in the experimental searches that the photons selection is based on photons depositing their energy within a cone of radius = 0.4, so to do this in MA, I used

select (a[1]) DELTAR(a[2]) = 0.4

is that right ?

Thanks,
Safinaz

Hi,

> # PS: btw: 'PT(a[1])/PT(a[2])' will not work. This must be a string (a
> keyword).
>
> So this way won't work for this variable ..
It will work but the left-hand side must be a string.

>
>
> Can I ensure from other thing please ..
>
> for photon isolation in p p > 2 a 2b , as I read in the experimental
> searches that the photons selection is based on photons depositing their
> energy within a cone of radius = 0.4, so to do this in MA, I used
>
> select (a[1]) DELTAR(a[2]) = 0.4
This means that you would keep events if you they feature two photons separated by R=0.4. I am not sure this is really what you want. Please double check.

Regards,

Benjamin