fh3PtHPtJDPhi(0),
fh3PtJet1VsMassVsHPtAllSel(0),
fh3PtJet1VsMassVsHPtAllSelMatch(0),
fh3PtHPtJDPhi(0),
fh3PtJet1VsMassVsHPtAllSel(0),
fh3PtJet1VsMassVsHPtAllSelMatch(0),
fh3PtHPtJDPhi = new TH3F*[fNcentBins];
fh3PtJet1VsMassVsHPtAllSel = new TH3F*[fNcentBins];
fh3PtJet1VsMassVsHPtAllSelMatch = new TH3F*[fNcentBins];
fh3PtHPtJDPhi = new TH3F*[fNcentBins];
fh3PtJet1VsMassVsHPtAllSel = new TH3F*[fNcentBins];
fh3PtJet1VsMassVsHPtAllSelMatch = new TH3F*[fNcentBins];
fh3PtHPtJDPhi(0),
fh3PtJet1VsMassVsHPtAllSel(0),
fh3PtJet1VsMassVsHPtAllSelMatch(0),
fh3PtHPtJDPhi(0),
fh3PtJet1VsMassVsHPtAllSel(0),
fh3PtJet1VsMassVsHPtAllSelMatch(0),
fh3PtHPtJDPhi = new TH3F*[fNcentBins];
fh3PtJet1VsMassVsHPtAllSel = new TH3F*[fNcentBins];
fh3PtJet1VsMassVsHPtAllSelMatch = new TH3F*[fNcentBins];
fh3PtHPtJDPhi = new TH3F*[fNcentBins];
fh3PtJet1VsMassVsHPtAllSel = new TH3F*[fNcentBins];
fh3PtJet1VsMassVsHPtAllSelMatch = new TH3F*[fNcentBins];
fh1PtHadron[i] = new TH1F(histName.Data(),histTitle.Data(),200.,0.,200.);
fOutput->Add(fh1PtHadron[i]);
fh1PtHadron[i] = new TH1F(histName.Data(),histTitle.Data(),200.,0.,200.);
fOutput->Add(fh1PtHadron[i]);
+ histName = TString::Format("fh1PtHadronMatch_%d",i);
+ histTitle = TString::Format("%s;#it{p}_{T,h}",histName.Data());
+ fh1PtHadronMatch[i] = new TH1F(histName.Data(),histTitle.Data(),200.,0.,200.);
+ fOutput->Add(fh1PtHadronMatch[i]);
+
histName = TString::Format("fh3PtHPtJDPhi_%d",i);
histTitle = TString::Format("%s;#it{p}_{T,h};#it{p}_{T,jet};#Delta#varphi_{h,jet}",histName.Data());
fh3PtHPtJDPhi[i] = new TH3F(histName.Data(),histTitle.Data(),nBinsPtH,minPtH,maxPtH,nBinsPt,minPt,maxPt,nBinsPhi,minPhi,maxPhi);
histName = TString::Format("fh3PtHPtJDPhi_%d",i);
histTitle = TString::Format("%s;#it{p}_{T,h};#it{p}_{T,jet};#Delta#varphi_{h,jet}",histName.Data());
fh3PtHPtJDPhi[i] = new TH3F(histName.Data(),histTitle.Data(),nBinsPtH,minPtH,maxPtH,nBinsPt,minPt,maxPt,nBinsPhi,minPhi,maxPhi);
vp = GetSingleInclusiveTT(pCont,fPtTTMin->At(it),fPtTTMax->At(it));
if(!vp) continue;
fh1PtHadron[fCentBin]->Fill(vp->Pt()); //all trigger tracks
vp = GetSingleInclusiveTT(pCont,fPtTTMin->At(it),fPtTTMax->At(it));
if(!vp) continue;
fh1PtHadron[fCentBin]->Fill(vp->Pt()); //all trigger tracks
- if(vp->TestBits(fEmbConstSel) != (Int_t)fEmbConstSel)
- AliDebug(11,Form("Skipping track with pT=%f because it does not match the bit mask (%d, %d)", vp->Pt(), fEmbConstSel, vp->TestBits(fEmbConstSel)));
- else {
+ if(fMarkMCLabel>0 && TMath::Abs(vp->GetLabel()) >= fMarkMCLabel ) {
+ // Printf("Accepting track with pT=%f because it does match the bit mask (%d, %d, %d). label: %d fMarkMCLabel: %d", vp->Pt(), fEmbConstSel, vp->TestBits(fEmbConstSel), vp->TestBit(fEmbConstSel),vp->GetLabel(),fMarkMCLabel);
fh3PtJet1VsMassVsHPtAllSelMatch[fCentBin]->Fill(ptJet,mJet,pt);
fh3PtJet1VsRatVsHPtAllSelMatch[fCentBin]->Fill(ptJet,rat,pt);
}
fh3PtJet1VsMassVsHPtAllSelMatch[fCentBin]->Fill(ptJet,mJet,pt);
fh3PtJet1VsRatVsHPtAllSelMatch[fCentBin]->Fill(ptJet,rat,pt);
}
fh3PtJet1VsRatVsHPtTagged[fCentBin]->Fill(ptJet,rat,pt);
if(fMinFractionShared>0. && fraction>fMinFractionShared) {
fh3PtJet1VsRatVsHPtTagged[fCentBin]->Fill(ptJet,rat,pt);
if(fMinFractionShared>0. && fraction>fMinFractionShared) {
- if(vp->TestBits(fEmbConstSel) != (Int_t)fEmbConstSel)
- AliDebug(11,Form("Skipping track with pT=%f because it does not match the bit mask (%d, %d)", vp->Pt(), fEmbConstSel, vp->TestBits(fEmbConstSel)));
- else {
+ if(fMarkMCLabel>0 && TMath::Abs(vp->GetLabel()) >= fMarkMCLabel ) {
fh3PtJet1VsMassVsHPtTaggedMatch[fCentBin]->Fill(ptJet,mJet,pt);
fh3PtJet1VsRatVsHPtTaggedMatch[fCentBin]->Fill(ptJet,rat,pt);
}
fh3PtJet1VsMassVsHPtTaggedMatch[fCentBin]->Fill(ptJet,mJet,pt);
fh3PtJet1VsRatVsHPtTaggedMatch[fCentBin]->Fill(ptJet,rat,pt);
}