Update Dpt Dpt and PID Corr: prabhat and debojit
authorsjena <sjena@cern.ch>
Fri, 7 Nov 2014 13:00:46 +0000 (14:00 +0100)
committersjena <sjena@cern.ch>
Fri, 7 Nov 2014 13:00:46 +0000 (14:00 +0100)
PWGCF/Correlations/DPhi/AliAnalysisTaskDptDptCorrelations.cxx
PWGCF/Correlations/DPhi/TriggerPID/AliTwoParticlePIDCorr.cxx
PWGCF/Correlations/macros/dptdptcorrelations/AddTaskDptDptCorr_dca1.C

index c03c7ba..e11f8d0 100644 (file)
@@ -908,7 +908,7 @@ void  AliAnalysisTaskDptDptCorrelations::createHistograms()
   name = "DCAz";    _dcaz     = createHisto1F(name,name, 500, -5.0, 5.0, "dcaZ","counts");
   name = "DCAxy";   _dcaxy    = createHisto1F(name,name, 500, -5.0, 5.0, "dcaXY","counts");
 
-  name = "Nclus1";   _Ncluster1    = createHisto1F(name,name, 200, 0, 200, "Ncluster1","counts");
+  // name = "Nclus1";   _Ncluster1    = createHisto1F(name,name, 200, 0, 200, "Ncluster1","counts");
   //name = "Nclus2";   _Ncluster2    = createHisto1F(name,name, 200, 0, 200, "Ncluster2","counts");
   
   if (_singlesOnly)
@@ -1153,7 +1153,7 @@ void  AliAnalysisTaskDptDptCorrelations::UserExec(Option_t */*option*/)
            }
        }
       
-      _vertexZ->Fill(vertexZ);
+      //_vertexZ->Fill(vertexZ);
       
       iVertex = int((vertexZ-_min_vertexZ)/_width_vertexZ);
       iVertexP1 = iVertex*_nBins_etaPhiPt_1;
@@ -1210,11 +1210,11 @@ void  AliAnalysisTaskDptDptCorrelations::UserExec(Option_t */*option*/)
          dedx   = t->GetTPCsignal();
          //dcaXY = t->DCA(); 
          //dcaZ  = t->ZAtDCA();  
-         nClus  = t->GetTPCNcls();       
+         //nClus  = t->GetTPCNcls();     
          
-          if ( nClus<_nClusterMin ) continue;
+         //if ( nClus<_nClusterMin ) continue;
          
-         _Ncluster1->Fill(nClus);
+         //_Ncluster1->Fill(nClus);
          
          /*
          //cuts on more than 0 shared cluster (suggested by Michael)
@@ -1266,7 +1266,7 @@ void  AliAnalysisTaskDptDptCorrelations::UserExec(Option_t */*option*/)
          //==== QA ===========================
          //_dcaz->Fill(DCAZ);
          //_dcaxy->Fill(DCAXY);
-         _etadis->Fill(eta);
+         //_etadis->Fill(eta);
          //_phidis->Fill(phi);
          //===================================
          //*************************************************
index 25a1855..4c30f13 100644 (file)
@@ -2502,6 +2502,40 @@ if (particletypeMC==SpProton)
   }
     }
 
+ //fill tracking efficiency
+ if(ffillefficiency)
+   {
+ if(particletypeMC==SpPion || particletypeMC==SpKaon)
+   {
+     if(TMath::Abs(pdgCode)==211 ||  TMath::Abs(pdgCode)==321) {
+       fTrackHistEfficiency[3]->Fill(allrecomatchedpid,4);//for mesons
+ if (((AliAODMCParticle*)recomatched)->IsPhysicalPrimary())  fTrackHistEfficiency[3]->Fill(allrecomatchedpid,3);//for mesons
+     }
+   }
+ if(particletypeMC==SpKaon || particletypeMC==SpProton)
+   {
+     if(TMath::Abs(pdgCode)==321 ||  TMath::Abs(pdgCode)==2212) {
+       fTrackHistEfficiency[4]->Fill(allrecomatchedpid,4);//for kaons+protons
+ if (((AliAODMCParticle*)recomatched)->IsPhysicalPrimary())  fTrackHistEfficiency[4]->Fill(allrecomatchedpid,3);
+     }
+   }
+ if(particletypeMC==SpPion && TMath::Abs(pdgCode)==211)  {
+   fTrackHistEfficiency[0]->Fill(allrecomatchedpid,4);//for pions
+ if (((AliAODMCParticle*)recomatched)->IsPhysicalPrimary())  fTrackHistEfficiency[0]->Fill(allrecomatchedpid,3);
+ } 
+ if(particletypeMC==SpKaon && TMath::Abs(pdgCode)==321) {
+   fTrackHistEfficiency[1]->Fill(allrecomatchedpid,4);//for kaons
+if (((AliAODMCParticle*)recomatched)->IsPhysicalPrimary())  fTrackHistEfficiency[1]->Fill(allrecomatchedpid,3);
+ }
+ if(particletypeMC==SpProton && TMath::Abs(pdgCode)==2212){
+   fTrackHistEfficiency[2]->Fill(allrecomatchedpid,4);//for protons
+if (((AliAODMCParticle*)recomatched)->IsPhysicalPrimary())  fTrackHistEfficiency[2]->Fill(allrecomatchedpid,3);
+ }
+   }
+
+ //remove the tracks which don't have proper TOF response-otherwise the misIDentification rate values will be wrong(do it after the efficiency filling is complete)
+if(fRequestTOFPID && track->Pt()>fPtTOFPIDmin && (!HasTOFPID(PIDtrack)) ) continue;
  //for misidentification fraction calculation(do it with tuneonPID)
  if(particletypeMC==SpPion )
    {
@@ -2524,7 +2558,7 @@ if(particletypeMC==SpKaon )
      if(TMath::Abs(pdgCode)==2212) fProtoncont->Fill(5.,track->Pt());
      if(TMath::Abs(pdgCode)!=211 && TMath::Abs(pdgCode)!=321 && TMath::Abs(pdgCode)!=2212) fProtoncont->Fill(7.,track->Pt());
    }
- if(particletypeMC==SpUndefined )
+ if(particletypeMC==SpUndefined )//these undefined are not due to absence of proper TOF response, rather due to the PID method only
    {
      if(TMath::Abs(pdgCode)==211) fUNIDcont->Fill(1.,track->Pt());
      if(TMath::Abs(pdgCode)==321) fUNIDcont->Fill(3.,track->Pt());
@@ -2538,37 +2572,6 @@ if(particletypeMC==SpKaon )
       FillPIDEventPlane(cent_v0,particletypeMC,track->Phi(),gReactionPlane);
     }
 
- //fill tracking efficiency
- if(ffillefficiency)
-   {
- if(particletypeMC==SpPion || particletypeMC==SpKaon)
-   {
-     if(TMath::Abs(pdgCode)==211 ||  TMath::Abs(pdgCode)==321) {
-       fTrackHistEfficiency[3]->Fill(allrecomatchedpid,4);//for mesons
- if (((AliAODMCParticle*)recomatched)->IsPhysicalPrimary())  fTrackHistEfficiency[3]->Fill(allrecomatchedpid,3);//for mesons
-     }
-   }
- if(particletypeMC==SpKaon || particletypeMC==SpProton)
-   {
-     if(TMath::Abs(pdgCode)==321 ||  TMath::Abs(pdgCode)==2212) {
-       fTrackHistEfficiency[4]->Fill(allrecomatchedpid,4);//for kaons+protons
- if (((AliAODMCParticle*)recomatched)->IsPhysicalPrimary())  fTrackHistEfficiency[4]->Fill(allrecomatchedpid,3);
-     }
-   }
- if(particletypeMC==SpPion && TMath::Abs(pdgCode)==211)  {
-   fTrackHistEfficiency[0]->Fill(allrecomatchedpid,4);//for pions
- if (((AliAODMCParticle*)recomatched)->IsPhysicalPrimary())  fTrackHistEfficiency[0]->Fill(allrecomatchedpid,3);
- } 
- if(particletypeMC==SpKaon && TMath::Abs(pdgCode)==321) {
-   fTrackHistEfficiency[1]->Fill(allrecomatchedpid,4);//for kaons
-if (((AliAODMCParticle*)recomatched)->IsPhysicalPrimary())  fTrackHistEfficiency[1]->Fill(allrecomatchedpid,3);
- }
- if(particletypeMC==SpProton && TMath::Abs(pdgCode)==2212){
-   fTrackHistEfficiency[2]->Fill(allrecomatchedpid,4);//for protons
-if (((AliAODMCParticle*)recomatched)->IsPhysicalPrimary())  fTrackHistEfficiency[2]->Fill(allrecomatchedpid,3);
- }
-   }
-
 if((track->Pt()>=fminPtAsso && track->Pt()<=fmaxPtAsso) || (track->Pt()>=fminPtTrig && track->Pt()<=fmaxPtTrig))//to reduce memory consumption in pool
   {
     Short_t chargeval=0;
@@ -4730,7 +4733,7 @@ Double_t AliTwoParticlePIDCorr::GetRefMultiOrCentrality(AliVEvent *mainevent, Bo
 
 if(fCentralityMethod=="V0M" || fCentralityMethod=="V0A" || fCentralityMethod=="V0C" || fCentralityMethod=="CL1" || fCentralityMethod=="ZNA" || fCentralityMethod=="V0AEq" || fCentralityMethod=="V0CEq" || fCentralityMethod=="V0MEq")//for PbPb, pPb, pp7TeV(still to be introduced)//data or RecoMC and also for TRUTH
     {
-      
+          
 if(fSampleType=="pp_7" && fPPVsMultUtils)
    {//for pp 7 TeV case only using Alianalysisutils class
        if(fAnalysisUtils) cent_v0 = fAnalysisUtils->GetMultiplicityPercentile((AliVEvent*)event,fCentralityMethod);
@@ -4742,7 +4745,7 @@ if(fSampleType=="pp_7" && fPPVsMultUtils)
   fHistCentStats->Fill(4.,fAnalysisUtils->GetMultiplicityPercentile((AliVEvent*)event,"V0CEq"));//only available for LHC10d at present (Quantile info)
   fHistCentStats->Fill(5.,fAnalysisUtils->GetMultiplicityPercentile((AliVEvent*)event,"V0MEq"));//only available for LHC10d at present (Quantile info)
       }
-      
+          
 else if(fSampleType=="pPb" || fSampleType=="PbPb")
   {
   AliCentrality *centralityObj=0;
index 5ed504e..f8654fc 100644 (file)
@@ -45,15 +45,11 @@ AliAnalysisTaskDptDptCorrelations *AddTaskDptDptCorr_dca1
     {
     if (centralityMethod == 4)
       {
-       minCentrality[0] = 0.0; maxCentrality[0] = 5.0;
-        minCentrality[1] = 5.0; maxCentrality[1] = 10.;
-        minCentrality[2] = 10.; maxCentrality[2] = 20.;
-        minCentrality[3] = 20.; maxCentrality[3] = 30.;
-        minCentrality[4] = 30.; maxCentrality[4] = 40.;
-        minCentrality[5] = 40.; maxCentrality[5] = 50.;
-        minCentrality[6] = 50.; maxCentrality[6] = 60.;
-        minCentrality[7] = 60.; maxCentrality[7] = 70.;
-       minCentrality[8] = 70.; maxCentrality[8] = 80.;
+       minCentrality[0] = 0.0;  maxCentrality[0]  = 5.0;
+        minCentrality[1] = 10.;  maxCentrality[1]  = 20.;
+        minCentrality[2] = 30.;  maxCentrality[2]  = 40.;
+        minCentrality[3] = 50.;  maxCentrality[3]  = 60.;
+        minCentrality[4] = 70.;  maxCentrality[4]  = 80.;
       }
     else
       {