1717#include < TH1F.h>
1818
1919#include < array>
20- #include < iostream>
2120// O2 includes
2221#include " Common/Core/TrackSelection.h"
2322#include " Common/DataModel/Centrality.h"
3635#include " TDatabasePDG.h"
3736#include < TPDGCode.h>
3837
38+ #include < string>
3939#include < unordered_set>
40+ #include < vector>
4041using std::array;
4142using namespace o2 ;
4243using namespace o2 ::framework;
4344using namespace o2 ::framework::expressions;
4445
46+ using namespace std ;
47+
4548TH1D* tmpFqErr[6 ][5 ][52 ];
4649
4750struct FactorialMoments {
@@ -406,10 +409,7 @@ struct FactorialMoments {
406409 fqEvent = {{{{{0 , 0 , 0 , 0 , 0 , 0 }}}}};
407410 binConEvent = {{{0 , 0 , 0 , 0 , 0 }}};
408411 for (auto const & track : tracks) {
409- if (track.hasTPC ())
410- // if (track.hasITS())
411- // if (track.isGlobalTrack())
412- {
412+ if (track.hasTPC ()) {
413413 histos.fill (HIST (" mCollID" ), track.collisionId ());
414414 histos.fill (HIST (" mEta" ), track.eta ());
415415 histos.fill (HIST (" mPt" ), track.pt ());
@@ -434,7 +434,6 @@ struct FactorialMoments {
434434 }
435435 }
436436 for (auto iPt = 0 ; iPt < confNumPt; ++iPt) {
437- // if (countTracks[iPt] > 0)
438437 if (countTracks[iPt] > 0 ) {
439438 mHistArrQA [iPt * 4 + 3 ]->Fill (countTracks[iPt]);
440439 }
@@ -486,9 +485,6 @@ struct FactorialMoments {
486485 fqEvent = {{{{{0 , 0 , 0 , 0 , 0 , 0 }}}}};
487486 binConEvent = {{{0 , 0 , 0 , 0 , 0 }}};
488487 for (auto const & track : colltracks) {
489- // if (track.hasITS())
490- // if (track.hasTPC())
491- // if (track.isGlobalTrack()) {
492488 histos.fill (HIST (" mCollID" ), track.collisionId ());
493489 histos.fill (HIST (" mEta" ), track.eta ());
494490 histos.fill (HIST (" mPt" ), track.pt ());
@@ -510,10 +506,8 @@ struct FactorialMoments {
510506 histos.fill (HIST (" mSharedClsvsPt" ), track.pt (), track.tpcNClsShared ());
511507 histos.fill (HIST (" mSharedClsProbvsPt" ), track.pt (), track.tpcFractionSharedCls () / track.tpcNClsCrossedRows ());
512508 checkpT (track);
513- // }
514509 }
515510 for (auto iPt = 0 ; iPt < confNumPt; ++iPt) {
516- // if (countTracks[iPt] > 0)countTracks = {0, 0, 0, 0, 0};
517511 if (countTracks[iPt] > 0 ) {
518512 mHistArrQA [iPt * 4 + 3 ]->Fill (countTracks[iPt]);
519513 }
@@ -571,7 +565,6 @@ struct FactorialMoments {
571565 histos1.fill (HIST (" mPrimariesPerEvent" ), count);
572566 }
573567 for (auto iPt = 0 ; iPt < confNumPt; ++iPt) {
574- // if (countTracks[iPt] > 0)
575568 if (countTracks[iPt] > 0 ) {
576569 mHistArrQA [iPt * 4 + 3 ]->Fill (countTracks[iPt]);
577570 }
@@ -703,9 +696,6 @@ struct FactorialMoments {
703696 countTracks = {0 , 0 , 0 , 0 , 0 };
704697 fqEvent = {{{{{0 , 0 , 0 , 0 , 0 , 0 }}}}};
705698 binConEvent = {{{0 , 0 , 0 , 0 , 0 }}};
706- // using McParticle = std::decay_t<decltype(*particles.begin())>;
707- // std::vector<McParticle> selectedParticles;
708- // if (counter > 1) return;
709699 for (auto const & track : mcParticles) {
710700 histos1.fill (HIST (" hRecoPtBefore" ), track.pt ());
711701 histos1.fill (HIST (" hRecoEtaBefore" ), track.eta ());
0 commit comments