99// granted to it by virtue of its status as an Intergovernmental Organization
1010// or submit itself to any jurisdiction.
1111// /
12- // / \file pmdqa .cxx
12+ // / \file pmdQa .cxx
1313// /
1414// / \brief QA task to check PMD info on Run 2 converted data
1515// / \author Abhi Modak (abhi.modak@cern.ch)
@@ -68,6 +68,7 @@ struct PmdQa {
6868 Configurable<int > fNcellCut {" fNcellCut" , 2 , " fNcellCut" };
6969 Configurable<float > fEtalow {" fEtalow" , 2.3 , " fEtalow" };
7070 Configurable<float > fEtahigh {" fEtalow" , 3.9 , " fEtahigh" };
71+ Configurable<float > fVtxCut {" fVtxCut" , 10.0 , " fVtxCut" };
7172
7273 void init (InitContext&)
7374 {
@@ -89,17 +90,17 @@ struct PmdQa {
8990 histos.add (" hclsncell" , " hclsncell" , kTH1F , {axisNcell});
9091 }
9192
92- using coltable = soa::Join<aod::Collisions, aod::PMDTracksIndex>;
93- using colevsel = soa::Join<coltable , aod::EvSels>;
93+ using ColTable = soa::Join<aod::Collisions, aod::PMDTracksIndex>;
94+ using ColevSel = soa::Join<ColTable , aod::EvSels>;
9495
95- void process (colevsel ::iterator const & collision, aod::Pmds const &)
96+ void process (ColevSel ::iterator const & collision, aod::Pmds const &)
9697 {
9798 histos.fill (HIST (" hEventHist" ), 1 );
9899 if (collision.sel7 ()) {
99100 return ;
100101 }
101102 histos.fill (HIST (" hEventHist" ), 2 );
102- if (std::abs (collision.posZ ()) >= 10 . ) {
103+ if (std::abs (collision.posZ ()) >= fVtxCut ) {
103104 return ;
104105 }
105106 histos.fill (HIST (" hEventHist" ), 3 );
@@ -121,9 +122,9 @@ struct PmdQa {
121122 }
122123 histos.fill (HIST (" hClusXY" ), track.pmdclsx (), track.pmdclsy ());
123124 histos.fill (HIST (" hClusAdc" ), track.pmdclsadc ());
124- float rdist = TMath::Sqrt (track.pmdclsx () * track.pmdclsx () + track.pmdclsy () * track.pmdclsy ());
125- float theta = TMath::ATan2 (rdist, track.pmdclsz ());
126- float etacls = -TMath::Log (TMath::Tan (0.5 * theta));
125+ float rdist = std::sqrt (track.pmdclsx () * track.pmdclsx () + track.pmdclsy () * track.pmdclsy ());
126+ float theta = std::atan2 (rdist, track.pmdclsz ());
127+ float etacls = -std::log (TMath::Tan (0.5 * theta));
127128 if (track.pmdclsadc () > fMipCut && track.pmdncell () > fNcellCut ) {
128129 if (etacls > fEtalow && etacls < fEtahigh ) {
129130 histos.fill (HIST (" hetacls" ), etacls);
0 commit comments