From 41bafcc6da77f86a8461624bb2d7b965f9a7d942 Mon Sep 17 00:00:00 2001 From: SusanneGlaessel Date: Tue, 22 Apr 2025 18:34:56 +0200 Subject: [PATCH 1/4] small change in parfile --- at_interface/example_parfiles/parfile2_add_omega.txt | 2 +- at_interface/example_parfiles/parfile2_add_xi.txt | 4 ++-- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/at_interface/example_parfiles/parfile2_add_omega.txt b/at_interface/example_parfiles/parfile2_add_omega.txt index dff542c..2d2e4f4 100644 --- a/at_interface/example_parfiles/parfile2_add_omega.txt +++ b/at_interface/example_parfiles/parfile2_add_omega.txt @@ -2,7 +2,7 @@ 3334, pdg mother Omega name mother -1, mass mother: = -1: pdg mass taken from KFParticleDatabase --1, sigma mother: = -1: pdg mass sigma taken from KFParticleDatabase +-1, mass sigma mother: = -1: pdg mass sigma taken from KFParticleDatabase 2, number of daughters -321, pdg daughter 1 0, optional: additional pdg to be considered for daughter 1 (=0: no additional pdg) diff --git a/at_interface/example_parfiles/parfile2_add_xi.txt b/at_interface/example_parfiles/parfile2_add_xi.txt index 08553a9..2b13752 100644 --- a/at_interface/example_parfiles/parfile2_add_xi.txt +++ b/at_interface/example_parfiles/parfile2_add_xi.txt @@ -1,8 +1,8 @@ ## Decay 3312, pdg mother -Xi name mother +Xi- name mother -1, mass mother: = -1: pdg mass taken from KFParticleDatabase --1, sigma mother: = -1: pdg mass sigma taken from KFParticleDatabase +-1, mass sigma mother: = -1: pdg mass sigma taken from KFParticleDatabase 2, number of daughters -211, pdg daughter 1 -1, optional: additional pdg to be considered for daughter 1 (=0: no additional pdg) From 9e1622e3ccacaebd6d3e1c11cb62eec676ac8002 Mon Sep 17 00:00:00 2001 From: SusanneGlaessel Date: Thu, 24 Apr 2025 11:35:08 +0200 Subject: [PATCH 2/4] change in parfile --- .../example_parfiles/parfile2_lambda_cascade_pidmode1.txt | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/at_interface/example_parfiles/parfile2_lambda_cascade_pidmode1.txt b/at_interface/example_parfiles/parfile2_lambda_cascade_pidmode1.txt index acd00b8..867f3f6 100644 --- a/at_interface/example_parfiles/parfile2_lambda_cascade_pidmode1.txt +++ b/at_interface/example_parfiles/parfile2_lambda_cascade_pidmode1.txt @@ -23,9 +23,9 @@ Lambda name mother -1, cut cosopen -------------------------------------------------------------------------- ## Geometrical cuts for mother --1, cut chi2topo +5.0, cut chi2topo -1, cut costopo -5.0, cut LdL +10.0, cut LdL -1, cut decaylength L -1, cut distance to PV line 3.0, cut invariant mass discrepancy in terms of expected width From b9cd9d5ea454112ba9d9da978d23614419db3840 Mon Sep 17 00:00:00 2001 From: SusanneGlaessel Date: Fri, 25 Apr 2025 19:47:23 +0200 Subject: [PATCH 3/4] bugfix in main.cpp --- README.md | 6 +++--- at_interface/ConverterIn.cpp | 2 +- at_interface/example_parfiles/parfile3_H3L_pidmode4.txt | 4 ++-- at_interface/main.cpp | 4 ++-- 4 files changed, 8 insertions(+), 8 deletions(-) diff --git a/README.md b/README.md index 95655e3..2773bc2 100644 --- a/README.md +++ b/README.md @@ -138,19 +138,19 @@ The user can select: - Output format: default output is analysistree format optional: plain tree: output is a simple root tree containing the candidates -For 3-body decays, the list of available cuts is extended, amoung others, to cuts +For **3-body decays**, the list of available cuts is extended, amoung others, to cuts on the secondary mothers (SM) of combinations of 2 daughters. E.g. for H3L->p + pi + d, cuts on the SM of p-pi can be applied. The cuts that test the mother against the PV, e.g. chi2topo, are inverted for the SM to exclude SMs that come from the PV. -If more than one decay should be reconstructed in the same run, additional +If **more than one decay** should be reconstructed in the same run, additional parfiles at_interface/parfile2_add.txt or at_interface/parfile3_add.txt can be added where 1.) - 3.) can be selected. 4.) - 6.) can not be changed in the same run. The decays can be independent or build a cascade decay. Multiple additional parfiles can be added. -The reconstruction of a cascade decay works the same way as adding an additional decay +The reconstruction of a **cascade decay** works the same way as adding an additional decay to the run as described above. The order of the parfiles need to be from last generation to first generation. Cascade decays with multiple stages as well as combinations of 2- and 3-body-decays can get reconstructed. diff --git a/at_interface/ConverterIn.cpp b/at_interface/ConverterIn.cpp index dde044d..6b12955 100644 --- a/at_interface/ConverterIn.cpp +++ b/at_interface/ConverterIn.cpp @@ -104,7 +104,7 @@ void ConverterIn::Init() { if (pid_mode_ > 1) rec_pdg_field_ = kf_tracks_.GetField("pid"); - + if (pid_mode_ > 2) { prob_p_field_ = kf_tracks_.GetField("prob_p"); prob_pi_field_ = kf_tracks_.GetField("prob_pi"); diff --git a/at_interface/example_parfiles/parfile3_H3L_pidmode4.txt b/at_interface/example_parfiles/parfile3_H3L_pidmode4.txt index e3ec2d0..e53c297 100644 --- a/at_interface/example_parfiles/parfile3_H3L_pidmode4.txt +++ b/at_interface/example_parfiles/parfile3_H3L_pidmode4.txt @@ -67,8 +67,8 @@ H3L name mother 0.7, specific purity tof for background (optional, if = -1: purity for all species will be used) -------------------------------------------------------------------------- ## Input / Output information -pTree treename in input analysistree -RecTracks branchname for reconstructed tracks in input analysistree +aTree treename in input analysistree +RecParticles branchname for reconstructed tracks in input analysistree 0, write detailed BG information -1, number of events to be processed (-1 = all events) 1, make plain tree: =1; make no plain tree: =0 diff --git a/at_interface/main.cpp b/at_interface/main.cpp index 0860796..57a617d 100644 --- a/at_interface/main.cpp +++ b/at_interface/main.cpp @@ -82,8 +82,8 @@ int main(int argc, char** argv) { fscanf(inputInfo, "%*[^\n]%*c"); fscanf(inputInfo, "%f %*[^\n]%*c", &dist); - fscanf(inputInfo, "%f %*[^\n]%*c", &distSV); if (ndaughters == 3) { + fscanf(inputInfo, "%f %*[^\n]%*c", &distSV); fscanf(inputInfo, "%f %*[^\n]%*c", &chi2geoSM.at(0)); fscanf(inputInfo, "%f %*[^\n]%*c", &chi2geoSM.at(1)); fscanf(inputInfo, "%f %*[^\n]%*c", &chi2geoSM.at(2)); @@ -148,7 +148,7 @@ int main(int argc, char** argv) { atree_name = atree_name_c; rec_tracks_name = rec_tracks_name_c; } - + fclose(inputInfo); std::vector daughters; From 4901bd7266f9c708274ab160e0dc2ad42fae3e10 Mon Sep 17 00:00:00 2001 From: SusanneGlaessel Date: Fri, 25 Apr 2025 20:07:25 +0200 Subject: [PATCH 4/4] change in readme --- README.md | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) diff --git a/README.md b/README.md index 2773bc2..8b79af5 100644 --- a/README.md +++ b/README.md @@ -138,19 +138,22 @@ The user can select: - Output format: default output is analysistree format optional: plain tree: output is a simple root tree containing the candidates -For **3-body decays**, the list of available cuts is extended, amoung others, to cuts +### 3-body-decays +For 3-body decays, the list of available cuts is extended, amoung others, to cuts on the secondary mothers (SM) of combinations of 2 daughters. E.g. for H3L->p + pi + d, cuts on the SM of p-pi can be applied. The cuts that test the mother against the PV, e.g. chi2topo, are inverted for the SM to exclude SMs that come from the PV. -If **more than one decay** should be reconstructed in the same run, additional +### Multiple decays +If more than one decay should be reconstructed in the same run, additional parfiles at_interface/parfile2_add.txt or at_interface/parfile3_add.txt can be added where 1.) - 3.) can be selected. 4.) - 6.) can not be changed in the same run. The decays can be independent or build a cascade decay. Multiple additional parfiles can be added. -The reconstruction of a **cascade decay** works the same way as adding an additional decay +### Cascade decays +The reconstruction of a cascade decay works the same way as adding an additional decay to the run as described above. The order of the parfiles need to be from last generation to first generation. Cascade decays with multiple stages as well as combinations of 2- and 3-body-decays can get reconstructed.