1515// / \author Alberto Caliva (alberto.caliva@cern.ch), Francesca Ercolessi (francesca.ercolessi@cern.ch), Nicolò Jacazio (nicolo.jacazio@cern.ch)
1616// / \since May 22, 2024
1717
18- #include < TLorentzVector.h>
1918#include < TMath.h>
2019#include < TObjArray.h>
2120#include < TPDGCode.h>
@@ -1681,6 +1680,7 @@ struct StrangenessInJets {
16811680 }
16821681 }
16831682
1683+ /*
16841684 // Reconstructed Tracks
16851685 for (const auto& track : tracksPerColl) {
16861686
@@ -1788,22 +1788,25 @@ struct StrangenessInJets {
17881788 registryMC.fill(HIST("pi_minus_rec_in_ue_tof"), multiplicity, track.pt(), wPiminusUe);
17891789 }
17901790 }
1791-
1791+ */
17921792 for (const auto & mcParticle : mcParticlesPerColl) {
17931793
17941794 if (mcParticle.eta () < etaMin || mcParticle.eta () > etaMax)
17951795 continue ;
17961796 if (!mcParticle.isPhysicalPrimary ())
17971797 continue ;
17981798
1799+ /*
17991800 double wPiplusJet(1.0), wPiplusUe(1.0);
18001801 double wPiminusJet(1.0), wPiminusUe(1.0);
18011802 double wKaplusJet(1.0), wKaplusUe(1.0);
18021803 double wKaminusJet(1.0), wKaminusUe(1.0);
18031804 double wPrplusJet(1.0), wPrplusUe(1.0);
18041805 double wPrminusJet(1.0), wPrminusUe(1.0);
1806+ */
18051807 double wK0jet (1.0 ), wK0Ue (1.0 ), wLambdaJet (1.0 ), wLambdaUe (1.0 ), wAntilambdaJet (1.0 ), wAntilambdaUe (1.0 );
18061808 if (applyReweighting) {
1809+ /*
18071810 auto getWeight = [&](TH2F* histo) {
18081811 if (!histo) {
18091812 return 1.0;
@@ -1822,7 +1825,7 @@ struct StrangenessInJets {
18221825 wPiplusUe = getWeight(twodWeightsPiplusUe);
18231826 wPiminusJet = getWeight(twodWeightsPiminusJet);
18241827 wPiminusUe = getWeight(twodWeightsPiminusUe);
1825-
1828+ */
18261829 int ix = twodWeightsK0Jet->GetXaxis ()->FindBin (mcParticle.pt ());
18271830 int iy = twodWeightsK0Jet->GetYaxis ()->FindBin (mcParticle.eta ());
18281831 wK0jet = twodWeightsK0Jet->GetBinContent (ix, iy);
@@ -1852,6 +1855,7 @@ struct StrangenessInJets {
18521855 }
18531856
18541857 switch (mcParticle.pdgCode ()) {
1858+ /*
18551859 case kPiPlus: // Pi+
18561860 registryMC.fill(HIST("mc_pi_plus/in_jet/gen"), multiplicity, mcParticle.pt(), wPiplusJet);
18571861 registryMC.fill(HIST("mc_pi_plus/in_ue/gen"), multiplicity, mcParticle.pt(), wPiplusUe);
@@ -1882,6 +1886,7 @@ struct StrangenessInJets {
18821886 registryMC.fill(HIST("mc_pr_minus/in_ue/gen"), multiplicity, mcParticle.pt(), wPrminusUe);
18831887 registryMC.fill(HIST("pr_minus_eta_pt_pythia"), mcParticle.pt(), mcParticle.eta());
18841888 break;
1889+ */
18851890 case kK0Short : // K0s
18861891 registryMC.fill (HIST (" K0s_generated_jet" ), multiplicity, mcParticle.pt (), wK0jet);
18871892 registryMC.fill (HIST (" K0s_generated_ue" ), multiplicity, mcParticle.pt (), wK0Ue);
@@ -2015,12 +2020,14 @@ struct StrangenessInJets {
20152020 // In jet
20162021 if (deltaRjet < rJet) {
20172022 switch (particle.pdgCode ()) {
2023+ /*
20182024 case kPiPlus:
20192025 registryMC.fill(HIST("pi_plus_eta_pt_jet"), particle.pt(), particle.eta());
20202026 break;
20212027 case kPiMinus:
20222028 registryMC.fill(HIST("pi_minus_eta_pt_jet"), particle.pt(), particle.eta());
20232029 break;
2030+ */
20242031 case kK0Short :
20252032 registryMC.fill (HIST (" K0s_eta_pt_jet" ), particle.pt (), particle.eta ());
20262033 break ;
@@ -2049,12 +2056,14 @@ struct StrangenessInJets {
20492056
20502057 if (deltaRue1 < rJet || deltaRue2 < rJet) {
20512058 switch (particle.pdgCode ()) {
2059+ /*
20522060 case kPiPlus:
20532061 registryMC.fill(HIST("pi_plus_eta_pt_ue"), particle.pt(), particle.eta());
20542062 break;
20552063 case kPiMinus:
20562064 registryMC.fill(HIST("pi_minus_eta_pt_ue"), particle.pt(), particle.eta());
20572065 break;
2066+ */
20582067 case kK0Short :
20592068 registryMC.fill (HIST (" K0s_eta_pt_ue" ), particle.pt (), particle.eta ());
20602069 break ;
0 commit comments