From 717bb45bd32612d71408843e2295cf11321f536e Mon Sep 17 00:00:00 2001 From: mcosenti Date: Thu, 19 Jul 2012 10:11:13 +0000 Subject: [PATCH] including isolation using cluster Et instead of the core --- .../AliAnalysisTaskEMCALIsoPhoton.cxx | 20 ++++++++++--------- 1 file changed, 11 insertions(+), 9 deletions(-) diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx index b6e8eb300a7..cba26c6a6dc 100644 --- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx +++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx @@ -124,13 +124,13 @@ void AliAnalysisTaskEMCALIsoPhoton::UserCreateOutputObjects() fNClusHighClusE = new TH2F("hNClusHighClusE","total number of clusters vs. highest clus energy in the event;E (GeV);NClus",200,0,100,301,-0.5,300.5); fOutputList->Add(fNClusHighClusE); - Int_t nEt=1000, nM02=400, nCeIso=1000, nTrIso=1000, nAllIso=1000, nTrClDphi=200, nTrClDeta=100, nClEta=140, nClPhi=128; - Int_t bins[] = {nEt, nM02, nCeIso, nTrIso, nAllIso, nTrClDphi, nTrClDeta,nClEta,nClPhi}; + Int_t nEt=1000, nM02=400, nCeIso=1000, nTrIso=1000, nAllIso=1000, nCeIsoNC=1000, nAllIsoNC=1000, nTrClDphi=200, nTrClDeta=100, nClEta=140, nClPhi=128; + Int_t bins[] = {nEt, nM02, nCeIso, nTrIso, nAllIso, nCeIsoNC, nAllIsoNC, nTrClDphi, nTrClDeta,nClEta,nClPhi}; fNDimensions = sizeof(bins)/sizeof(Int_t); const Int_t ndims = fNDimensions; - Double_t xmin[] = { 0., 0., -10., -10., -10., -0.1,-0.05, -0.7, 1.4}; - Double_t xmax[] = { 100., 4., 190., 190., 190., 0.1, 0.05, 0.7, 3.192}; - fHnOutput = new THnSparseF("fHnOutput","Output matrix: E_{T},M02,CeIso,TrIso,AllIso, d#phi_{trk},d#eta_{trk},#eta_{clus},#phi_{clus}", ndims, bins, xmin, xmax); + Double_t xmin[] = { 0., 0., -10., -10., -10., -10., -10., -0.1,-0.05, -0.7, 1.4}; + Double_t xmax[] = { 100., 4., 190., 190., 190., 190., 190., 0.1, 0.05, 0.7, 3.192}; + fHnOutput = new THnSparseF("fHnOutput","Output matrix: E_{T},M02,CeIso,TrIso,AllIso, CeIsoNoCore, AllIsoNoCore, d#phi_{trk},d#eta_{trk},#eta_{clus},#phi_{clus}", ndims, bins, xmin, xmax); fOutputList->Add(fHnOutput); @@ -260,10 +260,12 @@ void AliAnalysisTaskEMCALIsoPhoton::FillClusHists() outputValues[2] = ceiso-cecore-ceisoue; outputValues[3] = triso-trisoue; outputValues[4] = alliso-cecore-allisoue; - outputValues[5] = c->GetTrackDx(); - outputValues[6] = c->GetTrackDz(); - outputValues[7] = clsVec.Eta(); - outputValues[8] = clsVec.Phi(); + outputValues[5] = ceiso-Et-ceisoue; + outputValues[6] = alliso-Et-allisoue; + outputValues[7] = c->GetTrackDx(); + outputValues[8] = c->GetTrackDz(); + outputValues[9] = clsVec.Eta(); + outputValues[10] = clsVec.Phi(); fHnOutput->Fill(outputValues); if(c->E()>maxE) maxE = c->E(); -- 2.43.0