@@ -64,7 +64,6 @@ namespace
6464constexpr int kNpart = 2 ;
6565constexpr float kTrackSels [12 ]{/* 60, */ 80 , 100 , 2 , 3 , /* 4, */ 0.05 , 0.1 , /* 0.15, */ 0.5 , 1 , /* 1.5, */ 2 , 3 /* , 4 */ , 2 , 3 , /* , 4 */ };
6666constexpr float kDcaSelsParam [3 ][3 ]{{-1 .e32 , -1 .e32 , -1 .e32 }, {-1 .e32 , -1 .e32 , -1 .e32 }, {-1 .e32 , -1 .e32 , -1 .e32 }};
67- constexpr float kDcaSels [3 ]{10 ., 10 ., 10 .};
6867constexpr double kBetheBlochDefault [kNpart ][6 ]{{-1 .e32 , -1 .e32 , -1 .e32 , -1 .e32 , -1 .e32 , -1 .e32 }, {-1 .e32 , -1 .e32 , -1 .e32 , -1 .e32 , -1 .e32 , -1 .e32 }};
6968constexpr double kBetheBlochDefaultITS [6 ]{-1 .e32 , -1 .e32 , -1 .e32 , -1 .e32 , -1 .e32 , -1 .e32 };
7069constexpr double kEstimatorsCorrelationCoef [2 ]{-0.669108 , 1.04489 };
@@ -258,7 +257,6 @@ struct EbyeMaker {
258257 Configurable<float > lambdaMassCut{" lambdaMassCut" , 0 .02f , " maximum deviation from PDG mass (for QA histograms)" };
259258
260259 Configurable<LabeledArray<float >> cfgTrackSels{" cfgTrackSels" , {kTrackSels , 1 , 12 , particleName, trackSelsNames}, " Track selections" };
261- Configurable<LabeledArray<float >> cfgDcaSels{" cfgDcaSels" , {kDcaSels , 1 , 3 , particleName, dcaSelsNames}, " DCA selections" };
262260 Configurable<LabeledArray<float >> cfgDcaSelsParam{" cfgDcaSelsParam" , {kDcaSelsParam [0 ], 3 , 3 , dcaSelsNames, dcaParNames}, " DCA threshold settings" };
263261
264262 std::array<float , kNpart > ptMin;
@@ -509,8 +507,8 @@ struct EbyeMaker {
509507
510508 // tracking variables QA
511509 histos.add <TH2>(" QA/tpcCRvsCls" , " ;#it{N}_{TPCCR};#it{N}_{TPCcls}" , HistType::kTH2F , {nTpcAxis, nTpcAxis});
512- histos.add <TH1 >(" QA/dcaxyVsPt" , " ;#it{p}_{T} (GeV/#it{c});DCA_{#it{xy}} (cm)" , HistType::kTH2F , {momAxis, dcaAxis});
513- histos.add <TH1 >(" QA/dcazVsPt" , " ;#it{p}_{T} (GeV/#it{c});DCA_{#it{z}} (cm)" , HistType::kTH2F , {momAxis, dcaAxis});
510+ histos.add <TH2 >(" QA/dcaxyVsPt" , " ;#it{p}_{T} (GeV/#it{c});DCA_{#it{xy}} (cm)" , HistType::kTH2F , {momAxis, dcaAxis});
511+ histos.add <TH2 >(" QA/dcazVsPt" , " ;#it{p}_{T} (GeV/#it{c});DCA_{#it{z}} (cm)" , HistType::kTH2F , {momAxis, dcaAxis});
514512
515513 ptMin = std::array<float , kNpart >{antipPtMin, antidPtMin};
516514 ptMax = std::array<float , kNpart >{antipPtMax, antidPtMax};
@@ -563,7 +561,7 @@ struct EbyeMaker {
563561 auto trackEta = trackParCov.getEta ();
564562 histos.fill (HIST (" QA/dcaxyVsPt" ), track.pt (), dcaInfo[0 ]);
565563 histos.fill (HIST (" QA/dcazVsPt" ), track.pt (), dcaInfo[1 ]);
566- if (std::abs (dcaInfo[0 ]) > cfgDcaSels-> get ( " dcaxy " ) * dcaSigma (track.pt (), " dcaxy" ) || std::abs (dcaInfo[1 ]) > cfgDcaSels-> get ( " dcaz " ) * dcaSigma (track.pt (), " dcaz" ) || dca > cfgDcaSels-> get ( " dca " ) * dcaSigma (track.pt (), " dca" )) { // dcaxy
564+ if (std::abs (dcaInfo[0 ]) > dcaSigma (track.pt (), " dcaxy" ) || std::abs (dcaInfo[1 ]) > dcaSigma (track.pt (), " dcaz" ) || dca > dcaSigma (track.pt (), " dca" )) { // dcaxy
567565 continue ;
568566 }
569567 histos.fill (HIST (" QA/tpcSignal" ), track.tpcInnerParam (), track.tpcSignal ());
0 commit comments