@@ -36,21 +36,21 @@ using namespace o2::framework::expressions;
3636
3737namespace o2 ::aod
3838{
39- namespace pmdtrack
40- {
41- DECLARE_SOA_INDEX_COLUMN (Collision, collision);
42- DECLARE_SOA_ARRAY_INDEX_COLUMN (Collision, collisions);
43- DECLARE_SOA_INDEX_COLUMN (BC, bc);
44- DECLARE_SOA_SLICE_INDEX_COLUMN (Pmd,pmd);
45- } // namespace pmdtrack
46-
47- DECLARE_SOA_INDEX_TABLE_USER (PMDTracksIndex, BCs, " PMDTRKIDX" , pmdtrack::CollisionId, pmdtrack::BCId, pmdtrack::PmdIdSlice);
48- }
39+ namespace pmdtrack
40+ {
41+ DECLARE_SOA_INDEX_COLUMN (Collision, collision);
42+ DECLARE_SOA_ARRAY_INDEX_COLUMN (Collision, collisions);
43+ DECLARE_SOA_INDEX_COLUMN (BC, bc);
44+ DECLARE_SOA_SLICE_INDEX_COLUMN (Pmd, pmd);
45+ } // namespace pmdtrack
46+
47+ DECLARE_SOA_INDEX_TABLE_USER (PMDTracksIndex, BCs, " PMDTRKIDX" , pmdtrack::CollisionId, pmdtrack::BCId, pmdtrack::PmdIdSlice);
48+ } // namespace o2::aod
4949
5050struct BuiltPmdIndex {
5151 // build the index table PMDTracksIndex
5252 Builds<aod::PMDTracksIndex> idx;
53- void init (InitContext const &){};
53+ void init (InitContext const &) {};
5454};
5555
5656struct PmdQa {
@@ -59,12 +59,12 @@ struct PmdQa {
5959 ConfigurableAxis axisEventBin{" axisEventBin" , {4 , 0.5 , 4.5 }, " " };
6060 ConfigurableAxis axisVtxZBin{" axisVtxZBin" , {40 , -20 , 20 }, " " };
6161 ConfigurableAxis axisNPMDtracksBin{" axisNPMDtracksBin" , {500 , 0 , 500 }, " Number of pmdtracks" };
62- ConfigurableAxis axisClsxyBin{" axisClsxyBin" , {200 ,-100 ,100 }, " " };
63- ConfigurableAxis axisAdcBin{" axisAdcBin" , {200 ,0 , 2000 }, " " };
64- ConfigurableAxis axisEtaBin{" axisEtaBin" , {10 ,2.1 ,4.1 }, " " };
65- ConfigurableAxis axisNcellBin{" axisNcellBin" , {50 ,-0.5 ,49.5 }, " " };
62+ ConfigurableAxis axisClsxyBin{" axisClsxyBin" , {200 , -100 , 100 }, " " };
63+ ConfigurableAxis axisAdcBin{" axisAdcBin" , {200 , 0 , 2000 }, " " };
64+ ConfigurableAxis axisEtaBin{" axisEtaBin" , {10 , 2.1 , 4.1 }, " " };
65+ ConfigurableAxis axisNcellBin{" axisNcellBin" , {50 , -0.5 , 49.5 }, " " };
6666 Configurable<int > fMipCut {" fMipCut" , 432 , " fMipCut" };
67-
67+
6868 void init (InitContext&)
6969 {
7070
@@ -75,7 +75,7 @@ struct PmdQa {
7575 AxisSpec axisAdc = {axisAdcBin, " Adc" , " AdcAxis" };
7676 AxisSpec axisEta = {axisEtaBin, " Eta" , " EtaAxis" };
7777 AxisSpec axisNcell = {axisNcellBin, " Ncell" , " NcellAxis" };
78-
78+
7979 histos.add (" hEventHist" , " hEventHist" , kTH1F , {axisEvent});
8080 histos.add (" hVtxZHist" , " hVtxZHist" , kTH1F , {axisVtxZ});
8181 histos.add (" hNPMDtracks" , " Number of pmdtracks" , kTH1F , {axisNPMDtracks});
@@ -87,7 +87,7 @@ struct PmdQa {
8787
8888 using coltable = soa::Join<aod::Collisions, aod::PMDTracksIndex>;
8989 using colevsel = soa::Join<coltable, aod::EvSels>;
90-
90+
9191 void process (colevsel::iterator const & collision, aod::Pmds const &)
9292 {
9393 histos.fill (HIST (" hEventHist" ), 1 );
@@ -100,32 +100,32 @@ struct PmdQa {
100100 }
101101 histos.fill (HIST (" hEventHist" ), 3 );
102102 histos.fill (HIST (" hVtxZHist" ), collision.posZ ());
103-
103+
104104 if (collision.has_pmd ()) {
105105 histos.fill (HIST (" hEventHist" ), 4 );
106106 auto tracks = collision.pmd ();
107107 histos.fill (HIST (" hNPMDtracks" ), tracks.size ());
108108 for (const auto & track : tracks) {
109- if (track.pmddet () == 1 ) {
110- return ;
111- }
112- if (track.pmdclsz () == 0 ) {
113- return ;
114- }
115- if (!track.pmdmodule ()) {
116- return ;
117- }
118- histos.fill (HIST (" hClusXY" ), track.pmdclsx (), track.pmdclsy ());
119- histos.fill (HIST (" hClusAdc" ), track.pmdclsadc ());
120- float rdist = TMath::Sqrt (track.pmdclsx ()* track.pmdclsx () + track.pmdclsy ()* track.pmdclsy ());
121- float theta = TMath::ATan2 (rdist,track.pmdclsz ());
122- float etacls = -TMath::Log (TMath::Tan (0.5 * theta));
123- if (track.pmdclsadc () > fMipCut && track.pmdncell () > 2 ) {
124- if (etacls> 2.3 && etacls< 3.9 ){
125- histos.fill (HIST (" hetacls" ), etacls);
126- histos.fill (HIST (" hclsncell" ), track.pmdncell ());
127- }
128- }
109+ if (track.pmddet () == 1 ) {
110+ return ;
111+ }
112+ if (track.pmdclsz () == 0 ) {
113+ return ;
114+ }
115+ if (!track.pmdmodule ()) {
116+ return ;
117+ }
118+ histos.fill (HIST (" hClusXY" ), track.pmdclsx (), track.pmdclsy ());
119+ histos.fill (HIST (" hClusAdc" ), track.pmdclsadc ());
120+ float rdist = TMath::Sqrt (track.pmdclsx () * track.pmdclsx () + track.pmdclsy () * track.pmdclsy ());
121+ float theta = TMath::ATan2 (rdist, track.pmdclsz ());
122+ float etacls = -TMath::Log (TMath::Tan (0.5 * theta));
123+ if (track.pmdclsadc () > fMipCut && track.pmdncell () > 2 ) {
124+ if (etacls > 2.3 && etacls < 3.9 ) {
125+ histos.fill (HIST (" hetacls" ), etacls);
126+ histos.fill (HIST (" hclsncell" ), track.pmdncell ());
127+ }
128+ }
129129 }
130130 }
131131 }
0 commit comments