@@ -335,10 +335,10 @@ struct lnnRecoTask {
335335 hPtItsTpcTofTrStr = qaRegistry.add <TH2>(" MC/McTracks/hPtItsTpcTofTrStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
336336 hPtTrkItsTpcTofTrStr = qaRegistry.add <TH2>(" MC/McTracks/hPtTrkItsTpcTofTrStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
337337 hPtItsTpcPiStr = qaRegistry.add <TH2>(" MC/McTracks/hPtItsTpcPiStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
338- hPtTrkItsTpcPiStr = qaRegistry.add <TH2>(" MC/McTracks/hPtTrkItsTpcPiStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
339- hPtGeneratedTrStr = qaRegistry.add <TH2>(" MC/McGen/hPtGeneratedTrStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
340- hPtGeneratedPiStr = qaRegistry.add <TH2>(" MC/McGen/hPtGeneratedPiStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
341-
338+ hPtTrkItsTpcPiStr = qaRegistry.add <TH2>(" MC/McTracks/hPtTrkItsTpcPiStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
339+ hPtGeneratedTrStr = qaRegistry.add <TH2>(" MC/McGen/hPtGeneratedTrStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
340+ hPtGeneratedPiStr = qaRegistry.add <TH2>(" MC/McGen/hPtGeneratedPiStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
341+
342342 hEvents->GetXaxis ()->SetBinLabel (1 , " All" );
343343 hEvents->GetXaxis ()->SetBinLabel (2 , " sel8" );
344344 hEvents->GetXaxis ()->SetBinLabel (3 , " z_{vtx}" );
@@ -489,7 +489,7 @@ struct lnnRecoTask {
489489 h2FT0CnTPCNClsCrossedRowsPiBfSel->Fill (collision.centFT0C (), pitrack.tpcNClsCrossedRows ());
490490 h2FT0Cchi2NClTPCtoTrBfSel->Fill (collision.centFT0C (), h3track.tpcChi2NCl ());
491491 h2FT0Cchi2NClITStoTrBfSel->Fill (collision.centFT0C (), h3track.itsChi2NCl ());
492-
492+
493493 if (doTrackQA) {
494494 bool passedTrTrackITS = h3track.hasITS ();
495495 bool passedTrTrackTOF = h3track.hasTOF ();
@@ -697,9 +697,9 @@ struct lnnRecoTask {
697697
698698 bool isLnnDecay (aod::McParticles::iterator const & mcPart)
699699 {
700- if (!mcPart.has_mothers ())
700+ if (!mcPart.has_mothers ())
701701 return false ;
702-
702+
703703 if (mcPart.getProcess () != 4 )
704704 return false ;
705705
@@ -713,9 +713,10 @@ struct lnnRecoTask {
713713 }
714714 return motherIsAccepted;
715715 }
716-
716+
717717 template <class Tcoll >
718- void fillMcHistograms (TracksFullMC::iterator const & mcTrack, aod::McParticles::iterator const & mc, Tcoll const & collision) {
718+ void fillMcHistograms (TracksFullMC::iterator const & mcTrack, aod::McParticles::iterator const & mc, Tcoll const & collision)
719+ {
719720
720721 bool passedTrackITS = mcTrack.hasITS ();
721722 bool passedTrackTPC = mcTrack.hasTPC ();
@@ -743,7 +744,7 @@ struct lnnRecoTask {
743744 hPtTrkItsTpcPiStr->Fill (collision.centFT0C (), mcTrack.pt ());
744745 }
745746 }
746- }
747+ }
747748
748749 void processData (CollisionsFull const & collisions, aod::V0s const & V0s, TracksFull const & tracks, aod::BCsWithTimestamps const &)
749750 {
@@ -767,7 +768,6 @@ struct lnnRecoTask {
767768 }
768769 hEvents->Fill (3 .);
769770
770-
771771 hZvtx->Fill (collision.posZ ());
772772 hCentFT0A->Fill (collision.centFT0A ());
773773 hCentFT0C->Fill (collision.centFT0C ());
@@ -798,7 +798,7 @@ struct lnnRecoTask {
798798 PROCESS_SWITCH (lnnRecoTask, processData, " Data analysis" , true );
799799
800800 // MC process
801- void processMC (CollisionsFullMC const & collisions, aod::McCollisions const & mcCollisions, aod::V0s const & V0s, aod::BCsWithTimestamps const &, TracksFull const & tracks, TracksFullMC const & tracksMC, aod::McTrackLabels const & trackLabelsMC, aod::McParticles const & particlesMC)
801+ void processMC (CollisionsFullMC const & collisions, aod::McCollisions const & mcCollisions, aod::V0s const & V0s, aod::BCsWithTimestamps const &, TracksFull const & tracks, TracksFullMC const & tracksMC, aod::McTrackLabels const & trackLabelsMC, aod::McParticles const & particlesMC)
802802 {
803803 filledMothers.clear ();
804804 isGoodCollision.clear ();
@@ -811,11 +811,11 @@ struct lnnRecoTask {
811811 for (const auto & collision : collisions) {
812812 for (auto const & trackMC : tracksMC) {
813813
814- if (!trackMC.has_mcParticle ())
814+ if (!trackMC.has_mcParticle ())
815815 continue ;
816816
817817 auto mc = trackMC.mcParticle ();
818-
818+
819819 if (!isLnnDecay (mc))
820820 continue ;
821821
@@ -830,7 +830,7 @@ struct lnnRecoTask {
830830 if (collision.has_mcCollision ()) {
831831 recoCollisionIds[collision.mcCollisionId ()] = collision.globalIndex ();
832832 }
833-
833+
834834 if ((!collision.sel8 ())) {
835835 continue ;
836836 }
@@ -843,7 +843,7 @@ struct lnnRecoTask {
843843 continue ;
844844 }
845845 hEvents->Fill (3 .);
846-
846+
847847 hZvtx->Fill (collision.posZ ());
848848 hCentFT0A->Fill (collision.centFT0A ());
849849 hCentFT0C->Fill (collision.centFT0C ());
@@ -889,7 +889,7 @@ struct lnnRecoTask {
889889
890890 // now we fill only the signal candidates that were not reconstructed
891891 for (const auto & mcPart : particlesMC) {
892-
892+
893893 if (std::abs (mcPart.pdgCode ()) != lnnPdg) {
894894 continue ;
895895 }
@@ -958,7 +958,7 @@ struct lnnRecoTask {
958958 centFT0A = recoCollision.centFT0A ();
959959 centFT0C = recoCollision.centFT0C ();
960960 centFT0M = recoCollision.centFT0M ();
961- }
961+ }
962962 outputMCTable (centFT0A, centFT0C, centFT0M,
963963 -1 , -1 , -1 ,
964964 0 ,
@@ -972,7 +972,7 @@ struct lnnRecoTask {
972972 -1 , -1 , -1 ,
973973 chargeFactor * lnnCand.genPt (), lnnCand.genPhi (), lnnCand.genEta (), lnnCand.genPt3H (),
974974 lnnCand.gDecVtx [0 ], lnnCand.gDecVtx [1 ], lnnCand.gDecVtx [2 ], lnnCand.isReco , lnnCand.isSignal , lnnCand.recoMcColl , lnnCand.survEvSelection );
975- }
975+ }
976976 }
977977 }
978978 PROCESS_SWITCH (lnnRecoTask, processMC, " MC analysis" , false );
@@ -983,4 +983,5 @@ WorkflowSpec
983983{
984984 return WorkflowSpec{
985985 adaptAnalysisTask<lnnRecoTask>(cfgc)};
986- }
986+ }
987+
0 commit comments