@@ -657,7 +657,7 @@ struct FlowGfwV02 {
657657 }
658658
659659 template <typename TTrack>
660- double getEfficiency (TTrack track, const int & pid= PidCharged)
660+ double getEfficiency (TTrack track, const int & pid = PidCharged)
661661 {
662662 double eff = 1 .;
663663 if (cfg.mEfficiency [pid])
@@ -668,7 +668,6 @@ struct FlowGfwV02 {
668668 return 1 . / eff;
669669 }
670670
671-
672671 template <typename TCollision>
673672 bool eventSelected (TCollision collision, const int & multTrk, const float & centrality)
674673 {
@@ -1027,62 +1026,63 @@ struct FlowGfwV02 {
10271026 }
10281027
10291028 template <QAFillTime ft, typename TTrack>
1030- inline void fillPidQA (TTrack track, const int & pid) {
1029+ inline void fillPidQA (TTrack track, const int & pid)
1030+ {
10311031 // Fill Nsigma QA
10321032 if (!cfgUseItsPID) {
10331033 if (ft == kBefore ) {
1034- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" TofTpcNsigma_pions" ), track.tpcNSigmaPi (), track.tofNSigmaPi (), track.pt ());
1034+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" TofTpcNsigma_pions" ), track.tpcNSigmaPi (), track.tofNSigmaPi (), track.pt ());
10351035 if (cfgGetdEdx) {
10361036 double tpcExpSignalPi = track.tpcSignal () - (track.tpcNSigmaPi () * track.tpcExpSigmaPi ());
10371037
1038- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" TpcdEdx_ptwise_pions" ), track.pt (), track.tpcSignal (), track.tofNSigmaPi ());
1039- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" ExpTpcdEdx_ptwise_pions" ), track.pt (), tpcExpSignalPi, track.tofNSigmaPi ());
1040- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" ExpSigma_ptwise_pions" ), track.pt (), track.tpcExpSigmaPi (), track.tofNSigmaPi ());
1038+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" TpcdEdx_ptwise_pions" ), track.pt (), track.tpcSignal (), track.tofNSigmaPi ());
1039+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" ExpTpcdEdx_ptwise_pions" ), track.pt (), tpcExpSignalPi, track.tofNSigmaPi ());
1040+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" ExpSigma_ptwise_pions" ), track.pt (), track.tpcExpSigmaPi (), track.tofNSigmaPi ());
10411041 }
1042- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" TofTpcNsigma_kaons" ), track.tpcNSigmaKa (), track.tofNSigmaKa (), track.pt ());
1042+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" TofTpcNsigma_kaons" ), track.tpcNSigmaKa (), track.tofNSigmaKa (), track.pt ());
10431043 if (cfgGetdEdx) {
10441044 double tpcExpSignalKa = track.tpcSignal () - (track.tpcNSigmaKa () * track.tpcExpSigmaKa ());
10451045
1046- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" TpcdEdx_ptwise_kaons" ), track.pt (), track.tpcSignal (), track.tofNSigmaKa ());
1047- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" ExpTpcdEdx_ptwise_kaons" ), track.pt (), tpcExpSignalKa, track.tofNSigmaKa ());
1048- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" ExpSigma_ptwise_kaons" ), track.pt (), track.tpcExpSigmaKa (), track.tofNSigmaKa ());
1046+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" TpcdEdx_ptwise_kaons" ), track.pt (), track.tpcSignal (), track.tofNSigmaKa ());
1047+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" ExpTpcdEdx_ptwise_kaons" ), track.pt (), tpcExpSignalKa, track.tofNSigmaKa ());
1048+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" ExpSigma_ptwise_kaons" ), track.pt (), track.tpcExpSigmaKa (), track.tofNSigmaKa ());
10491049 }
1050- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" TofTpcNsigma_protons" ), track.tpcNSigmaPr (), track.tofNSigmaPr (), track.pt ());
1050+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" TofTpcNsigma_protons" ), track.tpcNSigmaPr (), track.tofNSigmaPr (), track.pt ());
10511051 if (cfgGetdEdx) {
10521052 double tpcExpSignalPr = track.tpcSignal () - (track.tpcNSigmaPr () * track.tpcExpSigmaPr ());
10531053
1054- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" TpcdEdx_ptwise_protons" ), track.pt (), track.tpcSignal (), track.tofNSigmaPr ());
1055- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" ExpTpcdEdx_ptwise_protons" ), track.pt (), tpcExpSignalPr, track.tofNSigmaPr ());
1056- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" ExpSigma_ptwise_protons" ), track.pt (), track.tpcExpSigmaPr (), track.tofNSigmaPr ());
1054+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" TpcdEdx_ptwise_protons" ), track.pt (), track.tpcSignal (), track.tofNSigmaPr ());
1055+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" ExpTpcdEdx_ptwise_protons" ), track.pt (), tpcExpSignalPr, track.tofNSigmaPr ());
1056+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" ExpSigma_ptwise_protons" ), track.pt (), track.tpcExpSigmaPr (), track.tofNSigmaPr ());
10571057 }
10581058 } else if (ft == kAfter ) {
10591059 if (pid == PidPions) {
1060- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" TofTpcNsigma_pions" ), track.tpcNSigmaPi (), track.tofNSigmaPi (), track.pt ());
1060+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" TofTpcNsigma_pions" ), track.tpcNSigmaPi (), track.tofNSigmaPi (), track.pt ());
10611061 if (cfgGetdEdx) {
10621062 double tpcExpSignalPi = track.tpcSignal () - (track.tpcNSigmaPi () * track.tpcExpSigmaPi ());
10631063
1064- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" TpcdEdx_ptwise_pions" ), track.pt (), track.tpcSignal (), track.tofNSigmaPi ());
1065- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" ExpTpcdEdx_ptwise_pions" ), track.pt (), tpcExpSignalPi, track.tofNSigmaPi ());
1066- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" ExpSigma_ptwise_pions" ), track.pt (), track.tpcExpSigmaPi (), track.tofNSigmaPi ());
1064+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" TpcdEdx_ptwise_pions" ), track.pt (), track.tpcSignal (), track.tofNSigmaPi ());
1065+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" ExpTpcdEdx_ptwise_pions" ), track.pt (), tpcExpSignalPi, track.tofNSigmaPi ());
1066+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" ExpSigma_ptwise_pions" ), track.pt (), track.tpcExpSigmaPi (), track.tofNSigmaPi ());
10671067 }
10681068 }
10691069 if (pid == PidKaons) {
1070- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" TofTpcNsigma_kaons" ), track.tpcNSigmaKa (), track.tofNSigmaKa (), track.pt ());
1071- if (cfgGetdEdx) {
1072- double tpcExpSignalKa = track.tpcSignal () - (track.tpcNSigmaKa () * track.tpcExpSigmaKa ());
1070+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" TofTpcNsigma_kaons" ), track.tpcNSigmaKa (), track.tofNSigmaKa (), track.pt ());
1071+ if (cfgGetdEdx) {
1072+ double tpcExpSignalKa = track.tpcSignal () - (track.tpcNSigmaKa () * track.tpcExpSigmaKa ());
10731073
1074- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" TpcdEdx_ptwise_kaons" ), track.pt (), track.tpcSignal (), track.tofNSigmaKa ());
1075- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" ExpTpcdEdx_ptwise_kaons" ), track.pt (), tpcExpSignalKa, track.tofNSigmaKa ());
1076- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" ExpSigma_ptwise_kaons" ), track.pt (), track.tpcExpSigmaKa (), track.tofNSigmaKa ());
1074+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" TpcdEdx_ptwise_kaons" ), track.pt (), track.tpcSignal (), track.tofNSigmaKa ());
1075+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" ExpTpcdEdx_ptwise_kaons" ), track.pt (), tpcExpSignalKa, track.tofNSigmaKa ());
1076+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" ExpSigma_ptwise_kaons" ), track.pt (), track.tpcExpSigmaKa (), track.tofNSigmaKa ());
10771077 }
10781078 }
10791079 if (pid == PidProtons) {
1080- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" TofTpcNsigma_protons" ), track.tpcNSigmaPr (), track.tofNSigmaPr (), track.pt ());
1080+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" TofTpcNsigma_protons" ), track.tpcNSigmaPr (), track.tofNSigmaPr (), track.pt ());
10811081 if (cfgGetdEdx) {
10821082 double tpcExpSignalPr = track.tpcSignal () - (track.tpcNSigmaPr () * track.tpcExpSigmaPr ());
1083- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" TpcdEdx_ptwise_protons" ), track.pt (), track.tpcSignal (), track.tofNSigmaPr ());
1084- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" ExpTpcdEdx_ptwise_protons" ), track.pt (), tpcExpSignalPr, track.tofNSigmaPr ());
1085- registry.fill (HIST (" QA_PID/" )+ HIST (FillTimeName[ft]) + HIST (" ExpSigma_ptwise_protons" ), track.pt (), track.tpcExpSigmaPr (), track.tofNSigmaPr ());
1083+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" TpcdEdx_ptwise_protons" ), track.pt (), track.tpcSignal (), track.tofNSigmaPr ());
1084+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" ExpTpcdEdx_ptwise_protons" ), track.pt (), tpcExpSignalPr, track.tofNSigmaPr ());
1085+ registry.fill (HIST (" QA_PID/" ) + HIST (FillTimeName[ft]) + HIST (" ExpSigma_ptwise_protons" ), track.pt (), track.tpcExpSigmaPr (), track.tofNSigmaPr ());
10861086 }
10871087 }
10881088 }
0 commit comments