@@ -1865,66 +1865,66 @@ struct Kstarqa {
18651865=======
18661866 hInvMass.fill (HIST (" CorrFactors/hMultiplicityVsMultMC" ), multiplicity, multiplicityNch);
18671867>>>>>>> c669bd852 (corrected event loss histograms)
1868- hInvMass.fill (HIST (" CorrFactors/hNrecInGen" ), collisions.size ());
1869- hInvMass.fill (HIST (" CorrFactors/MultiplicityGen" ), multiplicityNch);
1870- if (isSelectedEvent) {
1868+ hInvMass.fill (HIST (" CorrFactors/hNrecInGen" ), collisions.size ());
1869+ hInvMass.fill (HIST (" CorrFactors/MultiplicityGen" ), multiplicityNch);
1870+ if (isSelectedEvent) {
18711871<<<<<<< HEAD
1872- hInvMass.fill (HIST (" CorrFactors/MultiplicityRec2D" ), multiplicity, mcCollision.multMCNParticlesEta08 ());
1872+ hInvMass.fill (HIST (" CorrFactors/MultiplicityRec2D" ), multiplicity, mcCollision.multMCNParticlesEta08 ());
18731873>>>>>>> 2ddf06887 (minor fix to code)
18741874=======
18751875 hInvMass.fill (HIST (" CorrFactors/MultiplicityRec" ), multiplicityNch);
18761876>>>>>>> c669bd852 (corrected event loss histograms)
1877- }
1877+ }
18781878
1879- for (const auto & mcParticle : mcParticles) {
1879+ for (const auto & mcParticle : mcParticles) {
18801880
1881- if (std::abs (mcParticle.y ()) >= selectionConfig.rapidityMotherData )
1882- continue ;
1881+ if (std::abs (mcParticle.y ()) >= selectionConfig.rapidityMotherData )
1882+ continue ;
18831883
1884- if (std::abs (mcParticle.pdgCode ()) == o2::constants::physics::kK0Star892 ) {
1884+ if (std::abs (mcParticle.pdgCode ()) == o2::constants::physics::kK0Star892 ) {
18851885
1886- auto kDaughters = mcParticle.daughters_as <aod::McParticles>();
1887- if (kDaughters .size () != selectionConfig.noOfDaughters ) {
1888- continue ;
1889- }
1886+ auto kDaughters = mcParticle.daughters_as <aod::McParticles>();
1887+ if (kDaughters .size () != selectionConfig.noOfDaughters ) {
1888+ continue ;
1889+ }
18901890
1891- bool hasPos = false ;
1892- bool hasNeg = false ;
1891+ bool hasPos = false ;
1892+ bool hasNeg = false ;
18931893
1894- auto passkaon = false ;
1895- auto passpion = false ;
1896- for (const auto & kCurrentDaughter : kDaughters ) {
1897- // if (!kCurrentDaughter.isPhysicalPrimary())
1898- // continue;
1894+ auto passkaon = false ;
1895+ auto passpion = false ;
1896+ for (const auto & kCurrentDaughter : kDaughters ) {
1897+ // if (!kCurrentDaughter.isPhysicalPrimary())
1898+ // continue;
18991899
1900- int pdgDau = kCurrentDaughter .pdgCode ();
1901- int sign = (pdgDau > 0 ) - (pdgDau < 0 );
1900+ int pdgDau = kCurrentDaughter .pdgCode ();
1901+ int sign = (pdgDau > 0 ) - (pdgDau < 0 );
19021902
1903- if (sign > 0 )
1904- hasPos = true ;
1905- if (sign < 0 )
1906- hasNeg = true ;
1903+ if (sign > 0 )
1904+ hasPos = true ;
1905+ if (sign < 0 )
1906+ hasNeg = true ;
19071907
1908- if (std::abs (kCurrentDaughter .pdgCode ()) == PDG_t::kKPlus ) {
1909- passkaon = true ;
1910- daughter1 = ROOT::Math::PxPyPzMVector (kCurrentDaughter .px (), kCurrentDaughter .py (), kCurrentDaughter .pz (), massKa);
1908+ if (std::abs (kCurrentDaughter .pdgCode ()) == PDG_t::kKPlus ) {
1909+ passkaon = true ;
1910+ daughter1 = ROOT::Math::PxPyPzMVector (kCurrentDaughter .px (), kCurrentDaughter .py (), kCurrentDaughter .pz (), massKa);
19111911
1912- } else if (std::abs (kCurrentDaughter .pdgCode ()) == PDG_t::kPiPlus ) {
1913- passpion = true ;
1914- daughter2 = ROOT::Math::PxPyPzMVector (kCurrentDaughter .px (), kCurrentDaughter .py (), kCurrentDaughter .pz (), massPi);
1912+ } else if (std::abs (kCurrentDaughter .pdgCode ()) == PDG_t::kPiPlus ) {
1913+ passpion = true ;
1914+ daughter2 = ROOT::Math::PxPyPzMVector (kCurrentDaughter .px (), kCurrentDaughter .py (), kCurrentDaughter .pz (), massPi);
1915+ }
19151916 }
1916- }
19171917
1918- if ((passkaon && passpion) && (hasPos && hasNeg)) {
1919- mother = daughter1 + daughter2; // Kstar meson
1918+ if ((passkaon && passpion) && (hasPos && hasNeg)) {
1919+ mother = daughter1 + daughter2; // Kstar meson
19201920
1921- hInvMass.fill (HIST (" CorrFactors/h2dGenKstar" ), multiplicity, mother.Pt ());
1921+ hInvMass.fill (HIST (" CorrFactors/h2dGenKstar" ), multiplicity, mother.Pt ());
19221922<<<<<<< HEAD
19231923<<<<<<< HEAD
1924- hInvMass.fill (HIST (" CorrFactors/h3dGenKstarVsMultMCVsMultiplicity" ), multiplicityNch, multiplicity, mother.Pt ());
1925- hInvMass.fill (HIST (" CorrFactors/hSignalLossDenominator" ), mother.pt (), multiplicityNch);
1926- if (isSelectedEvent) {
1927- hInvMass.fill (HIST (" CorrFactors/hSignalLossNumerator" ), mother.pt (), multiplicityNch);
1924+ hInvMass.fill (HIST (" CorrFactors/h3dGenKstarVsMultMCVsMultiplicity" ), multiplicityNch, multiplicity, mother.Pt ());
1925+ hInvMass.fill (HIST (" CorrFactors/hSignalLossDenominator" ), mother.pt (), multiplicityNch);
1926+ if (isSelectedEvent) {
1927+ hInvMass.fill (HIST (" CorrFactors/hSignalLossNumerator" ), mother.pt (), multiplicityNch);
19281928=======
19291929 hInvMass.fill (HIST (" CorrFactors/h3dGenKstarVsMultMCVsMultiplicity" ), mcCollision.multMCNParticlesEta08 (), multiplicity, mother.Pt ());
19301930 hInvMass.fill (HIST (" CorrFactors/hSignalLossDenominator3D" ), mother.pt (), multiplicity, mcCollision.multMCNParticlesEta08 ());
@@ -1933,23 +1933,23 @@ struct Kstarqa {
19331933 hInvMass.fill (HIST (" CorrFactors/hSignalLossNumerator" ), mother.pt (), multiplicity);
19341934>>>>>>> 087e3b03d (different event/signal loss method)
19351935=======
1936- hInvMass.fill (HIST (" CorrFactors/hSignalLossNumerator3D" ), mother.pt (), multiplicity, mcCollision.multMCNParticlesEta08 ());
1936+ hInvMass.fill (HIST (" CorrFactors/hSignalLossNumerator3D" ), mother.pt (), multiplicity, mcCollision.multMCNParticlesEta08 ());
19371937>>>>>>> 2ddf06887 (minor fix to code)
19381938=======
19391939 hInvMass.fill (HIST (" CorrFactors/h3dGenKstarVsMultMCVsMultiplicity" ), multiplicityNch, multiplicity, mother.Pt ());
19401940 hInvMass.fill (HIST (" CorrFactors/hSignalLossDenominator" ), mother.pt (), multiplicityNch);
19411941 if (isSelectedEvent) {
19421942 hInvMass.fill (HIST (" CorrFactors/hSignalLossNumerator" ), mother.pt (), multiplicityNch);
19431943>>>>>>> c669bd852 (corrected event loss histograms)
1944+ }
19441945 }
19451946 }
19461947 }
1947- }
19481948
1949- if (collisions.size () == 0 )
1950- return ;
1949+ if (collisions.size () == 0 )
1950+ return ;
19511951
1952- hInvMass.fill (HIST (" CorrFactors/hGenEvents" ), multiplicityNch, 3.5 );
1952+ hInvMass.fill (HIST (" CorrFactors/hGenEvents" ), multiplicityNch, 3.5 );
19531953 }
19541954 PROCESS_SWITCH (Kstarqa, processEvSigLossFactors, " Process Event and Signal loss" , false );
19551955
0 commit comments