#include "AliRDHFCuts.h"
#include "AliKFVertex.h"
#include "AliESDVertex.h"
-#include "AliTOFPIDResponse.h"
-#include "AliAODpidUtil.h"
+//#include "AliAODpidUtil.h"
#include "AliAODPid.h"
#include "AliInputEventHandler.h"
+#include "AliPID.h"
ClassImp(AliAnalysisTaskSELambdac)
//________________________________________________________________________
-AliAnalysisTaskSELambdac::AliAnalysisTaskSELambdac():
-AliAnalysisTaskSE(),
-fOutput(0),
-fHistNEvents(0),
-fhChi2(0),
-fhMassPtGreater3(0),
-fhMassPtGreater3TC(0),
-fhMassPtGreater3Kp(0),
-fhMassPtGreater3KpTC(0),
-fhMassPtGreater3Lpi(0),
-fhMassPtGreater3LpiTC(0),
-fhMassPtGreater2(0),
-fhMassPtGreater2TC(0),
-fhMassPtGreater2Kp(0),
-fhMassPtGreater2KpTC(0),
-fhMassPtGreater2Lpi(0),
-fhMassPtGreater2LpiTC(0),
-fNtupleLambdac(0),
-fUpmasslimit(2.486),
-fLowmasslimit(2.086),
-fNPtBins(0),
-fRDCutsAnalysis(0),
-fRDCutsProduction(0),
-fListCuts(0),
-fFillNtuple(kFALSE),
-fReadMC(kFALSE),
-fMCPid(kFALSE),
-fRealPid(kFALSE),
-fResPid(kTRUE),
-fUseKF(kFALSE),
-fAnalysis(kFALSE),
-fVHF(0),
-fFillVarHists(kTRUE),
-fNentries(0),
-fOutputMC(0),
-fUtilPid(0)
+ AliAnalysisTaskSELambdac::AliAnalysisTaskSELambdac():
+ AliAnalysisTaskSE(),
+ fOutput(0),
+ fHistNEvents(0),
+ fhChi2(0),
+ fhMassPtGreater3(0),
+ fhMassPtGreater3TC(0),
+ fhMassPtGreater3Kp(0),
+ fhMassPtGreater3KpTC(0),
+ fhMassPtGreater3Lpi(0),
+ fhMassPtGreater3LpiTC(0),
+ fhMassPtGreater2(0),
+ fhMassPtGreater2TC(0),
+ fhMassPtGreater2Kp(0),
+ fhMassPtGreater2KpTC(0),
+ fhMassPtGreater2Lpi(0),
+ fhMassPtGreater2LpiTC(0),
+ fNtupleLambdac(0),
+ fUpmasslimit(2.486),
+ fLowmasslimit(2.086),
+ fNPtBins(0),
+ fRDCutsAnalysis(0),
+ fRDCutsProduction(0),
+ fListCuts(0),
+ fFillNtuple(kFALSE),
+ fReadMC(kFALSE),
+ fMCPid(kFALSE),
+ fRealPid(kFALSE),
+ fResPid(kTRUE),
+ fUseKF(kFALSE),
+ fAnalysis(kFALSE),
+ fVHF(0),
+ fFillVarHists(kFALSE),
+ fMultiplicityHists(kFALSE),
+ fPriorsHists(kFALSE),
+ fNentries(0),
+ fOutputMC(0),
+ fAPriori(0),
+ fMultiplicity(0),
+ //fUtilPid(0),
+ fPIDResponse(0)
{
- // Default constructor
+ // Default constructor
}
//________________________________________________________________________
AliAnalysisTaskSELambdac::AliAnalysisTaskSELambdac(const char *name,Bool_t fillNtuple,AliRDHFCutsLctopKpi *lccutsana,AliRDHFCutsLctopKpi *lccutsprod):
-AliAnalysisTaskSE(name),
-fOutput(0),
-fHistNEvents(0),
-fhChi2(0),
-fhMassPtGreater3(0),
-fhMassPtGreater3TC(0),
-fhMassPtGreater3Kp(0),
-fhMassPtGreater3KpTC(0),
-fhMassPtGreater3Lpi(0),
-fhMassPtGreater3LpiTC(0),
-fhMassPtGreater2(0),
-fhMassPtGreater2TC(0),
-fhMassPtGreater2Kp(0),
-fhMassPtGreater2KpTC(0),
-fhMassPtGreater2Lpi(0),
-fhMassPtGreater2LpiTC(0),
-fNtupleLambdac(0),
-fUpmasslimit(2.486),
-fLowmasslimit(2.086),
-fNPtBins(0),
-fRDCutsAnalysis(lccutsana),
-fRDCutsProduction(lccutsprod),
-fListCuts(0),
-fFillNtuple(fillNtuple),
-fReadMC(kFALSE),
-fMCPid(kFALSE),
-fRealPid(kTRUE),
-fResPid(kFALSE),
-fUseKF(kFALSE),
-fAnalysis(kFALSE),
-fVHF(0),
-fFillVarHists(kTRUE),
-fNentries(0),
-fOutputMC(0),
-fUtilPid(0)
+ AliAnalysisTaskSE(name),
+ fOutput(0),
+ fHistNEvents(0),
+ fhChi2(0),
+ fhMassPtGreater3(0),
+ fhMassPtGreater3TC(0),
+ fhMassPtGreater3Kp(0),
+ fhMassPtGreater3KpTC(0),
+ fhMassPtGreater3Lpi(0),
+ fhMassPtGreater3LpiTC(0),
+ fhMassPtGreater2(0),
+ fhMassPtGreater2TC(0),
+ fhMassPtGreater2Kp(0),
+ fhMassPtGreater2KpTC(0),
+ fhMassPtGreater2Lpi(0),
+ fhMassPtGreater2LpiTC(0),
+ fNtupleLambdac(0),
+ fUpmasslimit(2.486),
+ fLowmasslimit(2.086),
+ fNPtBins(0),
+ fRDCutsAnalysis(lccutsana),
+ fRDCutsProduction(lccutsprod),
+ fListCuts(0),
+ fFillNtuple(fillNtuple),
+ fReadMC(kFALSE),
+ fMCPid(kFALSE),
+ fRealPid(kTRUE),
+ fResPid(kFALSE),
+ fUseKF(kFALSE),
+ fAnalysis(kFALSE),
+ fVHF(0),
+ fFillVarHists(kFALSE),
+ fMultiplicityHists(kFALSE),
+ fPriorsHists(kFALSE),
+ fNentries(0),
+ fOutputMC(0),
+ fAPriori(0),
+ fMultiplicity(0),
+ //fUtilPid(0),
+ fPIDResponse(0)
{
- SetPtBinLimit(fRDCutsAnalysis->GetNPtBins()+1,fRDCutsAnalysis->GetPtBinLimits());
+ SetPtBinLimit(fRDCutsAnalysis->GetNPtBins()+1,
+ fRDCutsAnalysis->GetPtBinLimits());
// Default constructor
- // Output slot #1 writes into a TList container
+ // Output slot #1 writes into a TList container
DefineOutput(1,TList::Class()); //My private output
DefineOutput(2,TList::Class());
DefineOutput(3,TList::Class());
DefineOutput(4,TH1F::Class());
- if(fFillNtuple){
+ DefineOutput(5,TList::Class());
+ DefineOutput(6,TList::Class());
+ if (fFillNtuple) {
// Output slot #2 writes into a TNtuple container
- DefineOutput(5,TNtuple::Class()); //My private output
+ DefineOutput(7,TNtuple::Class()); //My private output
}
+
}
//________________________________________________________________________
delete fOutputMC;
fOutputMC = 0;
}
+ if (fAPriori) {
+ delete fAPriori;
+ fAPriori = 0;
+ }
+ if (fMultiplicity) {
+ delete fMultiplicity;
+ fMultiplicity = 0;
+ }
if (fVHF) {
delete fVHF;
fVHF = 0;
}
- if(fRDCutsAnalysis){
+ if(fRDCutsAnalysis){
delete fRDCutsAnalysis;
fRDCutsAnalysis = 0;
- }
- if(fRDCutsProduction){
+ }
+ if(fRDCutsProduction){
delete fRDCutsProduction;
fRDCutsProduction = 0;
- }
-
- if (fListCuts) {
- delete fListCuts;
- fListCuts = 0;
- }
- if (fNentries){
- delete fNentries;
- fNentries = 0;
- }
- if (fUtilPid){
- delete fUtilPid;
- fUtilPid = 0;
- }
+ }
+
+ if (fListCuts) {
+ delete fListCuts;
+ fListCuts = 0;
+ }
+ if (fNentries){
+ delete fNentries;
+ fNentries = 0;
+ }
+ /*
+ if (fUtilPid){
+ delete fUtilPid;
+ fUtilPid = 0;
+ }
+ */
+ if (fPIDResponse) {
+ delete fPIDResponse;
+ }
+
}
//_________________________________________________________________
void AliAnalysisTaskSELambdac::SetMassLimits(Float_t range){
{
// Initialization
- if(fDebug > 1) printf("AnalysisTaskSELambdac::Init() \n");
+ if (fDebug > 1) printf("AnalysisTaskSELambdac::Init() \n");
fListCuts=new TList();
fListCuts->SetOwner();
{
// Create the output container
//
- if(fDebug > 1) printf("AnalysisTaskSELambdac::UserCreateOutputObjects() \n");
+ if (fDebug > 1) printf("AnalysisTaskSELambdac::UserCreateOutputObjects() \n");
// Several histograms are more conveniently managed in a TList
fOutput = new TList();
hisname.Form("hMassPtKp%dTC",i);
fMassHistKpTC[index]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
fMassHistKpTC[index]->Sumw2();
-//signal
+ //signal
index=GetSignalHistoIndex(i);
hisname.Form("hSigPt%d",i);
fMassHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
hisname.Form("hBkgPtKp%dTC",i);
fMassHistKpTC[index]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
fMassHistKpTC[index]->Sumw2();
-}
+ }
for(Int_t i=0; i<3*fNPtBins; i++){
fOutput->Add(fMassHist[i]);
fOutputMC->SetOwner();
fOutputMC->SetName("QAMCHistos");
-// const char* nameoutput=GetOutputSlot(4)->GetContainer()->GetName();
+ // const char* nameoutput=GetOutputSlot(4)->GetContainer()->GetName();
fNentries=new TH1F("fNentries", "Integral(1,2) = number of AODs *** Integral(2,3) = number of candidates selected with cuts *** Integral(3,4) = number of Lc selected with cuts *** Integral(4,5) = events with good vertex *** Integral(5,6) = pt out of bounds", 11,-0.5,10.5);
hisname.Form("hbpIDTot");
TH1F *hBProtonIDTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
fOutputMC->Add(hBProtonIDTot);
- hisname.Form("hpnonIDTot");
- TH1F *hProtonnonIDTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
- fOutputMC->Add(hProtonnonIDTot);
- hisname.Form("hbpnonIDTot");
- TH1F *hBProtonnonIDTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
- fOutputMC->Add(hBProtonnonIDTot);
+ hisname.Form("hnopIDp");
+ TH1F *hnoProtonIDpTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
+ fOutputMC->Add(hnoProtonIDpTot);
+ hisname.Form("hbnopIDp");
+ TH1F *hBnoProtonIDpTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
+ fOutputMC->Add(hBnoProtonIDpTot);
hisname.Form("hpd0Prong");
TH1F *hProtond0Prong=new TH1F(hisname.Data(),hisname.Data(),100,-0.1,0.1);
TH2F *hBpSigmaVspTPC=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,-10.0,10.0);
fOutputMC->Add(hBpSigmaVspTPC);
- //kaon specific
+ //kaon specific
hisname.Form("hKTOFSignal");
TH1F *hKaonTOFSignal=new TH1F(hisname.Data(),hisname.Data(),100,12000.,50000.0);
fOutputMC->Add(hKaonTOFSignal);
hisname.Form("hbKTOFSignal");
TH1F *hBKaonTOFSignal=new TH1F(hisname.Data(),hisname.Data(),100,12000.,50000.0);
fOutputMC->Add(hBKaonTOFSignal);
-
hisname.Form("hKTPCSignal");
TH1F *hKaonTPCSignal=new TH1F(hisname.Data(),hisname.Data(),150,0.,150.0);
fOutputMC->Add(hKaonTPCSignal);
hisname.Form("hbKptProng");
TH1F *hBKaonPtProng=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.0);
fOutputMC->Add(hBKaonPtProng);
-
hisname.Form("hKRealTot");
TH1F *hKaonRealTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
fOutputMC->Add(hKaonRealTot);
hisname.Form("hbKIDTot");
TH1F *hBKaonIDTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
fOutputMC->Add(hBKaonIDTot);
- hisname.Form("hKnonIDTot");
- TH1F *hKaonnonIDTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
- fOutputMC->Add(hKaonnonIDTot);
- hisname.Form("hbKnonIDTot");
- TH1F *hBKaonnonIDTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
- fOutputMC->Add(hBKaonnonIDTot);
+ hisname.Form("hnokIDk");
+ TH1F *hnoKaonIDkTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
+ fOutputMC->Add(hnoKaonIDkTot);
+ hisname.Form("hbnokIDk");
+ TH1F *hBnoKaonIDkTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
+ fOutputMC->Add(hBnoKaonIDkTot);
hisname.Form("hKd0Prong");
hisname.Form("hKSignalVspTPC");
TH2F *hKSignalVspTPC=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,150,0.,150.0);
fOutputMC->Add(hKSignalVspTPC);
-
hisname.Form("hKSigmaVspTOF");
TH2F *hKSigmaVspTOF=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,-10.0,10.0);
fOutputMC->Add(hKSigmaVspTOF);
-
hisname.Form("hKSigmaVspTPC");
TH2F *hKSigmaVspTPC=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,-10.0,10.0);
fOutputMC->Add(hKSigmaVspTPC);
-
hisname.Form("hbKSigmaVspTOF");
TH2F *hBKSigmaVspTOF=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,-10.0,10.0);
fOutputMC->Add(hBKSigmaVspTOF);
TH2F *hBKSigmaVspTPC=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,-10.0,10.0);
fOutputMC->Add(hBKSigmaVspTPC);
+
// pion specific
hisname.Form("hpiTOFSignal");
TH1F *hPionTOFSignal=new TH1F(hisname.Data(),hisname.Data(),100,12000.,50000.0);
hisname.Form("hbpiTOFSignal");
TH1F *hBPionTOFSignal=new TH1F(hisname.Data(),hisname.Data(),100,12000.,50000.0);
fOutputMC->Add(hBPionTOFSignal);
-
hisname.Form("hpiTPCSignal");
TH1F *hPionTPCSignal=new TH1F(hisname.Data(),hisname.Data(),100,30.,100.0);
fOutputMC->Add(hPionTPCSignal);
hisname.Form("hbpiIDGood");
TH1F *hBPionIDGood=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
fOutputMC->Add(hBPionIDGood);
- hisname.Form("hpinonIDTot");
- TH1F *hPionnonIDTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
- fOutputMC->Add(hPionnonIDTot);
- hisname.Form("hbpinonIDTot");
- TH1F *hBPionnonIDTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
- fOutputMC->Add(hBPionnonIDTot);
-
+ hisname.Form("hnopiIDpi");
+ TH1F *hnoPionIDpiTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
+ fOutputMC->Add(hnoPionIDpiTot);
+ hisname.Form("hbnopiIDpi");
+ TH1F *hBnoPionIDpiTot=new TH1F(hisname.Data(),hisname.Data(),100,0.,10.0);
+ fOutputMC->Add(hBnoPionIDpiTot);
hisname.Form("hpid0Prong");
TH1F *hPiond0Prong=new TH1F(hisname.Data(),hisname.Data(),100,-0.1,-0.1);
hisname.Form("hbpid0Prong");
TH1F *hBPiond0Prong=new TH1F(hisname.Data(),hisname.Data(),100,-0.1,0.1);
fOutputMC->Add(hBPiond0Prong);
- hisname.Form("hbpiSignalVspTOF");
- TH2F *hbpiSignalVspTOF=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,12000.,50000.0);
- fOutputMC->Add(hbpiSignalVspTOF);
- hisname.Form("hbpiSignalVspTPC");
- TH2F *hbpiSignalVspTPC=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,150,0.,150.0);
- fOutputMC->Add(hbpiSignalVspTPC);
hisname.Form("hpiSignalVspTOF");
TH2F *hpiSignalVspTOF=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,12000.,50000.0);
hisname.Form("hpiSignalVspTPC");
TH2F *hpiSignalVspTPC=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,150,0.,150.0);
fOutputMC->Add(hpiSignalVspTPC);
-
- hisname.Form("hbpiSigmaVspTOF");
- TH2F *hBpiSigmaVspTOF=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,-10.0,10.0);
- fOutputMC->Add(hBpiSigmaVspTOF);
-
- hisname.Form("hbpiSigmaVspTPC");
- TH2F *hBpiSigmaVspTPC=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,-10.0,10.0);
- fOutputMC->Add(hBpiSigmaVspTPC);
+ hisname.Form("hbpiSignalVspTOF");
+ TH2F *hbpiSignalVspTOF=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,12000.,50000.0);
+ fOutputMC->Add(hbpiSignalVspTOF);
+ hisname.Form("hbpiSignalVspTPC");
+ TH2F *hbpiSignalVspTPC=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,150,0.,150.0);
+ fOutputMC->Add(hbpiSignalVspTPC);
hisname.Form("hpiSigmaVspTOF");
TH2F *hpiSigmaVspTOF=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,-10.0,10.0);
fOutputMC->Add(hpiSigmaVspTOF);
-
hisname.Form("hpiSigmaVspTPC");
TH2F *hpiSigmaVspTPC=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,-10.0,10.0);
fOutputMC->Add(hpiSigmaVspTPC);
+ hisname.Form("hbpiSigmaVspTOF");
+ TH2F *hBpiSigmaVspTOF=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,-10.0,10.0);
+ fOutputMC->Add(hBpiSigmaVspTOF);
+ hisname.Form("hbpiSigmaVspTPC");
+ TH2F *hBpiSigmaVspTPC=new TH2F(hisname.Data(),hisname.Data(),100,0.,5.0,100,-10.0,10.0);
+ fOutputMC->Add(hBpiSigmaVspTPC);
// other generic
hisname.Form("hLcpt");
+ fAPriori = new TList(); // AdC
+ fAPriori->SetOwner(); // AdC
+ fAPriori->SetName("APrioriMCHistos"); // AdC
+
+ hisname.Form("hElIn3Prong");
+ TH1F *hElIn3Prong=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hElIn3Prong);
+ hisname.Form("hMuIn3Prong");
+ TH1F *hMuIn3Prong=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hMuIn3Prong);
+ hisname.Form("hPiIn3Prong");
+ TH1F *hPiIn3Prong=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPiIn3Prong);
+ hisname.Form("hKaIn3Prong");
+ TH1F *hKaIn3Prong=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hKaIn3Prong);
+ hisname.Form("hPrIn3Prong");
+ TH1F *hPrIn3Prong=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPrIn3Prong);
+
+ hisname.Form("hElIn3Prong1");
+ TH1F *hElIn3Prong1=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hElIn3Prong1);
+ hisname.Form("hMuIn3Prong1");
+ TH1F *hMuIn3Prong1=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hMuIn3Prong1);
+ hisname.Form("hPiIn3Prong1");
+ TH1F *hPiIn3Prong1=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPiIn3Prong1);
+ hisname.Form("hKaIn3Prong1");
+ TH1F *hKaIn3Prong1=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hKaIn3Prong1);
+ hisname.Form("hPrIn3Prong1");
+ TH1F *hPrIn3Prong1=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPrIn3Prong1);
+
+ hisname.Form("hElIn3Prong2");
+ TH1F *hElIn3Prong2=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hElIn3Prong2);
+ hisname.Form("hMuIn3Prong2");
+ TH1F *hMuIn3Prong2=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hMuIn3Prong2);
+ hisname.Form("hPiIn3Prong2");
+ TH1F *hPiIn3Prong2=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPiIn3Prong2);
+ hisname.Form("hKaIn3Prong2");
+ TH1F *hKaIn3Prong2=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hKaIn3Prong2);
+ hisname.Form("hPrIn3Prong2");
+ TH1F *hPrIn3Prong2=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPrIn3Prong2);
+
+ hisname.Form("hElIn3Prong3");
+ TH1F *hElIn3Prong3=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hElIn3Prong3);
+ hisname.Form("hMuIn3Prong3");
+ TH1F *hMuIn3Prong3=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hMuIn3Prong3);
+ hisname.Form("hPiIn3Prong3");
+ TH1F *hPiIn3Prong3=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPiIn3Prong3);
+ hisname.Form("hKaIn3Prong3");
+ TH1F *hKaIn3Prong3=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hKaIn3Prong3);
+ hisname.Form("hPrIn3Prong3");
+ TH1F *hPrIn3Prong3=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPrIn3Prong3);
+
+ hisname.Form("hElIn3Prong4");
+ TH1F *hElIn3Prong4=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hElIn3Prong4);
+ hisname.Form("hMuIn3Prong4");
+ TH1F *hMuIn3Prong4=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hMuIn3Prong4);
+ hisname.Form("hPiIn3Prong4");
+ TH1F *hPiIn3Prong4=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPiIn3Prong4);
+ hisname.Form("hKaIn3Prong4");
+ TH1F *hKaIn3Prong4=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hKaIn3Prong4);
+ hisname.Form("hPrIn3Prong4");
+ TH1F *hPrIn3Prong4=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPrIn3Prong4);
+
+ hisname.Form("hElIn3Prong5");
+ TH1F *hElIn3Prong5=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hElIn3Prong5);
+ hisname.Form("hMuIn3Prong5");
+ TH1F *hMuIn3Prong5=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hMuIn3Prong5);
+ hisname.Form("hPiIn3Prong5");
+ TH1F *hPiIn3Prong5=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPiIn3Prong5);
+ hisname.Form("hKaIn3Prong5");
+ TH1F *hKaIn3Prong5=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hKaIn3Prong5);
+ hisname.Form("hPrIn3Prong5");
+ TH1F *hPrIn3Prong5=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPrIn3Prong5);
+
+ hisname.Form("hElIn3Prong6");
+ TH1F *hElIn3Prong6=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hElIn3Prong6);
+ hisname.Form("hMuIn3Prong6");
+ TH1F *hMuIn3Prong6=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hMuIn3Prong6);
+ hisname.Form("hPiIn3Prong6");
+ TH1F *hPiIn3Prong6=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPiIn3Prong6);
+ hisname.Form("hKaIn3Prong6");
+ TH1F *hKaIn3Prong6=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hKaIn3Prong6);
+ hisname.Form("hPrIn3Prong6");
+ TH1F *hPrIn3Prong6=new TH1F(hisname.Data(),hisname.Data(),100,0.,5.);
+ fAPriori->Add(hPrIn3Prong6);
+
+
+ fMultiplicity = new TList(); // AdC
+ fMultiplicity->SetOwner(); // AdC
+ fMultiplicity->SetName("MultiplicityMCHistos"); // AdC
+
+ hisname.Form("hLcinEvent");
+ TH1I*hLcinEvent=new TH1I(hisname.Data(),hisname.Data(),3,-1,2);
+ fMultiplicity->Add(hLcinEvent);
+
+ hisname.Form("hPrimariesvsAOD");
+ TH2I *hPrimariesvsAOD=new TH2I(hisname.Data(),hisname.Data(),300,0,300,300,0,300);
+ fMultiplicity->Add(hPrimariesvsAOD);
+
+ hisname.Form("hMultiplicityInLcEvent");
+ TH1I *hMultiplicityInLcEvent=new TH1I(hisname.Data(),hisname.Data(),300,0,300);
+ fMultiplicity->Add(hMultiplicityInLcEvent);
+ hisname.Form("h2MultiplicityInLcEvent");
+ TH1I *h2MultiplicityInLcEvent=new TH1I(hisname.Data(),hisname.Data(),300,0,300);
+ fMultiplicity->Add(h2MultiplicityInLcEvent);
+ hisname.Form("hAll2MultiplicityInEvent");
+ TH1I *hAll2MultiplicityInEvent=new TH1I(hisname.Data(),hisname.Data(),300,0,300);
+ fMultiplicity->Add(hAll2MultiplicityInEvent);
+ hisname.Form("hAllMultiplicityInEvent");
+ TH1I *hAllMultiplicityInEvent=new TH1I(hisname.Data(),hisname.Data(),300,0,300);
+ fMultiplicity->Add(hAllMultiplicityInEvent);
+ hisname.Form("hAllMultiplicityPrimaryInEvent");
+ TH1I *hAllMultiplicityPrimaryInEvent=new TH1I(hisname.Data(),hisname.Data(),300,0,300);
+ fMultiplicity->Add(hAllMultiplicityPrimaryInEvent);
+ hisname.Form("hAll2MultiplicityPrimaryInEvent");
+ TH1I *hAll2MultiplicityPrimaryInEvent=new TH1I(hisname.Data(),hisname.Data(),300,0,300);
+ fMultiplicity->Add(hAll2MultiplicityPrimaryInEvent);
+ hisname.Form("hMultiplicityInEvent");
+ TH1I *hMultiplicityInEvent=new TH1I(hisname.Data(),hisname.Data(),300,0.,300);
+ fMultiplicity->Add(hMultiplicityInEvent);
+ hisname.Form("hMultiplicityIn3ProngLC");
+ TH1I *hMultiplicityIn3ProngLC=new TH1I(hisname.Data(),hisname.Data(),300,0,300);
+ fMultiplicity->Add(hMultiplicityIn3ProngLC);
+ hisname.Form("hMultiplicityInLCpid");
+ TH1I *hMultiplicityInLCpid=new TH1I(hisname.Data(),hisname.Data(),300,0,300);
+ fMultiplicity->Add(hMultiplicityInLCpid);
+ hisname.Form("hMultiplicityInLCmc");
+ TH1I *hMultiplicityInLCmc=new TH1I(hisname.Data(),hisname.Data(),300,0,300);
+ fMultiplicity->Add(hMultiplicityInLCmc);
+ hisname.Form("hMultiplicityInLCNomc");
+ TH1I *hMultiplicityInLCNomc=new TH1I(hisname.Data(),hisname.Data(),300,0,300);
+ fMultiplicity->Add(hMultiplicityInLCNomc);
+ hisname.Form("hMultiplicityYesC");
+ TH1I *hMultiplicityYesC=new TH1I(hisname.Data(),hisname.Data(),300,0,300);
+ fMultiplicity->Add(hMultiplicityYesC);
+ hisname.Form("hMultiplicityYesB");
+ TH1I *hMultiplicityYesB=new TH1I(hisname.Data(),hisname.Data(),300,0,300);
+ fMultiplicity->Add(hMultiplicityYesB);
+ hisname.Form("hMultiplicityJPsi");
+ TH1I *hMultiplicityJPsi=new TH1I(hisname.Data(),hisname.Data(),300,0,300);
+ fMultiplicity->Add(hMultiplicityJPsi);
+
+
+ AliAnalysisManager *man=AliAnalysisManager::GetAnalysisManager();
+ AliInputEventHandler* inputHandler = (AliInputEventHandler*) (man->GetInputEventHandler());
+ fPIDResponse = inputHandler->GetPIDResponse();
+
if(fRDCutsProduction->GetIsUsePID()){
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- AliInputEventHandler *inputHandler=(AliInputEventHandler*)mgr->GetInputEventHandler();
- AliPIDResponse *pidResp=inputHandler->GetPIDResponse();
- fRDCutsProduction->GetPidHF()->SetPidResponse(pidResp);
- fRDCutsProduction->GetPidpion()->SetPidResponse(pidResp);
- fRDCutsProduction->GetPidprot()->SetPidResponse(pidResp);
- fUtilPid=new AliAODpidUtil(pidResp);
- fUtilPid=new AliAODpidUtil();
- }
+ fRDCutsProduction->GetPidHF()->SetPidResponse(fPIDResponse);
+ fRDCutsProduction->GetPidpion()->SetPidResponse(fPIDResponse);
+ fRDCutsProduction->GetPidprot()->SetPidResponse(fPIDResponse);
+ //fUtilPid=new AliAODpidUtil(fPIDResponse);
+ fRDCutsProduction->GetPidHF()->SetOldPid(kFALSE);
+ fRDCutsProduction->GetPidpion()->SetOldPid(kFALSE);
+ fRDCutsProduction->GetPidprot()->SetOldPid(kFALSE);
+ }
if(fRDCutsAnalysis->GetIsUsePID()){
- AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
- AliInputEventHandler *inputHandler=(AliInputEventHandler*)mgr->GetInputEventHandler();
- AliPIDResponse *pidResp=inputHandler->GetPIDResponse();
- fRDCutsAnalysis->GetPidHF()->SetPidResponse(pidResp);
- fRDCutsAnalysis->GetPidpion()->SetPidResponse(pidResp);
- fRDCutsAnalysis->GetPidprot()->SetPidResponse(pidResp);
+ fRDCutsAnalysis->GetPidHF()->SetPidResponse(fPIDResponse);
+ fRDCutsAnalysis->GetPidpion()->SetPidResponse(fPIDResponse);
+ fRDCutsAnalysis->GetPidprot()->SetPidResponse(fPIDResponse);
+ fRDCutsAnalysis->GetPidHF()->SetOldPid(kFALSE);
+ fRDCutsAnalysis->GetPidpion()->SetOldPid(kFALSE);
+ fRDCutsAnalysis->GetPidprot()->SetOldPid(kFALSE);
}
-
PostData(1,fOutput);
- if(fFillVarHists) PostData(3,fOutputMC);
+ if (fFillVarHists) PostData(3,fOutputMC);
PostData(4,fNentries);
- if(fFillNtuple){
+ if (fPriorsHists) PostData(5,fAPriori);
+ if (fMultiplicityHists) PostData(6,fMultiplicity);
+ if (fFillNtuple) {
//OpenFile(3); // 2 is the slot number of the ntuple
-
- fNtupleLambdac = new TNtuple("fNtupleLambdac","D +","pdg:Px:Py:Pz:PtTrue:VxTrue:VyTrue:VzTrue:Ptpi:PtK:Ptpi2:PtRec:PointingAngle:DecLeng:VxRec:VyRec:VzRec:InvMass:sigvert:d0Pi:d0K:d0Pi2:dca:d0square");
- PostData(5,fNtupleLambdac);
-
+ fNtupleLambdac = new TNtuple("fNtupleLambdac","D +",
+ "pdg:Px:Py:Pz:PtTrue:VxTrue:VyTrue:VzTrue:Ptpi:PtK:Ptpi2:PtRec:PointingAngle:DecLeng:VxRec:VyRec:VzRec:InvMass:sigvert:d0Pi:d0K:d0Pi2:dca:d0square");
+ PostData(7,fNtupleLambdac);
}
-
+
return;
}
// Execute analysis for current event:
// heavy flavor candidates association to MC truth
-
AliAODEvent *aod = dynamic_cast<AliAODEvent*> (InputEvent());
- //tmp
- fHistNEvents->Fill(0); // count event
+ AliAODEvent *aod = dynamic_cast<AliAODEvent*> (InputEvent());
+ //tmp
+ fHistNEvents->Fill(0); // count event
// Post the data already here
TClonesArray *array3Prong = 0;
arrayLikeSign=(TClonesArray*)aod->GetList()->FindObject("LikeSign3Prong");
}
- if(!array3Prong || !aod) {
- printf("AliAnalysisTaskSELambdac::UserExec: Charm3Prong branch not found!\n");
- return;
- }
- if(!arrayLikeSign) {
- printf("AliAnalysisTaskSELambdac::UserExec: LikeSign3Prong branch not found!\n");
- // return;
- }
+ if(!aod) return;
-
TClonesArray *arrayMC=0;
AliAODMCHeader *mcHeader=0;
- fNentries->Fill(0);
- TString trigclass=aod->GetFiredTriggerClasses();
- if(trigclass.Contains("C0SMH-B-NOPF-ALLNOTRD") || trigclass.Contains("C0SMH-B-NOPF-ALL")) fNentries->Fill(14);
- if(!fRDCutsProduction->IsEventSelected(aod)) {
- if(fRDCutsProduction->GetWhyRejection()==1) // rejected for pileup
- fNentries->Fill(13);
- return;
- }
-
- // AOD primary vertex
- AliAODVertex *vtx1 = (AliAODVertex*)aod->GetPrimaryVertex();
- if(!vtx1) return;
-
-
// load MC particles
if(fReadMC){
printf("AliAnalysisTaskSELambdac::UserExec: MC particles branch not found!\n");
return;
}
-
- // load MC header
+
+
+ // load MC header
mcHeader = (AliAODMCHeader*)aod->GetList()->FindObject(AliAODMCHeader::StdBranchName());
if(!mcHeader) {
- printf("AliAnalysisTaskSELambdac::UserExec: MC header branch not found!\n");
- return;
+ printf("AliAnalysisTaskSELambdac::UserExec: MC header branch not found!\n");
+ return;
}
}
-
+
+ TString fillthis="";
+ Int_t numberOfPrimaries= NumberPrimaries(aod);
+
+ if (fMultiplicityHists && fReadMC) {
+ fillthis="hPrimariesvsAOD";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(aod->GetNTracks(),numberOfPrimaries);
+
+ fillthis="hAll2MultiplicityInEvent";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(aod->GetNTracks());
+
+ fillthis="hAll2MultiplicityPrimaryInEvent";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(numberOfPrimaries);
+
+ if (IsThereAGeneratedLc(arrayMC)) {
+ fillthis="h2MultiplicityInLcEvent";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(numberOfPrimaries);
+ }
+
+ }
+
+ if(!array3Prong || !aod) {
+ printf("AliAnalysisTaskSELambdac::UserExec: Charm3Prong branch not found!\n");
+ return;
+ }
+ if(!arrayLikeSign) {
+ printf("AliAnalysisTaskSELambdac::UserExec: LikeSign3Prong branch not found!\n");
+ // return;
+ }
+
+ fNentries->Fill(0);
+ TString trigclass=aod->GetFiredTriggerClasses();
+ if(trigclass.Contains("C0SMH-B-NOPF-ALLNOTRD") || trigclass.Contains("C0SMH-B-NOPF-ALL")) fNentries->Fill(14);
+ if(!fRDCutsProduction->IsEventSelected(aod)) {
+ if(fRDCutsProduction->GetWhyRejection()==1) // rejected for pileup
+ fNentries->Fill(13);
+ return;
+ }
+
+ // AOD primary vertex
+ AliAODVertex *vtx1 = (AliAODVertex*)aod->GetPrimaryVertex();
+ if(!vtx1) return;
+
+ Bool_t isThereA3prongWithGoodTracks = kFALSE;
+ Bool_t isThereA3ProngLcKine = kFALSE;
+ Bool_t isThereA3ProngLcKineANDpid = kFALSE;
+ Bool_t isThereA3ProngLcMC = kFALSE;
+ Bool_t isThereA3ProngCyes = kFALSE;
+ Bool_t isThereA3ProngByes = kFALSE;
+ Bool_t isThereA3ProngJPsi = kFALSE;
+
Int_t n3Prong = array3Prong->GetEntriesFast();
-
-
for (Int_t i3Prong = 0; i3Prong < n3Prong; i3Prong++) {
AliAODRecoDecayHF3Prong *d = (AliAODRecoDecayHF3Prong*)array3Prong->UncheckedAt(i3Prong);
-
-
+
Bool_t unsetvtx=kFALSE;
if(!d->GetOwnPrimaryVtx()){
d->SetOwnPrimaryVtx(vtx1);
unsetvtx=kTRUE;
}
-// if(d->GetSelectionMap()) {if(!d->HasSelectionBit(AliRDHFCuts::kLcCuts)) continue;}
-// if(d->GetSelectionMap()) if(!d->HasSelectionBit(AliRDHFCuts::kLcPID)) continue;
+ //if(d->GetSelectionMap()) {if(!d->HasSelectionBit(AliRDHFCuts::kLcCuts)) continue;}
+ //if(d->GetSelectionMap()) {if(!d->HasSelectionBit(AliRDHFCuts::kLcPID)) continue;}
Int_t isSelectedTracks = fRDCutsProduction->IsSelected(d,AliRDHFCuts::kTracks,aod);
if(!isSelectedTracks) continue;
+ isThereA3prongWithGoodTracks=kTRUE;
+
if (fRDCutsProduction->IsInFiducialAcceptance(d->Pt(),d->Y(4122))) fNentries->Fill(6);
Int_t ptbin=fRDCutsProduction->PtBin(d->Pt());
if(ptbin==-1) {fNentries->Fill(4); continue;} //out of bounds
FillMassHists(aod,d,arrayMC,fRDCutsProduction);
- if(fFillVarHists) FillVarHists(d,arrayMC,fRDCutsProduction,fOutputMC,aod);
-
- /*
- //start OS analysis
- if(labDp<0)fHistOSbkg->Fill(d->InvMassDplus());
- fHistOS->Fill(d->InvMassDplus());
- */
+ if(fFillVarHists) FillVarHists(d,arrayMC,fRDCutsProduction,/*fOutputMC,*/aod);
+
+ if (fPriorsHists && fReadMC)
+ FillAPrioriConcentrations(d, fRDCutsProduction, aod, arrayMC);
+ if (fMultiplicityHists && fReadMC)
+ MultiplicityStudies(d, fRDCutsProduction, aod, arrayMC,
+ isThereA3ProngLcKine,isThereA3ProngLcKineANDpid,isThereA3ProngLcMC,
+ isThereA3ProngCyes,isThereA3ProngByes,isThereA3ProngJPsi);
+
+ /*
+ //start OS analysis
+ if(labDp<0)fHistOSbkg->Fill(d->InvMassDplus());
+ fHistOS->Fill(d->InvMassDplus());
+ */
if(unsetvtx) d->UnsetOwnPrimaryVtx();
}
+ if (fMultiplicityHists && fReadMC) {
+ fillthis="hAllMultiplicityInEvent";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(aod->GetNTracks());
+
+ fillthis="hAllMultiplicityPrimaryInEvent";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(numberOfPrimaries);
+
+ if (IsThereAGeneratedLc(arrayMC)) {
+ fillthis="hMultiplicityInLcEvent";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(numberOfPrimaries);
+ }
+
+ if (isThereA3prongWithGoodTracks) {
+ fillthis="hMultiplicityInEvent";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(numberOfPrimaries);
+ }
+ if (isThereA3ProngLcKine) {
+ fillthis="hMultiplicityIn3ProngLC";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(numberOfPrimaries);
+ }
+ if (isThereA3ProngLcKineANDpid) {
+ fillthis="hMultiplicityInLCpid";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(numberOfPrimaries);
+ }
+ if (isThereA3ProngLcMC) {
+ fillthis="hMultiplicityInLCmc";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(numberOfPrimaries);
+ }
+ if (isThereA3ProngLcKine && !isThereA3ProngLcMC) {
+ fillthis="hMultiplicityInLCNomc";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(numberOfPrimaries);
+ }
+
+ if (isThereA3ProngCyes) {
+ fillthis="hMultiplicityYesC";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(numberOfPrimaries);
+ }
+ if (isThereA3ProngByes) {
+ fillthis="hMultiplicityYesB";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(numberOfPrimaries);
+ }
+ if (isThereA3ProngJPsi) {
+ fillthis="hMultiplicityJPsi";
+ ((TH1I*)fMultiplicity->FindObject(fillthis))->Fill(numberOfPrimaries);
+ }
+
+ }
+
PostData(1,fOutput);
- if(fFillVarHists) PostData(3,fOutputMC);
+ if (fFillVarHists) PostData(3,fOutputMC);
+ if (fPriorsHists) PostData(5,fAPriori);
+ if (fMultiplicityHists) PostData(6,fMultiplicity);
+
PostData(4,fNentries);
return;
{
// Terminate analysis
//
- if(fDebug > 1) printf("AnalysisTaskSELambdac: Terminate() \n");
+ if (fDebug > 1) printf("AnalysisTaskSELambdac: Terminate() \n");
fOutput = dynamic_cast<TList*> (GetOutputData(1));
if (!fOutput) {
printf("ERROR: fOutput not available\n");
return;
}
- //fHistNEvents = dynamic_cast<TH1F*>(fOutput->FindObject("fHistNEvents"));
+ //fHistNEvents = dynamic_cast<TH1F*>(fOutput->FindObject("fHistNEvents"));
- TString hisname;
if(fFillNtuple){
fNtupleLambdac = dynamic_cast<TNtuple*>(GetOutputData(3));
}
- return;
+ return;
}
//________________________________________________________________________
Int_t lambdacLab[3]={0,0,0};
Int_t pdgs[3]={0,0,0};
for(Int_t i=0;i<3;i++){
- AliAODTrack *daugh=(AliAODTrack*)d->GetDaughter(i);
- Int_t lab=daugh->GetLabel();
- if(lab<0) return 0;
- AliAODMCParticle *part= (AliAODMCParticle*)arrayMC->At(lab);
- if(!part) continue;
- pdgs[i]=part->GetPdgCode();
- Int_t partPdgcode = TMath::Abs(part->GetPdgCode());
- if(partPdgcode==211 || partPdgcode==321 || partPdgcode==2212){
- Int_t motherLabel=part->GetMother();
- if(motherLabel<0) return 0;
- AliAODMCParticle *motherPart = (AliAODMCParticle*)arrayMC->At(motherLabel);
- if(!motherPart) continue;
- Int_t motherPdg = TMath::Abs(motherPart->GetPdgCode());
- if(motherPdg==4122) {
+ AliAODTrack *daugh=(AliAODTrack*)d->GetDaughter(i);
+ Int_t lab=daugh->GetLabel();
+ if(lab<0) return 0;
+ AliAODMCParticle *part= (AliAODMCParticle*)arrayMC->At(lab);
+ if(!part) continue;
+ pdgs[i]=part->GetPdgCode();
+ Int_t partPdgcode = TMath::Abs(part->GetPdgCode());
+ if(partPdgcode==211 || partPdgcode==321 || partPdgcode==2212){
+ Int_t motherLabel=part->GetMother();
+ if(motherLabel<0) return 0;
+ AliAODMCParticle *motherPart = (AliAODMCParticle*)arrayMC->At(motherLabel);
+ if(!motherPart) continue;
+ Int_t motherPdg = TMath::Abs(motherPart->GetPdgCode());
+ if(motherPdg==4122) {
if(GetLambdacDaugh(motherPart,arrayMC)){lambdacLab[i]=motherLabel;continue;}
- }
- if(motherPdg==313 || motherPdg==2224 || motherPdg==3124){
- Int_t GmotherLabel=motherPart->GetMother();
- if(GmotherLabel<0) return 0;
- AliAODMCParticle *GmotherPart = (AliAODMCParticle*)arrayMC->At(GmotherLabel);
- if(!GmotherPart) continue;
- Int_t GmotherPdg = TMath::Abs(GmotherPart->GetPdgCode());
- if(GmotherPdg==4122) {
- if(GetLambdacDaugh(GmotherPart,arrayMC)) {lambdacLab[i]=GmotherLabel;continue;}
- }
+ }
+ if(motherPdg==313 || motherPdg==2224 || motherPdg==3124){
+ Int_t GmotherLabel=motherPart->GetMother();
+ if(GmotherLabel<0) return 0;
+ AliAODMCParticle *GmotherPart = (AliAODMCParticle*)arrayMC->At(GmotherLabel);
+ if(!GmotherPart) continue;
+ Int_t GmotherPdg = TMath::Abs(GmotherPart->GetPdgCode());
+ if(GmotherPdg==4122) {
+ if(GetLambdacDaugh(GmotherPart,arrayMC)) {lambdacLab[i]=GmotherLabel;continue;}
}
- }
+ }
+ }
}
- if(lambdacLab[0]==lambdacLab[1] && lambdacLab[1]==lambdacLab[2]) {return lambdacLab[0];}
- return 0;
+ if(lambdacLab[0]==lambdacLab[1] && lambdacLab[1]==lambdacLab[2]) {return lambdacLab[0];}
+ return 0;
}
//------------------------
Bool_t AliAnalysisTaskSELambdac::GetLambdacDaugh(AliAODMCParticle *part,TClonesArray *arrayMC) const{
- // check if the particle is a lambdac and if its decay mode is the correct one
- Int_t numberOfLambdac=0;
- if(TMath::Abs(part->GetPdgCode())!=4122) return kFALSE;
- Int_t daugh_tmp[2];
- daugh_tmp[0]=part->GetDaughter(0);
- daugh_tmp[1]=part->GetDaughter(1);
- Int_t nDaugh = (Int_t)part->GetNDaughters();
- if(nDaugh<2) return kFALSE;
- if(nDaugh>3) return kFALSE;
- AliAODMCParticle* pdaugh1 = (AliAODMCParticle*)arrayMC->At(part->GetDaughter(0));
- if(!pdaugh1) {return kFALSE;}
- Int_t number1 = TMath::Abs(pdaugh1->GetPdgCode());
- AliAODMCParticle* pdaugh2 = (AliAODMCParticle*)arrayMC->At(part->GetDaughter(1));
- if(!pdaugh2) {return kFALSE;}
+ // check if the particle is a lambdac and if its decay mode is the correct one
+ Int_t numberOfLambdac=0;
+ if(TMath::Abs(part->GetPdgCode())!=4122) return kFALSE;
+ Int_t daugh_tmp[2];
+ daugh_tmp[0]=part->GetDaughter(0);
+ daugh_tmp[1]=part->GetDaughter(1);
+ Int_t nDaugh = (Int_t)part->GetNDaughters();
+ if(nDaugh<2) return kFALSE;
+ if(nDaugh>3) return kFALSE;
+ AliAODMCParticle* pdaugh1 = (AliAODMCParticle*)arrayMC->At(part->GetDaughter(0));
+ if(!pdaugh1) {return kFALSE;}
+ Int_t number1 = TMath::Abs(pdaugh1->GetPdgCode());
+ AliAODMCParticle* pdaugh2 = (AliAODMCParticle*)arrayMC->At(part->GetDaughter(1));
+ if(!pdaugh2) {return kFALSE;}
Int_t number2 = TMath::Abs(pdaugh2->GetPdgCode());
if(nDaugh==3){
- Int_t thirdDaugh=part->GetDaughter(1)-1;
- AliAODMCParticle* pdaugh3 = (AliAODMCParticle*)arrayMC->At(thirdDaugh);
- Int_t number3 = TMath::Abs(pdaugh3->GetPdgCode());
- if((number1==321 && number2==211 && number3==2212) || (number1==211 && number2==321 && number3==2212) || (number1==211 && number2==2212 && number3==321) || (number1==321 && number2==2212 && number3==211) || (number1==2212 && number2==321 && number3==211) || (number1==2212 && number2==211 && number3==321)) {
- numberOfLambdac++;
- }
+ Int_t thirdDaugh=part->GetDaughter(1)-1;
+ AliAODMCParticle* pdaugh3 = (AliAODMCParticle*)arrayMC->At(thirdDaugh);
+ Int_t number3 = TMath::Abs(pdaugh3->GetPdgCode());
+ if((number1==321 && number2==211 && number3==2212) ||
+ (number1==211 && number2==321 && number3==2212) ||
+ (number1==211 && number2==2212 && number3==321) ||
+ (number1==321 && number2==2212 && number3==211) ||
+ (number1==2212 && number2==321 && number3==211) ||
+ (number1==2212 && number2==211 && number3==321)) {
+ numberOfLambdac++;
+ }
}
- if(nDaugh==2){
+ if(nDaugh==2){
- //Lambda resonant
+ //Lambda resonant
- //Lambda -> p K*0
- //
- Int_t nfiglieK=0;
-
- if((number1==2212 && number2==313)){
- nfiglieK=pdaugh2->GetNDaughters();
- if(nfiglieK!=2) return kFALSE;
- AliAODMCParticle* pdaughK1 = (AliAODMCParticle*)arrayMC->At(pdaugh2->GetDaughter(0));
- AliAODMCParticle* pdaughK2 = (AliAODMCParticle*)arrayMC->At(pdaugh2->GetDaughter(1));
- if(!pdaughK1) return kFALSE;
- if(!pdaughK2) return kFALSE;
- if((TMath::Abs(pdaughK1->GetPdgCode())==211 && TMath::Abs(pdaughK2->GetPdgCode())==321) || (TMath::Abs(pdaughK1->GetPdgCode())==321 && TMath::Abs(pdaughK2->GetPdgCode())==211)) numberOfLambdac++;
+ //Lambda -> p K*0
+ //
+ Int_t nfiglieK=0;
+
+ if((number1==2212 && number2==313)){
+ nfiglieK=pdaugh2->GetNDaughters();
+ if(nfiglieK!=2) return kFALSE;
+ AliAODMCParticle* pdaughK1 = (AliAODMCParticle*)arrayMC->At(pdaugh2->GetDaughter(0));
+ AliAODMCParticle* pdaughK2 = (AliAODMCParticle*)arrayMC->At(pdaugh2->GetDaughter(1));
+ if(!pdaughK1) return kFALSE;
+ if(!pdaughK2) return kFALSE;
+ if((TMath::Abs(pdaughK1->GetPdgCode())==211 && TMath::Abs(pdaughK2->GetPdgCode())==321) || (TMath::Abs(pdaughK1->GetPdgCode())==321 && TMath::Abs(pdaughK2->GetPdgCode())==211)) numberOfLambdac++;
+ }
+
+ if((number1==313 && number2==2212)){
+ nfiglieK=pdaugh1->GetNDaughters();
+ if(nfiglieK!=2) return kFALSE;
+ AliAODMCParticle* pdaughK1 = (AliAODMCParticle*)arrayMC->At(pdaugh1->GetDaughter(0));
+ AliAODMCParticle* pdaughK2 = (AliAODMCParticle*)arrayMC->At(pdaugh1->GetDaughter(1));
+ if(!pdaughK1) return kFALSE;
+ if(!pdaughK2) return kFALSE;
+ if((TMath::Abs(pdaughK1->GetPdgCode())==211 && TMath::Abs(pdaughK2->GetPdgCode())==321) || (TMath::Abs(pdaughK1->GetPdgCode())==321 && TMath::Abs(pdaughK2->GetPdgCode())==211)) numberOfLambdac++;
}
- if((number1==313 && number2==2212)){
- nfiglieK=pdaugh1->GetNDaughters();
- if(nfiglieK!=2) return kFALSE;
- AliAODMCParticle* pdaughK1 = (AliAODMCParticle*)arrayMC->At(pdaugh1->GetDaughter(0));
- AliAODMCParticle* pdaughK2 = (AliAODMCParticle*)arrayMC->At(pdaugh1->GetDaughter(1));
- if(!pdaughK1) return kFALSE;
- if(!pdaughK2) return kFALSE;
- if((TMath::Abs(pdaughK1->GetPdgCode())==211 && TMath::Abs(pdaughK2->GetPdgCode())==321) || (TMath::Abs(pdaughK1->GetPdgCode())==321 && TMath::Abs(pdaughK2->GetPdgCode())==211)) numberOfLambdac++;
- }
-
- //Lambda -> Delta++ k
- Int_t nfiglieDelta=0;
- if(number1==321 && number2==2224){
- nfiglieDelta=pdaugh2->GetNDaughters();
- if(nfiglieDelta!=2) return kFALSE;
- AliAODMCParticle *pdaughD1=(AliAODMCParticle*)arrayMC->At(pdaugh2->GetDaughter(0));
- AliAODMCParticle *pdaughD2=(AliAODMCParticle*)arrayMC->At(pdaugh2->GetDaughter(1));
- if(!pdaughD1) return kFALSE;
- if(!pdaughD2) return kFALSE;
- if((TMath::Abs(pdaughD1->GetPdgCode())==211 && TMath::Abs(pdaughD2->GetPdgCode())==2212) || (TMath::Abs(pdaughD1->GetPdgCode())==2212 && TMath::Abs(pdaughD2->GetPdgCode())==211)) numberOfLambdac++;
- }
- if(number1==2224 && number2==321){
- nfiglieDelta=pdaugh1->GetNDaughters();
- if(nfiglieDelta!=2) return kFALSE;
- AliAODMCParticle* pdaughD1 = (AliAODMCParticle*)arrayMC->At(pdaugh1->GetDaughter(0));
- AliAODMCParticle* pdaughD2 = (AliAODMCParticle*)arrayMC->At(pdaugh1->GetDaughter(1));
- if(!pdaughD1) return kFALSE;
- if(!pdaughD2) return kFALSE;
- if((TMath::Abs(pdaughD1->GetPdgCode())==211 && TMath::Abs(pdaughD2->GetPdgCode())==2212) || (TMath::Abs(pdaughD1->GetPdgCode())==2212 && TMath::Abs(pdaughD2->GetPdgCode())==211)) numberOfLambdac++;
- }
+ //Lambda -> Delta++ k
+ Int_t nfiglieDelta=0;
+ if(number1==321 && number2==2224){
+ nfiglieDelta=pdaugh2->GetNDaughters();
+ if(nfiglieDelta!=2) return kFALSE;
+ AliAODMCParticle *pdaughD1=(AliAODMCParticle*)arrayMC->At(pdaugh2->GetDaughter(0));
+ AliAODMCParticle *pdaughD2=(AliAODMCParticle*)arrayMC->At(pdaugh2->GetDaughter(1));
+ if(!pdaughD1) return kFALSE;
+ if(!pdaughD2) return kFALSE;
+ if((TMath::Abs(pdaughD1->GetPdgCode())==211 && TMath::Abs(pdaughD2->GetPdgCode())==2212) || (TMath::Abs(pdaughD1->GetPdgCode())==2212 && TMath::Abs(pdaughD2->GetPdgCode())==211)) numberOfLambdac++;
+ }
+ if(number1==2224 && number2==321){
+ nfiglieDelta=pdaugh1->GetNDaughters();
+ if(nfiglieDelta!=2) return kFALSE;
+ AliAODMCParticle* pdaughD1 = (AliAODMCParticle*)arrayMC->At(pdaugh1->GetDaughter(0));
+ AliAODMCParticle* pdaughD2 = (AliAODMCParticle*)arrayMC->At(pdaugh1->GetDaughter(1));
+ if(!pdaughD1) return kFALSE;
+ if(!pdaughD2) return kFALSE;
+ if((TMath::Abs(pdaughD1->GetPdgCode())==211 && TMath::Abs(pdaughD2->GetPdgCode())==2212) || (TMath::Abs(pdaughD1->GetPdgCode())==2212 && TMath::Abs(pdaughD2->GetPdgCode())==211)) numberOfLambdac++;
+ }
- //Lambdac -> Lambda(1520) pi
- Int_t nfiglieLa=0;
- if(number1==3124 && number2==211){
- nfiglieLa=pdaugh1->GetNDaughters();
- if(nfiglieLa!=2) return kFALSE;
- AliAODMCParticle *pdaughL1=(AliAODMCParticle*)arrayMC->At(pdaugh1->GetDaughter(0));
- AliAODMCParticle *pdaughL2=(AliAODMCParticle*)arrayMC->At(pdaugh1->GetDaughter(1));
- if(!pdaughL1) return kFALSE;
- if(!pdaughL2) return kFALSE;
- if((TMath::Abs(pdaughL1->GetPdgCode())==321 && TMath::Abs(pdaughL2->GetPdgCode())==2212) || (TMath::Abs(pdaughL1->GetPdgCode())==2212 && TMath::Abs(pdaughL2->GetPdgCode())==321)) numberOfLambdac++;
- }
- if(number1==211 && number2==3124){
- nfiglieLa=pdaugh2->GetNDaughters();
- if(nfiglieLa!=2) return kFALSE;
- AliAODMCParticle *pdaughL1=(AliAODMCParticle*)arrayMC->At(pdaugh2->GetDaughter(0));
- AliAODMCParticle *pdaughL2=(AliAODMCParticle*)arrayMC->At(pdaugh2->GetDaughter(1));
- if(!pdaughL1) return kFALSE;
- if(!pdaughL2) return kFALSE;
- if((TMath::Abs(pdaughL1->GetPdgCode())==321 && TMath::Abs(pdaughL2->GetPdgCode())==2212) || (TMath::Abs(pdaughL1->GetPdgCode())==2212 && TMath::Abs(pdaughL2->GetPdgCode())==321)) numberOfLambdac++;
+ //Lambdac -> Lambda(1520) pi
+ Int_t nfiglieLa=0;
+ if(number1==3124 && number2==211){
+ nfiglieLa=pdaugh1->GetNDaughters();
+ if(nfiglieLa!=2) return kFALSE;
+ AliAODMCParticle *pdaughL1=(AliAODMCParticle*)arrayMC->At(pdaugh1->GetDaughter(0));
+ AliAODMCParticle *pdaughL2=(AliAODMCParticle*)arrayMC->At(pdaugh1->GetDaughter(1));
+ if(!pdaughL1) return kFALSE;
+ if(!pdaughL2) return kFALSE;
+ if((TMath::Abs(pdaughL1->GetPdgCode())==321 && TMath::Abs(pdaughL2->GetPdgCode())==2212) || (TMath::Abs(pdaughL1->GetPdgCode())==2212 && TMath::Abs(pdaughL2->GetPdgCode())==321)) numberOfLambdac++;
+ }
+ if(number1==211 && number2==3124){
+ nfiglieLa=pdaugh2->GetNDaughters();
+ if(nfiglieLa!=2) return kFALSE;
+ AliAODMCParticle *pdaughL1=(AliAODMCParticle*)arrayMC->At(pdaugh2->GetDaughter(0));
+ AliAODMCParticle *pdaughL2=(AliAODMCParticle*)arrayMC->At(pdaugh2->GetDaughter(1));
+ if(!pdaughL1) return kFALSE;
+ if(!pdaughL2) return kFALSE;
+ if((TMath::Abs(pdaughL1->GetPdgCode())==321 && TMath::Abs(pdaughL2->GetPdgCode())==2212) || (TMath::Abs(pdaughL1->GetPdgCode())==2212 && TMath::Abs(pdaughL2->GetPdgCode())==321)) numberOfLambdac++;
- }
- }
+ }
+ }
- if(numberOfLambdac>0) {return kTRUE;}
- return kFALSE;
+ if(numberOfLambdac>0) {return kTRUE;}
+ return kFALSE;
}
//-----------------------------
Bool_t AliAnalysisTaskSELambdac::IspKpiMC(AliAODRecoDecayHF3Prong *d,TClonesArray *arrayMC) const{
// Apply MC PID
- Int_t lab[3]={0,0,0},pdgs[3]={0,0,0};
- for(Int_t i=0;i<3;i++){
+ Int_t lab[3]={0,0,0},pdgs[3]={0,0,0};
+ for(Int_t i=0;i<3;i++){
AliAODTrack *daugh=(AliAODTrack*)d->GetDaughter(i);
lab[i]=daugh->GetLabel();
if(lab[i]<0) return kFALSE;
AliAODMCParticle *part= (AliAODMCParticle*)arrayMC->At(lab[i]);
if(!part) return kFALSE;
pdgs[i]=TMath::Abs(part->GetPdgCode());
- }
+ }
- if(pdgs[0]==2212 && pdgs[1]==321 && pdgs[2]==211) return kTRUE;
+ if(pdgs[0]==2212 && pdgs[1]==321 && pdgs[2]==211) return kTRUE;
- return kFALSE;
+ return kFALSE;
}
//-----------------------------
Bool_t AliAnalysisTaskSELambdac::IspiKpMC(AliAODRecoDecayHF3Prong *d,TClonesArray *arrayMC) const{
// Apply MC PID
- Int_t lab[3]={0,0,0},pdgs[3]={0,0,0};
- for(Int_t i=0;i<3;i++){
+ Int_t lab[3]={0,0,0},pdgs[3]={0,0,0};
+ for(Int_t i=0;i<3;i++){
AliAODTrack *daugh=(AliAODTrack*)d->GetDaughter(i);
lab[i]=daugh->GetLabel();
if(lab[i]<0) return kFALSE;
AliAODMCParticle *part= (AliAODMCParticle*)arrayMC->At(lab[i]);
if(!part) return kFALSE;
pdgs[i]=TMath::Abs(part->GetPdgCode());
- }
+ }
- if(pdgs[2]==2212 && pdgs[1]==321 && pdgs[0]==211) {return kTRUE;}
+ if(pdgs[2]==2212 && pdgs[1]==321 && pdgs[0]==211) {return kTRUE;}
- return kFALSE;
+ return kFALSE;
}
//--------------------------------------
Bool_t AliAnalysisTaskSELambdac::VertexingKF(AliAODRecoDecayHF3Prong *d,Int_t *pdgs,Double_t field) const{
- // apply vertexing KF
- Int_t iprongs[3]={0,1,2};
- Double_t mass[2]={0.,0.};
- //topological constr
- AliKFParticle *lambdac=d->ApplyVertexingKF(iprongs,3,pdgs,kTRUE,field,mass);
- if(!lambdac) return kFALSE;
-// Double_t probTot=TMath::Prob(Lambdac->GetChi2(),Lambdac->GetNDF());
-// if(probTot<fCutsKF[0]) return kFALSE;
+ // apply vertexing KF
+ Int_t iprongs[3]={0,1,2};
+ Double_t mass[2]={0.,0.};
+ //topological constr
+ AliKFParticle *lambdac=d->ApplyVertexingKF(iprongs,3,pdgs,kTRUE,field,mass);
+ if(!lambdac) return kFALSE;
+ // Double_t probTot=TMath::Prob(Lambdac->GetChi2(),Lambdac->GetNDF());
+ // if(probTot<fCutsKF[0]) return kFALSE;
if(lambdac->GetChi2()>fCutsKF[0]) return kFALSE;
- //mass constr for K*
- Int_t ipRes[2];
- Int_t pdgres[2];
- mass[0]=0.8961;mass[1]=0.03;
- if(TMath::Abs(pdgs[0])==211){
+ //mass constr for K*
+ Int_t ipRes[2];
+ Int_t pdgres[2];
+ mass[0]=0.8961;mass[1]=0.03;
+ if(TMath::Abs(pdgs[0])==211){
ipRes[0]=0;ipRes[1]=1;
pdgres[0]=pdgs[0];pdgres[1]=321;
- }
- if(TMath::Abs(pdgs[2])==211){
+ }
+ if(TMath::Abs(pdgs[2])==211){
ipRes[0]=2;ipRes[1]=1;
pdgres[0]=pdgs[2];pdgres[1]=321;
- }
- AliKFParticle *kappaStar=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
+ }
+ AliKFParticle *kappaStar=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
Double_t probKstar=TMath::Prob(kappaStar->GetChi2(),kappaStar->GetNDF());
if(probKstar>fCutsKF[1]) {
AliKFParticle prong2(*esdProng2,pdgres[1]);
if(kappaStar->GetPt()<fCutsKF[2] && prong1.GetAngle(prong2)>fCutsKF[3]) return kFALSE;
}
- //mass constr for Lambda
- mass[0]=1.520;mass[1]=0.005;
- if(TMath::Abs(pdgs[0])==2212){
+ //mass constr for Lambda
+ mass[0]=1.520;mass[1]=0.005;
+ if(TMath::Abs(pdgs[0])==2212){
ipRes[0]=0;ipRes[1]=1;
pdgres[0]=pdgs[0];pdgres[1]=pdgs[1];
- }
- if(TMath::Abs(pdgs[2])==2212){
+ }
+ if(TMath::Abs(pdgs[2])==2212){
ipRes[0]=2;ipRes[1]=1;
pdgres[0]=pdgs[2];pdgres[1]=pdgs[1];
- }
- AliKFParticle *lambda1520=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
+ }
+ AliKFParticle *lambda1520=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
Double_t probLa=TMath::Prob(lambda1520->GetChi2(),lambda1520->GetNDF());
if(probLa>fCutsKF[4]) {
AliAODTrack *esdProng1=(AliAODTrack*)d->GetDaughter(ipRes[0]);
AliKFParticle prong2(*esdProng2,pdgres[1]);
if(lambda1520->GetPt()<fCutsKF[5] && prong1.GetAngle(prong2)>fCutsKF[6]) return kFALSE;
}
- //mass constr for Delta
- mass[0]=1.2;mass[1]=0.15;
- ipRes[0]=0;ipRes[1]=2;
- pdgres[0]=pdgs[0];pdgres[1]=pdgs[2];
- AliKFParticle *delta=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
+ //mass constr for Delta
+ mass[0]=1.2;mass[1]=0.15;
+ ipRes[0]=0;ipRes[1]=2;
+ pdgres[0]=pdgs[0];pdgres[1]=pdgs[2];
+ AliKFParticle *delta=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
Double_t probDelta=TMath::Prob(delta->GetChi2(),delta->GetNDF());
if(probDelta>fCutsKF[7]) {
AliAODTrack *esdProng1=(AliAODTrack*)d->GetDaughter(ipRes[0]);
AliKFParticle prong2(*esdProng2,pdgres[1]);
if(delta->GetPt()<fCutsKF[8] && prong1.GetAngle(prong2)>fCutsKF[9]) return kFALSE;
}
- return kTRUE;
+ return kTRUE;
}
//-------------------------------------
Int_t AliAnalysisTaskSELambdac::IspiKpResonant(AliAODRecoDecayHF3Prong *d,Double_t field) const{
- // apply PID using the resonant channels
- //if lambda* -> pk
- Double_t mass[2]={1.520,0.005};
- Int_t ipRes[2]={1,2};
- Int_t pdgres[2]={321,2212};
- AliKFParticle *lambda1520=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
- Double_t probLa=TMath::Prob(lambda1520->GetChi2(),lambda1520->GetNDF());
- if(probLa>0.1) return 1;
- //K* -> kpi
- mass[0]=0.8961;mass[1]=0.03;
- ipRes[0]=0;ipRes[1]=1;
- pdgres[0]=211;pdgres[1]=321;
- AliKFParticle *kstar=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
- Double_t probKa=TMath::Prob(kstar->GetChi2(),kstar->GetNDF());
- if(probKa>0.1) return 2;
-
- return 0;
+ // apply PID using the resonant channels
+ //if lambda* -> pk
+ Double_t mass[2]={1.520,0.005};
+ Int_t ipRes[2]={1,2};
+ Int_t pdgres[2]={321,2212};
+ AliKFParticle *lambda1520=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
+ Double_t probLa=TMath::Prob(lambda1520->GetChi2(),lambda1520->GetNDF());
+ if(probLa>0.1) return 1;
+ //K* -> kpi
+ mass[0]=0.8961;mass[1]=0.03;
+ ipRes[0]=0;ipRes[1]=1;
+ pdgres[0]=211;pdgres[1]=321;
+ AliKFParticle *kstar=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
+ Double_t probKa=TMath::Prob(kstar->GetChi2(),kstar->GetNDF());
+ if(probKa>0.1) return 2;
+
+ return 0;
}
//-------------------------------------
Int_t AliAnalysisTaskSELambdac::IspKpiResonant(AliAODRecoDecayHF3Prong *d,Double_t field) const{
- // apply PID using the resonant channels
- //if lambda* -> pk
- Double_t mass[2]={1.520,0.005};
- Int_t ipRes[2]={0,1};
- Int_t pdgres[2]={2212,321};
- AliKFParticle *lambda1520=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
- Double_t probLa=TMath::Prob(lambda1520->GetChi2(),lambda1520->GetNDF());
- if(probLa>0.1) return 1;
- //K* -> kpi
- mass[0]=0.8961;mass[1]=0.03;
- ipRes[0]=1;ipRes[1]=2;
- pdgres[1]=211;pdgres[0]=321;
- AliKFParticle *kstar=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
- Double_t probKa=TMath::Prob(kstar->GetChi2(),kstar->GetNDF());
- if(probKa>0.1) return 2;
-
- return 0;
+ // apply PID using the resonant channels
+ //if lambda* -> pk
+ Double_t mass[2]={1.520,0.005};
+ Int_t ipRes[2]={0,1};
+ Int_t pdgres[2]={2212,321};
+ AliKFParticle *lambda1520=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
+ Double_t probLa=TMath::Prob(lambda1520->GetChi2(),lambda1520->GetNDF());
+ if(probLa>0.1) return 1;
+ //K* -> kpi
+ mass[0]=0.8961;mass[1]=0.03;
+ ipRes[0]=1;ipRes[1]=2;
+ pdgres[1]=211;pdgres[0]=321;
+ AliKFParticle *kstar=d->ApplyVertexingKF(ipRes,2,pdgres,kFALSE,field,mass);
+ Double_t probKa=TMath::Prob(kstar->GetChi2(),kstar->GetNDF());
+ if(probKa>0.1) return 2;
+
+ return 0;
}
//---------------------------
-void AliAnalysisTaskSELambdac::FillMassHists(AliAODEvent *aod,AliAODRecoDecayHF3Prong *part, TClonesArray *arrayMC, AliRDHFCutsLctopKpi *cuts){
-
- //if MC PID or no PID, unset PID
- if(!fRealPid) cuts->SetUsePID(kFALSE);
- Int_t selection=cuts->IsSelected(part,AliRDHFCuts::kCandidate,aod);
- if(selection>0){
- Int_t iPtBin = -1;
- Double_t ptCand = part->Pt();
- Int_t index=0;
-
- for(Int_t ibin=0;ibin<fNPtBins&&iPtBin<0&&ptCand>fArrayBinLimits[0]&&ptCand<fArrayBinLimits[fNPtBins];ibin++){
- if(ptCand<fArrayBinLimits[ibin+1])iPtBin=ibin;
- }
-
-
- Float_t pdgCode=-2;
- Float_t pdgCode1=-1;
- Float_t pdgCode2=-1;
- Int_t labDp=-1;
- Float_t deltaPx=0.;
- Float_t deltaPy=0.;
- Float_t deltaPz=0.;
- Float_t truePt=0.;
- Float_t xDecay=0.;
- Float_t yDecay=0.;
- Float_t zDecay=0.;
-
- if(fReadMC){
- labDp = MatchToMCLambdac(part,arrayMC);
- if(labDp>0){
- AliAODMCParticle *partDp = (AliAODMCParticle*)arrayMC->At(labDp);
- AliAODMCParticle *dg0 = (AliAODMCParticle*)arrayMC->At(partDp->GetDaughter(0));
- AliAODMCParticle *dg1 = (AliAODMCParticle*)arrayMC->At(partDp->GetDaughter(1));
- deltaPx=partDp->Px()-part->Px();
- deltaPy=partDp->Py()-part->Py();
- deltaPz=partDp->Pz()-part->Pz();
- truePt=partDp->Pt();
- xDecay=dg0->Xv();
- yDecay=dg0->Yv();
- zDecay=dg0->Zv();
- pdgCode=TMath::Abs(partDp->GetPdgCode());
- pdgCode1=TMath::Abs(dg0->GetPdgCode());
- pdgCode2=TMath::Abs(dg1->GetPdgCode());
+void AliAnalysisTaskSELambdac::FillMassHists(AliAODEvent *aod,AliAODRecoDecayHF3Prong *part,
+ TClonesArray *arrayMC, AliRDHFCutsLctopKpi *cuts)
+{
- }else{
- pdgCode=-1;
- }
- }
-
- Double_t invMasspKpi=-1.;
- Double_t invMasspiKp=-1.;
- Double_t invMasspKpiLpi=-1.;
- Double_t invMasspiKpLpi=-1.;
- Double_t invMasspKpiKp=-1.;
- Double_t invMasspiKpKp=-1.;
- Int_t pdgs[3]={0,0,0};
- Double_t field=aod->GetMagneticField();
-//apply MC PID
- if(fReadMC && fMCPid){
-
- if(IspKpiMC(part,arrayMC)) {
- invMasspKpi=part->InvMassLcpKpi();
- if(fUseKF){
- pdgs[0]=2212;pdgs[1]=321;pdgs[2]=211;
- if(!VertexingKF(part,pdgs,field)) invMasspKpi=-1.;
- }
- }
- if(IspiKpMC(part,arrayMC)) {
- invMasspiKp=part->InvMassLcpiKp();
- if(fUseKF){
- pdgs[0]=211;pdgs[1]=321;pdgs[2]=2212;
- if(!VertexingKF(part,pdgs,field)) invMasspiKp=-1.;
- }
- }
- }
-
- // apply realistic PID
- if(fRealPid){
- if(selection==1 || selection==3) {
- invMasspKpi=part->InvMassLcpKpi();
- pdgs[0]=2212;pdgs[1]=321;pdgs[2]=211;
- if(fUseKF){
- pdgs[0]=2212;pdgs[1]=321;pdgs[2]=211;
- if(!VertexingKF(part,pdgs,field)) invMasspKpi=-1.;
- }
- }
- if(selection>=2) {
- invMasspiKp=part->InvMassLcpiKp();
- pdgs[0]=211;pdgs[1]=321;pdgs[2]=2212;
- if(fUseKF){
- pdgs[0]=211;pdgs[1]=321;pdgs[2]=2212;
- if(!VertexingKF(part,pdgs,field)) invMasspiKp=-1.;
- }
- }
- }
- //apply PID using resonances
- if(fResPid && fRealPid){
- if((selection==3 || selection==1)) {
- if(IspKpiResonant(part,field)==1){
- invMasspKpiLpi=part->InvMassLcpKpi();
+ //if MC PID or no PID, unset PID
+ if(!fRealPid) cuts->SetUsePID(kFALSE);
+ Int_t selection=cuts->IsSelected(part,AliRDHFCuts::kCandidate,aod);
+ if(selection>0){
+ Int_t iPtBin = -1;
+ Double_t ptCand = part->Pt();
+ Int_t index=0;
+
+ for(Int_t ibin=0;ibin<fNPtBins&&iPtBin<0&&ptCand>fArrayBinLimits[0]&&ptCand<fArrayBinLimits[fNPtBins];ibin++){
+ if(ptCand<fArrayBinLimits[ibin+1])iPtBin=ibin;
}
- if(IspKpiResonant(part,field)==2){
- invMasspKpiKp=part->InvMassLcpKpi();
+
+ Float_t pdgCode=-2;
+ Float_t pdgCode1=-1;
+ Float_t pdgCode2=-1;
+ Int_t labDp=-1;
+ Float_t deltaPx=0.;
+ Float_t deltaPy=0.;
+ Float_t deltaPz=0.;
+ Float_t truePt=0.;
+ Float_t xDecay=0.;
+ Float_t yDecay=0.;
+ Float_t zDecay=0.;
+
+ if(fReadMC){
+ labDp = MatchToMCLambdac(part,arrayMC);
+ if(labDp>0){
+ AliAODMCParticle *partDp = (AliAODMCParticle*)arrayMC->At(labDp);
+ AliAODMCParticle *dg0 = (AliAODMCParticle*)arrayMC->At(partDp->GetDaughter(0));
+ AliAODMCParticle *dg1 = (AliAODMCParticle*)arrayMC->At(partDp->GetDaughter(1));
+ deltaPx=partDp->Px()-part->Px();
+ deltaPy=partDp->Py()-part->Py();
+ deltaPz=partDp->Pz()-part->Pz();
+ truePt=partDp->Pt();
+ xDecay=dg0->Xv();
+ yDecay=dg0->Yv();
+ zDecay=dg0->Zv();
+ pdgCode=TMath::Abs(partDp->GetPdgCode());
+ pdgCode1=TMath::Abs(dg0->GetPdgCode());
+ pdgCode2=TMath::Abs(dg1->GetPdgCode());
+
+ }else{
+ pdgCode=-1;
+ }
}
- }
- if(selection>=2) {
- if(IspiKpResonant(part,field)==1){
- invMasspiKpLpi=part->InvMassLcpiKp();
+
+ Double_t invMasspKpi=-1.;
+ Double_t invMasspiKp=-1.;
+ Double_t invMasspKpiLpi=-1.;
+ Double_t invMasspiKpLpi=-1.;
+ Double_t invMasspKpiKp=-1.;
+ Double_t invMasspiKpKp=-1.;
+ Int_t pdgs[3]={0,0,0};
+ Double_t field=aod->GetMagneticField();
+ //apply MC PID
+ if(fReadMC && fMCPid){
+
+ if(IspKpiMC(part,arrayMC)) {
+ invMasspKpi=part->InvMassLcpKpi();
+ if(fUseKF){
+ pdgs[0]=2212;pdgs[1]=321;pdgs[2]=211;
+ if(!VertexingKF(part,pdgs,field)) invMasspKpi=-1.;
+ }
+ }
+ if(IspiKpMC(part,arrayMC)) {
+ invMasspiKp=part->InvMassLcpiKp();
+ if(fUseKF){
+ pdgs[0]=211;pdgs[1]=321;pdgs[2]=2212;
+ if(!VertexingKF(part,pdgs,field)) invMasspiKp=-1.;
+ }
+ }
}
- if(IspiKpResonant(part,field)==2){
- invMasspiKpKp=part->InvMassLcpiKp();
+
+ // apply realistic PID
+ if(fRealPid){
+ if(selection==1 || selection==3) {
+ invMasspKpi=part->InvMassLcpKpi();
+ pdgs[0]=2212;pdgs[1]=321;pdgs[2]=211;
+ if(fUseKF){
+ pdgs[0]=2212;pdgs[1]=321;pdgs[2]=211;
+ if(!VertexingKF(part,pdgs,field)) invMasspKpi=-1.;
+ }
+ }
+ if(selection>=2) {
+ invMasspiKp=part->InvMassLcpiKp();
+ pdgs[0]=211;pdgs[1]=321;pdgs[2]=2212;
+ if(fUseKF){
+ pdgs[0]=211;pdgs[1]=321;pdgs[2]=2212;
+ if(!VertexingKF(part,pdgs,field)) invMasspiKp=-1.;
+ }
+ }
+ }
+ //apply PID using resonances
+ if (fResPid && fRealPid) {
+ if (selection==3 || selection==1) {
+ if (IspKpiResonant(part,field)==1) {
+ invMasspKpiLpi=part->InvMassLcpKpi();
+ }
+ if(IspKpiResonant(part,field)==2){
+ invMasspKpiKp=part->InvMassLcpKpi();
+ }
+ }
+ if(selection>=2) {
+ if(IspiKpResonant(part,field)==1){
+ invMasspiKpLpi=part->InvMassLcpiKp();
+ }
+ if(IspiKpResonant(part,field)==2){
+ invMasspiKpKp=part->InvMassLcpiKp();
+ }
+ }
}
- }
- }
- // no PID
- if(!fResPid && !fRealPid && !fMCPid){
- if(selection==2 || selection==3) {
- invMasspiKp=part->InvMassLcpiKp();
- if(fUseKF){
- pdgs[0]=211;pdgs[1]=321;pdgs[2]=2212;
- if(!VertexingKF(part,pdgs,field)) invMasspiKp=-1.;
- }
- }
- if(selection==1 || selection==3){
- invMasspKpi=part->InvMassLcpKpi();
- if(fUseKF){
- pdgs[0]=2212;pdgs[1]=321;pdgs[2]=211;
- if(!VertexingKF(part,pdgs,field)) invMasspKpi=-1.;
- }
+ // no PID
+ if(!fResPid && !fRealPid && !fMCPid){
+ if(selection==2 || selection==3) {
+ invMasspiKp=part->InvMassLcpiKp();
+ if(fUseKF){
+ pdgs[0]=211;pdgs[1]=321;pdgs[2]=2212;
+ if(!VertexingKF(part,pdgs,field)) invMasspiKp=-1.;
+ }
+ }
+ if(selection==1 || selection==3){
+ invMasspKpi=part->InvMassLcpKpi();
+ if(fUseKF){
+ pdgs[0]=2212;pdgs[1]=321;pdgs[2]=211;
+ if(!VertexingKF(part,pdgs,field)) invMasspKpi=-1.;
+ }
+ }
}
- }
-
- if(invMasspiKp<0. && invMasspKpi<0.) return;
-
- Int_t passTightCuts=0;
- if(fAnalysis) passTightCuts=fRDCutsAnalysis->IsSelected(part,AliRDHFCuts::kCandidate,aod);
+ if(invMasspiKp<0. && invMasspKpi<0.) return;
-
- Float_t tmp[24];
- if(fFillNtuple){
- tmp[0]=pdgCode;
- tmp[1]=deltaPx;
- tmp[2]=deltaPy;
- tmp[3]=deltaPz;
- tmp[4]=truePt;
- tmp[5]=xDecay;
- tmp[6]=yDecay;
- tmp[7]=zDecay;
- if(pdgCode1==2212) {tmp[8]=part->PtProng(0);}else{tmp[8]=0.;}
- if(pdgCode1==211) {tmp[10]=part->PtProng(0);}else{tmp[10]=0.;}
- tmp[9]=part->PtProng(1);
- if(pdgCode2==211) {tmp[10]=part->PtProng(2);}else{tmp[10]=0.;}
- tmp[11]=part->Pt();
- tmp[12]=part->CosPointingAngle();
- tmp[13]=part->DecayLength();
- tmp[14]=part->Xv();
- tmp[15]=part->Yv();
- tmp[16]=part->Zv();
- if(invMasspiKp>0.) tmp[17]=invMasspiKp;
- if(invMasspKpi>0.) tmp[17]=invMasspKpi;
- tmp[18]=part->GetSigmaVert();
- tmp[19]=part->Getd0Prong(0);
- tmp[20]=part->Getd0Prong(1);
- tmp[21]=part->Getd0Prong(2);
- tmp[22]=part->GetDCA();
- tmp[23]=part->Prodd0d0();
- fNtupleLambdac->Fill(tmp);
- PostData(5,fNtupleLambdac);
- }
+ Int_t passTightCuts=0;
+ if(fAnalysis) passTightCuts=fRDCutsAnalysis->IsSelected(part,AliRDHFCuts::kCandidate,aod);
+
+ Float_t tmp[24];
+ if (fFillNtuple) {
+ tmp[0]=pdgCode;
+ tmp[1]=deltaPx;
+ tmp[2]=deltaPy;
+ tmp[3]=deltaPz;
+ tmp[4]=truePt;
+ tmp[5]=xDecay;
+ tmp[6]=yDecay;
+ tmp[7]=zDecay;
+ if(pdgCode1==2212) {tmp[8]=part->PtProng(0);}else{tmp[8]=0.;}
+ if(pdgCode1==211) {tmp[10]=part->PtProng(0);}else{tmp[10]=0.;}
+ tmp[9]=part->PtProng(1);
+ if(pdgCode2==211) {tmp[10]=part->PtProng(2);}else{tmp[10]=0.;}
+ tmp[11]=part->Pt();
+ tmp[12]=part->CosPointingAngle();
+ tmp[13]=part->DecayLength();
+ tmp[14]=part->Xv();
+ tmp[15]=part->Yv();
+ tmp[16]=part->Zv();
+ if(invMasspiKp>0.) tmp[17]=invMasspiKp;
+ if(invMasspKpi>0.) tmp[17]=invMasspKpi;
+ tmp[18]=part->GetSigmaVert();
+ tmp[19]=part->Getd0Prong(0);
+ tmp[20]=part->Getd0Prong(1);
+ tmp[21]=part->Getd0Prong(2);
+ tmp[22]=part->GetDCA();
+ tmp[23]=part->Prodd0d0();
+ fNtupleLambdac->Fill(tmp);
+ PostData(7,fNtupleLambdac);
+ }
if(part->Pt()>3.){
- if(invMasspiKp>0. && invMasspKpi>0.){
- if(invMasspiKp>0.) fhMassPtGreater3->Fill(invMasspiKp,0.5);
- if(invMasspKpi>0.) fhMassPtGreater3->Fill(invMasspKpi,0.5);
- }else{
- if(invMasspiKp>0.) fhMassPtGreater3->Fill(invMasspiKp);
- if(invMasspKpi>0.) fhMassPtGreater3->Fill(invMasspKpi);
- }
- if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
- if(invMasspiKpLpi>0.) fhMassPtGreater3Lpi->Fill(invMasspiKpLpi,0.5);
- if(invMasspKpiLpi>0.) fhMassPtGreater3Lpi->Fill(invMasspKpiLpi,0.5);
- }else{
- if(invMasspiKpLpi>0.) fhMassPtGreater3Lpi->Fill(invMasspiKpLpi);
- if(invMasspKpiLpi>0.) fhMassPtGreater3Lpi->Fill(invMasspKpiLpi);
- }
- if(invMasspiKpKp>0. && invMasspKpiKp>0.){
- if(invMasspiKpKp>0.) fhMassPtGreater3Kp->Fill(invMasspiKpKp,0.5);
- if(invMasspKpiKp>0.) fhMassPtGreater3Kp->Fill(invMasspKpiKp,0.5);
- }else{
- if(invMasspiKpKp>0.) fhMassPtGreater3Kp->Fill(invMasspiKpKp);
- if(invMasspKpiKp>0.) fhMassPtGreater3Kp->Fill(invMasspKpiKp);
- }
- if(passTightCuts>0){
if(invMasspiKp>0. && invMasspKpi>0.){
- if(invMasspiKp>0.) fhMassPtGreater3TC->Fill(invMasspiKp,0.5);
- if(invMasspKpi>0.) fhMassPtGreater3TC->Fill(invMasspKpi,0.5);
+ if(invMasspiKp>0.) fhMassPtGreater3->Fill(invMasspiKp,0.5);
+ if(invMasspKpi>0.) fhMassPtGreater3->Fill(invMasspKpi,0.5);
+ }else{
+ if(invMasspiKp>0.) fhMassPtGreater3->Fill(invMasspiKp);
+ if(invMasspKpi>0.) fhMassPtGreater3->Fill(invMasspKpi);
+ }
+ if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
+ if(invMasspiKpLpi>0.) fhMassPtGreater3Lpi->Fill(invMasspiKpLpi,0.5);
+ if(invMasspKpiLpi>0.) fhMassPtGreater3Lpi->Fill(invMasspKpiLpi,0.5);
+ }else{
+ if(invMasspiKpLpi>0.) fhMassPtGreater3Lpi->Fill(invMasspiKpLpi);
+ if(invMasspKpiLpi>0.) fhMassPtGreater3Lpi->Fill(invMasspKpiLpi);
+ }
+ if(invMasspiKpKp>0. && invMasspKpiKp>0.){
+ if(invMasspiKpKp>0.) fhMassPtGreater3Kp->Fill(invMasspiKpKp,0.5);
+ if(invMasspKpiKp>0.) fhMassPtGreater3Kp->Fill(invMasspKpiKp,0.5);
}else{
- if(invMasspiKp>0.) fhMassPtGreater3TC->Fill(invMasspiKp);
- if(invMasspKpi>0.) fhMassPtGreater3TC->Fill(invMasspKpi);
+ if(invMasspiKpKp>0.) fhMassPtGreater3Kp->Fill(invMasspiKpKp);
+ if(invMasspKpiKp>0.) fhMassPtGreater3Kp->Fill(invMasspKpiKp);
+ }
+ if(passTightCuts>0){
+ if(invMasspiKp>0. && invMasspKpi>0.){
+ if(invMasspiKp>0.) fhMassPtGreater3TC->Fill(invMasspiKp,0.5);
+ if(invMasspKpi>0.) fhMassPtGreater3TC->Fill(invMasspKpi,0.5);
+ }else{
+ if(invMasspiKp>0.) fhMassPtGreater3TC->Fill(invMasspiKp);
+ if(invMasspKpi>0.) fhMassPtGreater3TC->Fill(invMasspKpi);
+ }
+ if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
+ if(invMasspiKpLpi>0.) fhMassPtGreater3LpiTC->Fill(invMasspiKpLpi,0.5);
+ if(invMasspKpiLpi>0.) fhMassPtGreater3LpiTC->Fill(invMasspKpiLpi,0.5);
+ }else{
+ if(invMasspiKpLpi>0.) fhMassPtGreater3LpiTC->Fill(invMasspiKpLpi);
+ if(invMasspKpiLpi>0.) fhMassPtGreater3LpiTC->Fill(invMasspKpiLpi);
+ }
+ if(invMasspiKpKp>0. && invMasspKpiKp>0.){
+ if(invMasspiKpKp>0.) fhMassPtGreater3KpTC->Fill(invMasspiKpKp,0.5);
+ if(invMasspKpiKp>0.) fhMassPtGreater3KpTC->Fill(invMasspKpiKp,0.5);
+ }else{
+ if(invMasspiKpKp>0.) fhMassPtGreater3KpTC->Fill(invMasspiKpKp);
+ if(invMasspKpiKp>0.) fhMassPtGreater3KpTC->Fill(invMasspKpiKp);
+ }
}
- if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
- if(invMasspiKpLpi>0.) fhMassPtGreater3LpiTC->Fill(invMasspiKpLpi,0.5);
- if(invMasspKpiLpi>0.) fhMassPtGreater3LpiTC->Fill(invMasspKpiLpi,0.5);
- }else{
- if(invMasspiKpLpi>0.) fhMassPtGreater3LpiTC->Fill(invMasspiKpLpi);
- if(invMasspKpiLpi>0.) fhMassPtGreater3LpiTC->Fill(invMasspKpiLpi);
- }
- if(invMasspiKpKp>0. && invMasspKpiKp>0.){
- if(invMasspiKpKp>0.) fhMassPtGreater3KpTC->Fill(invMasspiKpKp,0.5);
- if(invMasspKpiKp>0.) fhMassPtGreater3KpTC->Fill(invMasspKpiKp,0.5);
- }else{
- if(invMasspiKpKp>0.) fhMassPtGreater3KpTC->Fill(invMasspiKpKp);
- if(invMasspKpiKp>0.) fhMassPtGreater3KpTC->Fill(invMasspKpiKp);
- }
- }
}
- if(part->Pt()>2.){
- if(invMasspiKp>0. && invMasspKpi>0.){
- if(invMasspiKp>0.) fhMassPtGreater2->Fill(invMasspiKp,0.5);
- if(invMasspKpi>0.) fhMassPtGreater2->Fill(invMasspKpi,0.5);
- }else{
- if(invMasspiKp>0.) fhMassPtGreater2->Fill(invMasspiKp);
- if(invMasspKpi>0.) fhMassPtGreater2->Fill(invMasspKpi);
- }
- if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
- if(invMasspiKpLpi>0.) fhMassPtGreater2Lpi->Fill(invMasspiKpLpi,0.5);
- if(invMasspKpiLpi>0.) fhMassPtGreater2Lpi->Fill(invMasspKpiLpi,0.5);
- }else{
- if(invMasspiKpLpi>0.) fhMassPtGreater2Lpi->Fill(invMasspiKpLpi);
- if(invMasspKpiLpi>0.) fhMassPtGreater2Lpi->Fill(invMasspKpiLpi);
- }
- if(invMasspiKpKp>0. && invMasspKpiKp>0.){
- if(invMasspiKpKp>0.) fhMassPtGreater2Kp->Fill(invMasspiKpKp,0.5);
- if(invMasspKpiKp>0.) fhMassPtGreater2Kp->Fill(invMasspKpiKp,0.5);
- }else{
- if(invMasspiKpKp>0.) fhMassPtGreater2Kp->Fill(invMasspiKpKp);
- if(invMasspKpiKp>0.) fhMassPtGreater2Kp->Fill(invMasspKpiKp);
- }
- if(passTightCuts>0){
+ if(part->Pt()>2.){
if(invMasspiKp>0. && invMasspKpi>0.){
- if(invMasspiKp>0.) fhMassPtGreater2TC->Fill(invMasspiKp,0.5);
- if(invMasspKpi>0.) fhMassPtGreater2TC->Fill(invMasspKpi,0.5);
+ if(invMasspiKp>0.) fhMassPtGreater2->Fill(invMasspiKp,0.5);
+ if(invMasspKpi>0.) fhMassPtGreater2->Fill(invMasspKpi,0.5);
+ }else{
+ if(invMasspiKp>0.) fhMassPtGreater2->Fill(invMasspiKp);
+ if(invMasspKpi>0.) fhMassPtGreater2->Fill(invMasspKpi);
+ }
+ if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
+ if(invMasspiKpLpi>0.) fhMassPtGreater2Lpi->Fill(invMasspiKpLpi,0.5);
+ if(invMasspKpiLpi>0.) fhMassPtGreater2Lpi->Fill(invMasspKpiLpi,0.5);
+ }else{
+ if(invMasspiKpLpi>0.) fhMassPtGreater2Lpi->Fill(invMasspiKpLpi);
+ if(invMasspKpiLpi>0.) fhMassPtGreater2Lpi->Fill(invMasspKpiLpi);
+ }
+ if(invMasspiKpKp>0. && invMasspKpiKp>0.){
+ if(invMasspiKpKp>0.) fhMassPtGreater2Kp->Fill(invMasspiKpKp,0.5);
+ if(invMasspKpiKp>0.) fhMassPtGreater2Kp->Fill(invMasspKpiKp,0.5);
}else{
- if(invMasspiKp>0.) fhMassPtGreater2TC->Fill(invMasspiKp);
- if(invMasspKpi>0.) fhMassPtGreater2TC->Fill(invMasspKpi);
+ if(invMasspiKpKp>0.) fhMassPtGreater2Kp->Fill(invMasspiKpKp);
+ if(invMasspKpiKp>0.) fhMassPtGreater2Kp->Fill(invMasspKpiKp);
+ }
+ if(passTightCuts>0){
+ if(invMasspiKp>0. && invMasspKpi>0.){
+ if(invMasspiKp>0.) fhMassPtGreater2TC->Fill(invMasspiKp,0.5);
+ if(invMasspKpi>0.) fhMassPtGreater2TC->Fill(invMasspKpi,0.5);
+ }else{
+ if(invMasspiKp>0.) fhMassPtGreater2TC->Fill(invMasspiKp);
+ if(invMasspKpi>0.) fhMassPtGreater2TC->Fill(invMasspKpi);
+ }
+ if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
+ if(invMasspiKpLpi>0.) fhMassPtGreater2LpiTC->Fill(invMasspiKpLpi,0.5);
+ if(invMasspKpiLpi>0.) fhMassPtGreater2LpiTC->Fill(invMasspKpiLpi,0.5);
+ }else{
+ if(invMasspiKpLpi>0.) fhMassPtGreater2LpiTC->Fill(invMasspiKpLpi);
+ if(invMasspKpiLpi>0.) fhMassPtGreater2LpiTC->Fill(invMasspKpiLpi);
+ }
+ if(invMasspiKpKp>0. && invMasspKpiKp>0.){
+ if(invMasspiKpKp>0.) fhMassPtGreater2KpTC->Fill(invMasspiKpKp,0.5);
+ if(invMasspKpiKp>0.) fhMassPtGreater2KpTC->Fill(invMasspKpiKp,0.5);
+ }else{
+ if(invMasspiKpKp>0.) fhMassPtGreater2KpTC->Fill(invMasspiKpKp);
+ if(invMasspKpiKp>0.) fhMassPtGreater2KpTC->Fill(invMasspKpiKp);
+ }
}
- if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
- if(invMasspiKpLpi>0.) fhMassPtGreater2LpiTC->Fill(invMasspiKpLpi,0.5);
- if(invMasspKpiLpi>0.) fhMassPtGreater2LpiTC->Fill(invMasspKpiLpi,0.5);
- }else{
- if(invMasspiKpLpi>0.) fhMassPtGreater2LpiTC->Fill(invMasspiKpLpi);
- if(invMasspKpiLpi>0.) fhMassPtGreater2LpiTC->Fill(invMasspKpiLpi);
- }
- if(invMasspiKpKp>0. && invMasspKpiKp>0.){
- if(invMasspiKpKp>0.) fhMassPtGreater2KpTC->Fill(invMasspiKpKp,0.5);
- if(invMasspKpiKp>0.) fhMassPtGreater2KpTC->Fill(invMasspKpiKp,0.5);
- }else{
- if(invMasspiKpKp>0.) fhMassPtGreater2KpTC->Fill(invMasspiKpKp);
- if(invMasspKpiKp>0.) fhMassPtGreater2KpTC->Fill(invMasspKpiKp);
- }
- }
}
- if(iPtBin>=0){
- index=GetHistoIndex(iPtBin);
- if(invMasspiKp>0. && invMasspKpi>0.){
- if(invMasspiKp>0.) fMassHist[index]->Fill(invMasspiKp,0.5);
- if(invMasspKpi>0.) fMassHist[index]->Fill(invMasspKpi,0.5);
- }else{
- if(invMasspiKp>0.) fMassHist[index]->Fill(invMasspiKp);
- if(invMasspKpi>0.) fMassHist[index]->Fill(invMasspKpi);
- }
- if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
- if(invMasspiKpLpi>0.) fMassHistLpi[index]->Fill(invMasspiKpLpi,0.5);
- if(invMasspKpiLpi>0.) fMassHistLpi[index]->Fill(invMasspKpiLpi,0.5);
- }else{
- if(invMasspiKpLpi>0.) fMassHistLpi[index]->Fill(invMasspiKpLpi);
- if(invMasspKpiLpi>0.) fMassHistLpi[index]->Fill(invMasspKpiLpi);
- }
- if(invMasspiKpKp>0. && invMasspKpiKp>0.){
- if(invMasspiKpKp>0.) fMassHistKp[index]->Fill(invMasspiKpKp,0.5);
- if(invMasspKpiKp>0.) fMassHistKp[index]->Fill(invMasspKpiKp,0.5);
- }else{
- if(invMasspiKpKp>0.) fMassHistKp[index]->Fill(invMasspiKpKp);
- if(invMasspKpiKp>0.) fMassHistKp[index]->Fill(invMasspKpiKp);
- }
- if(passTightCuts>0){
- if(invMasspiKp>0. && invMasspKpi>0. && passTightCuts==3){
- if(invMasspiKp>0.) fMassHistTC[index]->Fill(invMasspiKp,0.5);
- if(invMasspKpi>0.) fMassHistTC[index]->Fill(invMasspKpi,0.5);
- }else{
- if(invMasspiKp>0. && passTightCuts==2) fMassHistTC[index]->Fill(invMasspiKp);
- if(invMasspKpi>0. && passTightCuts==1) fMassHistTC[index]->Fill(invMasspKpi);
- }
- if(invMasspiKpLpi>0. && invMasspKpiLpi>0. && passTightCuts==3){
- if(invMasspiKpLpi>0.) fMassHistLpiTC[index]->Fill(invMasspiKpLpi,0.5);
- if(invMasspKpiLpi>0.) fMassHistLpiTC[index]->Fill(invMasspKpiLpi,0.5);
- }else{
- if(invMasspiKpLpi>0. && passTightCuts==2) fMassHistLpiTC[index]->Fill(invMasspiKpLpi);
- if(invMasspKpiLpi>0.&& passTightCuts==1) fMassHistLpiTC[index]->Fill(invMasspKpiLpi);
- }
- if(invMasspiKpKp>0. && invMasspKpiKp>0. && passTightCuts==3){
- if(invMasspiKpKp>0.) fMassHistKpTC[index]->Fill(invMasspiKpKp,0.5);
- if(invMasspKpiKp>0.) fMassHistKpTC[index]->Fill(invMasspKpiKp,0.5);
- }else{
- if(invMasspiKpKp>0. && passTightCuts==2) fMassHistKpTC[index]->Fill(invMasspiKpKp);
- if(invMasspKpiKp>0.&& passTightCuts==1) fMassHistKpTC[index]->Fill(invMasspKpiKp);
- }
- }
- if(fReadMC){
- if(labDp>0) {
- index=GetSignalHistoIndex(iPtBin);
- if(invMasspiKp>0. && invMasspKpi>0.){
- if(invMasspiKp>0.) fMassHist[index]->Fill(invMasspiKp,0.5);
- if(invMasspKpi>0.) fMassHist[index]->Fill(invMasspKpi,0.5);
- }else{
- if(invMasspiKp>0.) fMassHist[index]->Fill(invMasspiKp);
- if(invMasspKpi>0.) fMassHist[index]->Fill(invMasspKpi);
- }
- if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
- if(invMasspiKpLpi>0.) fMassHistLpi[index]->Fill(invMasspiKpLpi,0.5);
- if(invMasspKpiLpi>0.) fMassHistLpi[index]->Fill(invMasspKpiLpi,0.5);
- }else{
- if(invMasspiKpLpi>0.) fMassHistLpi[index]->Fill(invMasspiKpLpi);
- if(invMasspKpiLpi>0.) fMassHistLpi[index]->Fill(invMasspKpiLpi);
- }
- if(invMasspiKpKp>0. && invMasspKpiKp>0.){
- if(invMasspiKpKp>0.) fMassHistKp[index]->Fill(invMasspiKpKp,0.5);
- if(invMasspKpiKp>0.) fMassHistKp[index]->Fill(invMasspKpiKp,0.5);
- }else{
- if(invMasspiKpKp>0.) fMassHistKp[index]->Fill(invMasspiKpKp);
- if(invMasspKpiKp>0.) fMassHistKp[index]->Fill(invMasspKpiKp);
- }
- if(passTightCuts>0){
- if(invMasspiKp>0. && invMasspKpi>0. && passTightCuts==3){
- if(invMasspiKp>0.) fMassHistTC[index]->Fill(invMasspiKp,0.5);
- if(invMasspKpi>0.) fMassHistTC[index]->Fill(invMasspKpi,0.5);
+ if(iPtBin>=0){
+ index=GetHistoIndex(iPtBin);
+ if(invMasspiKp>0. && invMasspKpi>0.){
+ if(invMasspiKp>0.) fMassHist[index]->Fill(invMasspiKp,0.5);
+ if(invMasspKpi>0.) fMassHist[index]->Fill(invMasspKpi,0.5);
}else{
- if(invMasspiKp>0. && passTightCuts==2) fMassHistTC[index]->Fill(invMasspiKp);
- if(invMasspKpi>0.&& passTightCuts==1) fMassHistTC[index]->Fill(invMasspKpi);
+ if(invMasspiKp>0.) fMassHist[index]->Fill(invMasspiKp);
+ if(invMasspKpi>0.) fMassHist[index]->Fill(invMasspKpi);
}
- if(invMasspiKpLpi>0. && invMasspKpiLpi>0. && passTightCuts==3){
- if(invMasspiKpLpi>0.) fMassHistLpiTC[index]->Fill(invMasspiKpLpi,0.5);
- if(invMasspKpiLpi>0.) fMassHistLpiTC[index]->Fill(invMasspKpiLpi,0.5);
- }else{
- if(invMasspiKpLpi>0. && passTightCuts==2) fMassHistLpiTC[index]->Fill(invMasspiKpLpi);
- if(invMasspKpiLpi>0.&& passTightCuts==1) fMassHistLpiTC[index]->Fill(invMasspKpiLpi);
- }
- if(invMasspiKpKp>0. && invMasspKpiKp>0. && passTightCuts==3){
- if(invMasspiKpKp>0.) fMassHistKpTC[index]->Fill(invMasspiKpKp,0.5);
- if(invMasspKpiKp>0.) fMassHistKpTC[index]->Fill(invMasspKpiKp,0.5);
- }else{
- if(invMasspiKpKp>0. && passTightCuts==2) fMassHistKpTC[index]->Fill(invMasspiKpKp);
- if(invMasspKpiKp>0.&& passTightCuts==1) fMassHistKpTC[index]->Fill(invMasspKpiKp);
- }
- }
+ if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
+ if(invMasspiKpLpi>0.) fMassHistLpi[index]->Fill(invMasspiKpLpi,0.5);
+ if(invMasspKpiLpi>0.) fMassHistLpi[index]->Fill(invMasspKpiLpi,0.5);
+ }else{
+ if(invMasspiKpLpi>0.) fMassHistLpi[index]->Fill(invMasspiKpLpi);
+ if(invMasspKpiLpi>0.) fMassHistLpi[index]->Fill(invMasspKpiLpi);
+ }
+ if(invMasspiKpKp>0. && invMasspKpiKp>0.){
+ if(invMasspiKpKp>0.) fMassHistKp[index]->Fill(invMasspiKpKp,0.5);
+ if(invMasspKpiKp>0.) fMassHistKp[index]->Fill(invMasspKpiKp,0.5);
+ }else{
+ if(invMasspiKpKp>0.) fMassHistKp[index]->Fill(invMasspiKpKp);
+ if(invMasspKpiKp>0.) fMassHistKp[index]->Fill(invMasspKpiKp);
+ }
+ if(passTightCuts>0){
+ if(invMasspiKp>0. && invMasspKpi>0. && passTightCuts==3){
+ if(invMasspiKp>0.) fMassHistTC[index]->Fill(invMasspiKp,0.5);
+ if(invMasspKpi>0.) fMassHistTC[index]->Fill(invMasspKpi,0.5);
+ }else{
+ if(invMasspiKp>0. && passTightCuts==2) fMassHistTC[index]->Fill(invMasspiKp);
+ if(invMasspKpi>0. && passTightCuts==1) fMassHistTC[index]->Fill(invMasspKpi);
+ }
+ if(invMasspiKpLpi>0. && invMasspKpiLpi>0. && passTightCuts==3){
+ if(invMasspiKpLpi>0.) fMassHistLpiTC[index]->Fill(invMasspiKpLpi,0.5);
+ if(invMasspKpiLpi>0.) fMassHistLpiTC[index]->Fill(invMasspKpiLpi,0.5);
+ }else{
+ if(invMasspiKpLpi>0. && passTightCuts==2) fMassHistLpiTC[index]->Fill(invMasspiKpLpi);
+ if(invMasspKpiLpi>0.&& passTightCuts==1) fMassHistLpiTC[index]->Fill(invMasspKpiLpi);
+ }
+ if(invMasspiKpKp>0. && invMasspKpiKp>0. && passTightCuts==3){
+ if(invMasspiKpKp>0.) fMassHistKpTC[index]->Fill(invMasspiKpKp,0.5);
+ if(invMasspKpiKp>0.) fMassHistKpTC[index]->Fill(invMasspKpiKp,0.5);
+ }else{
+ if(invMasspiKpKp>0. && passTightCuts==2) fMassHistKpTC[index]->Fill(invMasspiKpKp);
+ if(invMasspKpiKp>0.&& passTightCuts==1) fMassHistKpTC[index]->Fill(invMasspKpiKp);
+ }
+ }
+ if(fReadMC){
+ if(labDp>0) {
+ index=GetSignalHistoIndex(iPtBin);
+ if(invMasspiKp>0. && invMasspKpi>0.){
+ if(invMasspiKp>0.) fMassHist[index]->Fill(invMasspiKp,0.5);
+ if(invMasspKpi>0.) fMassHist[index]->Fill(invMasspKpi,0.5);
+ }else{
+ if(invMasspiKp>0.) fMassHist[index]->Fill(invMasspiKp);
+ if(invMasspKpi>0.) fMassHist[index]->Fill(invMasspKpi);
+ }
+ if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
+ if(invMasspiKpLpi>0.) fMassHistLpi[index]->Fill(invMasspiKpLpi,0.5);
+ if(invMasspKpiLpi>0.) fMassHistLpi[index]->Fill(invMasspKpiLpi,0.5);
+ }else{
+ if(invMasspiKpLpi>0.) fMassHistLpi[index]->Fill(invMasspiKpLpi);
+ if(invMasspKpiLpi>0.) fMassHistLpi[index]->Fill(invMasspKpiLpi);
+ }
+ if(invMasspiKpKp>0. && invMasspKpiKp>0.){
+ if(invMasspiKpKp>0.) fMassHistKp[index]->Fill(invMasspiKpKp,0.5);
+ if(invMasspKpiKp>0.) fMassHistKp[index]->Fill(invMasspKpiKp,0.5);
+ }else{
+ if(invMasspiKpKp>0.) fMassHistKp[index]->Fill(invMasspiKpKp);
+ if(invMasspKpiKp>0.) fMassHistKp[index]->Fill(invMasspKpiKp);
+ }
+ if(passTightCuts>0){
+ if(invMasspiKp>0. && invMasspKpi>0. && passTightCuts==3){
+ if(invMasspiKp>0.) fMassHistTC[index]->Fill(invMasspiKp,0.5);
+ if(invMasspKpi>0.) fMassHistTC[index]->Fill(invMasspKpi,0.5);
+ }else{
+ if(invMasspiKp>0. && passTightCuts==2) fMassHistTC[index]->Fill(invMasspiKp);
+ if(invMasspKpi>0.&& passTightCuts==1) fMassHistTC[index]->Fill(invMasspKpi);
+ }
+ if(invMasspiKpLpi>0. && invMasspKpiLpi>0. && passTightCuts==3){
+ if(invMasspiKpLpi>0.) fMassHistLpiTC[index]->Fill(invMasspiKpLpi,0.5);
+ if(invMasspKpiLpi>0.) fMassHistLpiTC[index]->Fill(invMasspKpiLpi,0.5);
+ }else{
+ if(invMasspiKpLpi>0. && passTightCuts==2) fMassHistLpiTC[index]->Fill(invMasspiKpLpi);
+ if(invMasspKpiLpi>0.&& passTightCuts==1) fMassHistLpiTC[index]->Fill(invMasspKpiLpi);
+ }
+ if(invMasspiKpKp>0. && invMasspKpiKp>0. && passTightCuts==3){
+ if(invMasspiKpKp>0.) fMassHistKpTC[index]->Fill(invMasspiKpKp,0.5);
+ if(invMasspKpiKp>0.) fMassHistKpTC[index]->Fill(invMasspKpiKp,0.5);
+ }else{
+ if(invMasspiKpKp>0. && passTightCuts==2) fMassHistKpTC[index]->Fill(invMasspiKpKp);
+ if(invMasspKpiKp>0.&& passTightCuts==1) fMassHistKpTC[index]->Fill(invMasspKpiKp);
+ }
+ }
- }else{
- index=GetBackgroundHistoIndex(iPtBin);
- if(invMasspiKp>0. && invMasspKpi>0.){
- fMassHist[index]->Fill(invMasspiKp,0.5);
- fMassHist[index]->Fill(invMasspKpi,0.5);
- }else{
- if(invMasspiKp>0.) fMassHist[index]->Fill(invMasspiKp);
- if(invMasspKpi>0.) fMassHist[index]->Fill(invMasspKpi);
- }
- if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
- if(invMasspiKpLpi>0.) fMassHistLpi[index]->Fill(invMasspiKpLpi,0.5);
- if(invMasspKpiLpi>0.) fMassHistLpi[index]->Fill(invMasspKpiLpi,0.5);
- }else{
- if(invMasspiKpLpi>0.) fMassHistLpi[index]->Fill(invMasspiKpLpi);
- if(invMasspKpiLpi>0.) fMassHistLpi[index]->Fill(invMasspKpiLpi);
+ }else{
+ index=GetBackgroundHistoIndex(iPtBin);
+ if(invMasspiKp>0. && invMasspKpi>0.){
+ fMassHist[index]->Fill(invMasspiKp,0.5);
+ fMassHist[index]->Fill(invMasspKpi,0.5);
+ }else{
+ if(invMasspiKp>0.) fMassHist[index]->Fill(invMasspiKp);
+ if(invMasspKpi>0.) fMassHist[index]->Fill(invMasspKpi);
+ }
+ if(invMasspiKpLpi>0. && invMasspKpiLpi>0.){
+ if(invMasspiKpLpi>0.) fMassHistLpi[index]->Fill(invMasspiKpLpi,0.5);
+ if(invMasspKpiLpi>0.) fMassHistLpi[index]->Fill(invMasspKpiLpi,0.5);
+ }else{
+ if(invMasspiKpLpi>0.) fMassHistLpi[index]->Fill(invMasspiKpLpi);
+ if(invMasspKpiLpi>0.) fMassHistLpi[index]->Fill(invMasspKpiLpi);
+ }
+ if(invMasspiKpKp>0. && invMasspKpiKp>0.){
+ if(invMasspiKpKp>0.) fMassHistKp[index]->Fill(invMasspiKpKp,0.5);
+ if(invMasspKpiKp>0.) fMassHistKp[index]->Fill(invMasspKpiKp,0.5);
+ }else{
+ if(invMasspiKpKp>0.) fMassHistKp[index]->Fill(invMasspiKpKp);
+ if(invMasspKpiKp>0.) fMassHistKp[index]->Fill(invMasspKpiKp);
+ }
+ if(invMasspiKp>0. && invMasspKpi>0. && passTightCuts==3){
+ fMassHistTC[index]->Fill(invMasspiKp,0.5);
+ fMassHistTC[index]->Fill(invMasspKpi,0.5);
+ }else{
+ if(invMasspiKp>0. && passTightCuts==2) fMassHistTC[index]->Fill(invMasspiKp);
+ if(invMasspKpi>0. && passTightCuts==1) fMassHistTC[index]->Fill(invMasspKpi);
+ }
+ if(invMasspiKpLpi>0. && invMasspKpiLpi>0. && passTightCuts==3){
+ if(invMasspiKpLpi>0.) fMassHistLpiTC[index]->Fill(invMasspiKpLpi,0.5);
+ if(invMasspKpiLpi>0.) fMassHistLpiTC[index]->Fill(invMasspKpiLpi,0.5);
+ }else{
+ if(invMasspiKpLpi>0. && passTightCuts==2) fMassHistLpiTC[index]->Fill(invMasspiKpLpi);
+ if(invMasspKpiLpi>0.&& passTightCuts==1) fMassHistLpiTC[index]->Fill(invMasspKpiLpi);
+ }
+ if(invMasspiKpKp>0. && invMasspKpiKp>0. && passTightCuts==3){
+ if(invMasspiKpKp>0.) fMassHistKpTC[index]->Fill(invMasspiKpKp,0.5);
+ if(invMasspKpiKp>0.) fMassHistKpTC[index]->Fill(invMasspKpiKp,0.5);
+ }else{
+ if(invMasspiKpKp>0. && passTightCuts==2) fMassHistKpTC[index]->Fill(invMasspiKpKp);
+ if(invMasspKpiKp>0.&& passTightCuts==1) fMassHistKpTC[index]->Fill(invMasspKpiKp);
+ }
+ }
+
+ }
}
- if(invMasspiKpKp>0. && invMasspKpiKp>0.){
- if(invMasspiKpKp>0.) fMassHistKp[index]->Fill(invMasspiKpKp,0.5);
- if(invMasspKpiKp>0.) fMassHistKp[index]->Fill(invMasspKpiKp,0.5);
- }else{
- if(invMasspiKpKp>0.) fMassHistKp[index]->Fill(invMasspiKpKp);
- if(invMasspKpiKp>0.) fMassHistKp[index]->Fill(invMasspKpiKp);
- }
- if(invMasspiKp>0. && invMasspKpi>0. && passTightCuts==3){
- fMassHistTC[index]->Fill(invMasspiKp,0.5);
- fMassHistTC[index]->Fill(invMasspKpi,0.5);
- }else{
- if(invMasspiKp>0. && passTightCuts==2) fMassHistTC[index]->Fill(invMasspiKp);
- if(invMasspKpi>0. && passTightCuts==1) fMassHistTC[index]->Fill(invMasspKpi);
- }
- if(invMasspiKpLpi>0. && invMasspKpiLpi>0. && passTightCuts==3){
- if(invMasspiKpLpi>0.) fMassHistLpiTC[index]->Fill(invMasspiKpLpi,0.5);
- if(invMasspKpiLpi>0.) fMassHistLpiTC[index]->Fill(invMasspKpiLpi,0.5);
- }else{
- if(invMasspiKpLpi>0. && passTightCuts==2) fMassHistLpiTC[index]->Fill(invMasspiKpLpi);
- if(invMasspKpiLpi>0.&& passTightCuts==1) fMassHistLpiTC[index]->Fill(invMasspKpiLpi);
- }
- if(invMasspiKpKp>0. && invMasspKpiKp>0. && passTightCuts==3){
- if(invMasspiKpKp>0.) fMassHistKpTC[index]->Fill(invMasspiKpKp,0.5);
- if(invMasspKpiKp>0.) fMassHistKpTC[index]->Fill(invMasspKpiKp,0.5);
- }else{
- if(invMasspiKpKp>0. && passTightCuts==2) fMassHistKpTC[index]->Fill(invMasspiKpKp);
- if(invMasspKpiKp>0.&& passTightCuts==1) fMassHistKpTC[index]->Fill(invMasspKpiKp);
- }
- }
-
}
- }
-}
- return;
+ return;
}
//-----------------------
-void AliAnalysisTaskSELambdac::FillVarHists(AliAODRecoDecayHF3Prong *part, TClonesArray *arrMC, AliRDHFCutsLctopKpi *cuts, TList *listout,AliAODEvent* aod){
+void AliAnalysisTaskSELambdac::FillVarHists(AliAODRecoDecayHF3Prong *part,
+ TClonesArray *arrMC,
+ AliRDHFCutsLctopKpi *cuts,
+ /* TList *listout,*/
+ AliAODEvent* aod) {
//
// function used in UserExec to fill variable histograms analysing MC
//
+ TString fillthis="";
+
+ Double_t mPDG=TDatabasePDG::Instance()->GetParticle(4122)->Mass();
+ Double_t invmasscut=0.05;
Int_t pdgDgLctopKpi[3]={2212,321,211};
Int_t lab=-9999;
- if(fReadMC) lab=part->MatchToMC(4122,arrMC,3,pdgDgLctopKpi); //return MC particle label if the array corresponds to a Lc, -1 if not (cf. AliAODRecoDecay.cxx)
- Int_t isSelectedPID=cuts->IsSelectedPID(part); //0 rejected,1 Lc -> p K- pi+,2 Lc -> pi+ K- p, (K at center because different sign is there)
- //3 both (it should never happen...)
+ if(fReadMC)
+ lab=part->MatchToMC(4122,arrMC,3,pdgDgLctopKpi); //return MC particle label if the array corresponds to a Lc, -1 if not (cf. AliAODRecoDecay.cxx)
+
+ Int_t isSelectedPID=cuts->IsSelectedPID(part); // 0 rejected, 1 Lc -> p K- pi+ (K at center because different sign is there),
+ // 2 Lc -> pi+ K- p (K at center because different sign is there),
+ // 3 both (it should never happen...)
if (isSelectedPID==0)fNentries->Fill(7);
if (isSelectedPID==1)fNentries->Fill(8);
if (isSelectedPID==2)fNentries->Fill(9);
if (isSelectedPID==3)fNentries->Fill(10);
-
- Double_t mPDG=TDatabasePDG::Instance()->GetParticle(4122)->Mass();
- Double_t minvLcpKpi = part->InvMassLcpKpi();
- Double_t minvLcpiKp = part->InvMassLcpiKp();
-
- Double_t invmasscut=0.05;
-
- TString fillthis="";
-
-
AliAODTrack *prong0=(AliAODTrack*)part->GetDaughter(0);
AliAODTrack *prong1=(AliAODTrack*)part->GetDaughter(1);
AliAODTrack *prong2=(AliAODTrack*)part->GetDaughter(2);
return;
}
+ Double_t minvLcpKpi = part->InvMassLcpKpi();
+ Double_t minvLcpiKp = part->InvMassLcpiKp();
+
//check pdg of the prongs
- Int_t labprong[3];
+ Int_t labprong[3]={-1,-1,-1};
if(fReadMC){
labprong[0]=prong0->GetLabel();
labprong[1]=prong1->GetLabel();
labprong[2]=prong2->GetLabel();
}
- AliAODMCParticle *mcprong=0;
- Int_t pdgProng[3]={0,0,0};
- Int_t pdgProngID[3]={0,0,0};
- if(fReadMC) {
- for (Int_t iprong=0;iprong<3;iprong++){
+ AliAODMCParticle *mcprong=0;
+ Int_t pdgProngMC[3]={-1,-1,-1};
+ if(fReadMC) {
+ for (Int_t iprong=0;iprong<3;iprong++){
if(labprong[iprong]<0) continue;
- mcprong= (AliAODMCParticle*)arrMC->At(labprong[iprong]);
- pdgProng[iprong]=TMath::Abs(mcprong->GetPdgCode());
- }
- if(isSelectedPID>0 && fReadMC){
- pdgProngID[1]=321;
- if(isSelectedPID==1) {pdgProngID[0]=2212;pdgProngID[2]=211;}
- if(isSelectedPID==2) {pdgProngID[0]=211;pdgProngID[2]=2212;}
- }
- } else {
- if (isSelectedPID>0) {
- pdgProng[1]=321;
- if(isSelectedPID==1) {pdgProng[0]=2212;pdgProng[2]=211;}
- if(isSelectedPID==2) {pdgProng[0]=211;pdgProng[2]=2212;}
+ mcprong = (AliAODMCParticle*)arrMC->At(labprong[iprong]);
+ pdgProngMC[iprong]=TMath::Abs(mcprong->GetPdgCode());
+ }
+ }
+
+ Int_t pdgProngPID[3]={-1,-1,-1};
+ if(isSelectedPID>0){
+ pdgProngPID[1]=321;
+ if(isSelectedPID==1) {pdgProngPID[0]=2212;pdgProngPID[2]=211;}
+ if(isSelectedPID==2) {pdgProngPID[0]=211;pdgProngPID[2]=2212;}
+ }
+
+ //fill hRealTot ---> PID efficiency denominators
+ cuts->SetUsePID(kFALSE); //PAOLA
+ Int_t selectionCandidate=cuts->IsSelected(part,AliRDHFCuts::kCandidate,aod);//PAOLA
+ if(fReadMC && selectionCandidate>0) { // 3prongs candidate x Lc (no PID)
+ if(lab>0){ // Signal
+ for (Int_t iprong=0; iprong<3; iprong++) {
+ switch (pdgProngMC[iprong]) {
+ case 2212:
+ fillthis="hpRealTot";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ break;
+ case 321:
+ fillthis="hKRealTot";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ break;
+ case 211:
+ fillthis="hpiRealTot";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ break;
+ default:
+ break;
+ }
+ }
+ } else { // bkg
+ for (Int_t iprong=0; iprong<3; iprong++) {
+ switch (pdgProngMC[iprong]) {
+ case 2212:
+ fillthis="hbpRealTot";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ break;
+ case 321:
+ fillthis="hbKRealTot";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ break;
+ case 211:
+ fillthis="hbpiRealTot";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ break;
+ default:
+ break;
+ }
}
}
+ }
- Int_t selection=cuts->IsSelected(part,AliRDHFCuts::kCandidate,aod);
- if((lab>=0 && fReadMC) || (!fReadMC && (isSelectedPID>0)) ){ //signal (from MC or PID)
+ cuts->SetUsePID(kTRUE); //PAOLA
+ Int_t selection=cuts->IsSelected(part,AliRDHFCuts::kCandidate,aod);
- fillthis="hMass";
+ if ( (lab>0 && fReadMC) || // signal X MC
+ (isSelectedPID>0 && !fReadMC) ) { // signal+bkg X real data
- if ((fReadMC && ((AliAODMCParticle*)arrMC->At(lab))->GetPdgCode() == 4122)
- || (!fReadMC && (isSelectedPID>0 && part->Charge()>0))){ //Lc
- if ( (pdgProng[0]==2212) && (pdgProng[1]==321) && (pdgProng[2]==211) ) ((TH1F*)listout->FindObject(fillthis))->Fill(minvLcpKpi);
- else if ( (pdgProng[0]==211) && (pdgProng[1]==321) && (pdgProng[2]==2212) ) ((TH1F*)listout->FindObject(fillthis))->Fill(minvLcpiKp);
+ fillthis="hMass";
+ if ( (fReadMC && ((AliAODMCParticle*)arrMC->At(lab))->GetPdgCode() == 4122) ||
+ (!fReadMC && (isSelectedPID>0 && part->Charge()>0)) ) { // Lc
+ if ( (pdgProngPID[0]==2212) && (pdgProngPID[1]==321) && (pdgProngPID[2]==211) )
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(minvLcpKpi); // signal+bkg X MC and real data
+ else if ( (pdgProngPID[0]==211) && (pdgProngPID[1]==321) && (pdgProngPID[2]==2212) )
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(minvLcpiKp); // signal+bkg X MC and real data
}
- else if ((fReadMC && ((AliAODMCParticle*)arrMC->At(lab))->GetPdgCode() == -4122)
- || (!fReadMC && (isSelectedPID>0 && part->Charge()<0))){ //anti-Lc
- if ( (pdgProng[0]==2212) && (pdgProng[1]==321) && (pdgProng[2]==211) ) ((TH1F*)listout->FindObject(fillthis))->Fill(minvLcpKpi);
- else if ( (pdgProng[0]==211) && (pdgProng[1]==321) && (pdgProng[2]==2212) ) ((TH1F*)listout->FindObject(fillthis))->Fill(minvLcpiKp);
+ else if ( (fReadMC && ((AliAODMCParticle*)arrMC->At(lab))->GetPdgCode() == -4122) ||
+ (!fReadMC && (isSelectedPID>0 && part->Charge()<0)) ) { // anti-Lc
+ if ( (pdgProngPID[0]==2212) && (pdgProngPID[1]==321) && (pdgProngPID[2]==211) )
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(minvLcpKpi); // signal+bkg X MC and real data
+ else if ( (pdgProngPID[0]==211) && (pdgProngPID[1]==321) && (pdgProngPID[2]==2212) )
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(minvLcpiKp); // signal+bkg X MC and real data
}
-
- //apply cut on invariant mass on the pair
- if(selection>0){
+ if (selection>0) { // 3prongs candidate x Lc (yes PID)
Double_t ptmax=0.;
+
for (Int_t iprong=0; iprong<3; iprong++) {
- AliAODTrack *prong=(AliAODTrack*)part->GetDaughter(iprong);
- if (fReadMC) {
- labprong[iprong]=prong->GetLabel();
- if(pdgProngID[iprong]==2212) {
- fillthis="hpIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- if(pdgProngID[iprong]==321) {
- fillthis="hKIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- if(pdgProngID[iprong]==211) {
- fillthis="hpiIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- }
- AliAODPid *pidObjtrk=prong->GetDetPid();
- AliAODPidHF *pidObj=new AliAODPidHF();
+ if (part->PtProng(iprong)>ptmax) ptmax=part->PtProng(iprong);
+
+ AliAODTrack *prong = (AliAODTrack*)part->GetDaughter(iprong);
+ AliAODPid *pidObjtrk = (AliAODPid*)prong->GetDetPid();
+ AliAODPidHF *pidObj = (AliAODPidHF*)cuts->GetPidHF();
Bool_t hasTOF=pidObj->CheckStatus(prong,"TOF");
Bool_t hasTPC=pidObj->CheckStatus(prong,"TPC");
- delete pidObj;
Double_t tofSignal=0.;
Double_t dedxTPC=0.;
Double_t momTOF=0.;
Double_t momTPC=0.;
if(hasTOF) {
- momTOF = prong->P();
- tofSignal=pidObjtrk->GetTOFsignal();
+ momTOF = prong->P();
+ tofSignal=pidObjtrk->GetTOFsignal();
}
if(hasTPC) {
- momTPC = pidObjtrk->GetTPCmomentum();
- dedxTPC=pidObjtrk->GetTPCsignal();
+ momTPC = pidObjtrk->GetTPCmomentum();
+ dedxTPC=pidObjtrk->GetTPCsignal();
}
- switch (pdgProng[iprong]) {
+ switch (pdgProngPID[iprong]) {
case 2212:
- fillthis="hpTOFSignal";
- ((TH1F*)listout->FindObject(fillthis))->Fill(tofSignal);
- fillthis="hpTPCSignal";
- ((TH1F*)listout->FindObject(fillthis))->Fill(dedxTPC);
- fillthis="hpptProng";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- fillthis="hpd0Prong";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->Getd0Prong(iprong));
- fillthis="hpSignalVspTPC";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTPC,dedxTPC);
- fillthis="hpSignalVspTOF";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTOF,tofSignal);
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->Getd0Prong(iprong));
- AliPID::EParticleType typep;
- typep=AliPID::EParticleType(4);
- if(hasTPC) {
- Double_t nsigmap = fUtilPid->NumberOfSigmasTPC(prong,typep);
- fillthis="hpSigmaVspTPC";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTPC,nsigmap);
- }
- if(hasTOF){
- Double_t nsigma=fUtilPid->NumberOfSigmasTOF(prong,typep);
- fillthis="hpSigmaVspTOF";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTOF,nsigma);
+ fillthis="hpTOFSignal";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(tofSignal);
+ fillthis="hpTPCSignal";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(dedxTPC);
+ fillthis="hpptProng";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ fillthis="hpd0Prong";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->Getd0Prong(iprong));
+ fillthis="hpSignalVspTPC";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTPC,dedxTPC);
+ fillthis="hpSignalVspTOF";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTOF,tofSignal);
+ AliPID::EParticleType typep;
+ typep=AliPID::EParticleType(4);
+ if(hasTPC) {
+ //Double_t nsigmap = fUtilPid->NumberOfSigmasTPC(prong,typep);
+ Double_t nsigmap = fPIDResponse->NumberOfSigmasTPC(prong,typep);
+ fillthis="hpSigmaVspTPC";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTPC,nsigmap);
+ }
+ if(hasTOF){
+ //Double_t nsigma=fUtilPid->NumberOfSigmasTOF(prong,typep);
+ Double_t nsigma=fPIDResponse->NumberOfSigmasTOF(prong,typep);
+ fillthis="hpSigmaVspTOF";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTOF,nsigma);
+ }
+
+ if (fReadMC) { // fill hpIDTot ---> PID contamination denominator
+ // hIDGood, hnoID ---> PID numerators
+ fillthis="hpIDTot";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ if(pdgProngMC[iprong]==2212) { // id protons
+ fillthis="hpIDGood";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
}
- if(fReadMC){
- // real protons
- fillthis="hpRealTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- // id protons
- if(pdgProngID[iprong]==2212) {
- fillthis="hpIDGood";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- // misidentified kaons, pions
- if(pdgProngID[iprong]==321) {
- fillthis="hKnonIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- if(pdgProngID[iprong]==211) {
- fillthis="hpinonIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
+ else { // misidentified electrons, muons, pions and kaons
+ fillthis="hnopIDp";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
}
+ }
+ break;
- break;
case 321:
- fillthis="hKTOFSignal";
- ((TH1F*)listout->FindObject(fillthis))->Fill(tofSignal);
- fillthis="hKTPCSignal";
- ((TH1F*)listout->FindObject(fillthis))->Fill(dedxTPC);
- fillthis="hKptProng";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- fillthis="hKd0Prong";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->Getd0Prong(iprong));
- fillthis="hKSignalVspTPC";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTPC,dedxTPC);
- fillthis="hKSignalVspTOF";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTOF,tofSignal);
- AliPID::EParticleType typek;
- typek=AliPID::EParticleType(3);
- if(hasTPC) {
- Double_t nsigmap = fUtilPid->NumberOfSigmasTPC(prong,typek);
- fillthis="hKSigmaVspTPC";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTPC,nsigmap);
+ fillthis="hKTOFSignal";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(tofSignal);
+ fillthis="hKTPCSignal";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(dedxTPC);
+ fillthis="hKptProng";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ fillthis="hKd0Prong";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->Getd0Prong(iprong));
+ fillthis="hKSignalVspTPC";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTPC,dedxTPC);
+ fillthis="hKSignalVspTOF";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTOF,tofSignal);
+ AliPID::EParticleType typek;
+ typek=AliPID::EParticleType(3);
+ if(hasTPC) {
+ //Double_t nsigmap = fUtilPid->NumberOfSigmasTPC(prong,typek);
+ Double_t nsigmap = fPIDResponse->NumberOfSigmasTPC(prong,typek);
+ fillthis="hKSigmaVspTPC";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTPC,nsigmap);
+ }
+ if(hasTOF){
+ //Double_t nsigma = fUtilPid->NumberOfSigmasTOF(prong,typek);
+ Double_t nsigma=fPIDResponse->NumberOfSigmasTOF(prong,typek); // protection against 'old' AODs...
+ fillthis="hKSigmaVspTOF";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTOF,nsigma);
+ }
+
+ if (fReadMC) { // fill hKIDTot ---> PID contamination denominator
+ // hIDGood, hnoID ---> PID numerators
+ fillthis="hKIDTot";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ if(pdgProngMC[iprong]==321) { // id kaons
+ fillthis="hKIDGood";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
}
- if(hasTOF){
- Double_t nsigma = fUtilPid->NumberOfSigmasTOF(prong,typek);
- fillthis="hKSigmaVspTOF";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTOF,nsigma);
+ else { // misidentified electrons, muons, pions and protons
+ fillthis="hnokIDk";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
}
+ }
+ break;
- if(fReadMC){
- // real kaons
- fillthis="hKRealTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- // id kaons
- if(pdgProngID[iprong]==321) {
- fillthis="hKIDGood";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- // misidentified protons, pions
- if(pdgProngID[iprong]==2212) {
- fillthis="hpnonIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- if(pdgProngID[iprong]==211) {
- fillthis="hpinonIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- }
- break;
case 211:
- fillthis="hpiTOFSignal";
- ((TH1F*)listout->FindObject(fillthis))->Fill(tofSignal);
- fillthis="hpiTPCSignal";
- ((TH1F*)listout->FindObject(fillthis))->Fill(dedxTPC);
- fillthis="hpiptProng";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- fillthis="hpid0Prong";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->Getd0Prong(iprong));
- fillthis="hpiSignalVspTPC";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTPC,dedxTPC);
- fillthis="hpiSignalVspTOF";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTOF,tofSignal);
- AliPID::EParticleType typepi;
- typepi=AliPID::EParticleType(2);
- if(hasTPC) {
- Double_t nsigmap = fUtilPid->NumberOfSigmasTPC(prong,typepi);
+ fillthis="hpiTOFSignal";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(tofSignal);
+ fillthis="hpiTPCSignal";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(dedxTPC);
+ fillthis="hpiptProng";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ fillthis="hpid0Prong";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->Getd0Prong(iprong));
+ fillthis="hpiSignalVspTPC";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTPC,dedxTPC);
+ fillthis="hpiSignalVspTOF";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTOF,tofSignal);
+ AliPID::EParticleType typepi;
+ typepi=AliPID::EParticleType(2);
+ if(hasTPC) {
+ //Double_t nsigmap = fUtilPid->NumberOfSigmasTPC(prong,typepi);
+ Double_t nsigmap = fPIDResponse->NumberOfSigmasTPC(prong,typepi);
fillthis="hpiSigmaVspTPC";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTPC,nsigmap);
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTPC,nsigmap);
+ }
+ if(hasTOF){
+ //Double_t nsigma = fUtilPid->NumberOfSigmasTOF(prong,typepi);
+ Double_t nsigma=fPIDResponse->NumberOfSigmasTOF(prong,typepi); // protection against 'old' AODs...
+ fillthis="hpiSigmaVspTOF";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTOF,nsigma);
+ }
+
+ if (fReadMC) { // fill hpiIDTot ---> PID contamination denominator
+ // hIDGood, hnoID ---> PID numerators
+ fillthis="hpiIDTot";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ if(pdgProngMC[iprong]==211) { // id pions
+ fillthis="hpiIDGood";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
}
- if(hasTOF){
- Double_t nsigma = fUtilPid->NumberOfSigmasTOF(prong,typepi);
- fillthis="hpiSigmaVspTOF";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTOF,nsigma);
- }
-
- if(fReadMC){
- // real pions
- fillthis="hpiRealTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- // id pions
- if(pdgProngID[iprong]==211) {
- fillthis="hpiIDGood";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- // misidentified protons, kaons
- if(pdgProngID[iprong]==2212) {
- fillthis="hpnonIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- if(pdgProngID[iprong]==321) {
- fillthis="hKnonIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
+ else { // misidentified electrons, muons, kaons and protons
+ fillthis="hnopiIDpi";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
}
- break;
+ }
+ break;
default:
break;
}
- if(part->PtProng(iprong)>ptmax)ptmax=part->PtProng(iprong);
} //end loop on prongs
- //now histograms where we don't need to check identity
+
+ //now histograms where we don't need to check identity
fillthis = "hDist12toPrim";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->GetDist12toPrim());
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->GetDist23toPrim());
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->GetDist12toPrim());
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->GetDist23toPrim());
fillthis = "hSigmaVert";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->GetSigmaVert());
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->GetSigmaVert());
fillthis = "hDCAs";
Double_t dcas[3];
part->GetDCAs(dcas);
- for (Int_t idca=0;idca<3;idca++) ((TH1F*)listout->FindObject(fillthis))->Fill(dcas[idca]);
+ for (Int_t idca=0;idca<3;idca++) ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(dcas[idca]);
fillthis = "hCosPointingAngle";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->CosPointingAngle());
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->CosPointingAngle());
fillthis = "hDecayLength";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->DecayLength());
- Double_t sum2=part->Getd0Prong(0)*part->Getd0Prong(0)+part->Getd0Prong(1)*part->Getd0Prong(1)+part->Getd0Prong(2)*part->Getd0Prong(2);
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->DecayLength());
+ Double_t sum2=part->Getd0Prong(0)*part->Getd0Prong(0)+
+ part->Getd0Prong(1)*part->Getd0Prong(1)+
+ part->Getd0Prong(2)*part->Getd0Prong(2);
fillthis = "hSum2";
- ((TH1F*)listout->FindObject(fillthis))->Fill(sum2);
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(sum2);
fillthis = "hptmax";
- ((TH1F*)listout->FindObject(fillthis))->Fill(ptmax);
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(ptmax);
fillthis="hLcpt";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->Pt());
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->Pt());
- } //end mass cut
+ } // end if (selection)
-
- } else if( (lab<0) && fReadMC) { // background **ONLY MC**
+ } else if( lab<=0 && fReadMC) { // bkg x MC
fillthis="hbMass";
- if (part->Charge()>0){ //Lc background
- if ( (pdgProng[0]==2212) && (pdgProng[1]==321) && (pdgProng[2]==211) ) ((TH1F*)listout->FindObject(fillthis))->Fill(minvLcpKpi);
- else if ( (pdgProng[0]==211) && (pdgProng[1]==321) && (pdgProng[2]==2212) ) ((TH1F*)listout->FindObject(fillthis))->Fill(minvLcpiKp);
- }
- else if (part->Charge()<0){ //anti-Lc background
- if ( (pdgProng[0]==2212) && (pdgProng[1]==321) && (pdgProng[2]==211) ) ((TH1F*)listout->FindObject(fillthis))->Fill(minvLcpKpi);
- else if ( (pdgProng[0]==211) && (pdgProng[1]==321) && (pdgProng[2]==2212) ) ((TH1F*)listout->FindObject(fillthis))->Fill(minvLcpiKp);
+ if (part->Charge()>0) { //Lc background
+ if ( (pdgProngPID[0]==2212) && (pdgProngPID[1]==321) && (pdgProngPID[2]==211) )
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(minvLcpKpi); // bkg X MC
+ else if ( (pdgProngPID[0]==211) && (pdgProngPID[1]==321) && (pdgProngPID[2]==2212) )
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(minvLcpiKp); // bkg X MC
+ }
+ else if (part->Charge()<0){ //anti-Lc background
+ if ( (pdgProngPID[0]==2212) && (pdgProngPID[1]==321) && (pdgProngPID[2]==211) )
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(minvLcpKpi); // bkg X MC
+ else if ( (pdgProngPID[0]==211) && (pdgProngPID[1]==321) && (pdgProngPID[2]==2212) )
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(minvLcpiKp); // bkg X MC
}
//apply cut on invariant mass on the pair
- if(TMath::Abs(minvLcpKpi-mPDG)<invmasscut || TMath::Abs(minvLcpiKp-mPDG)<invmasscut){
- if(selection>0){
-
-
- Double_t ptmax=0.;
-
- for (Int_t iprong=0; iprong<3; iprong++) {
- AliAODTrack *prong=(AliAODTrack*)part->GetDaughter(iprong);
- if (fReadMC) {
- labprong[iprong]=prong->GetLabel();
- if(pdgProngID[iprong]==2212) {
- fillthis="hbpIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- if(pdgProngID[iprong]==321) {
- fillthis="hbKIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- if(pdgProngID[iprong]==211) {
- fillthis="hbpiIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- }
- AliAODPid *pidObjtrk=prong->GetDetPid();
- AliAODPidHF *pidObj=new AliAODPidHF();
- Bool_t hasTOF=pidObj->CheckStatus(prong,"TOF");
- Bool_t hasTPC=pidObj->CheckStatus(prong,"TPC");
- delete pidObj;
- Double_t tofSignal=0.;
- Double_t dedxTPC=0.;
- Double_t momTOF=0.;
- Double_t momTPC=0.;
- if(hasTOF) {
- momTOF = prong->P();
- tofSignal=pidObjtrk->GetTOFsignal();
- }
- if(hasTPC) {
- momTPC = pidObjtrk->GetTPCmomentum();
- dedxTPC=pidObjtrk->GetTPCsignal();
- }
-
- switch (pdgProng[iprong]) {
- case 2212:
+ if (TMath::Abs(minvLcpKpi-mPDG)<invmasscut || TMath::Abs(minvLcpiKp-mPDG)<invmasscut) {
+
+ if (selection>0) { // 3prongs candidate x Lc (yes PID)
+
+ Double_t ptmax=0.;
+ for (Int_t iprong=0; iprong<3; iprong++) {
+ if(part->PtProng(iprong)>ptmax)ptmax=part->PtProng(iprong);
+
+ AliAODTrack *prong = (AliAODTrack*)part->GetDaughter(iprong);
+ AliAODPid *pidObjtrk = (AliAODPid*)prong->GetDetPid();
+ AliAODPidHF *pidObj = (AliAODPidHF*)cuts->GetPidHF();
+ Bool_t hasTOF=pidObj->CheckStatus(prong,"TOF");
+ Bool_t hasTPC=pidObj->CheckStatus(prong,"TPC");
+ Double_t tofSignal=0.;
+ Double_t dedxTPC=0.;
+ Double_t momTOF=0.;
+ Double_t momTPC=0.;
+ if(hasTOF) {
+ momTOF = prong->P();
+ tofSignal=pidObjtrk->GetTOFsignal();
+ }
+ if(hasTPC) {
+ momTPC = pidObjtrk->GetTPCmomentum();
+ dedxTPC=pidObjtrk->GetTPCsignal();
+ }
+
+ switch (pdgProngPID[iprong]) {
+ case 2212:
fillthis="hbpTOFSignal";
- ((TH1F*)listout->FindObject(fillthis))->Fill(tofSignal);
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(tofSignal);
fillthis="hbpTPCSignal";
- ((TH1F*)listout->FindObject(fillthis))->Fill(dedxTPC);
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(dedxTPC);
fillthis="hbpptProng";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
fillthis="hbpd0Prong";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->Getd0Prong(iprong));
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->Getd0Prong(iprong));
fillthis="hbpSignalVspTPC";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTPC,dedxTPC);
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTPC,dedxTPC);
fillthis="hbpSignalVspTOF";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTOF,tofSignal);
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTOF,tofSignal);
AliPID::EParticleType typep;
typep=AliPID::EParticleType(4);
if(hasTPC) {
- Double_t nsigmap = fUtilPid->NumberOfSigmasTPC(prong,typep);
- fillthis="hbpSigmaVspTPC";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTPC,nsigmap);
+ //Double_t nsigmap = fUtilPid->NumberOfSigmasTPC(prong,typep);
+ Double_t nsigmap = fPIDResponse->NumberOfSigmasTPC(prong,typep);
+ fillthis="hbpSigmaVspTPC";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTPC,nsigmap);
}
if(hasTOF){
- Double_t nsigma = fUtilPid->NumberOfSigmasTOF(prong,typep);
- fillthis="hbpSigmaVspTOF";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTOF,nsigma);
+ //Double_t nsigma = fUtilPid->NumberOfSigmasTOF(prong,typep);
+ Double_t nsigma=fPIDResponse->NumberOfSigmasTOF(prong,typep); // protection against 'old' AODs...
+ fillthis="hbpSigmaVspTOF";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTOF,nsigma);
}
- if(fReadMC){
- // real protons
- fillthis="hbpRealTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- // id protons
- if(pdgProngID[iprong]==2212) {
- fillthis="hbpIDGood";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- // misidentified kaons, pions
- if(pdgProngID[iprong]==321) {
- fillthis="hbKnonIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ if(fReadMC){ // fill hbpIDTot ---> PID contamination denominator
+ // hbIDGood, hbnoID ---> PID numerators
+ fillthis="hbpIDTot";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ if(pdgProngMC[iprong]==2212) { // id protons
+ fillthis="hbpIDGood";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
}
- if(pdgProngID[iprong]==211) {
- fillthis="hbpinonIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ else { // misidentified electrons, muons, pions and kaons
+ fillthis="hbnopIDp";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
}
}
break;
- case 321:
+
+ case 321:
fillthis="hbKTOFSignal";
- ((TH1F*)listout->FindObject(fillthis))->Fill(tofSignal);
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(tofSignal);
fillthis="hbKTPCSignal";
- ((TH1F*)listout->FindObject(fillthis))->Fill(dedxTPC);
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(dedxTPC);
fillthis="hbKptProng";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
fillthis="hbKd0Prong";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->Getd0Prong(iprong));
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->Getd0Prong(iprong));
fillthis="hbKSignalVspTPC";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTPC,dedxTPC);
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTPC,dedxTPC);
fillthis="hbKSignalVspTOF";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTOF,tofSignal);
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTOF,tofSignal);
AliPID::EParticleType typek;
typek=AliPID::EParticleType(3);
if(hasTPC) {
- Double_t nsigmap = fUtilPid->NumberOfSigmasTPC(prong,typek);
- fillthis="hbKSigmaVspTPC";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTPC,nsigmap);
+ //Double_t nsigmap = fUtilPid->NumberOfSigmasTPC(prong,typek);
+ Double_t nsigmap = fPIDResponse->NumberOfSigmasTPC(prong,typek);
+ fillthis="hbKSigmaVspTPC";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTPC,nsigmap);
}
if(hasTOF){
- Double_t nsigma = fUtilPid->NumberOfSigmasTOF(prong,typek);
- fillthis="hbKSigmaVspTOF";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTOF,nsigma);
+ //Double_t nsigma = fUtilPid->NumberOfSigmasTOF(prong,typek);
+ Double_t nsigma=fPIDResponse->NumberOfSigmasTOF(prong,typek); // protection against 'old' AODs...
+ fillthis="hbKSigmaVspTOF";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTOF,nsigma);
}
- if(fReadMC){
- // real kaons
- fillthis="hbKRealTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- // id kaons
- if(pdgProngID[iprong]==321) {
- fillthis="hbKIDGood";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ if (fReadMC) { // fill hbKIDTot ---> PID contamination denominator
+ // hIDGood, hnoID ---> PID numerators
+ fillthis="hbKIDTot";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ if(pdgProngMC[iprong]==321) { // id kaons
+ fillthis="hbKIDGood";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
}
- // misidentified protons, pions
- if(pdgProngID[iprong]==2212) {
- fillthis="hbpnonIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- if(pdgProngID[iprong]==211) {
- fillthis="hbpinonIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ else { // misidentified electrons, muons, pions and protons
+ fillthis="hbnokIDk";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
}
}
break;
- case 211:
+
+ case 211:
fillthis="hbpiTOFSignal";
- ((TH1F*)listout->FindObject(fillthis))->Fill(tofSignal);
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(tofSignal);
fillthis="hbpiTPCSignal";
- ((TH1F*)listout->FindObject(fillthis))->Fill(dedxTPC);
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(dedxTPC);
fillthis="hbpiptProng";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
fillthis="hbpid0Prong";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->Getd0Prong(iprong));
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->Getd0Prong(iprong));
fillthis="hbpiSignalVspTPC";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTPC,dedxTPC);
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTPC,dedxTPC);
fillthis="hbpiSignalVspTOF";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTOF,tofSignal);
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTOF,tofSignal);
AliPID::EParticleType typepi;
typepi=AliPID::EParticleType(2);
if(hasTPC) {
- Double_t nsigmap = fUtilPid->NumberOfSigmasTPC(prong,typepi);
- fillthis="hbpiSigmaVspTPC";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTPC,nsigmap);
+ //Double_t nsigmap = fUtilPid->NumberOfSigmasTPC(prong,typepi);
+ Double_t nsigmap = fPIDResponse->NumberOfSigmasTPC(prong,typepi);
+ fillthis="hbpiSigmaVspTPC";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTPC,nsigmap);
}
if(hasTOF){
- Double_t nsigma = fUtilPid->NumberOfSigmasTOF(prong,typepi);
- fillthis="hbpiSigmaVspTOF";
- ((TH2F*)listout->FindObject(fillthis))->Fill(momTOF,nsigma);
+ //Double_t nsigma = fUtilPid->NumberOfSigmasTOF(prong,typepi);
+ Double_t nsigma=fPIDResponse->NumberOfSigmasTOF(prong,typepi); // protection against 'old' AODs...
+ fillthis="hbpiSigmaVspTOF";
+ ((TH2F*)fOutputMC->FindObject(fillthis))->Fill(momTOF,nsigma);
}
- if(fReadMC){
- // real pions
- fillthis="hbpiRealTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- // id pions
- if(pdgProngID[iprong]==211) {
- fillthis="hbpiIDGood";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
- }
- // misidentified protons, kaons
- if(pdgProngID[iprong]==2212) {
- fillthis="hbpnonIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ if (fReadMC) { // fill hbpiIDTot ---> PID contamination denominator
+ // hIDGood, hnonIDTot ---> PID numerators
+ fillthis="hbpiIDTot";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ if(pdgProngMC[iprong]==211) { // id pions
+ fillthis="hbpiIDGood";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
}
- if(pdgProngID[iprong]==321) {
- fillthis="hbKnonIDTot";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->PtProng(iprong));
+ else { // misidentified electrons, muons, kaons and protons
+ fillthis="hbnopiIDpi";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->PtProng(iprong));
}
}
break;
- default:
- break;
- }
- if(part->PtProng(iprong)>ptmax)ptmax=part->PtProng(iprong);
- } //end loop on prongs
+ default:
+ break;
+ }
+
+ } //end loop on prongs
+
//now histograms where we don't need to check identity
fillthis="hbLcpt";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->Pt());
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->Pt());
+
fillthis = "hbDist12toPrim";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->GetDist12toPrim());
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->GetDist23toPrim());
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->GetDist12toPrim());
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->GetDist23toPrim());
fillthis = "hbSigmaVert";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->GetSigmaVert());
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->GetSigmaVert());
fillthis = "hbDCAs";
Double_t dcas[3];
part->GetDCAs(dcas);
- for (Int_t idca=0;idca<3;idca++) ((TH1F*)listout->FindObject(fillthis))->Fill(dcas[idca]);
+ for (Int_t idca=0;idca<3;idca++)
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(dcas[idca]);
fillthis = "hbCosPointingAngle";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->CosPointingAngle());
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->CosPointingAngle());
fillthis = "hbDecayLength";
- ((TH1F*)listout->FindObject(fillthis))->Fill(part->DecayLength());
- Double_t sum2=part->Getd0Prong(0)*part->Getd0Prong(0)+part->Getd0Prong(1)*part->Getd0Prong(1)+part->Getd0Prong(2)*part->Getd0Prong(2);
- fillthis = "hbSum2";
- ((TH1F*)listout->FindObject(fillthis))->Fill(sum2);
- fillthis = "hbptmax";
- ((TH1F*)listout->FindObject(fillthis))->Fill(ptmax);
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(part->DecayLength());
+ Double_t sum2=part->Getd0Prong(0)*part->Getd0Prong(0)+
+ part->Getd0Prong(1)*part->Getd0Prong(1)+
+ part->Getd0Prong(2)*part->Getd0Prong(2);
+ fillthis = "hbSum2";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(sum2);
+ fillthis = "hbptmax";
+ ((TH1F*)fOutputMC->FindObject(fillthis))->Fill(ptmax);
- }
+ } // end if (selection)
- } //end mass cut
+ } // end mass cut
} // end background case
+
return;
+
+}
+
+//-------------------------
+void AliAnalysisTaskSELambdac::FillAPrioriConcentrations(AliAODRecoDecayHF3Prong *part,
+ AliRDHFCutsLctopKpi *cuts,
+ AliAODEvent* aod,
+ TClonesArray *arrMC)
+{
+
+ cuts->SetUsePID(kFALSE); //Annalisa
+ Int_t isSelected3ProngByLc=cuts->IsSelected(part,AliRDHFCuts::kCandidate,aod);
+
+ TString fillthis="";
+
+ if(isSelected3ProngByLc>0 && fReadMC) {
+
+ for (Int_t ii=0; ii<3; ii++) {
+ AliAODTrack *prongTest=(AliAODTrack*)part->GetDaughter(ii);
+ if (!prongTest) continue;
+ Int_t labprongTest = prongTest->GetLabel();
+ if(labprongTest<0) continue;
+ AliAODMCParticle *mcprongTest = (AliAODMCParticle*)arrMC->At(labprongTest);
+
+ switch (TMath::Abs(mcprongTest->GetPdgCode())) {
+ case 11:
+ fillthis="hElIn3Prong";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 13:
+ fillthis="hMuIn3Prong";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 211:
+ fillthis="hPiIn3Prong";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 321:
+ fillthis="hKaIn3Prong";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 2212:
+ fillthis="hPrIn3Prong";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ default:
+ break;
+ }
+ }
+ }
+
+ cuts->SetUsePID(kTRUE); //Annalisa
+
+}
+
+//-----------------------
+Bool_t AliAnalysisTaskSELambdac::Is3ProngFromPDG(AliAODRecoDecayHF3Prong *d,
+ TClonesArray *arrayMC,
+ Int_t pdgToBeCompared)
+{
+ //
+ // Returs kTRUE if at least one of tracks in the current 3prong
+ // cames from particle with pdgCode=pdgToBeCompared
+ //
+
+ Bool_t localFlag = kFALSE;
+
+ for (Int_t ii=0;ii<3;ii++) {
+ AliAODTrack *daugh = (AliAODTrack*)d->GetDaughter(ii);
+
+ Int_t lab = daugh->GetLabel();
+ if(lab<0) continue;
+ AliAODMCParticle *part= (AliAODMCParticle*)arrayMC->At(lab);
+ //Int_t partPdg = part->GetPdgCode();
+ //printf(" -------------------------------------- partLab=%d partPdg=%d ---\n", lab, partPdg);
+
+ Int_t motherLabel = part->GetMother();
+ if(motherLabel<0) continue;
+
+ AliAODMCParticle *motherPart = 0;
+ Int_t motherPdg = 0;
+ while (!localFlag && motherLabel>=0) {
+ motherPart = (AliAODMCParticle*)arrayMC->At(motherLabel);
+ motherPdg = motherPart->GetPdgCode();
+ //printf(" ------------- motherLab=%d motherPdg=%d ---\n", motherLabel, motherPdg);
+ if (TMath::Abs(motherPdg)==pdgToBeCompared) {
+ //printf("-------------------------- trovato quark c/cbar\n");
+ localFlag = kTRUE;
+ }
+ else {
+ motherLabel = motherPart->GetMother();
+ }
+ }
+
+ }
+
+ return localFlag;
+
+}
+
+//-----------------------
+Bool_t AliAnalysisTaskSELambdac::IsTrackFromPDG(AliAODTrack *daugh,
+ TClonesArray *arrayMC,
+ Int_t pdgToBeCompared)
+{
+ //
+ // Returs kTRUE if at the tracks comes
+ // from particle with pdgCode=pdgToBeCompared
+ //
+
+ Bool_t localFlag = kFALSE;
+
+ Int_t lab = daugh->GetLabel();
+ if(lab<0) return localFlag;
+ AliAODMCParticle *part= (AliAODMCParticle*)arrayMC->At(lab);
+ //Int_t partPdg = part->GetPdgCode();
+ //printf(" -------------------------------------- partLab=%d partPdg=%d ---\n", lab, partPdg);
+
+ Int_t motherLabel = part->GetMother();
+ if(motherLabel<0) return localFlag;
+
+ AliAODMCParticle *motherPart = 0;
+ Int_t motherPdg = 0;
+ while (!localFlag && motherLabel>=0) {
+ motherPart = (AliAODMCParticle*)arrayMC->At(motherLabel);
+ motherPdg = motherPart->GetPdgCode();
+ //printf(" ------------- motherLab=%d motherPdg=%d ---\n", motherLabel, motherPdg);
+ if (TMath::Abs(motherPdg)==pdgToBeCompared) {
+ //printf("-------------------------- trovato quark c/cbar\n");
+ localFlag = kTRUE;
+ }
+ else {
+ motherLabel = motherPart->GetMother();
+ }
+ }
+
+ return localFlag;
+
}
+
+
+
+//-----------------------
+Bool_t AliAnalysisTaskSELambdac::IsThereAGeneratedLc(TClonesArray *arrayMC)
+{
+ //
+ // Returs kTRUE if there is a lepton related to the LambdaC
+ //
+
+ Bool_t localFlag = kFALSE;
+
+ AliAODMCParticle *searchLc;
+
+ for (Int_t iii=0; iii<arrayMC->GetEntries(); iii++) {
+ searchLc = (AliAODMCParticle*)arrayMC->At(iii);
+ Int_t searchLcpdg = searchLc->GetPdgCode();
+ if (TMath::Abs(searchLcpdg) == 4122){
+ localFlag = kTRUE;
+ break;
+ }
+ }
+
+
+ return localFlag;
+
+}
+//_________________________________
+///////////////estimate primaries
+Int_t AliAnalysisTaskSELambdac::NumberPrimaries(AliAODEvent *aods)
+{
+ Int_t counter=0;
+
+
+ TClonesArray *aodtracks=(TClonesArray *)aods->GetTracks();
+
+ // for(Int_t ji=0;ji<aods->GetNTracks();ji++)
+ for(Int_t ji=0;ji<aodtracks->GetEntriesFast();ji++)
+ {
+ AliAODTrack*aodTrack=(AliAODTrack*)aodtracks->UncheckedAt(ji);
+ if(aodTrack->IsPrimaryCandidate()) counter++;
+ }
+ return counter;
+
+}
+//_________________________________
+
+void AliAnalysisTaskSELambdac::MultiplicityStudies(AliAODRecoDecayHF3Prong *part,
+ AliRDHFCutsLctopKpi *cuts,
+ AliAODEvent* aod,
+ TClonesArray *arrMC,
+ Bool_t &flag1,Bool_t &flag2,Bool_t &flag3,
+ Bool_t &flag4, Bool_t &flag5, Bool_t &flag6)
+{
+
+ TString fillthis="";
+
+ Int_t pdgDgLctopKpi[3]={2212,321,211};
+ Int_t lab=-9999;
+ if(fReadMC)
+ lab=part->MatchToMC(4122,arrMC,3,pdgDgLctopKpi); //return MC particle label if the array corresponds to a Lc, -1 if not (cf. AliAODRecoDecay.cxx)
+
+ cuts->SetUsePID(kFALSE); //Annalisa
+ Int_t isSelected3ProngByLc=cuts->IsSelected(part,AliRDHFCuts::kCandidate,aod);
+
+ if(isSelected3ProngByLc>0 && fReadMC) {
+ flag1 = kTRUE;
+ if (lab>0) {
+ flag3 = kTRUE;
+ }
+
+ Bool_t is3ProngFromJPsi = Is3ProngFromPDG(part,arrMC,443);
+ if (is3ProngFromJPsi) flag6=is3ProngFromJPsi;
+
+ Bool_t is3ProngFromC = Is3ProngFromPDG(part,arrMC,4);
+ if (is3ProngFromC) flag4=is3ProngFromC;
+
+ Bool_t is3ProngFromB = Is3ProngFromPDG(part,arrMC,5);
+ if (is3ProngFromB) flag5=is3ProngFromB;
+
+ for (Int_t ii=0; ii<3; ii++) {
+ AliAODTrack *prongTest=(AliAODTrack*)part->GetDaughter(ii);
+ if (!prongTest) continue;
+ Int_t labprongTest = prongTest->GetLabel();
+ if(labprongTest<0) continue;
+ AliAODMCParticle *mcprongTest = (AliAODMCParticle*)arrMC->At(labprongTest);
+
+ switch (TMath::Abs(mcprongTest->GetPdgCode())) {
+ case 11:
+ fillthis="hElIn3Prong";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ if (IsTrackFromPDG(prongTest,arrMC,443)) {
+ fillthis="hElIn3Prong6";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ break;
+ case 13:
+ fillthis="hMuIn3Prong";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ if (IsTrackFromPDG(prongTest,arrMC,443)) {
+ fillthis="hMuIn3Prong6";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ break;
+ case 211:
+ fillthis="hPiIn3Prong";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ if (IsTrackFromPDG(prongTest,arrMC,443)) {
+ fillthis="hPiIn3Prong6";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ break;
+ case 321:
+ fillthis="hKaIn3Prong";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ if (IsTrackFromPDG(prongTest,arrMC,443)) {
+ fillthis="hKaIn3Prong6";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ break;
+ case 2212:
+ fillthis="hPrIn3Prong";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ if (IsTrackFromPDG(prongTest,arrMC,443)) {
+ fillthis="hPrIn3Prong6";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ break;
+ default:
+ break;
+ }
+
+ if (lab>0) {
+ switch (TMath::Abs(mcprongTest->GetPdgCode())) {
+ case 11:
+ fillthis="hElIn3Prong1";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 13:
+ fillthis="hMuIn3Prong1";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 211:
+ fillthis="hPiIn3Prong1";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 321:
+ fillthis="hKaIn3Prong1";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 2212:
+ fillthis="hPrIn3Prong1";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ default:
+ break;
+ }
+ } else {
+ switch (TMath::Abs(mcprongTest->GetPdgCode())) {
+ case 11:
+ fillthis="hElIn3Prong2";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ if (IsTrackFromPDG(prongTest,arrMC,4)) {
+ fillthis="hElIn3Prong3";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ if (IsTrackFromPDG(prongTest,arrMC,5)) {
+ fillthis="hElIn3Prong4";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ if (!IsTrackFromPDG(prongTest,arrMC,4) && !IsTrackFromPDG(prongTest,arrMC,5)) {
+ fillthis="hElIn3Prong5";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ break;
+ case 13:
+ fillthis="hMuIn3Prong2";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ if (IsTrackFromPDG(prongTest,arrMC,4)) {
+ fillthis="hMuIn3Prong3";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ if (IsTrackFromPDG(prongTest,arrMC,5)) {
+ fillthis="hMuIn3Prong4";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ if (!IsTrackFromPDG(prongTest,arrMC,4) && !IsTrackFromPDG(prongTest,arrMC,5)) {
+ fillthis="hMuIn3Prong5";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ break;
+ case 211:
+ fillthis="hPiIn3Prong2";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ if (IsTrackFromPDG(prongTest,arrMC,4)) {
+ fillthis="hPiIn3Prong3";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ if (IsTrackFromPDG(prongTest,arrMC,5)) {
+ fillthis="hPiIn3Prong4";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ if (!IsTrackFromPDG(prongTest,arrMC,4) && !IsTrackFromPDG(prongTest,arrMC,5)) {
+ fillthis="hPiIn3Prong5";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ break;
+ case 321:
+ fillthis="hKaIn3Prong2";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ if (IsTrackFromPDG(prongTest,arrMC,4)) {
+ fillthis="hKaIn3Prong3";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ if (IsTrackFromPDG(prongTest,arrMC,5)) {
+ fillthis="hKaIn3Prong4";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ if (!IsTrackFromPDG(prongTest,arrMC,4) && !IsTrackFromPDG(prongTest,arrMC,5)) {
+ fillthis="hKaIn3Prong5";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ break;
+ case 2212:
+ fillthis="hPrIn3Prong2";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ if (IsTrackFromPDG(prongTest,arrMC,4)) {
+ fillthis="hPrIn3Prong3";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ if (IsTrackFromPDG(prongTest,arrMC,5)) {
+ fillthis="hPrIn3Prong4";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ if (!IsTrackFromPDG(prongTest,arrMC,4) && !IsTrackFromPDG(prongTest,arrMC,5)) {
+ fillthis="hPrIn3Prong5";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ }
+ break;
+ default:
+ break;
+ }
+ }
+ /*
+ if (is3ProngFromC) {
+ switch (TMath::Abs(mcprongTest->GetPdgCode())) {
+ case 11:
+ fillthis="hElIn3Prong3";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 13:
+ fillthis="hMuIn3Prong3";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 211:
+ fillthis="hPiIn3Prong3";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 321:
+ fillthis="hKaIn3Prong3";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 2212:
+ fillthis="hPrIn3Prong3";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ default:
+ break;
+ }
+ } else { // !is3ProngFromC
+ if (is3ProngFromB) {
+ switch (TMath::Abs(mcprongTest->GetPdgCode())) {
+ case 11:
+ fillthis="hElIn3Prong4";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 13:
+ fillthis="hMuIn3Prong4";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 211:
+ fillthis="hPiIn3Prong4";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 321:
+ fillthis="hKaIn3Prong4";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 2212:
+ fillthis="hPrIn3Prong4";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ default:
+ break;
+ }
+ } else {//!is3ProngFromB
+ switch (TMath::Abs(mcprongTest->GetPdgCode())) {
+ case 11:
+ fillthis="hElIn3Prong5";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 13:
+ fillthis="hMuIn3Prong5";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 211:
+ fillthis="hPiIn3Prong5";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 321:
+ fillthis="hKaIn3Prong5";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ case 2212:
+ fillthis="hPrIn3Prong5";
+ ((TH1F*)fAPriori->FindObject(fillthis))->Fill(part->PtProng(ii));
+ break;
+ default:
+ break;
+ }
+ }
+ }
+ */
+ }
+ }
+
+ Double_t mPDG=TDatabasePDG::Instance()->GetParticle(4122)->Mass();
+ Double_t invmasscut=0.05;
+
+ Double_t minvLcpKpi = part->InvMassLcpKpi();
+ Double_t minvLcpiKp = part->InvMassLcpiKp();
+
+ cuts->SetUsePID(kTRUE); //Annalisa
+ Int_t isSelected3ProngByLcPID=cuts->IsSelected(part,AliRDHFCuts::kCandidate,aod);
+ if (isSelected3ProngByLcPID>0) {
+ if (TMath::Abs(minvLcpKpi-mPDG)<invmasscut || TMath::Abs(minvLcpiKp-mPDG)<invmasscut) {
+ flag2 = kTRUE;
+ }
+ }
+
+
+}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff