11int External ()
22{
33 std ::string path {"o2sim_Kine.root" };
4- std ::vector < int > checkPdgHadron {521 }; // Bmeson
5- std ::vector < int > nucleiDauPdg {1000010020 }; // d
4+ std ::vector < int > checkPdgHadron {521 };
5+ std ::vector < int > nucleiDauPdg {1000010020 };
66
77 TFile file (path .c_str (), "READ ");
88 if (file .IsZombie ())
@@ -24,10 +24,9 @@ int External()
2424 for (auto & track : * tracks )
2525 {
2626 auto pdg = track .GetPdgCode ();
27- if (std ::find (checkPdgHadron .begin (), checkPdgHadron .end (), std ::abs (pdg )) != checkPdgHadron .end ()) // found signal
27+ if (std ::find (checkPdgHadron .begin (), checkPdgHadron .end (), std ::abs (pdg )) != checkPdgHadron .end ())
2828 {
29- // count signal PDG
30- if (std ::abs (track .GetRapidity ()) > 1.5 ) continue ; // skip if outside rapidity window
29+ if (std ::abs (track .GetRapidity ()) > 1.5 ) continue ;
3130 nSignals ++ ;
3231 for (int j {track .getFirstDaughterTrackId ()}; j <= track .getLastDaughterTrackId (); ++ j )
3332 {
@@ -42,16 +41,8 @@ int External()
4241 }
4342 std ::cout << "--------------------------------\n" ;
4443 std ::cout << "# Events: " << nEvents << "\n" ;
45- std ::cout <<"# signal hadrons: " << nSignals << "\n" ;
46- std ::cout <<"# signal hadrons decaying into nuclei: " << nSignalGoodDecay << "\n" ;
47-
48- float fracForcedDecays = nSignals ? float (nSignalGoodDecay ) / nSignals : 0.0f ;
49- float uncFracForcedDecays = nSignals ? std ::sqrt (fracForcedDecays * (1 - fracForcedDecays ) / nSignals ) : 1.0f ;
50- if (1 - fracForcedDecays > 0.2 + uncFracForcedDecays ) // we put some tolerance (lambdaB in MB events do not coalesce)
51- {
52- std ::cerr << "Fraction of signals decaying into nuclei: " << fracForcedDecays << ", lower than expected\n" ;
53- return 1 ;
54- }
44+ std ::cout << "# signal hadrons: " << nSignals << "\n" ;
45+ std ::cout << "# signal hadrons decaying into nuclei: " << nSignalGoodDecay << "\n" ;
5546
5647 return 0 ;
57- }
48+ }
0 commit comments