From 1e3ad34df85013adbf01cd6c0abddbc5eb2fd363 Mon Sep 17 00:00:00 2001 From: fkrizek Date: Thu, 16 Oct 2014 17:03:08 +0200 Subject: [PATCH] multiple trigger counter histograms --- PWGJE/AliAnalysisTaskJetCorePP.cxx | 62 +++++++++++++++++++++++++++++- PWGJE/AliAnalysisTaskJetCorePP.h | 5 +++ 2 files changed, 66 insertions(+), 1 deletion(-) diff --git a/PWGJE/AliAnalysisTaskJetCorePP.cxx b/PWGJE/AliAnalysisTaskJetCorePP.cxx index 4171d2215fc..173e610550d 100644 --- a/PWGJE/AliAnalysisTaskJetCorePP.cxx +++ b/PWGJE/AliAnalysisTaskJetCorePP.cxx @@ -137,6 +137,8 @@ fhDphiTriggerJet(0x0), fhDphiTriggerJetAccept(0x0), fhCentrality(0x0), fhCentralityAccept(0x0), +fhNofMultipleTriggers(0x0), +fhDeltaRMultTriggers(0x0), //fHJetPtRaw(0x0), //fHLeadingJetPtRaw(0x0), //fHDphiVsJetPtAll(0x0), @@ -165,6 +167,8 @@ fhEntriesToMedian(0x0), fhEntriesToMedianGen(0x0), fhCellAreaToMedian(0x0), fhCellAreaToMedianGen(0x0), +fhNofMultipleTriggersGen(0x0), +fhDeltaRMultTriggersGen(0x0), fIsChargedMC(0), fIsKine(0), fIsFullMC(0), @@ -266,6 +270,8 @@ fhDphiTriggerJet(0x0), fhDphiTriggerJetAccept(0x0), fhCentrality(0x0), fhCentralityAccept(0x0), +fhNofMultipleTriggers(0x0), +fhDeltaRMultTriggers(0x0), //fHJetPtRaw(0x0), //fHLeadingJetPtRaw(0x0), //fHDphiVsJetPtAll(0x0), @@ -294,6 +300,8 @@ fhEntriesToMedian(0x0), fhEntriesToMedianGen(0x0), fhCellAreaToMedian(0x0), fhCellAreaToMedianGen(0x0), +fhNofMultipleTriggersGen(0x0), +fhDeltaRMultTriggersGen(0x0), fIsChargedMC(0), fIsKine(0), fIsFullMC(0), @@ -401,6 +409,8 @@ fhDphiTriggerJet(a.fhDphiTriggerJet), fhDphiTriggerJetAccept(a.fhDphiTriggerJetAccept), fhCentrality(a.fhCentrality), fhCentralityAccept(a.fhCentralityAccept), +fhNofMultipleTriggers(a.fhNofMultipleTriggers), +fhDeltaRMultTriggers(a.fhDeltaRMultTriggers), //fHJetPtRaw(a.fHJetPtRaw), //fHLeadingJetPtRaw(a.fHLeadingJetPtRaw), //fHDphiVsJetPtAll(a.fHDphiVsJetPtAll), @@ -429,6 +439,8 @@ fhEntriesToMedian(a.fhEntriesToMedian), fhEntriesToMedianGen(a.fhEntriesToMedianGen), fhCellAreaToMedian(a.fhCellAreaToMedian), fhCellAreaToMedianGen(a.fhCellAreaToMedianGen), +fhNofMultipleTriggersGen(a.fhNofMultipleTriggersGen), +fhDeltaRMultTriggersGen(a.fhDeltaRMultTriggersGen), fIsChargedMC(a.fIsChargedMC), fIsKine(a.fIsKine), fIsFullMC(a.fIsFullMC), @@ -711,6 +723,9 @@ void AliAnalysisTaskJetCorePP::UserCreateOutputObjects() fhEntriesToMedian = new TH1D("fhEntriesToMedian","fhEntriesToMedian",30,0,30); fhCellAreaToMedian = new TH1D("fhCellAreaToMedian", "fhCellAreaToMedian", 75,0,1.5); + fhNofMultipleTriggers = new TH1D("fhNofMultipleTriggers","fhNofMultipleTriggers",100,0,100); + fhDeltaRMultTriggers = new TH1D("fhDeltaRMultTriggers","fhDeltaRMultTriggers", 100,0,4); + if(!fIsKine){ fOutputList->Add(fhJetPhi); fOutputList->Add(fhTriggerPhi); @@ -728,6 +743,8 @@ void AliAnalysisTaskJetCorePP::UserCreateOutputObjects() fOutputList->Add(fhCentralityAccept); fOutputList->Add(fhEntriesToMedian); fOutputList->Add(fhCellAreaToMedian); + fOutputList->Add(fhNofMultipleTriggers); + fOutputList->Add(fhDeltaRMultTriggers); } // raw spectra of INCLUSIVE jets //Centrality, pTjet, A @@ -855,6 +872,13 @@ void AliAnalysisTaskJetCorePP::UserCreateOutputObjects() fhCellAreaToMedianGen = (TH1D*) fhCellAreaToMedian->Clone("fhCellAreaToMedianGen"); fhCellAreaToMedianGen->SetTitle(Form("%s Gen MC", fhCellAreaToMedian->GetTitle())); fOutputList->Add(fhCellAreaToMedianGen); + + fhNofMultipleTriggersGen = (TH1D*) fhNofMultipleTriggers->Clone("fhNofMultipleTriggersGen"); + fOutputList->Add(fhNofMultipleTriggersGen); + + fhDeltaRMultTriggersGen = (TH1D*) fhDeltaRMultTriggers->Clone("fhDeltaRMultTriggersGen"); + fOutputList->Add(fhDeltaRMultTriggersGen); + } //------------------------------------- // pythia histograms @@ -1134,7 +1158,6 @@ void AliAnalysisTaskJetCorePP::UserExec(Option_t *) //============== analyze generator level MC ================ TList particleListGen; //list of tracks in MC - if(fIsChargedMC || fIsKine){ if(fIsKine){ //pure kine @@ -1255,6 +1278,25 @@ void AliAnalysisTaskJetCorePP::UserExec(Option_t *) if(ntriggersMC>0){ //there is at least one trigger Int_t rnd = fRandom->Integer(ntriggersMC); //0 to ntriggers-1 indexTriggGen = triggersMC[rnd]; + + fhNofMultipleTriggersGen->Fill(ntriggersMC-1); + if(ntriggersMC>1){ + Double_t deltaPhi, deltaEta, deltaR; + for(Int_t ia=0; iaPhi(),tGenII->Phi()); + deltaEta = tGenI->Eta()-tGenII->Eta(); + deltaR = sqrt(deltaPhi*deltaPhi + deltaEta*deltaEta); + fhDeltaRMultTriggersGen->Fill(deltaR); + } + } + } + }else{ indexTriggGen = -1; //trigger not found } @@ -1742,6 +1784,7 @@ Int_t AliAnalysisTaskJetCorePP::GetListOfTracks(TList *list){ Double_t ptmax = -10; Int_t triggers[200]; Int_t ntriggers = 0; //index in triggers array + for(Int_t it = 0; it < aodevt->GetNumberOfTracks(); it++){ AliAODTrack *tr = aodevt->GetTrack(it); @@ -1774,6 +1817,23 @@ Int_t AliAnalysisTaskJetCorePP::GetListOfTracks(TList *list){ if(fHardest==0 && ntriggers>0){ //select random inclusive trigger Int_t rnd = fRandom->Integer(ntriggers); //0 to ntriggers-1 index = triggers[rnd]; + + fhNofMultipleTriggers->Fill(ntriggers-1); + if(ntriggers>1){ + Double_t deltaPhi, deltaEta, deltaR; + for(Int_t ia=0; iaAt(ia); + if(!tGeni) continue; + for(Int_t ib=ia+1; ibAt(ib); + if(!tGenii) continue; + deltaPhi = RelativePhi(tGeni->Phi(),tGenii->Phi()); + deltaEta = tGeni->Eta()-tGenii->Eta(); + deltaR = sqrt(deltaPhi*deltaPhi + deltaEta*deltaEta); + fhDeltaRMultTriggers->Fill(deltaR); + } + } + } } return index; diff --git a/PWGJE/AliAnalysisTaskJetCorePP.h b/PWGJE/AliAnalysisTaskJetCorePP.h index b1bd2fae550..92e6751e540 100644 --- a/PWGJE/AliAnalysisTaskJetCorePP.h +++ b/PWGJE/AliAnalysisTaskJetCorePP.h @@ -171,6 +171,9 @@ private: TH1D *fhDphiTriggerJetAccept; //Deltaphi between trigger and jet after cut TH1D *fhCentrality; //Deltaphi between trigger and jet TH1D *fhCentralityAccept; //Deltaphi between trigger and jet after cut + TH1D *fhNofMultipleTriggers; // The number of additional triggers in events with at least one trigger + TH1D *fhDeltaRMultTriggers; // Angular distributions of mult. triggers. + //THnSparse *fHJetPtRaw; //bg unsubtr. vs bg subtr. pT spectrum of jets vs jet area //THnSparse *fHLeadingJetPtRaw; //bg unsubtr. vs bg. subtr. leading jet pT vs area @@ -204,6 +207,8 @@ private: TH1D *fhCellAreaToMedian; //how many entries were used to calculate TH1D *fhCellAreaToMedianGen; //how many entries were used to calculate in MC + TH1D *fhNofMultipleTriggersGen; // The number of additional triggers in events with at least one trigger + TH1D *fhDeltaRMultTriggersGen; // Angular distributions of mult. triggers. Bool_t fIsChargedMC; //flag analysis on MC data with true and on the real+kine data false Bool_t fIsKine; //flag analysis on kine data with true and on the real+MC data false -- 2.43.0