@@ -192,10 +192,10 @@ struct FactorialMomentsTask {
192192 for (int iPt = 0 ; iPt < numPt; ++iPt) {
193193 mHistArrQA .push_back (std::get<std::shared_ptr<TH1>>(histos.add (Form (" bin%i/mEta" , iPt + 1 ), Form (" #eta for bin %.2f-%.2f;#eta" , ptCuts.value [2 * iPt], ptCuts.value [2 * iPt + 1 ]), HistType::kTH1F , {{1000 , -2 , 2 }})));
194194 mHistArrQA .push_back (std::get<std::shared_ptr<TH1>>(histos.add (Form (" bin%i/mPt" , iPt + 1 ), Form (" pT for bin %.2f-%.2f;pT" , ptCuts.value [2 * iPt], ptCuts.value [2 * iPt + 1 ]), HistType::kTH1F , {axisPt[iPt]})));
195- mHistArrQA .push_back (std::get<std::shared_ptr<TH1>>(histos.add (Form (" bin%i/mPhi" , iPt + 1 ), Form (" #phi for bin %.2f-%.2f;#phi" , ptCuts.value [2 * iPt], ptCuts.value [2 * iPt + 1 ]), HistType::kTH1F , {{1000 , 0 , o2::constants::math::TwoPI}})));
195+ mHistArrQA .push_back (std::get<std::shared_ptr<TH1>>(histos.add (Form (" bin%i/mPhi" , iPt + 1 ), Form (" #phi for bin %.2f-%.2f;#phi" , ptCuts.value [2 * iPt], ptCuts.value [2 * iPt + 1 ]), HistType::kTH1F , {{1000 , 0 , o2::constants::math::TwoPI}})));
196196 mHistArrQA .push_back (std::get<std::shared_ptr<TH1>>(histos.add (Form (" bin%i/mMultiplicity" , iPt + 1 ), Form (" Multiplicity for bin %.2f-%.2f;Multiplicity" , ptCuts.value [2 * iPt], ptCuts.value [2 * iPt + 1 ]), HistType::kTH1F , {{1000 , 0 , 15000 }})));
197197 for (int iM = 0 ; iM < nBins; ++iM) {
198- auto mHistsR = std::get<std::shared_ptr<TH2>>(histos.add (Form (" bin%i/Reset/mEtaPhi%i" , iPt + 1 , iM), Form (" #eta#phi_%i for bin %.2f-%.2f;#eta;#phi" , iM, ptCuts.value [2 * iPt], ptCuts.value [2 * iPt + 1 ]), HistType::kTH2F , {{binningM[iM], -0.8 , 0.8 }, {binningM[iM], 0 , o2::constants::math::TwoPI}}));
198+ auto mHistsR = std::get<std::shared_ptr<TH2>>(histos.add (Form (" bin%i/Reset/mEtaPhi%i" , iPt + 1 , iM), Form (" #eta#phi_%i for bin %.2f-%.2f;#eta;#phi" , iM, ptCuts.value [2 * iPt], ptCuts.value [2 * iPt + 1 ]), HistType::kTH2F , {{binningM[iM], -0.8 , 0.8 }, {binningM[iM], 0 , o2::constants::math::TwoPI}}));
199199 mHistArrReset .push_back (mHistsR );
200200 for (int iq = 0 ; iq < nfqOrder; ++iq) {
201201 tmpFqErr[iq][iPt][iM] = new TH1D (Form (" tmpFqErr%i%i%i" , iq, iPt, iM), Form (" tmpFqErr%i%i%i" , iq, iPt, iM), 100 , 0 , 10 );
@@ -403,29 +403,29 @@ struct FactorialMomentsTask {
403403 binConEvent = {{{0 , 0 , 0 , 0 , 0 }}};
404404 for (auto const & track : colltracks) {
405405 // if (track.hasITS())
406- // if (track.hasTPC())
407- // if (track.isGlobalTrack()) {
408- histos.fill (HIST (" mCollID" ), track.collisionId ());
409- histos.fill (HIST (" mEta" ), track.eta ());
410- histos.fill (HIST (" mPt" ), track.pt ());
411- histos.fill (HIST (" mPhi" ), track.phi ());
412- histos.fill (HIST (" mNFindableClsTPC" ), track.tpcNClsFindable ());
413- histos.fill (HIST (" mNClsTPC" ), track.tpcNClsFound ());
414- histos.fill (HIST (" mNClsITS" ), track.itsNCls ());
415- histos.fill (HIST (" mChi2TPC" ), track.tpcChi2NCl ());
416- histos.fill (HIST (" mChi2ITS" ), track.itsChi2NCl ());
417- histos.fill (HIST (" mChi2TRD" ), track.trdChi2 ());
418- histos.fill (HIST (" mDCAxy" ), track.dcaXY ());
419- histos.fill (HIST (" mDCAx" ), track.dcaZ ());
420- histos.fill (HIST (" mDCAxyPt" ), track.pt (), track.dcaXY ());
421- histos.fill (HIST (" mDCAzPt" ), track.pt (), track.dcaZ ());
422- histos.fill (HIST (" mNSharedClsTPC" ), track.tpcNClsShared ());
423- histos.fill (HIST (" mCrossedRowsTPC" ), track.tpcNClsCrossedRows ());
424- histos.fill (HIST (" mNFinClsminusCRows" ), track.tpcNClsFindableMinusCrossedRows ());
425- histos.fill (HIST (" mNFractionShClsTPC" ), track.tpcFractionSharedCls ());
426- histos.fill (HIST (" mSharedClsvsPt" ), track.pt (), track.tpcNClsShared ());
427- histos.fill (HIST (" mSharedClsProbvsPt" ), track.pt (), track.tpcFractionSharedCls () / track.tpcNClsCrossedRows ());
428- checkpT (track);
406+ // if (track.hasTPC())
407+ // if (track.isGlobalTrack()) {
408+ histos.fill (HIST (" mCollID" ), track.collisionId ());
409+ histos.fill (HIST (" mEta" ), track.eta ());
410+ histos.fill (HIST (" mPt" ), track.pt ());
411+ histos.fill (HIST (" mPhi" ), track.phi ());
412+ histos.fill (HIST (" mNFindableClsTPC" ), track.tpcNClsFindable ());
413+ histos.fill (HIST (" mNClsTPC" ), track.tpcNClsFound ());
414+ histos.fill (HIST (" mNClsITS" ), track.itsNCls ());
415+ histos.fill (HIST (" mChi2TPC" ), track.tpcChi2NCl ());
416+ histos.fill (HIST (" mChi2ITS" ), track.itsChi2NCl ());
417+ histos.fill (HIST (" mChi2TRD" ), track.trdChi2 ());
418+ histos.fill (HIST (" mDCAxy" ), track.dcaXY ());
419+ histos.fill (HIST (" mDCAx" ), track.dcaZ ());
420+ histos.fill (HIST (" mDCAxyPt" ), track.pt (), track.dcaXY ());
421+ histos.fill (HIST (" mDCAzPt" ), track.pt (), track.dcaZ ());
422+ histos.fill (HIST (" mNSharedClsTPC" ), track.tpcNClsShared ());
423+ histos.fill (HIST (" mCrossedRowsTPC" ), track.tpcNClsCrossedRows ());
424+ histos.fill (HIST (" mNFinClsminusCRows" ), track.tpcNClsFindableMinusCrossedRows ());
425+ histos.fill (HIST (" mNFractionShClsTPC" ), track.tpcFractionSharedCls ());
426+ histos.fill (HIST (" mSharedClsvsPt" ), track.pt (), track.tpcNClsShared ());
427+ histos.fill (HIST (" mSharedClsProbvsPt" ), track.pt (), track.tpcFractionSharedCls () / track.tpcNClsCrossedRows ());
428+ checkpT (track);
429429 // }
430430 }
431431 auto mcParts = mcParticles.sliceBy (perMcCollision, coll.mcCollision ().globalIndex ());
@@ -463,13 +463,13 @@ struct FactorialMomentsTask {
463463 using TracksRecSim = soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::McTrackLabels>;
464464 using CollisionRecSimRun2 = soa::Filtered<soa::Join<aod::Collisions, aod::McCollisionLabels, EventSelectionrun2>>::iterator;
465465 using BCsWithRun2Info = soa::Join<aod::BCs, aod::Run2BCInfos, aod::Timestamps>;
466- void processMcRun2 (CollisionRecSimRun2 const & coll,
467- aod::BCs const & ,
468- TracksRecSim const & tracks,
469- aod::McParticles const & mcParticles,
470- aod::McCollisions const &,
471- BCsWithRun2Info const & )
472- {
466+ void processMcRun2 (CollisionRecSimRun2 const & coll,
467+ aod::BCs const &,
468+ TracksRecSim const & tracks,
469+ aod::McParticles const & mcParticles,
470+ aod::McCollisions const &,
471+ BCsWithRun2Info const &)
472+ {
473473 auto bc = coll.bc_as <BCsWithRun2Info>();
474474 if (!(bc.eventCuts () & BIT (aod::Run2EventCuts::kAliEventCutsAccepted ))) {
475475 return ;
0 commit comments