[Halld-offline] No photons from hdparsim

[Will Levine]

Hi David,

So first I made sure that I was running an unmodified version of
hdparsim and then re-ran everything from scratch, using your
jana.conf. I still see an empty histogram and no photons. My call
graph (attached) looks more or less the same as yours.

I'll take another look tomorrow to see if I can find anything else
wrong, but I'm a little mystified that we would get different results.

Thanks for all the help.

Will

On Wed, Mar 24, 2010 at 9:06 AM, David Lawrence <davidl at jlab.org> wrote:
>
> Hi Will,
>
>     I'm CC-ing the Hall-offline list since this may be of interest to
> others.
>
>     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.
> 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.
>
> To run this using the output.hddm file produced using your recipe below and
> the attached jana.conf file, just do this:
>
> hd_root --config=jana.conf output.hddm
>
>
> For documentation purposes on the mail archives, here's a copy of the
> contents of jana.conf:
>
>
> # Attach a few plugins. The janadot plugin is optional, but may be
> # useful in profiling the job and showing the dependencies of the
> # factories so you know what is calling what and how often.
> PLUGINS hdparsim,invariant_mass_hists,janadot
>
> # Use HDParSim for both charged track and photon reconstruction
> DEFTAG:DParticle          HDParSim
> DEFTAG:DTrackTimeBased    HDParSim
> DEFTAG:DPhoton            HDParSim
>
> # Tell HDParSim not to apply efficiency cut for photons
> # The efficiency cuts includes the geometric acceptance
> # so setting this to zero means we'll get "reconstructed"
> # photons that are not detectable with the real detector
> HDPARSIM:APPLY_EFFICIENCY_PHOTON 0
>
> # Tell JANA to launch as many processing threads as there are
> # cores on the system (real cores + hyperthreads). This is
> # optional.
> NTHREADS Ncores
>
>
>
>
> Will Levine wrote:
>
> 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 <davidl at jlab.org> wrote:
>
>
> 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:
>
>
> 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
> Halld-offline at jlab.org
> https://mailman.jlab.org/mailman/listinfo/halld-offline
>
>
>
> --
>
> ------------------------------------------------------------------------
> David Lawrence Ph.D.
> Staff Scientist                 Office: (757)269-5567   [[[  [   [ [
> Jefferson Lab                   Pager:  (757)584-5567   [  [ [ [ [ [
> http://www.jlab.org/~davidl     davidl at jlab.org         [[[  [[ [[ [[[
> ------------------------------------------------------------------------
>
>
>
>
>
>
>
>
> --
>
> ------------------------------------------------------------------------
>  David Lawrence Ph.D.
>  Staff Scientist                 Office: (757)269-5567   [[[  [   [ [
>  Jefferson Lab                   Pager:  (757)584-5567   [  [ [ [ [ [
>
> http://www.jlab.org/~davidl     davidl at jlab.org         [[[  [[ [[ [[[
> ------------------------------------------------------------------------
>
>
>
> # Attach a few plugins. The janadot plugin is optional, but may be
> # useful in profiling the job and showing the dependencies of the
> # factories so you know what is calling what and how often.
> PLUGINS hdparsim,invariant_mass_hists,janadot
>
> # Use HDParSim for both charged track and photon reconstruction
> DEFTAG:DParticle          HDParSim
> DEFTAG:DTrackTimeBased    HDParSim
> DEFTAG:DPhoton            HDParSim
>
> # Tell HDParSim not to apply efficiency cut for photons
> # The efficiency cuts includes the geometric acceptance
> # so setting this to zero means we'll get "reconstructed"
> # photons that are not detectable with the real detector
> HDPARSIM:APPLY_EFFICIENCY_PHOTON 0
>
> # Tell JANA to launch as many processing threads as there are
> # cores on the system (real cores + hyperthreads). This is
> # optional.
> NTHREADS Ncores
>
>
>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: jana.pdf
Type: application/pdf
Size: 3831 bytes
Desc: not available
URL: <https://mailman.jlab.org/pipermail/halld-offline/attachments/20100324/d48e1cb7/attachment-0002.pdf>