issues with Delphes data at ma5

[Mouchemhou]

Question

Hi,
I am generating my signal and background within MG5+pythia8+delphes within cmd defaults card. but when I am tryin o analysis this data at ma5, I find some probleme like:

MA5-ERROR: The dataset '/home/mohamed/MG5_aMC_v3_5_5/BP1_Delphes/Events/runBPn_decayed_1/tag_1_delphes_events.root' has not been found or has a unsupported format.

MA5: ma5>set BP1.xsection = 0.01794491279
MA5-ERROR: no object called 'BP1' is found

MA5: ma5>set BP1.weighted_events = true
MA5-ERROR: no object called 'BP1' is found
Also, I have issues on how to apply Btagging in ma5
MA5: ma5>select N(Jet where Jet.BTag == 1) >= 2
MA5-ERROR: 'Jet' is not a defined (multi)particle.

this is the script used for run ma5 in reco mode
"set main.fastsim.package = delphesMA5tune
set main.fastsim.detector = cms
set main.stacking_method = normalize2one
set main.lumi = 300.
set main.normalize = lumi
set main.fom.formula = 6

import /home/mohamed/MG5_aMC_v3_5_5/BP1_Delphes/Events/runBPn_decayed_1/tag_1_delphes_events.root as BP1
set BP1.xsection = 0.01794491279
set BP1.weighted_events = true"

Any idea please?
Best regards, MO

[BFuks]

Hi Mouchemhou.

I see two problems:

• your event file does not exist or is not located where you think it is.
• the syntax is wrong (please check the manual or the reference cards).

Best regards,

Benjamin

[Mouchemhou]

Dear Benjamin,
The data is in this locations,
it work if i remove "set main.fastsim.package = delphesMA5tune ## or delphes
set main.fastsim.detector = cms"

best regards, MO

[BFuks]

Hi Mouchemhou,

I am a bit confused. Can you please share with me the exact script you have used, and the problems you got. Also, how do you start the code ? In reco mode (bin/ma5 - R) ?

Cheers,

Benjamin

[Mouchemhou]

Dear Benjamin,

this is the scripts that I used "#set nb_core 4
#set main.fastsim.package = delphes
#set main.fastsim.detector = cms
#set main.fastsim.package = fastjet
#set main.fastsim.algorithm = antikt
#set main.lumi = 300.
#set main.fastsim.radius = 0.4
#set main.outputfile = output.lhe

set main.fastsim.package = delphesMA5tune
set main.fastsim.detector = cms
set main.stacking_method = normalize2one
set main.lumi = 300.
set main.normalize = lumi
set main.fom.formula = 6

import /home/mohamed/MG5_aMC_v3_5_5/BP1_Delphes/Events/runBPn_decayed_1/tag_1_delphes_events.root as BP1

set BP1.xsection = 0.01794491279
set BP1.weighted_events = true

select N(Jet where Jet.BTag == 1) >= 2

Distributions Cuts

#select N(b)>= 4 and N(l)==1
#select PT(b)>= 15
#select PT(l)>= 15
#select PT(j)>= 10
#select MET >=15. and MHT>=15.
#select M (b[1] b[2]) >=5. and M (b[1] b[2])<=30.
#select M (b[3] b[4]) >=5. and M (b[3] b[4])<=20.
#select M (b[1] b[2] b[3] b[4])>=20. and M (b[1] b[2] b[3] b[4])<=55.
#select MT_MET(l)<=100.0
#select M (b[1] b[2] b[3] b[4] l[1]) <=150.0
#select MT (b[1] b[2] b[3] b[4] l[1]) <=200.

===== Distributions's plot

plot MET 25 0 400
plot MHT 25 0 400

plot TET 50 0 1000
plot THT 50 0 1000

plot PT(b[1]) 25 0 100 [logY]
plot PT(b[2]) 25 0 100 [logY]
plot PT(b[3]) 25 0 100 [logY]
plot PT(b[4]) 25 0 100 [logY]

plot PT(l[1]) 25 0 100 [logY]

plot PT(j[1]) 25 0 100 [logY]
plot PT(j[2]) 25 0 100 [logY]

plot ETA(b) 25 -3 3
plot ETA(l) 25 -3 3
plot ETA(j) 25 -3 3

plot M (b[1] b[2]) 25 0 120 [logY]
plot M (b[3] b[4]) 25 0 100 [logY]

plot M (b[1] b[2] b[3] b[4]) 25 0 140 [logY]
plot M (b[1] b[2] b[3] b[4] l[1]) 25 0 400 [logY]

plot MT_MET(l) 25 0 100 [logY]
plot MT(b[1] b[2] b[3] b[4] l[1]) 25 0 500 [logY]
plot MT_MET(b[1] b[2] b[3] b[4] l[1]) 25 0 500 [logY]

set BP1.type = signal

submit BP1_Optimised_wp2
exit

". but the problem I got, error like :
MA5: ma5>import /home/mohamed/MG5_aMC_v3_5_5/BP1_Delphes/Events/runBPn_decayed_1/tag_1_delphes_events.root as BP1
MA5-ERROR: The dataset '/home/mohamed/MG5_aMC_v3_5_5/BP1_Delphes/Events/runBPn_decayed_1/tag_1_delphes_events.root' has not been found or has a unsupported format.
MA5: ma5>set BP1.xsection = 0.01794491279
MA5-ERROR: no object called 'BP1' is found
MA5: ma5>set BP1.weighted_events = true
MA5-ERROR: no object called 'BP1' is found
MA5: ma5>select N(Jet where Jet.BTag == 1) >= 2
MA5-ERROR: 'Jet' is not a defined (multi)particle.

and Even that If I remove MA5: ma5>select N(Jet where Jet.BTag == 1) >= 2. the probleme persist.
I think the problem is on the
"set main.fastsim.package = delphesMA5tune
set main.fastsim.detector = cms
set main.stacking_method = normalize2one
set main.lumi = 300.
set main.normalize = lumi
set main.fom.formula = 6
"

best regards, MO

[BFuks]

Hi, The problem is that you are trying to run delphes on a sample that has already been proceesed by delphes. The input must be a hepmc file. Regards, Benjamin

…

On 23 Jun 2025, at 17:31, Mohamed Ouchemhou ***@***.***> wrote: Mouchemhou left a comment (MadAnalysis/madanalysis5#291) <#291 (comment)> Dear Benjamin, this is the scripts that I used "#set nb_core 4 #set main.fastsim.package = delphes #set main.fastsim.detector = cms #set main.fastsim.package = fastjet #set main.fastsim.algorithm = antikt #set main.lumi = 300. #set main.fastsim.radius = 0.4 #set main.outputfile = output.lhe set main.fastsim.package = delphesMA5tune set main.fastsim.detector = cms set main.stacking_method = normalize2one set main.lumi = 300. set main.normalize = lumi set main.fom.formula = 6 import /home/mohamed/MG5_aMC_v3_5_5/BP1_Delphes/Events/runBPn_decayed_1/tag_1_delphes_events.root as BP1 set BP1.xsection = 0.01794491279 set BP1.weighted_events = true select N(Jet where Jet.BTag == 1) >= 2 Distributions Cuts #select N(b)>= 4 and N(l)==1 #select PT(b)>= 15 #select PT(l)>= 15 #select PT(j)>= 10 #select MET >=15. and MHT>=15. #select M (b[1] b[2]) >=5. and M (b[1] b[2])<=30. #select M (b[3] b[4]) >=5. and M (b[3] b[4])<=20. #select M (b[1] b[2] b[3] b[4])>=20. and M (b[1] b[2] b[3] b[4])<=55. #select MT_MET(l)<=100.0 #select M (b[1] b[2] b[3] b[4] l[1]) <=150.0 #select MT (b[1] b[2] b[3] b[4] l[1]) <=200. ===== Distributions's plot plot MET 25 0 400 plot MHT 25 0 400 plot TET 50 0 1000 plot THT 50 0 1000 plot PT(b[1]) 25 0 100 [logY] plot PT(b[2]) 25 0 100 [logY] plot PT(b[3]) 25 0 100 [logY] plot PT(b[4]) 25 0 100 [logY] plot PT(l[1]) 25 0 100 [logY] plot PT(j[1]) 25 0 100 [logY] plot PT(j[2]) 25 0 100 [logY] plot ETA(b) 25 -3 3 plot ETA(l) 25 -3 3 plot ETA(j) 25 -3 3 plot M (b[1] b[2]) 25 0 120 [logY] plot M (b[3] b[4]) 25 0 100 [logY] plot M (b[1] b[2] b[3] b[4]) 25 0 140 [logY] plot M (b[1] b[2] b[3] b[4] l[1]) 25 0 400 [logY] plot MT_MET(l) 25 0 100 [logY] plot MT(b[1] b[2] b[3] b[4] l[1]) 25 0 500 [logY] plot MT_MET(b[1] b[2] b[3] b[4] l[1]) 25 0 500 [logY] set BP1.type = signal submit BP1_Optimised_wp2 exit ". but the problem I got, error like : MA5: ma5>import /home/mohamed/MG5_aMC_v3_5_5/BP1_Delphes/Events/runBPn_decayed_1/tag_1_delphes_events.root as BP1 MA5-ERROR: The dataset '/home/mohamed/MG5_aMC_v3_5_5/BP1_Delphes/Events/runBPn_decayed_1/tag_1_delphes_events.root' has not been found or has a unsupported format. MA5: ma5>set BP1.xsection = 0.01794491279 MA5-ERROR: no object called 'BP1' is found MA5: ma5>set BP1.weighted_events = true MA5-ERROR: no object called 'BP1' is found MA5: ma5>select N(Jet where Jet.BTag == 1) >= 2 MA5-ERROR: 'Jet' is not a defined (multi)particle. and Even that If I remove MA5: ma5>select N(Jet where Jet.BTag == 1) >= 2. the probleme persist. I think the problem is on the "set main.fastsim.package = delphesMA5tune set main.fastsim.detector = cms set main.stacking_method = normalize2one set main.lumi = 300. set main.normalize = lumi set main.fom.formula = 6 " best regards, MO — Reply to this email directly, view it on GitHub <#291 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AIWJNJ24RCRINR433OD4XC33FAMW5AVCNFSM6AAAAAB75GC5UGVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDSOJWHE2DENJVGA>. You are receiving this because you commented.

[Mouchemhou]

Thank you, dear Benjamin.
I had a doubt, but you’ve now confirmed it for me.

Best regards,
MO