fHistDcaV0DaughtersK0vsMassK0(0),
fHistCosPointAngleK0vsMassK0(0),
- fHistDcaPosToPrimVertexK0vsMassK0pt1(0),
- fHistDcaNegToPrimVertexK0vsMassK0pt1(0),
- fHistRadiusV0K0vsMassK0pt1(0),
- fHistDecayLengthV0K0vsMassK0pt1(0),
- fHistDcaV0DaughtersK0vsMassK0pt1(0),
- fHistCosPointAngleK0vsMassK0pt1(0),
-
- fHistDcaPosToPrimVertexK0vsMassK0pt2(0),
- fHistDcaNegToPrimVertexK0vsMassK0pt2(0),
- fHistRadiusV0K0vsMassK0pt2(0),
- fHistDecayLengthV0K0vsMassK0pt2(0),
- fHistDcaV0DaughtersK0vsMassK0pt2(0),
- fHistCosPointAngleK0vsMassK0pt2(0),
-
- fHistDcaPosToPrimVertexK0vsMassK0pt3(0),
- fHistDcaNegToPrimVertexK0vsMassK0pt3(0),
- fHistRadiusV0K0vsMassK0pt3(0),
- fHistDecayLengthV0K0vsMassK0pt3(0),
- fHistDcaV0DaughtersK0vsMassK0pt3(0),
- fHistCosPointAngleK0vsMassK0pt3(0),
+ fHistDcaPosToPrimVertexK0vsMassK0pt1(0),
+ fHistDcaNegToPrimVertexK0vsMassK0pt1(0),
+ fHistRadiusV0K0vsMassK0pt1(0),
+ fHistDecayLengthV0K0vsMassK0pt1(0),
+ fHistDcaV0DaughtersK0vsMassK0pt1(0),
+ fHistCosPointAngleK0vsMassK0pt1(0),
+
+ fHistDcaPosToPrimVertexK0vsMassK0pt2(0),
+ fHistDcaNegToPrimVertexK0vsMassK0pt2(0),
+ fHistRadiusV0K0vsMassK0pt2(0),
+ fHistDecayLengthV0K0vsMassK0pt2(0),
+ fHistDcaV0DaughtersK0vsMassK0pt2(0),
+ fHistCosPointAngleK0vsMassK0pt2(0),
+
+ fHistDcaPosToPrimVertexK0vsMassK0pt3(0),
+ fHistDcaNegToPrimVertexK0vsMassK0pt3(0),
+ fHistRadiusV0K0vsMassK0pt3(0),
+ fHistDecayLengthV0K0vsMassK0pt3(0),
+ fHistDcaV0DaughtersK0vsMassK0pt3(0),
+ fHistCosPointAngleK0vsMassK0pt3(0),
/////////// Lambda ///////////////////////////
- fHistDcaPosToPrimVertexL(0),
- fHistDcaNegToPrimVertexL(0),
- fHistRadiusV0L(0),
- fHistDecayLengthV0L(0),
- fHistDcaV0DaughtersL(0),
- fHistChi2L(0),
- fHistCosPointAngleL(0),
+ fHistDcaPosToPrimVertexL(0),
+ fHistDcaNegToPrimVertexL(0),
+ fHistRadiusV0L(0),
+ fHistDecayLengthV0L(0),
+ fHistDcaV0DaughtersL(0),
+ fHistChi2L(0),
+ fHistCosPointAngleL(0),
- fHistDcaPosToPrimVertexLvsMassL(0),
- fHistDcaNegToPrimVertexLvsMassL(0),
- fHistRadiusV0LvsMassL(0),
- fHistDecayLengthV0LvsMassL(0),
- fHistDcaV0DaughtersLvsMassL(0),
- fHistCosPointAngleLvsMassL(0),
+ fHistDcaPosToPrimVertexLvsMassL(0),
+ fHistDcaNegToPrimVertexLvsMassL(0),
+ fHistRadiusV0LvsMassL(0),
+ fHistDecayLengthV0LvsMassL(0),
+ fHistDcaV0DaughtersLvsMassL(0),
+ fHistCosPointAngleLvsMassL(0),
- fHistDcaPosToPrimVertexLambdaVsMasspt1(0),
- fHistDcaNegToPrimVertexLambdaVsMasspt1(0),
- fHistRadiusV0LambdaVsMasspt1(0),
- fHistDecayLengthV0LambdaVsMasspt1(0),
- fHistDcaV0DaughtersLambdaVsMasspt1(0),
- fHistCosPointAngleLambdaVsMasspt1(0),
+ fHistDcaPosToPrimVertexLambdaVsMasspt1(0),
+ fHistDcaNegToPrimVertexLambdaVsMasspt1(0),
+ fHistRadiusV0LambdaVsMasspt1(0),
+ fHistDecayLengthV0LambdaVsMasspt1(0),
+ fHistDcaV0DaughtersLambdaVsMasspt1(0),
+ fHistCosPointAngleLambdaVsMasspt1(0),
- fHistDcaPosToPrimVertexLambdaVsMasspt2(0),
- fHistDcaNegToPrimVertexLambdaVsMasspt2(0),
- fHistRadiusV0LambdaVsMasspt2(0),
- fHistDecayLengthV0LambdaVsMasspt2(0),
- fHistDcaV0DaughtersLambdaVsMasspt2(0),
- fHistCosPointAngleLambdaVsMasspt2(0),
+ fHistDcaPosToPrimVertexLambdaVsMasspt2(0),
+ fHistDcaNegToPrimVertexLambdaVsMasspt2(0),
+ fHistRadiusV0LambdaVsMasspt2(0),
+ fHistDecayLengthV0LambdaVsMasspt2(0),
+ fHistDcaV0DaughtersLambdaVsMasspt2(0),
+ fHistCosPointAngleLambdaVsMasspt2(0),
- fHistDcaPosToPrimVertexLambdaVsMasspt3(0),
- fHistDcaNegToPrimVertexLambdaVsMasspt3(0),
- fHistRadiusV0LambdaVsMasspt3(0),
- fHistDecayLengthV0LambdaVsMasspt3(0),
- fHistDcaV0DaughtersLambdaVsMasspt3(0),
- fHistCosPointAngleLambdaVsMasspt3(0),
+ fHistDcaPosToPrimVertexLambdaVsMasspt3(0),
+ fHistDcaNegToPrimVertexLambdaVsMasspt3(0),
+ fHistRadiusV0LambdaVsMasspt3(0),
+ fHistDecayLengthV0LambdaVsMasspt3(0),
+ fHistDcaV0DaughtersLambdaVsMasspt3(0),
+ fHistCosPointAngleLambdaVsMasspt3(0),
/////////Antilambda ///////////////////
- fHistDcaPosToPrimVertexAntiL(0),
- fHistDcaNegToPrimVertexAntiL(0),
- fHistRadiusV0AntiL(0),
- fHistDecayLengthV0AntiL(0),
- fHistDcaV0DaughtersAntiL(0),
- fHistChi2AntiL(0),
- fHistCosPointAngleAntiL(0),
-
- fHistDcaPosToPrimVertexAntiLvsMass(0),
- fHistDcaNegToPrimVertexAntiLvsMass(0),
- fHistRadiusV0AntiLvsMass(0),
- fHistDecayLengthV0AntiLvsMass(0),
- fHistDcaV0DaughtersAntiLvsMass(0),
- fHistCosPointAngleAntiLvsMass(0),
+ fHistDcaPosToPrimVertexAntiL(0),
+ fHistDcaNegToPrimVertexAntiL(0),
+ fHistRadiusV0AntiL(0),
+ fHistDecayLengthV0AntiL(0),
+ fHistDcaV0DaughtersAntiL(0),
+ fHistChi2AntiL(0),
+ fHistCosPointAngleAntiL(0),
+
+ fHistDcaPosToPrimVertexAntiLvsMass(0),
+ fHistDcaNegToPrimVertexAntiLvsMass(0),
+ fHistRadiusV0AntiLvsMass(0),
+ fHistDecayLengthV0AntiLvsMass(0),
+ fHistDcaV0DaughtersAntiLvsMass(0),
+ fHistCosPointAngleAntiLvsMass(0),
- fHistDcaPosToPrimVertexAntiLVsMasspt1(0),
- fHistDcaNegToPrimVertexAntiLVsMasspt1(0),
- fHistRadiusV0AntiLVsMasspt1(0),
- fHistDecayLengthV0AntiLVsMasspt1(0),
- fHistDcaV0DaughtersAntiLVsMasspt1(0),
- fHistCosPointAngleAntiLVsMasspt1(0),
-
- fHistDcaPosToPrimVertexAntiLVsMasspt2(0),
- fHistDcaNegToPrimVertexAntiLVsMasspt2(0),
- fHistRadiusV0AntiLVsMasspt2(0),
- fHistDecayLengthV0AntiLVsMasspt2(0),
- fHistDcaV0DaughtersAntiLVsMasspt2(0),
- fHistCosPointAngleAntiLVsMasspt2(0),
-
- fHistDcaPosToPrimVertexAntiLVsMasspt3(0),
- fHistDcaNegToPrimVertexAntiLVsMasspt3(0),
- fHistRadiusV0AntiLVsMasspt3(0),
- fHistDecayLengthV0AntiLVsMasspt3(0),
- fHistDcaV0DaughtersAntiLVsMasspt3(0),
- fHistCosPointAngleAntiLVsMasspt3(0),
+ fHistDcaPosToPrimVertexAntiLVsMasspt1(0),
+ fHistDcaNegToPrimVertexAntiLVsMasspt1(0),
+ fHistRadiusV0AntiLVsMasspt1(0),
+ fHistDecayLengthV0AntiLVsMasspt1(0),
+ fHistDcaV0DaughtersAntiLVsMasspt1(0),
+ fHistCosPointAngleAntiLVsMasspt1(0),
+
+ fHistDcaPosToPrimVertexAntiLVsMasspt2(0),
+ fHistDcaNegToPrimVertexAntiLVsMasspt2(0),
+ fHistRadiusV0AntiLVsMasspt2(0),
+ fHistDecayLengthV0AntiLVsMasspt2(0),
+ fHistDcaV0DaughtersAntiLVsMasspt2(0),
+ fHistCosPointAngleAntiLVsMasspt2(0),
+
+ fHistDcaPosToPrimVertexAntiLVsMasspt3(0),
+ fHistDcaNegToPrimVertexAntiLVsMasspt3(0),
+ fHistRadiusV0AntiLVsMasspt3(0),
+ fHistDecayLengthV0AntiLVsMasspt3(0),
+ fHistDcaV0DaughtersAntiLVsMasspt3(0),
+ fHistCosPointAngleAntiLVsMasspt3(0),
/////////////////////////////////////////
fHistV0Multiplicity(0),
///////////////////////////////////////
fHistPtVsMassK0(0),
fHistPtVsMassLambda(0),
- fHistPtVsMassAntiLambda(0),
+ fHistPtVsMassAntiLambda(0),
////////////////////////////////////////
fHistArmenterosPodolanski(0),
fHistK0sMassVsLambdaMass(0),
+ fHistTPCsignalPt(0),
+ fHistNSigmaProton(0),
fHistNsigmaPosPionAntiLambda(0),
fHistNsigmaNegProtonAntiLambda(0),
fHistNsigmaPosProtonLambda(0),
fHistAsMcSecondaryMotherPdgCodeLambda(0),
fHistAsMcSecondaryMotherPdgCodeAntiLambda(0),
fHistAsMcSecondaryPtLambdaFromSigma(0),
- fHistAsMcSecondaryPtAntiLambdaFromSigma(0),
- fHistTPCsignalPt(0)
+ fHistAsMcSecondaryPtAntiLambdaFromSigma(0)
+
{
// Constructor
}
fHistCosPointAngleK0vsMassK0(0),
- fHistDcaPosToPrimVertexK0vsMassK0pt1(0),
- fHistDcaNegToPrimVertexK0vsMassK0pt1(0),
- fHistRadiusV0K0vsMassK0pt1(0),
- fHistDecayLengthV0K0vsMassK0pt1(0),
- fHistDcaV0DaughtersK0vsMassK0pt1(0),
- fHistCosPointAngleK0vsMassK0pt1(0),
-
- fHistDcaPosToPrimVertexK0vsMassK0pt2(0),
- fHistDcaNegToPrimVertexK0vsMassK0pt2(0),
- fHistRadiusV0K0vsMassK0pt2(0),
- fHistDecayLengthV0K0vsMassK0pt2(0),
- fHistDcaV0DaughtersK0vsMassK0pt2(0),
- fHistCosPointAngleK0vsMassK0pt2(0),
-
- fHistDcaPosToPrimVertexK0vsMassK0pt3(0),
- fHistDcaNegToPrimVertexK0vsMassK0pt3(0),
- fHistRadiusV0K0vsMassK0pt3(0),
- fHistDecayLengthV0K0vsMassK0pt3(0),
- fHistDcaV0DaughtersK0vsMassK0pt3(0),
- fHistCosPointAngleK0vsMassK0pt3(0),
+ fHistDcaPosToPrimVertexK0vsMassK0pt1(0),
+ fHistDcaNegToPrimVertexK0vsMassK0pt1(0),
+ fHistRadiusV0K0vsMassK0pt1(0),
+ fHistDecayLengthV0K0vsMassK0pt1(0),
+ fHistDcaV0DaughtersK0vsMassK0pt1(0),
+ fHistCosPointAngleK0vsMassK0pt1(0),
+
+ fHistDcaPosToPrimVertexK0vsMassK0pt2(0),
+ fHistDcaNegToPrimVertexK0vsMassK0pt2(0),
+ fHistRadiusV0K0vsMassK0pt2(0),
+ fHistDecayLengthV0K0vsMassK0pt2(0),
+ fHistDcaV0DaughtersK0vsMassK0pt2(0),
+ fHistCosPointAngleK0vsMassK0pt2(0),
+
+ fHistDcaPosToPrimVertexK0vsMassK0pt3(0),
+ fHistDcaNegToPrimVertexK0vsMassK0pt3(0),
+ fHistRadiusV0K0vsMassK0pt3(0),
+ fHistDecayLengthV0K0vsMassK0pt3(0),
+ fHistDcaV0DaughtersK0vsMassK0pt3(0),
+ fHistCosPointAngleK0vsMassK0pt3(0),
/////////////////////////////////////////
- fHistDcaPosToPrimVertexLambdaVsMasspt1(0),
- fHistDcaNegToPrimVertexLambdaVsMasspt1(0),
- fHistRadiusV0LambdaVsMasspt1(0),
- fHistDecayLengthV0LambdaVsMasspt1(0),
- fHistDcaV0DaughtersLambdaVsMasspt1(0),
- fHistCosPointAngleLambdaVsMasspt1(0),
-
- fHistDcaPosToPrimVertexLambdaVsMasspt2(0),
- fHistDcaNegToPrimVertexLambdaVsMasspt2(0),
- fHistRadiusV0LambdaVsMasspt2(0),
- fHistDecayLengthV0LambdaVsMasspt2(0),
- fHistDcaV0DaughtersLambdaVsMasspt2(0),
- fHistCosPointAngleLambdaVsMasspt2(0),
-
- fHistDcaPosToPrimVertexLambdaVsMasspt3(0),
- fHistDcaNegToPrimVertexLambdaVsMasspt3(0),
- fHistRadiusV0LambdaVsMasspt3(0),
- fHistDecayLengthV0LambdaVsMasspt3(0),
- fHistDcaV0DaughtersLambdaVsMasspt3(0),
- fHistCosPointAngleLambdaVsMasspt3(0),
+ fHistDcaPosToPrimVertexLambdaVsMasspt1(0),
+ fHistDcaNegToPrimVertexLambdaVsMasspt1(0),
+ fHistRadiusV0LambdaVsMasspt1(0),
+ fHistDecayLengthV0LambdaVsMasspt1(0),
+ fHistDcaV0DaughtersLambdaVsMasspt1(0),
+ fHistCosPointAngleLambdaVsMasspt1(0),
+
+ fHistDcaPosToPrimVertexLambdaVsMasspt2(0),
+ fHistDcaNegToPrimVertexLambdaVsMasspt2(0),
+ fHistRadiusV0LambdaVsMasspt2(0),
+ fHistDecayLengthV0LambdaVsMasspt2(0),
+ fHistDcaV0DaughtersLambdaVsMasspt2(0),
+ fHistCosPointAngleLambdaVsMasspt2(0),
+
+ fHistDcaPosToPrimVertexLambdaVsMasspt3(0),
+ fHistDcaNegToPrimVertexLambdaVsMasspt3(0),
+ fHistRadiusV0LambdaVsMasspt3(0),
+ fHistDecayLengthV0LambdaVsMasspt3(0),
+ fHistDcaV0DaughtersLambdaVsMasspt3(0),
+ fHistCosPointAngleLambdaVsMasspt3(0),
////////////AntiLambda
fHistDcaPosToPrimVertexAntiL(0),
- fHistDcaPosToPrimVertexAntiLVsMasspt1(0),
- fHistDcaNegToPrimVertexAntiLVsMasspt1(0),
- fHistRadiusV0AntiLVsMasspt1(0),
- fHistDecayLengthV0AntiLVsMasspt1(0),
- fHistDcaV0DaughtersAntiLVsMasspt1(0),
- fHistCosPointAngleAntiLVsMasspt1(0),
-
- fHistDcaPosToPrimVertexAntiLVsMasspt2(0),
- fHistDcaNegToPrimVertexAntiLVsMasspt2(0),
- fHistRadiusV0AntiLVsMasspt2(0),
- fHistDecayLengthV0AntiLVsMasspt2(0),
- fHistDcaV0DaughtersAntiLVsMasspt2(0),
- fHistCosPointAngleAntiLVsMasspt2(0),
-
- fHistDcaPosToPrimVertexAntiLVsMasspt3(0),
- fHistDcaNegToPrimVertexAntiLVsMasspt3(0),
- fHistRadiusV0AntiLVsMasspt3(0),
- fHistDecayLengthV0AntiLVsMasspt3(0),
- fHistDcaV0DaughtersAntiLVsMasspt3(0),
- fHistCosPointAngleAntiLVsMasspt3(0),
+ fHistDcaPosToPrimVertexAntiLVsMasspt1(0),
+ fHistDcaNegToPrimVertexAntiLVsMasspt1(0),
+ fHistRadiusV0AntiLVsMasspt1(0),
+ fHistDecayLengthV0AntiLVsMasspt1(0),
+ fHistDcaV0DaughtersAntiLVsMasspt1(0),
+ fHistCosPointAngleAntiLVsMasspt1(0),
+
+ fHistDcaPosToPrimVertexAntiLVsMasspt2(0),
+ fHistDcaNegToPrimVertexAntiLVsMasspt2(0),
+ fHistRadiusV0AntiLVsMasspt2(0),
+ fHistDecayLengthV0AntiLVsMasspt2(0),
+ fHistDcaV0DaughtersAntiLVsMasspt2(0),
+ fHistCosPointAngleAntiLVsMasspt2(0),
+
+ fHistDcaPosToPrimVertexAntiLVsMasspt3(0),
+ fHistDcaNegToPrimVertexAntiLVsMasspt3(0),
+ fHistRadiusV0AntiLVsMasspt3(0),
+ fHistDecayLengthV0AntiLVsMasspt3(0),
+ fHistDcaV0DaughtersAntiLVsMasspt3(0),
+ fHistCosPointAngleAntiLVsMasspt3(0),
/////////////////////////////////////////
///////////////////////////////////////////////////
fHistArmenterosPodolanski(0),
fHistK0sMassVsLambdaMass(0),
+ fHistTPCsignalPt(0),
+ fHistNSigmaProton(0),
fHistNsigmaPosPionAntiLambda(0),
fHistNsigmaNegProtonAntiLambda(0),
fHistNsigmaPosProtonLambda(0),
fHistAsMcSecondaryMotherPdgCodeLambda(0),
fHistAsMcSecondaryMotherPdgCodeAntiLambda(0),
fHistAsMcSecondaryPtLambdaFromSigma(0),
- fHistAsMcSecondaryPtAntiLambdaFromSigma(0),
- fHistTPCsignalPt(0)
+ fHistAsMcSecondaryPtAntiLambdaFromSigma(0)
{
// Constructor
fHistMCPtLambda = new TH1F("h1MCPtLambda", "#Lambda^{0};p_{t} (GeV/c)",240,0,12);
fListHist->Add(fHistMCPtLambda);
- fHistMCPtAntiLambda = new TH1F("h1MCPtAntiLambda", "#AntiLambda^{0};p_{t} (GeV/c)",240,0,12);
+fHistMCPtAntiLambda = new TH1F("h1MCPtAntiLambda", "#AntiLambda^{0};p_{t} (GeV/c)",240,0,12);
fListHist->Add(fHistMCPtAntiLambda);
// Pt distribution of Lambda coming from Sigma decay
fHistMCDaughterTrack = new TH1F("h1MCDaughterTrack","Distribution of mc id for daughters;id tags;Counts",15,0,15);
fListHist->Add(fHistMCDaughterTrack);
- fHistTPCTracks = new TH1F("h1TPCTracks","Distribution of TPC tracks;Number of TPC tracks:Number of events",1000,0,10000);
+ fHistTPCTracks = new TH1F("h1TPCTracks","Distribution of TPC tracks;Number of TPC tracks:Number of events",1000,0,10000);
fListHist->Add(fHistTPCTracks);
// Primary Vertex:
fHistCosPointAngleK0vsMassK0 = new TH2F("h2CosPointAngleK0vsMassK0", "Cosine of V0's pointing angle", 200,0.997,1.007,200,0.4,0.6);
fListHist->Add(fHistCosPointAngleK0vsMassK0);
- //// pt1
- fHistDcaPosToPrimVertexK0vsMassK0pt1 = new TH2F("h2DcaPosToPrimVertexK0vsMassK0pt1", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
- fListHist->Add(fHistDcaPosToPrimVertexK0vsMassK0pt1);
+ //// pt1
+ fHistDcaPosToPrimVertexK0vsMassK0pt1 = new TH2F("h2DcaPosToPrimVertexK0vsMassK0pt1", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
+ fListHist->Add(fHistDcaPosToPrimVertexK0vsMassK0pt1);
- fHistDcaNegToPrimVertexK0vsMassK0pt1 = new TH2F("h2DcaNegToPrimVertexK0vsMassK0pt1", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
- fListHist->Add(fHistDcaNegToPrimVertexK0vsMassK0pt1);
+ fHistDcaNegToPrimVertexK0vsMassK0pt1 = new TH2F("h2DcaNegToPrimVertexK0vsMassK0pt1", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
+ fListHist->Add(fHistDcaNegToPrimVertexK0vsMassK0pt1);
- fHistRadiusV0K0vsMassK0pt1 = new TH2F("h2RadiusV0K0vsMassK0pt1", "Radius;Radius(cm);K0s inv. mass",110,0,110,200,0.4,0.6);
- fListHist->Add(fHistRadiusV0K0vsMassK0pt1);
+ fHistRadiusV0K0vsMassK0pt1 = new TH2F("h2RadiusV0K0vsMassK0pt1", "Radius;Radius(cm);K0s inv. mass",110,0,110,200,0.4,0.6);
+ fListHist->Add(fHistRadiusV0K0vsMassK0pt1);
- fHistDecayLengthV0K0vsMassK0pt1 = new TH2F("h2DecayLengthV0K0vsMassK0pt1", "V0s decay Length;decay length(cm);K0s inv. mass",100,0,100,200,0.4,0.6);
- fListHist->Add(fHistDecayLengthV0K0vsMassK0pt1);
+ fHistDecayLengthV0K0vsMassK0pt1 = new TH2F("h2DecayLengthV0K0vsMassK0pt1", "V0s decay Length;decay length(cm);K0s inv. mass",100,0,100,200,0.4,0.6);
+ fListHist->Add(fHistDecayLengthV0K0vsMassK0pt1);
- fHistDcaV0DaughtersK0vsMassK0pt1 = new TH2F("h2DcaV0DaughtersK0vsMassK0pt1", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,200,0.4,0.6);
- fListHist->Add(fHistDcaV0DaughtersK0vsMassK0pt1);
+ fHistDcaV0DaughtersK0vsMassK0pt1 = new TH2F("h2DcaV0DaughtersK0vsMassK0pt1", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,200,0.4,0.6);
+ fListHist->Add(fHistDcaV0DaughtersK0vsMassK0pt1);
- fHistCosPointAngleK0vsMassK0pt1 = new TH2F("h2CosPointAngleK0vsMassK0pt1", "Cosine of V0's pointing angle", 200,0.997,1.007,200,0.4,0.6);
- fListHist->Add(fHistCosPointAngleK0vsMassK0pt1);
+ fHistCosPointAngleK0vsMassK0pt1 = new TH2F("h2CosPointAngleK0vsMassK0pt1", "Cosine of V0's pointing angle", 200,0.997,1.007,200,0.4,0.6);
+ fListHist->Add(fHistCosPointAngleK0vsMassK0pt1);
- /// pt2
- fHistDcaPosToPrimVertexK0vsMassK0pt2 = new TH2F("h2DcaPosToPrimVertexK0vsMassK0pt2", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
- fListHist->Add(fHistDcaPosToPrimVertexK0vsMassK0pt2);
+ /// pt2
+ fHistDcaPosToPrimVertexK0vsMassK0pt2 = new TH2F("h2DcaPosToPrimVertexK0vsMassK0pt2", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
+ fListHist->Add(fHistDcaPosToPrimVertexK0vsMassK0pt2);
- fHistDcaNegToPrimVertexK0vsMassK0pt2 = new TH2F("h2DcaNegToPrimVertexK0vsMassK0pt2", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
- fListHist->Add(fHistDcaNegToPrimVertexK0vsMassK0pt2);
+ fHistDcaNegToPrimVertexK0vsMassK0pt2 = new TH2F("h2DcaNegToPrimVertexK0vsMassK0pt2", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
+ fListHist->Add(fHistDcaNegToPrimVertexK0vsMassK0pt2);
- fHistRadiusV0K0vsMassK0pt2 = new TH2F("h2RadiusV0K0vsMassK0pt2", "Radius;Radius(cm);K0s inv. mass",110,0,110,200,0.4,0.6);
- fListHist->Add(fHistRadiusV0K0vsMassK0pt2);
+ fHistRadiusV0K0vsMassK0pt2 = new TH2F("h2RadiusV0K0vsMassK0pt2", "Radius;Radius(cm);K0s inv. mass",110,0,110,200,0.4,0.6);
+ fListHist->Add(fHistRadiusV0K0vsMassK0pt2);
- fHistDecayLengthV0K0vsMassK0pt2 = new TH2F("h2DecayLengthV0K0vsMassK0pt2", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,200,0.4,0.6);
- fListHist->Add(fHistDecayLengthV0K0vsMassK0pt2);
+ fHistDecayLengthV0K0vsMassK0pt2 = new TH2F("h2DecayLengthV0K0vsMassK0pt2", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,200,0.4,0.6);
+ fListHist->Add(fHistDecayLengthV0K0vsMassK0pt2);
- fHistDcaV0DaughtersK0vsMassK0pt2 = new TH2F("h2DcaV0DaughtersK0vsMassK0pt2", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,200,0.4,0.6);
- fListHist->Add(fHistDcaV0DaughtersK0vsMassK0pt2);
+ fHistDcaV0DaughtersK0vsMassK0pt2 = new TH2F("h2DcaV0DaughtersK0vsMassK0pt2", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,200,0.4,0.6);
+ fListHist->Add(fHistDcaV0DaughtersK0vsMassK0pt2);
- fHistCosPointAngleK0vsMassK0pt2 = new TH2F("h2CosPointAngleK0vsMassK0pt2", "Cosine of V0's pointing angle", 200,0.997,1.007,200,0.4,0.6);
- fListHist->Add(fHistCosPointAngleK0vsMassK0pt2);
+ fHistCosPointAngleK0vsMassK0pt2 = new TH2F("h2CosPointAngleK0vsMassK0pt2", "Cosine of V0's pointing angle", 200,0.997,1.007,200,0.4,0.6);
+ fListHist->Add(fHistCosPointAngleK0vsMassK0pt2);
- /// pt3
- fHistDcaPosToPrimVertexK0vsMassK0pt3 = new TH2F("h2DcaPosToPrimVertexK0vsMassK0pt3", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
- fListHist->Add(fHistDcaPosToPrimVertexK0vsMassK0pt3);
+ /// pt3
+ fHistDcaPosToPrimVertexK0vsMassK0pt3 = new TH2F("h2DcaPosToPrimVertexK0vsMassK0pt3", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
+ fListHist->Add(fHistDcaPosToPrimVertexK0vsMassK0pt3);
- fHistDcaNegToPrimVertexK0vsMassK0pt3 = new TH2F("h2DcaNegToPrimVertexK0vsMassK0pt3", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
- fListHist->Add(fHistDcaNegToPrimVertexK0vsMassK0pt3);
+ fHistDcaNegToPrimVertexK0vsMassK0pt3 = new TH2F("h2DcaNegToPrimVertexK0vsMassK0pt3", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,200,0.4,0.6);
+ fListHist->Add(fHistDcaNegToPrimVertexK0vsMassK0pt3);
- fHistRadiusV0K0vsMassK0pt3 = new TH2F("h2RadiusV0K0vsMassK0pt3", "Radius;Radius(cm);K0s inv. mass",110,0,110,200,0.4,0.6);
- fListHist->Add(fHistRadiusV0K0vsMassK0pt3);
+ fHistRadiusV0K0vsMassK0pt3 = new TH2F("h2RadiusV0K0vsMassK0pt3", "Radius;Radius(cm);K0s inv. mass",110,0,110,200,0.4,0.6);
+ fListHist->Add(fHistRadiusV0K0vsMassK0pt3);
- fHistDecayLengthV0K0vsMassK0pt3 = new TH2F("h2DecayLengthV0K0vsMassK0pt3", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,200,0.4,0.6);
- fListHist->Add(fHistDecayLengthV0K0vsMassK0pt3);
+ fHistDecayLengthV0K0vsMassK0pt3 = new TH2F("h2DecayLengthV0K0vsMassK0pt3", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,200,0.4,0.6);
+ fListHist->Add(fHistDecayLengthV0K0vsMassK0pt3);
- fHistDcaV0DaughtersK0vsMassK0pt3 = new TH2F("h2DcaV0DaughtersK0vsMassK0pt3", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,200,0.4,0.6);
- fListHist->Add(fHistDcaV0DaughtersK0vsMassK0pt3);
+ fHistDcaV0DaughtersK0vsMassK0pt3 = new TH2F("h2DcaV0DaughtersK0vsMassK0pt3", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,200,0.4,0.6);
+ fListHist->Add(fHistDcaV0DaughtersK0vsMassK0pt3);
- fHistCosPointAngleK0vsMassK0pt3 = new TH2F("h2CosPointAngleK0vsMassK0pt3", "Cosine of V0's pointing angle", 200,0.997,1.007,200,0.4,0.6);
- fListHist->Add(fHistCosPointAngleK0vsMassK0pt3);
+ fHistCosPointAngleK0vsMassK0pt3 = new TH2F("h2CosPointAngleK0vsMassK0pt3", "Cosine of V0's pointing angle", 200,0.997,1.007,200,0.4,0.6);
+ fListHist->Add(fHistCosPointAngleK0vsMassK0pt3);
//////////Lambda////////////// 2D histos: cut vs on fly status////
fHistCosPointAngleLvsMassL = new TH2F("h2CosPointAngleLvsMassL", "Cosine of V0's pointing angle", 200,0.997,1.007,140, 1.06, 1.2);
fListHist->Add(fHistCosPointAngleLvsMassL);
- //// pt1
- fHistDcaPosToPrimVertexLambdaVsMasspt1 = new TH2F("h2DcaPosToPrimVertexLambdaVsMasspt1", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
- fListHist->Add(fHistDcaPosToPrimVertexLambdaVsMasspt1);
+ //// pt1
+ fHistDcaPosToPrimVertexLambdaVsMasspt1 = new TH2F("h2DcaPosToPrimVertexLambdaVsMasspt1", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
+ fListHist->Add(fHistDcaPosToPrimVertexLambdaVsMasspt1);
- fHistDcaNegToPrimVertexLambdaVsMasspt1 = new TH2F("h2DcaNegToPrimVertexLambdaVsMasspt1", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
- fListHist->Add(fHistDcaNegToPrimVertexLambdaVsMasspt1);
+ fHistDcaNegToPrimVertexLambdaVsMasspt1 = new TH2F("h2DcaNegToPrimVertexLambdaVsMasspt1", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
+ fListHist->Add(fHistDcaNegToPrimVertexLambdaVsMasspt1);
- fHistRadiusV0LambdaVsMasspt1 = new TH2F("h2RadiusV0LambdaVsMasspt1", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
- fListHist->Add(fHistRadiusV0LambdaVsMasspt1);
+ fHistRadiusV0LambdaVsMasspt1 = new TH2F("h2RadiusV0LambdaVsMasspt1", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
+ fListHist->Add(fHistRadiusV0LambdaVsMasspt1);
- fHistDecayLengthV0LambdaVsMasspt1 = new TH2F("h2DecayLengthV0LambdaVsMasspt1", "V0s decay Length;decay length(cm);K0s inv. mass",100,0,100,140,1.06,1.2);
- fListHist->Add(fHistDecayLengthV0LambdaVsMasspt1);
+ fHistDecayLengthV0LambdaVsMasspt1 = new TH2F("h2DecayLengthV0LambdaVsMasspt1", "V0s decay Length;decay length(cm);K0s inv. mass",100,0,100,140,1.06,1.2);
+ fListHist->Add(fHistDecayLengthV0LambdaVsMasspt1);
- fHistDcaV0DaughtersLambdaVsMasspt1 = new TH2F("h2DcaV0DaughtersLambdaVsMasspt1", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
- fListHist->Add(fHistDcaV0DaughtersLambdaVsMasspt1);
+ fHistDcaV0DaughtersLambdaVsMasspt1 = new TH2F("h2DcaV0DaughtersLambdaVsMasspt1", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
+ fListHist->Add(fHistDcaV0DaughtersLambdaVsMasspt1);
- fHistCosPointAngleLambdaVsMasspt1 = new TH2F("h2CosPointAngleLambdaVsMasspt1", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
- fListHist->Add(fHistCosPointAngleLambdaVsMasspt1);
+ fHistCosPointAngleLambdaVsMasspt1 = new TH2F("h2CosPointAngleLambdaVsMasspt1", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
+ fListHist->Add(fHistCosPointAngleLambdaVsMasspt1);
- /// pt2
- fHistDcaPosToPrimVertexLambdaVsMasspt2 = new TH2F("h2DcaPosToPrimVertexLambdaVsMasspt2", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
- fListHist->Add(fHistDcaPosToPrimVertexLambdaVsMasspt2);
+ /// pt2
+ fHistDcaPosToPrimVertexLambdaVsMasspt2 = new TH2F("h2DcaPosToPrimVertexLambdaVsMasspt2", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
+ fListHist->Add(fHistDcaPosToPrimVertexLambdaVsMasspt2);
- fHistDcaNegToPrimVertexLambdaVsMasspt2 = new TH2F("h2DcaNegToPrimVertexLambdaVsMasspt2", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
- fListHist->Add(fHistDcaNegToPrimVertexLambdaVsMasspt2);
+ fHistDcaNegToPrimVertexLambdaVsMasspt2 = new TH2F("h2DcaNegToPrimVertexLambdaVsMasspt2", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
+ fListHist->Add(fHistDcaNegToPrimVertexLambdaVsMasspt2);
- fHistRadiusV0LambdaVsMasspt2 = new TH2F("h2RadiusV0LambdaVsMasspt2", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
- fListHist->Add(fHistRadiusV0LambdaVsMasspt2);
+ fHistRadiusV0LambdaVsMasspt2 = new TH2F("h2RadiusV0LambdaVsMasspt2", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
+ fListHist->Add(fHistRadiusV0LambdaVsMasspt2);
- fHistDecayLengthV0LambdaVsMasspt2 = new TH2F("h2DecayLengthV0LambdaVsMasspt2", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,140,1.06,1.2);
- fListHist->Add(fHistDecayLengthV0LambdaVsMasspt2);
+ fHistDecayLengthV0LambdaVsMasspt2 = new TH2F("h2DecayLengthV0LambdaVsMasspt2", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,140,1.06,1.2);
+ fListHist->Add(fHistDecayLengthV0LambdaVsMasspt2);
- fHistDcaV0DaughtersLambdaVsMasspt2 = new TH2F("h2DcaV0DaughtersLambdaVsMasspt2", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
- fListHist->Add(fHistDcaV0DaughtersLambdaVsMasspt2);
+ fHistDcaV0DaughtersLambdaVsMasspt2 = new TH2F("h2DcaV0DaughtersLambdaVsMasspt2", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
+ fListHist->Add(fHistDcaV0DaughtersLambdaVsMasspt2);
- fHistCosPointAngleLambdaVsMasspt2 = new TH2F("h2CosPointAngleLambdaVsMasspt2", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
- fListHist->Add(fHistCosPointAngleLambdaVsMasspt2);
+ fHistCosPointAngleLambdaVsMasspt2 = new TH2F("h2CosPointAngleLambdaVsMasspt2", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
+ fListHist->Add(fHistCosPointAngleLambdaVsMasspt2);
- /// pt3
- fHistDcaPosToPrimVertexLambdaVsMasspt3 = new TH2F("h2DcaPosToPrimVertexLambdaVsMasspt3", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
- fListHist->Add(fHistDcaPosToPrimVertexLambdaVsMasspt3);
+ /// pt3
+ fHistDcaPosToPrimVertexLambdaVsMasspt3 = new TH2F("h2DcaPosToPrimVertexLambdaVsMasspt3", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
+ fListHist->Add(fHistDcaPosToPrimVertexLambdaVsMasspt3);
- fHistDcaNegToPrimVertexLambdaVsMasspt3 = new TH2F("h2DcaNegToPrimVertexLambdaVsMasspt3", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
- fListHist->Add(fHistDcaNegToPrimVertexLambdaVsMasspt3);
+ fHistDcaNegToPrimVertexLambdaVsMasspt3 = new TH2F("h2DcaNegToPrimVertexLambdaVsMasspt3", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
+ fListHist->Add(fHistDcaNegToPrimVertexLambdaVsMasspt3);
- fHistRadiusV0LambdaVsMasspt3 = new TH2F("h2RadiusV0LambdaVsMasspt3", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
- fListHist->Add(fHistRadiusV0LambdaVsMasspt3);
+ fHistRadiusV0LambdaVsMasspt3 = new TH2F("h2RadiusV0LambdaVsMasspt3", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
+ fListHist->Add(fHistRadiusV0LambdaVsMasspt3);
- fHistDecayLengthV0LambdaVsMasspt3 = new TH2F("h2DecayLengthV0LambdaVsMasspt3", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,140,1.06,1.2);
- fListHist->Add(fHistDecayLengthV0LambdaVsMasspt3);
+ fHistDecayLengthV0LambdaVsMasspt3 = new TH2F("h2DecayLengthV0LambdaVsMasspt3", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,140,1.06,1.2);
+ fListHist->Add(fHistDecayLengthV0LambdaVsMasspt3);
- fHistDcaV0DaughtersLambdaVsMasspt3 = new TH2F("h2DcaV0DaughtersLambdaVsMasspt3", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
- fListHist->Add(fHistDcaV0DaughtersLambdaVsMasspt3);
+ fHistDcaV0DaughtersLambdaVsMasspt3 = new TH2F("h2DcaV0DaughtersLambdaVsMasspt3", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
+ fListHist->Add(fHistDcaV0DaughtersLambdaVsMasspt3);
- fHistCosPointAngleLambdaVsMasspt3 = new TH2F("h2CosPointAngleLambdaVsMasspt3", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
- fListHist->Add(fHistCosPointAngleLambdaVsMasspt3);
+ fHistCosPointAngleLambdaVsMasspt3 = new TH2F("h2CosPointAngleLambdaVsMasspt3", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
+ fListHist->Add(fHistCosPointAngleLambdaVsMasspt3);
//////////AntiLambda////////////// 2D histos: cut vs on fly status////
fHistCosPointAngleAntiLvsMass = new TH2F("h2CosPointAngleAntiLvsMass", "Cosine of V0's pointing angle", 200,0.997,1.007,140, 1.06, 1.2);
fListHist->Add(fHistCosPointAngleAntiLvsMass);
- //// pt1
- fHistDcaPosToPrimVertexAntiLVsMasspt1 = new TH2F("h2DcaPosToPrimVertexAntiLVsMasspt1", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
- fListHist->Add(fHistDcaPosToPrimVertexAntiLVsMasspt1);
+ //// pt1
+ fHistDcaPosToPrimVertexAntiLVsMasspt1 = new TH2F("h2DcaPosToPrimVertexAntiLVsMasspt1", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
+ fListHist->Add(fHistDcaPosToPrimVertexAntiLVsMasspt1);
- fHistDcaNegToPrimVertexAntiLVsMasspt1 = new TH2F("h2DcaNegToPrimVertexAntiLVsMasspt1", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
- fListHist->Add(fHistDcaNegToPrimVertexAntiLVsMasspt1);
+ fHistDcaNegToPrimVertexAntiLVsMasspt1 = new TH2F("h2DcaNegToPrimVertexAntiLVsMasspt1", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
+ fListHist->Add(fHistDcaNegToPrimVertexAntiLVsMasspt1);
- fHistRadiusV0AntiLVsMasspt1 = new TH2F("h2RadiusV0AntiLVsMasspt1", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
- fListHist->Add(fHistRadiusV0AntiLVsMasspt1);
+ fHistRadiusV0AntiLVsMasspt1 = new TH2F("h2RadiusV0AntiLVsMasspt1", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
+ fListHist->Add(fHistRadiusV0AntiLVsMasspt1);
- fHistDecayLengthV0AntiLVsMasspt1 = new TH2F("h2DecayLengthV0AntiLVsMasspt1", "V0s decay Length;decay length(cm);K0s inv. mass",100,0,100,140,1.06,1.2);
- fListHist->Add(fHistDecayLengthV0AntiLVsMasspt1);
+ fHistDecayLengthV0AntiLVsMasspt1 = new TH2F("h2DecayLengthV0AntiLVsMasspt1", "V0s decay Length;decay length(cm);K0s inv. mass",100,0,100,140,1.06,1.2);
+ fListHist->Add(fHistDecayLengthV0AntiLVsMasspt1);
- fHistDcaV0DaughtersAntiLVsMasspt1 = new TH2F("h2DcaV0DaughtersAntiLVsMasspt1", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
- fListHist->Add(fHistDcaV0DaughtersAntiLVsMasspt1);
+ fHistDcaV0DaughtersAntiLVsMasspt1 = new TH2F("h2DcaV0DaughtersAntiLVsMasspt1", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
+ fListHist->Add(fHistDcaV0DaughtersAntiLVsMasspt1);
- fHistCosPointAngleAntiLVsMasspt1 = new TH2F("h2CosPointAngleAntiLVsMasspt1", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
- fListHist->Add(fHistCosPointAngleAntiLVsMasspt1);
+ fHistCosPointAngleAntiLVsMasspt1 = new TH2F("h2CosPointAngleAntiLVsMasspt1", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
+ fListHist->Add(fHistCosPointAngleAntiLVsMasspt1);
- /// pt2
- fHistDcaPosToPrimVertexAntiLVsMasspt2 = new TH2F("h2DcaPosToPrimVertexAntiLVsMasspt2", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
- fListHist->Add(fHistDcaPosToPrimVertexAntiLVsMasspt2);
+ /// pt2
+ fHistDcaPosToPrimVertexAntiLVsMasspt2 = new TH2F("h2DcaPosToPrimVertexAntiLVsMasspt2", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
+ fListHist->Add(fHistDcaPosToPrimVertexAntiLVsMasspt2);
- fHistDcaNegToPrimVertexAntiLVsMasspt2 = new TH2F("h2DcaNegToPrimVertexAntiLVsMasspt2", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
- fListHist->Add(fHistDcaNegToPrimVertexAntiLVsMasspt2);
+ fHistDcaNegToPrimVertexAntiLVsMasspt2 = new TH2F("h2DcaNegToPrimVertexAntiLVsMasspt2", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
+ fListHist->Add(fHistDcaNegToPrimVertexAntiLVsMasspt2);
- fHistRadiusV0AntiLVsMasspt2 = new TH2F("h2RadiusV0AntiLVsMasspt2", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
- fListHist->Add(fHistRadiusV0AntiLVsMasspt2);
+ fHistRadiusV0AntiLVsMasspt2 = new TH2F("h2RadiusV0AntiLVsMasspt2", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
+ fListHist->Add(fHistRadiusV0AntiLVsMasspt2);
- fHistDecayLengthV0AntiLVsMasspt2 = new TH2F("h2DecayLengthV0AntiLVsMasspt2", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,140,1.06,1.2);
- fListHist->Add(fHistDecayLengthV0AntiLVsMasspt2);
+ fHistDecayLengthV0AntiLVsMasspt2 = new TH2F("h2DecayLengthV0AntiLVsMasspt2", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,140,1.06,1.2);
+ fListHist->Add(fHistDecayLengthV0AntiLVsMasspt2);
- fHistDcaV0DaughtersAntiLVsMasspt2 = new TH2F("h2DcaV0DaughtersAntiLVsMasspt2", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
- fListHist->Add(fHistDcaV0DaughtersAntiLVsMasspt2);
+ fHistDcaV0DaughtersAntiLVsMasspt2 = new TH2F("h2DcaV0DaughtersAntiLVsMasspt2", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
+ fListHist->Add(fHistDcaV0DaughtersAntiLVsMasspt2);
- fHistCosPointAngleAntiLVsMasspt2 = new TH2F("h2CosPointAngleAntiLVsMasspt2", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
- fListHist->Add(fHistCosPointAngleAntiLVsMasspt2);
+ fHistCosPointAngleAntiLVsMasspt2 = new TH2F("h2CosPointAngleAntiLVsMasspt2", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
+ fListHist->Add(fHistCosPointAngleAntiLVsMasspt2);
- /// pt3
- fHistDcaPosToPrimVertexAntiLVsMasspt3 = new TH2F("h2DcaPosToPrimVertexAntiLVsMasspt3", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
- fListHist->Add(fHistDcaPosToPrimVertexAntiLVsMasspt3);
+ /// pt3
+ fHistDcaPosToPrimVertexAntiLVsMasspt3 = new TH2F("h2DcaPosToPrimVertexAntiLVsMasspt3", "Positive V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
+ fListHist->Add(fHistDcaPosToPrimVertexAntiLVsMasspt3);
- fHistDcaNegToPrimVertexAntiLVsMasspt3 = new TH2F("h2DcaNegToPrimVertexAntiLVsMasspt3", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
- fListHist->Add(fHistDcaNegToPrimVertexAntiLVsMasspt3);
+ fHistDcaNegToPrimVertexAntiLVsMasspt3 = new TH2F("h2DcaNegToPrimVertexAntiLVsMasspt3", "Negative V0 daughter;dca(cm);K0s inv. mass",500,0,10,140,1.06,1.2);
+ fListHist->Add(fHistDcaNegToPrimVertexAntiLVsMasspt3);
- fHistRadiusV0AntiLVsMasspt3 = new TH2F("h2RadiusV0AntiLVsMasspt3", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
- fListHist->Add(fHistRadiusV0AntiLVsMasspt3);
+ fHistRadiusV0AntiLVsMasspt3 = new TH2F("h2RadiusV0AntiLVsMasspt3", "Radius;Radius(cm);K0s inv. mass",110,0,110,140,1.06,1.2);
+ fListHist->Add(fHistRadiusV0AntiLVsMasspt3);
- fHistDecayLengthV0AntiLVsMasspt3 = new TH2F("h2DecayLengthV0AntiLVsMasspt3", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,140,1.06,1.2);
- fListHist->Add(fHistDecayLengthV0AntiLVsMasspt3);
+ fHistDecayLengthV0AntiLVsMasspt3 = new TH2F("h2DecayLengthV0AntiLVsMasspt3", "V0s decay Length;decay length(cm);K0s inv. mass", 100, 0, 100,140,1.06,1.2);
+ fListHist->Add(fHistDecayLengthV0AntiLVsMasspt3);
- fHistDcaV0DaughtersAntiLVsMasspt3 = new TH2F("h2DcaV0DaughtersAntiLVsMasspt3", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
- fListHist->Add(fHistDcaV0DaughtersAntiLVsMasspt3);
+ fHistDcaV0DaughtersAntiLVsMasspt3 = new TH2F("h2DcaV0DaughtersAntiLVsMasspt3", "DCA between daughters;dca(cm);K0s inv. mass", 110, 0, 1.1,140,1.06,1.2);
+ fListHist->Add(fHistDcaV0DaughtersAntiLVsMasspt3);
- fHistCosPointAngleAntiLVsMasspt3 = new TH2F("h2CosPointAngleAntiLVsMasspt3", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
- fListHist->Add(fHistCosPointAngleAntiLVsMasspt3);
+ fHistCosPointAngleAntiLVsMasspt3 = new TH2F("h2CosPointAngleAntiLVsMasspt3", "Cosine of V0's pointing angle", 200,0.997,1.007,140,1.06,1.2);
+ fListHist->Add(fHistCosPointAngleAntiLVsMasspt3);
// V0 Multiplicity
fHistK0sMassVsLambdaMass = new TH2F("h2HistK0sMassVsLambdaMass","K^{0} vs #Lambda^{0} candidates; M(#pi^{+}#pi^{-}) (GeV/c^{2}); M(p#pi^{-}) (GeV/c^{2})",200, 0.4, 0.6,140, 1.06, 1.2);
fListHist->Add(fHistK0sMassVsLambdaMass);
+ //dE/dx vs P daughters
+ fHistTPCsignalPt = new TH2F("h2TPCsignalVsPdaughters","TPC signal Vs p_{t} daughters; p (GeV/c);TPC signal",1000,0,2,1000,0,1000);
+ fListHist->Add(fHistTPCsignalPt);
+ fHistNSigmaProton =new TH1F("h1NSigmaProton","Number of sigmas for proton;;Count",600,0,6);
+ fListHist->Add(fHistNSigmaProton);
+
+
//PID
fHistNsigmaPosPionAntiLambda = new TH1F("h1NsigmaPosPionAntiLambda", "Positive daughter of Antilambda;NsigmaPion;Counts",25,0,5);
fListHist->Add(fHistNsigmaPosPionAntiLambda);
fHistAsMcPtLambda = new TH1F("h1AsMcPtLambda", "#Lambda^{0} associated;p_{t} (GeV/c);Counts", 240,0,12);
fListHist->Add(fHistAsMcPtLambda);
- fHistAsMcPtAntiLambda = new TH1F("h1AsMcPtAntiLambda", "#AntiLambda^{0} associated;p_{t} (GeV/c);Counts", 240,0,12);
+ fHistAsMcPtAntiLambda = new TH1F("h1AsMcPtAntiLambda", "#AntiLambda^{0} associated;p_{t} (GeV/c);Counts", 240,0,12);
fListHist->Add(fHistAsMcPtAntiLambda);
fHistAsMcPtZoomLambda = new TH1F("h1AsMcPtZoomLambda", "#Lambda^{0} candidates in -1 <y<1;p_{t} (GeV/c);Counts",20,0,1);
fListHist->Add(fHistAsMcPtZoomLambda);
- fHistAsMcPtZoomAntiLambda = new TH1F("h1AsMcPtZoomAntiLambda", "#AntiLambda^{0} candidates in -1 <y<1;p_{t} (GeV/c);Counts",20,0,1);
+ fHistAsMcPtZoomAntiLambda = new TH1F("h1AsMcPtZoomAntiLambda", "#AntiLambda^{0} candidates in -1 <y<1;p_{t} (GeV/c);Counts",20,0,1);
fListHist->Add(fHistAsMcPtZoomAntiLambda);
// Radius distribution
PostData(1, fListHist);
PostData(2, fCentrSelector);
PostData(3, fTracksCuts);
- //dE/dx vs P daughters
- fHistTPCsignalPt = new TH2F("h2TPCsignalVsPdaughters","TPC signal Vs p_{t} daughters; p_{t} (GeV/c);TPC signal",1000,0,2,1000,0,1000);
- fListHist->Add(fHistTPCsignalPt);
}
//________________________________________________________________________
return;
}
// PID
- if (fUsePID.Contains("withPID")) {
- AliAnalysisManager *man=AliAnalysisManager::GetAnalysisManager();
- AliInputEventHandler* inputHandler = (AliInputEventHandler*) (man->GetInputEventHandler());
- fPIDResponse = inputHandler->GetPIDResponse();
- }
+ if (fUsePID.Contains("withPID")) {
+ AliAnalysisManager *man=AliAnalysisManager::GetAnalysisManager();
+ AliInputEventHandler* inputHandler = (AliInputEventHandler*) (man->GetInputEventHandler());
+ fPIDResponse = inputHandler->GetPIDResponse();
+ }
//******************
//printf("found Lambda MC Plabel=%d PPDGcode=%d Nlabel=%d NPDGcode=%d\n\n",id0,lPdgCurrentDaughter0,id1,lPdgCurrentDaughter1);
}
- else
- if (lPdgcodeCurrentPart==-3122) {
- fHistMCProdRadiusAntiLambda->Fill(mcPosR);
+ else
+ if (lPdgcodeCurrentPart==-3122) {
+ fHistMCProdRadiusAntiLambda->Fill(mcPosR);
- fHistMCPtAntiLambda->Fill(lPtCurrentPart);
+ fHistMCPtAntiLambda->Fill(lPtCurrentPart);
- fHistNTimesRecAntiLambda->Fill(lNtimesReconstructedAntiLambda);
- fHistNTimesRecAntiLambdaVsPt->Fill(lPtCurrentPart,lNtimesReconstructedAntiLambda);
- if (lComeFromSigma) fHistMCPtAntiLambdaFromSigma->Fill(lPtCurrentPart);
+ fHistNTimesRecAntiLambda->Fill(lNtimesReconstructedAntiLambda);
+ fHistNTimesRecAntiLambdaVsPt->Fill(lPtCurrentPart,lNtimesReconstructedAntiLambda);
+ if (lComeFromSigma) fHistMCPtAntiLambdaFromSigma->Fill(lPtCurrentPart);
- //printf("found Lambda MC pT=%e\n",lPtCurrentPart);
- //printf("found Lambda MC Plabel=%d PPDGcode=%d Nlabel=%d NPDGcode=%d\n\n",id0,lPdgCurrentDaughter0,id1,lPdgCurrentDaughter1);
+ //printf("found Lambda MC pT=%e\n",lPtCurrentPart);
+ //printf("found Lambda MC Plabel=%d PPDGcode=%d Nlabel=%d NPDGcode=%d\n\n",id0,lPdgCurrentDaughter0,id1,lPdgCurrentDaughter1);
- }
+ }
} // end loop ESD MC
Double_t rcPosRK0s = 0;
Double_t rcPosXLambda = 0, rcPosYLambda = 0, rcPosZLambda = 0;
Double_t rcPosRLambda = 0;
- Double_t rcPosXAntiLambda = 0, rcPosYAntiLambda = 0, rcPosZAntiLambda = 0;
- Double_t rcPosRAntiLambda = 0;
+ Double_t rcPosXAntiLambda = 0, rcPosYAntiLambda = 0, rcPosZAntiLambda = 0;
+ Double_t rcPosRAntiLambda = 0;
// Pt resolution
- Double_t deltaPtK0s = 0, deltaPtLambda = 0, deltaPtAntiLambda = 0;
+ Double_t deltaPtK0s = 0, deltaPtLambda = 0, deltaPtAntiLambda = 0;
// Daughters
AliESDtrack *myTrackPos = NULL;
lMagneticField = ((AliESDEvent*)fESD)->GetMagneticField();
- fHistTPCTracks->Fill(AliESDtrackCuts::GetReferenceMultiplicity((AliESDEvent*)fESD, kTRUE));
+ fHistTPCTracks->Fill(AliESDtrackCuts::GetReferenceMultiplicity((AliESDEvent*)fESD, kTRUE));
////////////////////////////////////////////////////////////////////////////////////////
Int_t i =0;
if (fUsePID.Contains("withPID")) {
// nSigmaPosPion = TMath::Abs(fESDpid->NumberOfSigmasTPC(myTrackPos,AliPID::kPion));
nSigmaPosPion = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(myTrackPos, AliPID::kPion));
- // nSigmaNegPion = TMath::Abs(fESDpid->NumberOfSigmasTPC(myTrackNeg,AliPID::kPion));
+ // nSigmaNegPion = TMath::Abs(fESDpid->NumberOfSigmasTPC(myTrackNeg,AliPID::kPion));
nSigmaNegPion = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(myTrackNeg, AliPID::kPion));
// nSigmaPosProton = TMath::Abs(fESDpid->NumberOfSigmasTPC(myTrackPos,AliPID::kProton));
nSigmaPosProton = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(myTrackPos, AliPID::kProton));
/*
- if ( lCheckPIDAntiLambdaNegDaughter==1 || lCheckPIDAntiLambdaPosDaughter==1 || lCheckPIDLambdaPosDaughter!=1 || lCheckPIDLambdaNegDaughter!=1)
- {
- if (negPiKF) delete negPiKF; negPiKF=NULL;
- if (posPiKF) delete posPiKF; posPiKF=NULL;
- if (posPKF) delete posPKF; posPKF=NULL;
- if (negAPKF) delete negAPKF; negAPKF=NULL;
+ if ( lCheckPIDAntiLambdaNegDaughter==1 || lCheckPIDAntiLambdaPosDaughter==1 || lCheckPIDLambdaPosDaughter!=1 || lCheckPIDLambdaNegDaughter!=1)
+ {
+ if (negPiKF) delete negPiKF; negPiKF=NULL;
+ if (posPiKF) delete posPiKF; posPiKF=NULL;
+ if (posPKF) delete posPKF; posPKF=NULL;
+ if (negAPKF) delete negAPKF; negAPKF=NULL;
- continue;
- }
+ continue;
+ }
*/
//*****************************
// filling histograms
//*****************************
/*
-///////////////values for cuts/////////////////////////////////////////////////////////////////////////////////////////
-if ((lDcaPosToPrimVertex<=0.1) || (lDcaNegToPrimVertex<=0.1) || (lDcaV0Daughters>=1.00) ||
-(lV0cosPointAngle<=0.998) || (lV0Radius<=0.9) || (lV0Radius>=100) )
+ ///////////////values for cuts/////////////////////////////////////////////////////////////////////////////////////////
+ if ((lDcaPosToPrimVertex<=0.1) || (lDcaNegToPrimVertex<=0.1) || (lDcaV0Daughters>=1.00) ||
+ (lV0cosPointAngle<=0.998) || (lV0Radius<=0.9) || (lV0Radius>=100) )
-{
-if (negPiKF) delete negPiKF; negPiKF=NULL;
-if (posPiKF) delete posPiKF; posPiKF=NULL;
-if (posPKF) delete posPKF; posPKF=NULL;
-if (negAPKF) delete negAPKF; negAPKF=NULL;
-
-continue;}
-//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
-*/
+ {
+ if (negPiKF) delete negPiKF; negPiKF=NULL;
+ if (posPiKF) delete posPiKF; posPiKF=NULL;
+ if (posPKF) delete posPKF; posPKF=NULL;
+ if (negAPKF) delete negAPKF; negAPKF=NULL;
+
+ continue;}
+ //////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
+ */
- if(fUsePID.Contains("withPID") && (lCheckPIDAntiLambdaNegDaughter==0) && lCheckPIDLambdaPosDaughter==1) LambdaPID = 1;
- else LambdaPID =0;
- if(fUsePID.Contains("withPID") && (lCheckPIDLambdaPosDaughter==0) && lCheckPIDAntiLambdaNegDaughter==1) AntiLambdaPID = 1;
- else AntiLambdaPID =0;
+ if(fUsePID.Contains("withPID") && (lCheckPIDAntiLambdaNegDaughter==0) && lCheckPIDLambdaPosDaughter==1) LambdaPID = 1;
+ else LambdaPID =0;
+ if(fUsePID.Contains("withPID") && (lCheckPIDLambdaPosDaughter==0) && lCheckPIDAntiLambdaNegDaughter==1) AntiLambdaPID = 1;
+ else AntiLambdaPID =0;
- V0momentum = TMath::Sqrt(lPzLambda*lPzLambda + lPtLambda*lPtLambda);
+ V0momentum = TMath::Sqrt(lPzLambda*lPzLambda + lPtLambda*lPtLambda);
- if ((LambdaPID==1 && V0momentum <=1) || (V0momentum>1) || !(fUsePID.Contains("withPID"))){
- if ((TMath::Abs(lRapK0s) < lCutRap) && lOnFlyStatus==0) {
- fHistTPCsignalPt->Fill(V0momentum,myTrackPos->GetTPCsignal());
- fHistTPCsignalPt->Fill(V0momentum,myTrackNeg->GetTPCsignal());
- }
- }
+ if ((LambdaPID==1 && V0momentum <=1) || (V0momentum>1) || !(fUsePID.Contains("withPID"))){
+ if ((TMath::Abs(lRapK0s) < lCutRap) && lOnFlyStatus==0) {
+ fHistTPCsignalPt->Fill(V0momentum,myTrackPos->GetTPCsignal());
+ // fHistTPCsignalPt->Fill(V0momentum,myTrackNeg->GetTPCsignal());
+ }
+ }
- if ((V0momentum)<1 && lOnFlyStatus==0 ){
- fHistArmenterosPodolanski->Fill(lAlphaV0,lPtArmV0);
- }
+ if ((V0momentum)<1 && lOnFlyStatus==0 ){
+ fHistArmenterosPodolanski->Fill(lAlphaV0,lPtArmV0);
+ }
- // insert PID condition for K0s instead of "kTRUE" value
- if ((fUsePID.Contains("withPID") && kTRUE) || !(fUsePID.Contains("withPID"))){
- if (TMath::Abs(lRapK0s) < lCutRap ) {
+ // insert PID condition for K0s instead of "kTRUE" value
+ if ((fUsePID.Contains("withPID") && kTRUE) || !(fUsePID.Contains("withPID"))){
+ if (TMath::Abs(lRapK0s) < lCutRap ) {
- //////2D histos: cut vs on fly status/////////////////////
+ //////2D histos: cut vs on fly status/////////////////////
- fHistDcaPosToPrimVertexK0->Fill(lDcaPosToPrimVertex,lOnFlyStatus);
- fHistDcaNegToPrimVertexK0->Fill(lDcaNegToPrimVertex,lOnFlyStatus);
- fHistRadiusV0K0->Fill(lV0Radius,lOnFlyStatus);
- fHistDecayLengthV0K0->Fill(lV0DecayLength,lOnFlyStatus);
- fHistDcaV0DaughtersK0->Fill(lDcaV0Daughters,lOnFlyStatus);
- fHistChi2K0->Fill(lChi2V0,lOnFlyStatus);
- fHistCosPointAngleK0->Fill(lV0cosPointAngle,lOnFlyStatus);
+ fHistDcaPosToPrimVertexK0->Fill(lDcaPosToPrimVertex,lOnFlyStatus);
+ fHistDcaNegToPrimVertexK0->Fill(lDcaNegToPrimVertex,lOnFlyStatus);
+ fHistRadiusV0K0->Fill(lV0Radius,lOnFlyStatus);
+ fHistDecayLengthV0K0->Fill(lV0DecayLength,lOnFlyStatus);
+ fHistDcaV0DaughtersK0->Fill(lDcaV0Daughters,lOnFlyStatus);
+ fHistChi2K0->Fill(lChi2V0,lOnFlyStatus);
+ fHistCosPointAngleK0->Fill(lV0cosPointAngle,lOnFlyStatus);
- //////2D histos: cut vs mass/////////////////////
+ //////2D histos: cut vs mass/////////////////////
- if (lOnFlyStatus==0){
+ if (lOnFlyStatus==0){
- fHistMassK0->Fill(lInvMassK0s);
- fHistMassVsRadiusK0->Fill(rcPosRK0s,lInvMassK0s);
- fHistPtVsMassK0->Fill(lInvMassK0s,lPtK0s);
-
-
- fHistDcaPosToPrimVertexK0vsMassK0->Fill(lDcaPosToPrimVertex,lInvMassK0s);
- fHistDcaNegToPrimVertexK0vsMassK0->Fill(lDcaNegToPrimVertex,lInvMassK0s);
- fHistRadiusV0K0vsMassK0->Fill(lV0Radius,lInvMassK0s);
- fHistDecayLengthV0K0vsMassK0->Fill(lV0DecayLength,lInvMassK0s);
- fHistDcaV0DaughtersK0vsMassK0->Fill(lDcaV0Daughters,lInvMassK0s);
- fHistCosPointAngleK0vsMassK0->Fill(lV0cosPointAngle,lInvMassK0s);
-
- if (lPtK0s>0 && lPtK0s <3){
- fHistDcaPosToPrimVertexK0vsMassK0pt1->Fill(lDcaPosToPrimVertex,lInvMassK0s);
- fHistDcaNegToPrimVertexK0vsMassK0pt1->Fill(lDcaNegToPrimVertex,lInvMassK0s);
- fHistRadiusV0K0vsMassK0pt1->Fill(lV0Radius,lInvMassK0s);
- fHistDecayLengthV0K0vsMassK0pt1->Fill(lV0DecayLength,lInvMassK0s);
- fHistDcaV0DaughtersK0vsMassK0pt1->Fill(lDcaV0Daughters,lInvMassK0s);
- fHistCosPointAngleK0vsMassK0pt1->Fill(lV0cosPointAngle,lInvMassK0s);
- }
- if (lPtK0s > 3 && lPtK0s < 6){
- fHistDcaPosToPrimVertexK0vsMassK0pt2->Fill(lDcaPosToPrimVertex,lInvMassK0s);
- fHistDcaNegToPrimVertexK0vsMassK0pt2->Fill(lDcaNegToPrimVertex,lInvMassK0s);
- fHistRadiusV0K0vsMassK0pt2->Fill(lV0Radius,lInvMassK0s);
- fHistDecayLengthV0K0vsMassK0pt2->Fill(lV0DecayLength,lInvMassK0s);
- fHistDcaV0DaughtersK0vsMassK0pt2->Fill(lDcaV0Daughters,lInvMassK0s);
- fHistCosPointAngleK0vsMassK0pt2->Fill(lV0cosPointAngle,lInvMassK0s);
- }
- if (lPtK0s > 6 && lPtK0s < 10){
- fHistDcaPosToPrimVertexK0vsMassK0pt3->Fill(lDcaPosToPrimVertex,lInvMassK0s);
- fHistDcaNegToPrimVertexK0vsMassK0pt3->Fill(lDcaNegToPrimVertex,lInvMassK0s);
- fHistRadiusV0K0vsMassK0pt3->Fill(lV0Radius,lInvMassK0s);
- fHistDecayLengthV0K0vsMassK0pt3->Fill(lV0DecayLength,lInvMassK0s);
- fHistDcaV0DaughtersK0vsMassK0pt3->Fill(lDcaV0Daughters,lInvMassK0s);
- fHistCosPointAngleK0vsMassK0pt3->Fill(lV0cosPointAngle,lInvMassK0s);
- }
- }
- } // if rap. condition
- } // end if withPID condition
+ fHistMassK0->Fill(lInvMassK0s);
+ fHistMassVsRadiusK0->Fill(rcPosRK0s,lInvMassK0s);
+ fHistPtVsMassK0->Fill(lInvMassK0s,lPtK0s);
+
+
+ fHistDcaPosToPrimVertexK0vsMassK0->Fill(lDcaPosToPrimVertex,lInvMassK0s);
+ fHistDcaNegToPrimVertexK0vsMassK0->Fill(lDcaNegToPrimVertex,lInvMassK0s);
+ fHistRadiusV0K0vsMassK0->Fill(lV0Radius,lInvMassK0s);
+ fHistDecayLengthV0K0vsMassK0->Fill(lV0DecayLength,lInvMassK0s);
+ fHistDcaV0DaughtersK0vsMassK0->Fill(lDcaV0Daughters,lInvMassK0s);
+ fHistCosPointAngleK0vsMassK0->Fill(lV0cosPointAngle,lInvMassK0s);
+
+ if (lPtK0s>0 && lPtK0s <3){
+ fHistDcaPosToPrimVertexK0vsMassK0pt1->Fill(lDcaPosToPrimVertex,lInvMassK0s);
+ fHistDcaNegToPrimVertexK0vsMassK0pt1->Fill(lDcaNegToPrimVertex,lInvMassK0s);
+ fHistRadiusV0K0vsMassK0pt1->Fill(lV0Radius,lInvMassK0s);
+ fHistDecayLengthV0K0vsMassK0pt1->Fill(lV0DecayLength,lInvMassK0s);
+ fHistDcaV0DaughtersK0vsMassK0pt1->Fill(lDcaV0Daughters,lInvMassK0s);
+ fHistCosPointAngleK0vsMassK0pt1->Fill(lV0cosPointAngle,lInvMassK0s);
+ }
+ if (lPtK0s > 3 && lPtK0s < 6){
+ fHistDcaPosToPrimVertexK0vsMassK0pt2->Fill(lDcaPosToPrimVertex,lInvMassK0s);
+ fHistDcaNegToPrimVertexK0vsMassK0pt2->Fill(lDcaNegToPrimVertex,lInvMassK0s);
+ fHistRadiusV0K0vsMassK0pt2->Fill(lV0Radius,lInvMassK0s);
+ fHistDecayLengthV0K0vsMassK0pt2->Fill(lV0DecayLength,lInvMassK0s);
+ fHistDcaV0DaughtersK0vsMassK0pt2->Fill(lDcaV0Daughters,lInvMassK0s);
+ fHistCosPointAngleK0vsMassK0pt2->Fill(lV0cosPointAngle,lInvMassK0s);
+ }
+ if (lPtK0s > 6 && lPtK0s < 10){
+ fHistDcaPosToPrimVertexK0vsMassK0pt3->Fill(lDcaPosToPrimVertex,lInvMassK0s);
+ fHistDcaNegToPrimVertexK0vsMassK0pt3->Fill(lDcaNegToPrimVertex,lInvMassK0s);
+ fHistRadiusV0K0vsMassK0pt3->Fill(lV0Radius,lInvMassK0s);
+ fHistDecayLengthV0K0vsMassK0pt3->Fill(lV0DecayLength,lInvMassK0s);
+ fHistDcaV0DaughtersK0vsMassK0pt3->Fill(lDcaV0Daughters,lInvMassK0s);
+ fHistCosPointAngleK0vsMassK0pt3->Fill(lV0cosPointAngle,lInvMassK0s);
+ }
+ }
+ } // if rap. condition
+ } // end if withPID condition
- // insert PID condition for Lambda instead of "kTRUE" value
- // if ((fUsePID.Contains("withPID") && kTRUE )|| !(fUsePID.Contains("withPID"))){
- if ((LambdaPID==1 && V0momentum <=1) || (V0momentum>1) || !(fUsePID.Contains("withPID"))){
+ // insert PID condition for Lambda instead of "kTRUE" value
+ // if ((fUsePID.Contains("withPID") && kTRUE )|| !(fUsePID.Contains("withPID"))){
+ if ((LambdaPID==1 && V0momentum <=1) || (V0momentum>1) || !(fUsePID.Contains("withPID"))){
- if (TMath::Abs(lRapLambda) < lCutRap) {
+ if (TMath::Abs(lRapLambda) < lCutRap) {
- //////2D histos: cut vs on fly status/////////////////////
+ //////2D histos: cut vs on fly status/////////////////////
- fHistDcaPosToPrimVertexL->Fill(lDcaPosToPrimVertex,lOnFlyStatus);
- fHistDcaNegToPrimVertexL->Fill(lDcaNegToPrimVertex,lOnFlyStatus);
- fHistRadiusV0L->Fill(lV0Radius,lOnFlyStatus);
- fHistDecayLengthV0L->Fill(lV0DecayLength,lOnFlyStatus);
- fHistDcaV0DaughtersL->Fill(lDcaV0Daughters,lOnFlyStatus);
- fHistChi2L->Fill(lChi2V0,lOnFlyStatus);
- fHistCosPointAngleL->Fill(lV0cosPointAngle,lOnFlyStatus);
+ fHistDcaPosToPrimVertexL->Fill(lDcaPosToPrimVertex,lOnFlyStatus);
+ fHistDcaNegToPrimVertexL->Fill(lDcaNegToPrimVertex,lOnFlyStatus);
+ fHistRadiusV0L->Fill(lV0Radius,lOnFlyStatus);
+ fHistDecayLengthV0L->Fill(lV0DecayLength,lOnFlyStatus);
+ fHistDcaV0DaughtersL->Fill(lDcaV0Daughters,lOnFlyStatus);
+ fHistChi2L->Fill(lChi2V0,lOnFlyStatus);
+ fHistCosPointAngleL->Fill(lV0cosPointAngle,lOnFlyStatus);
- //////2D histos: cut vs mass/////////////////////
+ //////2D histos: cut vs mass/////////////////////
- if (lOnFlyStatus==0){
- fHistMassLambda->Fill(lInvMassLambda);
- fHistMassVsRadiusLambda->Fill(rcPosRLambda,lInvMassLambda);
- fHistPtVsMassLambda->Fill(lInvMassLambda,lPtLambda);
+ if (lOnFlyStatus==0){
+ fHistMassLambda->Fill(lInvMassLambda);
+ fHistMassVsRadiusLambda->Fill(rcPosRLambda,lInvMassLambda);
+ fHistPtVsMassLambda->Fill(lInvMassLambda,lPtLambda);
+ if(V0momentum <=1) fHistNSigmaProton->Fill(nSigmaPosProton);
-
- fHistDcaPosToPrimVertexLvsMassL->Fill(lDcaPosToPrimVertex,lInvMassLambda);
- fHistDcaNegToPrimVertexLvsMassL->Fill(lDcaNegToPrimVertex,lInvMassLambda);
- fHistRadiusV0LvsMassL->Fill(lV0Radius,lInvMassLambda);
- fHistDecayLengthV0LvsMassL->Fill(lV0DecayLength,lInvMassLambda);
- fHistDcaV0DaughtersLvsMassL->Fill(lDcaV0Daughters,lInvMassLambda);
- fHistCosPointAngleLvsMassL->Fill(lV0cosPointAngle,lInvMassLambda);
+ fHistDcaPosToPrimVertexLvsMassL->Fill(lDcaPosToPrimVertex,lInvMassLambda);
+ fHistDcaNegToPrimVertexLvsMassL->Fill(lDcaNegToPrimVertex,lInvMassLambda);
+ fHistRadiusV0LvsMassL->Fill(lV0Radius,lInvMassLambda);
+ fHistDecayLengthV0LvsMassL->Fill(lV0DecayLength,lInvMassLambda);
+ fHistDcaV0DaughtersLvsMassL->Fill(lDcaV0Daughters,lInvMassLambda);
+ fHistCosPointAngleLvsMassL->Fill(lV0cosPointAngle,lInvMassLambda);
- if (lPtLambda>0 && lPtLambda <3){
- fHistDcaPosToPrimVertexLambdaVsMasspt1->Fill(lDcaPosToPrimVertex,lInvMassLambda);
- fHistDcaNegToPrimVertexLambdaVsMasspt1->Fill(lDcaNegToPrimVertex,lInvMassLambda);
- fHistRadiusV0LambdaVsMasspt1->Fill(lV0Radius,lInvMassLambda);
- fHistDecayLengthV0LambdaVsMasspt1->Fill(lV0DecayLength,lInvMassLambda);
- fHistDcaV0DaughtersLambdaVsMasspt1->Fill(lDcaV0Daughters,lInvMassLambda);
- fHistCosPointAngleLambdaVsMasspt1->Fill(lV0cosPointAngle,lInvMassLambda);
- }
- if (lPtLambda > 3 && lPtLambda < 6){
- fHistDcaPosToPrimVertexLambdaVsMasspt2->Fill(lDcaPosToPrimVertex,lInvMassLambda);
- fHistDcaNegToPrimVertexLambdaVsMasspt2->Fill(lDcaNegToPrimVertex,lInvMassLambda);
- fHistRadiusV0LambdaVsMasspt2->Fill(lV0Radius,lInvMassLambda);
- fHistDecayLengthV0LambdaVsMasspt2->Fill(lV0DecayLength,lInvMassLambda);
- fHistDcaV0DaughtersLambdaVsMasspt2->Fill(lDcaV0Daughters,lInvMassLambda);
- fHistCosPointAngleLambdaVsMasspt2->Fill(lV0cosPointAngle,lInvMassLambda);
- }
- if (lPtLambda > 6 && lPtLambda < 10){
- fHistDcaPosToPrimVertexLambdaVsMasspt3->Fill(lDcaPosToPrimVertex,lInvMassLambda);
- fHistDcaNegToPrimVertexLambdaVsMasspt3->Fill(lDcaNegToPrimVertex,lInvMassLambda);
- fHistRadiusV0LambdaVsMasspt3->Fill(lV0Radius,lInvMassLambda);
- fHistDecayLengthV0LambdaVsMasspt3->Fill(lV0DecayLength,lInvMassLambda);
- fHistDcaV0DaughtersLambdaVsMasspt3->Fill(lDcaV0Daughters,lInvMassLambda);
- fHistCosPointAngleLambdaVsMasspt3->Fill(lV0cosPointAngle,lInvMassLambda);
- }
- }
- } //end of Rap condition
- } // end of PID condition
-
- /////////// Anti Lambda ///////////////
- if ((AntiLambdaPID==1 && lPtAntiLambda <=1) || (lPtAntiLambda>1) || !(fUsePID.Contains("withPID"))){
- if (TMath::Abs(lRapAntiLambda) < lCutRap) {
-
- //////2D histos: cut vs on fly status/////////////////////
-
- fHistDcaPosToPrimVertexAntiL->Fill(lDcaPosToPrimVertex,lOnFlyStatus);
- fHistDcaNegToPrimVertexAntiL->Fill(lDcaNegToPrimVertex,lOnFlyStatus);
- fHistRadiusV0AntiL->Fill(lV0Radius,lOnFlyStatus);
- fHistDecayLengthV0AntiL->Fill(lV0DecayLength,lOnFlyStatus);
- fHistDcaV0DaughtersAntiL->Fill(lDcaV0Daughters,lOnFlyStatus);
- fHistChi2AntiL->Fill(lChi2V0,lOnFlyStatus);
- fHistCosPointAngleAntiL->Fill(lV0cosPointAngle,lOnFlyStatus);
-
- //////2D histos: cut vs mass/////////////////////
-
- if (lOnFlyStatus==0){
-
- fHistMassAntiLambda->Fill(lInvMassAntiLambda);
- fHistMassVsRadiusAntiLambda->Fill(rcPosRAntiLambda,lInvMassAntiLambda);
- fHistPtVsMassAntiLambda->Fill(lInvMassAntiLambda,lPtAntiLambda);
-
-
- fHistDcaPosToPrimVertexAntiLvsMass->Fill(lDcaPosToPrimVertex,lInvMassAntiLambda);
- fHistDcaNegToPrimVertexAntiLvsMass->Fill(lDcaNegToPrimVertex,lInvMassAntiLambda);
- fHistRadiusV0AntiLvsMass->Fill(lV0Radius,lInvMassAntiLambda);
- fHistDecayLengthV0AntiLvsMass->Fill(lV0DecayLength,lInvMassAntiLambda);
- fHistDcaV0DaughtersAntiLvsMass->Fill(lDcaV0Daughters,lInvMassAntiLambda);
- fHistCosPointAngleAntiLvsMass->Fill(lV0cosPointAngle,lInvMassAntiLambda);
+ if (lPtLambda>0 && lPtLambda <3){
+ fHistDcaPosToPrimVertexLambdaVsMasspt1->Fill(lDcaPosToPrimVertex,lInvMassLambda);
+ fHistDcaNegToPrimVertexLambdaVsMasspt1->Fill(lDcaNegToPrimVertex,lInvMassLambda);
+ fHistRadiusV0LambdaVsMasspt1->Fill(lV0Radius,lInvMassLambda);
+ fHistDecayLengthV0LambdaVsMasspt1->Fill(lV0DecayLength,lInvMassLambda);
+ fHistDcaV0DaughtersLambdaVsMasspt1->Fill(lDcaV0Daughters,lInvMassLambda);
+ fHistCosPointAngleLambdaVsMasspt1->Fill(lV0cosPointAngle,lInvMassLambda);
+ }
+ if (lPtLambda > 3 && lPtLambda < 6){
+ fHistDcaPosToPrimVertexLambdaVsMasspt2->Fill(lDcaPosToPrimVertex,lInvMassLambda);
+ fHistDcaNegToPrimVertexLambdaVsMasspt2->Fill(lDcaNegToPrimVertex,lInvMassLambda);
+ fHistRadiusV0LambdaVsMasspt2->Fill(lV0Radius,lInvMassLambda);
+ fHistDecayLengthV0LambdaVsMasspt2->Fill(lV0DecayLength,lInvMassLambda);
+ fHistDcaV0DaughtersLambdaVsMasspt2->Fill(lDcaV0Daughters,lInvMassLambda);
+ fHistCosPointAngleLambdaVsMasspt2->Fill(lV0cosPointAngle,lInvMassLambda);
+ }
+ if (lPtLambda > 6 && lPtLambda < 10){
+ fHistDcaPosToPrimVertexLambdaVsMasspt3->Fill(lDcaPosToPrimVertex,lInvMassLambda);
+ fHistDcaNegToPrimVertexLambdaVsMasspt3->Fill(lDcaNegToPrimVertex,lInvMassLambda);
+ fHistRadiusV0LambdaVsMasspt3->Fill(lV0Radius,lInvMassLambda);
+ fHistDecayLengthV0LambdaVsMasspt3->Fill(lV0DecayLength,lInvMassLambda);
+ fHistDcaV0DaughtersLambdaVsMasspt3->Fill(lDcaV0Daughters,lInvMassLambda);
+ fHistCosPointAngleLambdaVsMasspt3->Fill(lV0cosPointAngle,lInvMassLambda);
+ }
+ }
+ } //end of Rap condition
+ } // end of PID condition
+
+ /////////// Anti Lambda ///////////////
+ if ((AntiLambdaPID==1 && lPtAntiLambda <=1) || (lPtAntiLambda>1) || !(fUsePID.Contains("withPID"))){
+ if (TMath::Abs(lRapAntiLambda) < lCutRap) {
+
+ //////2D histos: cut vs on fly status/////////////////////
+
+ fHistDcaPosToPrimVertexAntiL->Fill(lDcaPosToPrimVertex,lOnFlyStatus);
+ fHistDcaNegToPrimVertexAntiL->Fill(lDcaNegToPrimVertex,lOnFlyStatus);
+ fHistRadiusV0AntiL->Fill(lV0Radius,lOnFlyStatus);
+ fHistDecayLengthV0AntiL->Fill(lV0DecayLength,lOnFlyStatus);
+ fHistDcaV0DaughtersAntiL->Fill(lDcaV0Daughters,lOnFlyStatus);
+ fHistChi2AntiL->Fill(lChi2V0,lOnFlyStatus);
+ fHistCosPointAngleAntiL->Fill(lV0cosPointAngle,lOnFlyStatus);
+
+ //////2D histos: cut vs mass/////////////////////
+
+ if (lOnFlyStatus==0){
+
+ fHistMassAntiLambda->Fill(lInvMassAntiLambda);
+ fHistMassVsRadiusAntiLambda->Fill(rcPosRAntiLambda,lInvMassAntiLambda);
+ fHistPtVsMassAntiLambda->Fill(lInvMassAntiLambda,lPtAntiLambda);
+
+
+ fHistDcaPosToPrimVertexAntiLvsMass->Fill(lDcaPosToPrimVertex,lInvMassAntiLambda);
+ fHistDcaNegToPrimVertexAntiLvsMass->Fill(lDcaNegToPrimVertex,lInvMassAntiLambda);
+ fHistRadiusV0AntiLvsMass->Fill(lV0Radius,lInvMassAntiLambda);
+ fHistDecayLengthV0AntiLvsMass->Fill(lV0DecayLength,lInvMassAntiLambda);
+ fHistDcaV0DaughtersAntiLvsMass->Fill(lDcaV0Daughters,lInvMassAntiLambda);
+ fHistCosPointAngleAntiLvsMass->Fill(lV0cosPointAngle,lInvMassAntiLambda);
- if (lPtAntiLambda>0 && lPtAntiLambda <3){
- fHistDcaPosToPrimVertexAntiLVsMasspt1->Fill(lDcaPosToPrimVertex,lInvMassAntiLambda);
- fHistDcaNegToPrimVertexAntiLVsMasspt1->Fill(lDcaNegToPrimVertex,lInvMassAntiLambda);
- fHistRadiusV0AntiLVsMasspt1->Fill(lV0Radius,lInvMassAntiLambda);
- fHistDecayLengthV0AntiLVsMasspt1->Fill(lV0DecayLength,lInvMassAntiLambda);
- fHistDcaV0DaughtersAntiLVsMasspt1->Fill(lDcaV0Daughters,lInvMassAntiLambda);
- fHistCosPointAngleAntiLVsMasspt1->Fill(lV0cosPointAngle,lInvMassAntiLambda);
- }
- if (lPtAntiLambda > 3 && lPtAntiLambda < 6){
- fHistDcaPosToPrimVertexAntiLVsMasspt2->Fill(lDcaPosToPrimVertex,lInvMassAntiLambda);
- fHistDcaNegToPrimVertexAntiLVsMasspt2->Fill(lDcaNegToPrimVertex,lInvMassAntiLambda);
- fHistRadiusV0AntiLVsMasspt2->Fill(lV0Radius,lInvMassAntiLambda);
- fHistDecayLengthV0AntiLVsMasspt2->Fill(lV0DecayLength,lInvMassAntiLambda);
- fHistDcaV0DaughtersAntiLVsMasspt2->Fill(lDcaV0Daughters,lInvMassAntiLambda);
- fHistCosPointAngleAntiLVsMasspt2->Fill(lV0cosPointAngle,lInvMassAntiLambda);
- }
- if (lPtAntiLambda > 6 && lPtAntiLambda < 10){
- fHistDcaPosToPrimVertexAntiLVsMasspt3->Fill(lDcaPosToPrimVertex,lInvMassAntiLambda);
- fHistDcaNegToPrimVertexAntiLVsMasspt3->Fill(lDcaNegToPrimVertex,lInvMassAntiLambda);
- fHistRadiusV0AntiLVsMasspt3->Fill(lV0Radius,lInvMassAntiLambda);
- fHistDecayLengthV0AntiLVsMasspt3->Fill(lV0DecayLength,lInvMassAntiLambda);
- fHistDcaV0DaughtersAntiLVsMasspt3->Fill(lDcaV0Daughters,lInvMassAntiLambda);
- fHistCosPointAngleAntiLVsMasspt3->Fill(lV0cosPointAngle,lInvMassAntiLambda);
- }
- }
- } //end of Rap condition
- } // end of PID condition
-
-
- ///////////////values for cuts end////////////////////////////////////////////////////////////////////////
-
-
- // Histo versus Rap and armenteros plot
- if (!lOnFlyStatus){
- if (lCheckMcK0Short) fHistAsMcRapK0->Fill(lRapK0s);
- if (lCheckMcLambda) fHistAsMcRapLambda->Fill(lRapLambda);
- if (lCheckMcAntiLambda) fHistAsMcRapLambda->Fill(lRapAntiLambda);
- // fHistArmenterosPodolanski->Fill(lAlphaV0,lPtArmV0);
- if ((TMath::Abs(lRapK0s) < lCutRap)&&(TMath::Abs(lRapLambda) < lCutRap)) fHistK0sMassVsLambdaMass->Fill(lInvMassK0s,lInvMassLambda);
- }
+ if (lPtAntiLambda>0 && lPtAntiLambda <3){
+ fHistDcaPosToPrimVertexAntiLVsMasspt1->Fill(lDcaPosToPrimVertex,lInvMassAntiLambda);
+ fHistDcaNegToPrimVertexAntiLVsMasspt1->Fill(lDcaNegToPrimVertex,lInvMassAntiLambda);
+ fHistRadiusV0AntiLVsMasspt1->Fill(lV0Radius,lInvMassAntiLambda);
+ fHistDecayLengthV0AntiLVsMasspt1->Fill(lV0DecayLength,lInvMassAntiLambda);
+ fHistDcaV0DaughtersAntiLVsMasspt1->Fill(lDcaV0Daughters,lInvMassAntiLambda);
+ fHistCosPointAngleAntiLVsMasspt1->Fill(lV0cosPointAngle,lInvMassAntiLambda);
+ }
+ if (lPtAntiLambda > 3 && lPtAntiLambda < 6){
+ fHistDcaPosToPrimVertexAntiLVsMasspt2->Fill(lDcaPosToPrimVertex,lInvMassAntiLambda);
+ fHistDcaNegToPrimVertexAntiLVsMasspt2->Fill(lDcaNegToPrimVertex,lInvMassAntiLambda);
+ fHistRadiusV0AntiLVsMasspt2->Fill(lV0Radius,lInvMassAntiLambda);
+ fHistDecayLengthV0AntiLVsMasspt2->Fill(lV0DecayLength,lInvMassAntiLambda);
+ fHistDcaV0DaughtersAntiLVsMasspt2->Fill(lDcaV0Daughters,lInvMassAntiLambda);
+ fHistCosPointAngleAntiLVsMasspt2->Fill(lV0cosPointAngle,lInvMassAntiLambda);
+ }
+ if (lPtAntiLambda > 6 && lPtAntiLambda < 10){
+ fHistDcaPosToPrimVertexAntiLVsMasspt3->Fill(lDcaPosToPrimVertex,lInvMassAntiLambda);
+ fHistDcaNegToPrimVertexAntiLVsMasspt3->Fill(lDcaNegToPrimVertex,lInvMassAntiLambda);
+ fHistRadiusV0AntiLVsMasspt3->Fill(lV0Radius,lInvMassAntiLambda);
+ fHistDecayLengthV0AntiLVsMasspt3->Fill(lV0DecayLength,lInvMassAntiLambda);
+ fHistDcaV0DaughtersAntiLVsMasspt3->Fill(lDcaV0Daughters,lInvMassAntiLambda);
+ fHistCosPointAngleAntiLVsMasspt3->Fill(lV0cosPointAngle,lInvMassAntiLambda);
+ }
+ }
+ } //end of Rap condition
+ } // end of PID condition
+
+
+ ///////////////values for cuts end////////////////////////////////////////////////////////////////////////
+
+
+ // Histo versus Rap and armenteros plot
+ if (!lOnFlyStatus){
+ if (lCheckMcK0Short) fHistAsMcRapK0->Fill(lRapK0s);
+ if (lCheckMcLambda) fHistAsMcRapLambda->Fill(lRapLambda);
+ if (lCheckMcAntiLambda) fHistAsMcRapLambda->Fill(lRapAntiLambda);
+ // fHistArmenterosPodolanski->Fill(lAlphaV0,lPtArmV0);
+ if ((TMath::Abs(lRapK0s) < lCutRap)&&(TMath::Abs(lRapLambda) < lCutRap)) fHistK0sMassVsLambdaMass->Fill(lInvMassK0s,lInvMassLambda);
+ }
- // K0s associated histograms in |rap| < lCutRap:
+ // K0s associated histograms in |rap| < lCutRap:
- if (TMath::Abs(lRapK0s) < lCutRap) {
+ if (TMath::Abs(lRapK0s) < lCutRap) {
- fHistNsigmaPosPionK0->Fill(nSigmaPosPion);
- fHistNsigmaNegPionK0->Fill(nSigmaNegPion);
+ fHistNsigmaPosPionK0->Fill(nSigmaPosPion);
+ fHistNsigmaNegPionK0->Fill(nSigmaNegPion);
- switch (lOnFlyStatus){
- case 0 :
-
-
- // fHistMultVsPtVsMassK0->Fill(multiplicity ,lInvMassK0s,lPtK0s);
- if(lCheckPIdK0Short) fHistPidMcMassK0->Fill(lInvMassK0s);
- if(lCheckMcK0Short) {
- fHistAsMcMassK0->Fill(lInvMassK0s);
- fHistAsMcPtK0->Fill(lPtK0s);
-
-
- fHistAsMcPtVsMassK0->Fill(lInvMassK0s,lPtK0s);
- if (lPtK0s <= 1) fHistAsMcPtZoomK0->Fill(lPtK0s);
- fHistAsMcMassVsRadiusK0->Fill(rcPosRK0s,lInvMassK0s);
- fHistAsMcResxK0->Fill(rcPosXK0s-mcPosX);
- fHistAsMcResyK0->Fill(rcPosYK0s-mcPosY);
- fHistAsMcReszK0->Fill(rcPosZK0s-mcPosZ);
- fHistAsMcResrVsRadiusK0->Fill(rcPosRK0s,rcPosRK0s-mcPosR);
- fHistAsMcReszVsRadiusK0->Fill(rcPosZK0s,rcPosZK0s-mcPosZ);
- fHistAsMcProdRadiusK0->Fill(mcPosMotherR);
- fHistAsMcProdRadiusXvsYK0s->Fill(mcPosMotherX,mcPosMotherY);
- fHistAsMcResPtK0->Fill(deltaPtK0s);
- fHistAsMcResPtVsRapK0->Fill(deltaPtK0s,lRapK0s);
- fHistAsMcResPtVsPtK0->Fill(deltaPtK0s,lPtK0s);
- }
- else if (lCheckSecondaryK0s) {
- fHistAsMcSecondaryPtVsRapK0s->Fill(lPtK0s,lRapK0s);
- fHistAsMcSecondaryProdRadiusK0s->Fill(mcPosMotherR);
- fHistAsMcSecondaryProdRadiusXvsYK0s->Fill(mcPosMotherX,mcPosMotherY);
- switch (lPdgcodeMotherOfMother) {
- case 130 : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(0.5);break; // K0L
- case 321 : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(1.5);break; // K+
- case -321 : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(2.5);break; // K-
- case -3122 : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(3.5);break; //AntiLambda
- default : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(6.5);break;
- }
- }
- break;
+ switch (lOnFlyStatus){
+ case 0 :
+
+
+ // fHistMultVsPtVsMassK0->Fill(multiplicity ,lInvMassK0s,lPtK0s);
+ if(lCheckPIdK0Short) fHistPidMcMassK0->Fill(lInvMassK0s);
+ if(lCheckMcK0Short) {
+ fHistAsMcMassK0->Fill(lInvMassK0s);
+ fHistAsMcPtK0->Fill(lPtK0s);
+
+
+ fHistAsMcPtVsMassK0->Fill(lInvMassK0s,lPtK0s);
+ if (lPtK0s <= 1) fHistAsMcPtZoomK0->Fill(lPtK0s);
+ fHistAsMcMassVsRadiusK0->Fill(rcPosRK0s,lInvMassK0s);
+ fHistAsMcResxK0->Fill(rcPosXK0s-mcPosX);
+ fHistAsMcResyK0->Fill(rcPosYK0s-mcPosY);
+ fHistAsMcReszK0->Fill(rcPosZK0s-mcPosZ);
+ fHistAsMcResrVsRadiusK0->Fill(rcPosRK0s,rcPosRK0s-mcPosR);
+ fHistAsMcReszVsRadiusK0->Fill(rcPosZK0s,rcPosZK0s-mcPosZ);
+ fHistAsMcProdRadiusK0->Fill(mcPosMotherR);
+ fHistAsMcProdRadiusXvsYK0s->Fill(mcPosMotherX,mcPosMotherY);
+ fHistAsMcResPtK0->Fill(deltaPtK0s);
+ fHistAsMcResPtVsRapK0->Fill(deltaPtK0s,lRapK0s);
+ fHistAsMcResPtVsPtK0->Fill(deltaPtK0s,lPtK0s);
+ }
+ else if (lCheckSecondaryK0s) {
+ fHistAsMcSecondaryPtVsRapK0s->Fill(lPtK0s,lRapK0s);
+ fHistAsMcSecondaryProdRadiusK0s->Fill(mcPosMotherR);
+ fHistAsMcSecondaryProdRadiusXvsYK0s->Fill(mcPosMotherX,mcPosMotherY);
+ switch (lPdgcodeMotherOfMother) {
+ case 130 : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(0.5);break; // K0L
+ case 321 : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(1.5);break; // K+
+ case -321 : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(2.5);break; // K-
+ case -3122 : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(3.5);break; //AntiLambda
+ default : fHistAsMcSecondaryMotherPdgCodeK0s->Fill(6.5);break;
+ }
+ }
+ break;
- }
- } // end rapidity condition
+ }
+ } // end rapidity condition
- // Associated Lambda histograms in |rap| < lCutRap
- if ((LambdaPID==1 && lPtLambda <=1) || (lPtLambda>1) || !(fUsePID.Contains("withPID"))){
+ // Associated Lambda histograms in |rap| < lCutRap
+ if ((LambdaPID==1 && V0momentum <=1) || (V0momentum>1) || !(fUsePID.Contains("withPID"))){
- if (TMath::Abs(lRapLambda) < lCutRap) {
+ if (TMath::Abs(lRapLambda) < lCutRap) {
- fHistNsigmaPosProtonLambda->Fill(nSigmaPosProton);
- fHistNsigmaNegPionLambda->Fill(nSigmaNegPion);
- switch (lOnFlyStatus){
- case 0 :
+ fHistNsigmaPosProtonLambda->Fill(nSigmaPosProton);
+ fHistNsigmaNegPionLambda->Fill(nSigmaNegPion);
+ switch (lOnFlyStatus){
+ case 0 :
- // fHistMultVsPtVsMassLambda->Fill(multiplicity ,lInvMassLambda,lPtLambda);
- if(lCheckPIdLambda) fHistPidMcMassLambda->Fill(lInvMassLambda);
+ // fHistMultVsPtVsMassLambda->Fill(multiplicity ,lInvMassLambda,lPtLambda);
+ if(lCheckPIdLambda) fHistPidMcMassLambda->Fill(lInvMassLambda);
- if(lCheckMcLambda) {
- fHistAsMcMassLambda->Fill(lInvMassLambda);
- fHistAsMcPtLambda->Fill(lPtLambda);
-
- fHistAsMcPtVsMassLambda->Fill(lInvMassLambda,lPtLambda);
- if (lPtLambda <= 1) fHistAsMcPtZoomLambda->Fill(lPtLambda);
- fHistAsMcMassVsRadiusLambda->Fill(rcPosRLambda,lInvMassLambda);
- fHistAsMcResxLambda->Fill(rcPosXLambda-mcPosX);
- fHistAsMcResyLambda->Fill(rcPosYLambda-mcPosY);
- fHistAsMcReszLambda->Fill(rcPosZLambda-mcPosZ);
- fHistAsMcResrVsRadiusLambda->Fill(rcPosRLambda,rcPosRLambda-mcPosR);
- fHistAsMcReszVsRadiusLambda->Fill(rcPosZLambda,rcPosZLambda-mcPosZ);
- fHistAsMcProdRadiusLambda->Fill(mcPosMotherR);
- fHistAsMcProdRadiusXvsYLambda->Fill(mcPosMotherX,mcPosMotherY);
- fHistAsMcResPtLambda->Fill(deltaPtLambda);
- fHistAsMcResPtVsRapLambda->Fill(deltaPtLambda,lRapLambda);
- fHistAsMcResPtVsPtLambda->Fill(deltaPtLambda,lPtLambda);
- if (lComeFromSigma) fHistAsMcPtLambdaFromSigma->Fill(lPtLambda);
- switch (lPdgcodeMotherOfMother) {
- case 3222 : fHistAsMcMotherPdgCodeLambda->Fill(0.5); break; // Sigma +
- case 3212 : fHistAsMcMotherPdgCodeLambda->Fill(1.5); break; // Sigma 0
- case 3112 : fHistAsMcMotherPdgCodeLambda->Fill(2.5); break;// Sigma -
- case 3224 : fHistAsMcMotherPdgCodeLambda->Fill(3.5); break;// Sigma(1385) +
- case 3214 : fHistAsMcMotherPdgCodeLambda->Fill(4.5); break;// Sigma(1385) 0
- case 3114 : fHistAsMcMotherPdgCodeLambda->Fill(5.5); break;// Sigma(1385) -
- case 3322 : fHistAsMcMotherPdgCodeLambda->Fill(6.5); break; // Xi 0
- case 3312 : fHistAsMcMotherPdgCodeLambda->Fill(7.5); break; // Xi -
- case 3334 : fHistAsMcMotherPdgCodeLambda->Fill(8.5); break; // Omega
- case -1 : fHistAsMcMotherPdgCodeLambda->Fill(9.5); break;
- default : fHistAsMcMotherPdgCodeLambda->Fill(10.5);break;
- }
+ if(lCheckMcLambda) {
+ fHistAsMcMassLambda->Fill(lInvMassLambda);
+ fHistAsMcPtLambda->Fill(lPtLambda);
+
+ fHistAsMcPtVsMassLambda->Fill(lInvMassLambda,lPtLambda);
+ if (lPtLambda <= 1) fHistAsMcPtZoomLambda->Fill(lPtLambda);
+ fHistAsMcMassVsRadiusLambda->Fill(rcPosRLambda,lInvMassLambda);
+ fHistAsMcResxLambda->Fill(rcPosXLambda-mcPosX);
+ fHistAsMcResyLambda->Fill(rcPosYLambda-mcPosY);
+ fHistAsMcReszLambda->Fill(rcPosZLambda-mcPosZ);
+ fHistAsMcResrVsRadiusLambda->Fill(rcPosRLambda,rcPosRLambda-mcPosR);
+ fHistAsMcReszVsRadiusLambda->Fill(rcPosZLambda,rcPosZLambda-mcPosZ);
+ fHistAsMcProdRadiusLambda->Fill(mcPosMotherR);
+ fHistAsMcProdRadiusXvsYLambda->Fill(mcPosMotherX,mcPosMotherY);
+ fHistAsMcResPtLambda->Fill(deltaPtLambda);
+ fHistAsMcResPtVsRapLambda->Fill(deltaPtLambda,lRapLambda);
+ fHistAsMcResPtVsPtLambda->Fill(deltaPtLambda,lPtLambda);
+ if (lComeFromSigma) fHistAsMcPtLambdaFromSigma->Fill(lPtLambda);
+ switch (lPdgcodeMotherOfMother) {
+ case 3222 : fHistAsMcMotherPdgCodeLambda->Fill(0.5); break; // Sigma +
+ case 3212 : fHistAsMcMotherPdgCodeLambda->Fill(1.5); break; // Sigma 0
+ case 3112 : fHistAsMcMotherPdgCodeLambda->Fill(2.5); break;// Sigma -
+ case 3224 : fHistAsMcMotherPdgCodeLambda->Fill(3.5); break;// Sigma(1385) +
+ case 3214 : fHistAsMcMotherPdgCodeLambda->Fill(4.5); break;// Sigma(1385) 0
+ case 3114 : fHistAsMcMotherPdgCodeLambda->Fill(5.5); break;// Sigma(1385) -
+ case 3322 : fHistAsMcMotherPdgCodeLambda->Fill(6.5); break; // Xi 0
+ case 3312 : fHistAsMcMotherPdgCodeLambda->Fill(7.5); break; // Xi -
+ case 3334 : fHistAsMcMotherPdgCodeLambda->Fill(8.5); break; // Omega
+ case -1 : fHistAsMcMotherPdgCodeLambda->Fill(9.5); break;
+ default : fHistAsMcMotherPdgCodeLambda->Fill(10.5);break;
+ }
- //printf("found Lambda RC dcaPos=%e dcaNeg=%e dcaDau=%e cosP=%e pT=%e mass=%e\n",lDcaPosToPrimVertex ,lDcaNegToPrimVertex ,lDcaV0Daughters,lV0cosPointAngle,lPtLambda,lInvMassLambda);
- //printf("found Lambda RC Pindex=%d Nindex=%d Plabel=%d Nlabel=%d\n\n",lIndexTrackPos,lIndexTrackNeg,lLabelTrackPos,lLabelTrackNeg);
+ //printf("found Lambda RC dcaPos=%e dcaNeg=%e dcaDau=%e cosP=%e pT=%e mass=%e\n",lDcaPosToPrimVertex ,lDcaNegToPrimVertex ,lDcaV0Daughters,lV0cosPointAngle,lPtLambda,lInvMassLambda);
+ //printf("found Lambda RC Pindex=%d Nindex=%d Plabel=%d Nlabel=%d\n\n",lIndexTrackPos,lIndexTrackNeg,lLabelTrackPos,lLabelTrackNeg);
- }
+ }
- else if (lCheckSecondaryLambda) {
- fHistAsMcSecondaryPtVsRapLambda->Fill(lPtLambda,lRapLambda);
- fHistAsMcSecondaryProdRadiusLambda->Fill(mcPosMotherR);
- fHistAsMcSecondaryProdRadiusXvsYLambda->Fill(mcPosMotherX,mcPosMotherY);
- if (lComeFromSigma) fHistAsMcSecondaryPtLambdaFromSigma->Fill(lPtLambda);
- printf(" lPdgcodeMotherOfMother= %d",lPdgcodeMotherOfMother);
- switch (lPdgcodeMotherOfMother) {
- case 3222 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(0.5); break;// Sigma +
- case 3212 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(1.5); break;// Sigma 0
- case 3112 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(2.5); break;// Sigma -
- case 3224 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(3.5); break;// Sigma(1385) +
- case 3214 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(4.5); break;// Sigma(1385) 0
- case 3114 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(5.5); break;// Sigma(1385) -
- case 3322 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(6.5); break; // Xi 0
- case 3312 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(7.5); break; // Xi -
- case 3334 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(8.5); break; // Omega
- case -1 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(9.5); break;
- default : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(10.5);break;
- }
- }
- break;
+ else if (lCheckSecondaryLambda) {
+ fHistAsMcSecondaryPtVsRapLambda->Fill(lPtLambda,lRapLambda);
+ fHistAsMcSecondaryProdRadiusLambda->Fill(mcPosMotherR);
+ fHistAsMcSecondaryProdRadiusXvsYLambda->Fill(mcPosMotherX,mcPosMotherY);
+ if (lComeFromSigma) fHistAsMcSecondaryPtLambdaFromSigma->Fill(lPtLambda);
+ printf(" lPdgcodeMotherOfMother= %d",lPdgcodeMotherOfMother);
+ switch (lPdgcodeMotherOfMother) {
+ case 3222 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(0.5); break;// Sigma +
+ case 3212 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(1.5); break;// Sigma 0
+ case 3112 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(2.5); break;// Sigma -
+ case 3224 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(3.5); break;// Sigma(1385) +
+ case 3214 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(4.5); break;// Sigma(1385) 0
+ case 3114 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(5.5); break;// Sigma(1385) -
+ case 3322 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(6.5); break; // Xi 0
+ case 3312 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(7.5); break; // Xi -
+ case 3334 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(8.5); break; // Omega
+ case -1 : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(9.5); break;
+ default : fHistAsMcSecondaryMotherPdgCodeLambda->Fill(10.5);break;
+ }
+ }
+ break;
- }
- } // end rapidity condition
- }// end PID condition
- // Associated AntiLambda histograms in |rap| < lCutRap
- if ((AntiLambdaPID==1 && lPtAntiLambda <=1) || (lPtAntiLambda>1) || !(fUsePID.Contains("withPID"))){
+ }
+ } // end rapidity condition
+ }// end PID condition
+ // Associated AntiLambda histograms in |rap| < lCutRap
+ if ((AntiLambdaPID==1 && lPtAntiLambda <=1) || (lPtAntiLambda>1) || !(fUsePID.Contains("withPID"))){
- if (TMath::Abs(lRapAntiLambda) < lCutRap) {
+ if (TMath::Abs(lRapAntiLambda) < lCutRap) {
- fHistNsigmaPosProtonAntiLambda->Fill(nSigmaPosProton);
- fHistNsigmaNegPionAntiLambda->Fill(nSigmaNegPion);
- switch (lOnFlyStatus){
- case 0 :
+ fHistNsigmaPosProtonAntiLambda->Fill(nSigmaPosProton);
+ fHistNsigmaNegPionAntiLambda->Fill(nSigmaNegPion);
+ switch (lOnFlyStatus){
+ case 0 :
- // fHistMultVsPtVsMassAntiLambda->Fill(multiplicity ,lInvMassAntiLambda,lPtAntiLambda);
- if(lCheckPIdAntiLambda) fHistPidMcMassAntiLambda->Fill(lInvMassAntiLambda);
+ // fHistMultVsPtVsMassAntiLambda->Fill(multiplicity ,lInvMassAntiLambda,lPtAntiLambda);
+ if(lCheckPIdAntiLambda) fHistPidMcMassAntiLambda->Fill(lInvMassAntiLambda);
- if(lCheckMcAntiLambda) {
- fHistAsMcMassAntiLambda->Fill(lInvMassAntiLambda);
- fHistAsMcPtAntiLambda->Fill(lPtAntiLambda);
-
-
- fHistAsMcPtVsMassAntiLambda->Fill(lInvMassAntiLambda,lPtAntiLambda);
- if (lPtAntiLambda <= 1) fHistAsMcPtZoomAntiLambda->Fill(lPtAntiLambda);
- fHistAsMcMassVsRadiusAntiLambda->Fill(rcPosRAntiLambda,lInvMassAntiLambda);
- fHistAsMcResxAntiLambda->Fill(rcPosXAntiLambda-mcPosX);
- fHistAsMcResyAntiLambda->Fill(rcPosYAntiLambda-mcPosY);
- fHistAsMcReszAntiLambda->Fill(rcPosZAntiLambda-mcPosZ);
- fHistAsMcResrVsRadiusAntiLambda->Fill(rcPosRAntiLambda,rcPosRAntiLambda-mcPosR);
- fHistAsMcReszVsRadiusAntiLambda->Fill(rcPosZAntiLambda,rcPosZAntiLambda-mcPosZ);
- fHistAsMcProdRadiusAntiLambda->Fill(mcPosMotherR);
- fHistAsMcProdRadiusXvsYAntiLambda->Fill(mcPosMotherX,mcPosMotherY);
- fHistAsMcResPtAntiLambda->Fill(deltaPtAntiLambda);
- fHistAsMcResPtVsRapAntiLambda->Fill(deltaPtAntiLambda,lRapAntiLambda);
- fHistAsMcResPtVsPtAntiLambda->Fill(deltaPtAntiLambda,lPtAntiLambda);
- if (lComeFromSigma) fHistAsMcPtAntiLambdaFromSigma->Fill(lPtAntiLambda);
- switch (lPdgcodeMotherOfMother) {
- case 3222 : fHistAsMcMotherPdgCodeAntiLambda->Fill(0.5); break; // Sigma +
- case 3212 : fHistAsMcMotherPdgCodeAntiLambda->Fill(1.5); break; // Sigma 0
- case 3112 : fHistAsMcMotherPdgCodeAntiLambda->Fill(2.5); break;// Sigma -
- case 3224 : fHistAsMcMotherPdgCodeAntiLambda->Fill(3.5); break;// Sigma(1385) +
- case 3214 : fHistAsMcMotherPdgCodeAntiLambda->Fill(4.5); break;// Sigma(1385) 0
- case 3114 : fHistAsMcMotherPdgCodeAntiLambda->Fill(5.5); break;// Sigma(1385) -
- case 3322 : fHistAsMcMotherPdgCodeAntiLambda->Fill(6.5); break; // Xi 0
- case 3312 : fHistAsMcMotherPdgCodeAntiLambda->Fill(7.5); break; // Xi -
- case 3334 : fHistAsMcMotherPdgCodeAntiLambda->Fill(8.5); break; // Omega
-
- case -1 : fHistAsMcMotherPdgCodeAntiLambda->Fill(9.5); break;
- default : fHistAsMcMotherPdgCodeAntiLambda->Fill(10.5);break;
- }
+ if(lCheckMcAntiLambda) {
+ fHistAsMcMassAntiLambda->Fill(lInvMassAntiLambda);
+ fHistAsMcPtAntiLambda->Fill(lPtAntiLambda);
+
+
+ fHistAsMcPtVsMassAntiLambda->Fill(lInvMassAntiLambda,lPtAntiLambda);
+ if (lPtAntiLambda <= 1) fHistAsMcPtZoomAntiLambda->Fill(lPtAntiLambda);
+ fHistAsMcMassVsRadiusAntiLambda->Fill(rcPosRAntiLambda,lInvMassAntiLambda);
+ fHistAsMcResxAntiLambda->Fill(rcPosXAntiLambda-mcPosX);
+ fHistAsMcResyAntiLambda->Fill(rcPosYAntiLambda-mcPosY);
+ fHistAsMcReszAntiLambda->Fill(rcPosZAntiLambda-mcPosZ);
+ fHistAsMcResrVsRadiusAntiLambda->Fill(rcPosRAntiLambda,rcPosRAntiLambda-mcPosR);
+ fHistAsMcReszVsRadiusAntiLambda->Fill(rcPosZAntiLambda,rcPosZAntiLambda-mcPosZ);
+ fHistAsMcProdRadiusAntiLambda->Fill(mcPosMotherR);
+ fHistAsMcProdRadiusXvsYAntiLambda->Fill(mcPosMotherX,mcPosMotherY);
+ fHistAsMcResPtAntiLambda->Fill(deltaPtAntiLambda);
+ fHistAsMcResPtVsRapAntiLambda->Fill(deltaPtAntiLambda,lRapAntiLambda);
+ fHistAsMcResPtVsPtAntiLambda->Fill(deltaPtAntiLambda,lPtAntiLambda);
+ if (lComeFromSigma) fHistAsMcPtAntiLambdaFromSigma->Fill(lPtAntiLambda);
+ switch (lPdgcodeMotherOfMother) {
+ case 3222 : fHistAsMcMotherPdgCodeAntiLambda->Fill(0.5); break; // Sigma +
+ case 3212 : fHistAsMcMotherPdgCodeAntiLambda->Fill(1.5); break; // Sigma 0
+ case 3112 : fHistAsMcMotherPdgCodeAntiLambda->Fill(2.5); break;// Sigma -
+ case 3224 : fHistAsMcMotherPdgCodeAntiLambda->Fill(3.5); break;// Sigma(1385) +
+ case 3214 : fHistAsMcMotherPdgCodeAntiLambda->Fill(4.5); break;// Sigma(1385) 0
+ case 3114 : fHistAsMcMotherPdgCodeAntiLambda->Fill(5.5); break;// Sigma(1385) -
+ case 3322 : fHistAsMcMotherPdgCodeAntiLambda->Fill(6.5); break; // Xi 0
+ case 3312 : fHistAsMcMotherPdgCodeAntiLambda->Fill(7.5); break; // Xi -
+ case 3334 : fHistAsMcMotherPdgCodeAntiLambda->Fill(8.5); break; // Omega
+
+ case -1 : fHistAsMcMotherPdgCodeAntiLambda->Fill(9.5); break;
+ default : fHistAsMcMotherPdgCodeAntiLambda->Fill(10.5);break;
+ }
- //printf("found AntiLambda RC dcaPos=%e dcaNeg=%e dcaDau=%e cosP=%e pT=%e mass=%e\n",lDcaPosToPrimVertex ,lDcaNegToPrimVertex ,lDcaV0Daughters,lV0cosPointAngle,lPtAntiLambda,lInvMassAntiLambda);
- //printf("found AntiLambda RC Pindex=%d Nindex=%d Plabel=%d Nlabel=%d\n\n",lIndexTrackPos,lIndexTrackNeg,lLabelTrackPos,lLabelTrackNeg);
+ //printf("found AntiLambda RC dcaPos=%e dcaNeg=%e dcaDau=%e cosP=%e pT=%e mass=%e\n",lDcaPosToPrimVertex ,lDcaNegToPrimVertex ,lDcaV0Daughters,lV0cosPointAngle,lPtAntiLambda,lInvMassAntiLambda);
+ //printf("found AntiLambda RC Pindex=%d Nindex=%d Plabel=%d Nlabel=%d\n\n",lIndexTrackPos,lIndexTrackNeg,lLabelTrackPos,lLabelTrackNeg);
- }
+ }
- else if (lCheckSecondaryAntiLambda) {
- fHistAsMcSecondaryPtVsRapAntiLambda->Fill(lPtAntiLambda,lRapAntiLambda);
- fHistAsMcSecondaryProdRadiusAntiLambda->Fill(mcPosMotherR);
- fHistAsMcSecondaryProdRadiusXvsYAntiLambda->Fill(mcPosMotherX,mcPosMotherY);
- if (lComeFromSigma) fHistAsMcSecondaryPtAntiLambdaFromSigma->Fill(lPtAntiLambda);
- printf(" lPdgcodeMotherOfMother= %d",lPdgcodeMotherOfMother);
- switch (lPdgcodeMotherOfMother) {
- case 3222 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(0.5); break;// Sigma +
- case 3212 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(1.5); break;// Sigma 0
- case 3112 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(2.5); break;// Sigma -
- case 3224 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(3.5); break;// Sigma(1385) +
- case 3214 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(4.5); break;// Sigma(1385) 0
- case 3114 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(5.5); break;// Sigma(1385) -
- case 3322 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(6.5); break; // Xi 0
- case 3312 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(7.5); break; // Xi -
- case 3334 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(8.5); break; // Omega
- case -1 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(9.5); break;
- default : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(10.5);break;
- }
- }
- break;
+ else if (lCheckSecondaryAntiLambda) {
+ fHistAsMcSecondaryPtVsRapAntiLambda->Fill(lPtAntiLambda,lRapAntiLambda);
+ fHistAsMcSecondaryProdRadiusAntiLambda->Fill(mcPosMotherR);
+ fHistAsMcSecondaryProdRadiusXvsYAntiLambda->Fill(mcPosMotherX,mcPosMotherY);
+ if (lComeFromSigma) fHistAsMcSecondaryPtAntiLambdaFromSigma->Fill(lPtAntiLambda);
+ printf(" lPdgcodeMotherOfMother= %d",lPdgcodeMotherOfMother);
+ switch (lPdgcodeMotherOfMother) {
+ case 3222 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(0.5); break;// Sigma +
+ case 3212 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(1.5); break;// Sigma 0
+ case 3112 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(2.5); break;// Sigma -
+ case 3224 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(3.5); break;// Sigma(1385) +
+ case 3214 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(4.5); break;// Sigma(1385) 0
+ case 3114 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(5.5); break;// Sigma(1385) -
+ case 3322 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(6.5); break; // Xi 0
+ case 3312 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(7.5); break; // Xi -
+ case 3334 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(8.5); break; // Omega
+ case -1 : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(9.5); break;
+ default : fHistAsMcSecondaryMotherPdgCodeAntiLambda->Fill(10.5);break;
+ }
+ }
+ break;
- }
- } // end rapidity condition
- }// end PID condition
-
- if (negPiKF) delete negPiKF; negPiKF= NULL;
- if (posPiKF) delete posPiKF; posPiKF= NULL;
- if (posPKF) delete posPKF; posPKF = NULL;
- if (negAPKF) delete negAPKF; negAPKF= NULL;
+ }
+ } // end rapidity condition
+ }// end PID condition
+
+ if (negPiKF) delete negPiKF; negPiKF= NULL;
+ if (posPiKF) delete posPiKF; posPiKF= NULL;
+ if (posPKF) delete posPKF; posPKF = NULL;
+ if (negAPKF) delete negAPKF; negAPKF= NULL;
} // end V0 loop