#include "AliKFParticle.h"
#include "AliKFVertex.h"
#include "AliLog.h"
+#include "AliExternalTrackParam.h"
#include "AliHFEcollection.h"
fQA->CreateTH1Fvector1(2, "h_cut_Gamma_DCA", "DCA between the gamma daughters; dca (cm); counts", 100, 0, 2);
fQA->CreateTH1Fvector1(2, "h_cut_Gamma_VtxR_old", "*old* Radius of the gamma conversion vertex; r (cm); counts", 1000, 0, 100);
fQA->CreateTH1Fvector1(2, "h_cut_Gamma_VtxR", "Radius of the gamma conversion vertex; r (cm); counts", 1000, 0, 100);
- fQA->CreateTH1Fvector1(2, "h_cut_Gamma_OA", "opening angle of the gamma products; opening angle (rad); counts", 100, 0, 1);
fQA->CreateTH1Fvector1(2, "h_cut_Gamma_PP", "gamma psi pair angle; psi pairangle (rad); counts", 100, 0, 2);
fQA->CreateTH1Fvector1(2, "h_cut_Gamma_Chi2", "gamma Chi2/NDF; Chi2/NDF; counts", 100, 0, 50);
+ fQA->CreateTH1Fvector1(2, "h_cut_Gamma_Sep", "gamma separation dist at TPC inned wall", 100, 0, 50);
fQA->CreateTH1Fvector1(7, "h_Gamma_Mass", "Invariant mass of gammas; mass (GeV/c^{2}); counts", 100, 0, 0.2);
-
+
// kaons
fQA->CreateTH1Fvector1(2, "h_cut_K0_CosPoint", "K0 Cosine pointing angle; cos point. angle; counts", 100, 0, 0.1);
fQA->CreateTH2F("h_L_checks", "Lambda candidate check[0] -v- check[1]; check[0]; check[1]", 5, -0.75, 1.75, 6, -0.75, 1.75 );
// electrons
- fQA->CreateTH1Fvector1(9, "h_Electron_P", "Momenta of conversion electrons -cuts-; P (GeV/c); counts", 50, 0.1, 20, 0);
+ fQA->CreateTH1Fvector1(7, "h_Electron_P", "Momenta of conversion electrons -cuts-; P (GeV/c); counts", 50, 0.1, 20, 0);
// K0 pions
fQA->CreateTH1Fvector1(8, "h_PionK0_P", "Momenta of K0 pions -cuts-; P (GeV/c) counts;", 50, 0.1, 20, 0);
//
// possibly new cuts
//
-
- // Gamma
- fQA->CreateTH2Fvector1(2, "h_cut_Gamma_OAvP", "open. ang. of the Gamma daughters versus Gamma mom; Gamma p (GeV/c); opening angle (pions) (rad)", 100, 0.1, 10, 200, 0., 0.2);
- // K0
- fQA->CreateTH2Fvector1(2, "h_cut_K0_OAvP", "open. ang. of the K0 daughters versus K0 momentum; K0 p (GeV/c); opening angle (pions) (rad)", 100, 0.1, 10, 100, 0, 3.5);
- // Lambda
- fQA->CreateTH2Fvector1(2, "h_cut_L_OAvP", "open. ang. of the L daughters versus L momentum; Lambda p (GeV/c); openeing angle pion-proton (rad)", 100, 0.1, 10, 100, 0, 3.5);
- fQA->CreateTH2Fvector1(2, "h_cut_L_rdp_v_mp", "relative L daughter mom -v- mother mom; L mom (GeV/c); relative daughter mom p2/p1", 100, 0.1, 10, 100, 0, 1);
-
+ fQA->CreateTH2Fvector1(2, "h_cut_L_rdp_v_mp", "relative L daughter mom -v- mother mom; L mom (GeV/c); relative daughter mom p2/p1", 100, 0.1, 10, 100, 0, 1);
// THnSparse histograms
fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_CosPoint_S", "S - Gamma Cosine pointing angle; mom (GeV/c); cos point. angle", pN, pMin, pMax, 50, 0, 0.1, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_DCA_S", "S - DCA between the gamma daughters; mom (GeV/c); dca (cm)", pN, pMin, pMax, 50, 0, 2, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_VtxR_S", "S - Radius of the gamma conversion vertex; mom (GeV/c); r (cm)", pN, pMin, pMax, 100, 0, 100, 0);
- fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_OA_S", "S - opening angle of the gamma products; mom (GeV/c); opening angle (rad)", pN, pMin, pMax, 50, 0, 0.3, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_PP_S", "S - gamma psi pair angle; mom (GeV/c); psi pairangle (rad)", pN, pMin, pMax, 50, 0, 0.5, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_Chi2_S", "S - gamma Chi2/NDF; mom (GeV/c); Chi2/NDF", pN, pMin, pMax, 50, 0, 100, 0);
+ fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_Sep_S", "S - gamma separation TPC-inner; mom (GeV/c); tracks separatin (cm)", pN, pMin, pMax, 100, 0, 50, 0);
- fQAmc->CreateTH1Fvector1(8, "h_Gamma_Mass_S", "S - Invariant mass of gammas; mass (GeV/c^{2}); counts", 100, 0, 0.2);
+ fQAmc->CreateTH1Fvector1(9, "h_Gamma_Mass_S", "S - Invariant mass of gammas; mass (GeV/c^{2}); counts", 100, 0, 0.2);
// gamma background
fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_CosPoint_B", "B - Gamma Cosine pointing angle; mom (GeV/c); cos point. angle", pN, pMin, pMax, 50, 0, 0.1, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_DCA_B", "B - DCA between the gamma daughters; mom (GeV/c); dca (cm)", pN, pMin, pMax, 50, 0, 2, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_VtxR_B", "B - Radius of the gamma conversion vertex; mom (GeV/c); r (cm)", pN, pMin, pMax, 100, 0, 100, 0);
- fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_OA_B", "B - opening angle of the gamma products; mom (GeV/c); opening angle (rad)", pN, pMin, pMax, 50, 0, 0.3, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_PP_B", "B - gamma psi pair angle; mom (GeV/c); psi pairangle (rad)", pN, pMin, pMax, 50, 0, 0.5, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_Chi2_B", "B - gamma Chi2/NDF; mom (GeV/c); Chi2/NDF", pN, pMin, pMax, 50, 0, 100, 0);
+ fQAmc->CreateTH2Fvector1(2, "h_cut_Gamma_Sep_B", "B - gamma separation TPC-inner; mom (GeV/c); tracks separatin (cm)", pN, pMin, pMax, 100, 0, 50, 0);
- fQAmc->CreateTH1Fvector1(8, "h_Gamma_Mass_B", "B - Invariant mass of gammas; mass (GeV/c^{2}); counts", 100, 0, 0.2);
+ fQAmc->CreateTH1Fvector1(9, "h_Gamma_Mass_B", "B - Invariant mass of gammas; mass (GeV/c^{2}); counts", 100, 0, 0.2);
// kaons signal
fQAmc->CreateTH2Fvector1(2, "h_cut_K0_CosPoint_S", "S - K0 Cosine pointing angle; mom (GeV/c); cos point. angle", pN, pMin, pMax, 50, 0, 0.1, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_K0_DCA_S", "S - DCA between the K0 daughters; mom (GeV/c); dca (cm)", pN, pMin, pMax, 50, 0, 2, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_K0_VtxR_S", "S - Radius of the K0 decay vertex; mom (GeV/c); r (cm)", pN, pMin, pMax, 50, 0, 100, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_K0_Chi2_S", "S - K0 Chi2/NDF; mom (GeV/c); Chi2/NDF", pN, pMin, pMax, 50, 0, 100, 0);
- fQAmc->CreateTH2Fvector1(2, "h_cut_K0_OA_S", "S - opening angle of the K0 pions; mom (GeV/c); opening angle (rad)", pN, pMin, pMax, 100, 0, 1, 0);
fQAmc->CreateTH1Fvector1(5, "h_K0_Mass_S", "S - Invariant mass of K0; mass (GeV/c^{2}); counts", 125, 0.45, 0.55);
// kaons background
fQAmc->CreateTH2Fvector1(2, "h_cut_K0_DCA_B", "B - DCA between the K0 daughters; mom (GeV/c); dca (cm)", pN, pMin, pMax, 50, 0, 2, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_K0_VtxR_B", "B - Radius of the K0 decay vertex; mom (GeV/c); r (cm)", pN, pMin, pMax, 50, 0, 100, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_K0_Chi2_B", "B - K0 Chi2/NDF; mom (GeV/c); Chi2/NDF", pN, pMin, pMax, 50, 0, 100, 0);
- fQAmc->CreateTH2Fvector1(2, "h_cut_K0_OA_B", "B - opening angle of the K0 pions; mom (GeV/c); opening angle (rad)", pN, pMin, pMax, 100, 0, 1, 0);
fQAmc->CreateTH1Fvector1(5, "h_K0_Mass_B", "B - Invariant mass of K0; mass (GeV/c^{2}); counts", 125, 0.45, 0.55);
fQAmc->CreateTH2Fvector1(2, "h_cut_L_DCA_S", "S - DCA between the L daughters; mom (GeV/c); dca (cm)", pN, pMin, pMax, 50, 0, 2, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_L_VtxR_S", "S - Radius of the L decay vertex; mom (GeV/c); r (cm)", pN, pMin, pMax, 50, 0, 100, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_L_Chi2_S", "S - L Chi2/NDF; mom (GeV/c); Chi2/NDF", pN, pMin, pMax, 50, 0, 100, 0);
- fQAmc->CreateTH2Fvector1(2, "h_cut_L_OA_S", "S - opening angle of the L p-p; mom (GeV/c); opening angle (rad)", pN, pMin, pMax, 100, 0, 1, 0);
fQAmc->CreateTH1Fvector1(5, "h_L_Mass_S", "S - Invariant mass of L; mass (GeV/c^{2}); counts", 60, 1.1, 1.13);
fQAmc->CreateTH1Fvector1(5, "h_AL_Mass_S", "S - Invariant mass of anti L; mass (GeV/c^{2}); counts", 60, 1.1, 1.13);
fQAmc->CreateTH2Fvector1(2, "h_cut_L_DCA_B", "B - DCA between the L daughters; mom (GeV/c); dca (cm)", pN, pMin, pMax, 50, 0, 2, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_L_VtxR_B", "B - Radius of the L decay vertex; mom (GeV/c); r (cm)", pN, pMin, pMax, 50, 0, 100, 0);
fQAmc->CreateTH2Fvector1(2, "h_cut_L_Chi2_B", "B - L Chi2/NDF; mom (GeV/c); Chi2/NDF", pN, pMin, pMax, 50, 0, 100, 0);
- fQAmc->CreateTH2Fvector1(2, "h_cut_L_OA_B", "B - opening angle of the L p-p; mom (GeV/c); opening angle (rad)", pN, pMin, pMax, 100, 0, 1, 0);
-
fQAmc->CreateTH2Fvector1(2, "h_cut_L_rdp_v_mp_S", "S - relative L daughter mom -v- mother mom; L mom (GeV/c); relative daughter mom p2/p1", 100, 0.1, 10, 100, 0, 1);
fQAmc->CreateTH2Fvector1(2, "h_cut_L_rdp_v_mp_B", "B - relative L daughter mom -v- mother mom; L mom (GeV/c); relative daughter mom p2/p1", 100, 0.1, 10, 100, 0, 1);
fQAmc->CreateTH1Fvector1(5, "h_LAL_Mass_B", "B - Invariant mass of anti L; mass (GeV/c^{2}); counts", 60, 1.1, 1.13);
// No. of TPC clusters
fQA->Fill("h_ST_NclsTPC", track->GetTPCNcls());
- if(track->GetTPCNcls() < 80) return kFALSE; //
+ if(track->GetTPCNcls() < 1) return kFALSE; //
// TPC refit
if((status & AliESDtrack::kTPCrefit)){
Int_t nTPCclusters = track->GetTPCclusters(0);
Float_t chi2perTPCcluster = track->GetTPCchi2()/Float_t(nTPCclusters);
fQA->Fill("h_ST_chi2TPCcls", chi2perTPCcluster);
- if(chi2perTPCcluster > 3.5) return kFALSE; // 4.0
+ if(chi2perTPCcluster > 4.0) return kFALSE; // 4.0
// TPC cluster ratio
Float_t cRatioTPC = track->GetTPCNclsF() > 0. ? static_cast<Float_t>(track->GetTPCNcls())/static_cast<Float_t> (track->GetTPCNclsF()) : 1.;
// pt
fQA->Fill("h_ST_pt",track->Pt());
- if(track->Pt() < 0.1 || track->Pt() > 100) return kFALSE; //
+ //if(track->Pt() < 0.1 || track->Pt() > 100) return kFALSE; //
// eta
fQA->Fill("h_ST_eta", track->Eta());
Float_t p[2] = {d[0]->GetP(), d[1]->GetP()};
// Cut values
- const Double_t cutChi2NDF = 40.; // ORG [7.]
+ const Double_t cutChi2NDF = 10.; // ORG [7.]
const Double_t cutCosPoint[2] = {0., 0.02}; // ORG [0., 0.03]
const Double_t cutDCA[2] = {0., 0.25}; // ORG [0., 0.25]
- const Double_t cutProdVtxR[2] = {8., 90.}; // ORG [6., 9999]
+ const Double_t cutProdVtxR[2] = {3., 90.}; // ORG [6., 9999]
const Double_t cutPsiPair[2] = {0., 0.05}; // ORG [0. 0.05]
- const Double_t cutOAngle[2] = {0, 0.1}; // ORG [0., 0.1]
// mass cut
const Double_t cutMass = 0.05; // ORG [0.05]
+
+ //
+ // possible new cuts
+ //
+ // separation cut at the entrance to the TPC
+ const Double_t cutSeparation = 0.; // ORG 3.0 cm
+
+
+
// Values
// cos pointing angle
r2 = TMath::Sqrt(xy[0]*xy[0] + xy[1]*xy[1]);
}
-
- // Opening angle
- Double_t oAngle = OpenAngle(v0);
-
// psi pair
Double_t psiPair = PsiPair(v0);
// V0 chi2/ndf
Double_t chi2ndf = kfMother->GetChi2()/kfMother->GetNDF();
-
if(kfMother) delete kfMother;
+
+ // Separation
+ AliExternalTrackParam const *param[2];
+ param[0] = d[0]->GetInnerParam();
+ param[1] = d[1]->GetInnerParam();
+ Double_t sep = 999.;
+ if(param[0] && param[1]){
+ TVector3 xyz[3];
+ xyz[0].SetXYZ(param[0]->GetX(), param[0]->GetY(), param[0]->GetZ());
+ xyz[1].SetXYZ(param[1]->GetX(), param[1]->GetY(), param[1]->GetZ());
+ xyz[2] = xyz[0] - xyz[1];
+ sep = xyz[2].Mag();
+ }
+
+
//
// Apply the cuts, produce QA plots (with mass cut)
fQA->Fill("h_cut_Gamma_DCA", 0, dca);
fQA->Fill("h_cut_Gamma_VtxR_old", 0, r);
fQA->Fill("h_cut_Gamma_VtxR", 0, r2);
- fQA->Fill("h_cut_Gamma_OA", 0, oAngle);
fQA->Fill("h_cut_Gamma_PP", 0, psiPair);
fQA->Fill("h_cut_Gamma_Chi2", 0, chi2ndf);
fQA->Fill("h_cut_Gamma_Chi2", 1, chi2ndf, iP);
- fQA->Fill("h_cut_Gamma_OAvP", 0, iP, oAngle);
+ fQA->Fill("h_cut_Gamma_Sep", 0, iP, sep);
+
if(fMCEvent){
// MC signal
fQAmc->Fill("h_cut_Gamma_CosPoint_S", 0, iP, cosPoint);
fQAmc->Fill("h_cut_Gamma_DCA_S", 0, iP, dca);
fQAmc->Fill("h_cut_Gamma_VtxR_S", 0, iP, r2);
- fQAmc->Fill("h_cut_Gamma_OA_S", 0, iP, oAngle);
fQAmc->Fill("h_cut_Gamma_PP_S", 0, iP, psiPair);
fQAmc->Fill("h_cut_Gamma_Chi2_S", 0, iP, chi2ndf);
fQAmc->Fill("h_cut_Gamma_Chi2_S", 1, iP, chi2ndf);
+ fQAmc->Fill("h_cut_Gamma_Sep_S", 0, iP, sep);
fQAmc->Fill("h_Electron_P_S", 0, p[0]);
fQAmc->Fill("h_Electron_P_S", 0, p[1]);
}
fQAmc->Fill("h_cut_Gamma_CosPoint_B", 0, iP, cosPoint);
fQAmc->Fill("h_cut_Gamma_DCA_B", 0, iP, dca);
fQAmc->Fill("h_cut_Gamma_VtxR_B", 0, iP, r2);
- fQAmc->Fill("h_cut_Gamma_OA_B", 0, iP, oAngle);
fQAmc->Fill("h_cut_Gamma_PP_B", 0, iP, psiPair);
fQAmc->Fill("h_cut_Gamma_Chi2_B", 0, iP, chi2ndf);
fQAmc->Fill("h_cut_Gamma_Chi2_B", 1, iP, chi2ndf);
+ fQAmc->Fill("h_cut_Gamma_Sep_B", 0, iP, sep);
fQAmc->Fill("h_Electron_P_B", 0, p[0]);
fQAmc->Fill("h_Electron_P_B", 0, p[1]);
}
if(psiPair < cutPsiPair[0] || psiPair > cutPsiPair[1]) return kFALSE;
fQA->Fill("h_Gamma_Mass", 5, iMass);
if(iMass < cutMass){
- fQA->Fill("h_cut_Gamma_OA", 1, oAngle);
- fQA->Fill("h_cut_Gamma_OAvP", 1, iP, oAngle);
+ fQA->Fill("h_cut_Gamma_Sep", 1, iP, sep);
fQA->Fill("h_Electron_P", 5, p[0]);
fQA->Fill("h_Electron_P", 5, p[1]);
}
if(iMass < cutMass){
if(1 == fCurrentV0id){
- fQAmc->Fill("h_cut_Gamma_OA_S", 1, iP, oAngle);
+ fQAmc->Fill("h_cut_Gamma_Sep_S", 1, iP, sep);
fQAmc->Fill("h_Electron_P_S", 5, p[0]);
fQAmc->Fill("h_Electron_P_S", 5, p[1]);
}
else if(-2 != fCurrentV0id){
- fQAmc->Fill("h_cut_Gamma_OA_B", 1, iP, oAngle);
+ fQAmc->Fill("h_cut_Gamma_Sep_B", 1, iP, sep);
fQAmc->Fill("h_Electron_P_B", 5, p[0]);
fQAmc->Fill("h_Electron_P_B", 5, p[1]);
}
}
}
+
+ // TESTING NEW CUT
//
- // Opening angle cut (obsolete?)
+ // distance of the tracks at the entrance of the TPC
//
- if(oAngle < cutOAngle[0] || oAngle > cutOAngle[1]) return kFALSE;
+ if(sep < cutSeparation) return kFALSE;
fQA->Fill("h_Gamma_Mass", 6, iMass);
if(iMass < cutMass){
fQA->Fill("h_Electron_P", 6, p[0]);
}
if(fMCEvent){
if(1 == fCurrentV0id) fQAmc->Fill("h_Gamma_Mass_S", 6, iMass);
- else if(-2 != fCurrentV0id) fQAmc->Fill("h_Gamma_Mass_B", 6, iMass);
+ else if(-2 != fCurrentV0id)fQAmc->Fill("h_Gamma_Mass_B", 6, iMass);
+
if(iMass < cutMass){
if(1 == fCurrentV0id){
fQAmc->Fill("h_Electron_P_S", 6, p[0]);
}
}
+ // .. test
+
if(iMass > cutMass) return kFALSE;
Double_t data[4] = {0., 0., 0., 0.};
// Cut values
- const Double_t cutChi2NDF = 40.; // ORG [7.]
+ const Double_t cutChi2NDF = 10.; // ORG [7.]
const Double_t cutCosPoint[2] = {0., 0.02}; // ORG [0., 0.03]
const Double_t cutDCA[2] = {0., 0.2}; // ORG [0., 0.1]
const Double_t cutProdVtxR[2] = {2.0, 30.}; // ORG [0., 8.1]
- const Double_t cutMass[2] = {0.49, 0.51}; // ORG [0.485, 0.51]
- //const Double_t cutOAngleP = (1.0/(iP + 0.3) - 0.1); // momentum dependent min. OAngle ~ 1/x
+ const Double_t cutMass[2] = {0.486, 0.508}; // ORG [0.485, 0.51]
// Values
// cos pointing angle
if(kfMother) delete kfMother;
- // Opening angle
- Double_t oAngle = OpenAngle(v0);
-
//
// Apply the cuts, produce QA plots (with mass cut)
//
fQAmc->Fill("h_cut_K0_VtxR_S", 0, iP, r);
fQAmc->Fill("h_cut_K0_Chi2_S", 0, iP, chi2ndf);
fQAmc->Fill("h_cut_K0_Chi2_S", 1, iP, chi2ndf);
- fQAmc->Fill("h_cut_K0_OA_S", 0, iP, oAngle);
fQAmc->Fill("h_PionK0_P_S", 0, p[0]);
fQAmc->Fill("h_PionK0_P_S", 0, p[1]);
}
fQAmc->Fill("h_cut_K0_VtxR_B", 0, iP, r);
fQAmc->Fill("h_cut_K0_Chi2_B", 0, iP, chi2ndf);
fQAmc->Fill("h_cut_K0_Chi2_B", 1, iP, chi2ndf);
- fQAmc->Fill("h_cut_K0_OA_B", 0, iP, oAngle);
fQAmc->Fill("h_PionK0_P_B", 0, p[0]);
fQAmc->Fill("h_PionK0_P_B", 0, p[1]);
}
if(iMass > cutMass[0] && iMass < cutMass[1]){
fQA->Fill("h_PionK0_P", 4, p[0]);
fQA->Fill("h_PionK0_P", 4, p[1]);
- fQA->Fill("h_cut_K0_OAvP", 1, iP, oAngle);
}
if(fMCEvent){
if(2 == fCurrentV0id) fQAmc->Fill("h_K0_Mass_S", 4, iMass);
else if(-2 != fCurrentV0id) fQAmc->Fill("h_K0_Mass_B", 4, iMass);
if(iMass > cutMass[0] && iMass < cutMass[1]){
if(2 == fCurrentV0id){
- fQAmc->Fill("h_cut_K0_OA_S", 1, iP, oAngle);
fQAmc->Fill("h_PionK0_P_S", 4, p[0]);
fQAmc->Fill("h_PionK0_P_S", 4, p[1]);
}
else if(-2 != fCurrentV0id){
- fQAmc->Fill("h_cut_K0_OA_B", 1, iP, oAngle);
fQAmc->Fill("h_PionK0_P_B", 4, p[0]);
fQAmc->Fill("h_PionK0_P_B", 4, p[1]);
}
Float_t iP = v0->P();
// Cuts
- const Double_t cutChi2NDF = 40.; // ORG [5.]
+ const Double_t cutChi2NDF = 10.; // ORG [5.]
const Double_t cutCosPoint[2] = {0., 0.02}; // ORG [0., 0.03]
const Double_t cutDCA[2] = {0., 0.2}; // ORG [0., 0.2]
const Double_t cutProdVtxR[2] = {2., 40.}; // ORG [0., 24.]
const Double_t cutMass[2] = {1.11, 1.12}; // ORG [1.11, 1.12]
// cundidate cuts
- // opening angle as a function of L momentum
- //const Double_t cutOAngleP = 0.3 - 0.2*iP; // momentum dependent min. OAngle linear cut
// relative daughter momentum versusu mother momentum
// compute the cut values
if(kfMother[0]) delete kfMother[0];
if(kfMother[1]) delete kfMother[1];
- // Opening angle
- Double_t oAngle = OpenAngle(v0);
-
// Relative daughter momentum
Double_t rP = (0 == check[0]) ? p[1]/p[0] : p[0]/p[1];
fQA->Fill("h_cut_L_CosPoint", 0, cosPoint);
fQA->Fill("h_cut_L_DCA", 0, dca);
fQA->Fill("h_cut_L_VtxR", 0, r);
- fQA->Fill("h_cut_L_OAvP", 0, iP, oAngle);
fQA->Fill("h_cut_L_rdp_v_mp", 0, iP, rP);
}
if(fMCEvent){
fQAmc->Fill("h_cut_L_CosPoint_S", 0, iP, cosPoint);
fQAmc->Fill("h_cut_L_DCA_S", 0, iP, dca);
fQAmc->Fill("h_cut_L_VtxR_S", 0, iP, r);
- fQAmc->Fill("h_cut_L_OA_S", 0, iP, oAngle);
fQAmc->Fill("h_cut_L_rdp_v_mp_S", 0, iP, rP);
fQAmc->Fill("h_ProtonL_P_S", 0, p[ixMC[0]]);
fQAmc->Fill("h_PionL_P_S", 0, p[ixMC[1]]);
fQAmc->Fill("h_cut_L_CosPoint_B", 0, iP, cosPoint);
fQAmc->Fill("h_cut_L_DCA_B", 0, iP, dca);
fQAmc->Fill("h_cut_L_VtxR_B", 0, iP, r);
- fQAmc->Fill("h_cut_L_OA_B", 0, iP, oAngle);
fQAmc->Fill("h_cut_L_rdp_v_mp_B", 0, iP, rP);
fQAmc->Fill("h_ProtonL_P_B", 0, p[ixMC[0]]);
fQAmc->Fill("h_PionL_P_B", 0, p[ixMC[1]]);
if(r < cutProdVtxR[0] || r > cutProdVtxR[1]) return kFALSE;
(type == 0) ? fQA->Fill("h_L_Mass", 4, iMass) : fQA->Fill("h_AL_Mass", 4, iMass);
if(iMass > cutMass[0] && iMass < cutMass[1]){
- fQA->Fill("h_cut_L_OAvP", 1, iP, oAngle);
fQA->Fill("h_ProtonL_P", 4, p[ix[0]]);
fQA->Fill("h_PionL_P", 4, p[ix[1]]);
}
else if(-2 != fCurrentV0id) fQAmc->Fill("h_LAL_Mass_B", 4, iMass);
if(iMass > cutMass[0] && iMass < cutMass[1]){
if(4 == TMath::Abs(fCurrentV0id)){
- fQAmc->Fill("h_cut_L_OA_S", 1, iP, oAngle);
fQAmc->Fill("h_ProtonL_P_S", 4, p[ixMC[0]]);
fQAmc->Fill("h_PionL_P_S", 4, p[ixMC[1]]);
}
else if(-2 != fCurrentV0id){
- fQAmc->Fill("h_cut_L_OA_B", 1, iP, oAngle);
fQAmc->Fill("h_ProtonL_P_B", 4, p[ixMC[0]]);
fQAmc->Fill("h_PionL_P_B", 4, p[ixMC[1]]);
}