come back origianl match
authorssakai <Shingo.Sakai@lnf.infn.it>
Sat, 2 Aug 2014 06:58:07 +0000 (08:58 +0200)
committerssakai <Shingo.Sakai@lnf.infn.it>
Tue, 5 Aug 2014 12:24:51 +0000 (14:24 +0200)
PWGHF/hfe/AliAnalysisTaskHFECal.cxx
PWGHF/hfe/AliAnalysisTaskHFECal.h

index 8d2d265..af6603d 100644 (file)
@@ -824,9 +824,7 @@ void AliAnalysisTaskHFECal::UserExec(Option_t*)
         double emctof = 0.0;
         Bool_t MaxEmatch = kFALSE;
 
-    Int_t clsId_ESD = track->GetEMCALcluster();
-    Int_t clsId = -9999;
-    if(phi>1.0 && phi<3.5 && fabs(eta)<0.75)clsId = TrMatch(track);
+    Int_t clsId = track->GetEMCALcluster();
 
     //cout << "match ID" << clsId << " ; " << clsId_Tender << endl;
 
@@ -840,15 +838,6 @@ void AliAnalysisTaskHFECal::UserExec(Option_t*)
                 eop = clustE/fabs(mom);
 
 
-                   MatchTrCheck->Fill(clsId_ESD,clsId);
-                   double eop_ESD = -10.0;
-                   if(clsId_ESD>=0)
-                     {
-                       AliESDCaloCluster *clust_ESD = fESD->GetCaloCluster(clsId_ESD);
-                       eop_ESD = clust_ESD->E()/fabs(mom);
-                     }
-                   if(fTPCnSigma>-1 && fTPCnSigma<3 && pt>3.0)MatchTrEop->Fill(eop_ESD,eop);
-
                 //double clustT = clust->GetTOF();
                 ncells = clust->GetNCells();
                 m02 = clust->GetM02();
@@ -2314,65 +2303,3 @@ double AliAnalysisTaskHFECal::NsigmaCorrection(double tmpeta, float central)
 
 }
 
-//________ matching
-int AliAnalysisTaskHFECal::TrMatch(AliESDtrack *tr)
-{
-  
-  Float_t  clsPos[3];
-  Double_t trkPos[3];
-  int matchID = -1;
-  Double_t MaxR = 0.05;
-  Double_t MaxEta = 0.025;
-  Double_t MaxPhi = 0.05;
-
-  Double_t  magF = fESD->GetMagneticField();
-  Double_t magSign = 1.0;
-  if(magF<0)magSign = -1.0;
-  //printf("magF ; %g ; %g \n", magF,magSign);
-  if (!TGeoGlobalMagField::Instance()->GetField()) {
-         printf("Loading field map...\n");
-         //AliMagF* field = new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG);
-         AliMagF* field = new AliMagF("Maps","Maps", magSign, magSign, AliMagF::k5kG); // for 10d
-         TGeoGlobalMagField::Instance()->SetField(field);
-                   }
-
-  AliEMCALTrack *emctrack = new AliEMCALTrack(*tr);
-  Double_t fieldB[3]; 
-  emctrack->GetBxByBz(fieldB);
-  //printf("%g %g %g \n", fieldB[0], fieldB[1], fieldB[2]);
-
-  for(Int_t icl=0; icl<fESD->GetNumberOfCaloClusters(); icl++)
-   {
-
-   AliVCluster *cluster = (AliVCluster*) fESD->GetCaloCluster(icl);
-   if(!cluster->IsEMCAL()) continue;
-
-   cluster->GetPosition(clsPos);
-   if(!emctrack->PropagateToGlobal(clsPos[0],clsPos[1],clsPos[2],0.,0.) )  continue;
-   emctrack->GetXYZ(trkPos);
-
-   TVector3 clsPosVec(clsPos[0],clsPos[1],clsPos[2]);
-   TVector3 trkPosVec(trkPos[0],trkPos[1],trkPos[2]);
-
-   Double_t delEmcphi = clsPosVec.Phi()-trkPosVec.Phi();  // track cluster matching
-   Double_t delEmceta = clsPosVec.Eta()-trkPosVec.Eta();  // track cluster matching
-
-   if(fabs(delEmcphi)>0.05)continue;
-   if(fabs(delEmceta)>0.025)continue;
-
-   double rmatch = sqrt(pow(delEmcphi,2)+pow(delEmceta,2));
-   
-   //if(rmatch<MaxR)
-   if(fabs(delEmcphi)<MaxPhi && fabs(delEmceta)<MaxEta )
-     {
-      matchID = icl;
-      //MaxR = rmatch;
-      MaxPhi = fabs(delEmcphi);
-      MaxEta = fabs(delEmceta);
-     }
-  }
- return matchID;
-}
-
index cd9d58c..447f217 100644 (file)
@@ -68,7 +68,6 @@ class AliAnalysisTaskHFECal : public AliAnalysisTaskSE {
   void FindTriggerClusters();
   double MCEopMeanCorrection(double pTmc, float central);
   double NsigmaCorrection(double tmpeta, float central);
-  int TrMatch(AliESDtrack *tr);
  private:
   
   Bool_t ProcessCutStep(Int_t cutStep, AliVParticle *track);