@@ -31,59 +31,64 @@ using namespace o2::framework::expressions;
3131struct propagatorQa {
3232 Configurable<float > windowDCA{" windowDCA" , 50 , " windowDCA" };
3333 Configurable<int > Nbins{" Nbins" , 10000 , " Nbins" };
34-
35- // Momentum distribution
34+
35+ // Momentum distribution
3636 OutputObj<TH1F> hPtBeforeXcut{TH1F (" hPtBeforeXcut" , " hPtBeforeXcut" , Nbins, 0 , 10 )};
3737 OutputObj<TH1F> hPt{TH1F (" hPt" , " hPt" , Nbins, 0 , 10 )};
3838 OutputObj<TH1F> hPtusedInSVertexer{TH1F (" hPtusedInSVertexer" , " hPtusedInSVertexer" , Nbins, 0 , 10 )};
39-
40- // IU radii
39+
40+ // IU radii
4141 OutputObj<TH1F> hUpdateRadii{TH1F (" hUpdateRadii" , " hUpdateRadii" , 5000 , 0 , 100 )};
4242 OutputObj<TH1F> hUpdateRadiiusedInSVertexer{TH1F (" hUpdateRadiiusedInSVertexer" , " hUpdateRadii" , 5000 , 0 , 100 )};
43-
44- // DCA
43+
44+ // DCA
4545 OutputObj<TH1F> hdcaXYall{TH1F (" hdcaXYall" , " hdcaXYall" , Nbins, -windowDCA, windowDCA)};
4646 OutputObj<TH1F> hdcaXYusedInSVertexer{TH1F (" hdcaXYusedInSVertexer" , " hdcaXYusedInSVertexer" , Nbins, -windowDCA, windowDCA)};
47-
47+
4848 o2::track::TrackPar lTrackParametrization;
49-
49+
5050 void process (aod::Collision const & collision, aod::V0s const & V0s, aod::Cascades const & cascades, soa::Join<aod::TracksIU, aod::TracksExtra, aod::TracksDCA> const & tracks)
5151 {
5252 for (auto & track : tracks) {
5353 hPtBeforeXcut->Fill (track.pt ());
54-
54+
5555 float maxXtoConsider = o2::constants::geom::XTPCInnerRef + 0.1 ;
56- if ( track.trackType () != aod::track::TrackIU && track.x () > maxXtoConsider ) continue ;
57-
56+ if (track.trackType () != aod::track::TrackIU && track.x () > maxXtoConsider)
57+ continue ;
58+
5859 std::array<float , 3 > pos;
5960 lTrackParametrization = getTrackPar (track);
6061 lTrackParametrization.getXYZGlo (pos);
61- float lRadiusOfLastUpdate = TMath::Sqrt (pos[0 ]* pos[0 ]+ pos[1 ]* pos[1 ]);
62-
62+ float lRadiusOfLastUpdate = TMath::Sqrt (pos[0 ] * pos[0 ] + pos[1 ] * pos[1 ]);
63+
6364 hUpdateRadii->Fill (lRadiusOfLastUpdate);
64-
65+
6566 // fill kinematic variables
6667 hPt->Fill (track.pt ());
67-
68+
6869 float lDCA = track.dcaXY ();
69-
70+
7071 hdcaXYall->Fill (lDCA);
71-
72- // determine if track was used in svertexer
72+
73+ // determine if track was used in svertexer
7374 bool usedInSVertexer = false ;
7475 bool lPosUsed, lNegUsed, lBachUsed;
7576 for (auto & V0 : V0s) {
76- lPosUsed = ( V0.posTrackId () == track.globalIndex () );
77- lNegUsed = ( V0.negTrackId () == track.globalIndex () );
77+ lPosUsed = (V0.posTrackId () == track.globalIndex ());
78+ lNegUsed = (V0.negTrackId () == track.globalIndex ());
7879 }
7980 for (auto & cascade : cascades) {
80- lBachUsed = ( cascade.bachelorId () == track.globalIndex () );
81+ lBachUsed = (cascade.bachelorId () == track.globalIndex ());
8182 }
82- if ( lPosUsed || lNegUsed || lBachUsed ) usedInSVertexer = true ;
83-
84- if ( usedInSVertexer ) hUpdateRadiiusedInSVertexer->Fill (lRadiusOfLastUpdate);
85- if ( usedInSVertexer ) hdcaXYusedInSVertexer->Fill (lDCA);
86- if ( usedInSVertexer ) hPtusedInSVertexer->Fill (track.pt ());
83+ if (lPosUsed || lNegUsed || lBachUsed)
84+ usedInSVertexer = true ;
85+
86+ if (usedInSVertexer)
87+ hUpdateRadiiusedInSVertexer->Fill (lRadiusOfLastUpdate);
88+ if (usedInSVertexer)
89+ hdcaXYusedInSVertexer->Fill (lDCA);
90+ if (usedInSVertexer)
91+ hPtusedInSVertexer->Fill (track.pt ());
8792 }
8893 }
8994};
0 commit comments