@@ -36,10 +36,11 @@ struct sgPIDAnalyzer {
3636
3737 ConfigurableAxis ptAxis{
3838 " ptAxis" ,
39- {198 , 0.1 , 10.0 },
39+ {200 , 0.0 , 10.0 },
4040 " Pt binning" };
4141
42- ConfigurableAxis sigmaAxis{" sigmaAxis" , {100 , -50 , 50 }, " nSigma TPC binning" };
42+ ConfigurableAxis sigmaAxis{" sigmaAxis" , {1000 , -20 , 180 }, " nSigma TPC binning" };
43+ ConfigurableAxis tofAxis{" tofAxis" , {200 , -10 , 10 }, " nSigma TOF binning" };
4344 Configurable<float > eta_min{" eta_min" , -0.9 , " Track Pseudorapidity" };
4445 Configurable<float > eta_max{" eta_max" , 0.9 , " Track Pseudorapidity" };
4546
@@ -48,6 +49,8 @@ struct sgPIDAnalyzer {
4849
4950 const AxisSpec ptBins{ptAxis, " p_{T} axis" };
5051 const AxisSpec nSigmaBins{sigmaAxis, " pseudo rapidity axis" };
52+ const AxisSpec ntofBins{tofAxis, " pseudo rapidity axis" };
53+ histos.add (" Events" , " Selected Events" , {HistType::kTH1F , {{3 , -.5 , 2.5 }}});
5154 histos.add (" TPC/pTPC_Pi" , " Positive TPC Pi Tracks" , {HistType::kTH2F , {ptBins, nSigmaBins}});
5255 histos.add (" TPC/nTPC_Pi" , " Negative TPC Pi Tracks" , {HistType::kTH2F , {ptBins, nSigmaBins}});
5356 histos.add (" TPC/pTPC_Ka" , " Positive TPC Ka Tracks" , {HistType::kTH2F , {ptBins, nSigmaBins}});
@@ -93,28 +96,29 @@ struct sgPIDAnalyzer {
9396 histos.add (" TPC/nTPC_El_Ka" , " Positive TPC Pr vs Ka" , {HistType::kTH2F , {ptBins, nSigmaBins}});
9497 histos.add (" TPC/nTPC_El_Pr" , " Positive TPC Pr vs El" , {HistType::kTH2F , {ptBins, nSigmaBins}});
9598
96- histos.add (" TOF/pPi" , " Positive TPC Pi vs TOF Pi vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
97- histos.add (" TOF/nPi" , " Negative TPC Pi vs TOF Pi vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
98- histos.add (" TOF/pKa" , " Positive TPC Ka vs TOF Ka vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
99- histos.add (" TOF/nKa" , " Negative TPC Ka vs TOF Ka vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
100- histos.add (" TOF/pPr" , " Positive TPC Pr vs TOF Pr vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
101- histos.add (" TOF/nPr" , " Negative TPC Pr vs TOF Pr vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
102- histos.add (" TOF/pEl" , " Positive TPC El vs TOF El vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
103- histos.add (" TOF/nEl" , " Negative TPC El vs TOF El vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
104- histos.add (" TOF/pDe" , " Positive TPC De vs TOF Pi vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
105- histos.add (" TOF/nDe" , " Negative TPC De vs TOF Pi vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
106- histos.add (" TOF/pTr" , " Positive TPC Tr vs TOF Ka vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
107- histos.add (" TOF/nTr" , " Negative TPC Tr vs TOF Ka vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
108- histos.add (" TOF/pHe" , " Positive TPC He vs TOF Pr vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
109- histos.add (" TOF/nHe" , " Negative TPC He vs TOF Pr vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
110- histos.add (" TOF/pAl" , " Positive TPC Al vs TOF El vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
111- histos.add (" TOF/nAl" , " Negative TPC Al vs TOF El vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
112- histos.add (" TOF/pMu" , " Positive TPC Mu vs TOF El vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
113- histos.add (" TOF/nMu" , " Negative TPC Mu vs TOF El vs pt" , {HistType::kTH3F , {ptBins, nSigmaBins, nSigmaBins }});
99+ histos.add (" TOF/pPi" , " Positive TPC Pi vs TOF Pi vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
100+ histos.add (" TOF/nPi" , " Negative TPC Pi vs TOF Pi vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
101+ histos.add (" TOF/pKa" , " Positive TPC Ka vs TOF Ka vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
102+ histos.add (" TOF/nKa" , " Negative TPC Ka vs TOF Ka vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
103+ histos.add (" TOF/pPr" , " Positive TPC Pr vs TOF Pr vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
104+ histos.add (" TOF/nPr" , " Negative TPC Pr vs TOF Pr vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
105+ histos.add (" TOF/pEl" , " Positive TPC El vs TOF El vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
106+ histos.add (" TOF/nEl" , " Negative TPC El vs TOF El vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
107+ histos.add (" TOF/pDe" , " Positive TPC De vs TOF Pi vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
108+ histos.add (" TOF/nDe" , " Negative TPC De vs TOF Pi vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
109+ histos.add (" TOF/pTr" , " Positive TPC Tr vs TOF Ka vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
110+ histos.add (" TOF/nTr" , " Negative TPC Tr vs TOF Ka vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
111+ histos.add (" TOF/pHe" , " Positive TPC He vs TOF Pr vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
112+ histos.add (" TOF/nHe" , " Negative TPC He vs TOF Pr vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
113+ histos.add (" TOF/pAl" , " Positive TPC Al vs TOF El vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
114+ histos.add (" TOF/nAl" , " Negative TPC Al vs TOF El vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
115+ histos.add (" TOF/pMu" , " Positive TPC Mu vs TOF El vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
116+ histos.add (" TOF/nMu" , " Negative TPC Mu vs TOF El vs pt" , {HistType::kTH3F , {ptBins, ntofBins, ntofBins }});
114117 }
115-
116- void process (aod::SGEvents const &, aod::SGTracks const & tracks)
118+ using SGEvent = aod::SGEvents::iterator;
119+ void process (SGEvent const & event , aod::SGTracks const & tracks)
117120 {
121+ histos.fill (HIST (" Events" ), event.gs ());
118122 for (const auto & track : tracks) {
119123 if (track.eta () < eta_min || track.eta () > eta_max)
120124 continue ;
0 commit comments