[Halld-offline] No photons from hdparsim

[Will Levine]

Ah okay, I figured out what the problem was. I forgot to change my
PATH to point to the new bin directory in sim-recon. Sorry for the
trouble.

Will

On Wed, Mar 24, 2010 at 9:50 PM, Will Levine <wilevine at andrew.cmu.edu> wrote:
> 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
>>
>>
>>
>