Crashes in BeamProperties

[Elton Smith]

Dear Software Help,

I am trying to access the energy-dependent polarization to add to my 2pi0 simulation. Looking at other examples, it seems like they should be easy to access using the BeamProperties class. I have copied code that I would think should work, but the program crashes hard each time. I don’t know if there are initializations that are needed that I am not doing. Any suggestions are welcome.

I am working in the npp_ana branch of halld_sim in the following directories

Code: /w/halld-scshelf2101/halld2/home/elton/git/halld_sim/src/programs/Simulation/gen_2pi0_primakoff

Running program: /w/halld-scshelf2101/halld2/home/elton/gen_2pi0_primakoff

Below are files with the program output and also the crash log.  

The crashing occurs in TwoPiAngles_primakoff.cc at about line 33.

 https://github.com/JeffersonLab/halld_sim/compare/master...npp_sim#diff-e10cbeb08afd68c703a894142d70044c6de03506e70d34e4c223fa647b20df2d

 BeamProperties beamProp(beamConfigFile);

 polFrac_vs_E = (TH1D*)beamProp.GetPolFrac();

 polAngle = beamProp.GetPolAngle();

Thanks, Elton.

------------

Elton Smith

Elton Smith

el...@jlab.org

Cell: 757-749-5262

[Alexander Austregesilo]

Hi Elton,

The BeamProperties class is not supported within the amptools libraries. Instead, you should directly read in the TPOL histograms from the root file. You can see an example in the SDME amplitude:

https://github.com/JeffersonLab/halld_sim/blob/ff41f5fe12c264880b557fb1b66f22bcfc1d8b58/src/libraries/AMPTOOLS_AMPS/TwoPiAngles.cc#L60

Best regards,

Alex

--
You received this message because you are subscribed to the Google Groups "GlueX Software Help" group.
To unsubscribe from this group and stop receiving emails from it, send an email to gluex-softwar...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/gluex-software/30BC3F2E-C21B-419E-AE67-4FA0F5E26F3E%40jlab.org.

-- 
Alexander Austregesilo

Staff Scientist - Experimental Nuclear Physics
Thomas Jefferson National Accelerator Facility
Newport News, VA
aaus...@jlab.org
(757) 269-6982

[Elton Smith]

Well that explains it….. Thanks for your prompt response.  I will look at your example.

Elton Smith

el...@jlab.org

Cell: 757-749-5262

To view this discussion visit https://groups.google.com/d/msgid/gluex-software/719cdc35-e043-46fa-9342-680be18539cb%40jlab.org.

[Justin Stevens]

Also, FYI the TPOL polarization information is documented here https://halldweb.jlab.org/wiki-private/index.php/TPOL_Polarization but I’m not sure if it yet includes the CPP/NPP datasets to produce histograms of this polarization vs beam energy.  For the GlueX data, I know there are people using it in AmpTools fits, for example Hao Li (ha...@jlab.org).

-Justin

To view this discussion visit https://groups.google.com/d/msgid/gluex-software/F9D6CA07-426F-4F87-A987-CBE9133E290E%40jlab.org.

[Alexander Austregesilo]

The TPOL analysis for CPP is documented here:

 https://userweb.jlab.org/~dugger/triPol/TPOL_2022-05.pdf

I'll update the wiki page.

- Alex

[Elton Smith]

Hi Alex and Justin,

Thanks for your comments and updating the web page and to Mike Dugger with clear and easy to follow instructions. Worked like a charm!

I will work with Alex’s example to use it in my generator.

Cheers, Elton.

Elton Smith

el...@jlab.org

Cell: 757-749-5262

[Elton Smith]

Hi Alex,

With your example, I now have a test program that reads in the ROOT file with TPOL polarizations. I copied the ROOT file into the directory where MCWrapper was running locally. This was ok for testing, but clearly just a kludge.

What is the best way of making the ROOT file accessible to MCWrapper? I guess would might be able  use the ANA_ENVIRONMENT_FILE option for use in batch mode, but I don’t think the file is sourced interactively. How do you do it for the rho analysis?

Thanks, Elton

Elton Smith

el...@jlab.org

Cell: 757-749-5262