remove use of EMC info from AliAODPID (M. Estienne)
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 3 May 2012 09:50:31 +0000 (09:50 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 3 May 2012 09:50:31 +0000 (09:50 +0000)
JETAN/AliJetAODFillUnitArrayEMCalDigits.cxx
JETAN/AliJetAODFillUnitArrayEMCalDigits.h
JETAN/AliJetAODFillUnitArrayTracks.cxx
JETAN/AliJetHadronCorrectionv1.cxx
JETAN/AliJetHadronCorrectionv1.h

index 43fb930..8efda07 100755 (executable)
@@ -186,8 +186,9 @@ void AliJetAODFillUnitArrayEMCalDigits::Exec(Option_t* const /*option*/)
 
       // Hadronic Correction
       double correction=0;
+// Temporarily commented - will be updated removing the AliAODpid dependence
+/*
       if (fApplyFractionHadronicCorrection) {
-      
         TArrayS* matched = new TArrayS();
        GetTracksPointingToCell(matched, etaD, phiD,0.02);
 
@@ -213,7 +214,7 @@ void AliJetAODFillUnitArrayEMCalDigits::Exec(Option_t* const /*option*/)
        delete matched;
        
       }//end hadronic correction
-
+*/
       double enerCorr = digitEn;
       if (correction >= enerCorr) enerCorr = 0;
       else enerCorr -= correction;
@@ -297,8 +298,9 @@ void AliJetAODFillUnitArrayEMCalDigits::Exec(Option_t* const /*option*/)
 
        // Hadronic Correction
        double correction=0;
+// Temporarily commented - will be updated removing the AliAODpid dependence
+/*
        if (fApplyFractionHadronicCorrection) {
-      
          TArrayS* matched = new TArrayS();
          GetTracksPointingToCell(matched, etaD, phiD,0.02);
          
@@ -324,7 +326,7 @@ void AliJetAODFillUnitArrayEMCalDigits::Exec(Option_t* const /*option*/)
          delete matched;
        
        }//end hadronic correction
-       
+*/     
        double enerCorr = digitEn;
        if (correction >= enerCorr) enerCorr = 0;
        else enerCorr -= correction;
@@ -383,42 +385,42 @@ void AliJetAODFillUnitArrayEMCalDigits::Exec(Option_t* const /*option*/)
     }
 }
 
-//____________________________________________________________________________
-void AliJetAODFillUnitArrayEMCalDigits::GetTracksPointingToCell(TArrayS* array,Double_t eta, Double_t phi, Double_t cut)
-{
-// Get all tracks pointing to cell 
-  int size=0;
-  
-  for (Int_t itrk =  0; itrk <  fAOD->GetNumberOfTracks() ; itrk++) { //track loop
-  
-    AliAODTrack * track = (AliAODTrack*) fAOD->GetTrack(itrk) ;  
-    AliAODPid*    pid   = (AliAODPid*) track->GetDetPid();
-    
-    if(pid) {
-      Double_t emcpos[3];
-      pid->GetEMCALPosition(emcpos);      
-      TVector3 tpos(emcpos[0],emcpos[1],emcpos[2]);
-
-      Double_t deta = tpos.Eta() - eta;
-      Double_t dphi = tpos.Phi() - phi;
-
-      if(dphi > TMath::Pi()) dphi -= 2*TMath::Pi();
-      if(dphi < -TMath::Pi()) dphi += 2*TMath::Pi();
-
-      Double_t res = sqrt(dphi*dphi + deta*deta);
-      
-      if (res< cut) {
-        //add this track-index
-        size++;
-        array->Set( size );
-        array->AddAt( itrk, (size-1) ); 
-       if(fDebug>1) printf("MTH:: track %d matched cell at eta=%f , phi=%f \n", itrk, eta, phi);
-                     
-      }
-    }
-  }
-
-}
+////____________________________________________________________________________
+//void AliJetAODFillUnitArrayEMCalDigits::GetTracksPointingToCell(TArrayS* array,Double_t eta, Double_t phi, Double_t cut)
+//{ // Temporarily commented -> will be corrected removing the dependence to AliAODPid
+//// Get all tracks pointing to cell 
+//  int size=0;
+//  
+//  for (Int_t itrk =  0; itrk <  fAOD->GetNumberOfTracks() ; itrk++) { //track loop
+//  
+//    AliAODTrack * track = (AliAODTrack*) fAOD->GetTrack(itrk) ;  
+//    AliAODPid*    pid   = (AliAODPid*) track->GetDetPid();
+//    
+//    if(pid) {
+//      Double_t emcpos[3];
+//      pid->GetEMCALPosition(emcpos);      
+//      TVector3 tpos(emcpos[0],emcpos[1],emcpos[2]);
+//
+//      Double_t deta = tpos.Eta() - eta;
+//      Double_t dphi = tpos.Phi() - phi;
+//
+//      if(dphi > TMath::Pi()) dphi -= 2*TMath::Pi();
+//      if(dphi < -TMath::Pi()) dphi += 2*TMath::Pi();
+//
+//      Double_t res = sqrt(dphi*dphi + deta*deta);
+//      
+//      if (res< cut) {
+//        //add this track-index
+//        size++;
+//        array->Set( size );
+//        array->AddAt( itrk, (size-1) ); 
+//     if(fDebug>1) printf("MTH:: track %d matched cell at eta=%f , phi=%f \n", itrk, eta, phi);
+//                   
+//      }
+//    }
+//  }
+//
+//}
 
 /*
 //____________________________________________________________________________
index f529a40..4277abc 100755 (executable)
@@ -54,7 +54,8 @@ class AliJetAODFillUnitArrayEMCalDigits : public AliJetFillUnitArray
   Double_t     fFractionHadronicCorrection;       // Fraction hadronic correction 
 
   //Track-matching (mth)
-  void GetTracksPointingToCell(TArrayS *arr, Double_t eta, Double_t phi, Double_t res);
+// Temporarily commented - will be updated removing the AliAODpid dependence
+//  void GetTracksPointingToCell(TArrayS *arr, Double_t eta, Double_t phi, Double_t res);
 
   // geometry info
   AliAODCaloCluster *fClus;           //! 
index 91fe93a..8628a29 100644 (file)
@@ -516,8 +516,10 @@ void AliJetAODFillUnitArrayTracks::Exec(Option_t* const /*option*/)
              // For the moment I apply MIP correction if p >= 0.5 GeV/c
              // and only for tracks inside EMCal acceptance
              // End of if(fGrid==1) -> hadron correction for all tracks
+// Temporarily commented - will be updated removing the AliAODpid dependence
+/*
              if (fApplyMIPCorrection != 0 && p3.Mag() >= 0.5) 
-               { 
+               {
                  ((AliJetHadronCorrectionv1*)fHadCorr)->SetGeometry("EMCAL_COMPLETEV1",1.);
 
                  Double_t etaOut = 0.;
@@ -533,6 +535,7 @@ void AliJetAODFillUnitArrayTracks::Exec(Option_t* const /*option*/)
                      unitEnergy -= hCEnergy*TMath::Sin(2.0*TMath::ATan(TMath::Exp(-eta)));
                    } 
                } //end Hadron Correction loop
+*/
 
              if((unitEnergy + pt) > 0.) uArray->SetUnitEnergy(unitEnergy + pt);
              else uArray->SetUnitEnergy(0.);
index c0f8874..a67f7eb 100644 (file)
@@ -168,19 +168,19 @@ Double_t AliJetHadronCorrectionv1::GetEnergy(Double_t pmom, Double_t eta, Int_t
    
 }
 
-void AliJetHadronCorrectionv1::TrackPositionEMCal(const AliAODTrack* track,Double_t &eta, Double_t &phi)
-{
-// Return track position on EMCal
-  AliAODPid*    pid   = (AliAODPid*) track->GetDetPid();
-    
-  if(pid) {
-    Double_t emcpos[3];
-    pid->GetEMCALPosition(emcpos);      
-    TVector3 tpos(emcpos[0],emcpos[1],emcpos[2]);
-    
-    eta = tpos.Eta();
-    phi = tpos.Phi();
-
-  }
-
-}
+// void AliJetHadronCorrectionv1::TrackPositionEMCal(const AliAODTrack* track,Double_t &eta, Double_t &phi)
+// { // Temporarily commented - will be updated removing the AliAODpid dependence
+// // Return track position on EMCal
+//   AliAODPid*    pid   = (AliAODPid*) track->GetDetPid();
+//    
+//   if(pid) {
+//     Double_t emcpos[3];
+//     pid->GetEMCALPosition(emcpos);      
+//     TVector3 tpos(emcpos[0],emcpos[1],emcpos[2]);
+//     
+//     eta = tpos.Eta();
+//     phi = tpos.Phi();
+// 
+//   }
+// 
+// }
index b210685..c4c6002 100644 (file)
@@ -24,7 +24,8 @@ class AliJetHadronCorrectionv1 : public AliJetHadronCorrection
   
   void SetGeometry(TString name, Double_t fs = 1.); 
   void SetGeometry2(const AliJetDummyGeo *geometry);
-  void TrackPositionEMCal(const AliAODTrack* track,Double_t &eta, Double_t &phi);
+// Temporarily commented - will be updated removing the AliAODpid dependence
+//  void TrackPositionEMCal(const AliAODTrack* track,Double_t &eta, Double_t &phi);
 
  protected:
   AliJetHadronCorrectionv1():fSamplingFraction(0) {for (Int_t i = 0; i < 8; i++) fPar[i] = 0.;}