@@ -9,7 +9,7 @@ int External() {
99 std ::array < float , 6 > freqRepl = {0.1 , 0.1 , 0.1 , 0.1 , 0.5 , 0.5 };
1010 std ::map < int , int > sumOrigReplacedParticles = {{10433 , 0 }, {435 , 0 }, {425 , 0 }};
1111
12- std ::array < int , 11 > checkPdgHadron {411 , 421 , 10433 , 30433 , 435 , 437 , 4325 , 4326 , 4315 , 4316 , 531 };
12+ std ::array < int , 12 > checkPdgHadron {411 , 421 , 10433 , 30433 , 435 , 437 , 4325 , 4326 , 4315 , 4316 , 531 , 425 };
1313 std ::map < int , std ::vector < std ::vector < int >>> checkHadronDecays { // sorted pdg of daughters
1414 {411 , {{-321 , 211 , 211 }, {-313 , 211 }, {211 , 311 }, {211 , 333 }}}, // D+
1515 {421 , {{-321 , 211 }, {-321 , 211 , 111 }}}, // D0
@@ -21,7 +21,8 @@ int External() {
2121 {4326 , {{411 , 3122 }}}, // Xic(3080)+
2222 {4315 , {{421 , 3122 }}}, // Xic(3055)+
2323 {4316 , {{421 , 3122 }}}, // Xic(3080)+
24- {531 , {{-435 , -11 , 12 }, {-10433 , -11 , 12 }, {-435 , -13 , 14 }, {-10433 , -13 , 14 }, {-435 , -15 , 16 }, {-10433 , -15 , 16 }, {-435 , 211 }}}// Bs0
24+ {531 , {{-435 , -11 , 12 }, {-10433 , -11 , 12 }, {-435 , -13 , 14 }, {-10433 , -13 , 14 }, {-435 , -15 , 16 }, {-10433 , -15 , 16 }, {-435 , 211 }}}, // Bs0
25+ {425 , {{413 , -211 }, {423 , 111 }, {411 , -211 }, {421 , 111 }, {413 , -211 , 111 }, {423 , 211 , -211 }}}
2526 };
2627
2728 TFile file (path .c_str (), "READ ");
@@ -81,7 +82,7 @@ int External() {
8182 for (int j {track .getFirstDaughterTrackId ()}; j <= track .getLastDaughterTrackId (); ++ j ) {
8283 auto pdgDau = tracks -> at (j ).GetPdgCode ();
8384 pdgsDecay .push_back (pdgDau );
84- if (pdgDau != 333 ) { // phi is antiparticle of itself
85+ if (pdgDau != 333 && pdgDau != 111 ) { // phi and pi0 are antiparticles of themselves
8586 pdgsDecayAntiPart .push_back (- pdgDau );
8687 } else {
8788 pdgsDecayAntiPart .push_back (pdgDau );
0 commit comments