<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
<br>
Hi Will,<br>
<br>
I'm CC-ing the Hall-offline list since this may be of interest to
others.<br>
<br>
I've run things as you've described, and I am seeing a significant
number of photons get reconstructed. There are some that are "lost" due
to efficiency cuts HDParSim applies using tables from the same root
files that supply the resolutions. With recent versions of HDParSim
(rev6122 or greater circa 2010-03-10) you can use the
HDPARSIM:APPLY_EFFICIENCY_PHOTON configuration parameter to turn off
the efficiency cut.<br>
I've attached a configuration file that uses this along with the
resulting plot of the mass_2gamma histogram to see if your installation
reproduces this. Also, find attached a call graph made using the
janadot plugin which shows which factories are being called and by
whom. Note that I did uncomment the line in
DEventProcessor_invariant_mass_hists.cc in order to get the photons as
you pointed out. I went ahead and checked that change into the
repository.<br>
<br>
To run this using the output.hddm file produced using your recipe below
and the attached jana.conf file, just do this:<br>
<br>
hd_root --config=jana.conf output.hddm<br>
<br>
<br>
For documentation purposes on the mail archives, here's a copy of the
contents of jana.conf:<br>
<br>
<br>
<tt># Attach a few plugins. The janadot plugin is optional, but may be <br>
# useful in profiling the job and showing the dependencies of the<br>
# factories so you know what is calling what and how often.<br>
PLUGINS hdparsim,invariant_mass_hists,janadot<br>
<br>
# Use HDParSim for both charged track and photon reconstruction<br>
DEFTAG:DParticle HDParSim<br>
DEFTAG:DTrackTimeBased HDParSim<br>
DEFTAG:DPhoton HDParSim<br>
<br>
# Tell HDParSim not to apply efficiency cut for photons<br>
# The efficiency cuts includes the geometric acceptance<br>
# so setting this to zero means we'll get "reconstructed"<br>
# photons that are not detectable with the real detector<br>
HDPARSIM:APPLY_EFFICIENCY_PHOTON 0<br>
<br>
# Tell JANA to launch as many processing threads as there are<br>
# cores on the system (real cores + hyperthreads). This is<br>
# optional.<br>
NTHREADS Ncores<br>
<br>
</tt><br>
<br>
<br>
Will Levine wrote:
<blockquote
cite="mid:d1cf954b1003231900u5b1ea07n6352f74cbb052dad@mail.gmail.com"
type="cite">
<pre wrap="">Hi David,
Yes, I am using hdparsim to produce DPhotons.
For genr8, I am just using the b1_pi.input file from
sim-recon/src/programs/Simulation/genr8/InputFiles . I'm using my own
plugin for analysis, but I see the same thing with the
invariant_mass_hists plugin (although I need to be careful here, since
for some reason the current version of this file has the line
loop->Get(photons);
commented out, which, if left commented, would of course leave me
without photons)
Anyway, you should be able to see the problem with something simple like this:
genr8 -Ab1_pi.asc -M8200 < b1_pi.input
genr8_2_hddm b1_pi.asc
hd_root --plugin=hdparsim --plugin=invariant_mass_hists
-PDEFTAG:DParticle=HDParSim -PDEFTAG:DPhoton=HDParSim output.hddm
Each event should have two photons associated with it, but the
histograms related to photons are unfilled, and furthermore, there
don't seem to be any DPhotons from hdparsim at all.
Thanks,
Will
On Tue, Mar 23, 2010 at 4:50 PM, David Lawrence <a class="moz-txt-link-rfc2396E" href="mailto:davidl@jlab.org"><davidl@jlab.org></a> wrote:
</pre>
<blockquote type="cite">
<pre wrap="">Hi Will,
Can you give me a little more information. Specifically:
1.) Are you using hdparsim to produce the "reconstructed" DPhotons (i.e.
using -PDEFTAG:DPhoton=HDParSim)
or are you reconstructing the DPhotons and using HDParSim just for charged
particles.
2.) Please send me your configuration file for genr8 and the exact commands
you're using so I can try and reproduce the error.
Regards,
-David
Will Levine wrote:
</pre>
<blockquote type="cite">
<pre wrap="">Hi all,
I just upgraded my Hall-D software to the latest version and now in
some cases, I do not get any DPhotons from hdparsim. Things work as
expected when I run hdparsim on an hddm file generated by bggen, but
not with an hddm file from genr8, which had worked previously. I know
that Blake posted about a similar-sounding problem earlier, but after
reading through that thread I'm not sure if this is the same issue.
Will
_______________________________________________
Halld-offline mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Halld-offline@jlab.org">Halld-offline@jlab.org</a>
<a class="moz-txt-link-freetext" href="https://mailman.jlab.org/mailman/listinfo/halld-offline">https://mailman.jlab.org/mailman/listinfo/halld-offline</a>
</pre>
</blockquote>
<pre wrap="">--
------------------------------------------------------------------------
David Lawrence Ph.D.
Staff Scientist Office: (757)269-5567 [[[ [ [ [
Jefferson Lab Pager: (757)584-5567 [ [ [ [ [ [
<a class="moz-txt-link-freetext" href="http://www.jlab.org/~davidl">http://www.jlab.org/~davidl</a> <a class="moz-txt-link-abbreviated" href="mailto:davidl@jlab.org">davidl@jlab.org</a> [[[ [[ [[ [[[
------------------------------------------------------------------------
</pre>
</blockquote>
<pre wrap=""><!---->
</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
------------------------------------------------------------------------
David Lawrence Ph.D.
Staff Scientist Office: (757)269-5567 [[[ [ [ [
Jefferson Lab Pager: (757)584-5567 [ [ [ [ [ [
<a class="moz-txt-link-freetext" href="http://www.jlab.org/~davidl">http://www.jlab.org/~davidl</a> <a class="moz-txt-link-abbreviated" href="mailto:davidl@jlab.org">davidl@jlab.org</a> [[[ [[ [[ [[[
------------------------------------------------------------------------
</pre>
</body>
</html>