@@ -102,7 +102,6 @@ struct PiDecandidate {
102102 std::array<float , 3 > momDe = {99 .f , 99 .f , 99 .f };
103103 std::array<float , 3 > momPi = {99 .f , 99 .f , 99 .f };
104104
105-
106105 float signDe = 1 .f;
107106 float signPi = 1 .f;
108107 float invMass = -10 .f;
@@ -289,7 +288,7 @@ struct PiDeuteronFemto {
289288 for (int i = 0 ; i < numParticles; i++) {
290289 mBBparamsDe [i] = settingBetheBlochParams->get (" De" , Form (" p%i" , i));
291290 }
292- mBBparamsDe [5 ] = settingBetheBlochParams->get (" De" , " resolution" );
291+ mBBparamsDe [5 ] = settingBetheBlochParams->get (" De" , " resolution" );
293292
294293 std::vector<std::string> selectionLabels = {" All" , " Track selection" , " PID" };
295294 for (int i = 0 ; i < Selections::kAll ; i++) {
@@ -404,7 +403,7 @@ struct PiDeuteronFemto {
404403 auto tpcNSigmaPi = candidate.tpcNSigmaPi ();
405404 mQaRegistry .fill (HIST (" h2NsigmaPiTPC_preselection" ), candidate.tpcInnerParam (), tpcNSigmaPi);
406405 if (std::abs (candidate.pt ()) < settingCutPiptMin || std::abs (candidate.pt ()) > settingCutPiptMax)
407- return false ;
406+ return false ;
408407 if (candidate.hasTOF () && candidate.tpcInnerParam () > settingCutPinMinTOFPi) {
409408 auto tofNSigmaPi = candidate.tofNSigmaPi ();
410409 auto combNsigma = std::sqrt (tofNSigmaPi * tofNSigmaPi + tpcNSigmaPi * tpcNSigmaPi);
@@ -416,12 +415,12 @@ struct PiDeuteronFemto {
416415 if (combNsigma > settingCutNsigmaTOFTPCPi) {
417416 return false ;
418417 }
419- mQaRegistry .fill (HIST (" h2NsigmaPiTPC" ), candidate.sign () *candidate.pt (), tpcNSigmaPi);
420- mQaRegistry .fill (HIST (" h2NsigmaPiTOF" ), candidate.sign () *candidate.pt (), tofNSigmaPi);
421- mQaRegistry .fill (HIST (" h2NsigmaPiComb" ), candidate.sign () *candidate.pt (), combNsigma);
418+ mQaRegistry .fill (HIST (" h2NsigmaPiTPC" ), candidate.sign () * candidate.pt (), tpcNSigmaPi);
419+ mQaRegistry .fill (HIST (" h2NsigmaPiTOF" ), candidate.sign () * candidate.pt (), tofNSigmaPi);
420+ mQaRegistry .fill (HIST (" h2NsigmaPiComb" ), candidate.sign () * candidate.pt (), combNsigma);
422421 return true ;
423422 } else if (std::abs (tpcNSigmaPi) < settingCutNsigmaTPCPi) {
424- mQaRegistry .fill (HIST (" h2NsigmaPiTPC" ), candidate.sign () *candidate.pt (), tpcNSigmaPi);
423+ mQaRegistry .fill (HIST (" h2NsigmaPiTPC" ), candidate.sign () * candidate.pt (), tpcNSigmaPi);
425424 return true ;
426425 }
427426 return false ;
@@ -445,16 +444,16 @@ struct PiDeuteronFemto {
445444 return false ;
446445 }
447446 float tpcNSigmaDe;
448- if (settingUseBBcomputeDeNsigma) {
449- tpcNSigmaDe = computeNSigmaDe (candidate);
450- }else {
451- tpcNSigmaDe = candidate.tpcNSigmaDe ();
447+ if (settingUseBBcomputeDeNsigma) {
448+ tpcNSigmaDe = computeNSigmaDe (candidate);
449+ } else {
450+ tpcNSigmaDe = candidate.tpcNSigmaDe ();
452451 }
453-
452+
454453 mQaRegistry .fill (HIST (" h2NsigmaDeTPC_preselection" ), candidate.sign () * candidate.pt (), tpcNSigmaDe);
455454 mQaRegistry .fill (HIST (" h2NsigmaDeTPC_preselecComp" ), candidate.sign () * candidate.pt (), candidate.tpcNSigmaDe ());
456455 if (std::abs (candidate.pt ()) < settingCutDeptMin || std::abs (candidate.pt ()) > settingCutDeptMax)
457- return false ;
456+ return false ;
458457 if (candidate.hasTOF () && candidate.tpcInnerParam () > settingCutPinMinTOFITSDe) {
459458 auto tofNSigmaDe = candidate.tofNSigmaDe ();
460459 auto combNsigma = std::sqrt (tofNSigmaDe * tofNSigmaDe + tpcNSigmaDe * tpcNSigmaDe);
@@ -478,7 +477,7 @@ struct PiDeuteronFemto {
478477 }
479478 mQaRegistry .fill (HIST (" h2NsigmaDeTPC" ), candidate.sign () * candidate.pt (), tpcNSigmaDe);
480479 mQaRegistry .fill (HIST (" h2NSigmaDeITS" ), candidate.sign () * candidate.pt (), itsnSigmaDe);
481- // mQaRegistry.fill(HIST("h2NsigmaDeComb"), candidate.sign() * candidate.pt(), combNsigma);
480+ // mQaRegistry.fill(HIST("h2NsigmaDeComb"), candidate.sign() * candidate.pt(), combNsigma);
482481 mQaRegistry .fill (HIST (" h2dEdxDecandidates" ), candidate.sign () * tpcInnerParam, candidate.tpcSignal ());
483482 return true ;
484483 }
0 commit comments