@@ -466,6 +466,13 @@ struct lambdapolsp {
466466 return false ;
467467 }
468468
469+ if (pid == 0 && (TMath::Abs (v0.dcapostopv ()) < cMinV0DCAPr || TMath::Abs (v0.dcanegtopv ()) < cMinV0DCAPi)) {
470+ return false ;
471+ }
472+ if (pid == 1 && (TMath::Abs (v0.dcapostopv ()) < cMinV0DCAPi || TMath::Abs (v0.dcanegtopv ()) < cMinV0DCAPr)) {
473+ return false ;
474+ }
475+
469476 // if we made it this far, it's good
470477 return true ;
471478 }
@@ -919,14 +926,30 @@ struct lambdapolsp {
919926
920927 // ___________________________________________________________________________________________________
921928 // retrieve further info provided by StraZDCSP
922- /* auto qxZDCA = collision.qxZDCA();
929+ auto qxZDCA = collision.qxZDCA ();
923930 auto qxZDCC = collision.qxZDCC ();
924931 auto qyZDCA = collision.qyZDCA ();
925- auto qyZDCC = collision.qyZDCC();*/
932+ auto qyZDCC = collision.qyZDCC ();
926933 auto psiZDCC = collision.psiZDCC ();
927934 auto psiZDCA = collision.psiZDCA ();
928- // auto psiZDC=TMath::ATan2((qyZDCC-qyZDCA),(qxZDCC-qxZDCA)); //full event plane
929- auto psiZDC = psiZDCC - psiZDCA;
935+ double modqxZDCA;
936+ double modqyZDCA;
937+ double modqxZDCC;
938+ double modqyZDCC;
939+
940+ if (cqvas) {
941+ modqxZDCA = TMath::Sqrt ((qxZDCA * qxZDCA) + (qyZDCA * qyZDCA)) * TMath::Cos (psiZDCA);
942+ modqyZDCA = TMath::Sqrt ((qxZDCA * qxZDCA) + (qyZDCA * qyZDCA)) * TMath::Sin (psiZDCA);
943+ modqxZDCC = TMath::Sqrt ((qxZDCC * qxZDCC) + (qyZDCC * qyZDCC)) * TMath::Cos (psiZDCC);
944+ modqyZDCC = TMath::Sqrt ((qxZDCC * qxZDCC) + (qyZDCC * qyZDCC)) * TMath::Sin (psiZDCC);
945+ } else {
946+ modqxZDCA = qxZDCA;
947+ modqyZDCA = qyZDCA;
948+ modqxZDCC = qxZDCC;
949+ modqyZDCC = qyZDCC;
950+ }
951+
952+ auto psiZDC = TMath::ATan2 ((modqyZDCC - modqyZDCA), (modqxZDCC - modqxZDCA)); // full event plane
930953
931954 // fill histograms
932955 histos.fill (HIST (" hCentrality" ), centrality);
@@ -970,6 +993,9 @@ struct lambdapolsp {
970993 continue ;
971994 }
972995
996+ if (TMath::Abs (v0.eta ()) > 0.8 )
997+ continue ;
998+
973999 int taga = LambdaTag;
9741000 int tagb = aLambdaTag;
9751001
0 commit comments