ClassImp(AliHFEmcQA)
//_______________________________________________________________________________________________
- AliHFEmcQA::AliHFEmcQA() :
- fMCEvent(NULL)
- ,fMCHeader(NULL)
- ,fMCArray(NULL)
- ,fQAhistos(NULL)
- ,fMCQACollection(NULL)
- ,fNparents(0)
- ,fCentrality(0)
- ,fIsPbPb(kFALSE)
- ,fIsppMultiBin(kFALSE)
- ,fContainerStep(0)
- ,fIsDebugStreamerON(kFALSE)
- ,fRecPt(-999)
- ,fRecEta(-999)
- ,fRecPhi(-999)
- ,fLyrhit(0)
- ,fLyrstat(0)
- ,fHfeImpactR(-999)
- ,fHfeImpactnsigmaR(-999)
- ,fTreeStream(NULL)
+AliHFEmcQA::AliHFEmcQA() :
+fMCEvent(NULL)
+ ,fMCHeader(NULL)
+ ,fMCArray(NULL)
+ ,fQAhistos(NULL)
+ ,fMCQACollection(NULL)
+ ,fNparents(0)
+ ,fCentrality(0)
+ ,fPerCentrality(-1)
+ ,fIsPbPb(kFALSE)
+ ,fIsppMultiBin(kFALSE)
+ ,fContainerStep(0)
+ ,fIsDebugStreamerON(kFALSE)
+ ,fRecPt(-999)
+ ,fRecEta(-999)
+ ,fRecPhi(-999)
+ ,fLyrhit(0)
+ ,fLyrstat(0)
+ ,fHfeImpactR(-999)
+ ,fHfeImpactnsigmaR(-999)
+ ,fTreeStream(NULL)
{
- // Default constructor
+ // Default constructor
for(Int_t mom = 0; mom < 9; mom++){
fhD[mom] = NULL;
}
//_______________________________________________________________________________________________
AliHFEmcQA::AliHFEmcQA(const AliHFEmcQA&p):
-TObject(p)
- ,fMCEvent(NULL)
- ,fMCHeader(NULL)
- ,fMCArray(NULL)
- ,fQAhistos(p.fQAhistos)
- ,fMCQACollection(p.fMCQACollection)
- ,fNparents(p.fNparents)
- ,fCentrality(0)
- ,fIsPbPb(kFALSE)
- ,fIsppMultiBin(kFALSE)
- ,fContainerStep(0)
- ,fIsDebugStreamerON(kFALSE)
- ,fRecPt(-999)
- ,fRecEta(-999)
- ,fRecPhi(-999)
- ,fLyrhit(0)
- ,fLyrstat(0)
- ,fHfeImpactR(0)
- ,fHfeImpactnsigmaR(0)
- ,fTreeStream(NULL)
+ TObject(p)
+ ,fMCEvent(NULL)
+ ,fMCHeader(NULL)
+ ,fMCArray(NULL)
+ ,fQAhistos(p.fQAhistos)
+ ,fMCQACollection(p.fMCQACollection)
+ ,fNparents(p.fNparents)
+ ,fCentrality(0)
+ ,fPerCentrality(-1)
+ ,fIsPbPb(kFALSE)
+ ,fIsppMultiBin(kFALSE)
+ ,fContainerStep(0)
+ ,fIsDebugStreamerON(kFALSE)
+ ,fRecPt(-999)
+ ,fRecEta(-999)
+ ,fRecPhi(-999)
+ ,fLyrhit(0)
+ ,fLyrstat(0)
+ ,fHfeImpactR(0)
+ ,fHfeImpactnsigmaR(0)
+ ,fTreeStream(NULL)
{
- // Copy constructor
+ // Copy constructor
for(Int_t mom = 0; mom < 9; mom++){
fhD[mom] = NULL;
}
memset(fElecBackgroundFactor, 0, sizeof(Double_t) * kElecBgSpecies * kBgPtBins * kCentBins * kBgLevels);
memset(fBinLimit, 0, sizeof(Double_t) * (kBgPtBins+1));
}
-
//_______________________________________________________________________________________________
AliHFEmcQA&
AliHFEmcQA::operator=(const AliHFEmcQA &)
{
// Assignment operator
-
+
AliInfo("Not yet implemented.");
return *this;
}
//_______________________________________________________________________________________________
AliHFEmcQA::~AliHFEmcQA()
{
- // Destructor
-
- if(fTreeStream && fIsDebugStreamerON) delete fTreeStream;
- AliInfo("Analysis Done.");
+ // Destructor
+
+ if(fTreeStream && fIsDebugStreamerON) delete fTreeStream;
+ AliInfo("Analysis Done.");
}
-
//_______________________________________________________________________________________________
void AliHFEmcQA::PostAnalyze() const
{
- //
- // Post analysis
- //
+ //
+ // Post analysis
+ //
}
-
//_______________________________________________________________________________________________
void AliHFEmcQA::SetBackgroundWeightFactor(Double_t *elecBackgroundFactor, Double_t *binLimit)
{
- //
- // copy background weighting factors into data member
- //
+ //
+ // copy background weighting factors into data member
+ //
memcpy(fElecBackgroundFactor,elecBackgroundFactor,sizeof(Double_t) * kElecBgSpecies * kBgPtBins * kCentBins * kBgLevels);
memcpy(fBinLimit,binLimit,sizeof(Double_t) * (kBgPtBins+1));
}
-
//__________________________________________
void AliHFEmcQA::CreatDefaultHistograms(TList * const qaList)
{
//
if(!qaList) return;
-
+
fQAhistos = qaList;
fQAhistos->SetName("MCqa");
-
+
CreateHistograms(AliHFEmcQA::kCharm); // create histograms for charm
CreateHistograms(AliHFEmcQA::kBeauty); // create histograms for beauty
CreateHistograms(AliHFEmcQA::kOthers); // create histograms for beauty
-
-// prepare 2D(pt vs Y) histogram for D spectra, we consider following 9 particles
+
+ // prepare 2D(pt vs Y) histogram for D spectra, we consider following 9 particles
const Int_t nbspecies = 9;
TString kDspecies[nbspecies];
kDspecies[0]="411"; //D+
fMCQACollection->CreateTH1Farray(Form("etapspectra_centrbin%i",centbin), "etap yields: MC p_{t} ", iBin[1],kPtRange);
fMCQACollection->CreateTH1Farray(Form("rhospectra_centrbin%i",centbin), "rho yields: MC p_{t} ", iBin[1],kPtRange);
+ fMCQACollection->CreateTH1Farray(Form("pionspectrapri_centrbin%i",centbin), "pion yields: MC p_{t} ", iBin[1],kPtRange);
+ fMCQACollection->CreateTH1Farray(Form("etaspectrapri_centrbin%i",centbin), "eta yields: MC p_{t} ", iBin[1],kPtRange);
+ fMCQACollection->CreateTH1Farray(Form("pionspectrasec_centrbin%i",centbin), "pion yields: MC p_{t} ", iBin[1],kPtRange);
+ fMCQACollection->CreateTH1Farray(Form("etaspectrasec_centrbin%i",centbin), "eta yields: MC p_{t} ", iBin[1],kPtRange);
+ fMCQACollection->CreateTH1Farray(Form("pionspectraLogpri_centrbin%i",centbin), "pion yields: MC p_{t} ", iBin[1],kPtRange);
+ fMCQACollection->CreateTH1Farray(Form("etaspectraLogpri_centrbin%i",centbin), "eta yields: MC p_{t} ", iBin[1],kPtRange);
+ fMCQACollection->CreateTH1Farray(Form("pionspectraLogsec_centrbin%i",centbin), "pion yields: MC p_{t} ", iBin[1],kPtRange);
+ fMCQACollection->CreateTH1Farray(Form("etaspectraLogsec_centrbin%i",centbin), "eta yields: MC p_{t} ", iBin[1],kPtRange);
+
fMCQACollection->CreateTH1F(Form("pionspectraLog_centrbin%i",centbin), "pion yields: MC p_{t} ", iBin[0],kPtbound[0], kPtbound[1], 1);
fMCQACollection->CreateTH1F(Form("etaspectraLog_centrbin%i",centbin), "eta yields: MC p_{t} ", iBin[0],kPtbound[0], kPtbound[1], 1);
fMCQACollection->CreateTH1F(Form("omegaspectraLog_centrbin%i",centbin), "omega yields: MC p_{t} ", iBin[0],kPtbound[0], kPtbound[1], 1);
xcorrbin[icorrbin]=icorrbin*0.1;
}
+ Int_t fCentrmax = 0;
+ if(!fIsPbPb&&!fIsppMultiBin) fCentrmax=0;
+ else fCentrmax = kCentBins;
TString hname;
if(kquark == kOthers){
for (Int_t icut = 0; icut < fgkEtaRanges; icut++ ){
- for (Int_t iqType = 0; iqType < 4; iqType++ ){
- hname = kqEtaRangeLabel[icut]+"Pt_"+kqTypeLabel[iqType];
- fHist[iq][iqType][icut].fPt = new TH1F(hname,hname+";p_{T} (GeV/c)",60,0.25,30.25);
- hname = kqEtaRangeLabel[icut]+"Y_"+kqTypeLabel[iqType];
- fHist[iq][iqType][icut].fY = new TH1F(hname,hname,150,-7.5,7.5);
- hname = kqEtaRangeLabel[icut]+"Eta_"+kqTypeLabel[iqType];
- fHist[iq][iqType][icut].fEta = new TH1F(hname,hname,150,-7.5,7.5);
- // Fill List
- if(fQAhistos) fHist[iq][iqType][icut].FillList(fQAhistos);
- }
+ for (Int_t iqType = 0; iqType < 4; iqType++ ){
+ for(Int_t icentr = 0; icentr < (fCentrmax+1); icentr++)
+ {
+ hname = kqEtaRangeLabel[icut]+"Pt_"+kqTypeLabel[iqType];
+ fHist[iq][iqType][icut][icentr].fPt = new TH1F(Form("%sCentr_%i",hname.Data(),icentr),Form("%sCentr_%i;p_{T} (GeV/c)",hname.Data(),icentr),60,0.25,30.25);
+ hname = kqEtaRangeLabel[icut]+"Y_"+kqTypeLabel[iqType];
+ fHist[iq][iqType][icut][icentr].fY = new TH1F(Form("%sCentr_%i",hname.Data(),icentr),Form("%sCentr_%i;y",hname.Data(),icentr),150,-7.5,7.5);
+ hname = kqEtaRangeLabel[icut]+"Eta_"+kqTypeLabel[iqType];
+ fHist[iq][iqType][icut][icentr].fEta = new TH1F(Form("%sCentr_%i",hname.Data(),icentr),Form("%sCentr_%i;eta",hname.Data(),icentr),150,-7.5,7.5);
+ // Fill List
+ if(fQAhistos) fHist[iq][iqType][icut][icentr].FillList(fQAhistos);
+ }
+ }
}
return;
}
for (Int_t icut = 0; icut < fgkEtaRanges; icut++ ){
for (Int_t iqType = 0; iqType < fgkqType; iqType++ ){
- if (iqType < keHadron && icut > 0) continue; // don't duplicate histogram for quark and hadron
- hname = kqEtaRangeLabel[icut]+"PdgCode_"+kqTypeLabel[iqType];
- fHist[iq][iqType][icut].fPdgCode = new TH1F(hname,hname,20001,-10000.5,10000.5);
- hname = kqEtaRangeLabel[icut]+"Pt_"+kqTypeLabel[iqType];
- fHist[iq][iqType][icut].fPt = new TH1F(hname,hname+";p_{T} (GeV/c)",iBin[1],kPtbinning1); // new binning
- hname = kqEtaRangeLabel[icut]+"Y_"+kqTypeLabel[iqType];
- fHist[iq][iqType][icut].fY = new TH1F(hname,hname,150,-7.5,7.5);
- hname = kqEtaRangeLabel[icut]+"Eta_"+kqTypeLabel[iqType];
- fHist[iq][iqType][icut].fEta = new TH1F(hname,hname,150,-7.5,7.5);
- // Fill List
- if(fQAhistos) fHist[iq][iqType][icut].FillList(fQAhistos);
+ if (iqType < keHadron && icut > 0) continue; // don't duplicate histogram for quark and hadron
+ for(Int_t icentr = 0; icentr<(fCentrmax+1); icentr++)
+ {
+ hname = kqEtaRangeLabel[icut]+"PdgCode_"+kqTypeLabel[iqType];
+ fHist[iq][iqType][icut][icentr].fPdgCode = new TH1F(Form("%sCentr_%i",hname.Data(),icentr),Form("%sCentr_%i;PdgCode",hname.Data(),icentr),20001,-10000.5,10000.5);
+ hname = kqEtaRangeLabel[icut]+"Pt_"+kqTypeLabel[iqType];
+ fHist[iq][iqType][icut][icentr].fPt = new TH1F(Form("%sCentr_%i",hname.Data(),icentr),Form("%sCentr_%i;p_{T} (GeV/c)",hname.Data(),icentr),iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"Y_"+kqTypeLabel[iqType];
+ fHist[iq][iqType][icut][icentr].fY = new TH1F(Form("%sCentr_%i",hname.Data(),icentr),Form("%sCentr_%i;y",hname.Data(),icentr),150,-7.5,7.5);
+ hname = kqEtaRangeLabel[icut]+"Eta_"+kqTypeLabel[iqType];
+ fHist[iq][iqType][icut][icentr].fEta = new TH1F(Form("%sCentr_%i",hname.Data(),icentr),Form("%sCentr_%i;eta",hname.Data(),icentr),150,-7.5,7.5);
+ // Fill List
+ if(fQAhistos) fHist[iq][iqType][icut][icentr].FillList(fQAhistos);
+ }
}
}
fHistComm[iq][icut].fPtCorrD0 = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1);
hname = kqEtaRangeLabel[icut]+"PtCorrDrest_"+kqTypeLabel[kQuark];
fHistComm[iq][icut].fPtCorrDrest = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1);
-
+
hname = kqEtaRangeLabel[icut]+"ePtRatio_"+kqTypeLabel[kQuark];
fHistComm[iq][icut].fePtRatio = new TH2F(hname,hname+";p_{T} (GeV/c);momentum fraction",200,0,20,100,0,1);
hname = kqEtaRangeLabel[icut]+"DePtRatio_"+kqTypeLabel[kQuark];
fHistComm[iq][icut].feDistance= new TH2F(hname,hname+";p_{T} (GeV/c);distance (cm)",100,0,20,200,0,2);
hname = kqEtaRangeLabel[icut]+"DeDistance_"+kqTypeLabel[kQuark];
fHistComm[iq][icut].fDeDistance= new TH2F(hname,hname+";p_{T} (GeV/c);distance (cm)",100,0,20,200,0,2);
+
+ if(icut <1){
+ hname = kqEtaRangeLabel[icut]+"PtCorrDinein_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrDinein = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrDineout_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrDineout = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrDoutein_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrDoutein = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrDouteout_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrDouteout = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+
+ hname = kqEtaRangeLabel[icut]+"PtCorrDpDinein_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrDpDinein = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrDpDineout_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrDpDineout = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrDpDoutein_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrDpDoutein = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrDpDouteout_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrDpDouteout = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+
+ hname = kqEtaRangeLabel[icut]+"PtCorrD0Dinein_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrD0Dinein = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrD0Dineout_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrD0Dineout = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrD0Doutein_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrD0Doutein = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrD0Douteout_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrD0Douteout = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+
+ hname = kqEtaRangeLabel[icut]+"PtCorrDrestDinein_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrDrestDinein = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrDrestDineout_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrDrestDineout = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrDrestDoutein_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrDrestDoutein = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrDrestDouteout_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrDrestDouteout = new TH2F(hname,hname+";D p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+
+ hname = kqEtaRangeLabel[icut]+"fEtaCorrD_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fEtaCorrD = new TH2F(hname,hname+";D Y;e eta",200,-10,10,200,-10,10);
+ hname = kqEtaRangeLabel[icut]+"fEtaCorrDp_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fEtaCorrDp = new TH2F(hname,hname+";D Y;e eta",200,-10,10,200,-10,10);
+ hname = kqEtaRangeLabel[icut]+"fEtaCorrD0_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fEtaCorrD0 = new TH2F(hname,hname+";D Y;e eta",200,-10,10,200,-10,10);
+ hname = kqEtaRangeLabel[icut]+"fEtaCorrDrest_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fEtaCorrDrest = new TH2F(hname,hname+";D Y;e eta",200,-10,10,200,-10,10);
+
+ hname = kqEtaRangeLabel[icut]+"fEtaCorrGD_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fEtaCorrGD = new TH2F(hname,hname+";D Y;e eta",200,-10,10,200,-10,10);
+ hname = kqEtaRangeLabel[icut]+"fEtaCorrGDp_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fEtaCorrGDp = new TH2F(hname,hname+";D Y;e eta",200,-10,10,200,-10,10);
+ hname = kqEtaRangeLabel[icut]+"fEtaCorrGD0_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fEtaCorrGD0 = new TH2F(hname,hname+";D Y;e eta",200,-10,10,200,-10,10);
+ hname = kqEtaRangeLabel[icut]+"fEtaCorrGDrest_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fEtaCorrGDrest = new TH2F(hname,hname+";D Y;e eta",200,-10,10,200,-10,10);
+
+ hname = kqEtaRangeLabel[icut]+"fEtaCorrB_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fEtaCorrB = new TH2F(hname,hname+";B Y;e eta",200,-10,10,200,-10,10);
+ hname = kqEtaRangeLabel[icut]+"fEtaCorrGB_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fEtaCorrGB = new TH2F(hname,hname+";B Y;e eta",200,-10,10,200,-10,10);
+
+ hname = kqEtaRangeLabel[icut]+"PtCorrBinein_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrBinein = new TH2F(hname,hname+";B p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrBineout_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrBineout = new TH2F(hname,hname+";B p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrBoutein_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrBoutein = new TH2F(hname,hname+";B p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ hname = kqEtaRangeLabel[icut]+"PtCorrBouteout_"+kqTypeLabel[kQuark];
+ fHistComm[iq][icut].fPtCorrBouteout = new TH2F(hname,hname+";B p_{T} (GeV/c);e p_{T} (GeV/c)",ndptbins,xcorrbin,iBin[1],kPtbinning1); // new binning
+ }
if(fQAhistos) fHistComm[iq][icut].FillList(fQAhistos);
}
+
hname = kqEtaRangeLabel[0]+"Nq_"+kqTypeLabel[kQuark];
fHistComm[iq][0].fNq = new TH1F(hname,hname,50,-0.5,49.5);
hname = kqEtaRangeLabel[0]+"ProcessID_"+kqTypeLabel[kQuark];
fParentSelect[1][5] = 5232; //Ksib0
fParentSelect[1][6] = 5332; //Omegab-
+
}
//__________________________________________
Int_t id1=0, id2=0;
- if(fCentrality>11) printf("warning centrality out of histogram array limits \n");
-
+ if(fCentrality>=11) {
+ AliWarning(Form("Centrality out of histogram array limits: %d", fCentrality));
+ return;
+ }
+ if(!fIsPbPb&&!fIsppMultiBin) fCentrality=0;
for(Int_t imc = 0; imc <fMCEvent->GetNumberOfPrimaries(); imc++){
if(!(mctrack2 = fMCEvent->GetTrack(imc))) continue;
mctrack0 = dynamic_cast<AliMCParticle *>(mctrack2);
if(!mctrack0) continue;
- if(!fIsPbPb&&!fIsppMultiBin) fCentrality=0;
+// if(!fIsPbPb&&!fIsppMultiBin) fCentrality=0;
- if(abs(mctrack0->PdgCode()) == 111) // pi0
+ if(TMath::Abs(mctrack0->PdgCode()) == 111) // pi0
{
if(TMath::Abs(AliHFEtools::GetRapidity(mcpart0))<0.8) {
fMCQACollection->Fill(Form("pionspectra_centrbin%i",fCentrality),mctrack0->Pt());
fMCQACollection->Fill(Form("pionspectraLog_centrbin%i",fCentrality),mctrack0->Pt());
+ if(imc>fMCEvent->GetNumberOfPrimaries()) {
+ fMCQACollection->Fill(Form("pionspectrasec_centrbin%i",fCentrality),mctrack0->Pt());
+ fMCQACollection->Fill(Form("pionspectraLogsec_centrbin%i",fCentrality),mctrack0->Pt());
+ }
+ else {
+ fMCQACollection->Fill(Form("pionspectrapri_centrbin%i",fCentrality),mctrack0->Pt());
+ fMCQACollection->Fill(Form("pionspectraLogpri_centrbin%i",fCentrality),mctrack0->Pt());
+ }
}
id1=mctrack0->GetFirstDaughter();
id2=mctrack0->GetLastDaughter();
for(int idx=id1; idx<=id2; idx++){
if(!(mctrackdaugt = fMCEvent->GetTrack(idx))) continue;
if(!(mctrackd = dynamic_cast<AliMCParticle *>(mctrackdaugt))) continue;
- if(abs(mctrackd->PdgCode()) == 11 && TMath::Abs(mctrackd->Eta())<0.8)
+ if(TMath::Abs(mctrackd->PdgCode()) == 11 && TMath::Abs(mctrackd->Eta())<0.8)
fMCQACollection->Fill(Form("piondaughters_centrbin%i",fCentrality),mctrackd->Pt());
}
}
- else if(abs(mctrack0->PdgCode()) == 221) // eta
+ else if(TMath::Abs(mctrack0->PdgCode()) == 221) // eta
{
if(TMath::Abs(AliHFEtools::GetRapidity(mcpart0))<0.8) {
fMCQACollection->Fill(Form("etaspectra_centrbin%i",fCentrality),mctrack0->Pt());
fMCQACollection->Fill(Form("etaspectraLog_centrbin%i",fCentrality),mctrack0->Pt());
+ if(imc>fMCEvent->GetNumberOfPrimaries()) {
+ fMCQACollection->Fill(Form("etaspectrasec_centrbin%i",fCentrality),mctrack0->Pt());
+ fMCQACollection->Fill(Form("etaspectraLogsec_centrbin%i",fCentrality),mctrack0->Pt());
+ }
+ else {
+ fMCQACollection->Fill(Form("etaspectrapri_centrbin%i",fCentrality),mctrack0->Pt());
+ fMCQACollection->Fill(Form("etaspectraLogpri_centrbin%i",fCentrality),mctrack0->Pt());
+ }
}
id1=mctrack0->GetFirstDaughter();
id2=mctrack0->GetLastDaughter();
for(int idx=id1; idx<=id2; idx++){
if(!(mctrackdaugt = fMCEvent->GetTrack(idx))) continue;
if(!(mctrackd = dynamic_cast<AliMCParticle *>(mctrackdaugt))) continue;
- if(abs(mctrackd->PdgCode()) == 11 && TMath::Abs(mctrackd->Eta())<0.8)
+ if(TMath::Abs(mctrackd->PdgCode()) == 11 && TMath::Abs(mctrackd->Eta())<0.8)
fMCQACollection->Fill(Form("etadaughters_centrbin%i",fCentrality),mctrackd->Pt());
}
}
- else if(abs(mctrack0->PdgCode()) == 223) // omega
+ else if(TMath::Abs(mctrack0->PdgCode()) == 223) // omega
{
if(TMath::Abs(AliHFEtools::GetRapidity(mcpart0))<0.8) {
fMCQACollection->Fill(Form("omegaspectra_centrbin%i",fCentrality),mctrack0->Pt());
for(int idx=id1; idx<=id2; idx++){
if(!(mctrackdaugt = fMCEvent->GetTrack(idx))) continue;
if(!(mctrackd = dynamic_cast<AliMCParticle *>(mctrackdaugt))) continue;
- if(abs(mctrackd->PdgCode()) == 11 && TMath::Abs(mctrackd->Eta())<0.8)
+ if(TMath::Abs(mctrackd->PdgCode()) == 11 && TMath::Abs(mctrackd->Eta())<0.8)
fMCQACollection->Fill(Form("omegadaughters_centrbin%i",fCentrality),mctrackd->Pt());
}
}
- else if(abs(mctrack0->PdgCode()) == 333) // phi
+ else if(TMath::Abs(mctrack0->PdgCode()) == 333) // phi
{
if(TMath::Abs(AliHFEtools::GetRapidity(mcpart0))<0.8) {
fMCQACollection->Fill(Form("phispectra_centrbin%i",fCentrality),mctrack0->Pt());
for(int idx=id1; idx<=id2; idx++){
if(!(mctrackdaugt = fMCEvent->GetTrack(idx))) continue;
if(!(mctrackd = dynamic_cast<AliMCParticle *>(mctrackdaugt))) continue;
- if(abs(mctrackd->PdgCode()) == 11 && TMath::Abs(mctrackd->Eta())<0.8)
+ if(TMath::Abs(mctrackd->PdgCode()) == 11 && TMath::Abs(mctrackd->Eta())<0.8)
fMCQACollection->Fill(Form("phidaughters_centrbin%i",fCentrality),mctrackd->Pt());
}
}
- else if(abs(mctrack0->PdgCode()) == 331) // eta prime
+ else if(TMath::Abs(mctrack0->PdgCode()) == 331) // eta prime
{
if(TMath::Abs(AliHFEtools::GetRapidity(mcpart0))<0.8) {
fMCQACollection->Fill(Form("etapspectra_centrbin%i",fCentrality),mctrack0->Pt());
for(int idx=id1; idx<=id2; idx++){
if(!(mctrackdaugt = fMCEvent->GetTrack(idx))) continue;
if(!(mctrackd = dynamic_cast<AliMCParticle *>(mctrackdaugt))) continue;
- if(abs(mctrackd->PdgCode()) == 11 && TMath::Abs(mctrackd->Eta())<0.8)
+ if(TMath::Abs(mctrackd->PdgCode()) == 11 && TMath::Abs(mctrackd->Eta())<0.8)
fMCQACollection->Fill(Form("etapdaughters_centrbin%i",fCentrality),mctrackd->Pt());
}
}
- else if(abs(mctrack0->PdgCode()) == 113) // rho
+ else if(TMath::Abs(mctrack0->PdgCode()) == 113) // rho
{
if(TMath::Abs(AliHFEtools::GetRapidity(mcpart0))<0.8) {
fMCQACollection->Fill(Form("rhospectra_centrbin%i",fCentrality),mctrack0->Pt());
for(int idx=id1; idx<=id2; idx++){
if(!(mctrackdaugt = fMCEvent->GetTrack(idx))) continue;
if(!(mctrackd = dynamic_cast<AliMCParticle *>(mctrackdaugt))) continue;
- if(abs(mctrackd->PdgCode()) == 11 && TMath::Abs(mctrackd->Eta())<0.8)
+ if(TMath::Abs(mctrackd->PdgCode()) == 11 && TMath::Abs(mctrackd->Eta())<0.8)
fMCQACollection->Fill(Form("rhodaughters_centrbin%i",fCentrality),mctrackd->Pt());
}
}
- else if(abs(mctrack0->PdgCode()) == 321) // kaon+-
+ else if(TMath::Abs(mctrack0->PdgCode()) == 321) // kaon+-
{
if(TMath::Abs(AliHFEtools::GetRapidity(mcpart0))<0.8) {
fMCQACollection->Fill(Form("kaonspectraLog_centrbin%i",fCentrality),mctrack0->Pt());
}
}
- else if(abs(mctrack0->PdgCode()) == 130) // k0L
+ else if(TMath::Abs(mctrack0->PdgCode()) == 130) // k0L
{
if(TMath::Abs(AliHFEtools::GetRapidity(mcpart0))<0.8) {
fMCQACollection->Fill(Form("k0LspectraLog_centrbin%i",fCentrality),mctrack0->Pt());
return;
}
+ if(!fIsPbPb&&!fIsppMultiBin) fCentrality=0;
+
AliMCParticle *mctrack = NULL;
Int_t partPdgcode = TMath::Abs(part->GetPdgCode());
// if the heavy particle comes from string which is denoted as particle status 12|12|12...12|11,[PYTHIA p.60]
// in this case, the mother of heavy particle can be one of the fragmented parton of the string
// should I make a condition that partMother should be quark or diquark? -> not necessary
- if ( abs(partMother->GetPdgCode()) == kquark || (partMother->GetStatusCode() == 11) ){
- //if ( abs(partMother->GetPdgCode()) == kquark || (partMother->GetStatusCode() == 11 || partMother->GetStatusCode() == 12) ){
+ if ( TMath::Abs(partMother->GetPdgCode()) == kquark || (partMother->GetStatusCode() == 11) ){
+ //if ( TMath::Abs(partMother->GetPdgCode()) == kquark || (partMother->GetStatusCode() == 11 || partMother->GetStatusCode() == 12) ){
- if ( abs(partMother->GetPdgCode()) != kquark ){
+ if ( TMath::Abs(partMother->GetPdgCode()) != kquark ){
// search fragmented partons in the same string
Bool_t isSameString = kTRUE;
for (Int_t i=1; i<fgkMaxIter; i++){
AliDebug(1, "Stack label is negative, return\n");
return;
}
- if ( abs(partMother->GetPdgCode()) == kquark ) break;
+ if ( TMath::Abs(partMother->GetPdgCode()) == kquark ) break;
if ( partMother->GetStatusCode() != 12 ) isSameString = kFALSE;
if (!isSameString) return;
}
}
AliDebug(1, "Can not find heavy parton of this heavy hadron in the string, return\n");
- if (abs(partMother->GetPdgCode()) != kquark) return;
+ if (TMath::Abs(partMother->GetPdgCode()) != kquark) return;
if (fIsHeavy[iq] >= 50) return;
fHeavyQuark[fIsHeavy[iq]] = partMother;
// fill kinematics for heavy parton
if (partMother->GetPdgCode() > 0) { // quark
- fHist[iq][kQuark][0].fPdgCode->Fill(partMother->GetPdgCode());
- fHist[iq][kQuark][0].fPt->Fill(partMother->Pt());
- fHist[iq][kQuark][0].fY->Fill(AliHFEtools::GetRapidity(partMother));
- fHist[iq][kQuark][0].fEta->Fill(partMother->Eta());
+ fHist[iq][kQuark][0][fCentrality].fPdgCode->Fill(partMother->GetPdgCode());
+ fHist[iq][kQuark][0][fCentrality].fPt->Fill(partMother->Pt());
+ fHist[iq][kQuark][0][fCentrality].fY->Fill(AliHFEtools::GetRapidity(partMother));
+ fHist[iq][kQuark][0][fCentrality].fEta->Fill(partMother->Eta());
} else{ // antiquark
- fHist[iq][kantiQuark][0].fPdgCode->Fill(partMother->GetPdgCode());
- fHist[iq][kantiQuark][0].fPt->Fill(partMother->Pt());
- fHist[iq][kantiQuark][0].fY->Fill(AliHFEtools::GetRapidity(partMother));
- fHist[iq][kantiQuark][0].fEta->Fill(partMother->Eta());
+ fHist[iq][kantiQuark][0][fCentrality].fPdgCode->Fill(partMother->GetPdgCode());
+ fHist[iq][kantiQuark][0][fCentrality].fPt->Fill(partMother->Pt());
+ fHist[iq][kantiQuark][0][fCentrality].fY->Fill(AliHFEtools::GetRapidity(partMother));
+ fHist[iq][kantiQuark][0][fCentrality].fEta->Fill(partMother->Eta());
}
} // end of heavy parton slection loop
return;
}
+ if(!fIsPbPb&&!fIsppMultiBin) fCentrality=0;
+
TParticle *partCopy = mcpart;
Int_t pdgcode = mcpart->GetPdgCode();
Int_t pdgcodeCopy = pdgcode;
// if the mother is charmed hadron
Bool_t isDirectCharm = kFALSE;
- if ( int(abs(pdgcode)/100.) == kCharm || int(abs(pdgcode)/1000.) == kCharm ) {
+ if ( int(TMath::Abs(pdgcode)/100.) == kCharm || int(TMath::Abs(pdgcode)/1000.) == kCharm ) {
// iterate until you find B hadron as a mother or become top ancester
for (Int_t i=1; i<fgkMaxIter; i++){
Int_t motherPDG = mother->GetPdgCode();
for (Int_t j=0; j<fNparents; j++){
- if (abs(motherPDG)==fParentSelect[1][j]) return; // return if this hadron is originated from b
+ if (TMath::Abs(motherPDG)==fParentSelect[1][j]) return; // return if this hadron is originated from b
}
mcpart = mother;
} // end of if
if((isDirectCharm == kTRUE && kquark == kCharm) || kquark == kBeauty) {
for (Int_t i=0; i<fNparents; i++){
- if (abs(pdgcodeCopy)==fParentSelect[iq][i]){
+ if (TMath::Abs(pdgcodeCopy)==fParentSelect[iq][i]){
// fill hadron kinematics
- fHist[iq][kHadron][0].fPdgCode->Fill(pdgcodeCopy);
- fHist[iq][kHadron][0].fPt->Fill(partCopy->Pt());
- fHist[iq][kHadron][0].fY->Fill(AliHFEtools::GetRapidity(partCopy));
- fHist[iq][kHadron][0].fEta->Fill(partCopy->Eta());
+ fHist[iq][kHadron][0][fCentrality].fPdgCode->Fill(pdgcodeCopy);
+ fHist[iq][kHadron][0][fCentrality].fPt->Fill(partCopy->Pt());
+ fHist[iq][kHadron][0][fCentrality].fY->Fill(AliHFEtools::GetRapidity(partCopy));
+ fHist[iq][kHadron][0][fCentrality].fEta->Fill(partCopy->Eta());
if(iq==0) {
fhD[i]->Fill(partCopy->Pt(),AliHFEtools::GetRapidity(partCopy));
}
}
// I also want to store D* info to compare with D* measurement
- if (abs(pdgcodeCopy)==413 && iq==0) { //D*+
+ if (TMath::Abs(pdgcodeCopy)==413 && iq==0) { //D*+
fhD[7]->Fill(partCopy->Pt(),AliHFEtools::GetRapidity(partCopy));
}
- if (abs(pdgcodeCopy)==423 && iq==0) { //D*0
+ if (TMath::Abs(pdgcodeCopy)==423 && iq==0) { //D*0
fhD[8]->Fill(partCopy->Pt(),AliHFEtools::GetRapidity(partCopy));
}
} // end of if
return;
}
+ if(!fIsPbPb&&!fIsppMultiBin) fCentrality=0;
+
Double_t eabsEta = TMath::Abs(mcpart->Eta());
Double_t eabsY = TMath::Abs(AliHFEtools::GetRapidity(mcpart));
if(kquark==kOthers){
Int_t esource = -1;
- if ( abs(mcpart->GetPdgCode()) != kdecayed ) esource = kMisID-4;
+ if ( TMath::Abs(mcpart->GetPdgCode()) != kdecayed ) esource = kMisID-4;
else esource =GetSource(mcpart)-4; // return for the cases kGamma=4, kPi0=5, kElse=6
if(esource==0|| esource==1 || esource==2 || esource==3){
- fHist[iq][esource][0].fPt->Fill(mcpart->Pt());
- fHist[iq][esource][0].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][esource][0].fEta->Fill(mcpart->Eta());
+ fHist[iq][esource][0][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][esource][0][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][esource][0][fCentrality].fEta->Fill(mcpart->Eta());
if(eabsEta<0.9){
- fHist[iq][esource][1].fPt->Fill(mcpart->Pt());
- fHist[iq][esource][1].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][esource][1].fEta->Fill(mcpart->Eta());
+ fHist[iq][esource][1][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][esource][1][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][esource][1][fCentrality].fEta->Fill(mcpart->Eta());
}
if(eabsY<0.5){
- fHist[iq][esource][2].fPt->Fill(mcpart->Pt());
- fHist[iq][esource][2].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][esource][2].fEta->Fill(mcpart->Eta());
+ fHist[iq][esource][2][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][esource][2][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][esource][2][fCentrality].fEta->Fill(mcpart->Eta());
}
return;
}
}
}
- if ( abs(mcpart->GetPdgCode()) != kdecayed ) return;
+ if ( TMath::Abs(mcpart->GetPdgCode()) != kdecayed ) return;
Int_t iLabel = mcpart->GetFirstMother();
if (iLabel<0){
// if the mother is charmed hadron
Bool_t isMotherDirectCharm = kFALSE;
- if ( int(abs(maPdgcode)/100.) == kCharm || int(abs(maPdgcode)/1000.) == kCharm ) {
+ if ( int(TMath::Abs(maPdgcode)/100.) == kCharm || int(TMath::Abs(maPdgcode)/1000.) == kCharm ) {
for (Int_t i=0; i<fNparents; i++){
- if (abs(maPdgcode)==fParentSelect[0][i]){
+ if (TMath::Abs(maPdgcode)==fParentSelect[0][i]){
isFinalOpenCharm = kTRUE;
}
}
Int_t grandMaPDG = grandMa->GetPdgCode();
for (Int_t j=0; j<fNparents; j++){
- if (abs(grandMaPDG)==fParentSelect[1][j]){
+ if (TMath::Abs(grandMaPDG)==fParentSelect[1][j]){
if (kquark == kCharm) return;
// fill electron kinematics
- fHist[iq][kElectron2nd][0].fPdgCode->Fill(mcpart->GetPdgCode());
- fHist[iq][kElectron2nd][0].fPt->Fill(mcpart->Pt());
- fHist[iq][kElectron2nd][0].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][kElectron2nd][0].fEta->Fill(mcpart->Eta());
+ fHist[iq][kElectron2nd][0][fCentrality].fPdgCode->Fill(mcpart->GetPdgCode());
+ fHist[iq][kElectron2nd][0][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][kElectron2nd][0][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][kElectron2nd][0][fCentrality].fEta->Fill(mcpart->Eta());
// fill mother hadron kinematics
- fHist[iq][kDeHadron][0].fPdgCode->Fill(grandMaPDG);
- fHist[iq][kDeHadron][0].fPt->Fill(grandMa->Pt());
- fHist[iq][kDeHadron][0].fY->Fill(AliHFEtools::GetRapidity(grandMa));
- fHist[iq][kDeHadron][0].fEta->Fill(grandMa->Eta());
+ fHist[iq][kDeHadron][0][fCentrality].fPdgCode->Fill(grandMaPDG);
+ fHist[iq][kDeHadron][0][fCentrality].fPt->Fill(grandMa->Pt());
+ fHist[iq][kDeHadron][0][fCentrality].fY->Fill(AliHFEtools::GetRapidity(grandMa));
+ fHist[iq][kDeHadron][0][fCentrality].fEta->Fill(grandMa->Eta());
if(eabsEta<0.9){
- fHist[iq][kElectron2nd][1].fPdgCode->Fill(mcpart->GetPdgCode());
- fHist[iq][kElectron2nd][1].fPt->Fill(mcpart->Pt());
- fHist[iq][kElectron2nd][1].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][kElectron2nd][1].fEta->Fill(mcpart->Eta());
+ fHist[iq][kElectron2nd][1][fCentrality].fPdgCode->Fill(mcpart->GetPdgCode());
+ fHist[iq][kElectron2nd][1][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][kElectron2nd][1][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][kElectron2nd][1][fCentrality].fEta->Fill(mcpart->Eta());
// fill mother hadron kinematics
- fHist[iq][kDeHadron][1].fPdgCode->Fill(grandMaPDG);
- fHist[iq][kDeHadron][1].fPt->Fill(grandMa->Pt());
- fHist[iq][kDeHadron][1].fY->Fill(AliHFEtools::GetRapidity(grandMa));
- fHist[iq][kDeHadron][1].fEta->Fill(grandMa->Eta());
+ fHist[iq][kDeHadron][1][fCentrality].fPdgCode->Fill(grandMaPDG);
+ fHist[iq][kDeHadron][1][fCentrality].fPt->Fill(grandMa->Pt());
+ fHist[iq][kDeHadron][1][fCentrality].fY->Fill(AliHFEtools::GetRapidity(grandMa));
+ fHist[iq][kDeHadron][1][fCentrality].fEta->Fill(grandMa->Eta());
}
if(eabsY<0.5){
- fHist[iq][kElectron2nd][2].fPdgCode->Fill(mcpart->GetPdgCode());
- fHist[iq][kElectron2nd][2].fPt->Fill(mcpart->Pt());
- fHist[iq][kElectron2nd][2].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][kElectron2nd][2].fEta->Fill(mcpart->Eta());
+ fHist[iq][kElectron2nd][2][fCentrality].fPdgCode->Fill(mcpart->GetPdgCode());
+ fHist[iq][kElectron2nd][2][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][kElectron2nd][2][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][kElectron2nd][2][fCentrality].fEta->Fill(mcpart->Eta());
// fill mother hadron kinematics
- fHist[iq][kDeHadron][2].fPdgCode->Fill(grandMaPDG);
- fHist[iq][kDeHadron][2].fPt->Fill(grandMa->Pt());
- fHist[iq][kDeHadron][2].fY->Fill(AliHFEtools::GetRapidity(grandMa));
- fHist[iq][kDeHadron][2].fEta->Fill(grandMa->Eta());
+ fHist[iq][kDeHadron][2][fCentrality].fPdgCode->Fill(grandMaPDG);
+ fHist[iq][kDeHadron][2][fCentrality].fPt->Fill(grandMa->Pt());
+ fHist[iq][kDeHadron][2][fCentrality].fY->Fill(AliHFEtools::GetRapidity(grandMa));
+ fHist[iq][kDeHadron][2][fCentrality].fEta->Fill(grandMa->Eta());
}
+ //mj: to calculate B to e eta correlation to calculate total heavy quark cross section
+ Int_t kLabel0 = grandMa->GetFirstMother();
+ Bool_t isGGrandmaYes = kFALSE;
+ Double_t ggmrapidwstmp=0;
+ if (!(kLabel0 < 0)){ // safety protection
+ if((mctrack = dynamic_cast<AliMCParticle *>(fMCEvent->GetTrack(TMath::Abs(kLabel0))))){
+ TParticle* ggrandMatmp = mctrack->Particle();
+ Int_t ggrandMaPDGtmp = ggrandMatmp->GetPdgCode();
+ if ( int(TMath::Abs(ggrandMaPDGtmp)/100.) == kBeauty || int(TMath::Abs(ggrandMaPDGtmp)/1000.) == kBeauty) isGGrandmaYes = kTRUE;
+ ggmrapidwstmp = AliHFEtools::GetRapidity(ggrandMatmp);
+ }
+ }
+
+ Double_t gmrapidwstmp0 = AliHFEtools::GetRapidity(grandMa);
+ Double_t eetawstmp0 = mcpart->Eta();
+
+ Double_t gmrapidtmp0 = TMath::Abs(gmrapidwstmp0);
+ Double_t eetatmp0 = TMath::Abs(eetawstmp0);
+
+ fHistComm[iq][0].fEtaCorrB->Fill(gmrapidwstmp0,eetawstmp0);
+ if(isGGrandmaYes) fHistComm[iq][0].fEtaCorrGB->Fill(ggmrapidwstmp,eetawstmp0);
+ else fHistComm[iq][0].fEtaCorrGB->Fill(gmrapidwstmp0,eetawstmp0);
+
+ if(gmrapidtmp0<0.5 && eetatmp0<0.5 ) fHistComm[iq][0].fPtCorrBinein->Fill(grandMa->Pt(),mcpart->Pt());
+ else if(gmrapidtmp0<0.5 && eetatmp0>0.5 ) fHistComm[iq][0].fPtCorrBineout->Fill(grandMa->Pt(),mcpart->Pt());
+ else if(gmrapidtmp0>0.5 && eetatmp0<0.5 ) fHistComm[iq][0].fPtCorrBoutein->Fill(grandMa->Pt(),mcpart->Pt());
+ else if(gmrapidtmp0>0.5 && eetatmp0>0.5 ) fHistComm[iq][0].fPtCorrBouteout->Fill(grandMa->Pt(),mcpart->Pt());
+ //======================================================================================
+
// ratio between pT of electron and pT of mother B hadron
if(grandMa->Pt()) {
fHistComm[iq][0].fDePtRatio->Fill(grandMa->Pt(),mcpart->Pt()/grandMa->Pt());
} // end of if
if((isMotherDirectCharm == kTRUE && kquark == kCharm) || kquark == kBeauty) {
for (Int_t i=0; i<fNparents; i++){
- if (abs(maPdgcodeCopy)==fParentSelect[iq][i]){
+ if (TMath::Abs(maPdgcodeCopy)==fParentSelect[iq][i]){
- fHist[iq][kElectron][0].fPdgCode->Fill(mcpart->GetPdgCode());
- fHist[iq][kElectron][0].fPt->Fill(mcpart->Pt());
- fHist[iq][kElectron][0].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][kElectron][0].fEta->Fill(mcpart->Eta());
+ fHist[iq][kElectron][0][fCentrality].fPdgCode->Fill(mcpart->GetPdgCode());
+ fHist[iq][kElectron][0][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][kElectron][0][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][kElectron][0][fCentrality].fEta->Fill(mcpart->Eta());
// fill mother hadron kinematics
- fHist[iq][keHadron][0].fPdgCode->Fill(maPdgcodeCopy);
- fHist[iq][keHadron][0].fPt->Fill(partMotherCopy->Pt());
- fHist[iq][keHadron][0].fY->Fill(AliHFEtools::GetRapidity(partMotherCopy));
- fHist[iq][keHadron][0].fEta->Fill(partMotherCopy->Eta());
+ fHist[iq][keHadron][0][fCentrality].fPdgCode->Fill(maPdgcodeCopy);
+ fHist[iq][keHadron][0][fCentrality].fPt->Fill(partMotherCopy->Pt());
+ fHist[iq][keHadron][0][fCentrality].fY->Fill(AliHFEtools::GetRapidity(partMotherCopy));
+ fHist[iq][keHadron][0][fCentrality].fEta->Fill(partMotherCopy->Eta());
if(eabsEta<0.9){
- fHist[iq][kElectron][1].fPdgCode->Fill(mcpart->GetPdgCode());
- fHist[iq][kElectron][1].fPt->Fill(mcpart->Pt());
- fHist[iq][kElectron][1].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][kElectron][1].fEta->Fill(mcpart->Eta());
+ fHist[iq][kElectron][1][fCentrality].fPdgCode->Fill(mcpart->GetPdgCode());
+ fHist[iq][kElectron][1][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][kElectron][1][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][kElectron][1][fCentrality].fEta->Fill(mcpart->Eta());
// fill mother hadron kinematics
- fHist[iq][keHadron][1].fPdgCode->Fill(maPdgcodeCopy);
- fHist[iq][keHadron][1].fPt->Fill(partMotherCopy->Pt());
- fHist[iq][keHadron][1].fY->Fill(AliHFEtools::GetRapidity(partMotherCopy));
- fHist[iq][keHadron][1].fEta->Fill(partMotherCopy->Eta());
+ fHist[iq][keHadron][1][fCentrality].fPdgCode->Fill(maPdgcodeCopy);
+ fHist[iq][keHadron][1][fCentrality].fPt->Fill(partMotherCopy->Pt());
+ fHist[iq][keHadron][1][fCentrality].fY->Fill(AliHFEtools::GetRapidity(partMotherCopy));
+ fHist[iq][keHadron][1][fCentrality].fEta->Fill(partMotherCopy->Eta());
}
if(eabsY<0.5){
- fHist[iq][kElectron][2].fPdgCode->Fill(mcpart->GetPdgCode());
- fHist[iq][kElectron][2].fPt->Fill(mcpart->Pt());
- fHist[iq][kElectron][2].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][kElectron][2].fEta->Fill(mcpart->Eta());
+ fHist[iq][kElectron][2][fCentrality].fPdgCode->Fill(mcpart->GetPdgCode());
+ fHist[iq][kElectron][2][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][kElectron][2][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][kElectron][2][fCentrality].fEta->Fill(mcpart->Eta());
// fill mother hadron kinematics
- fHist[iq][keHadron][2].fPdgCode->Fill(maPdgcodeCopy);
- fHist[iq][keHadron][2].fPt->Fill(partMotherCopy->Pt());
- fHist[iq][keHadron][2].fY->Fill(AliHFEtools::GetRapidity(partMotherCopy));
- fHist[iq][keHadron][2].fEta->Fill(partMotherCopy->Eta());
+ fHist[iq][keHadron][2][fCentrality].fPdgCode->Fill(maPdgcodeCopy);
+ fHist[iq][keHadron][2][fCentrality].fPt->Fill(partMotherCopy->Pt());
+ fHist[iq][keHadron][2][fCentrality].fY->Fill(AliHFEtools::GetRapidity(partMotherCopy));
+ fHist[iq][keHadron][2][fCentrality].fEta->Fill(partMotherCopy->Eta());
}
// ratio between pT of electron and pT of mother B or direct D hadron
}
}
+ //mj: to calculate D to e eta correlation to calculate total heavy quark cross section
+ Int_t kLabel = partMotherCopy->GetFirstMother();
+ Bool_t isGrandmaYes = kFALSE;
+ Double_t gmrapidwstmp =0;
+ if (!(kLabel < 0)){ // safety protection
+ if((mctrack = dynamic_cast<AliMCParticle *>(fMCEvent->GetTrack(TMath::Abs(kLabel))))){
+ TParticle* grandMatmp = mctrack->Particle();
+ Int_t grandMaPDGtmp = grandMatmp->GetPdgCode();
+ if ( int(TMath::Abs(grandMaPDGtmp)/100.) == kCharm || int(TMath::Abs(grandMaPDGtmp)/1000.) == kCharm ) isGrandmaYes = kTRUE;
+ if ( int(TMath::Abs(grandMaPDGtmp)/100.) == kBeauty || int(TMath::Abs(grandMaPDGtmp)/1000.) == kBeauty) isGrandmaYes = kTRUE;
+ gmrapidwstmp = AliHFEtools::GetRapidity(grandMatmp);
+ }
+ }
+
+ Double_t mrapidwstmp = AliHFEtools::GetRapidity(partMotherCopy);
+ Double_t eetawstmp = mcpart->Eta();
+
+ Double_t mrapidtmp = TMath::Abs(mrapidwstmp);
+ Double_t eetatmp = TMath::Abs(eetawstmp);
+
+ fHistComm[iq][0].fEtaCorrD->Fill(mrapidwstmp,eetawstmp);
+ if(isGrandmaYes) fHistComm[iq][0].fEtaCorrGD->Fill(gmrapidwstmp,eetawstmp);
+ else fHistComm[iq][0].fEtaCorrGD->Fill(mrapidwstmp,eetawstmp);
+
+ if(mrapidtmp<0.5 && eetatmp<0.5 ) fHistComm[iq][0].fPtCorrDinein->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ else if(mrapidtmp<0.5 && eetatmp>0.5 ) fHistComm[iq][0].fPtCorrDineout->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ else if(mrapidtmp>0.5 && eetatmp<0.5 ) fHistComm[iq][0].fPtCorrDoutein->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ else if(mrapidtmp>0.5 && eetatmp>0.5 ) fHistComm[iq][0].fPtCorrDouteout->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ if(TMath::Abs(partMotherCopy->GetPdgCode())==411) {
+ fHistComm[iq][0].fEtaCorrDp->Fill(mrapidwstmp,eetawstmp);
+ if(isGrandmaYes) fHistComm[iq][0].fEtaCorrGDp->Fill(gmrapidwstmp,eetawstmp);
+ else fHistComm[iq][0].fEtaCorrGDp->Fill(mrapidwstmp,eetawstmp);
+ if(mrapidtmp<0.5 && eetatmp<0.5 ) fHistComm[iq][0].fPtCorrDpDinein->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ else if(mrapidtmp<0.5 && eetatmp>0.5 ) fHistComm[iq][0].fPtCorrDpDineout->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ else if(mrapidtmp>0.5 && eetatmp<0.5 ) fHistComm[iq][0].fPtCorrDpDoutein->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ else if(mrapidtmp>0.5 && eetatmp>0.5 ) fHistComm[iq][0].fPtCorrDpDouteout->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ }
+ else if(TMath::Abs(partMotherCopy->GetPdgCode())==421) {
+ fHistComm[iq][0].fEtaCorrD0->Fill(mrapidwstmp,eetawstmp);
+ if(isGrandmaYes) fHistComm[iq][0].fEtaCorrGD0->Fill(gmrapidwstmp,eetawstmp);
+ else fHistComm[iq][0].fEtaCorrGD0->Fill(mrapidwstmp,eetawstmp);
+ if(mrapidtmp<0.5 && eetatmp<0.5 ) fHistComm[iq][0].fPtCorrD0Dinein->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ else if(mrapidtmp<0.5 && eetatmp>0.5 ) fHistComm[iq][0].fPtCorrD0Dineout->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ else if(mrapidtmp>0.5 && eetatmp<0.5 ) fHistComm[iq][0].fPtCorrD0Doutein->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ else if(mrapidtmp>0.5 && eetatmp>0.5 ) fHistComm[iq][0].fPtCorrD0Douteout->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ }
+ else {
+ fHistComm[iq][0].fEtaCorrDrest->Fill(mrapidwstmp,eetawstmp);
+ if(isGrandmaYes) fHistComm[iq][0].fEtaCorrGDrest->Fill(gmrapidwstmp,eetawstmp);
+ else fHistComm[iq][0].fEtaCorrGDrest->Fill(mrapidwstmp,eetawstmp);
+ if(mrapidtmp<0.5 && eetatmp<0.5 ) fHistComm[iq][0].fPtCorrDrestDinein->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ else if(mrapidtmp<0.5 && eetatmp>0.5 ) fHistComm[iq][0].fPtCorrDrestDineout->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ else if(mrapidtmp>0.5 && eetatmp<0.5 ) fHistComm[iq][0].fPtCorrDrestDoutein->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ else if(mrapidtmp>0.5 && eetatmp>0.5 ) fHistComm[iq][0].fPtCorrDrestDouteout->Fill(partMotherCopy->Pt(),mcpart->Pt());
+ }
+
// distance between electron production point and primary vertex
fHistComm[iq][0].feDistance->Fill(partMotherCopy->Pt(),decayLxy);
if(eabsEta<0.9){
return;
}
- if ( abs(mcpart->GetPdgCode()) != kdecayed ) return;
+ if ( TMath::Abs(mcpart->GetPdgCode()) != kdecayed ) return;
Double_t eabsEta = TMath::Abs(mcpart->Eta());
Double_t eabsY = TMath::Abs(AliHFEtools::GetRapidity(mcpart));
return;
}
+ if(!fIsPbPb&&!fIsppMultiBin) fCentrality=0;
+
AliAODMCParticle *partMother = (AliAODMCParticle*)fMCArray->At(iLabel);
AliAODMCParticle *partMotherCopy = partMother;
Int_t maPdgcode = partMother->GetPdgCode();
Int_t maPdgcodeCopy = maPdgcode;
Bool_t isMotherDirectCharm = kFALSE;
- if ( int(abs(maPdgcode)/100.) == kCharm || int(abs(maPdgcode)/1000.) == kCharm ) {
+ if ( int(TMath::Abs(maPdgcode)/100.) == kCharm || int(TMath::Abs(maPdgcode)/1000.) == kCharm ) {
for (Int_t i=0; i<fNparents; i++){
- if (abs(maPdgcode)==fParentSelect[0][i]){
+ if (TMath::Abs(maPdgcode)==fParentSelect[0][i]){
isFinalOpenCharm = kTRUE;
}
}
Int_t grandMaPDG = grandMa->GetPdgCode();
for (Int_t j=0; j<fNparents; j++){
- if (abs(grandMaPDG)==fParentSelect[1][j]){
+ if (TMath::Abs(grandMaPDG)==fParentSelect[1][j]){
if (kquark == kCharm) return;
// fill electron kinematics
- fHist[iq][kElectron2nd][0].fPdgCode->Fill(mcpart->GetPdgCode());
- fHist[iq][kElectron2nd][0].fPt->Fill(mcpart->Pt());
- fHist[iq][kElectron2nd][0].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][kElectron2nd][0].fEta->Fill(mcpart->Eta());
+ fHist[iq][kElectron2nd][0][fCentrality].fPdgCode->Fill(mcpart->GetPdgCode());
+ fHist[iq][kElectron2nd][0][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][kElectron2nd][0][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][kElectron2nd][0][fCentrality].fEta->Fill(mcpart->Eta());
// fill mother hadron kinematics
- fHist[iq][kDeHadron][0].fPdgCode->Fill(grandMaPDG);
- fHist[iq][kDeHadron][0].fPt->Fill(grandMa->Pt());
- fHist[iq][kDeHadron][0].fY->Fill(AliHFEtools::GetRapidity(grandMa));
- fHist[iq][kDeHadron][0].fEta->Fill(grandMa->Eta());
+ fHist[iq][kDeHadron][0][fCentrality].fPdgCode->Fill(grandMaPDG);
+ fHist[iq][kDeHadron][0][fCentrality].fPt->Fill(grandMa->Pt());
+ fHist[iq][kDeHadron][0][fCentrality].fY->Fill(AliHFEtools::GetRapidity(grandMa));
+ fHist[iq][kDeHadron][0][fCentrality].fEta->Fill(grandMa->Eta());
if(eabsEta<0.9){
// fill electron kinematics
- fHist[iq][kElectron2nd][1].fPdgCode->Fill(mcpart->GetPdgCode());
- fHist[iq][kElectron2nd][1].fPt->Fill(mcpart->Pt());
- fHist[iq][kElectron2nd][1].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][kElectron2nd][1].fEta->Fill(mcpart->Eta());
+ fHist[iq][kElectron2nd][1][fCentrality].fPdgCode->Fill(mcpart->GetPdgCode());
+ fHist[iq][kElectron2nd][1][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][kElectron2nd][1][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][kElectron2nd][1][fCentrality].fEta->Fill(mcpart->Eta());
// fill mother hadron kinematics
- fHist[iq][kDeHadron][1].fPdgCode->Fill(grandMaPDG);
- fHist[iq][kDeHadron][1].fPt->Fill(grandMa->Pt());
- fHist[iq][kDeHadron][1].fY->Fill(AliHFEtools::GetRapidity(grandMa));
- fHist[iq][kDeHadron][1].fEta->Fill(grandMa->Eta());
+ fHist[iq][kDeHadron][1][fCentrality].fPdgCode->Fill(grandMaPDG);
+ fHist[iq][kDeHadron][1][fCentrality].fPt->Fill(grandMa->Pt());
+ fHist[iq][kDeHadron][1][fCentrality].fY->Fill(AliHFEtools::GetRapidity(grandMa));
+ fHist[iq][kDeHadron][1][fCentrality].fEta->Fill(grandMa->Eta());
}
if(eabsY<0.5){
// fill electron kinematics
- fHist[iq][kElectron2nd][2].fPdgCode->Fill(mcpart->GetPdgCode());
- fHist[iq][kElectron2nd][2].fPt->Fill(mcpart->Pt());
- fHist[iq][kElectron2nd][2].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][kElectron2nd][2].fEta->Fill(mcpart->Eta());
+ fHist[iq][kElectron2nd][2][fCentrality].fPdgCode->Fill(mcpart->GetPdgCode());
+ fHist[iq][kElectron2nd][2][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][kElectron2nd][2][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][kElectron2nd][2][fCentrality].fEta->Fill(mcpart->Eta());
// fill mother hadron kinematics
- fHist[iq][kDeHadron][2].fPdgCode->Fill(grandMaPDG);
- fHist[iq][kDeHadron][2].fPt->Fill(grandMa->Pt());
- fHist[iq][kDeHadron][2].fY->Fill(AliHFEtools::GetRapidity(grandMa));
- fHist[iq][kDeHadron][2].fEta->Fill(grandMa->Eta());
+ fHist[iq][kDeHadron][2][fCentrality].fPdgCode->Fill(grandMaPDG);
+ fHist[iq][kDeHadron][2][fCentrality].fPt->Fill(grandMa->Pt());
+ fHist[iq][kDeHadron][2][fCentrality].fY->Fill(AliHFEtools::GetRapidity(grandMa));
+ fHist[iq][kDeHadron][2][fCentrality].fEta->Fill(grandMa->Eta());
}
return;
} // end of if
if ((isMotherDirectCharm == kTRUE && kquark == kCharm) || kquark == kBeauty) {
for (Int_t i=0; i<fNparents; i++){
- if (abs(maPdgcodeCopy)==fParentSelect[iq][i]){
+ if (TMath::Abs(maPdgcodeCopy)==fParentSelect[iq][i]){
// fill electron kinematics
- fHist[iq][kElectron][0].fPdgCode->Fill(mcpart->GetPdgCode());
- fHist[iq][kElectron][0].fPt->Fill(mcpart->Pt());
- fHist[iq][kElectron][0].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][kElectron][0].fEta->Fill(mcpart->Eta());
+ fHist[iq][kElectron][0][fCentrality].fPdgCode->Fill(mcpart->GetPdgCode());
+ fHist[iq][kElectron][0][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][kElectron][0][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][kElectron][0][fCentrality].fEta->Fill(mcpart->Eta());
// fill mother hadron kinematics
- fHist[iq][keHadron][0].fPdgCode->Fill(maPdgcodeCopy);
- fHist[iq][keHadron][0].fPt->Fill(partMotherCopy->Pt());
- fHist[iq][keHadron][0].fY->Fill(AliHFEtools::GetRapidity(partMotherCopy));
- fHist[iq][keHadron][0].fEta->Fill(partMotherCopy->Eta());
+ fHist[iq][keHadron][0][fCentrality].fPdgCode->Fill(maPdgcodeCopy);
+ fHist[iq][keHadron][0][fCentrality].fPt->Fill(partMotherCopy->Pt());
+ fHist[iq][keHadron][0][fCentrality].fY->Fill(AliHFEtools::GetRapidity(partMotherCopy));
+ fHist[iq][keHadron][0][fCentrality].fEta->Fill(partMotherCopy->Eta());
if(eabsEta<0.9){
// fill electron kinematics
- fHist[iq][kElectron][1].fPdgCode->Fill(mcpart->GetPdgCode());
- fHist[iq][kElectron][1].fPt->Fill(mcpart->Pt());
- fHist[iq][kElectron][1].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][kElectron][1].fEta->Fill(mcpart->Eta());
+ fHist[iq][kElectron][1][fCentrality].fPdgCode->Fill(mcpart->GetPdgCode());
+ fHist[iq][kElectron][1][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][kElectron][1][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][kElectron][1][fCentrality].fEta->Fill(mcpart->Eta());
// fill mother hadron kinematics
- fHist[iq][keHadron][1].fPdgCode->Fill(maPdgcodeCopy);
- fHist[iq][keHadron][1].fPt->Fill(partMotherCopy->Pt());
- fHist[iq][keHadron][1].fY->Fill(AliHFEtools::GetRapidity(partMotherCopy));
- fHist[iq][keHadron][1].fEta->Fill(partMotherCopy->Eta());
+ fHist[iq][keHadron][1][fCentrality].fPdgCode->Fill(maPdgcodeCopy);
+ fHist[iq][keHadron][1][fCentrality].fPt->Fill(partMotherCopy->Pt());
+ fHist[iq][keHadron][1][fCentrality].fY->Fill(AliHFEtools::GetRapidity(partMotherCopy));
+ fHist[iq][keHadron][1][fCentrality].fEta->Fill(partMotherCopy->Eta());
}
if(eabsY<0.5){
// fill electron kinematics
- fHist[iq][kElectron][2].fPdgCode->Fill(mcpart->GetPdgCode());
- fHist[iq][kElectron][2].fPt->Fill(mcpart->Pt());
- fHist[iq][kElectron][2].fY->Fill(AliHFEtools::GetRapidity(mcpart));
- fHist[iq][kElectron][2].fEta->Fill(mcpart->Eta());
+ fHist[iq][kElectron][2][fCentrality].fPdgCode->Fill(mcpart->GetPdgCode());
+ fHist[iq][kElectron][2][fCentrality].fPt->Fill(mcpart->Pt());
+ fHist[iq][kElectron][2][fCentrality].fY->Fill(AliHFEtools::GetRapidity(mcpart));
+ fHist[iq][kElectron][2][fCentrality].fEta->Fill(mcpart->Eta());
// fill mother hadron kinematics
- fHist[iq][keHadron][2].fPdgCode->Fill(maPdgcodeCopy);
- fHist[iq][keHadron][2].fPt->Fill(partMotherCopy->Pt());
- fHist[iq][keHadron][2].fY->Fill(AliHFEtools::GetRapidity(partMotherCopy));
- fHist[iq][keHadron][2].fEta->Fill(partMotherCopy->Eta());
+ fHist[iq][keHadron][2][fCentrality].fPdgCode->Fill(maPdgcodeCopy);
+ fHist[iq][keHadron][2][fCentrality].fPt->Fill(partMotherCopy->Pt());
+ fHist[iq][keHadron][2][fCentrality].fY->Fill(AliHFEtools::GetRapidity(partMotherCopy));
+ fHist[iq][keHadron][2][fCentrality].fEta->Fill(partMotherCopy->Eta());
}
}
motherlabel = -1;
- if (abs(afterinitialrad1->GetPdgCode()) == fgkGluon && abs(afterinitialrad2->GetPdgCode()) == fgkGluon){
+ if (TMath::Abs(afterinitialrad1->GetPdgCode()) == fgkGluon && TMath::Abs(afterinitialrad2->GetPdgCode()) == fgkGluon){
AliDebug(1,"heavy from gluon gluon pair creation!\n");
mothertype = -1;
motherID = fgkGluon;
}
- else if (abs(afterinitialrad1->GetPdgCode()) == kquark || abs(afterinitialrad2->GetPdgCode()) == kquark){ // one from Q and the other from g
+ else if (TMath::Abs(afterinitialrad1->GetPdgCode()) == kquark || TMath::Abs(afterinitialrad2->GetPdgCode()) == kquark){ // one from Q and the other from g
AliDebug(1,"heavy from flavor exitation!\n");
mothertype = -1;
motherID = kquark;
}
- else if (abs(afterinitialrad1->GetPdgCode()) == abs(afterinitialrad2->GetPdgCode())){
+ else if (TMath::Abs(afterinitialrad1->GetPdgCode()) == TMath::Abs(afterinitialrad2->GetPdgCode())){
AliDebug(1,"heavy from q-qbar pair creation!\n");
mothertype = -1;
motherID = 1;
}
//__________________________________________
-Int_t AliHFEmcQA::GetSource(const AliAODMCParticle * const mcpart)
+Int_t AliHFEmcQA::GetSource(const AliVParticle* const mcpart) const
{
// decay particle's origin
- //if ( abs(mcpart->GetPdgCode()) != AliHFEmcQA::kElectronPDG ) return -1;
+ //if ( TMath::Abs(mcpart->GetPdgCode()) != AliHFEmcQA::kElectronPDG ) return -1;
Int_t origin = -1;
Bool_t isFinalOpenCharm = kFALSE;
}
// mother
- Int_t iLabel = mcpart->GetMother();
+ // Information not in the base class, cast necessary
+ Int_t iLabel = GetMother(mcpart);
if (iLabel<0){
AliDebug(1, "Stack label is negative, return\n");
return -1;
}
- AliAODMCParticle *partMother = (AliAODMCParticle*)fMCArray->At(iLabel);
- Int_t maPdgcode = partMother->GetPdgCode();
+ const AliVParticle *partMother = fMCEvent->GetTrack(iLabel);
+ Int_t maPdgcode = partMother->PdgCode();
// if the mother is charmed hadron
- if ( int(abs(maPdgcode)/100.) == kCharm || int(abs(maPdgcode)/1000.) == kCharm ) {
+ if ( int(TMath::Abs(maPdgcode)/100.) == kCharm || int(TMath::Abs(maPdgcode)/1000.) == kCharm ) {
for (Int_t i=0; i<fNparents; i++){
- if (abs(maPdgcode)==fParentSelect[0][i]){
+ if (TMath::Abs(maPdgcode)==fParentSelect[0][i]){
isFinalOpenCharm = kTRUE;
}
}
// iterate until you find B hadron as a mother or become top ancester
for (Int_t i=1; i<fgkMaxIter; i++){
-
- Int_t jLabel = partMother->GetMother();
+
+ Int_t jLabel = GetMother(partMother);
if (jLabel == -1){
origin = kDirectCharm;
return origin;
}
// if there is an ancester
- AliAODMCParticle* grandMa = (AliAODMCParticle*)fMCArray->At(jLabel);
- Int_t grandMaPDG = grandMa->GetPdgCode();
+ const AliVParticle* grandMa = fMCEvent->GetTrack(jLabel);
+ Int_t grandMaPDG = grandMa->PdgCode();
for (Int_t j=0; j<fNparents; j++){
- if (abs(grandMaPDG)==fParentSelect[1][j]){
+ if (TMath::Abs(grandMaPDG)==fParentSelect[1][j]){
origin = kBeautyCharm;
return origin;
}
partMother = grandMa;
} // end of iteration
} // end of if
- else if ( int(abs(maPdgcode)/100.) == kBeauty || int(abs(maPdgcode)/1000.) == kBeauty ) {
+ else if ( int(TMath::Abs(maPdgcode)/100.) == kBeauty || int(TMath::Abs(maPdgcode)/1000.) == kBeauty ) {
for (Int_t i=0; i<fNparents; i++){
- if (abs(maPdgcode)==fParentSelect[1][i]){
+ if (TMath::Abs(maPdgcode)==fParentSelect[1][i]){
origin = kDirectBeauty;
return origin;
}
}
} // end of if
- else if ( abs(maPdgcode) == 22 ) {
+ else if ( TMath::Abs(maPdgcode) == 22 ) {
origin = kGamma;
return origin;
} // end of if
- else if ( abs(maPdgcode) == 111 ) {
+ else if ( TMath::Abs(maPdgcode) == 111 ) {
origin = kPi0;
return origin;
} // end of if
}
//__________________________________________
-Int_t AliHFEmcQA::GetSource(const TParticle * const mcpart)
+Int_t AliHFEmcQA::GetSource(const TParticle * const mcpart) const
{
// decay particle's origin
- //if ( abs(mcpart->GetPdgCode()) != AliHFEmcQA::kElectronPDG ) return -1;
+ //if ( TMath::Abs(mcpart->GetPdgCode()) != AliHFEmcQA::kElectronPDG ) return -1;
Int_t origin = -1;
Bool_t isFinalOpenCharm = kFALSE;
Int_t maPdgcode = partMother->GetPdgCode();
// if the mother is charmed hadron
- if ( int(abs(maPdgcode)/100.) == kCharm || int(abs(maPdgcode)/1000.) == kCharm ) {
+ if ( int(TMath::Abs(maPdgcode)/100.) == kCharm || int(TMath::Abs(maPdgcode)/1000.) == kCharm ) {
for (Int_t i=0; i<fNparents; i++){
- if (abs(maPdgcode)==fParentSelect[0][i]){
+ if (TMath::Abs(maPdgcode)==fParentSelect[0][i]){
isFinalOpenCharm = kTRUE;
}
}
Int_t grandMaPDG = grandMa->GetPdgCode();
for (Int_t j=0; j<fNparents; j++){
- if (abs(grandMaPDG)==fParentSelect[1][j]){
+ if (TMath::Abs(grandMaPDG)==fParentSelect[1][j]){
origin = kBeautyCharm;
return origin;
}
partMother = grandMa;
} // end of iteration
} // end of if
- else if ( int(abs(maPdgcode)/100.) == kBeauty || int(abs(maPdgcode)/1000.) == kBeauty ) {
+ else if ( int(TMath::Abs(maPdgcode)/100.) == kBeauty || int(TMath::Abs(maPdgcode)/1000.) == kBeauty ) {
for (Int_t i=0; i<fNparents; i++){
- if (abs(maPdgcode)==fParentSelect[1][i]){
+ if (TMath::Abs(maPdgcode)==fParentSelect[1][i]){
origin = kDirectBeauty;
return origin;
}
}
} // end of if
- else if ( abs(maPdgcode) == 22 ) {
+ else if ( TMath::Abs(maPdgcode) == 22 ) {
origin = kGamma;
return origin;
} // end of if
- else if ( abs(maPdgcode) == 111 ) {
+ else if ( TMath::Abs(maPdgcode) == 111 ) {
origin = kPi0;
return origin;
} // end of if
}
//__________________________________________
-Int_t AliHFEmcQA::GetElecSource(TParticle * const mcpart)
+Int_t AliHFEmcQA::GetElecSource(TParticle * const mcpart) const
{
// decay particle's origin
return -1;
}
- if ( abs(mcpart->GetPdgCode()) != AliHFEmcQA::kElectronPDG ) return kMisID;
+ if ( TMath::Abs(mcpart->GetPdgCode()) != AliHFEmcQA::kElectronPDG ) return kMisID;
Int_t origin = -1;
Bool_t isFinalOpenCharm = kFALSE;
TParticle *partMother = mctrack->Particle();
TParticle *partMotherCopy = mctrack->Particle();
Int_t maPdgcode = partMother->GetPdgCode();
+ Int_t grmaPdgcode;
+ Int_t ggrmaPdgcode;
// if the mother is charmed hadron
- if ( (int(abs(maPdgcode)/100.)%10) == kCharm || (int(abs(maPdgcode)/1000.)%10) == kCharm ) {
+
+ if(TMath::Abs(maPdgcode)==443){ // J/spi
+ Int_t jLabel = partMother->GetFirstMother();
+ if((mctrack = dynamic_cast<AliMCParticle *>(fMCEvent->GetTrack(TMath::Abs(jLabel))))){
+ TParticle* grandMa = mctrack->Particle();
+ Int_t grandMaPDG = grandMa->GetPdgCode();
+ if((int(TMath::Abs(grandMaPDG)/100.)%10) == kBeauty || (int(TMath::Abs(grandMaPDG)/1000.)%10) == kBeauty) return kB2Jpsi;
+ }
+ return kJpsi;
+ }
+ else if ( (int(TMath::Abs(maPdgcode)/100.)%10) == kCharm || (int(TMath::Abs(maPdgcode)/1000.)%10) == kCharm ) {
for (Int_t i=0; i<fNparents; i++){
- if (abs(maPdgcode)==fParentSelect[0][i]){
+ if (TMath::Abs(maPdgcode)==fParentSelect[0][i]){
isFinalOpenCharm = kTRUE;
}
}
Int_t jLabel = partMother->GetFirstMother();
if (jLabel == -1){
- origin = kDirectCharm;
- return origin;
+ return kDirectCharm;
}
if (jLabel < 0){ // safety protection
AliDebug(1, "Stack label is negative, return\n");
Int_t grandMaPDG = grandMa->GetPdgCode();
for (Int_t j=0; j<fNparents; j++){
- if (abs(grandMaPDG)==fParentSelect[1][j]){
- origin = kBeautyCharm;
- return origin;
+ if (TMath::Abs(grandMaPDG)==fParentSelect[1][j]){
+ return kBeautyCharm;
}
}
partMother = grandMa;
} // end of iteration
} // end of if
- else if ( (int(abs(maPdgcode)/100.)%10) == kBeauty || (int(abs(maPdgcode)/1000.)%10) == kBeauty ) {
+ else if ( (int(TMath::Abs(maPdgcode)/100.)%10) == kBeauty || (int(TMath::Abs(maPdgcode)/1000.)%10) == kBeauty ) {
for (Int_t i=0; i<fNparents; i++){
- if (abs(maPdgcode)==fParentSelect[1][i]){
- origin = kDirectBeauty;
- return origin;
+ if (TMath::Abs(maPdgcode)==fParentSelect[1][i]){
+ return kDirectBeauty;
}
}
} // end of if
- else if ( abs(maPdgcode) == 22 ) { //conversion
+ else if ( TMath::Abs(maPdgcode) == 22 ) { //conversion
tmpMomLabel = partMotherCopy->GetFirstMother();
if(!(mctrack = dynamic_cast<AliMCParticle *>(fMCEvent->GetTrack(TMath::Abs(tmpMomLabel))))) return -1;
partMother = mctrack->Particle();
maPdgcode = partMother->GetPdgCode();
- if ( abs(maPdgcode) == 111 ) {
- origin = kGammaPi0;
- return origin;
+
+ // check if the ligth meson is the decay product of heavy mesons
+ tmpMomLabel = partMother->GetFirstMother();
+ if((mctrack = dynamic_cast<AliMCParticle *>(fMCEvent->GetTrack(TMath::Abs(tmpMomLabel))))) {
+ partMother = mctrack->Particle();
+ grmaPdgcode = partMother->GetPdgCode();
+
+ if ( (int(TMath::Abs(grmaPdgcode)/100.)%10) == kBeauty || (int(TMath::Abs(grmaPdgcode)/1000.)%10) == kBeauty ) return kGammaB2M;
+ if ( (int(TMath::Abs(grmaPdgcode)/100.)%10) == kCharm || (int(TMath::Abs(grmaPdgcode)/1000.)%10) == kCharm ) return kGammaD2M;
+ if ( TMath::Abs(grmaPdgcode) == 221 ) return kGammaEta2Pi0; //mj to remove secondary pi0 decay electrons from eta decays, mainly to eta signal enhance samples
+ //if ( TMath::Abs(grmaPdgcode) == 221 ) return kElse; //mj to remove secondary pi0 decay electrons from eta decays, mainly to eta signal enhance samples
+
+ tmpMomLabel = partMother->GetFirstMother();
+ if((mctrack = dynamic_cast<AliMCParticle *>(fMCEvent->GetTrack(TMath::Abs(tmpMomLabel))))) {
+ partMother = mctrack->Particle();
+ ggrmaPdgcode = partMother->GetPdgCode();
+
+ if ( (int(TMath::Abs(ggrmaPdgcode)/100.)%10) == kBeauty || (int(TMath::Abs(ggrmaPdgcode)/1000.)%10) == kBeauty ) return kGammaB2M;
+ if ( (int(TMath::Abs(ggrmaPdgcode)/100.)%10) == kCharm || (int(TMath::Abs(ggrmaPdgcode)/1000.)%10) == kCharm ) return kGammaD2M;
+ }
+ }
+
+ if ( TMath::Abs(maPdgcode) == 111 ) {
+ return kGammaPi0;
}
- else if ( abs(maPdgcode) == 221 ) {
- origin = kGammaEta;
- return origin;
+ else if ( TMath::Abs(maPdgcode) == 221 ) {
+ return kGammaEta;
}
- else if ( abs(maPdgcode) == 223 ) {
- origin = kGammaOmega;
- return origin;
+ else if ( TMath::Abs(maPdgcode) == 223 ) {
+ return kGammaOmega;
}
- else if ( abs(maPdgcode) == 333 ) {
- origin = kGammaPhi;
- return origin;
+ else if ( TMath::Abs(maPdgcode) == 333 ) {
+ return kGammaPhi;
}
- else if ( abs(maPdgcode) == 331 ) {
- origin = kGammaEtaPrime;
- return origin;
+ else if ( TMath::Abs(maPdgcode) == 331 ) {
+ return kGammaEtaPrime;
}
- else if ( abs(maPdgcode) == 113 ) {
- origin = kGammaRho0;
- return origin;
+ else if ( TMath::Abs(maPdgcode) == 113 ) {
+ return kGammaRho0;
}
else origin = kElse;
- //origin = kGamma; // finer category above
return origin;
- } // end of if
- else if ( abs(maPdgcode) == 111 ) {
- origin = kPi0;
- return origin;
- } // end of if
- else if ( abs(maPdgcode) == 221 ) {
- origin = kEta;
- return origin;
- } // end of if
- else if ( abs(maPdgcode) == 223 ) {
- origin = kOmega;
- return origin;
- } // end of if
- else if ( abs(maPdgcode) == 333 ) {
- origin = kPhi;
- return origin;
- } // end of if
- else if ( abs(maPdgcode) == 331 ) {
- origin = kEtaPrime;
- return origin;
- } // end of if
- else if ( abs(maPdgcode) == 113 ) {
- origin = kRho0;
- return origin;
- } // end of if
- else{
- origin = kElse;
+ }
+ else {
+
+ // check if the ligth meson is the decay product of heavy mesons
+ tmpMomLabel = partMotherCopy->GetFirstMother();
+ if((mctrack = dynamic_cast<AliMCParticle *>(fMCEvent->GetTrack(TMath::Abs(tmpMomLabel))))) {
+ partMother = mctrack->Particle();
+ grmaPdgcode = partMother->GetPdgCode();
+
+ if ( (int(TMath::Abs(grmaPdgcode)/100.)%10) == kBeauty || (int(TMath::Abs(grmaPdgcode)/1000.)%10) == kBeauty ) return kB2M;
+ if ( (int(TMath::Abs(grmaPdgcode)/100.)%10) == kCharm || (int(TMath::Abs(grmaPdgcode)/1000.)%10) == kCharm ) return kD2M;
+ if ( TMath::Abs(grmaPdgcode) == 221 ) return kEta2Pi0; //mj to remove secondary pi0 decay electrons from eta decays, mainly to eta signal enhance samples
+
+ tmpMomLabel = partMother->GetFirstMother();
+ if((mctrack = dynamic_cast<AliMCParticle *>(fMCEvent->GetTrack(TMath::Abs(tmpMomLabel))))) {
+ partMother = mctrack->Particle();
+ ggrmaPdgcode = partMother->GetPdgCode();
+
+ if ( (int(TMath::Abs(ggrmaPdgcode)/100.)%10) == kBeauty || (int(TMath::Abs(ggrmaPdgcode)/1000.)%10) == kBeauty ) return kB2M;
+ if ( (int(TMath::Abs(ggrmaPdgcode)/100.)%10) == kCharm || (int(TMath::Abs(ggrmaPdgcode)/1000.)%10) == kCharm ) return kD2M;
+ }
+ }
+
+ if ( TMath::Abs(maPdgcode) == 111 ) {
+ return kPi0;
+ }
+ else if ( TMath::Abs(maPdgcode) == 221 ) {
+ return kEta;
+ }
+ else if ( TMath::Abs(maPdgcode) == 223 ) {
+ return kOmega;
+ }
+ else if ( TMath::Abs(maPdgcode) == 333 ) {
+ return kPhi;
+ }
+ else if ( TMath::Abs(maPdgcode) == 331 ) {
+ return kEtaPrime;
+ }
+ else if ( TMath::Abs(maPdgcode) == 113 ) {
+ return kRho0;
+ }
+ else if ( TMath::Abs(maPdgcode) == 321 || TMath::Abs(maPdgcode) == 130 ) {
+ return kKe3;
+ }
+ else{
+ origin = kElse;
+ }
+
}
return origin;
}
+
+//__________________________________________
+Int_t AliHFEmcQA::GetElecSource(const AliVParticle * const mctrack) const
+{
+ //
+ // decay particle's origin
+ //
+
+ if(!mctrack){
+ AliDebug(1, "no mcparticle, return\n");
+ return -1;
+ }
+
+ TClass *type = mctrack->IsA();
+
+ if(type == AliMCParticle::Class()) {
+ const AliMCParticle *esdmc = dynamic_cast<const AliMCParticle *>(mctrack);
+ if(esdmc){
+ TParticle *mcpart = esdmc->Particle();
+ return GetElecSource(mcpart);
+ }
+ else return -1;
+ }
+ if(type == AliAODMCParticle::Class()) {
+ const AliAODMCParticle *aodmc = dynamic_cast<const AliAODMCParticle *>(mctrack);
+ if(aodmc){
+ return GetElecSource(aodmc);
+ }
+ else return -1;
+ }
+ return -1;
+}
+//__________________________________________
+Int_t AliHFEmcQA::GetElecSource(const AliAODMCParticle * const mcpart) const
+{
+ //
+ // Function for AliAODMCParticle
+ //
+
+ if (!mcpart) return -1;
+ if (!fMCArray) return -1;
+
+ if ( TMath::Abs(mcpart->GetPdgCode()) != AliHFEmcQA::kElectronPDG ) return kMisID;
+
+ Int_t origin = -1;
+ Bool_t isFinalOpenCharm = kFALSE;
+
+ Int_t iLabel = mcpart->GetMother();
+ if ((iLabel<0) || (iLabel>=fMCArray->GetEntriesFast())){
+ AliDebug(1, "label is out of range, return\n");
+ return -1;
+ }
+
+ AliAODMCParticle *mctrack = NULL; // will change all the time
+ Int_t tmpMomLabel=0;
+ if(!(mctrack = dynamic_cast<AliAODMCParticle *>(fMCArray->At(TMath::Abs(iLabel))))) return -1;
+ AliAODMCParticle *partMother = mctrack;
+ AliAODMCParticle *partMotherCopy = mctrack;
+ Int_t maPdgcode = mctrack->GetPdgCode();
+ Int_t grmaPdgcode;
+ Int_t ggrmaPdgcode;
+
+ // if the mother is charmed hadron
+
+ if(TMath::Abs(maPdgcode)==443){
+ //
+ // J/spi
+ //
+ Int_t jLabel = partMother->GetMother();
+ if ((jLabel>=0) && (jLabel<fMCArray->GetEntriesFast())){
+ if((mctrack = dynamic_cast<AliAODMCParticle *>(fMCArray->At(TMath::Abs(jLabel))))){
+ Int_t grandMaPDG = mctrack->GetPdgCode();
+ if((int(TMath::Abs(grandMaPDG)/100.)%10) == kBeauty || (int(TMath::Abs(grandMaPDG)/1000.)%10) == kBeauty) {
+ return kB2Jpsi;
+ }
+ }
+ }
+ return kJpsi;
+ }
+ else if ( (int(TMath::Abs(maPdgcode)/100.)%10) == kCharm || (int(TMath::Abs(maPdgcode)/1000.)%10) == kCharm ) {
+ //
+ // charm
+ //
+ for (Int_t i=0; i<fNparents; i++){
+ if (TMath::Abs(maPdgcode)==fParentSelect[0][i]){
+ isFinalOpenCharm = kTRUE;
+ }
+ }
+ if (!isFinalOpenCharm) {
+ return -1;
+ }
+
+ // iterate until you find B hadron as a mother or become top ancester
+ for (Int_t i=1; i<fgkMaxIter; i++){
+
+ Int_t jLabel = partMother->GetMother();
+ if (jLabel == -1){
+ return kDirectCharm;
+ }
+ if ((jLabel<0) || (jLabel>=fMCArray->GetEntriesFast())){
+ AliDebug(1, "Stack label is negative, return\n");
+ return -1;
+ }
+
+ // if there is an ancester
+ if(!(mctrack = dynamic_cast<AliAODMCParticle *>(fMCArray->At(TMath::Abs(jLabel))))) {
+ return -1;
+ }
+ Int_t grandMaPDG = mctrack->GetPdgCode();
+ for (Int_t j=0; j<fNparents; j++){
+ if (TMath::Abs(grandMaPDG)==fParentSelect[1][j]){
+ return kBeautyCharm;
+ }
+ }
+ partMother = mctrack;
+ } // end of iteration
+
+ } // end of if
+ else if ( (int(TMath::Abs(maPdgcode)/100.)%10) == kBeauty || (int(TMath::Abs(maPdgcode)/1000.)%10) == kBeauty ) {
+ //
+ // beauty
+ //
+ for (Int_t i=0; i<fNparents; i++){
+ if (TMath::Abs(maPdgcode)==fParentSelect[1][i]){
+ return kDirectBeauty;
+ }
+ }
+ } // end of if
+ else if ( TMath::Abs(maPdgcode) == 22 ) {
+ //
+ //conversion
+ //
+ tmpMomLabel = partMotherCopy->GetMother();
+ if((tmpMomLabel<0) || (tmpMomLabel>=fMCArray->GetEntriesFast())) {
+ return -1;
+ }
+ if(!(mctrack = dynamic_cast<AliAODMCParticle *>(fMCArray->At(TMath::Abs(tmpMomLabel))))) {
+ return -1;
+ }
+ partMother = mctrack;
+ maPdgcode = partMother->GetPdgCode();
+
+ // check if the ligth meson is the decay product of heavy mesons
+ tmpMomLabel = partMother->GetMother();
+ if((tmpMomLabel>=0) && (tmpMomLabel<fMCArray->GetEntriesFast())) {
+ if((mctrack = dynamic_cast<AliAODMCParticle *>(fMCArray->At(TMath::Abs(tmpMomLabel))))) {
+ partMother = mctrack;
+ grmaPdgcode = partMother->GetPdgCode();
+
+ if ( (int(TMath::Abs(grmaPdgcode)/100.)%10) == kBeauty || (int(TMath::Abs(grmaPdgcode)/1000.)%10) == kBeauty ) {
+ return kGammaB2M;
+ }
+ if ( (int(TMath::Abs(grmaPdgcode)/100.)%10) == kCharm || (int(TMath::Abs(grmaPdgcode)/1000.)%10) == kCharm ) {
+ return kGammaD2M;
+ }
+ if ( TMath::Abs(grmaPdgcode) == 221 ) {
+ return kGammaEta2Pi0; //mj to remove secondary pi0 decay electrons from eta decays, mainly to eta signal enhance samples
+ //return kElse; //mj to remove secondary pi0 decay electrons from eta decays, mainly to eta signal enhance samples
+ }
+
+ tmpMomLabel = partMother->GetMother();
+ if((tmpMomLabel>=0) && (tmpMomLabel<fMCArray->GetEntriesFast())) {
+ if((mctrack = dynamic_cast<AliAODMCParticle *>(fMCArray->At(TMath::Abs(tmpMomLabel))))) {
+ partMother = mctrack;
+ ggrmaPdgcode = partMother->GetPdgCode();
+
+ if ( (int(TMath::Abs(ggrmaPdgcode)/100.)%10) == kBeauty || (int(TMath::Abs(ggrmaPdgcode)/1000.)%10) == kBeauty ) {
+ return kGammaB2M;
+ }
+ if ( (int(TMath::Abs(ggrmaPdgcode)/100.)%10) == kCharm || (int(TMath::Abs(ggrmaPdgcode)/1000.)%10) == kCharm ) {
+ return kGammaD2M;
+ }
+ }
+ }
+ }
+ }
+
+ if ( TMath::Abs(maPdgcode) == 111 ) {
+ return kGammaPi0;
+ }
+ else if ( TMath::Abs(maPdgcode) == 221 ) {
+ return kGammaEta;
+ }
+ else if ( TMath::Abs(maPdgcode) == 223 ) {
+ return kGammaOmega;
+ }
+ else if ( TMath::Abs(maPdgcode) == 333 ) {
+ return kGammaPhi;
+ }
+ else if ( TMath::Abs(maPdgcode) == 331 ) {
+ return kGammaEtaPrime;
+ }
+ else if ( TMath::Abs(maPdgcode) == 113 ) {
+ return kGammaRho0;
+ }
+ else origin = kElse;
+
+ return origin;
+
+ }
+ else {
+ //
+ // check if the ligth meson is the decay product of heavy mesons
+ //
+ tmpMomLabel = partMotherCopy->GetMother();
+ if((tmpMomLabel>=0) && (tmpMomLabel<fMCArray->GetEntriesFast())) {
+ if((mctrack = dynamic_cast<AliAODMCParticle *>(fMCArray->At(TMath::Abs(tmpMomLabel))))) {
+ partMother = mctrack;
+ grmaPdgcode = partMother->GetPdgCode();
+
+ if ( (int(TMath::Abs(grmaPdgcode)/100.)%10) == kBeauty || (int(TMath::Abs(grmaPdgcode)/1000.)%10) == kBeauty ) {
+ return kB2M;
+ }
+ if ( (int(TMath::Abs(grmaPdgcode)/100.)%10) == kCharm || (int(TMath::Abs(grmaPdgcode)/1000.)%10) == kCharm ) {
+ return kD2M;
+ }
+ if ( TMath::Abs(grmaPdgcode) == 221 ) {
+ return kEta2Pi0; //mj to remove secondary pi0 decay electrons from eta decays, mainly to eta signal enhance samples
+ }
+
+ tmpMomLabel = partMother->GetMother();
+ if((tmpMomLabel>=0) && (tmpMomLabel<fMCArray->GetEntriesFast())) {
+ if((mctrack = dynamic_cast<AliAODMCParticle *>(fMCArray->At(TMath::Abs(tmpMomLabel))))) {
+ partMother = mctrack;
+ ggrmaPdgcode = partMother->GetPdgCode();
+
+ if ( (int(TMath::Abs(ggrmaPdgcode)/100.)%10) == kBeauty || (int(TMath::Abs(ggrmaPdgcode)/1000.)%10) == kBeauty ) {
+ return kB2M;
+ }
+ if ( (int(TMath::Abs(ggrmaPdgcode)/100.)%10) == kCharm || (int(TMath::Abs(ggrmaPdgcode)/1000.)%10) == kCharm ) {
+ return kD2M;
+ }
+ }
+ }
+ }
+
+ if ( TMath::Abs(maPdgcode) == 111 ) {
+ return kPi0;
+ }
+ else if ( TMath::Abs(maPdgcode) == 221 ) {
+ return kEta;
+ }
+ else if ( TMath::Abs(maPdgcode) == 223 ) {
+ return kOmega;
+ }
+ else if ( TMath::Abs(maPdgcode) == 333 ) {
+ return kPhi;
+ }
+ else if ( TMath::Abs(maPdgcode) == 331 ) {
+ return kEtaPrime;
+ }
+ else if ( TMath::Abs(maPdgcode) == 113 ) {
+ return kRho0;
+ }
+ else if ( TMath::Abs(maPdgcode) == 321 || TMath::Abs(maPdgcode) == 130 ) {
+ return kKe3;
+ }
+ else{
+ origin = kElse;
+ }
+ }
+ }
+ return origin;
+}
//__________________________________________
Double_t AliHFEmcQA::GetWeightFactor(AliMCParticle *mctrack, const Int_t iBgLevel){
//
else if(mesonID==kGammaEtaPrime || mesonID==kEtaPrime) mArr=4; //etaprime
else if(mesonID==kGammaRho0 || mesonID==kRho0) mArr=5; //rho
- Double_t datamc[24]={-999,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999, -999, -999, -999, -999, -999, -999, -999, -999};
+ Double_t datamc[25]={-999,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999,-999, -999, -999, -999, -999, -999, -999, -999, -999, -999};
Double_t xr[3]={-999,-999,-999};
datamc[0] = mesonID;
datamc[17] = mctrack->Pt(); //electron pt
Int_t glabel=TMath::Abs(mctrack->GetMother());
if((mctrackmother = dynamic_cast<AliMCParticle *>(fMCEvent->GetTrack(glabel)))){
mesonPt=mctrackmother->Pt(); //meson pt
- bgcategory = -1.;
+ if(mesonID==kEta) bgcategory = -1.41; // to consider new branching ratio for the eta Dalitz decay
+ else bgcategory = -1.;
datamc[1] = bgcategory;
datamc[2] = mesonPt;
mctrackmother->XvYvZv(xr);
datamc[15] = mcpart->GetUniqueID();
}
if(glabel>fMCEvent->GetNumberOfPrimaries()) {
- bgcategory = -2.;
+ if(mesonID==kEta) bgcategory = -2.82; // to consider new branching ratio for the eta Dalitz decay (1.41)
+ else bgcategory = -2.;
datamc[1] = bgcategory;
glabel=TMath::Abs(mctrackmother->GetMother()); // gamma's mother's mother
if((mctrackmother = dynamic_cast<AliMCParticle *>(fMCEvent->GetTrack(glabel)))){
if(fIsPbPb){
- if(fCentrality < 0)return 0.;
- weightElecBg=fElecBackgroundFactor[iBgLevel][fCentrality][mArr][kBgPtBins-1];
+ Int_t centBin = GetWeightCentralityBin(fPerCentrality);
+ if(centBin < 0)return 0.;
+ weightElecBg=fElecBackgroundFactor[iBgLevel][centBin][mArr][kBgPtBins-1];
for(int ii=0; ii<kBgPtBins; ii++){
if((mesonPt > fBinLimit[ii]) && (mesonPt < fBinLimit[ii+1])){
- weightElecBg = fElecBackgroundFactor[iBgLevel][fCentrality][mArr][ii];
+ weightElecBg = fElecBackgroundFactor[iBgLevel][centBin][mArr][ii];
break;
}
}
datamc[21] = fRecPhi;
datamc[22] = fLyrhit;
datamc[23] = fLyrstat;
+ datamc[24] = fCentrality;
if(fIsDebugStreamerON && fTreeStream){
if(!iBgLevel){
"ereceta="<<datamc[20]<<
"erecphi="<<datamc[21]<<
"itshit="<<datamc[22]<<
- "itsstat="<<datamc[23]
+ "itsstat="<<datamc[23]<<
+ "centrality="<<datamc[24]
<< "\n";
}
}
return returnval;
}
+//__________________________________________
+Double_t AliHFEmcQA::GetWeightFactor(const AliAODMCParticle * const mcpart, const Int_t iBgLevel){
+ //
+ // Get weighting factor for the realistic background estimation, for three possible background yield levels, indicated by the argument "iLevel": the best estimate (0), the lower uncertainty level (1), and the upper uncertainty level (2)
+ //
+
+ if (!mcpart) return 0;
+ if (!fMCArray) return 0;
+
+ Double_t weightElecBg = 0.;
+ Double_t mesonPt = 0.;
+ Double_t bgcategory = 0.;
+ Int_t mArr = -1;
+ Int_t mesonID = GetElecSource(mcpart);
+ if(mesonID==kGammaPi0 || mesonID==kPi0) mArr=0; //pion
+ else if(mesonID==kGammaEta || mesonID==kEta) mArr=1; //eta
+ else if(mesonID==kGammaOmega || mesonID==kOmega) mArr=2; //omega
+ else if(mesonID==kGammaPhi || mesonID==kPhi) mArr=3; //phi
+ else if(mesonID==kGammaEtaPrime || mesonID==kEtaPrime) mArr=4; //etaprime
+ else if(mesonID==kGammaRho0 || mesonID==kRho0) mArr=5; //rho
+
+ if(!(mArr<0)){
+
+ AliAODMCParticle *mctrackmother = NULL; // will change all the time
+
+ if(mesonID>=kGammaPi0) { // conversion electron, be careful with the enum odering
+ Int_t iLabel = mcpart->GetMother(); //gamma label
+ if(!(mctrackmother = dynamic_cast<AliAODMCParticle *>(fMCArray->At(TMath::Abs(iLabel))))) return 0;
+ iLabel = mctrackmother->GetMother(); //gamma's mother's label
+ if(!(mctrackmother= dynamic_cast<AliAODMCParticle *>(fMCArray->At(TMath::Abs(iLabel))))) return 0;
+ mesonPt = mctrackmother->Pt(); //meson pt
+ bgcategory = 1.;
+
+ //if(glabel>fMCEvent->GetNumberOfPrimaries()) bgcategory = 2.;
+ }
+ else{ // nonHFE except for the conversion electron
+ Int_t iLabel = mcpart->GetMother(); //meson label
+ if(!(mctrackmother = dynamic_cast<AliAODMCParticle *>(fMCArray->At(TMath::Abs(iLabel))))) return 0;
+ mesonPt = mctrackmother->Pt(); //meson pt
+ if(mesonID==kEta) bgcategory = -1.41; // to consider new branching ratio for the eta Dalitz decay
+ else bgcategory = -1.;
+
+ //if(glabel>fMCEvent->GetNumberOfPrimaries()) bgcategory = -2.;
+ }
+
+ if(fIsPbPb){
+ Int_t centBin = GetWeightCentralityBin(fPerCentrality);
+ if(centBin < 0)return 0.;
+ weightElecBg=fElecBackgroundFactor[iBgLevel][centBin][mArr][kBgPtBins-1];
+ for(int ii=0; ii<kBgPtBins; ii++){
+ if((mesonPt > fBinLimit[ii]) && (mesonPt < fBinLimit[ii+1])){
+ weightElecBg = fElecBackgroundFactor[iBgLevel][centBin][mArr][ii];
+ break;
+ }
+ }
+ }
+ else{
+ weightElecBg=fElecBackgroundFactor[iBgLevel][0][mArr][kBgPtBins-1];
+ for(int ii=0; ii<kBgPtBins; ii++){
+ if((mesonPt > fBinLimit[ii]) && (mesonPt < fBinLimit[ii+1])){
+ weightElecBg = fElecBackgroundFactor[iBgLevel][0][mArr][ii];
+ break;
+ }
+ }
+ }
+ }
+
+ Double_t returnval = bgcategory*weightElecBg;
+ if(TMath::Abs(bgcategory)>1) returnval = bgcategory/2.;
+
+ return returnval;
+}
+
+//__________________________________________
+Int_t AliHFEmcQA::GetMother(const AliVParticle * const mcpart) const {
+ //
+ // Wrapper to get the mother label
+ //
+ Int_t label = -1;
+ TClass *type = mcpart->IsA();
+ if(type == AliMCParticle::Class()){
+ // ESD analysis
+ const AliMCParticle *emcpart = static_cast<const AliMCParticle *>(mcpart);
+ label = emcpart->GetMother();
+ } else if(type == AliAODMCParticle::Class()){
+ // AOD analysis
+ const AliAODMCParticle *amcpart = static_cast<const AliAODMCParticle *>(mcpart);
+ label = amcpart->GetMother();
+ }
+ return label;
+}
+//__________________________________________
+Int_t AliHFEmcQA::GetWeightCentralityBin(const Float_t percentile) const {
+ //
+ //translate the centrality percentile into the centrality bin of the reference weighting histograms for electron background
+ //
+
+ Float_t centralityLimits[12]= {0.,5.,10., 20., 30., 40., 50., 60.,70.,80., 90., 100.};
+ Int_t bin = -1;
+ for(Int_t ibin = 0; ibin < 11; ibin++){
+ if(percentile >= centralityLimits[ibin] && percentile < centralityLimits[ibin+1]){
+ bin = ibin;
+ break;
+ }
+ }
+ return bin;
+}
//__________________________________________
void AliHFEmcQA::AliHists::FillList(TList *l) const {
//
if(fPtCorrDp) l->Add(fPtCorrDp);
if(fPtCorrD0) l->Add(fPtCorrD0);
if(fPtCorrDrest) l->Add(fPtCorrDrest);
+
+ if(fPtCorrDinein) l->Add(fPtCorrDinein);
+ if(fPtCorrDineout) l->Add(fPtCorrDineout);
+ if(fPtCorrDoutein) l->Add(fPtCorrDoutein);
+ if(fPtCorrDouteout) l->Add(fPtCorrDouteout);
+ if(fPtCorrDpDinein) l->Add(fPtCorrDpDinein);
+ if(fPtCorrDpDineout) l->Add(fPtCorrDpDineout);
+ if(fPtCorrDpDoutein) l->Add(fPtCorrDpDoutein);
+ if(fPtCorrDpDouteout) l->Add(fPtCorrDpDouteout);
+ if(fPtCorrD0Dinein) l->Add(fPtCorrD0Dinein);
+ if(fPtCorrD0Dineout) l->Add(fPtCorrD0Dineout);
+ if(fPtCorrD0Doutein) l->Add(fPtCorrD0Doutein);
+ if(fPtCorrD0Douteout) l->Add(fPtCorrD0Douteout);
+ if(fPtCorrDrestDinein) l->Add(fPtCorrDrestDinein);
+ if(fPtCorrDrestDineout) l->Add(fPtCorrDrestDineout);
+ if(fPtCorrDrestDoutein) l->Add(fPtCorrDrestDoutein);
+ if(fPtCorrDrestDouteout) l->Add(fPtCorrDrestDouteout);
+
+ if(fEtaCorrD) l->Add(fEtaCorrD);
+ if(fEtaCorrDp) l->Add(fEtaCorrDp);
+ if(fEtaCorrD0) l->Add(fEtaCorrD0);
+ if(fEtaCorrDrest) l->Add(fEtaCorrDrest);
+
+ if(fEtaCorrGD) l->Add(fEtaCorrGD);
+ if(fEtaCorrGDp) l->Add(fEtaCorrGDp);
+ if(fEtaCorrGD0) l->Add(fEtaCorrGD0);
+ if(fEtaCorrGDrest) l->Add(fEtaCorrGDrest);
+
+ if(fEtaCorrB) l->Add(fEtaCorrB);
+ if(fEtaCorrGB) l->Add(fEtaCorrGB);
+ if(fPtCorrBinein) l->Add(fPtCorrBinein);
+ if(fPtCorrBineout) l->Add(fPtCorrBineout);
+ if(fPtCorrBoutein) l->Add(fPtCorrBoutein);
+ if(fPtCorrBouteout) l->Add(fPtCorrBouteout);
+
if(fDePtRatio) l->Add(fDePtRatio);
if(feDistance) l->Add(feDistance);
if(fDeDistance) l->Add(fDeDistance);
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff