including isolation using cluster Et instead of the core
authormcosenti <mcosenti@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 19 Jul 2012 10:11:13 +0000 (10:11 +0000)
committermcosenti <mcosenti@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 19 Jul 2012 10:11:13 +0000 (10:11 +0000)
PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx

index b6e8eb3..cba26c6 100644 (file)
@@ -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();