@@ -1721,7 +1721,7 @@ struct Derivedupcanalysis {
17211721
17221722 int selGapSide = collision.isUPC () ? getGapSide (collision, false ) : -1 ;
17231723 for (const auto & neutron : neutrons) {
1724- if (!neutron.has_straMCCollision () || !collision.has_straMCCollision ())
1724+ if (!neutron.has_straMCCollision () || !collision.has_straMCCollision ())
17251725 continue ;
17261726
17271727 const auto & mcCollisionNeutron = neutron.straMCCollision_as <StraMCCollisionsFull>(); // take gen. collision associated to the neutron
@@ -1733,14 +1733,18 @@ struct Derivedupcanalysis {
17331733 const float eta = neutron.eta ();
17341734 switch (selGapSide) {
17351735 case 0 : // SGA
1736- if (eta > 8.8 ) selGapSide = -1 ;
1736+ if (eta > 8.8 )
1737+ selGapSide = -1 ;
17371738 break ;
17381739 case 1 : // SGC
1739- if (eta < -8.8 ) selGapSide = -1 ;
1740+ if (eta < -8.8 )
1741+ selGapSide = -1 ;
17401742 break ;
17411743 case 2 : // DG
1742- if (eta > 8.8 ) selGapSide = 1 ;
1743- else if (eta < -8.8 ) selGapSide = 0 ;
1744+ if (eta > 8.8 )
1745+ selGapSide = 1 ;
1746+ else if (eta < -8.8 )
1747+ selGapSide = 0 ;
17441748 break ;
17451749 }
17461750 }
@@ -1759,7 +1763,7 @@ struct Derivedupcanalysis {
17591763 return listBestCollisionIds;
17601764 }
17611765
1762- void fillGenMCHistogramsQA (StraMCCollisionsFull const & mcCollisions,
1766+ void fillGenMCHistogramsQA (StraMCCollisionsFull const & mcCollisions,
17631767 StraCollisonsFullMC const & collisions,
17641768 NeutronsMC const & neutrons)
17651769 {
@@ -1796,7 +1800,7 @@ struct Derivedupcanalysis {
17961800
17971801 int selGapSide = collision.isUPC () ? getGapSide (collision, false ) : -1 ;
17981802 for (const auto & neutron : neutrons) {
1799- if (!neutron.has_straMCCollision () || !collision.has_straMCCollision ())
1803+ if (!neutron.has_straMCCollision () || !collision.has_straMCCollision ())
18001804 continue ;
18011805
18021806 const auto & mcCollisionNeutron = neutron.straMCCollision_as <StraMCCollisionsFull>(); // take gen. collision associated to the neutron
@@ -1808,14 +1812,18 @@ struct Derivedupcanalysis {
18081812 const float eta = neutron.eta ();
18091813 switch (selGapSide) {
18101814 case 0 : // SGA
1811- if (eta > 8.8 ) selGapSide = -1 ;
1815+ if (eta > 8.8 )
1816+ selGapSide = -1 ;
18121817 break ;
18131818 case 1 : // SGC
1814- if (eta < -8.8 ) selGapSide = -1 ;
1819+ if (eta < -8.8 )
1820+ selGapSide = -1 ;
18151821 break ;
18161822 case 2 : // DG
1817- if (eta > 8.8 ) selGapSide = 1 ;
1818- else if (eta < -8.8 ) selGapSide = 0 ;
1823+ if (eta > 8.8 )
1824+ selGapSide = 1 ;
1825+ else if (eta < -8.8 )
1826+ selGapSide = 0 ;
18191827 break ;
18201828 }
18211829 }
@@ -1934,11 +1942,11 @@ struct Derivedupcanalysis {
19341942
19351943 int selGapSide = collision.isUPC () ? getGapSide (collision, true ) : -1 ;
19361944 for (const auto & neutron : neutrons) {
1937- if (!neutron.has_straMCCollision () || !collision.has_straMCCollision ())
1945+ if (!neutron.has_straMCCollision () || !collision.has_straMCCollision ())
19381946 continue ;
19391947
19401948 const auto & mcCollisionNeutron = neutron.straMCCollision_as <StraMCCollisionsFull>(); // take gen. collision associated to the neutron
1941- const auto & mcCollision = collision.straMCCollision_as <StraMCCollisionsFull>(); // take gen. collision associated to the rec. collision
1949+ const auto & mcCollision = collision.straMCCollision_as <StraMCCollisionsFull>(); // take gen. collision associated to the rec. collision
19421950
19431951 // Consider neutrons from the same collision
19441952 if (mcCollisionNeutron.globalIndex () != mcCollision.globalIndex ())
@@ -1947,14 +1955,18 @@ struct Derivedupcanalysis {
19471955 const float eta = neutron.eta ();
19481956 switch (selGapSide) {
19491957 case 0 : // SGA
1950- if (eta > 8.8 ) selGapSide = -1 ;
1958+ if (eta > 8.8 )
1959+ selGapSide = -1 ;
19511960 break ;
19521961 case 1 : // SGC
1953- if (eta < -8.8 ) selGapSide = -1 ;
1962+ if (eta < -8.8 )
1963+ selGapSide = -1 ;
19541964 break ;
19551965 case 2 : // DG
1956- if (eta > 8.8 ) selGapSide = 1 ;
1957- else if (eta < -8.8 ) selGapSide = 0 ;
1966+ if (eta > 8.8 )
1967+ selGapSide = 1 ;
1968+ else if (eta < -8.8 )
1969+ selGapSide = 0 ;
19581970 break ;
19591971 }
19601972 }
@@ -2047,11 +2059,11 @@ struct Derivedupcanalysis {
20472059
20482060 int selGapSide = collision.isUPC () ? getGapSide (collision, true ) : -1 ;
20492061 for (const auto & neutron : neutrons) {
2050- if (!neutron.has_straMCCollision () || !collision.has_straMCCollision ())
2062+ if (!neutron.has_straMCCollision () || !collision.has_straMCCollision ())
20512063 continue ;
2052-
2064+
20532065 const auto & mcCollisionNeutron = neutron.straMCCollision_as <StraMCCollisionsFull>(); // take gen. collision associated to the neutron
2054- const auto & mcCollision = collision.straMCCollision_as <StraMCCollisionsFull>(); // take gen. collision associated to the rec. collision
2066+ const auto & mcCollision = collision.straMCCollision_as <StraMCCollisionsFull>(); // take gen. collision associated to the rec. collision
20552067
20562068 // Consider neutrons from the same collision
20572069 if (mcCollisionNeutron.globalIndex () != mcCollision.globalIndex ())
@@ -2060,14 +2072,18 @@ struct Derivedupcanalysis {
20602072 const float eta = neutron.eta ();
20612073 switch (selGapSide) {
20622074 case 0 : // SGA
2063- if (eta > 8.8 ) selGapSide = -1 ;
2075+ if (eta > 8.8 )
2076+ selGapSide = -1 ;
20642077 break ;
20652078 case 1 : // SGC
2066- if (eta < -8.8 ) selGapSide = -1 ;
2079+ if (eta < -8.8 )
2080+ selGapSide = -1 ;
20672081 break ;
20682082 case 2 : // DG
2069- if (eta > 8.8 ) selGapSide = 1 ;
2070- else if (eta < -8.8 ) selGapSide = 0 ;
2083+ if (eta > 8.8 )
2084+ selGapSide = 1 ;
2085+ else if (eta < -8.8 )
2086+ selGapSide = 0 ;
20712087 break ;
20722088 }
20732089 }
@@ -2162,13 +2178,13 @@ struct Derivedupcanalysis {
21622178 const bool isNegativePion = (pdgNeg == PDG_t::kPiMinus ) || (doTreatPiToMuon && pdgNeg == PDG_t::kMuonMinus );
21632179
21642180 // Fill histograms
2165- if ( (pdgV0 == PDG_t::kK0Short ) && isPositivePion && isNegativePion) {
2181+ if ((pdgV0 == PDG_t::kK0Short ) && isPositivePion && isNegativePion) {
21662182 histos.fill (HIST (kParticlenames [0 ]) + HIST (" /mc/h7dGen" ), ft0ampl, nTracksGlobal, mcCollision.multMCNParticlesEta08 (), pTmc, static_cast <int >(upcCuts.genGapSide ), ymc, mcCollision.generatorsID ());
21672183 }
2168- if ( (pdgV0 == PDG_t::kLambda0 ) && isPositiveProton && isNegativePion) {
2184+ if ((pdgV0 == PDG_t::kLambda0 ) && isPositiveProton && isNegativePion) {
21692185 histos.fill (HIST (kParticlenames [1 ]) + HIST (" /mc/h7dGen" ), ft0ampl, nTracksGlobal, mcCollision.multMCNParticlesEta08 (), pTmc, static_cast <int >(upcCuts.genGapSide ), ymc, mcCollision.generatorsID ());
21702186 }
2171- if ( (pdgV0 == PDG_t::kLambda0Bar ) && isPositivePion && isNegativeProton) {
2187+ if ((pdgV0 == PDG_t::kLambda0Bar ) && isPositivePion && isNegativeProton) {
21722188 histos.fill (HIST (kParticlenames [2 ]) + HIST (" /mc/h7dGen" ), ft0ampl, nTracksGlobal, mcCollision.multMCNParticlesEta08 (), pTmc, static_cast <int >(upcCuts.genGapSide ), ymc, mcCollision.generatorsID ());
21732189 }
21742190 } // V0 end
@@ -2181,10 +2197,9 @@ struct Derivedupcanalysis {
21812197 // Kinematics (|y| < rapidityCut)
21822198 float pTmc = cascMC.ptMC ();
21832199 float ymc = 1e3 ;
2184- if ( (cascMC.pdgCode () == PDG_t::kXiMinus ) || (cascMC.pdgCode () == PDG_t::kXiPlusBar ) ) {
2200+ if ((cascMC.pdgCode () == PDG_t::kXiMinus ) || (cascMC.pdgCode () == PDG_t::kXiPlusBar )) {
21852201 ymc = cascMC.rapidityMC (0 );
2186- }
2187- else if ( (cascMC.pdgCode () == PDG_t::kOmegaMinus ) || (cascMC.pdgCode () == PDG_t::kOmegaPlusBar ) ) {
2202+ } else if ((cascMC.pdgCode () == PDG_t::kOmegaMinus ) || (cascMC.pdgCode () == PDG_t::kOmegaPlusBar )) {
21882203 ymc = cascMC.rapidityMC (2 );
21892204 }
21902205 if (std::abs (ymc) > rapidityCut)
@@ -2227,16 +2242,16 @@ struct Derivedupcanalysis {
22272242 const bool isBachelorNegativeKaon = (pdgBach == PDG_t::kKMinus );
22282243
22292244 // Fill histograms
2230- if ( (pdgCasc == PDG_t::kXiMinus ) && isPositiveProton && isNegativePion && isBachelorNegativePion) {
2245+ if ((pdgCasc == PDG_t::kXiMinus ) && isPositiveProton && isNegativePion && isBachelorNegativePion) {
22312246 histos.fill (HIST (kParticlenames [3 ]) + HIST (" /mc/h7dGen" ), ft0ampl, nTracksGlobal, mcCollision.multMCNParticlesEta08 (), pTmc, static_cast <int >(upcCuts.genGapSide ), ymc, mcCollision.generatorsID ());
22322247 }
2233- if ( (pdgCasc == PDG_t::kXiPlusBar ) && isNegativeProton && isPositivePion && isBachelorPositivePion) {
2248+ if ((pdgCasc == PDG_t::kXiPlusBar ) && isNegativeProton && isPositivePion && isBachelorPositivePion) {
22342249 histos.fill (HIST (kParticlenames [4 ]) + HIST (" /mc/h7dGen" ), ft0ampl, nTracksGlobal, mcCollision.multMCNParticlesEta08 (), pTmc, static_cast <int >(upcCuts.genGapSide ), ymc, mcCollision.generatorsID ());
22352250 }
2236- if ( (pdgCasc == PDG_t::kOmegaMinus ) && isPositiveProton && isNegativePion && isBachelorNegativeKaon) {
2251+ if ((pdgCasc == PDG_t::kOmegaMinus ) && isPositiveProton && isNegativePion && isBachelorNegativeKaon) {
22372252 histos.fill (HIST (kParticlenames [5 ]) + HIST (" /mc/h7dGen" ), ft0ampl, nTracksGlobal, mcCollision.multMCNParticlesEta08 (), pTmc, static_cast <int >(upcCuts.genGapSide ), ymc, mcCollision.generatorsID ());
22382253 }
2239- if ( (pdgCasc == PDG_t::kOmegaPlusBar ) && isNegativeProton && isPositivePion && isBachelorPositiveKaon) {
2254+ if ((pdgCasc == PDG_t::kOmegaPlusBar ) && isNegativeProton && isPositivePion && isBachelorPositiveKaon) {
22402255 histos.fill (HIST (kParticlenames [6 ]) + HIST (" /mc/h7dGen" ), ft0ampl, nTracksGlobal, mcCollision.multMCNParticlesEta08 (), pTmc, static_cast <int >(upcCuts.genGapSide ), ymc, mcCollision.generatorsID ());
22412256 }
22422257 } // Cascade end
0 commit comments