@@ -35,16 +35,16 @@ struct HfBsToDsPiCandidateSelector {
3535 Configurable<double > pTCandMin{" pTCandMin" , 0 ., " Lower bound of candidate pT" };
3636 Configurable<double > pTCandMax{" pTCandMax" , 50 ., " Upper bound of candidate pT" };
3737
38- // Track quality
38+ // Track quality
3939 Configurable<double > TPCNClsFindablePIDCut{" TPCNClsFindablePIDCut" , 70 ., " Lower bound of TPC findable clusters for good PID" };
4040
41- // TPC PID
41+ // TPC PID
4242 Configurable<double > pidTPCMinpT{" pidTPCMinpT" , 0.15 , " Lower bound of track pT for TPC PID" };
4343 Configurable<double > pidTPCMaxpT{" pidTPCMaxpT" , 10 ., " Upper bound of track pT for TPC PID" };
4444 Configurable<double > nSigmaTPC{" nSigmaTPC" , 5 ., " Nsigma cut on TPC only" };
4545 Configurable<double > nSigmaTPCCombined{" nSigmaTPCCombined" , 5 ., " Nsigma cut on TPC combined with TOF" };
4646
47- // TOF PID
47+ // TOF PID
4848 Configurable<double > pidTOFMinpT{" pidTOFMinpT" , 0.15 , " Lower bound of track pT for TOF PID" };
4949 Configurable<double > pidTOFMaxpT{" pidTOFMaxpT" , 10 ., " Upper bound of track pT for TOF PID" };
5050 Configurable<double > nSigmaTOF{" nSigmaTOF" , 5 ., " Nsigma cut on TOF only" };
@@ -65,56 +65,56 @@ struct HfBsToDsPiCandidateSelector {
6565 return false ;
6666 }
6767
68- // Bs mass cut
68+ // Bs mass cut
6969 if (std::abs (InvMassBsToDsPi (hfCandBs) - RecoDecay::getMassPDG (pdg::Code::kBs )) > cuts->get (pTBin, " m" )) {
70- // Printf("Bs topol selection failed at mass diff check");
70+ // Printf("Bs topol selection failed at mass diff check");
7171 return false ;
7272 }
7373
74- // pion pt
74+ // pion pt
7575 if (trackPi.pt () < cuts->get (pTBin, " pT Pi" )) {
7676 return false ;
7777 }
7878
79- // Ds+ pt
79+ // Ds+ pt
8080 if (hfCandDs.pt () < cuts->get (pTBin, " pT Ds+" )) {
8181 return false ;
8282 }
8383
84- // Ds mass
85- // if (trackPi.sign() < 0) {
86- // if (std::abs(InvMassDspKpi(hfCandDs) - RecoDecay::getMassPDG(pdg::Code::kDs)) > cuts->get(pTBin, "DeltaMDs")) {
87- // return false;
88- // }
89- // }
84+ // Ds mass
85+ // if (trackPi.sign() < 0) {
86+ // if (std::abs(InvMassDspKpi(hfCandDs) - RecoDecay::getMassPDG(pdg::Code::kDs)) > cuts->get(pTBin, "DeltaMDs")) {
87+ // return false;
88+ // }
89+ // }
9090
91- // Bs Decay length
91+ // Bs Decay length
9292 if (hfCandBs.decayLength () < cuts->get (pTBin, " Bs decLen" )) {
9393 return false ;
9494 }
9595
96- // Bs Decay length XY
96+ // Bs Decay length XY
9797 if (hfCandBs.decayLengthXY () < cuts->get (pTBin, " Bs decLenXY" )) {
9898 return false ;
9999 }
100100
101- // Bs chi2PCA cut
101+ // Bs chi2PCA cut
102102 if (hfCandBs.chi2PCA () > cuts->get (pTBin, " Chi2PCA" )) {
103- // Printf("Bs selection failed at chi2PCA");
103+ // Printf("Bs selection failed at chi2PCA");
104104 return false ;
105105 }
106106
107- // Bs CPA cut
107+ // Bs CPA cut
108108 if (hfCandBs.cpa () < cuts->get (pTBin, " CPA" )) {
109109 return false ;
110110 }
111111
112- // d0 of pi
112+ // d0 of pi
113113 if (std::abs (hfCandBs.impactParameter1 ()) < cuts->get (pTBin, " d0 Pi" )) {
114114 return false ;
115115 }
116116
117- // d0 of Ds+
117+ // d0 of Ds+
118118 if (std::abs (hfCandBs.impactParameter0 ()) < cuts->get (pTBin, " d0 Ds+" )) {
119119 return false ;
120120 }
@@ -124,30 +124,30 @@ struct HfBsToDsPiCandidateSelector {
124124
125125 void process (aod::HfCandBs const & hfCandBss, soa::Join<aod::HfCandProng3, aod::HFSelDsCandidate>, aod::BigTracksPID const &)
126126 {
127- for (auto & hfCandBs : hfCandBss) { // looping over Bs candidates
127+ for (auto & hfCandBs : hfCandBss) { // looping over Bs candidates
128128
129129 int statusBs = 0 ;
130130
131131 // check if flagged as Bs --> Ds+ π-
132132 if (!(hfCandBs.hfflag () & 1 << hf_cand_bs::DecayType::BsToDsPi)) {
133133 hfSelBsToDsPiCandidate (statusBs);
134- // Printf("Bs candidate selection failed at hfflag check");
134+ // Printf("Bs candidate selection failed at hfflag check");
135135 continue ;
136136 }
137137
138138 // Ds is always index0 and pi is index1 by default
139139 auto candDs = hfCandBs.index0_as <soa::Join<aod::HfCandProng3, aod::HFSelDsCandidate>>();
140140 auto trackPi = hfCandBs.index1_as <aod::BigTracksPID>();
141141
142- // topological cuts
142+ // topological cuts
143143 if (!selectionTopol (hfCandBs, candDs, trackPi)) {
144144 hfSelBsToDsPiCandidate (statusBs);
145145 // Printf("Bs candidate selection failed at selection topology");
146146 continue ;
147147 }
148148
149149 hfSelBsToDsPiCandidate (1 );
150- // Printf("Bs candidate selection successful, candidate should be selected");
150+ // Printf("Bs candidate selection successful, candidate should be selected");
151151 }
152152 }
153153};
0 commit comments