@@ -199,71 +199,71 @@ struct StrangeCascTrack {
199199 return true ;
200200 }
201201 // checks TOF PID of dau tracks
202- template <typename TCascade>
203- bool passesTOF (TCascade cascade, TString particle)
204- {
205- return true ;
206- // const auto& bachDaughterTrackCasc = cascade.bachTrackExtra_as<DauTracks>();
207- // const auto& posDaughterTrackCasc = cascade.posTrackExtra_as<DauTracks>();
208- // const auto& negDaughterTrackCasc = cascade.negTrackExtra_as<DauTracks>();
209- // bool xiPassTOFSelection = true;
210- // bool omegaPassTOFSelection = true;
211- // if (cascade.sign() < 0) {
212- // if (posDaughterTrackCasc.hasTOF()) {
213- // if (std::abs(cascade.tofNSigmaXiLaPr()) > selCuts.NSigmaTOFProton) {
214- // xiPassTOFSelection &= false;
215- // }
216- // if (std::abs(cascade.tofNSigmaOmLaPr()) > selCuts.NSigmaTOFProton) {
217- // omegaPassTOFSelection &= false;
218- // }
219- // }
220- // if (negDaughterTrackCasc.hasTOF()) {
221- // if (std::abs(cascade.tofNSigmaXiLaPi()) > selCuts.NSigmaTOFPion) {
222- // xiPassTOFSelection &= false;
223- // }
224- // if (std::abs(cascade.tofNSigmaOmLaPi()) > selCuts.NSigmaTOFPion) {
225- // omegaPassTOFSelection &= false;
226- // }
227- // }
228- // } else {
229- // if (posDaughterTrackCasc.hasTOF()) {
230- // if (std::abs(cascade.tofNSigmaXiLaPi()) > selCuts.NSigmaTOFPion) {
231- // xiPassTOFSelection &= false;
232- // }
233- // if (std::abs(cascade.tofNSigmaOmLaPi()) > selCuts.NSigmaTOFPion) {
234- // omegaPassTOFSelection &= false;
235- // }
236- // }
237- // if (negDaughterTrackCasc.hasTOF()) {
238- // if (std::abs(cascade.tofNSigmaXiLaPr()) > selCuts.NSigmaTOFProton) {
239- // xiPassTOFSelection &= false;
240- // }
241- // if (std::abs(cascade.tofNSigmaOmLaPr()) > selCuts.NSigmaTOFProton) {
242- // omegaPassTOFSelection &= false;
243- // }
244- // }
245- // }
246-
247- // if (bachDaughterTrackCasc.hasTOF()) {
248- // if (std::abs(cascade.tofNSigmaXiPi()) > selCuts.NSigmaTOFPion) {
249- // xiPassTOFSelection &= false;
250- // }
251- // if (std::abs(cascade.tofNSigmaOmKa()) > selCuts.NSigmaTOFKaon) {
252- // omegaPassTOFSelection &= false;
253- // }
254- // }
255-
256- // if (bachDaughterTrackCasc.hasTOF()) {
257- // if (std::abs(cascade.tofNSigmaXiPi()) > selCuts.NSigmaTOFPion) {
258- // xiPassTOFSelection &= false;
259- // }
260- // if (std::abs(cascade.tofNSigmaOmKa()) > selCuts.NSigmaTOFKaon) {
261- // omegaPassTOFSelection &= false;
262- // }
263- // }
264-
265- // if (particle == "xi") {return xiPassTOFSelection;} else {return omegaPassTOFSelection;}
266- }
202+ // template <typename TCascade>
203+ // bool passesTOF(TCascade cascade, TString particle)
204+ // {
205+ // return true;
206+ // // const auto& bachDaughterTrackCasc = cascade.bachTrackExtra_as<DauTracks>();
207+ // // const auto& posDaughterTrackCasc = cascade.posTrackExtra_as<DauTracks>();
208+ // // const auto& negDaughterTrackCasc = cascade.negTrackExtra_as<DauTracks>();
209+ // // bool xiPassTOFSelection = true;
210+ // // bool omegaPassTOFSelection = true;
211+ // // if (cascade.sign() < 0) {
212+ // // if (posDaughterTrackCasc.hasTOF()) {
213+ // // if (std::abs(cascade.tofNSigmaXiLaPr()) > selCuts.NSigmaTOFProton) {
214+ // // xiPassTOFSelection &= false;
215+ // // }
216+ // // if (std::abs(cascade.tofNSigmaOmLaPr()) > selCuts.NSigmaTOFProton) {
217+ // // omegaPassTOFSelection &= false;
218+ // // }
219+ // // }
220+ // // if (negDaughterTrackCasc.hasTOF()) {
221+ // // if (std::abs(cascade.tofNSigmaXiLaPi()) > selCuts.NSigmaTOFPion) {
222+ // // xiPassTOFSelection &= false;
223+ // // }
224+ // // if (std::abs(cascade.tofNSigmaOmLaPi()) > selCuts.NSigmaTOFPion) {
225+ // // omegaPassTOFSelection &= false;
226+ // // }
227+ // // }
228+ // // } else {
229+ // // if (posDaughterTrackCasc.hasTOF()) {
230+ // // if (std::abs(cascade.tofNSigmaXiLaPi()) > selCuts.NSigmaTOFPion) {
231+ // // xiPassTOFSelection &= false;
232+ // // }
233+ // // if (std::abs(cascade.tofNSigmaOmLaPi()) > selCuts.NSigmaTOFPion) {
234+ // // omegaPassTOFSelection &= false;
235+ // // }
236+ // // }
237+ // // if (negDaughterTrackCasc.hasTOF()) {
238+ // // if (std::abs(cascade.tofNSigmaXiLaPr()) > selCuts.NSigmaTOFProton) {
239+ // // xiPassTOFSelection &= false;
240+ // // }
241+ // // if (std::abs(cascade.tofNSigmaOmLaPr()) > selCuts.NSigmaTOFProton) {
242+ // // omegaPassTOFSelection &= false;
243+ // // }
244+ // // }
245+ // // }
246+ //
247+ // // if (bachDaughterTrackCasc.hasTOF()) {
248+ // // if (std::abs(cascade.tofNSigmaXiPi()) > selCuts.NSigmaTOFPion) {
249+ // // xiPassTOFSelection &= false;
250+ // // }
251+ // // if (std::abs(cascade.tofNSigmaOmKa()) > selCuts.NSigmaTOFKaon) {
252+ // // omegaPassTOFSelection &= false;
253+ // // }
254+ // // }
255+ //
256+ // // if (bachDaughterTrackCasc.hasTOF()) {
257+ // // if (std::abs(cascade.tofNSigmaXiPi()) > selCuts.NSigmaTOFPion) {
258+ // // xiPassTOFSelection &= false;
259+ // // }
260+ // // if (std::abs(cascade.tofNSigmaOmKa()) > selCuts.NSigmaTOFKaon) {
261+ // // omegaPassTOFSelection &= false;
262+ // // }
263+ // // }
264+ //
265+ // // if (particle == "xi") {return xiPassTOFSelection;} else {return omegaPassTOFSelection;}
266+ // }
267267 // checks whether gen cascade corresponds to PDG code
268268 template <typename TCascade>
269269 bool isValidPDG (TCascade cascade, TString particle)
@@ -403,18 +403,18 @@ struct StrangeCascTrack {
403403 // apply tof pid
404404 bool passedAllSelsXi = passedAllSels;
405405 bool passedAllSelsOmega = passedAllSels;
406- if (doApplyTOFPID) {
407- if (passesTOF (cascade, " xi" )) {
408- histos.fill (HIST (kTypeNames [type]) + HIST (" /MassXiTOFPID" ), massXi);
409- } else {
410- passedAllSelsXi = false ;
411- }
412- if (passesTOF (cascade, " omega" )) {
413- histos.fill (HIST (kTypeNames [type]) + HIST (" /MassOmegaTOFPID" ), massOmega);
414- } else {
415- passedAllSelsOmega = false ;
416- }
417- }
406+ // if (doApplyTOFPID) {
407+ // if (passesTOF(cascade, "xi")) {
408+ // histos.fill(HIST(kTypeNames[type]) + HIST("/MassXiTOFPID"), massXi);
409+ // } else {
410+ // passedAllSelsXi = false;
411+ // }
412+ // if (passesTOF(cascade, "omega")) {
413+ // histos.fill(HIST(kTypeNames[type]) + HIST("/MassOmegaTOFPID"), massOmega);
414+ // } else {
415+ // passedAllSelsOmega = false;
416+ // }
417+ // }
418418 // apply competing mass rej
419419 if (doCompetingMassRej) {
420420 if (std::abs (massXi - pdgDB->Mass (3312 )) > selCuts.compMassRej ) {
0 commit comments