Update from Prabhat DptDpt Corr
authorsjena <sjena@cern.ch>
Wed, 2 Jul 2014 07:45:27 +0000 (09:45 +0200)
committersjena <sjena@cern.ch>
Wed, 2 Jul 2014 07:45:27 +0000 (09:45 +0200)
PWGCF/Correlations/DPhi/AliAnalysisTaskDptDptCorrelations.cxx

index 17a287e..ca805f3 100644 (file)
@@ -1166,7 +1166,7 @@ void  AliAnalysisTaskDptDptCorrelations::UserExec(Option_t */*option*/)
       //====================== 
       
       //*********************************************************
-      TExMap *trackMap = new TExMap();//Mapping matrix----                                            
+       TExMap *trackMap = new TExMap();//Mapping matrix----                                            
 
       //1st loop track for Global tracks                                                                                
       for(Int_t i = 0; i < _nTracks; i++)
@@ -1190,10 +1190,10 @@ void  AliAnalysisTaskDptDptCorrelations::UserExec(Option_t */*option*/)
            AliError(Form("Could not receive track %d", iTrack));
            continue;
          }
-
+         
          bitOK  = t->TestFilterBit(_trackFilterBit);
          if (!bitOK) continue; //128bit or 272bit
-
+         
          Int_t gID = t->GetID();
          newAodTrack = gID >= 0 ?t : fAODEvent->GetTrack(trackMap->GetValue(-1-gID));
          
@@ -1209,7 +1209,7 @@ void  AliAnalysisTaskDptDptCorrelations::UserExec(Option_t */*option*/)
          //dcaZ  = t->ZAtDCA();  
          
          //for Global tracks
-         Double_t nsigmaelectron = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kElectron));
+          Double_t nsigmaelectron = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kElectron));
          Double_t nsigmapion = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kPion));
          Double_t nsigmakaon = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kKaon));
          Double_t nsigmaproton = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,(AliPID::EParticleType)AliPID::kProton));
@@ -1221,12 +1221,13 @@ void  AliAnalysisTaskDptDptCorrelations::UserExec(Option_t */*option*/)
             && nsigmakaon   > fNSigmaCut
             && nsigmaproton > fNSigmaCut ) continue;
          
+
          if(charge == 0) continue;
          // Kinematics cuts used                                                                                        
          if( pt < _min_pt_1 || pt > _max_pt_1) continue;
          if( eta < _min_eta_1 || eta > _max_eta_1) continue;
          
-         Double_t pos[3];
+         /*      Double_t pos[3];
          newAodTrack->GetXYZ(pos);
 
          Double_t DCAX = pos[0] - vertexX;
@@ -1238,14 +1239,14 @@ void  AliAnalysisTaskDptDptCorrelations::UserExec(Option_t */*option*/)
          if (DCAZ     <  _dcaZMin || 
              DCAZ     >  _dcaZMax ||
              DCAXY    >  _dcaXYMax ) continue; 
-
+         */
 
            //------- Eff. test---------- //just for checking
-         Double_t yy = (1 - 0.7)/1.8;
-         Double_t zz = (pt - 0.2);
-         Double_t effValue = 0.7 + yy*zz;
-         Double_t R = gRandom->Rndm();
-          if(R > effValue) continue;
+         //Double_t yy = (1 - 0.7)/1.8;
+         //Double_t zz = (pt - 0.2);
+         //Double_t effValue = 0.7 + yy*zz;
+         //Double_t R = gRandom->Rndm();
+          //if(R > effValue) continue;
          //---------------------------   
          
          //==== QA ===========================
@@ -1257,11 +1258,8 @@ void  AliAnalysisTaskDptDptCorrelations::UserExec(Option_t */*option*/)
          //*************************************************
                  
          //Particle 1
-         if (_requestedCharge_1 == charge &&
-             eta >= _min_eta_1 &&
-             eta < _max_eta_1) 
+         if (_requestedCharge_1 == charge)
            {
-             
              iPhi   = int( phi/_width_phi_1);
              
              if (iPhi<0 || iPhi>=_nBins_phi_1 ) 
@@ -1330,9 +1328,8 @@ void  AliAnalysisTaskDptDptCorrelations::UserExec(Option_t */*option*/)
                }
            }
          
-         if (!_sameFilter && _requestedCharge_2 == charge &&
-             eta >= _min_eta_2 &&
-             eta < _max_eta_2)  
+         if (!_sameFilter && _requestedCharge_2 == charge)
+              
            {
              
              iPhi   = int( phi/_width_phi_2);