From: cvetan Date: Mon, 18 Feb 2013 14:45:03 +0000 (+0000) Subject: Adding histos with muon eta. X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=b897e6f84ee4a4bd8df66a33fa7e05e6424651a8 Adding histos with muon eta. --- diff --git a/PWGCF/Correlations/DPhi/MuonHadron/AddAnalysisTaskMuonHadronCorrelations.C b/PWGCF/Correlations/DPhi/MuonHadron/AddAnalysisTaskMuonHadronCorrelations.C index 9ecad21d3f4..249b68a7b09 100644 --- a/PWGCF/Correlations/DPhi/MuonHadron/AddAnalysisTaskMuonHadronCorrelations.C +++ b/PWGCF/Correlations/DPhi/MuonHadron/AddAnalysisTaskMuonHadronCorrelations.C @@ -23,6 +23,10 @@ AliAnalysisTaskMuonHadronCorrelations *AddAnalysisTaskMuonHadronCorrelations(con Double_t ptLimits[nBinPt+1] = {0., 1., 2., 4.}; task->SetPtBinning(nBinPt, ptLimits); + const Int_t nBinEta = 3; + Double_t etaLimits[nBinEta+1] = {-4., -3.6, -3.2, -2.5}; + task->SetEtaBinning(nBinEta, etaLimits); + mgr->AddTask(task); // create output container diff --git a/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskMuonHadronCorrelations.cxx b/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskMuonHadronCorrelations.cxx index 8881b5cbea7..869eddd74c5 100644 --- a/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskMuonHadronCorrelations.cxx +++ b/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskMuonHadronCorrelations.cxx @@ -33,6 +33,7 @@ AliAnalysisTaskMuonHadronCorrelations::AliAnalysisTaskMuonHadronCorrelations() : fNbinsPt(1), fCentAxis(0x0), fPtAxis(0x0), + fEtaAxis(0x0), fHistV0Multiplicity(0x0), fHistITSMultiplicity(0x0), fHistCentrality(0x0), @@ -48,6 +49,8 @@ AliAnalysisTaskMuonHadronCorrelations::AliAnalysisTaskMuonHadronCorrelations() : for (Int_t iPtBinMA=0; iPtBinMAIsProofMode()) delete fOutputList; @@ -164,6 +175,39 @@ void AliAnalysisTaskMuonHadronCorrelations::UserCreateOutputObjects() { fOutputList -> Add(fHistDeltaPhi[iCent][iPtBinCB][iPtBinMA]); fOutputList -> Add(fHistDeltaPhiMix[iCent][iPtBinCB][iPtBinMA]); + fHistEtaDeltaPhi[iCent][iPtBinCB][iPtBinMA] = new TH2D(Form("fHistEtaDeltaPhi_Cent%02d_PtBin%02d_%02d",iCent,iPtBinCB,iPtBinMA), + Form("%d-%d %%, %3.1fGetBinLowEdge(iCent+1)), + Int_t(fCentAxis->GetBinUpEdge(iCent+1)), + fPtAxis->GetBinLowEdge(iPtBinCB+1), + fPtAxis->GetBinUpEdge(iPtBinCB+1), + fPtAxis->GetBinLowEdge(iPtBinMA+1), + fPtAxis->GetBinUpEdge(iPtBinMA+1)), + 100, -0.5*TMath::RadToDeg()*TMath::Pi(), 1.5*TMath::RadToDeg()*TMath::Pi(), + fEtaAxis->GetNbins(),(Double_t*)fEtaAxis->GetXbins()->GetArray()); + + fHistEtaDeltaPhiMix[iCent][iPtBinCB][iPtBinMA] = new TH2D(Form("fHistEtaDeltaPhiMix_Cent%02d_PtBin%02d_%02d",iCent,iPtBinCB,iPtBinMA), + Form("%d-%d %%, %3.1fGetBinLowEdge(iCent+1)), + Int_t(fCentAxis->GetBinUpEdge(iCent+1)), + fPtAxis->GetBinLowEdge(iPtBinCB+1), + fPtAxis->GetBinUpEdge(iPtBinCB+1), + fPtAxis->GetBinLowEdge(iPtBinMA+1), + fPtAxis->GetBinUpEdge(iPtBinMA+1)), + 100, -0.5*TMath::RadToDeg()*TMath::Pi(), 1.5*TMath::RadToDeg()*TMath::Pi(), + fEtaAxis->GetNbins(),(Double_t*)fEtaAxis->GetXbins()->GetArray()); + + fHistEtaDeltaPhi[iCent][iPtBinCB][iPtBinMA] -> SetXTitle("#Delta#varphi [degrees]"); + fHistEtaDeltaPhiMix[iCent][iPtBinCB][iPtBinMA] -> SetXTitle("#Delta#varphi [degrees]"); + fHistEtaDeltaPhi[iCent][iPtBinCB][iPtBinMA] -> SetYTitle("#eta"); + fHistEtaDeltaPhiMix[iCent][iPtBinCB][iPtBinMA] -> SetYTitle("#eta"); + + fHistEtaDeltaPhi[iCent][iPtBinCB][iPtBinMA] -> Sumw2(); + fHistEtaDeltaPhiMix[iCent][iPtBinCB][iPtBinMA] -> Sumw2(); + + fOutputList -> Add(fHistEtaDeltaPhi[iCent][iPtBinCB][iPtBinMA]); + fOutputList -> Add(fHistEtaDeltaPhiMix[iCent][iPtBinCB][iPtBinMA]); + } } @@ -202,6 +246,17 @@ void AliAnalysisTaskMuonHadronCorrelations::UserCreateOutputObjects() { fOutputList -> Add(fHistSingleMuonsPt[iCent]); fOutputList -> Add(fHistSingleMuonsPtmixed[iCent]); + fHistSingleMuonsEtaPt[iCent] = new TH2D(Form("fHistSingleMuonEtaPt_Cent%02d",iCent), + "#eta vs p_{T} for single muons", + fNbinsPt, (Double_t*)fPtAxis->GetXbins()->GetArray(), + fEtaAxis->GetNbins(),(Double_t*)fEtaAxis->GetXbins()->GetArray()); + fHistSingleMuonsEtaPtmixed[iCent] = new TH2D(Form("fHistSingleMuonEtaPtmixed_Cent%02d",iCent), + "#eta vs p_{T} for single muons", + fNbinsPt, (Double_t*)fPtAxis->GetXbins()->GetArray(), + fEtaAxis->GetNbins(),(Double_t*)fEtaAxis->GetXbins()->GetArray()); + fOutputList -> Add(fHistSingleMuonsEtaPt[iCent]); + fOutputList -> Add(fHistSingleMuonsEtaPtmixed[iCent]); + fHistSingleMuonsEtaVsPt[iCent] = new TH2D(Form("fHistSingleMuonsEtaVsPt_Cent%02d",iCent), "#eta vs p_{T} for single muons", 100, -4.5, -2., 100, 0., 10.); @@ -303,6 +358,7 @@ void AliAnalysisTaskMuonHadronCorrelations::UserExec(Option_t *) { for (Int_t iTrMA=0; iTrMAGetEntriesFast(); iTrMA++) { fTrackMA = (AliAODTrack*) tracksMuonArm->At(iTrMA); fHistSingleMuonsPt[centBin]->Fill(fTrackMA->Pt()); + fHistSingleMuonsEtaPt[centBin]->Fill(fTrackMA->Pt(),fTrackMA->Eta()); for (Int_t iTrCB=0; iTrCBGetEntriesFast(); iTrCB++) { fTrackCB = (AliAODTrack*) tracksCentralBarrel -> At(iTrCB); FillHistograms(centBin, kSingleEvent); @@ -320,6 +376,7 @@ void AliAnalysisTaskMuonHadronCorrelations::UserExec(Option_t *) { for (Int_t iTrMA=0; iTrMAGetEntriesFast(); iTrMA++) { fTrackMA = (AliAODTrack*) tracksMuonArm->At(iTrMA); fHistSingleMuonsPtmixed[centBin]->Fill(fTrackMA->Pt()); + fHistSingleMuonsEtaPtmixed[centBin]->Fill(fTrackMA->Pt(),fTrackMA->Eta()); for (Int_t iTrCB=0; iTrCBGetEntriesFast(); iTrCB++) { fTrackCB = (AliAODTrack*) mixedTracks -> At(iTrCB); FillHistograms(centBin, kMixedEvent); @@ -351,6 +408,9 @@ void AliAnalysisTaskMuonHadronCorrelations::FillHistograms(Int_t centrality, Int if (option==kSingleEvent) fHistDeltaPhi[centrality][ptBinTrackCB-1][ptBinTrackMA-1] -> Fill(TMath::RadToDeg()*deltaPhi); else if (option==kMixedEvent) fHistDeltaPhiMix[centrality][ptBinTrackCB-1][ptBinTrackMA-1] -> Fill(TMath::RadToDeg()*deltaPhi); + if (option==kSingleEvent) fHistEtaDeltaPhi[centrality][ptBinTrackCB-1][ptBinTrackMA-1] -> Fill(TMath::RadToDeg()*deltaPhi,fTrackMA->Eta()); + else if (option==kMixedEvent) fHistEtaDeltaPhiMix[centrality][ptBinTrackCB-1][ptBinTrackMA-1] -> Fill(TMath::RadToDeg()*deltaPhi,fTrackMA->Eta()); + if (option==kSingleEvent) fHistTracksEtaMAvsEtaCB[centrality]->Fill(fTrackMA->Eta(),fTrackCB->Eta()); else if (option==kMixedEvent) fHistTracksEtaMAvsEtaCBmixed[centrality]->Fill(fTrackMA->Eta(),fTrackCB->Eta()); } @@ -470,6 +530,23 @@ void AliAnalysisTaskMuonHadronCorrelations::SetPtBinning(Int_t nBins, Double_t * //==================================================================================================================================================== +void AliAnalysisTaskMuonHadronCorrelations::SetEtaBinning(Int_t nBins, Double_t *limits) { + + if (nBins>fNMaxBinsEta) { + AliInfo(Form("WARNING : only %d pt bins (out of the %d proposed) will be considered",fNMaxBinsEta,nBins)); + nBins = fNMaxBinsEta; + } + if (nBins<=0) { + AliInfo("WARNING : at least one pt bin must be considered"); + nBins = 1; + } + + fEtaAxis = new TAxis(nBins, limits); + +} + +//==================================================================================================================================================== + Int_t AliAnalysisTaskMuonHadronCorrelations::GetCentBin() { Double_t percentile = fAOD->GetCentrality()->GetCentralityPercentile(fCentMethod.Data()); diff --git a/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskMuonHadronCorrelations.h b/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskMuonHadronCorrelations.h index 00d9ffd6257..4f6cd737bc5 100644 --- a/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskMuonHadronCorrelations.h +++ b/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskMuonHadronCorrelations.h @@ -45,6 +45,7 @@ class AliAnalysisTaskMuonHadronCorrelations : public AliAnalysisTaskSE { TObjArray* GetAcceptedTracksMuonArm(AliAODEvent *aodEvent, Int_t centBin); void SetPtBinning(Int_t nBins, Double_t *limits); void SetCentBinning(Int_t nBins, Double_t *limits); + void SetEtaBinning(Int_t nBins, Double_t *limits); void SetCentMethod(const Char_t *method) { fCentMethod = method; } void FillHistograms(Int_t centrality, Int_t option); Int_t GetCentBin(); @@ -53,6 +54,7 @@ class AliAnalysisTaskMuonHadronCorrelations : public AliAnalysisTaskSE { static const Int_t fNMaxBinsCentrality = 20; static const Int_t fNMaxBinsPt = 10; + static const Int_t fNMaxBinsEta = 10; AliAODEvent *fAOD; //! AliEventPoolManager *fPoolMgr; //! event pool manager @@ -69,15 +71,20 @@ class AliAnalysisTaskMuonHadronCorrelations : public AliAnalysisTaskSE { TAxis *fCentAxis; TAxis *fPtAxis; + TAxis *fEtaAxis; TH1D *fHistDeltaPhi[fNMaxBinsCentrality][fNMaxBinsPt][fNMaxBinsPt]; //! TH1D *fHistDeltaPhiMix[fNMaxBinsCentrality][fNMaxBinsPt][fNMaxBinsPt]; //! + TH2D *fHistEtaDeltaPhi[fNMaxBinsCentrality][fNMaxBinsPt][fNMaxBinsPt]; //! + TH2D *fHistEtaDeltaPhiMix[fNMaxBinsCentrality][fNMaxBinsPt][fNMaxBinsPt]; //! TH2D *fHistNTracksCB_vs_NTracksMA[fNMaxBinsCentrality]; //! TH2D *fHistNTracksCB_vs_NTracksMAmixed[fNMaxBinsCentrality]; //! TH2D *fHistTracksEtaMAvsEtaCB[fNMaxBinsCentrality]; //! TH2D *fHistTracksEtaMAvsEtaCBmixed[fNMaxBinsCentrality]; //! TH1D *fHistSingleMuonsPt[fNMaxBinsCentrality]; //! TH1D *fHistSingleMuonsPtmixed[fNMaxBinsCentrality]; //! + TH2D *fHistSingleMuonsEtaPt[fNMaxBinsCentrality]; //! + TH2D *fHistSingleMuonsEtaPtmixed[fNMaxBinsCentrality]; //! TH1D *fHistV0Multiplicity; //! TH1D *fHistITSMultiplicity; //! @@ -94,7 +101,7 @@ class AliAnalysisTaskMuonHadronCorrelations : public AliAnalysisTaskSE { AliAnalysisTaskMuonHadronCorrelations(const AliAnalysisTaskMuonHadronCorrelations&);//not implimented AliAnalysisTaskMuonHadronCorrelations& operator=(const AliAnalysisTaskMuonHadronCorrelations&);//not implimnted - ClassDef(AliAnalysisTaskMuonHadronCorrelations, 1) // example of analysis + ClassDef(AliAnalysisTaskMuonHadronCorrelations, 2) // example of analysis };