@@ -331,10 +331,10 @@ struct lnnRecoTask {
331331 hPtItsTpcTofTrStr = qaRegistry.add <TH2>(" MC/McTracks/hPtItsTpcTofTrStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
332332 hPtTrkItsTpcTofTrStr = qaRegistry.add <TH2>(" MC/McTracks/hPtTrkItsTpcTofTrStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
333333 hPtItsTpcPiStr = qaRegistry.add <TH2>(" MC/McTracks/hPtItsTpcPiStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
334- hPtTrkItsTpcPiStr = qaRegistry.add <TH2>(" MC/McTracks/hPtTrkItsTpcPiStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
335- hPtGeneratedTrStr = qaRegistry.add <TH2>(" MC/McGen/hPtGeneratedTrStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
336- hPtGeneratedPiStr = qaRegistry.add <TH2>(" MC/McGen/hPtGeneratedPiStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
337-
334+ hPtTrkItsTpcPiStr = qaRegistry.add <TH2>(" MC/McTracks/hPtTrkItsTpcPiStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
335+ hPtGeneratedTrStr = qaRegistry.add <TH2>(" MC/McGen/hPtGeneratedTrStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
336+ hPtGeneratedPiStr = qaRegistry.add <TH2>(" MC/McGen/hPtGeneratedPiStr" , " ;FT0C (%);#it{p}_{T} (GeV/#it{c})" , HistType::kTH2F , {{10 , 0 , 100 }, {160 , 0 , 8 }});
337+
338338 hEvents->GetXaxis ()->SetBinLabel (1 , " All" );
339339 hEvents->GetXaxis ()->SetBinLabel (2 , " sel8" );
340340 hEvents->GetXaxis ()->SetBinLabel (3 , " z_{vtx}" );
@@ -486,7 +486,7 @@ struct lnnRecoTask {
486486 h2FT0CnClusTPCtoPiBfSel->Fill (collision.centFT0C (), pitrack.tpcNClsFound ());
487487 h2FT0Cchi2NClTPCtoTrBfSel->Fill (collision.centFT0C (), h3track.tpcChi2NCl ());
488488 h2FT0Cchi2NClITStoTrBfSel->Fill (collision.centFT0C (), h3track.itsChi2NCl ());
489-
489+
490490 if (doTrackQA) {
491491 bool passedTrTrackITS = h3track.hasITS ();
492492 bool passedTrTrackTOF = h3track.hasTOF ();
@@ -694,25 +694,26 @@ struct lnnRecoTask {
694694
695695 bool isLnnDecay (aod::McParticles::iterator const & mcPart)
696696 {
697- if (!mcPart.has_mothers ())
697+ if (!mcPart.has_mothers ())
698698 return false ;
699699
700700 bool motherIsAccepted = true ;
701701 auto mothers = mcPart.mothers_as <o2::aod::McParticles>();
702702 for (const auto & mother : mothers) {
703703 if (mcPart.getProcess () == 4 ) {
704- motherIsAccepted = false ;
705- // só aceita se a mãe for um Lnn
706- if (std::abs (mother.pdgCode ()) == lnnPdg) {
707- motherIsAccepted = true ;
708- }
704+ motherIsAccepted = false ;
705+ // só aceita se a mãe for um Lnn
706+ if (std::abs (mother.pdgCode ()) == lnnPdg) {
707+ motherIsAccepted = true ;
709708 }
709+ }
710710 }
711711 return motherIsAccepted;
712712 }
713-
713+
714714 template <class Tcoll >
715- void fillMcHistograms (TracksFullMC::iterator const & mcTrack, aod::McParticles::iterator const & mc, bool motherIsAccepted, Tcoll const & collision) {
715+ void fillMcHistograms (TracksFullMC::iterator const & mcTrack, aod::McParticles::iterator const & mc, bool motherIsAccepted, Tcoll const & collision)
716+ {
716717
717718 bool passedTrackITS = mcTrack.hasITS ();
718719 bool passedTrackTPC = mcTrack.hasTPC ();
@@ -741,7 +742,7 @@ struct lnnRecoTask {
741742 hPtTrkItsTpcPiStr->Fill (collision.centFT0C (), mcTrack.pt ());
742743 }
743744 }
744- }
745+ }
745746
746747 void processData (CollisionsFull const & collisions, aod::V0s const & V0s, TracksFull const & tracks, aod::BCsWithTimestamps const &)
747748 {
@@ -791,7 +792,7 @@ struct lnnRecoTask {
791792 PROCESS_SWITCH (lnnRecoTask, processData, " Data analysis" , true );
792793
793794 // MC process
794- 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)
795+ 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)
795796 {
796797 filledMothers.clear ();
797798 isGoodCollision.clear ();
@@ -804,14 +805,14 @@ struct lnnRecoTask {
804805 for (const auto & collision : collisions) {
805806 for (auto const & trackMC : tracksMC) {
806807
807- if (!trackMC.has_mcParticle ())
808- continue ;
808+ if (!trackMC.has_mcParticle ())
809+ continue ;
809810
810- auto mc = trackMC.mcParticle ();
811- bool motherIsAccepted = isLnnDecay (mc);
811+ auto mc = trackMC.mcParticle ();
812+ bool motherIsAccepted = isLnnDecay (mc);
812813
813- fillMcHistograms (trackMC, mc, motherIsAccepted, collision);
814- }
814+ fillMcHistograms (trackMC, mc, motherIsAccepted, collision);
815+ }
815816
816817 lnnCandidates.clear ();
817818 auto bc = collision.bc_as <aod::BCsWithTimestamps>();
@@ -829,7 +830,7 @@ struct lnnRecoTask {
829830 continue ;
830831 }
831832 hEvents->Fill (2 .);
832-
833+
833834 hZvtx->Fill (collision.posZ ());
834835 hCentFT0A->Fill (collision.centFT0A ());
835836 hCentFT0C->Fill (collision.centFT0C ());
@@ -884,7 +885,7 @@ struct lnnRecoTask {
884885
885886 // now we fill only the signal candidates that were not reconstructed
886887 for (const auto & mcPart : particlesMC) {
887-
888+
888889 if (std::abs (mcPart.pdgCode ()) != lnnPdg) {
889890 continue ;
890891 }
@@ -940,7 +941,7 @@ struct lnnRecoTask {
940941 lnnCand.posTrackID = -1 ;
941942 lnnCand.negTrackID = -1 ;
942943 lnnCand.isSignal = true ;
943- if (lnnCand.isSignal ){
944+ if (lnnCand.isSignal ) {
944945 h2FT0CPtGenColGenCandMC->Fill (cent, chargeFactor * lnnCand.genPt ());
945946 h2FT0CPtGenColGenTrStrMC->Fill (cent, chargeFactor * lnnCand.genPt3H ());
946947 h2FT0CPtGenColGenPiStrMC->Fill (cent, chargeFactor * lnnCand.genPtPi ());
@@ -950,7 +951,7 @@ struct lnnRecoTask {
950951 if (lnnCand.recoMcColl ) {
951952 auto recoCollision = collisions.rawIteratorAt (recoCollisionIds[mcPart.mcCollisionId ()]);
952953 centFT0C = recoCollision.centFT0C ();
953- }
954+ }
954955 outputMCTable (-1 , centFT0C, -1 ,
955956 -1 , -1 , -1 ,
956957 0 ,
@@ -964,7 +965,7 @@ struct lnnRecoTask {
964965 -1 , -1 , -1 ,
965966 chargeFactor * lnnCand.genPt (), lnnCand.genPhi (), lnnCand.genEta (), lnnCand.genPt3H (),
966967 lnnCand.gDecVtx [0 ], lnnCand.gDecVtx [1 ], lnnCand.gDecVtx [2 ], lnnCand.isReco , lnnCand.isSignal , lnnCand.recoMcColl , lnnCand.survEvSelection );
967- }
968+ }
968969 }
969970 }
970971 PROCESS_SWITCH (lnnRecoTask, processMC, " MC analysis" , false );
0 commit comments