* provided "as is" without express or implied warranty. *
**************************************************************************/
-/////////////////////////////////////////////////////////////
-//
-// AliAnalysisTaskSE for the extraction of signal(e.g D+) of heavy flavor
-// decay candidates with the MC truth.
+//*************************************************************************
+// Class AliAnalysisTaskSEDplus
+// AliAnalysisTaskSE for the D+ candidates Invariant Mass Histogram and
+//comparison of heavy-flavour decay candidates
+// to MC truth (kinematics stored in the AOD)
// Authors: Renu Bala, bala@to.infn.it
// F. Prino, prino@to.infn.it
// G. Ortona, ortona@to.infn.it
//________________________________________________________________________
AliAnalysisTaskSEDplus::AliAnalysisTaskSEDplus():
AliAnalysisTaskSE(),
-fOutput(0),
-fHistNEvents(0),
-fNtupleDplus(0),
-fUpmasslimit(1.965),
-fLowmasslimit(1.765),
-fNPtBins(0),
-fListCuts(0),
-fRDCutsProduction(0),
-fRDCutsAnalysis(0),
-fFillNtuple(kFALSE),
-fReadMC(kFALSE),
-fDoLS(kFALSE)
+ fOutput(0),
+ fHistNEvents(0),
+ fPtVsMass(0),
+ fPtVsMassTC(0),
+ fYVsPt(0),
+ fYVsPtTC(0),
+ fYVsPtSig(0),
+ fYVsPtSigTC(0),
+ fNtupleDplus(0),
+ fUpmasslimit(1.965),
+ fLowmasslimit(1.765),
+ fNPtBins(0),
+ fBinWidth(0.002),
+ fListCuts(0),
+ fRDCutsProduction(0),
+ fRDCutsAnalysis(0),
+ fFillNtuple(kFALSE),
+ fReadMC(kFALSE),
+ fDoLS(kFALSE)
{
// Default constructor
}
AliAnalysisTaskSE(name),
fOutput(0),
fHistNEvents(0),
+fPtVsMass(0),
+fPtVsMassTC(0),
+fYVsPt(0),
+fYVsPtTC(0),
+fYVsPtSig(0),
+fYVsPtSigTC(0),
fNtupleDplus(0),
fUpmasslimit(1.965),
fLowmasslimit(1.765),
fNPtBins(0),
+fBinWidth(0.002),
fListCuts(0),
fRDCutsProduction(dpluscutsprod),
fRDCutsAnalysis(dpluscutsana),
delete fOutput;
fOutput = 0;
}
-
+ if(fHistNEvents){
+ delete fHistNEvents;
+ fHistNEvents=0;
+ }
+ for(Int_t i=0;i<3*fNPtBins;i++){
+ if(fMassHist[i]){ delete fMassHist[i]; fMassHist[i]=0;}
+ if(fCosPHist[i]){ delete fCosPHist[i]; fCosPHist[i]=0;}
+ if(fDLenHist[i]){ delete fDLenHist[i]; fDLenHist[i]=0;}
+ if(fSumd02Hist[i]){ delete fSumd02Hist[i]; fSumd02Hist[i]=0;}
+ if(fSigVertHist[i]){ delete fSigVertHist[i]; fSigVertHist[i]=0;}
+ if(fPtMaxHist[i]){ delete fPtMaxHist[i]; fPtMaxHist[i]=0;}
+ if(fDCAHist[i]){ delete fDCAHist[i]; fDCAHist[i]=0;}
+ if(fMassHistTC[i]){ delete fMassHistTC[i]; fMassHistTC[i]=0;}
+ if(fMassHistTCPlus[i]){ delete fMassHistTCPlus[i]; fMassHistTCPlus[i]=0;}
+ if(fMassHistTCMinus[i]){ delete fMassHistTCMinus[i]; fMassHistTCMinus[i]=0;}
+
+ if(fMassHistLS[i]){ delete fMassHistLS[i]; fMassHistLS[i]=0;}
+ if(fCosPHistLS[i]){ delete fCosPHistLS[i]; fCosPHistLS[i]=0;}
+ if(fDLenHistLS[i]){ delete fDLenHistLS[i]; fDLenHistLS[i]=0;}
+ if(fSumd02HistLS[i]){ delete fSumd02HistLS[i]; fSumd02HistLS[i]=0;}
+ if(fSigVertHistLS[i]){ delete fSigVertHistLS[i]; fSigVertHistLS[i]=0;}
+ if(fPtMaxHistLS[i]){ delete fPtMaxHistLS[i]; fPtMaxHistLS[i]=0;}
+ if(fDCAHistLS[i]){ delete fDCAHistLS[i]; fDCAHistLS[i]=0;}
+ if(fMassHistLSTC[i]){ delete fMassHistLSTC[i]; fMassHistLSTC[i]=0;}
+ }
+ if(fPtVsMass){
+ delete fPtVsMass;
+ fPtVsMass=0;
+ }
+ if(fPtVsMassTC){
+ delete fPtVsMassTC;
+ fPtVsMassTC=0;
+ }
+ if(fYVsPt){
+ delete fYVsPt;
+ fYVsPt=0;
+ }
+ if(fYVsPtTC){
+ delete fYVsPtTC;
+ fYVsPtTC=0;
+ }
+ if(fYVsPtSig){
+ delete fYVsPtSig;
+ fYVsPtSig=0;
+ }
+ if(fYVsPtSigTC){
+ delete fYVsPtSigTC;
+ fYVsPtSigTC=0;
+ }
+ if(fNtupleDplus){
+ delete fNtupleDplus;
+ fNtupleDplus=0;
+ }
if (fListCuts) {
delete fListCuts;
fListCuts = 0;
}
-
if(fRDCutsProduction){
delete fRDCutsProduction;
fRDCutsProduction = 0;
}
-
- if(fRDCutsAnalysis){
+ if(fRDCutsAnalysis){
delete fRDCutsAnalysis;
fRDCutsAnalysis = 0;
}
+
+
}
//_________________________________________________________________
void AliAnalysisTaskSEDplus::SetMassLimits(Float_t range){
// set invariant mass limits
+ Float_t bw=GetBinWidth();
fUpmasslimit = 1.865+range;
fLowmasslimit = 1.865-range;
+ SetBinWidth(bw);
}
//_________________________________________________________________
void AliAnalysisTaskSEDplus::SetMassLimits(Float_t lowlimit, Float_t uplimit){
// set invariant mass limits
if(uplimit>lowlimit)
{
+ Float_t bw=GetBinWidth();
fUpmasslimit = lowlimit;
fLowmasslimit = uplimit;
+ SetBinWidth(bw);
}
}
//________________________________________________________________________
for(Int_t i=0; i<fNPtBins+1; i++) printf(" Bin%d = %8.2f-%8.2f\n",i,fArrayBinLimits[i],fArrayBinLimits[i+1]);
}
}
+//________________________________________________________________
+void AliAnalysisTaskSEDplus::SetBinWidth(Float_t w){
+ Float_t width=w;
+ Int_t nbins=(Int_t)((fUpmasslimit-fLowmasslimit)/width+0.5);
+ Int_t missingbins=4-nbins%4;
+ nbins=nbins+missingbins;
+ width=(fUpmasslimit-fLowmasslimit)/nbins;
+ if(missingbins!=0){
+ printf("AliAnalysisTaskSEDplus::SetBinWidth: W-bin width of %f will produce histograms not rebinnable by 4. New width set to %f\n",w,width);
+ }
+ else{
+ if(fDebug>1) printf("AliAnalysisTaskSEDplus::SetBinWidth: width set to %f\n",width);
+ }
+ fBinWidth=width;
+}
//_________________________________________________________________
Double_t AliAnalysisTaskSEDplus::GetPtBinLimit(Int_t ibin){
// get pt bin limit
if(ibin>fNPtBins)return -1;
return fArrayBinLimits[ibin];
}
-
+//_________________________________________________________________
+Int_t AliAnalysisTaskSEDplus::GetNBinsHistos(){
+ return (Int_t)((fUpmasslimit-fLowmasslimit)/fBinWidth+0.5);
+}
//_________________________________________________________________
void AliAnalysisTaskSEDplus::LSAnalysis(TClonesArray *arrayOppositeSign,TClonesArray *arrayLikeSign,AliAODEvent *aod,AliAODVertex *vtx1, Int_t nDplusOS){
//
Double_t cosp=d->CosPointingAngle();
Double_t sumD02=d->Getd0Prong(0)*d->Getd0Prong(0)+d->Getd0Prong(1)*d->Getd0Prong(1)+d->Getd0Prong(2)*d->Getd0Prong(2);
Double_t dca=d->GetDCA();
+ Double_t sigvert=d->GetSigmaVert();
Double_t ptmax=0;
for(Int_t i=0;i<3;i++){
if(d->PtProng(i)>ptmax)ptmax=d->PtProng(i);
fCosPHistLS[indexcut]->Fill(cosp);
fDLenHistLS[indexcut]->Fill(dlen);
fSumd02HistLS[indexcut]->Fill(sumD02);
+ fSigVertHistLS[indexcut]->Fill(sigvert);
fPtMaxHistLS[indexcut]->Fill(ptmax);
fDCAHistLS[indexcut]->Fill(dca);
- if(passTightCuts==1){
+ if(passTightCuts){
fMassHistLSTC[index]->Fill(invMass,wei);
fMassHistLSTC[index+1]->Fill(invMass);
fMassHistLSTC[index+2]->Fill(invMass,wei2);
TString hisname;
Int_t index=0;
Int_t indexLS=0;
+ Int_t nbins=GetNBinsHistos();
for(Int_t i=0;i<fNPtBins;i++){
index=GetHistoIndex(i);
indexLS=GetLSHistoIndex(i);
hisname.Form("hMassPt%d",i);
- fMassHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHist[index]->Sumw2();
hisname.Form("hCosPAllPt%d",i);
- fCosPHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.5,1.);
+ fCosPHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.5,1.);
fCosPHist[index]->Sumw2();
hisname.Form("hDLenAllPt%d",i);
- fDLenHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.5);
+ fDLenHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.5);
fDLenHist[index]->Sumw2();
hisname.Form("hSumd02AllPt%d",i);
- fSumd02Hist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,1.);
+ fSumd02Hist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,1.);
fSumd02Hist[index]->Sumw2();
hisname.Form("hSigVertAllPt%d",i);
- fSigVertHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.1);
+ fSigVertHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.1);
fSigVertHist[index]->Sumw2();
hisname.Form("hPtMaxAllPt%d",i);
- fPtMaxHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.5,5.);
+ fPtMaxHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.5,5.);
fPtMaxHist[index]->Sumw2();
hisname.Form("hDCAAllPt%d",i);
- fDCAHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.1);
+ fDCAHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.1);
fDCAHist[index]->Sumw2();
hisname.Form("hMassPt%dTC",i);
- fMassHistTC[index]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHistTC[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHistTC[index]->Sumw2();
-
+ hisname.Form("hMassPt%dTCPlus",i);
+ fMassHistTCPlus[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
+ fMassHistTCPlus[index]->Sumw2();
+ hisname.Form("hMassPt%dTCMinus",i);
+ fMassHistTCMinus[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
+ fMassHistTCMinus[index]->Sumw2();
hisname.Form("hCosPAllPt%dLS",i);
- fCosPHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.5,1.);
+ fCosPHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.5,1.);
fCosPHistLS[index]->Sumw2();
hisname.Form("hDLenAllPt%dLS",i);
- fDLenHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.5);
+ fDLenHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.5);
fDLenHistLS[index]->Sumw2();
hisname.Form("hSumd02AllPt%dLS",i);
- fSumd02HistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,1.);
+ fSumd02HistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,1.);
fSumd02HistLS[index]->Sumw2();
hisname.Form("hSigVertAllPt%dLS",i);
- fSigVertHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.1);
+ fSigVertHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.1);
fSigVertHistLS[index]->Sumw2();
hisname.Form("hPtMaxAllPt%dLS",i);
- fPtMaxHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.5,5.);
+ fPtMaxHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.5,5.);
fPtMaxHistLS[index]->Sumw2();
hisname.Form("hDCAAllPt%dLS",i);
- fDCAHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.1);
+ fDCAHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.1);
fDCAHistLS[index]->Sumw2();
hisname.Form("hLSPt%dLC",i);
- fMassHistLS[indexLS] = new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHistLS[indexLS] = new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHistLS[indexLS]->Sumw2();
hisname.Form("hLSPt%dTC",i);
- fMassHistLSTC[indexLS] = new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHistLSTC[indexLS] = new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHistLSTC[indexLS]->Sumw2();
index=GetSignalHistoIndex(i);
indexLS++;
hisname.Form("hSigPt%d",i);
- fMassHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHist[index]->Sumw2();
hisname.Form("hCosPSigPt%d",i);
- fCosPHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.5,1.);
+ fCosPHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.5,1.);
fCosPHist[index]->Sumw2();
hisname.Form("hDLenSigPt%d",i);
- fDLenHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.5);
+ fDLenHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.5);
fDLenHist[index]->Sumw2();
hisname.Form("hSumd02SigPt%d",i);
- fSumd02Hist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,1.);
+ fSumd02Hist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,1.);
fSumd02Hist[index]->Sumw2();
hisname.Form("hSigVertSigPt%d",i);
- fSigVertHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.1);
+ fSigVertHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.1);
fSigVertHist[index]->Sumw2();
hisname.Form("hPtMaxSigPt%d",i);
- fPtMaxHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.5,5.);
+ fPtMaxHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.5,5.);
fPtMaxHist[index]->Sumw2();
hisname.Form("hDCASigPt%d",i);
- fDCAHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.1);
+ fDCAHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.1);
fDCAHist[index]->Sumw2();
hisname.Form("hSigPt%dTC",i);
- fMassHistTC[index]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHistTC[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHistTC[index]->Sumw2();
+ hisname.Form("hSigPt%dTCPlus",i);
+ fMassHistTCPlus[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
+ fMassHistTCPlus[index]->Sumw2();
+ hisname.Form("hSigPt%dTCMinus",i);
+ fMassHistTCMinus[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
+ fMassHistTCMinus[index]->Sumw2();
hisname.Form("hLSPt%dLCnw",i);
- fMassHistLS[indexLS]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHistLS[indexLS]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHistLS[indexLS]->Sumw2();
hisname.Form("hLSPt%dTCnw",i);
- fMassHistLSTC[indexLS]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHistLSTC[indexLS]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHistLSTC[indexLS]->Sumw2();
hisname.Form("hCosPSigPt%dLS",i);
- fCosPHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.5,1.);
+ fCosPHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.5,1.);
fCosPHistLS[index]->Sumw2();
hisname.Form("hDLenSigPt%dLS",i);
- fDLenHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.5);
+ fDLenHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.5);
fDLenHistLS[index]->Sumw2();
hisname.Form("hSumd02SigPt%dLS",i);
- fSumd02HistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,1.);
+ fSumd02HistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,1.);
fSumd02HistLS[index]->Sumw2();
hisname.Form("hSigVertSigPt%dLS",i);
- fSigVertHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.1);
+ fSigVertHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.1);
fSigVertHistLS[index]->Sumw2();
hisname.Form("hPtMaxSigPt%dLS",i);
- fPtMaxHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.5,5.);
+ fPtMaxHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.5,5.);
fPtMaxHistLS[index]->Sumw2();
hisname.Form("hDCASigPt%dLS",i);
- fDCAHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.1);
+ fDCAHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.1);
fDCAHistLS[index]->Sumw2();
index=GetBackgroundHistoIndex(i);
indexLS++;
hisname.Form("hBkgPt%d",i);
- fMassHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHist[index]->Sumw2();
hisname.Form("hCosPBkgPt%d",i);
- fCosPHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.5,1.);
+ fCosPHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.5,1.);
fCosPHist[index]->Sumw2();
hisname.Form("hDLenBkgPt%d",i);
- fDLenHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.5);
+ fDLenHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.5);
fDLenHist[index]->Sumw2();
hisname.Form("hSumd02BkgPt%d",i);
- fSumd02Hist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,1.);
+ fSumd02Hist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,1.);
fSumd02Hist[index]->Sumw2();
hisname.Form("hSigVertBkgPt%d",i);
- fSigVertHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.1);
+ fSigVertHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.1);
fSigVertHist[index]->Sumw2();
hisname.Form("hPtMaxBkgPt%d",i);
- fPtMaxHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.5,5.);
+ fPtMaxHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.5,5.);
fPtMaxHist[index]->Sumw2();
hisname.Form("hDCABkgPt%d",i);
- fDCAHist[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.1);
+ fDCAHist[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.1);
fDCAHist[index]->Sumw2();
hisname.Form("hBkgPt%dTC",i);
- fMassHistTC[index]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHistTC[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHistTC[index]->Sumw2();
+ hisname.Form("hBkgPt%dTCPlus",i);
+ fMassHistTCPlus[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
+ fMassHistTCPlus[index]->Sumw2();
+ hisname.Form("hBkgPt%dTCMinus",i);
+ fMassHistTCMinus[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
+ fMassHistTCMinus[index]->Sumw2();
hisname.Form("hLSPt%dLCntrip",i);
- fMassHistLS[indexLS]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHistLS[indexLS]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHistLS[indexLS]->Sumw2();
hisname.Form("hLSPt%dTCntrip",i);
- fMassHistLSTC[indexLS]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHistLSTC[indexLS]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHistLSTC[indexLS]->Sumw2();
hisname.Form("hCosPBkgPt%dLS",i);
- fCosPHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.5,1.);
+ fCosPHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.5,1.);
fCosPHistLS[index]->Sumw2();
hisname.Form("hDLenBkgPt%dLS",i);
- fDLenHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.5);
+ fDLenHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.5);
fDLenHistLS[index]->Sumw2();
hisname.Form("hSumd02BkgPt%dLS",i);
- fSumd02HistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,1.);
+ fSumd02HistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,1.);
fSumd02HistLS[index]->Sumw2();
hisname.Form("hSigVertBkgPt%dLS",i);
- fSigVertHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.1);
+ fSigVertHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.1);
fSigVertHistLS[index]->Sumw2();
hisname.Form("hPtMaxBkgPt%dLS",i);
- fPtMaxHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.5,5.);
+ fPtMaxHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.5,5.);
fPtMaxHistLS[index]->Sumw2();
hisname.Form("hDCABkgPt%dLS",i);
- fDCAHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),100,0.,0.1);
+ fDCAHistLS[index]=new TH1F(hisname.Data(),hisname.Data(),nbins,0.,0.1);
fDCAHistLS[index]->Sumw2();
indexLS++;
hisname.Form("hLSPt%dLCntripsinglecut",i);
- fMassHistLS[indexLS]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHistLS[indexLS]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHistLS[indexLS]->Sumw2();
hisname.Form("hLSPt%dTCntripsinglecut",i);
- fMassHistLSTC[indexLS]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHistLSTC[indexLS]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHistLSTC[indexLS]->Sumw2();
indexLS++;
hisname.Form("hLSPt%dLCspc",i);
- fMassHistLS[indexLS]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHistLS[indexLS]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHistLS[indexLS]->Sumw2();
hisname.Form("hLSPt%dTCspc",i);
- fMassHistLSTC[indexLS]=new TH1F(hisname.Data(),hisname.Data(),100,fLowmasslimit,fUpmasslimit);
+ fMassHistLSTC[indexLS]=new TH1F(hisname.Data(),hisname.Data(),nbins,fLowmasslimit,fUpmasslimit);
fMassHistLSTC[indexLS]->Sumw2();
}
fOutput->Add(fPtMaxHist[i]);
fOutput->Add(fDCAHist[i]);
fOutput->Add(fMassHistTC[i]);
+ fOutput->Add(fMassHistTCPlus[i]);
+ fOutput->Add(fMassHistTCMinus[i]);
}
for(Int_t i=0; i<3*fNPtBins&&fDoLS; i++){
fOutput->Add(fCosPHistLS[i]);
fHistNEvents->Sumw2();
fHistNEvents->SetMinimum(0);
fOutput->Add(fHistNEvents);
-
-
+
+ fPtVsMass=new TH2F("hPtVsMass","PtVsMass (prod. cuts)",nbins,fLowmasslimit,fUpmasslimit,40,0.,20.);
+ fPtVsMassTC=new TH2F("hPtVsMassTC","PtVsMass (analysis cuts)",nbins,fLowmasslimit,fUpmasslimit,40,0.,20.);
+ fYVsPt=new TH2F("hYVsPt","YvsPt (prod. cuts)",40,0.,20.,80,-2.,2.);
+ fYVsPtTC=new TH2F("hYVsPtTC","YvsPt (analysis cuts)",40,0.,20.,80,-2.,2.);
+ fYVsPtSig=new TH2F("hYVsPtSig","YvsPt (MC, only sig., prod. cuts)",40,0.,20.,80,-2.,2.);
+ fYVsPtSigTC=new TH2F("hYVsPtSigTC","YvsPt (MC, only Sig, analysis cuts)",40,0.,20.,80,-2.,2.);
+
+ fOutput->Add(fPtVsMass);
+ fOutput->Add(fPtVsMassTC);
+ fOutput->Add(fYVsPt);
+ fOutput->Add(fYVsPtTC);
+ fOutput->Add(fYVsPtSig);
+ fOutput->Add(fYVsPtSigTC);
if(fFillNtuple){
OpenFile(2); // 2 is the slot number of the ntuple
// heavy flavor candidates association to MC truth
AliAODEvent *aod = dynamic_cast<AliAODEvent*> (InputEvent());
- fHistNEvents->Fill(0); // count event
- // Post the data already here
- PostData(1,fOutput);
TClonesArray *array3Prong = 0;
TClonesArray *arrayLikeSign =0;
return;
}
+
+ // fix for temporary bug in ESDfilter
+ // the AODs with null vertex pointer didn't pass the PhysSel
+ if(!aod->GetPrimaryVertex() || TMath::Abs(aod->GetMagneticField())<0.001) return;
+
+ fHistNEvents->Fill(0); // count event
+ // Post the data already here
+ PostData(1,fOutput);
TClonesArray *arrayMC=0;
AliAODMCHeader *mcHeader=0;
}
}
Double_t invMass=d->InvMassDplus();
-
+ Double_t rapid=d->YDplus();
+ fYVsPt->Fill(ptCand,rapid);
+ if(passTightCuts) fYVsPtTC->Fill(ptCand,rapid);
+ Bool_t isFidAcc=fRDCutsAnalysis->IsInFiducialAcceptance(ptCand,rapid);
+ if(isFidAcc){
+ fPtVsMass->Fill(invMass,ptCand);
+ if(passTightCuts) fPtVsMassTC->Fill(invMass,ptCand);
+ }
Float_t tmp[24];
if(fFillNtuple){
tmp[0]=pdgCode;
Double_t dlen=d->DecayLength();
Double_t cosp=d->CosPointingAngle();
Double_t sumD02=d->Getd0Prong(0)*d->Getd0Prong(0)+d->Getd0Prong(1)*d->Getd0Prong(1)+d->Getd0Prong(2)*d->Getd0Prong(2);
- Double_t dca=d->GetDCA();
-Double_t ptmax=0;
+ Double_t dca=d->GetDCA();
+ Double_t sigvert=d->GetSigmaVert();
+ Double_t ptmax=0;
for(Int_t i=0;i<3;i++){
if(d->PtProng(i)>ptmax)ptmax=d->PtProng(i);
}
if(iPtBin>=0){
index=GetHistoIndex(iPtBin);
- fMassHist[index]->Fill(invMass);
- fCosPHist[index]->Fill(cosp);
- fDLenHist[index]->Fill(dlen);
- fSumd02Hist[index]->Fill(sumD02);
- fPtMaxHist[index]->Fill(ptmax);
- fDCAHist[index]->Fill(dca);
-
- if(passTightCuts==1){
- fMassHistTC[index]->Fill(invMass);
+ if(isFidAcc){
+ fMassHist[index]->Fill(invMass);
+ fCosPHist[index]->Fill(cosp);
+ fDLenHist[index]->Fill(dlen);
+ fSumd02Hist[index]->Fill(sumD02);
+ fSigVertHist[index]->Fill(sigvert);
+ fPtMaxHist[index]->Fill(ptmax);
+ fDCAHist[index]->Fill(dca);
+
+ if(passTightCuts){
+ fMassHistTC[index]->Fill(invMass);
+ if(d->GetCharge()>0) fMassHistTCPlus[index]->Fill(invMass);
+ else if(d->GetCharge()<0) fMassHistTCMinus[index]->Fill(invMass);
+ }
}
if(fReadMC){
if(labDp>=0) {
index=GetSignalHistoIndex(iPtBin);
- fMassHist[index]->Fill(invMass);
- fCosPHist[index]->Fill(cosp);
- fDLenHist[index]->Fill(dlen);
- fSumd02Hist[index]->Fill(sumD02);
- fPtMaxHist[index]->Fill(ptmax);
- fDCAHist[index]->Fill(dca);
- if(passTightCuts==1){
- fMassHistTC[index]->Fill(invMass);
-
- }
-
+ if(isFidAcc){
+ fMassHist[index]->Fill(invMass);
+ fCosPHist[index]->Fill(cosp);
+ fDLenHist[index]->Fill(dlen);
+ fSumd02Hist[index]->Fill(sumD02);
+ fSigVertHist[index]->Fill(sigvert);
+ fPtMaxHist[index]->Fill(ptmax);
+ fDCAHist[index]->Fill(dca);
+ if(passTightCuts){
+ fMassHistTC[index]->Fill(invMass);
+ if(d->GetCharge()>0) fMassHistTCPlus[index]->Fill(invMass);
+ else if(d->GetCharge()<0) fMassHistTCMinus[index]->Fill(invMass);
+ }
+ }
+ fYVsPtSig->Fill(ptCand,rapid);
+ if(passTightCuts) fYVsPtSigTC->Fill(ptCand,rapid);
}else{
index=GetBackgroundHistoIndex(iPtBin);
- fMassHist[index]->Fill(invMass);
- fCosPHist[index]->Fill(cosp);
- fDLenHist[index]->Fill(dlen);
- fSumd02Hist[index]->Fill(sumD02);
- fPtMaxHist[index]->Fill(ptmax);
- fDCAHist[index]->Fill(dca);
- if(passTightCuts==1){
- fMassHistTC[index]->Fill(invMass);
-
+ if(isFidAcc){
+ fMassHist[index]->Fill(invMass);
+ fCosPHist[index]->Fill(cosp);
+ fDLenHist[index]->Fill(dlen);
+ fSumd02Hist[index]->Fill(sumD02);
+ fSigVertHist[index]->Fill(sigvert);
+ fPtMaxHist[index]->Fill(ptmax);
+ fDCAHist[index]->Fill(dca);
+ if(passTightCuts){
+ fMassHistTC[index]->Fill(invMass);
+ if(d->GetCharge()>0) fMassHistTCPlus[index]->Fill(invMass);
+ else if(d->GetCharge()<0) fMassHistTCMinus[index]->Fill(invMass);
+ }
}
}
}
- }
+ }
}
if(unsetvtx) d->UnsetOwnPrimaryVtx();
}
printf("ERROR: fOutput not available\n");
return;
}
- fHistNEvents = dynamic_cast<TH1F*>(fOutput->FindObject("fHistNEvents"));
+ fHistNEvents = dynamic_cast<TH1F*>(fOutput->FindObject("fHistNEvents"));
+ fYVsPt = dynamic_cast<TH2F*>(fOutput->FindObject("hYVsPt"));
+ fYVsPtTC = dynamic_cast<TH2F*>(fOutput->FindObject("hYVsPtTC"));
+ fYVsPtSig = dynamic_cast<TH2F*>(fOutput->FindObject("hYVsPtSig"));
+ fYVsPtSigTC = dynamic_cast<TH2F*>(fOutput->FindObject("hYVsPtSigTC"));
+ fPtVsMass = dynamic_cast<TH2F*>(fOutput->FindObject("hPtVsMass"));
+ fPtVsMassTC = dynamic_cast<TH2F*>(fOutput->FindObject("hPtVsMassTC"));
- TString hisname;
- Int_t index=0;
+ TString hisname;
+ Int_t index=0;
- Int_t indexLS=0;
- for(Int_t i=0;i<fNPtBins;i++){
+ Int_t indexLS=0;
+ for(Int_t i=0;i<fNPtBins;i++){
index=GetHistoIndex(i);
if(fDoLS)indexLS=GetLSHistoIndex(i);
hisname.Form("hMassPt%d",i);
fMassHist[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
- hisname.Form("hCosPAllPt%d",i);
+ hisname.Form("hCosPAllPt%d",i);
fCosPHist[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
- hisname.Form("hDLenAllPt%d",i);
+ hisname.Form("hDLenAllPt%d",i);
fDLenHist[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
- hisname.Form("hSumd02AllPt%d",i);
- fSumd02Hist[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
- hisname.Form("hSigVertAllPt%d",i);
- fSigVertHist[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
- hisname.Form("hPtMaxAllPt%d",i);
- fPtMaxHist[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
- hisname.Form("hDCAAllPt%d",i);
- fDCAHist[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
- hisname.Form("hMassPt%dTC",i);
+ hisname.Form("hSumd02AllPt%d",i);
+ fSumd02Hist[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
+ hisname.Form("hSigVertAllPt%d",i);
+ fSigVertHist[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
+ hisname.Form("hPtMaxAllPt%d",i);
+ fPtMaxHist[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
+ hisname.Form("hDCAAllPt%d",i);
+ fDCAHist[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
+ hisname.Form("hMassPt%dTC",i);
fMassHistTC[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
+ hisname.Form("hMassPt%dTCPlus",i);
+ fMassHistTCPlus[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
+ hisname.Form("hMassPt%dTCMinus",i);
+ fMassHistTCMinus[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
if(fDoLS){
hisname.Form("hLSPt%dLC",i);
fMassHistLS[indexLS]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
hisname.Form("hSigPt%dTC",i);
fMassHistTC[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
+ hisname.Form("hSigPt%dTCPlus",i);
+ fMassHistTCPlus[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
+ hisname.Form("hSigPt%dTCMinus",i);
+ fMassHistTCMinus[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
if(fDoLS){
hisname.Form("hLSPt%dLCnw",i);
fMassHistLS[indexLS]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
hisname.Form("hLSPt%dTCnw",i);
fMassHistLSTC[indexLS]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
- }
+ }
index=GetBackgroundHistoIndex(i);
if(fDoLS)indexLS++;
fDCAHist[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
hisname.Form("hBkgPt%dTC",i);
fMassHistTC[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
+ hisname.Form("hBkgPt%dTCPlus",i);
+ fMassHistTCPlus[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
+ hisname.Form("hBkgPt%dTCMinus",i);
+ fMassHistTCMinus[index]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
if(fDoLS){
hisname.Form("hLSPt%dLCntrip",i);
fMassHistLS[indexLS]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
fMassHistLSTC[indexLS]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
}
- }
+ }
if(fFillNtuple){
fNtupleDplus = dynamic_cast<TNtuple*>(GetOutputData(3));
hMassPt->SetXTitle("M[GeV/c^{2}]");
hMassPt->Draw();
- return;
+ return;
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff