5656#include < TPDGCode.h>
5757#include < TF1.h>
5858
59- namespace o2 ::aod
60- {
61-
62- namespace gen_ebyecolltable
63- {
64- DECLARE_SOA_COLUMN (CentralityGen, centralityGen, float );
65- DECLARE_SOA_COLUMN (NetProtNoGen, netProtNoGen, float ); // ! net proton no. in an event
66- DECLARE_SOA_COLUMN (ProtNoGen, protNoGen, float ); // ! proton no. in an event
67- DECLARE_SOA_COLUMN (AntiProtNoGen, antiProtNoGen, float ); // ! antiproton no. in an event
68- } // namespace gen_ebyecolltable
69-
70- DECLARE_SOA_TABLE (ProtGenCollEbyeTables, " AOD" , " PROTGENCOLLEBYETABLE" ,
71- gen_ebyecolltable::CentralityGen,
72- gen_ebyecolltable::NetProtNoGen,
73- gen_ebyecolltable::ProtNoGen,
74- gen_ebyecolltable::AntiProtNoGen);
75- using ProtGenCollEbyeTable = ProtGenCollEbyeTables::iterator;
76-
77- namespace rec_ebyecolltable
78- {
79- DECLARE_SOA_COLUMN (CentralityRec, centralityRec, float );
80- DECLARE_SOA_COLUMN (NetProtNoRec, netProtNoRec, float ); // ! net proton no. in an event
81- DECLARE_SOA_COLUMN (ProtNoRec, protNoRec, float ); // ! proton no. in an event
82- DECLARE_SOA_COLUMN (AntiProtNoRec, antiProtNoRec, float ); // ! antiproton no. in an event
83- } // namespace rec_ebyecolltable
84-
85- DECLARE_SOA_TABLE (ProtRecCollEbyeTables, " AOD" , " PROTRECCOLLEBYETABLE" ,
86- rec_ebyecolltable::CentralityRec,
87- rec_ebyecolltable::NetProtNoRec,
88- rec_ebyecolltable::ProtNoRec,
89- rec_ebyecolltable::AntiProtNoRec);
90- using ProtRecCollEbyeTable = ProtRecCollEbyeTables::iterator;
91-
92- DECLARE_SOA_TABLE (RecCollTables, " AOD" , " RECCOLLTABLE" ,
93- o2::soa::Index<>,
94- rec_ebyecolltable::CentralityRec);
95- using RecCollTable = RecCollTables::iterator;
96-
97- namespace rec_ebyetracktable
98- {
99- DECLARE_SOA_INDEX_COLUMN (RecCollTable, recCollTable);
100- DECLARE_SOA_COLUMN (Pt, pt, float );
101- DECLARE_SOA_COLUMN (Eta, eta, float );
102- DECLARE_SOA_COLUMN (Charge, charge, int );
103- } // namespace rec_ebyetracktable
104-
105- DECLARE_SOA_TABLE (ProtRecCompleteEbyeTables, " AOD" , " PROTRECCOMPLETEEBYETABLE" ,
106- o2::soa::Index<>,
107- rec_ebyetracktable::RecCollTableId,
108- rec_ebyetracktable::Pt,
109- rec_ebyetracktable::Eta,
110- rec_ebyetracktable::Charge);
111- using ProtRecCompleteEbyeTable = ProtRecCompleteEbyeTables::iterator;
112-
113- } // namespace o2::aod
114-
11559using namespace o2 ;
11660using namespace o2 ::framework;
11761using namespace o2 ::framework::expressions;
@@ -991,8 +935,6 @@ struct NetprotonCumulantsMc {
991935 }
992936 }
993937
994- Produces<aod::ProtGenCollEbyeTables> genEbyeCollisions; // ! MC Gen table creation
995-
996938 void processMCGen (aod::McCollision const & mcCollision, aod::McParticles const & mcParticles, const soa::SmallGroups<EventCandidatesMC>& collisions)
997939 {
998940 histos.fill (HIST (" hMC" ), 0.5 );
@@ -1085,7 +1027,6 @@ struct NetprotonCumulantsMc {
10851027 histos.fill (HIST (" hgenProfileTotalProton" ), cent, (nProt + nAntiprot));
10861028 histos.fill (HIST (" hgenProfileProton" ), cent, nProt);
10871029 histos.fill (HIST (" hgenProfileAntiproton" ), cent, nAntiprot);
1088- genEbyeCollisions (cent, netProt, nProt, nAntiprot);
10891030
10901031 // Profiles for generated level cumulants
10911032 // -------------------------------------------------------------------------------------------
@@ -1119,10 +1060,6 @@ struct NetprotonCumulantsMc {
11191060 }
11201061 PROCESS_SWITCH (NetprotonCumulantsMc, processMCGen, " Process Generated" , true );
11211062
1122- Produces<aod::ProtRecCollEbyeTables> recEbyeCollisions; // ! MC Rec table creation
1123- Produces<aod::RecCollTables> recCollisions; // ! MC Rec table creation
1124- Produces<aod::ProtRecCompleteEbyeTables> recEbyeCompleteCollisions; // ! MC Rec table creation with tracks
1125-
11261063 void processMCRec (MyMCRecCollision const & collision, MyMCTracks const & tracks, aod::McCollisions const &, aod::McParticles const &)
11271064 {
11281065 if (!collision.has_mcCollision ()) {
@@ -1147,7 +1084,6 @@ struct NetprotonCumulantsMc {
11471084 histos.fill (HIST (" hCentrec" ), cent);
11481085 histos.fill (HIST (" hMC" ), 5.5 );
11491086 histos.fill (HIST (" hZvtx_after_sel" ), collision.posZ ());
1150- recCollisions (cent);
11511087
11521088 float nProt = 0.0 ;
11531089 float nAntiprot = 0.0 ;
@@ -1207,7 +1143,6 @@ struct NetprotonCumulantsMc {
12071143 trackSelected = selectionPIDold (track);
12081144
12091145 if (trackSelected) {
1210- recEbyeCompleteCollisions (recCollisions.lastIndex (), particle.pt (), particle.eta (), track.sign ());
12111146 // filling nSigma distribution
12121147 histos.fill (HIST (" h2DnsigmaTpcVsPt" ), track.pt (), track.tpcNSigmaPr ());
12131148 histos.fill (HIST (" h2DnsigmaTofVsPt" ), track.pt (), track.tofNSigmaPr ());
@@ -1271,7 +1206,6 @@ struct NetprotonCumulantsMc {
12711206 histos.fill (HIST (" hCorrProfileTotalProton" ), cent, (powerEffProt[1 ] + powerEffAntiprot[1 ]));
12721207 histos.fill (HIST (" hCorrProfileProton" ), cent, powerEffProt[1 ]);
12731208 histos.fill (HIST (" hCorrProfileAntiproton" ), cent, powerEffAntiprot[1 ]);
1274- recEbyeCollisions (cent, netProt, nProt, nAntiprot);
12751209
12761210 // Calculating q_{r,s} as required
12771211 for (int i = 1 ; i < 7 ; i++) {
@@ -2091,7 +2025,6 @@ struct NetprotonCumulantsMc {
20912025 // variables
20922026 auto cent = coll.centFT0M ();
20932027 histos.fill (HIST (" hCentrec" ), cent);
2094- recCollisions (cent);
20952028
20962029 float nProt = 0.0 ;
20972030 float nAntiprot = 0.0 ;
@@ -2143,7 +2076,6 @@ struct NetprotonCumulantsMc {
21432076 trackSelected = selectionPIDold (track);
21442077
21452078 if (trackSelected) {
2146- recEbyeCompleteCollisions (recCollisions.lastIndex (), track.pt (), track.eta (), track.sign ());
21472079 // filling nSigma distribution
21482080 histos.fill (HIST (" h2DnsigmaTpcVsPt" ), track.pt (), track.tpcNSigmaPr ());
21492081 histos.fill (HIST (" h2DnsigmaTofVsPt" ), track.pt (), track.tofNSigmaPr ());
@@ -2200,7 +2132,6 @@ struct NetprotonCumulantsMc {
22002132 histos.fill (HIST (" hCorrProfileTotalProton" ), cent, (powerEffProt[1 ] + powerEffAntiprot[1 ]));
22012133 histos.fill (HIST (" hCorrProfileProton" ), cent, powerEffProt[1 ]);
22022134 histos.fill (HIST (" hCorrProfileAntiproton" ), cent, powerEffAntiprot[1 ]);
2203- recEbyeCollisions (cent, netProt, nProt, nAntiprot);
22042135
22052136 // Calculating q_{r,s} as required
22062137 for (int i = 1 ; i < 7 ; i++) {
0 commit comments