-
Notifications
You must be signed in to change notification settings - Fork 14
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Merge pull request #1151 from JeffersonLab/elton_2pi_primakoff
Elton 2pi primakoff
- Loading branch information
Showing
6 changed files
with
141 additions
and
58 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -1,33 +1,2 @@ | ||
| # Elton Smith. Updated June 7, 2017. Mod for Z2pi_trees | ||
| # Here is the sequence of instructions to generate, fit and view generated distributions | ||
|
|
||
| gen_2pi_primakoff> gen_2pi_primakoff -c gen_2pi_primakoff_signal.cfg -o tree_gen_2pi_primakoff_signal.root -hd gen_2pi_primakoff_signal.hddm -a 5.5 -b 6.0 -p 6.0 -m 11.6 -n 100000 -r 31000 | ||
| gen_2pi_primakoff> gen_2pi_primakoff -c gen_2pi_primakoff_flat.cfg -o tree_gen_2pi_primakoff_flat.root -hd gen_2pi_primakoff_flat.hddm -a 5.5 -b 6.0 -p 6.0 -m 11.6 -n 100000 -r 31000 | ||
|
|
||
| # Here are instructions for processing MC smeared output files / or data | ||
|
|
||
| gen_2pi_primakoff> hd_root -PPLUGINS=monitoring_hists,Z2pi_trees -PNTHREADS=4 -PEVENTS_TO_KEEP=10000 ../hddm/dana_rest_gen_2pi_primakoff_030000_00*.hddm -o hd_root_Z2pi_trees_signal.root | ||
| gen_2pi_primakoff> mv tree_Z2pi_trees.root tree_hd_root_Z2pi_trees_signal.root | ||
|
|
||
| gen_2pi_primakoff> root -l -b tree_hd_root_Z2pi_trees_signal.root | ||
| root [0] .x $ROOT_ANALYSIS_HOME/scripts/Load_DSelector.C | ||
| root [0] Z2pi_trees_Tree->Process("DSelector_Z2pi_trees.C+"); | ||
| root [0] .q; | ||
| gen_2pi_primakoff> mv DSelector_Z2pi_trees.root DSelector_Z2pi_trees_signal.root | ||
| gen_2pi_primakoff> mv tree_DSelector_Z2pi_trees.root tree_DSelector_Z2pi_trees_signal.root | ||
| gen_2pi_primakoff> root -l | ||
| root [0].x plot_Z2pi_trees.C("DSelector_Z2pi_trees_signal") | ||
| root [0].q | ||
|
|
||
| gen_2pi_primakoff> tree_to_amptools tree_DSelector_Z2pi_trees_signal.root Z2pi_trees_Tree | ||
| gen_2pi_primakoff> mv AmpToolsInputTree.root tree_DSelector_Z2pi_trees_signal_amptools.root | ||
|
|
||
| gen_2pi_primakoff> fit -c fit_2pi_primakoff.cfg | ||
| gen_2pi_primakoff> cp twopi_primakoff.fit twopi_primakoff_DSelect.fit | ||
| gen_2pi_primakoff> twopi_plotter_primakoff twopi_primakoff_DSelect.fit -o twopi_primakoff_DSelect.root | ||
| gen_2pi_primakoff> mv twopi_fitPars.txt twopi_primakoff_DSelect.fit2 | ||
| gen_2pi_primakoff> root -l | ||
| root [0] .x twopi_primakoff.C | ||
|
|
||
|
|
||
|
|
||
| # Elton Smith. Updated June 14, 2018. Mod for Z2pi_trees | ||
| # Use the script streamline2.csh to execute commands to process CPP MC through Amplitude Analysis |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
114 changes: 114 additions & 0 deletions
114
src/plugins/Analysis/Z2pi_trees/scripts/fit_2pi_primakoff_Intime_1000000.cfg
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,114 @@ | ||
| # | ||
| ##################################### | ||
| #### THIS IS A CONFIG FILE #### | ||
| ##################################### | ||
| ## | ||
| ## Blank lines or lines beginning with a "#" are ignored. | ||
| ## | ||
| ## Double colons (::) are treated like a space. | ||
| ## This is sometimes useful for grouping (for example, | ||
| ## grouping strings like "reaction::sum::amplitudeName") | ||
| ## | ||
| ## All non-comment lines must begin with one of the following keywords. | ||
| ## | ||
| ## (note: <word> means necessary | ||
| ## (word) means optional) | ||
| ## | ||
| ## include <file> | ||
| ## define <word> (defn1) (defn2) (defn3) ... | ||
| ## fit <fitname> | ||
| ## keyword <keyword> <min arguments> <max arguments> | ||
| ## reaction <reaction> <particle1> <particle2> (particle3) ... | ||
| ## data <reaction> <class> (arg1) (arg2) (arg3) ... | ||
| ## genmc <reaction> <class> (arg1) (arg2) (arg3) ... | ||
| ## accmc <reaction> <class> (arg1) (arg2) (arg3) ... | ||
| ## normintfile <reaction> <file> | ||
| ## sum <reaction> <sum> (sum2) (sum3) ... | ||
| ## amplitude <reaction> <sum> <amp> <class> (arg1) (arg2) ([par]) ... | ||
| ## initialize <reaction> <sum> <amp> <"events"/"polar"/"cartesian"> | ||
| ## <value1> <value2> ("fixed"/"real") | ||
| ## scale <reaction> <sum> <amp> <value or [parameter]> | ||
| ## constrain <reaction1> <sum1> <amp1> <reaction2> <sum2> <amp2> ... | ||
| ## permute <reaction> <sum> <amp> <index1> <index2> ... | ||
| ## parameter <par> <value> ("fixed"/"bounded"/"gaussian") | ||
| ## (lower/central) (upper/error) | ||
| ## DEPRECATED: | ||
| ## datafile <reaction> <file> (file2) (file3) ... | ||
| ## genmcfile <reaction> <file> (file2) (file3) ... | ||
| ## accmcfile <reaction> <file> (file2) (file3) ... | ||
| ## | ||
| ##################################### | ||
|
|
||
| #Take useful definitions from gen_3pi.cfg | ||
| define rho 0.775 0.146 | ||
| # some definitions for adjusting the beam polarization. In the equations beamX: Polarization>0, beamY: Polarization<0. | ||
| # Phi is the azimuthal angle of the polarization vector in degrees (in the lab coordinates) | ||
| define phipol 0 | ||
| define polFrac 0.7 | ||
| define beamX 0 polFrac | ||
| define beamY 90 polFrac | ||
|
|
||
| # Uniform angles: flat=1; YLMs: flat=0; | ||
| define flat 0 | ||
|
|
||
| fit twopi_primakoff | ||
|
|
||
| reaction Primakoff gamma Pi+ Pi- Pb208 | ||
|
|
||
| normintfile Primakoff twopi_primakoff_ni.txt | ||
|
|
||
| # sum contains two terms for s-wave production of pi+pi- | ||
| sum Primakoff Aplus | ||
| sum Primakoff Aminus | ||
|
|
||
| # genmc Primakoff ROOTDataReader tree_gen_2pi_primakoff_flat_100000.root | ||
| # accmc Primakoff ROOTDataReader tree_gen_2pi_primakoff_flat_100000.root | ||
| # data Primakoff ROOTDataReader tree_gen_2pi_primakoff_signal_100000.root | ||
| genmc Primakoff ROOTDataReader treeFlat_gen_2pi_primakoff_sw1pw1000_NOTAG_flat_1000000_amptools.root | ||
| accmc Primakoff ROOTDataReader treeFlat_DSelector_Z2pi_trees_sw1pw1000_NOTAG_flat_1000000_amptools.root | ||
| data Primakoff ROOTDataReader treeFlat_DSelector_Z2pi_trees_sw1pw1000_TAG10_signal_1000000_amptools_InTimeW.root | ||
|
|
||
| # amplitude Primakoff::Aplus TwoPiAngles_primakoff (phipol, pol fraction, m_rho, PhaseFactor, flat) | ||
| # phipol is the lab azimuthal angle of the polarization vector. | ||
| # pol fraction is the linear polarization of the beam | ||
| # m_rho Jz component of rho | ||
| # Phasefactor determines prefix factor to amplitudes in computation | ||
| # flat=1 generates distribution uniform in angles. flat=0 use YLMs | ||
|
|
||
| amplitude Primakoff::Aplus::g1V00 TwoPiAngles_primakoff phipol polFrac 0 0 flat | ||
| amplitude Primakoff::Aplus::g1V00 TwoPiWt_primakoff 1.29 0. 1 2 | ||
| amplitude Primakoff::Aplus::g1V11 TwoPiAngles_primakoff phipol polFrac 1 2 flat | ||
| amplitude Primakoff::Aplus::g1V11 BreitWigner rho 1 1 2 | ||
| amplitude Primakoff::Aplus::g1V10 TwoPiAngles_primakoff phipol polFrac 0 2 flat | ||
| amplitude Primakoff::Aplus::g1V10 BreitWigner rho 1 1 2 | ||
| amplitude Primakoff::Aplus::g1V1-1 TwoPiAngles_primakoff phipol polFrac -1 2 flat | ||
| amplitude Primakoff::Aplus::g1V1-1 BreitWigner rho 1 1 2 | ||
|
|
||
| amplitude Primakoff::Aminus::g1V00 TwoPiAngles_primakoff phipol polFrac 0 1 flat | ||
| amplitude Primakoff::Aminus::g1V00 TwoPiWt_primakoff 1.29 0. 1 2 | ||
| amplitude Primakoff::Aminus::g1V11 TwoPiAngles_primakoff phipol polFrac 1 3 flat | ||
| amplitude Primakoff::Aminus::g1V11 BreitWigner rho 1 1 2 | ||
| amplitude Primakoff::Aminus::g1V10 TwoPiAngles_primakoff phipol polFrac 0 3 flat | ||
| amplitude Primakoff::Aminus::g1V10 BreitWigner rho 1 1 2 | ||
| amplitude Primakoff::Aminus::g1V1-1 TwoPiAngles_primakoff phipol polFrac -1 3 flat | ||
| amplitude Primakoff::Aminus::g1V1-1 BreitWigner rho 1 1 2 | ||
|
|
||
| initialize Primakoff::Aplus::g1V00 cartesian 10.0 0.0 | ||
| initialize Primakoff::Aplus::g1V11 cartesian 2000.0 0.0 real | ||
| initialize Primakoff::Aplus::g1V10 cartesian 0.0 0.0 | ||
| initialize Primakoff::Aplus::g1V1-1 cartesian 0.0 0.0 | ||
|
|
||
| initialize Primakoff::Aminus::g1V00 cartesian 10.0 0.0 | ||
| initialize Primakoff::Aminus::g1V11 cartesian 2000.0 0.0 | ||
| initialize Primakoff::Aminus::g1V10 cartesian 0.0 0.0 | ||
| initialize Primakoff::Aminus::g1V1-1 cartesian 0.0 0.0 | ||
|
|
||
| constrain Primakoff::Aplus::g1V00 Primakoff::Aminus::g1V00 | ||
| constrain Primakoff::Aplus::g1V11 Primakoff::Aminus::g1V11 | ||
| constrain Primakoff::Aplus::g1V10 Primakoff::Aminus::g1V10 | ||
| constrain Primakoff::Aplus::g1V1-1 Primakoff::Aminus::g1V1-1 | ||
|
|
||
|
|
||
|
|
||
|
|
||
|
|
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters