@@ -60,7 +60,7 @@ namespace o2::aod
6060
6161namespace gen_ebyecolltable
6262{
63- DECLARE_SOA_COLUMN (CentralityGen, centralityGen, uint8_t );
63+ DECLARE_SOA_COLUMN (CentralityGen, centralityGen, float );
6464DECLARE_SOA_COLUMN (NetProtNoGen, netProtNoGen, float ); // ! net proton no. in an event
6565DECLARE_SOA_COLUMN (ProtNoGen, protNoGen, float ); // ! proton no. in an event
6666DECLARE_SOA_COLUMN (AntiProtNoGen, antiProtNoGen, float ); // ! antiproton no. in an event
@@ -75,31 +75,35 @@ using ProtGenCollEbyeTable = ProtGenCollEbyeTables::iterator;
7575
7676namespace rec_ebyecolltable
7777{
78- DECLARE_SOA_COLUMN (CentralityRec, centralityRec, uint8_t );
78+ DECLARE_SOA_COLUMN (CentralityRec, centralityRec, float );
7979DECLARE_SOA_COLUMN (NetProtNoRec, netProtNoRec, float ); // ! net proton no. in an event
8080DECLARE_SOA_COLUMN (ProtNoRec, protNoRec, float ); // ! proton no. in an event
8181DECLARE_SOA_COLUMN (AntiProtNoRec, antiProtNoRec, float ); // ! antiproton no. in an event
8282} // namespace rec_ebyecolltable
8383
8484DECLARE_SOA_TABLE (ProtRecCollEbyeTables, " AOD" , " PROTRECCOLLEBYETABLE" ,
85- o2::soa::Index<>,
8685 rec_ebyecolltable::CentralityRec,
8786 rec_ebyecolltable::NetProtNoRec,
8887 rec_ebyecolltable::ProtNoRec,
8988 rec_ebyecolltable::AntiProtNoRec);
9089using ProtRecCollEbyeTable = ProtRecCollEbyeTables::iterator;
9190
91+ DECLARE_SOA_TABLE (ProtRecCollTables, " AOD" , " PROTRECCOLLTABLE" ,
92+ o2::soa::Index<>,
93+ rec_ebyecolltable::CentralityRec);
94+ using ProtRecCollTable = ProtRecCollTables::iterator;
95+
9296namespace rec_ebyetracktable
9397{
94- DECLARE_SOA_INDEX_COLUMN (ProtRecCollEbyeTable, protRecCollEbyeTable );
98+ DECLARE_SOA_INDEX_COLUMN (ProtRecCollTable, protRecCollTable );
9599DECLARE_SOA_COLUMN (Pt, pt, float );
96100DECLARE_SOA_COLUMN (Eta, eta, float );
97101DECLARE_SOA_COLUMN (Charge, charge, int );
98102} // namespace rec_ebyetracktable
99103
100104DECLARE_SOA_TABLE (ProtRecCompleteEbyeTables, " AOD" , " PROTRECCOMPLETEEBYETABLE" ,
101105 o2::soa::Index<>,
102- rec_ebyetracktable::ProtRecCollEbyeTableId ,
106+ rec_ebyetracktable::ProtRecCollTableId ,
103107 rec_ebyetracktable::Pt,
104108 rec_ebyetracktable::Eta,
105109 rec_ebyetracktable::Charge);
@@ -370,6 +374,7 @@ struct NetprotonCumulantsMc {
370374 PROCESS_SWITCH (NetprotonCumulantsMc, processMCGen, " Process Generated" , true );
371375
372376 Produces<aod::ProtRecCollEbyeTables> recEbyeCollisions; // ! MC Rec table creation
377+ Produces<aod::ProtRecCollTables> recCollisions; // ! MC Rec table creation
373378 Produces<aod::ProtRecCompleteEbyeTables> recEbyeCompleteCollisions; // ! MC Rec table creation with tracks
374379
375380 void processMCRec (MyMCRecCollision const & collision, MyMCTracks const & tracks, aod::McCollisions const &, aod::McParticles const &)
@@ -381,6 +386,7 @@ struct NetprotonCumulantsMc {
381386 histos.fill (HIST (" Centrec" ), cent);
382387 histos.fill (HIST (" hMC" ), 5.5 );
383388 histos.fill (HIST (" hZvtx_after_sel" ), collision.posZ ());
389+ recCollisions (cent);
384390
385391 float nProt = 0.0 ;
386392 float nAntiprot = 0.0 ;
@@ -409,7 +415,7 @@ struct NetprotonCumulantsMc {
409415 trackSelected = selectionPIDnew (track);
410416
411417 if (trackSelected) {
412- recEbyeCompleteCollisions (recEbyeCollisions .lastIndex (), particle.pt (), particle.eta (), track.sign ());
418+ recEbyeCompleteCollisions (recCollisions .lastIndex (), particle.pt (), particle.eta (), track.sign ());
413419 if (track.sign () > 0 ) {
414420 histos.fill (HIST (" hrecPtProton" ), particle.pt ()); // ! hist for p rec
415421 histos.fill (HIST (" hrecPtDistProtonVsCentrality" ), particle.pt (), cent);
0 commit comments