SetCurrentOutputList(0);
- // Cuts
+ TH1* tmpHisto = AddHistogram1D<TH1D>("hVertexAcceptance", "Accepted vertices for different conditions", "", 4, 0, 4, "stage","N^{Events}/cut");
+ tmpHisto->GetXaxis()->SetBinLabel(1, "Triggered all");
+ tmpHisto->GetXaxis()->SetBinLabel(2, "Triggered w/ vertex");
+ tmpHisto->GetXaxis()->SetBinLabel(3, "Pile-up corrected all");
+ tmpHisto->GetXaxis()->SetBinLabel(4, "Pile-up corrected w vertex");
+
TH2* tmpHisto2D = AddHistogram2D<TH2D>("hCentrality", Form("Accepted events in centrality (%s)", fCentralityType.Data()), "COLZ", 102, 0., 102., 4, 0,4,"Centrality","Cut stage","dN^{Events}");
tmpHisto2D->GetYaxis()->SetBinLabel(1, "Before cuts");
tmpHisto2D->GetYaxis()->SetBinLabel(2, "After pile up");
tmpHisto2D->GetYaxis()->SetBinLabel(3, "After vertex demand");
tmpHisto2D->GetYaxis()->SetBinLabel(4, "After vertex cuts");
- TH1* tmpHisto = AddHistogram1D<TH1D>("hTrackAcceptance", "Accepted tracks (0 = before cuts, 1 = after eta, 2 = after pT)", "", 3, 0, 3, "stage","N^{Tracks}/cut");
+ tmpHisto = AddHistogram1D<TH1D>("hTrackAcceptance", "Accepted tracks (0 = before cuts, 1 = after eta, 2 = after pT)", "", 3, 0, 3, "stage","N^{Tracks}/cut");
tmpHisto->GetXaxis()->SetBinLabel(1, "Before cuts");
tmpHisto->GetXaxis()->SetBinLabel(2, "After eta");
tmpHisto->GetXaxis()->SetBinLabel(3, "After p_{T}");
AddHistogram1D<TH1D>("hLeadingJetDeltaPhi", "1st and 2nd leading jet #Delta #phi", "", 250, 0., TMath::Pi(), "#Delta #phi","dN^{Jets}/d(#Delta #phi)");
// Background distributions
+
AddHistogram2D<TH2D>("hKTBackgroundExternal", "KT background density (External task)", "LEGO2", 400, 0., 40., fNumberOfCentralityBins, 0, 100, "#rho (GeV/c)","Centrality", "dN^{Events}/d#rho");
+ AddHistogram2D<TH2D>("hKTBackgroundExternalVsPt", "KT background density (External task)", "LEGO2", 400, 0., 40., 200, 0, 200, "#rho (GeV/c)","Raw jet p_{T}", "dN^{Events}/d#rho");
AddHistogram2D<TH2D>("hKTBackgroundExternal20GeV", "KT background density (External task, jet p_{T} > 20 GeV)", "LEGO2", 400, 0., 40., fNumberOfCentralityBins, 0, 100, "#rho (GeV/c)","Centrality", "dN^{Events}/d#rho");
AddHistogram2D<TH2D>("hKTBackgroundImprovedCMS", "KT background density (Improved CMS approach)", "LEGO2", 400, 0., 40., fNumberOfCentralityBins, 0, 100, "#rho (GeV/c)","Centrality", "dN^{Events}/d#rho");
AddHistogram2D<TH2D>("hPPBackground", "PP background density (Michals approach)", "LEGO2", 400, 0., 40., fNumberOfCentralityBins, 0, 100, "#rho (GeV/c)","Centrality", "dN^{Events}/d#rho");
AddHistogram2D<TH2D>("hDeltaPtPP", "Background fluctuations #delta p_{T} (PP approach)", "", 1801, -40.0, 80.0, fNumberOfCentralityBins, 0, 100, "#delta p_{T} (GeV/c)","Centrality","dN^{Jets}/d#delta p_{T}");
AddHistogram2D<TH2D>("hDeltaPtExternalBgrd", "Background fluctuations #delta p_{T} (KT, External)", "", 1801, -40.0, 80.0, fNumberOfCentralityBins, 0, 100, "#delta p_{T} (GeV/c)","Centrality","dN^{Jets}/d#delta p_{T}");
AddHistogram2D<TH2D>("hDeltaPtExternalBgrdVsPt", "Background fluctuations #delta p_{T} (KT, External, in p_{T} bins)", "", 1801, -40.0, 80.0, 200, 0, 200, "#delta p_{T} (GeV/c)","Raw jet p_{T}","dN^{Jets}/d#delta p_{T}");
+ AddHistogram2D<TH2D>("hDeltaPtExternalBgrdPartialExclusion", "Background fluctuations #delta p_{T} (KT, External, partial jet exclusion)", "", 1801, -40.0, 80.0, fNumberOfCentralityBins, 0, 100, "#delta p_{T} (GeV/c)","Centrality","dN^{Jets}/d#delta p_{T}");
AddHistogram2D<TH2D>("hDeltaPtKTImprovedCMS", "Background fluctuations #delta p_{T} (KT, Improved CMS-like)", "", 1801, -40.0, 80.0, fNumberOfCentralityBins, 0, 100, "#delta p_{T} (GeV/c)","Centrality","dN^{Jets}/d#delta p_{T}");
- AddHistogram2D<TH2D>("hDeltaPtKTImprovedCMSFullExclusion", "Background fluctuations #delta p_{T} (KT, Improved CMS-like, full leading jet exclusion)", "", 1801, -40.0, 80.0, fNumberOfCentralityBins, 0, 100, "#delta p_{T} (GeV/c)","Centrality","dN^{Jets}/d#delta p_{T}");
AddHistogram2D<TH2D>("hDeltaPtNoBackground", "Background fluctuations #delta p_{T} (No background)", "", 1801, -40.0, 80.0, fNumberOfCentralityBins, 0, 100, "#delta p_{T} (GeV/c)","Centrality","dN^{Jets}/d#delta p_{T}");
AddHistogram2D<TH2D>("hDeltaPtKTPbPb", "Background fluctuations #delta p_{T} (KT, PbPb w/o ghosts)", "", 1801, -40.0, 80.0, fNumberOfCentralityBins, 0, 100, "#delta p_{T} (GeV/c)","Centrality","dN^{Jets}/d#delta p_{T}");
AddHistogram2D<TH2D>("hDeltaPtKTPbPbWithGhosts", "Background fluctuations #delta p_{T} (KT, PbPb w/ ghosts)", "", 1801, -40.0, 80.0, fNumberOfCentralityBins, 0, 100, "#delta p_{T} (GeV/c)","Centrality","dN^{Jets}/d#delta p_{T}");
}
//________________________________________________________________________
-AliAnalysisTaskChargedJetsPA::AliAnalysisTaskChargedJetsPA(const char *name, const char* trackArrayName, const char* jetArrayName, const char* backgroundJetArrayName, Bool_t analyzeJetProfile, Bool_t analyzeTrackcuts) : AliAnalysisTaskSE(name), fOutputLists(), fCurrentOutputList(0), fDoJetAnalysis(1), fAnalyzeJetProfile(0), fAnalyzeTrackcuts(0), fAnalyzeJetConstituents(1), fParticleLevel(0), fUseDefaultVertexCut(1), fUsePileUpCut(1), fSetCentralityToOne(0), fNoExternalBackground(0), fBackgroundForJetProfile(0), fPartialAnalysisNParts(1), fPartialAnalysisIndex(0), fJetArray(0), fTrackArray(0), fBackgroundJetArray(0), fJetArrayName(), fTrackArrayName(), fBackgroundJetArrayName(), fRhoTaskName(), fRandConeRadius(0.4), fSignalJetRadius(0.4), fBackgroundJetRadius(0.4), fNumberExcludedJets(-1), fMinEta(-0.9), fMaxEta(0.9), fMinJetEta(-0.5), fMaxJetEta(0.5), fMinTrackPt(0.150), fMinJetPt(5.0), fMinJetArea(0.5), fMinBackgroundJetPt(0.0), fMinNCrossedRows(70), fUsePtDepCrossedRowsCut(0), fNumberOfCentralityBins(20), fCentralityType("V0A"), fMatchTr(), fMatchChi(), fPrimaryVertex(0), fFirstLeadingJet(0), fSecondLeadingJet(0), fFirstLeadingKTJet(0), fSecondLeadingKTJet(0), fNumberSignalJets(0), fNumberSignalJetsAbove5GeV(0), fRandom(0), fHelperClass(0), fInitialized(0), fTaskInstanceCounter(0), fIsDEBUG(0), fIsPA(1), fNoTerminate(1), fEventCounter(0), fHybridESDtrackCuts(0), fHybridESDtrackCuts_variedPtDep(0), fHybridESDtrackCuts_variedPtDep2(0)
+AliAnalysisTaskChargedJetsPA::AliAnalysisTaskChargedJetsPA(const char *name, const char* trackArrayName, const char* jetArrayName, const char* backgroundJetArrayName, Bool_t analyzeJetProfile, Bool_t analyzeTrackcuts) : AliAnalysisTaskSE(name), fOutputLists(), fCurrentOutputList(0), fDoJetAnalysis(1), fAnalyzeJetProfile(0), fAnalyzeTrackcuts(0), fAnalyzeJetConstituents(1), fParticleLevel(0), fUseDefaultVertexCut(1), fUsePileUpCut(1), fSetCentralityToOne(0), fNoExternalBackground(0), fBackgroundForJetProfile(0), fPartialAnalysisNParts(1), fPartialAnalysisIndex(0), fJetArray(0), fTrackArray(0), fBackgroundJetArray(0), fJetArrayName(), fTrackArrayName(), fBackgroundJetArrayName(), fRhoTaskName(), fRandConeRadius(0.4), fRandConeNumber(10), fSignalJetRadius(0.4), fBackgroundJetRadius(0.4), fNumberExcludedJets(-1), fMinEta(-0.9), fMaxEta(0.9), fMinJetEta(-0.5), fMaxJetEta(0.5), fMinTrackPt(0.150), fMinJetPt(5.0), fMinJetArea(0.5), fMinBackgroundJetPt(0.0), fMinNCrossedRows(70), fUsePtDepCrossedRowsCut(0), fNumberOfCentralityBins(20), fCentralityType("V0A"), fMatchTr(), fMatchChi(), fPrimaryVertex(0), fFirstLeadingJet(0), fSecondLeadingJet(0), fFirstLeadingKTJet(0), fSecondLeadingKTJet(0), fNumberSignalJets(0), fNumberSignalJetsAbove5GeV(0), fRandom(0), fHelperClass(0), fInitialized(0), fTaskInstanceCounter(0), fIsDEBUG(0), fIsPA(1), fNoTerminate(1), fEventCounter(0), fTempExcludedRCs(0), fTempAllRCs(1), fTempOverlapCounter(0), fTempMeanExclusionProbability(0), fHybridESDtrackCuts(0), fHybridESDtrackCuts_variedPtDep(0), fHybridESDtrackCuts_variedPtDep2(0)
{
#ifdef DEBUGMODE
AliInfo("Calling constructor.");
if (!event)
return kFALSE;
- // ### Get centrality values
+ fPrimaryVertex = event->GetPrimaryVertex();
+
+ // ### Create plot for vertex acceptance
+ fHelperClass->SetMaxVtxZ(10000.);
+ // Check vertex existance
+ Bool_t hasVertex = kFALSE;
+ if(fUseDefaultVertexCut)
+ {
+ if(fHelperClass->IsVertexSelected2013pA(event))
+ hasVertex = kTRUE;
+ }
+ else
+ {
+ if(!fPrimaryVertex || (fPrimaryVertex->GetNContributors()<2) || (TMath::Sqrt(fPrimaryVertex->GetX()*fPrimaryVertex->GetX() + fPrimaryVertex->GetY()*fPrimaryVertex->GetY()) > 1.0))
+ hasVertex = kTRUE;
+ }
+
+ // All triggered events
+ FillHistogram("hVertexAcceptance", 0.5); // all triggered events all
+ if(hasVertex)
+ FillHistogram("hVertexAcceptance", 1.5); // all triggered events w/ vertex
+
+ // Pile-up corrected events
+ if(!fHelperClass->IsPileUpEvent(event))
+ {
+ FillHistogram("hVertexAcceptance", 2.5); // pile-up corr. events all
+ if(hasVertex)
+ FillHistogram("hVertexAcceptance", 3.5); // pile-up corr. events w/ vertex
+ }
+
+ fHelperClass->SetMaxVtxZ(10.);
+
+ // ### Get centrality values
AliCentrality* tmpCentrality = event->GetCentrality();
Double_t centralityPercentile = -1.0;
if (tmpCentrality != NULL)
FillHistogram("hCentrality", centralityPercentile, 1.5); // after pileup cut
// ### CUT STAGE 2: Existence of primary vertex
- fPrimaryVertex = event->GetPrimaryVertex();
- if(fUseDefaultVertexCut)
- {
- fHelperClass->SetMaxVtxZ(10000.);
- if(!fHelperClass->IsVertexSelected2013pA(event))
- {
- fHelperClass->SetMaxVtxZ(10.);
- return kFALSE;
- }
- fHelperClass->SetMaxVtxZ(10.);
- }
- else // Vertex cut for pp
- {
- if(!fPrimaryVertex || (fPrimaryVertex->GetNContributors()<2) || (TMath::Sqrt(fPrimaryVertex->GetX()*fPrimaryVertex->GetX() + fPrimaryVertex->GetY()*fPrimaryVertex->GetY()) > 1.0))
- return kFALSE;
- }
+ if(!hasVertex)
+ return kFALSE;
FillHistogram("hVertexZ",fPrimaryVertex->GetZ());
FillHistogram("hCentrality", centralityPercentile, 2.5); // after vertex existance cut
//________________________________________________________________________
-Double_t AliAnalysisTaskChargedJetsPA::GetDeltaPt(Double_t rho, Double_t leadingJetExclusionProbability)
+Double_t AliAnalysisTaskChargedJetsPA::GetDeltaPt(Double_t rho, Double_t overlappingJetExclusionProbability)
{
#ifdef DEBUGMODE
AliInfo("Getting Delta Pt.");
// Define an invalid delta pt
Double_t deltaPt = -10000.0;
+ // Define the ratio of excluded RCs over all RCs
+ Double_t ratioExcludedRCs = static_cast<Double_t>(fTempExcludedRCs)/fTempAllRCs;
+
// Define eta range
Double_t etaMin, etaMax;
etaMin = fMinEta+fRandConeRadius;
Double_t tmpRandConeEta = etaMin + fRandom->Rndm()*(etaMax-etaMin);
Double_t tmpRandConePhi = fRandom->Rndm()*TMath::TwoPi();
- // if there is a jet, check for overlap if demanded
- if(leadingJetExclusionProbability)
+ // Check if a signal jet is overlapping with random cone
+ if(overlappingJetExclusionProbability)
{
- AliEmcalJet* tmpLeading = dynamic_cast<AliEmcalJet*>(fJetArray->At(0));
- // Get leading jet (regardless of pT)
- for (Int_t i = 1; i<fJetArray->GetEntries(); i++)
+ // Calculate the mean exclusion probability
+ fTempOverlapCounter++;
+ fTempMeanExclusionProbability += overlappingJetExclusionProbability;
+ // For all jets, check overlap
+ for (Int_t i = 0; i<fJetArray->GetEntries(); i++)
{
AliEmcalJet* tmpJet = static_cast<AliEmcalJet*>(fJetArray->At(i));
- // if jet is in acceptance and higher, take as new leading
- if (tmpJet)
- if ( ((tmpJet->Eta() >= fMinJetEta) && (tmpJet->Eta() < fMaxJetEta)) && (tmpJet->Area() >= fMinJetArea))
- if((!tmpLeading) || (tmpJet->Pt() > tmpLeading->Pt()))
- tmpLeading = tmpJet;
- }
- if(tmpLeading)
- {
- Double_t excludedJetPhi = tmpLeading->Phi();
- Double_t excludedJetEta = tmpLeading->Eta();
- Double_t tmpDeltaPhi = GetDeltaPhi(tmpRandConePhi, excludedJetPhi);
+ if (!IsSignalJetInAcceptance(tmpJet, kTRUE)) continue;
- // Check, if cone has overlap with jet
- if ( tmpDeltaPhi*tmpDeltaPhi + TMath::Abs(tmpRandConeEta-excludedJetEta)*TMath::Abs(tmpRandConeEta-excludedJetEta) <= fRandConeRadius*fRandConeRadius)
+ // Check overlap
+ Double_t tmpDeltaPhi = GetDeltaPhi(tmpRandConePhi, tmpJet->Phi());
+ if ( tmpDeltaPhi*tmpDeltaPhi + (tmpRandConeEta-tmpJet->Eta())*(tmpRandConeEta-tmpJet->Eta()) <= fRandConeRadius*fRandConeRadius )
{
- // Define probability to exclude the RC
- Double_t probability = leadingJetExclusionProbability;
-
- // Only exclude cone with a given probability
- if (fRandom->Rndm()<=probability)
+ fTempAllRCs++;
+ // If an overlap is given, discard or accept it according to the exclusion prob.
+ if(ratioExcludedRCs < fTempMeanExclusionProbability/fTempOverlapCounter) // to few RCs excluded -> exclude this one
+ {
coneValid = kFALSE;
+ fTempExcludedRCs++;
+ }
+ else // to many RCs excluded -> take this one
+ coneValid = kTRUE;
}
}
}
+
// Get the cones' pt and calculate delta pt
if (coneValid)
deltaPt = GetConePt(tmpRandConeEta,tmpRandConePhi,fRandConeRadius) - (rho*fRandConeRadius*fRandConeRadius*TMath::Pi());
FillHistogram("hJetPtSubtractedRhoExternal", tmpJet->Pt(), centralityPercentile, tmpJet->Pt() - GetCorrectedJetPt(tmpJet, backgroundExternal));
FillHistogram("hJetPtSubtractedRhoKTImprovedCMS", tmpJet->Pt(), centralityPercentile, tmpJet->Pt() - GetCorrectedJetPt(tmpJet, backgroundKTImprovedCMS));
FillHistogram("hJetPtSubtractedRhoPP", tmpJet->Pt(), centralityPercentile, tmpJet->Pt() - GetCorrectedJetPt(tmpJet, backgroundPP));
- FillHistogram("hDeltaPtExternalBgrdVsPt", GetDeltaPt(backgroundExternal), GetCorrectedJetPt(tmpJet, backgroundExternal));
+ for(Int_t j=0; j<fRandConeNumber; j++)
+ FillHistogram("hDeltaPtExternalBgrdVsPt", GetDeltaPt(backgroundExternal), GetCorrectedJetPt(tmpJet, backgroundExternal));
+ FillHistogram("hKTBackgroundExternalVsPt", backgroundExternal, GetCorrectedJetPt(tmpJet, backgroundExternal));
// ###### CONSTITUENT ANALYSIS
FillHistogram("hTRBackgroundExact", backgroundTRExact, centralityPercentile);
// Calculate the delta pt
- Double_t tmpDeltaPtNoBackground = GetDeltaPt(0.0);
- Double_t tmpDeltaPtExternalBgrd = GetDeltaPt(backgroundExternal);
- Double_t tmpDeltaPtPP = GetDeltaPt(backgroundPP);
- Double_t tmpDeltaPtKTImprovedCMS = GetDeltaPt(backgroundKTImprovedCMS);
- Double_t tmpDeltaPtKTImprovedCMSFullExclusion = GetDeltaPt(backgroundKTImprovedCMS, 1.0);
-
- Double_t tmpDeltaPtKTPbPb = 0;
- Double_t tmpDeltaPtKTPbPbWithGhosts = 0;
- Double_t tmpDeltaPtKTCMS = 0;
- Double_t tmpDeltaPtKTMean = 0;
- Double_t tmpDeltaPtKTTrackLike = 0;
- Double_t tmpDeltaPtTR = 0;
-
- tmpDeltaPtKTPbPb = GetDeltaPt(backgroundKTPbPb);
- tmpDeltaPtKTPbPbWithGhosts = GetDeltaPt(backgroundKTPbPbWithGhosts);
- tmpDeltaPtKTCMS = GetDeltaPt(backgroundKTCMS);
- tmpDeltaPtKTMean = GetDeltaPt(backgroundKTMean);
- tmpDeltaPtKTTrackLike = GetDeltaPt(backgroundKTTrackLike);
- tmpDeltaPtTR = GetDeltaPt(backgroundTRCone06);
-
-
- // If valid, fill the delta pt histograms
-
- if(tmpDeltaPtExternalBgrd > -10000.0)
- FillHistogram("hDeltaPtExternalBgrd", tmpDeltaPtExternalBgrd, centralityPercentile);
- if(tmpDeltaPtKTImprovedCMS > -10000.0)
- FillHistogram("hDeltaPtKTImprovedCMS", tmpDeltaPtKTImprovedCMS, centralityPercentile);
- if(tmpDeltaPtKTImprovedCMSFullExclusion > -10000.0)
- FillHistogram("hDeltaPtKTImprovedCMSFullExclusion", tmpDeltaPtKTImprovedCMSFullExclusion, centralityPercentile);
- if(tmpDeltaPtPP > -10000.0)
- FillHistogram("hDeltaPtPP", tmpDeltaPtPP, centralityPercentile);
-
- if(tmpDeltaPtNoBackground > -10000.0)
- FillHistogram("hDeltaPtNoBackground", tmpDeltaPtNoBackground, centralityPercentile);
-
- if(tmpDeltaPtKTPbPb > -10000.0)
- FillHistogram("hDeltaPtKTPbPb", tmpDeltaPtKTPbPb, centralityPercentile);
- if(tmpDeltaPtKTPbPbWithGhosts > -10000.0)
- FillHistogram("hDeltaPtKTPbPbWithGhosts", tmpDeltaPtKTPbPbWithGhosts, centralityPercentile);
- if(tmpDeltaPtKTCMS > -10000.0)
- FillHistogram("hDeltaPtKTCMS", tmpDeltaPtKTCMS, centralityPercentile);
- if(tmpDeltaPtKTMean > -10000.0)
- FillHistogram("hDeltaPtKTMean", tmpDeltaPtKTMean, centralityPercentile);
- if(tmpDeltaPtKTTrackLike > -10000.0)
- FillHistogram("hDeltaPtKTTrackLike", tmpDeltaPtKTTrackLike, centralityPercentile);
- if(tmpDeltaPtTR > -10000.0)
- FillHistogram("hDeltaPtTR", tmpDeltaPtTR, centralityPercentile);
+ for(Int_t i=0; i<fRandConeNumber; i++)
+ {
+ Double_t tmpRatio =1./10.;
+ if(fNumberSignalJets)
+ tmpRatio =1./fNumberSignalJets;
+
+ Double_t tmpDeltaPtNoBackground = GetDeltaPt(0.0);
+ Double_t tmpDeltaPtExternalBgrd = GetDeltaPt(backgroundExternal);
+ Double_t tmpDeltaPtExternalBgrdPartialExclusion = GetDeltaPt(backgroundExternal, tmpRatio);
+ Double_t tmpDeltaPtPP = GetDeltaPt(backgroundPP);
+ Double_t tmpDeltaPtKTImprovedCMS = GetDeltaPt(backgroundKTImprovedCMS);
+
+ Double_t tmpDeltaPtKTPbPb = 0;
+ Double_t tmpDeltaPtKTPbPbWithGhosts = 0;
+ Double_t tmpDeltaPtKTCMS = 0;
+ Double_t tmpDeltaPtKTMean = 0;
+ Double_t tmpDeltaPtKTTrackLike = 0;
+ Double_t tmpDeltaPtTR = 0;
+
+ tmpDeltaPtKTPbPb = GetDeltaPt(backgroundKTPbPb);
+ tmpDeltaPtKTPbPbWithGhosts = GetDeltaPt(backgroundKTPbPbWithGhosts);
+ tmpDeltaPtKTCMS = GetDeltaPt(backgroundKTCMS);
+ tmpDeltaPtKTMean = GetDeltaPt(backgroundKTMean);
+ tmpDeltaPtKTTrackLike = GetDeltaPt(backgroundKTTrackLike);
+ tmpDeltaPtTR = GetDeltaPt(backgroundTRCone06);
+
+
+ // If valid, fill the delta pt histograms
+
+ if(tmpDeltaPtExternalBgrd > -10000.0)
+ FillHistogram("hDeltaPtExternalBgrd", tmpDeltaPtExternalBgrd, centralityPercentile);
+ if(tmpDeltaPtKTImprovedCMS > -10000.0)
+ FillHistogram("hDeltaPtKTImprovedCMS", tmpDeltaPtKTImprovedCMS, centralityPercentile);
+ if(tmpDeltaPtExternalBgrdPartialExclusion > -10000.0)
+ FillHistogram("hDeltaPtExternalBgrdPartialExclusion", tmpDeltaPtExternalBgrdPartialExclusion, centralityPercentile);
+ if(tmpDeltaPtPP > -10000.0)
+ FillHistogram("hDeltaPtPP", tmpDeltaPtPP, centralityPercentile);
+
+ if(tmpDeltaPtNoBackground > -10000.0)
+ FillHistogram("hDeltaPtNoBackground", tmpDeltaPtNoBackground, centralityPercentile);
+
+ if(tmpDeltaPtKTPbPb > -10000.0)
+ FillHistogram("hDeltaPtKTPbPb", tmpDeltaPtKTPbPb, centralityPercentile);
+ if(tmpDeltaPtKTPbPbWithGhosts > -10000.0)
+ FillHistogram("hDeltaPtKTPbPbWithGhosts", tmpDeltaPtKTPbPbWithGhosts, centralityPercentile);
+ if(tmpDeltaPtKTCMS > -10000.0)
+ FillHistogram("hDeltaPtKTCMS", tmpDeltaPtKTCMS, centralityPercentile);
+ if(tmpDeltaPtKTMean > -10000.0)
+ FillHistogram("hDeltaPtKTMean", tmpDeltaPtKTMean, centralityPercentile);
+ if(tmpDeltaPtKTTrackLike > -10000.0)
+ FillHistogram("hDeltaPtKTTrackLike", tmpDeltaPtKTTrackLike, centralityPercentile);
+ if(tmpDeltaPtTR > -10000.0)
+ FillHistogram("hDeltaPtTR", tmpDeltaPtTR, centralityPercentile);
+ }
#ifdef DEBUGMODE
AliInfo("Calculate()::Background done.");