Error message while defining the multiparticle

Hi Shiba,

I assume that you are trying to make use of partons in the reco mode of the program. This cannot work as we only have (possibly b-tagged) jets at that level.

Cheers,

Benjamin

On 17 Jan 2016, at 11:32 , shibasipu <email address hidden> wrote:

> New question #281085 on MadAnalysis 5:
> https://answers.launchpad.net/madanalysis5/+question/281085
>
> Hello Benjamin,
> Why the multiparticle definition is not working for g ?
>
> ma5>define j = g u c d s u~ c~ d~ s~
> ** ERROR: The (multi)particle 'g' is not defined.
>
>
>
>
>
>
>
> Regards,
> Shiba
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

ok. Thanks. Suppose I want various distributions for jets. Which mode should I use for .hep files ?

Regards,
Shiba

Hi,
In general, stdhep files contain hadrons as final state particles. You should then: - use the reco mode and its fasjet interface to create a reconstructed lhe file; - analyze this lhe file with the default mode.
Cheers,
Benjamin

-------- Original message --------
From: shibasipu <email address hidden>
Date: 18/01/2016 06:17 (GMT+01:00)
To: <email address hidden>
Subject: Re: [Question #281085]: Error message while defining the multiparticle

Question #281085 on MadAnalysis 5 changed:
https://answers.launchpad.net/madanalysis5/+question/281085

    Status: Answered => Open

shibasipu is still having a problem:
ok. Thanks. Suppose I want various distributions for jets. Which mode
should I use for  .hep files  ?

Regards,
Shiba

--
You received this question notification because you are an answer
contact for MadAnalysis 5.

Hi Benjamin,
                                 I am providing the following way to do my analysis in MA5.

1. First I am generating .hep files from Madgraph interfacing with Pythia.

2. Then I am using the following script in reco mode in MA5.
set main.fastsim.package = delphes
set main.fastsim.detector = atlas
import sample/IDM/IDM200.hep.gz as sig
import sample/SM/Wplusjets.hep as Wplusjets
import sample/SM/Zplusjets.hep.gz as Zplusjets
set sig.type = signal
set Wplusjets.type = background
set Zplusjets.type = background
set sig.backcolor = red
set Wplusjets.backcolor = yellow
set Zplusjets.backcolor = green
set sig.xsection = 0.008283
set Wplusjets.xsection = 240.4
set Zplusjets.xsection = 175.5
set main.stacking_method = superimpose
set main.stacking_method = normalize2one
set main.lumi = 3000
set main.normalize = lumi
define j = g u c d s u~ c~ d~ s~
select N(j)=4
select MET > 160
select PT(j[1])>130
select PT(j[2])>60
select PT(j[3])>60
select PT(j[4])>60
plot MET
plot PT(j[1])
plot PT(j[2])
plot PT(j[3])
plot PT(j[4])
plot MT(j[1])
plot MT(j[2])
plot MT(j[3])
plot MT(j[4])
set main.SBratio = "S/sqrt(S+B)"
submit

Is this not the right way ? still I am getting the error message.
ma5>define j = g u c d s u~ c~ d~ s~
 ** ERROR: The (multi)particle 'g' is not defined.

Regards,
Shiba

Hi Shiba,

First, you need to proceed in two steps.
 1. generate the delphes root files for all your samples;
 2. analyze these delphes files.
Second, it still makes no sense to define a jet as you do at the delphes level. After hadronization, you don't have any quark and gluon left in your sample.

Cheers,

Benjamin

On 20 Jan 2016, at 06:12 , shibasipu <email address hidden> wrote:

> Question #281085 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/281085
>
> Status: Answered => Open
>
> shibasipu is still having a problem:
> Hi Benjamin,
> I am providing the following way to do my analysis in MA5.
>
> 1. First I am generating .hep files from Madgraph interfacing with
> Pythia.
>
> 2. Then I am using the following script in reco mode in MA5.
> set main.fastsim.package = delphes
> set main.fastsim.detector = atlas
> import sample/IDM/IDM200.hep.gz as sig
> import sample/SM/Wplusjets.hep as Wplusjets
> import sample/SM/Zplusjets.hep.gz as Zplusjets
> set sig.type = signal
> set Wplusjets.type = background
> set Zplusjets.type = background
> set sig.backcolor = red
> set Wplusjets.backcolor = yellow
> set Zplusjets.backcolor = green
> set sig.xsection = 0.008283
> set Wplusjets.xsection = 240.4
> set Zplusjets.xsection = 175.5
> set main.stacking_method = superimpose
> set main.stacking_method = normalize2one
> set main.lumi = 3000
> set main.normalize = lumi
> define j = g u c d s u~ c~ d~ s~
> select N(j)=4
> select MET > 160
> select PT(j[1])>130
> select PT(j[2])>60
> select PT(j[3])>60
> select PT(j[4])>60
> plot MET
> plot PT(j[1])
> plot PT(j[2])
> plot PT(j[3])
> plot PT(j[4])
> plot MT(j[1])
> plot MT(j[2])
> plot MT(j[3])
> plot MT(j[4])
> set main.SBratio = "S/sqrt(S+B)"
> submit
>
> Is this not the right way ? still I am getting the error message.
> ma5>define j = g u c d s u~ c~ d~ s~
> ** ERROR: The (multi)particle 'g' is not defined.
>
>
> Regards,
> Shiba
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Thank you very much. I need another clarification.

>>>Second, it still makes no sense to define a jet as you do at the delphes level. After hadronization, you don't have any quark and gluon left in your sample.

Suppose I have generated .hep files with only SHOWERING (without hadronization). Will it make sense for the procedures that I was following previously ?

Regards,
Shiba

Not to me. The code is not supopsed to work that way at all.

Cheers,

Benjamin

On 20 Jan 2016, at 07:47 , shibasipu <email address hidden> wrote:

> Question #281085 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/281085
>
> Status: Answered => Open
>
> shibasipu is still having a problem:
> Thank you very much. I need another clarification.
>
>>>> Second, it still makes no sense to define a jet as you do at the
> delphes level. After hadronization, you don't have any quark and gluon
> left in your sample.
>
> Suppose I have generated .hep files with only SHOWERING (without
> hadronization). Will it make sense for the procedures that I was
> following previously ?
>
>
>
> Regards,
> Shiba
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Hi Shiba,

Could you please describe exactly what you want to do, so that we could try to find a solution?

Regards,

Benjamin

On 20 Jan 2016, at 07:57 , Benjamin Fuks <email address hidden> wrote:

> Question #281085 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/281085
>
> Status: Open => Answered
>
> Benjamin Fuks proposed the following answer:
> Not to me. The code is not supopsed to work that way at all.
>
> Cheers,
>
> Benjamin
>
> On 20 Jan 2016, at 07:47 , shibasipu
> <email address hidden> wrote:
>
>> Question #281085 on MadAnalysis 5 changed:
>> https://answers.launchpad.net/madanalysis5/+question/281085
>>
>> Status: Answered => Open
>>
>> shibasipu is still having a problem:
>> Thank you very much. I need another clarification.
>>
>>>>> Second, it still makes no sense to define a jet as you do at the
>> delphes level. After hadronization, you don't have any quark and gluon
>> left in your sample.
>>
>> Suppose I have generated .hep files with only SHOWERING (without
>> hadronization). Will it make sense for the procedures that I was
>> following previously ?
>>
>>
>>
>> Regards,
>> Shiba
>>
>> --
>> You received this question notification because you are an answer
>> contact for MadAnalysis 5.
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Hi Benjamin,
                          Basically, I want to do a signal vs background simulation for 4j+MET final states. . For this, I am generating the signal and background events in Madgraph(with Pythia6 on) i.e. I am getting .hep files for signal and background.
My dominant SM backgrounds are
1. p p > W+ j j j j , ( w+ > l+ vl )
2. p p > z j j j j , ( z > vl vl~ )
3. p p > t t~ , ( t > w+ b , w+ > j j ) , (t~ > w- b~ , w- > l- vl~ )

At the end, I want a signal significance which will helpful for experimenalists. Please tell the details which I way I should follow.

Hope, you understand my point.

Regards,
Shiba

Hi,

In this case, you should not switch off hadronization in pythia. Then, you start from the stdhep files and you make use of delphes for the simulation of the detector. You get root files for each sample. Finally, you can analyze those root files and implement your selection strategy on them. I hoep this helps. You btw do not need to define any new symbol at any step.

I hope this helps.

Regards,

Benjamin

On 20 Jan 2016, at 09:57 , shibasipu <email address hidden> wrote:

> Question #281085 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/281085
>
> Status: Answered => Open
>
> shibasipu is still having a problem:
> Hi Benjamin,
> Basically, I want to do a signal vs background simulation for 4j+MET final states. . For this, I am generating the signal and background events in Madgraph(with Pythia6 on) i.e. I am getting .hep files for signal and background.
> My dominant SM backgrounds are
> 1. p p > W+ j j j j , ( w+ > l+ vl )
> 2. p p > z j j j j , ( z > vl vl~ )
> 3. p p > t t~ , ( t > w+ b , w+ > j j ) , (t~ > w- b~ , w- > l- vl~ )
>
>
> At the end, I want a signal significance which will helpful for experimenalists. Please tell the details which I way I should follow.
>
>
> Hope, you understand my point.
>
>
> Regards,
> Shiba
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Sorry. I need one clarification regarding your comment "In this case, you should not switch off hadronization in pythia". Should I switch on hadronization or switch off hadronization ? As earlier, you have mentioned that after hadronization, I do not have any jet. To be more specific, I am using the following card for Pythia.

!...Parton showering on or off
      MSTP(61)=1
      MSTP(71)=1

!...Fragmentation/hadronization on or off
      MSTJ(1)=1

!...Multiple interactions on or off
      MSTP(81)=20

!...Don't stop execution after 10 errors
      MSTU(21)=1

!...PDFset if MG set not supported by pythia-pgs package (set in lhapdf5 or higher)
      LHAID= 10041

Hope, you understand my confusion.

Regards,
Shiba

Hi Shiba,

That's very good. After pythia, you will get an event file woith tons of hadrons whose interactions with the detector will be taken care of by delphes. Then, delphes is internally calling fastet to cluster tracks and deposists into jets. Those jets are already available when you import the root file in ma5 (in the reco mode), and you can study their properties.

Cheers,

Benj

On 20 Jan 2016, at 10:22 , shibasipu <email address hidden> wrote:

> Question #281085 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/281085
>
> Status: Answered => Open
>
> shibasipu is still having a problem:
> Sorry. I need one clarification regarding your comment "In this case,
> you should not switch off hadronization in pythia". Should I switch
> on hadronization or switch off hadronization ? As earlier, you have
> mentioned that after hadronization, I do not have any jet. To be
> more specific, I am using the following card for Pythia.
>
> !...Parton showering on or off
> MSTP(61)=1
> MSTP(71)=1
>
> !...Fragmentation/hadronization on or off
> MSTJ(1)=1
>
> !...Multiple interactions on or off
> MSTP(81)=20
>
> !...Don't stop execution after 10 errors
> MSTU(21)=1
>
> !...PDFset if MG set not supported by pythia-pgs package (set in lhapdf5 or higher)
> LHAID= 10041
>
>
> Hope, you understand my confusion.
>
>
>
> Regards,
> Shiba
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Thanks for the detailed explanation. That means after importing the root file, my new script(earlier script without the definition,define j = g u c d s u~ c~ d~ s~ ) will work fine without any error message. Is not it ?

Regards,
Shiba

That's correct.

Cheers,

Benjamin

On 20 Jan 2016, at 10:42 , shibasipu <email address hidden> wrote:

> Question #281085 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/281085
>
> Status: Answered => Open
>
> shibasipu is still having a problem:
> Thanks for the detailed explanation. That means after importing the
> root file, my new script(earlier script without the definition,define j
> = g u c d s u~ c~ d~ s~ ) will work fine without any error message. Is
> not it ?
>
>
>
>
>
> Regards,
> Shiba
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Thanks Benjamin Fuks, that solved my question.

Hi Benjamin,
                               As you suggested to work with root files, I got the root files from .hep files. But the below script is not working for the root file.
set main.fastsim.package = delphes
set main.fastsim.detector = atlas
import sample/IDM/IDM200.root as sig
import sample/SM/Wplusjets.root as Wplusjets
import sample/SM/Zplusjets.root as Zplusjets
import sample/SM/ttbar_Whad_Wlep.root as ttbarWhadWlep
set sig.type = signal
set Wplusjets.type = background
set Zplusjets.type = background
set ttbarWhadWlep.type = background
set sig.backcolor = red
set Wplusjets.backcolor = yellow
set Zplusjets.backcolor = green
set ttbarWhadWlep.backcolor = purple
set sig.xsection = 0.008283
set Wplusjets.xsection = 240.4
set Zplusjets.xsection = 175.5
set ttbarWhadWlep.xsection = 0.8675
set main.stacking_method = superimpose
set main.stacking_method = normalize2one
set main.lumi = 3000
set main.normalize = lumi
#define j = g u c d s u~ c~ d~ s~
plot MET
plot THT
plot PT(j[1])
plot PT(j[2])
plot PT(j[3])
plot PT(j[4])
plot E(j[1])
plot E(j[2])
plot E(j[3])
plot E(j[4])
plot M(j[1] j[2])
plot M(j[1] j[3])
plot M(j[1] j[4])
plot M(j[2] j[3])
plot M(j[2] j[3])
plot M(j[2] j[4])
plot M(j[3] j[4])
plot DELTAR(j[1],j[2])
plot DELTAR(j[1],j[3])
plot DELTAR(j[1],j[4])
plot DELTAR(j[2],j[3])
plot DELTAR(j[2],j[4])
plot DELTAR(j[3],j[4])
plot MT(j[1])
plot MT(j[2])
plot MT(j[3])
plot MT(j[4])
plot ETA(j[1])
plot ETA(j[2])
plot ETA(j[3])
plot ETA(j[4])
set main.SBratio = "S/sqrt(S+B)"
submit

ERROR MESSAGE
=================
ma5>import sample/IDM/IDM200.root as sig
 ** ERROR: The dataset 'sample/IDM/IDM200.root' has not been found or has a unsupported format.
ma5>import sample/SM/Wplusjets.root as Wplusjets
 ** ERROR: The dataset 'sample/SM/Wplusjets.root' has not been found or has a unsupported format.
ma5>import sample/SM/Zplusjets.root as Zplusjets
 ** ERROR: The dataset 'sample/SM/Zplusjets.root' has not been found or has a unsupported format.
ma5>import sample/SM/ttbar_Whad_Wlep.root as ttbarWhadWlep
 ** ERROR: The dataset 'sample/SM/ttbar_Whad_Wlep.root' has not been found or has a unsupported format.

NOTE : Another silly question. Suppose in the above script, I want to comment a line. Which symbol should work, ! or # or % ?

Hi Shiba,

You need to proceed following two independent steps

1) you create the root files as in the example below:
   set main.fastsim.package = delphes
   set main.fastsim.detector = atlas
   import <my-hep-sample> as blablabla
   submit
The output root file will be located in the ANALYSIS_X/Output/ directory.

2) Then, you start ma5 in the reco mode and proceed with the second half of your analysis
   import sample/SM/Wplusjets.root as Wplusjets
   import sample/SM/Zplusjets.root as Zplusjets
   import sample/SM/ttbar_Whad_Wlep.root as ttbarWhadWlep
   etc…
   submit

> NOTE : Another silly question. Suppose in the above script, I want to
> comment a line. Which symbol should work, ! or # or % ?
You need to use '#' (this is python)

Cheers,

Benjamin

Hi Benjamin,
                         Thanks for the reply. But the first one is not clear to me.
1) you create the root files as in the example below:
   set main.fastsim.package = delphes
   set main.fastsim.detector = atlas
   import <my-hep-sample> as blablabla
   submit

As per my understanding, I have to write the above portion in a script file. Then where to run that script file for generating root file from .hep file ?

Regards,
Shiba

Hi,

Either in a script file, or directly within the interpreter.

Regards,

Benjamin

On 21 Jan 2016, at 18:13 , shibasipu <email address hidden> wrote:

> Question #281085 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/281085
>
> Status: Answered => Open
>
> shibasipu is still having a problem:
> Hi Benjamin,
> Thanks for the reply. But the first one is not clear to me.
> 1) you create the root files as in the example below:
> set main.fastsim.package = delphes
> set main.fastsim.detector = atlas
> import <my-hep-sample> as blablabla
> submit
>
> As per my understanding, I have to write the above portion in a script
> file. Then where to run that script file for generating root file from
> .hep file ?
>
>
> Regards,
> Shiba
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Hi Benjamin,
                            Apologizing for lengthing this thread. Actually, I was asking which mode should I use for the script.

1. To get the root file, I am using above script in reco mode of madanalysis 5. Also, I am getting the root file in Output/directory i.e. TheMouth.root. Is it correct ?

2. Another thing, Ealier I was creating the .root file from .hep file by using Delphes as a standalone via command

./DelphesSTDHEP cards/delphes_card_CMS.tcl XXXXX.root XXXXX.hep

When I am trying to use .root file obtained from Delphes (standalone), Ma5 is complaing about the format of the .root file.
ERROR MESSAGE
=================
ma5>import sample/IDM/IDM200.root as sig
 ** ERROR: The dataset 'sample/IDM/IDM200.root' has not been found or has a unsupported format.
ma5>import sample/SM/Wplusjets.root as Wplusjets
 ** ERROR: The dataset 'sample/SM/Wplusjets.root' has not been found or has a unsupported format.
ma5>import sample/SM/Zplusjets.root as Zplusjets
 ** ERROR: The dataset 'sample/SM/Zplusjets.root' has not been found or has a unsupported format.
ma5>import sample/SM/ttbar_Whad_Wlep.root as ttbarWhadWlep
 ** ERROR: The dataset 'sample/SM/ttbar_Whad_Wlep.root' has not been found or has a unsupported format.

I am not getting where is the problem for Method 2.

Regards,
Shiba

Hi Shiba,

> 1. To get the root file, I am using above script in reco mode of
> madanalysis 5. Also, I am getting the root file in Output/directory
> i.e. TheMouth.root. Is it correct ?
Yes it is.

> 2. Another thing, Ealier I was creating the .root file from .hep
> file by using Delphes as a standalone via command
>
> ./DelphesSTDHEP cards/delphes_card_CMS.tcl XXXXX.root XXXXX.hep
>
> When I am trying to use .root file obtained from Delphes (standalone), Ma5 is complaing about the format of the .root file.
> ERROR MESSAGE
> =================
> ma5>import sample/IDM/IDM200.root as sig
> ** ERROR: The dataset 'sample/IDM/IDM200.root' has not been found or has a unsupported format.
> ma5>import sample/SM/Wplusjets.root as Wplusjets
> ** ERROR: The dataset 'sample/SM/Wplusjets.root' has not been found or has a unsupported format.
> ma5>import sample/SM/Zplusjets.root as Zplusjets
> ** ERROR: The dataset 'sample/SM/Zplusjets.root' has not been found or has a unsupported format.
> ma5>import sample/SM/ttbar_Whad_Wlep.root as ttbarWhadWlep
> ** ERROR: The dataset 'sample/SM/ttbar_Whad_Wlep.root' has not been found or has a unsupported format.
>
> I am not getting where is the problem for Method 2.
Are you running madanalysis5 in the reco mode here? Could you please also try with the delphes version installed from madanalysis5?

Cheers,

Benjamin

Hi Benjamin,
                              Thanks.
1. >>Are you running madanalysis5 in the reco mode here?
Yes.

2. >>Could you please also try with the delphes version installed from madanalysis5?
The below message is coming.
[shibananda@localhost delphes]$ ./DelphesSTDHEP cards/delphes_card_CMS.tcl Wplusjets.root Wplusjets.hep
bash: ./DelphesSTDHEP: No such file or directory

Regards,
Shiba

Hi Shiba,

Well, is the executable present?

Cheers,

Benjamin

On 22 Jan 2016, at 09:32 , shibasipu <email address hidden> wrote:

> Question #281085 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/281085
>
> Status: Answered => Open
>
> shibasipu is still having a problem:
> Hi Benjamin,
> Thanks.
> 1. >>Are you running madanalysis5 in the reco mode here?
> Yes.
>
>
> 2. >>Could you please also try with the delphes version installed from madanalysis5?
> The below message is coming.
> [shibananda@localhost delphes]$ ./DelphesSTDHEP cards/delphes_card_CMS.tcl Wplusjets.root Wplusjets.hep
> bash: ./DelphesSTDHEP: No such file or directory
>
>
> Regards,
> Shiba
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Hi Benjamin,
                             None of the executables are present. But in Delphes-3.3.2, all the executables are present.

Regards,
Shiba

That means that delphes has not been properly compiled. The reason is probably due to these pythia8 libraries that do not seem to be correctly installed and prevent delphes form correctly compiling. Maybe could Eric (in cc) be more precise here?

Cheers,

Benjamin

On 22 Jan 2016, at 13:47 , shibasipu <email address hidden> wrote:

> Question #281085 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/281085
>
> Status: Answered => Open
>
> shibasipu is still having a problem:
> Hi Benjamin,
> None of the executables are present. But in Delphes-3.3.2, all the executables are present.
>
>
>
> Regards,
> Shiba
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.