MadAnalysis 5

bracket is wrong for the complex conditions

Asked by lily

Hi,
Using the user-friendly framework of Madanalysis5 (version: Date : 2014/08/01, 1.1.12.1), I analyze the p p > A H, A > ta+ ta-, H > ta+ ta- with m_A = 50 GeV and m_H = 300 GeV.
I want to make a cut on the invariant mass of two pair taus.

I input a command,
select (40 < M(ta+[1] ta-[1]) < 60 and 290 < M(ta+[2] ta-[2]) < 310) or (290 < M(ta+[1] ta-[1]) < 310 and 40 < M(ta+[2] ta-[2]) < 60) or (40 < M(ta+[1] ta-[2]) < 60 and 290 < M(ta+[2] ta-[1]) < 310) or (290 < M(ta+[1] ta-[2]) < 310 and 40 < M(ta+[2] ta-[1]) < 60)

Madanalysis5 report: "** ERROR: ')' must be a float value."

However, If I delete the final a pair "( )", the Madanalysis5 can work. The command is,
select (40 < M(ta+[1] ta-[1]) < 60 and 290 < M(ta+[2] ta-[2]) < 310) or (290 < M(ta+[1] ta-[1]) < 310 and 40 < M(ta+[2] ta-[2]) < 60) or (40 < M(ta+[1] ta-[2]) < 60 and 290 < M(ta+[2] ta-[1]) < 310) or 290 < M(ta+[1] ta-[2]) < 310 and 40 < M(ta+[2] ta-[1]) < 60

It looks strange. What is wrong?
Many thanks.

Question information

Language:
English Edit question
Status:
Answered
For:
MadAnalysis 5 Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Benjamin Fuks (fuks) said : 		#1

Hi,

This was in fact a bug. I have fixed the problem. Could you please update your version (via bazaar) of madanalysis5 and let me know whether everything is now fine?

Cheers,

Benjamin

Revision history for this message
lily (lilymag) said : 		#2

Hi Benjamin,

For "select 60 < PT(mu+[1]) < 100 and 60 < PT(mu-[1]) < 100", Ma5 still selects the events follow "60 < PT(mu+[1]) < 100". Where "and 60 < PT(mu-[1]) < 100" is not valid.

Also the problem appears for the "select 60 < PT(mu+[1]) < 100 or 60 < PT(mu-[1]) < 100".

This should be a bug.

Many thanks.

Revision history for this message
lily (lilymag) said : 		#3

Hi,

In addition, for the command " select (PT(mu[2]) > 10 and PT(j[2]) > 10) or (PT(mu[2]) > 10 and PT(j[1]) > 10)", Ma5 still selects the events follow "(PT(mu[2]) > 10 and PT(j[2]) > 10) ". Where "and PT(j[2]) > 10) or (PT(mu[2]) > 10 and PT(j[1]) > 10)" is not valid.

Many thanks

Revision history for this message
Benjamin Fuks (fuks) said : 		#4

Hi Lily,

I cannot reproduce the issue: I have double-checked and MadAnalysis5 works perfectly in this case. Have you tried with the latest version of the program?

Cheers,

Benjamin

On 03 Oct 2015, at 14:42 , lily <email address hidden> wrote:

> Question #271824 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/271824
>
> Status: Answered => Open
>
> lily is still having a problem:
> Hi Benjamin,
>
> For "select 60 < PT(mu+[1]) < 100 and 60 < PT(mu-[1]) < 100", Ma5 still
> selects the events follow "60 < PT(mu+[1]) < 100". Where "and 60 <
> PT(mu-[1]) < 100" is not valid.
>
> Also the problem appears for the "select 60 < PT(mu+[1]) < 100 or 60 <
> PT(mu-[1]) < 100".
>
> This should be a bug.
>
> Many thanks.
>
> --
> You received this question notification because you are subscribed to
> the question.

Revision history for this message
lily (lilymag) said : 		#5

Hi Benjamin,
Thank you for your reply. I checked it again, but the problem appeared.

My operations are:

1. bzr branch lp:madanalysis5

Although I obtain a report "You have not informed bzr of your Launchpad ID, and you must do this to
write to Launchpad or access private data. See "bzr help launchpad-login".
Branched 96 revisions.", but madanalysis5 (up to 96-1.3-(2015/10/01) )is downloaded. I install the fastjet and delphes in the Ma5.

2. I generate 100000 p p > z z, z > mu+ mu-, z > mu+ mu- using madgraph5, and the .root file is automatically produced by the madgraph5.

3. I use the user-friendly framework of Madanalysis5

bin/ma5 -R
import .root

select 60 < PT(mu+[1]) < 100
Events kept:30583

select 60 < PT(mu-[1]) < 100
Events kept:9116

4. I again run the Ma5 and import the same .root file.
select 60 < PT(mu+[1]) < 100 and 60 < PT(mu-[1]) < 100
Events kept:30583

Obviously, here "and 60 < PT(mu-[1]) < 100" is not valid.

Also, I download "MadAnalysis5_v1.2_patch2.tar.gz" directly from the "https://launchpad.net/madanalysis5/+download".
However, the installation failed. The report is:

" Component 6/9 - library: interface to Delphes
     - Cleaning the project before building the library ...
     - Compiling the source files ...
     - Linking the library ...
 ** ERROR: impossible to link the project. For more details, see the log file:
 ** ERROR: /home/mylocal/madanalysis5/tools/SampleAnalyzer/Interfaces/linking_delphes.log
 ** ERROR: The library building aborted."

linking_delphes.log gives:
"g++ -shared -o ../Lib/libdelphes_for_ma5.so delphes/DelphesDataFormat2.o delphes/DelphesDataFormat.o delphes/DelphesMemoryInterface.o delphes/DelphesTreeReader.o delphes/DetectorDelphes.o -L/home/mylocal/madanalysis5/tools/SampleAnalyzer/Lib -L/home/mylocal/madanalysis5/tools/delphes/ -lDelphes -lcommons_for_ma5 -L/home/mylocal/root/lib -lCore -lCint -lRIO -lNet -lHist -lGraf -lGraf3d -lGpad -lTree -lRint -lPostscript -lMatrix -lPhysics -lMathCore -lThread -pthread -lm -ldl -rdynamic -lEG
/home/mylocal/madanalysis5/tools/delphes//libDelphes.so: file not recognized: File format not recognized
collect2: ld returned 1 exit status
make: *** [link] Error 1

Could you help me fix the problem of installation? I can use it to check the problem again. My system is ubuntu 12.04 x86_64.

Many thanks.

Revision history for this message
Benjamin Fuks (fuks) said : 		#6

Hi Lily,

Thanks for your feedback. I created a patch3 version where the conditional bug is fixed. In principle, you can use the latest version that was working for you and replace the madanalysis/job/job_*_cut.py files by those of the patch3 version.

For the other problem, I cannot reproduce the bug, but you are not the only one suffering from it. Could you please check whether the file tools/delphes//libDelphes.so is present or not? Which libraries are present?

Cheers,

Benjamin

Can you help with this problem?

Provide an answer of your own, or ask lily for more information if necessary.

To post a message you must log in.