@@ -308,7 +308,7 @@ struct decay3bodyBuilder {
308308 Configurable<bool > applySVertexerV0Cuts{" kfparticleConfigurations.applySVertexerV0Cuts" , false , " Apply virtual V0 cuts applied in SVertexer in case of proton mixing" };
309309 ConfigurableAxis bins3BodyRadius{" kfparticleConfigurations.bins3BodyRadius" , {VARIABLE_WIDTH , 0 .0f , 0 .5f , 1 .0f , 1 .5f , 2 .0f , 3 .0f , 4 .0f , 6 .0f , 8 .0f , 10 .0f , 12 .0f , 14 .0f , 16 .0f , 18 .0f , 20 .0f , 30 .0f , 1000.0 }, " Mixing bins - 3body radius" };
310310 // ConfigurableAxis bins3BodyPhi{"kfparticleConfigurations.bins3BodyPhi", {VARIABLE_WIDTH, -180.0f*TMath::Pi()/180, -170.0f*TMath::Pi()/180, -160.0f*TMath::Pi()/180, -150.0f*TMath::Pi()/180, -140.0f*TMath::Pi()/180, -130.0f*TMath::Pi()/180, -120.0f*TMath::Pi()/180, -110.0f*TMath::Pi()/180, -100.0f*TMath::Pi()/180, -90.0f*TMath::Pi()/180, -80.0f*TMath::Pi()/180, -70.0f*TMath::Pi()/180, -60.0f*TMath::Pi()/180, -50.0f*TMath::Pi()/180, -40.0f*TMath::Pi()/180, -30.0f*TMath::Pi()/180, -20.0f*TMath::Pi()/180, -10.0f*TMath::Pi()/180, 0.0f, 10.0f*TMath::Pi()/180, 20.0f*TMath::Pi()/180, 30.0f*TMath::Pi()/180, 40.0f*TMath::Pi()/180, 50.0f*TMath::Pi()/180, 60.0f*TMath::Pi()/180, 70.0f*TMath::Pi()/180, 80.0f*TMath::Pi()/180, 90.0f*TMath::Pi()/180, 100.0f*TMath::Pi()/180, 110.0f*TMath::Pi()/180, 120.0f*TMath::Pi()/180, 130.0f*TMath::Pi()/180, 140.0f*TMath::Pi()/180, 150.0f*TMath::Pi()/180, 160.0f*TMath::Pi()/180, 170.0f*TMath::Pi()/180, 180.0f*TMath::Pi()/180}, "Mixing bins - 3body phi"};
311- ConfigurableAxis bins3BodyPhi{" kfparticleConfigurations.bins3BodyPhi" , {VARIABLE_WIDTH , -180 .0f * TMath::Pi ()/ 180 , -160 .0f * TMath::Pi ()/ 180 , -140 .0f * TMath::Pi ()/ 180 , -120 .0f * TMath::Pi ()/ 180 , -100 .0f * TMath::Pi ()/ 180 , -80 .0f * TMath::Pi ()/ 180 , -60 .0f * TMath::Pi ()/ 180 , -40 .0f * TMath::Pi ()/ 180 , -20 .0f * TMath::Pi ()/ 180 , 0 .0f , 20 .0f * TMath::Pi ()/ 180 , 40 .0f * TMath::Pi ()/ 180 , 60 .0f * TMath::Pi ()/ 180 , 80 .0f * TMath::Pi ()/ 180 , 100 .0f * TMath::Pi ()/ 180 , 120 .0f * TMath::Pi ()/ 180 , 140 .0f * TMath::Pi ()/ 180 , 160 .0f * TMath::Pi ()/ 180 , 180 .0f * TMath::Pi ()/ 180 }, " Mixing bins - 3body phi" };
311+ ConfigurableAxis bins3BodyPhi{" kfparticleConfigurations.bins3BodyPhi" , {VARIABLE_WIDTH , -180 .0f * TMath::Pi () / 180 , -160 .0f * TMath::Pi () / 180 , -140 .0f * TMath::Pi () / 180 , -120 .0f * TMath::Pi () / 180 , -100 .0f * TMath::Pi () / 180 , -80 .0f * TMath::Pi () / 180 , -60 .0f * TMath::Pi () / 180 , -40 .0f * TMath::Pi () / 180 , -20 .0f * TMath::Pi () / 180 , 0 .0f , 20 .0f * TMath::Pi () / 180 , 40 .0f * TMath::Pi () / 180 , 60 .0f * TMath::Pi () / 180 , 80 .0f * TMath::Pi () / 180 , 100 .0f * TMath::Pi () / 180 , 120 .0f * TMath::Pi () / 180 , 140 .0f * TMath::Pi () / 180 , 160 .0f * TMath::Pi () / 180 , 180 .0f * TMath::Pi () / 180 }, " Mixing bins - 3body phi" };
312312 ConfigurableAxis bins3BodyPosZ{" kfparticleConfigurations.bins3BodyPosZ" , {VARIABLE_WIDTH , -300 .0f , -42 .0f , -13 .0f , -6 .0f , -4 .0f , -2 .0f , 0 .0f , 2 .0f , 4 .0f , 6 .0f , 13 .0f , 42 .0f , 300 .0f }, " Mixing bins - 3body z position" };
313313 Configurable<bool > selectVtxZ3bodyMixing{" kfparticleConfigurations.selectVtxZ3bodyMixing" , true , " Select same VtxZ events in case of 3body mixing" };
314314 Configurable<float > VtxZBin3bodyMixing{" kfparticleConfigurations.VtxZBin3bodyMixing" , 1 ., " Bin width for event vtx z position in case of 3body mixing" };
@@ -1497,15 +1497,15 @@ struct decay3bodyBuilder {
14971497 // apply virtual V0 cuts used in SVertexer in case of 3body mixing with proton track
14981498 if (kfparticleConfigurations.mixingType == 1 && kfparticleConfigurations.applySVertexerV0Cuts ) {
14991499 // V0 radius
1500- if (std::sqrt (KFV0 .GetX ()* KFV0 .GetX () + KFV0 .GetY ()* KFV0 .GetY ()) <= 0.5 ) {
1500+ if (std::sqrt (KFV0 .GetX () * KFV0 .GetX () + KFV0 .GetY () * KFV0 .GetY ()) <= 0.5 ) {
15011501 return ;
15021502 }
15031503 // pT
15041504 if (KFV0 .GetPt () <= 0.01 ) {
15051505 return ;
15061506 }
15071507 // pz/pT
1508- if (KFV0 .GetPz ()/ KFV0 .GetPt () >= 2 ) {
1508+ if (KFV0 .GetPz () / KFV0 .GetPt () >= 2 ) {
15091509 return ;
15101510 }
15111511 // cos(PA)
@@ -2205,9 +2205,9 @@ struct decay3bodyBuilder {
22052205 // Function to find bin index (returns -1 if out of range)
22062206 auto findBin = [](float value, const std::vector<float >& binEdges) -> int {
22072207 for (size_t i = 0 ; i < binEdges.size () - 1 ; ++i) {
2208- if (value > binEdges[i] && value <= binEdges[i + 1 ]) {
2209- return i;
2210- }
2208+ if (value > binEdges[i] && value <= binEdges[i + 1 ]) {
2209+ return i;
2210+ }
22112211 }
22122212 return -1 ; // Out of range
22132213 };
@@ -2228,8 +2228,8 @@ struct decay3bodyBuilder {
22282228
22292229 // Determine bin indices
22302230 int radiusBin = findBin (radius, {0 .0f , 0 .5f , 1 .0f , 1 .5f , 2 .0f , 3 .0f , 4 .0f , 6 .0f , 8 .0f , 10 .0f , 12 .0f , 14 .0f , 16 .0f , 18 .0f , 20 .0f , 30 .0f , 1000.0 });
2231- int phiBin = findBin (phi, {-180 .0f * TMath::Pi ()/ 180 , -160 .0f * TMath::Pi ()/ 180 , -140 .0f * TMath::Pi ()/ 180 , -120 .0f * TMath::Pi ()/ 180 , -100 .0f * TMath::Pi ()/ 180 , -80 .0f * TMath::Pi ()/ 180 , -60 .0f * TMath::Pi ()/ 180 , -40 .0f * TMath::Pi ()/ 180 , -20 .0f * TMath::Pi ()/ 180 , 0 .0f , 20 .0f * TMath::Pi ()/ 180 , 40 .0f * TMath::Pi ()/ 180 , 60 .0f * TMath::Pi ()/ 180 , 80 .0f * TMath::Pi ()/ 180 , 100 .0f * TMath::Pi ()/ 180 , 120 .0f * TMath::Pi ()/ 180 , 140 .0f * TMath::Pi ()/ 180 , 160 .0f * TMath::Pi ()/ 180 , 180 .0f * TMath::Pi ()/ 180 });
2232- if (radiusBin >= 0 && phiBin >= 0 ) { // && posZBin >= 0) {
2231+ int phiBin = findBin (phi, {-180 .0f * TMath::Pi () / 180 , -160 .0f * TMath::Pi () / 180 , -140 .0f * TMath::Pi () / 180 , -120 .0f * TMath::Pi () / 180 , -100 .0f * TMath::Pi () / 180 , -80 .0f * TMath::Pi () / 180 , -60 .0f * TMath::Pi () / 180 , -40 .0f * TMath::Pi () / 180 , -20 .0f * TMath::Pi () / 180 , 0 .0f , 20 .0f * TMath::Pi () / 180 , 40 .0f * TMath::Pi () / 180 , 60 .0f * TMath::Pi () / 180 , 80 .0f * TMath::Pi () / 180 , 100 .0f * TMath::Pi () / 180 , 120 .0f * TMath::Pi () / 180 , 140 .0f * TMath::Pi () / 180 , 160 .0f * TMath::Pi () / 180 , 180 .0f * TMath::Pi () / 180 });
2232+ if (radiusBin >= 0 && phiBin >= 0 ) { // && posZBin >= 0) {
22332233 bin3bodyCounts[radiusBin][phiBin]++; // [posZBin]++;
22342234 }
22352235 }
@@ -2285,7 +2285,7 @@ struct decay3bodyBuilder {
22852285 continue ;
22862286 }
22872287 registry.fill (HIST (" QA/EM/h3bodyCombinationCounter" ), 3.5 );
2288-
2288+
22892289 // ---------- do candidate analysis ----------
22902290 bool isMatter1 = false ;
22912291 if (trackBach1.sign () > 0 ) {
@@ -2365,7 +2365,7 @@ struct kfdecay3bodyDataLinkBuilder {
23652365 kfvtxdataLink (lIndices[ii]);
23662366 }
23672367 }
2368-
2368+
23692369 void processStandard (aod::Decay3Bodys const & decay3bodytable, aod::KFVtx3BodyDatas const & vtxdatatable)
23702370 {
23712371 buildDataLink (decay3bodytable, vtxdatatable); // build Decay3Body -> KFDecay3BodyData link table
0 commit comments