]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
improving cell time cuts in a way that MC productions also get the same cut.
authormcosenti <mcosenti@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Dec 2012 09:42:43 +0000 (09:42 +0000)
committermcosenti <mcosenti@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Dec 2012 09:42:43 +0000 (09:42 +0000)
PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.cxx
PWGGA/EMCALTasks/AliAnalysisTaskEMCALIsoPhoton.h

index d4e0bdf14e64eb3443a317a6d50358482d165734..72af76b39415f6ab119d371ac34578eefdce184e 100644 (file)
@@ -392,7 +392,7 @@ void AliAnalysisTaskEMCALIsoPhoton::FillClusHists()
     Float_t alliso, allphiband, allcore;
     Float_t phibandArea = (1.4 - 2*fIsoConeR)*2*fIsoConeR;
     Float_t netConeArea = TMath::Pi()*(fIsoConeR*fIsoConeR - 0.04*0.04);
-    GetCeIso(clsVec, ceiso, cephiband, cecore);
+    GetCeIso(clsVec, id, ceiso, cephiband, cecore);
     GetTrIso(clsVec, triso, trphiband, trcore);
     Double_t dr = TMath::Sqrt(c->GetTrackDx()*c->GetTrackDx() + c->GetTrackDz()*c->GetTrackDz());
     if(Et>10 && Et<15 && dr>0.025){
@@ -444,7 +444,7 @@ void AliAnalysisTaskEMCALIsoPhoton::FillClusHists()
 } 
 
 //________________________________________________________________________
-void AliAnalysisTaskEMCALIsoPhoton::GetCeIso(TVector3 vec, Float_t &iso, Float_t &phiband, Float_t &core)
+void AliAnalysisTaskEMCALIsoPhoton::GetCeIso(TVector3 vec, Int_t maxid, Float_t &iso, Float_t &phiband, Float_t &core)
 {
   // Get cell isolation.
 
@@ -457,6 +457,7 @@ void AliAnalysisTaskEMCALIsoPhoton::GetCeIso(TVector3 vec, Float_t &iso, Float_t
   Float_t etacl = vec.Eta();
   Float_t phicl = vec.Phi();
   Float_t thetacl = vec.Theta();
+  Float_t maxtcl = fEMCalCells->GetCellTime(maxid);
   if(phicl<0)
     phicl+=TMath::TwoPi();
   Int_t absid = -1;
@@ -464,7 +465,8 @@ void AliAnalysisTaskEMCALIsoPhoton::GetCeIso(TVector3 vec, Float_t &iso, Float_t
   for(int icell=0;icell<ncells;icell++){
     absid = TMath::Abs(fEMCalCells->GetCellNumber(icell));
     Float_t celltime = fEMCalCells->GetCellTime(absid);
-    if((celltime>2e-8 || celltime<-2e-8) && (!fIsMc))
+    //if(TMath::Abs(celltime)>2e-8 && (!fIsMc))
+    if(TMath::Abs(celltime-maxtcl)>2e-8 )
       continue;
     if(!fGeom)
       return;
index a737350201b07e8975ef9dc473c4c20d170736c4..a357ca292680ab3960099dd3996881300d637b7e 100644 (file)
@@ -31,7 +31,7 @@ class AliAnalysisTaskEMCALIsoPhoton : public AliAnalysisTaskSE {
   void                   UserExec(Option_t *option);
   void                   Terminate(Option_t *);
 
-  void                   GetCeIso(TVector3 vec, Float_t &iso, Float_t &phiband, Float_t &core);
+  void                   GetCeIso(TVector3 vec, Int_t maxid, Float_t &iso, Float_t &phiband, Float_t &core);
   Double_t               GetCrossEnergy(const AliVCluster *cluster, Short_t &idmax);
   Double_t               GetMaxCellEnergy(const AliVCluster *cluster, Short_t &id) const; 
   void                   GetTrIso(TVector3 vec, Float_t &iso, Float_t &phiband, Float_t &core);