On interfacing with other packages

Asked by SAUMYEN KUNDU

Hi Juergen
As mentioned in ticket number 693212, I wanted to know how can I link my Pythia8 package with Whizard. How do I re-configure this with Pythia8 and HepMC? Also if possible with Delphes (althou, this is not critically important). Do I have to reinstall this ?
Another thing, I was trying to get plots of a couple of variables with the following commands:
n_events = 1000
analysis = record e_photon (eval E [A]) ;
           record pt_photon (eval Pt [A])
histogram e_photon (0 GeV, 500 GeV, 10 GeV)
histogram pt_photon (0 GeV, 500 GeV, 10 GeV)
plot e_photon
sample_format = lhef
simulate (Bhaba)
compile_analysis {$out_file = "test.dat"}

But no plots are coming. What might be the probable reason? Is some
package missing in my system? I am though getting the .dat file and some
files in .dvi, .ps, .pdf format. Also .mp file is there. But in the PDF
or PS files no histogram is drawn there. Only the following things are
written:
Data within bounds:
hObservablei = 18.3 ± 0.90 [n entries = 1000]
All data:
hObservablei = 18.3 ± 0.90 [n entries = 1000]

Would have been good if could see the Histos.

With regards,
Saumyen

Question information

Language:
English Edit question
Status:
Solved
For:
WHIZARD Edit question
Assignee:
Juergen Reuter Edit question
Last query:
Last reply:
Whiteboard:
These are again different things, let me try to answer: yes, to include a previously not linked package, you have to reconfigure, recompile, and reinstall. Linking Pythia8 and HepMC (v3 recommended over v2) is described in the manual, Sec. 2.2.7 and 2.2.8. The internal analysis is not so trivial to use, without a complete SINDARIN it is hard to tell anything. Do you have Metapost on your system? Does the .dat file show a non-trivial histogram?
Revision history for this message
Juergen Reuter (j.r.reuter) said :
#1

These are again different things, let me try to answer: yes, to include a previously not linked package, you have to reconfigure, recompile, and reinstall. Linking Pythia8 and HepMC (v3 recommended over v2) is described in the manual, Sec. 2.2.7 and 2.2.8. The internal analysis is not so trivial to use, without a complete SINDARIN it is hard to tell anything. Do you have Metapost on your system? Does the .dat file show a non-trivial histogram?

Revision history for this message
SAUMYEN KUNDU (saumyen.k) said :
#2

Thanks a lot, Juergen, for the help.
Indeed, the Metapost package was not there in my system. So now the plots are coming.
But, I tried to re-configure with Pythia8 and HepMC but failed. I will try this again later, though.

In te meantime, (if this is the right place to ask) could you just tell me if it's possible to update this version of Whizard to a newer version. And if not, how can I uninstall this so as not to keep both the version the system.

I have another question regarding the results (errors, accuracy etc). Should I open a new ticket?

Thanks again.

With regards,
Saumyen

Revision history for this message
Juergen Reuter (j.r.reuter) said :
#3

make uninstall will uninstall WHIZARD. If you are compiling from source tar distributions, I'd recommend to unpack and build in a completely fresh build folder. Yes, different topic -> different ticket.

Revision history for this message
SAUMYEN KUNDU (saumyen.k) said :
#4

Okay, thanks a lot. I'll install in a new build folder.
Just the last question in this thread. Is it possible to plot histograms of the non-ISR and ISR photons separately like the way we put the cuts? I did this:
analysis = record e_photon (eval E [extract index 1 [sort by - Pt [A]]]) ;

But like in the cuts here I couldn't do "index > 2" to identify the non-ISR photons.

Regards,
Saumyen

Revision history for this message
Juergen Reuter (j.r.reuter) said :
#5

If you don't switch off the ISR handler you can separate matrix element photons always by Pt. Otherwise it is not possible to separate them, and also not physical.

Revision history for this message
SAUMYEN KUNDU (saumyen.k) said :
#6

Okay. So I guess I can separate withis ISR configuration that I'm using.
?keep_beams = true
!isr_order = 3
?isr_handler = true
$isr_handler_mode = "recoil"
isr_alpha = 0.0072993
isr_mass = me

Thanks a lot again.

Regards,
Saumyen

Revision history for this message
Juergen Reuter (j.r.reuter) said :
#7

Closing.