fESDTrackCutsLowPt(0),
fOutputContainer(0),
//thnf
- fTHnFK0s(0),
- fTHnFL(0),
- fTHnFAL(0),
+// fTHnFK0s(0),
+// fTHnFL(0),
+// fTHnFAL(0),
fTHnFK0sDauEta(0),
fTHnFLDauEta(0),
fTHnFALDauEta(0),
fESDTrackCutsLowPt(0),
fOutputContainer(0),
//thnf
- fTHnFK0s(0),
- fTHnFL(0),
- fTHnFAL(0),
+ // fTHnFK0s(0),
+ // fTHnFL(0),
+ // fTHnFAL(0),
fTHnFK0sDauEta(0),
fTHnFLDauEta(0),
fTHnFALDauEta(0),
// ------------------- add always ---------------------------//
//THnF
Double_t piForAx = 2.0*TMath::Pi();
- Int_t binsTHnV0K0s[4] = {150,100,80,18};
- Int_t binsTHnV0L[4] = {100,100,80,18};
+ // Int_t binsTHnV0K0s[4] = {150,100,80,18};
+ // Int_t binsTHnV0L[4] = {100,100,80,18};
Int_t binsTHnV0DauEtaK0s[4] = {150,100,40,18};
Int_t binsTHnV0DauEtaL[4] = {100,100,40,18};
Int_t binsTHnV0DauPhiK0s[4] = {150, 100,40,18};
Int_t binsTHnV0DauPhiL[4] = {100, 100,40,18};
- Double_t minK0s[4] = {0.35,0.0,0.0,0.0};
- Double_t maxK0s[4] = {0.65,50.0,40.0,piForAx};
+ //Double_t minK0s[4] = {0.35,0.0,0.0,0.0};
+ //Double_t maxK0s[4] = {0.65,50.0,40.0,piForAx};
Double_t minK0sDauEta[4] = {0.35, 0.0,-0.8,0.0};
Double_t maxK0sDauEta[4] = {0.65,50.0, 0.8,piForAx};
Double_t minK0sDauPhi[4] = {0.35,0.0,-0.8,0.0};
Double_t maxK0sDauPhi[4] = {0.65,50.0,0.8,piForAx};
- Double_t minL[4] = {1.07, 0.0,0.0,0.0};
- Double_t maxL[4] = {1.17,50.0, 40.0,piForAx};
+ //Double_t minL[4] = {1.07, 0.0,0.0,0.0};
+ // Double_t maxL[4] = {1.17,50.0, 40.0,piForAx};
Double_t minLDauEta[4] = {1.07, 0.0,-0.8,0.0};
Double_t maxLDauEta[4] = {1.17,50.0, 0.8, piForAx};
Double_t minLDauPhi[4] = {1.07,0.0,-0.8,0.0};
Double_t maxLDauPhi[4] = {1.17,50.0,0.8,piForAx};
-
+ /*
char histTitK0s[255];
snprintf(histTitK0s,255,"fTHnFK0s");
char histTitL[255];
snprintf(histTitL,255,"fTHnFL");
char histTitAL[255];
snprintf(histTitAL,255,"fTHnFAL");
-
+ */
char histTitK0sDauEta[255];
snprintf(histTitK0sDauEta,255,"fTHnFK0sDauEta");
snprintf(histTitLDauPhi,255,"fTHnFLDauPhi");
char histTitALDauPhi[255];
snprintf(histTitALDauPhi,255,"fTHnFALDauPhi");
-
+ /*
char axTitK0s[255];
snprintf(axTitK0s,255,"K^{0}_{s};m_{inv} (GeV/c^{2});p_{T} (GeV/c);c#tau(V0);#phi(V0)");
char axTitL[255];
snprintf(axTitL,255,"#Lambda;m_{inv} (GeV/c^{2});p_{T} (GeV/c);c#tau(V0);#phi(V0)");
char axTitAL[255];
snprintf(axTitAL,255,"#bar{#Lambda};m_{inv} (GeV/c^{2});p_{T} (GeV/c);c#tau(V0);#phi(V0)");
-
+ */
char axTitK0sDauEta[255];
snprintf(axTitK0sDauEta,255,"K^{0}_{s} daughter;m_{inv} (GeV/c^{2});p_{T} (Gev/c);#eta_{pos};#phi_{pos}");
char axTitALDauPhi[255];
snprintf(axTitALDauPhi,255,"#bar{#Lambda} daughter;m_{inv} (GeV/c^{2});p_{T}(GeV/c);#eta_{neg};#phi_{neg}");
-
+ /*
fTHnFK0s = new THnF(histTitK0s,axTitK0s,4,binsTHnV0K0s,minK0s,maxK0s);
// fTHnFK0s->Sumw2();
fTHnFL = new THnF(histTitL ,axTitL ,4,binsTHnV0L,minL ,maxL);
// fTHnFL->Sumw2();
fTHnFAL = new THnF(histTitAL ,axTitAL ,4,binsTHnV0L,minL ,maxL);
// fTHnFAL->Sumw2();
-
+ */
fTHnFK0sDauEta = new THnF(histTitK0sDauEta,axTitK0sDauEta,4,binsTHnV0DauEtaK0s,minK0sDauEta,maxK0sDauEta);
// fTHnFK0sDauEta->Sumw2();
fOutputContainer->Add(fHistPiPiMonitorMCCuts);
fOutputContainer->Add(fHistPiPiK0sVsLambdaMass);
fOutputContainer->Add(fHistPiPiK0sVsALambdaMass);
- fOutputContainer->Add(fTHnFK0s);
+ // fOutputContainer->Add(fTHnFK0s);
fOutputContainer->Add(fTHnFK0sDauEta);
fOutputContainer->Add(fTHnFK0sDauPhi);
//fOutputContainer->Add(fHistPiPiPhiPosVsPtPosVsMass);//xxx
fOutputContainer->Add(fHistPiPK0sVsLambdaMass);
fOutputContainer->Add(fHistPiPALambdaVsLambdaMass);
// fOutputContainer->Add(fHistPiPPhiPosVsPtPosVsMass);//xxx
- fOutputContainer->Add(fTHnFL);
+ //fOutputContainer->Add(fTHnFL);
fOutputContainer->Add(fTHnFLDauEta);
fOutputContainer->Add(fTHnFLDauPhi);
fOutputContainer->Add(fHistPiAPK0sVsALambdaMass);
fOutputContainer->Add(fHistPiAPLambdaVsALambdaMass);
// if(fMCTruthMode) fOutputContainer->Add(fHistPiAPPhiPosVsPtPosVsMass);//xxx
- fOutputContainer->Add(fTHnFAL);
+ //fOutputContainer->Add(fTHnFAL);
fOutputContainer->Add(fTHnFALDauEta);
fOutputContainer->Add(fTHnFALDauPhi);
fHistPiAPMonitorMCCuts[isSecd]->Fill(16*fillFlagAL);
- Double_t phiMC = p0->Phi();
+ // Double_t phiMC = p0->Phi();
//Double_t etaMC = p0->Eta();
- Double_t valTHnMC[4] = {massV0MC,ptV0MC,0.0,phiMC};
+ // Double_t valTHnMC[4] = {massV0MC,ptV0MC,0.0,phiMC};
Double_t valTHnMCDauEta[4] = {massV0MC,ptV0MC,etaMCPos,phiMCPos};
Double_t valTHnMCDauPhi[4] = {massV0MC,ptV0MC,etaMCNeg,phiMCNeg};
// fHistPiPiPhiPosVsPtPosVsMass->Fill(massV0MC,ctTK0s,ptV0MC);//,ctK0s);//phiPosMC);//xxx
fHistPiPiK0sVsLambdaMass->Fill(calcLambdamass,calcK0smass);
fHistPiPiK0sVsALambdaMass->Fill(calcALambdamass,calcK0smass);
- valTHnMC[2] = ctTK0s;
- fTHnFK0s->Fill(valTHnMC);
+ // valTHnMC[2] = ctTK0s;
+ // fTHnFK0s->Fill(valTHnMC);
fTHnFK0sDauEta->Fill(valTHnMCDauEta);
fTHnFK0sDauPhi->Fill(valTHnMCDauPhi);
}
fHistPiPK0sVsLambdaMass->Fill(calcLambdamass,calcK0smass);
- valTHnMC[2] = ctTL;
-fTHnFL->Fill(valTHnMC);
+ //valTHnMC[2] = ctTL;
+ //fTHnFL->Fill(valTHnMC);
fTHnFLDauEta->Fill(valTHnMCDauEta);
fTHnFLDauPhi->Fill(valTHnMCDauPhi);
}
// if(isSecd <1) fHistPiPPhiPosVsPtPosVsMass->Fill(massV0MC,ctTL,ptV0MC);//,ctK0s);//phiPosMC);//xxx
//else fHistPiAPPhiPosVsPtPosVsMass->Fill(massV0MC,ctTL,ptV0MC);//,ctK0s);//phiPosMC);//xxx
//fHistPiAPPhiPosVsPtPosVsMass->Fill(massV0MC,ctTL,ptV0MC);//,ctK0s);//phiPosMC);//xxx
- valTHnMC[2] = ctTAL;
- fTHnFAL->Fill(valTHnMC);
+ // valTHnMC[2] = ctTAL;
+ //fTHnFAL->Fill(valTHnMC);
fTHnFALDauEta->Fill(valTHnMCDauEta);
fTHnFALDauPhi->Fill(valTHnMCDauPhi);
}
// if( ppTrack.Angle(pmTrack)<0.001) continue;
// if( ppTrack.Angle(pmTrack)<0.004) continue;
- Double_t px = v0K0.Px();
- Double_t py = v0K0.Py();
- Double_t phi = TMath::Pi()+TMath::ATan2(-py, -px);
+ // Double_t px = v0K0.Px();
+ // Double_t py = v0K0.Py();
+ // Double_t phi = TMath::Pi()+TMath::ATan2(-py, -px);
// Double_t eta = v0K0.Eta();
/*
fHistPiPiPtVSY->Fill(rapK0s,ptK0s);
fHistPiPiDecayLengthVsMass->Fill(massK0s,dim2V0Radius);//decayLength);
// fHistPiPiPhiPosVsPtPosVsMass->Fill(massK0s,ctTK0,ptV0MC);//,ctK0);//posDaughterPhi);//xxx
- Double_t valTHnK0s[4]= {massK0s,ptV0MC,ctTK0,phi};
+ //Double_t valTHnK0s[4]= {massK0s,ptV0MC,ctTK0,phi};
Double_t valTHnK0sDauEta[4]= {massK0s,ptV0MC,posDaughterEta,posDaughterPhi};
Double_t valTHnK0sDauPhi[4]= {massK0s,ptV0MC,negDaughterEta,negDaughterPhi};
if(massK0s >=0.35 && massK0s <= 0.65){
- fTHnFK0s->Fill(valTHnK0s);
+ // fTHnFK0s->Fill(valTHnK0s);
fTHnFK0sDauEta->Fill(valTHnK0sDauEta);
fTHnFK0sDauPhi->Fill(valTHnK0sDauPhi);
}
//Double_t valTHnLDauEta[4]= {massLambda,ptV0MC,posDaughterEta,negDaughterEta};
// Double_t valTHnLDauPhi[5]= {massLambda,posDaughterPhi,negDaughterPhi,double(nclsITSPos),double(nclsITSNeg)};
- Double_t valTHnL[4]= {massLambda,ptV0MC,ctTL,phi};
+ // Double_t valTHnL[4]= {massLambda,ptV0MC,ctTL,phi};
Double_t valTHnLDauEta[4]= {massLambda,ptV0MC,posDaughterEta,posDaughterPhi};
Double_t valTHnLDauPhi[4]= {massLambda,ptV0MC,negDaughterEta,negDaughterPhi};
if(massLambda >=1.07 && massLambda <= 1.17){
- fTHnFL->Fill(valTHnL);
+ // fTHnFL->Fill(valTHnL);
fTHnFLDauEta->Fill(valTHnLDauEta);
fTHnFLDauPhi->Fill(valTHnLDauPhi);
}
Double_t valTHnALDauEta[4]={massALambda,ptV0MC,posDaughterEta,negDaughterEta};
Double_t valTHnALDauPhi[5]={massALambda,posDaughterPhi,negDaughterPhi,double(nclsITSPos),double(nclsITSNeg)};
*/
- Double_t valTHnAL[4]= {massALambda,ptV0MC,ctTAL,phi};
+ //Double_t valTHnAL[4]= {massALambda,ptV0MC,ctTAL,phi};
Double_t valTHnALDauEta[4]= {massALambda,ptV0MC,posDaughterEta,posDaughterPhi};
Double_t valTHnALDauPhi[4]= {massALambda,ptV0MC,negDaughterEta,negDaughterPhi};
if(massALambda >=1.07 && massALambda <= 1.17){
fTHnFALDauEta->Fill(valTHnALDauEta);
fTHnFALDauPhi->Fill(valTHnALDauPhi);
- fTHnFAL->Fill(valTHnAL);
+ //fTHnFAL->Fill(valTHnAL);
}
if(fMCMode && !fMCTruthMode) fHistPiAPPDGCode->Fill(pdgBG);
if( massALambda>1.108 && massALambda<1.123 ) fHistPiAPDecayLengthVsPt[isSecd]->Fill(ptV0MC,dim2V0Radius);//decayLength);