]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGLF/STRANGENESS/Hypernuclei/AliAnalysisTaskNucleiv2SP.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
Changed TOF PID
[u/mrichter/AliRoot.git] / PWGLF / STRANGENESS / Hypernuclei / AliAnalysisTaskNucleiv2SP.cxx
diff --git a/PWGLF/STRANGENESS/Hypernuclei/AliAnalysisTaskNucleiv2SP.cxx b/PWGLF/STRANGENESS/Hypernuclei/AliAnalysisTaskNucleiv2SP.cxx
index f0bac1db5633c70ad9cc4a8afda2decd0445b495..54a3fee8d5da523e3cc2fb4eabebf3a894425f07 100644 (file)
--- a/PWGLF/STRANGENESS/Hypernuclei/AliAnalysisTaskNucleiv2SP.cxx
+++ b/PWGLF/STRANGENESS/Hypernuclei/AliAnalysisTaskNucleiv2SP.cxx
@@ -98,6 +98,12 @@ AliAnalysisTaskNucleiv2SP::AliAnalysisTaskNucleiv2SP()
   hCos2DeltaVzMTPCnvsCentrality(0),
   hCos2DeltaTPCpTPCnvsCentrality(0),
   hQVzAQVzCvsCentrality(0),
+  hQxVzAvsCentrality(0),
+  hQyVzAvsCentrality(0),
+  hQxVzCvsCentrality(0),
+  hQyVzCvsCentrality(0),
+  hQxVzMvsCentrality(0),
+  hQyVzMvsCentrality(0),
   ftree(0),           
   tCentrality(0),     
   tType(0),  
@@ -165,6 +171,12 @@ AliAnalysisTaskNucleiv2SP::AliAnalysisTaskNucleiv2SP(const char *name)
   hCos2DeltaVzMTPCnvsCentrality(0),
   hCos2DeltaTPCpTPCnvsCentrality(0),
   hQVzAQVzCvsCentrality(0),
+  hQxVzAvsCentrality(0),
+  hQyVzAvsCentrality(0),
+  hQxVzCvsCentrality(0),
+  hQyVzCvsCentrality(0),
+  hQxVzMvsCentrality(0),
+  hQyVzMvsCentrality(0),
   ftree(0),           
   tCentrality(0),     
   tType(0),  
@@ -311,17 +323,17 @@ void AliAnalysisTaskNucleiv2SP::UserCreateOutputObjects()
   } 
  
   if(! fhBB ){
-    fhBB = new TH2F( "fhBB" , "BetheBlochTPC" , 240,-6,6,250,0,1000);
+    fhBB = new TH2F( "fhBB" , "BetheBlochTPC" , 240,-10,10,250,0,1000);
     fListHist->Add(fhBB);
   }
   
   if(! fhBBDeu ){
-    fhBBDeu = new TH2F( "fhBBDeu" , "BetheBlochTPC - Deuteron" , 240,-6,6,250,0,1000);
+    fhBBDeu = new TH2F( "fhBBDeu" , "BetheBlochTPC - Deuteron" , 240,-10,10,250,0,1000);
     fListHist->Add(fhBBDeu);
   }
  
   if(! fhTOF ){
-    fhTOF = new TH2F( "fhTOF" , "Scatter Plot TOF" , 240,-6,6,500,0,1.2);
+    fhTOF = new TH2F( "fhTOF" , "Scatter Plot TOF" , 240,-10,10,500,0,1.2);
     fListHist->Add(fhTOF);
   }
   if(! fhMassTOF){
@@ -385,6 +397,19 @@ void AliAnalysisTaskNucleiv2SP::UserCreateOutputObjects()
   hQVzAQVzCvsCentrality = new TH2F("hQVzAQVzCvsCentrality","hQVzAQVzCvsCentrality",1000,-5,5,105,-0.5,105.5);
   fListHist->Add(hQVzAQVzCvsCentrality);
  
+  hQxVzAvsCentrality = new TH2F("hQxVzAvsCentrality","hQxVzAvsCentrality",100,-5,5,105,-0.5,105.5);
+  hQyVzAvsCentrality = new TH2F("hQyVzAvsCentrality","hQyVzAvsCentrality",100,-5,5,105,-0.5,105.5);
+  hQxVzCvsCentrality = new TH2F("hQxVzCvsCentrality","hQxVzCvsCentrality",100,-5,5,105,-0.5,105.5);
+  hQyVzCvsCentrality = new TH2F("hQyVzCvsCentrality","hQyVzCvsCentrality",100,-5,5,105,-0.5,105.5);
+  hQxVzMvsCentrality = new TH2F("hQxVzMvsCentrality","hQxVzMvsCentrality",100,-5,5,105,-0.5,105.5);
+  hQyVzMvsCentrality = new TH2F("hQyVzMvsCentrality","hQyVzMvsCentrality",100,-5,5,105,-0.5,105.5);
+
+  fListHist->Add(hQxVzAvsCentrality);
+  fListHist->Add(hQyVzAvsCentrality);
+  fListHist->Add(hQxVzCvsCentrality);
+  fListHist->Add(hQyVzCvsCentrality);
+  fListHist->Add(hQxVzMvsCentrality);
+  fListHist->Add(hQyVzMvsCentrality);
 
   if(!ftree){
    
@@ -461,7 +486,7 @@ void AliAnalysisTaskNucleiv2SP::UserExec(Option_t *)
   ULong_t  status=0;
  
   Double_t pmax  = 10;
-  Double_t ptmax = 7;
+  Double_t ptmax = 10;
   // Primary vertex cut
   
   const AliESDVertex *vtx = lESDevent->GetPrimaryVertexTracks();
@@ -609,19 +634,34 @@ void AliAnalysisTaskNucleiv2SP::UserExec(Option_t *)
   Double_t  QV0AQV0C = qxEPa * qxEPc + qyEPa*qyEPc;
   hQVzAQVzCvsCentrality->Fill(QV0AQV0C,percentile);
   
+  //NUA correction
+
+  hQxVzAvsCentrality->Fill(qxEPa,percentile);
+  hQyVzAvsCentrality->Fill(qyEPa,percentile);
+  hQxVzCvsCentrality->Fill(qxEPc,percentile);
+  hQyVzCvsCentrality->Fill(qyEPc,percentile);
+  hQxVzMvsCentrality->Fill(qxEP ,percentile);
+  hQyVzMvsCentrality->Fill(qyEP ,percentile);
+
   //====================================================================================================================
   
   // To remove auto-correlation
   TVector2 *q = 0x0;
   q = pl->GetQVector();
 
-  Float_t ptcExp  = -999;
+  Double_t ptcExp  = -999;
   Double_t pullTPC = -999;
+  Double_t expbeta = -999;
+  Double_t pullTOF = -999;
   Float_t deltaphiTPC = -3;
   Float_t deltaphiV0  = -3;
   Float_t deltaphiV0A = -3;
   Float_t deltaphiV0C = -3;
 
+  Double_t massd   = 1.875612859;
+  Double_t masst   = 2.808939;
+  Double_t mass3he = 2.80892;
+
   Float_t  uqV0A = -999;
   Float_t  uqV0C = -999; 
 
@@ -646,7 +686,8 @@ void AliAnalysisTaskNucleiv2SP::UserExec(Option_t *)
     
     Double_t ptot = esdtrack->GetInnerParam()->GetP(); // momentum for dEdx determination
  
-    if(ptot<0.2)continue;
+    //   if(ptot<0.2)continue;
+    if(ptot<0.6)continue;
     fhBB->Fill(ptot*esdtrack->GetSign(),TPCSignal);
     esdtrack->GetImpactParameters(impactXY, impactZ);
               
@@ -668,37 +709,59 @@ void AliAnalysisTaskNucleiv2SP::UserExec(Option_t *)
     Float_t  gamma = 0;
     Float_t  mass  = -99;
     Double_t pt  = esdtrack->Pt();
+ 
+    if(fptc==1)
+      expbeta = TMath::Sqrt(1-((massd*massd)/(ptot*ptot+massd*massd))); 
+    if(fptc==2)
+      expbeta = TMath::Sqrt(1-((masst*masst)/(ptot*ptot+masst*masst))); 
+    if(fptc==3)
+      expbeta = TMath::Sqrt(1-((mass3he*mass3he)/(ptot*ptot+mass3he*mass3he))); 
+     
 
     if(fptc==3)
       pt = 2*pt;
 
-    if(TMath::Abs(ptot) < pmax && TMath::Abs(pullTPC) <= 3 && TMath::Abs(pt) < ptmax){
-
-      fhBBDeu->Fill(ptot*esdtrack->GetSign(),TPCSignal);
+    if(TMath::Abs(ptot) < pmax  && TMath::Abs(pt) < ptmax){
+      // if(TMath::Abs(pullTPC) <= 3)
       
       //
       // Process TOF information
       //
-      if (hasTOF) {
-       beta = length / (2.99792457999999984e-02 * tof);
-       gamma = 1/TMath::Sqrt(1 - beta*beta);
-       mass = ptot/TMath::Sqrt(gamma*gamma - 1); // using inner TPC mom. as approx.
-       
-       fhTOF->Fill(ptot*esdtrack->GetSign(),beta);
-       if(fptc==1){
-         if(TMath::Abs(mass) > 2.7)continue;
-         if(TMath::Abs(mass) < 1. )continue;
-       }
-       if(fptc==2){
-         if(TMath::Abs(mass) > 5.0)continue;
-         if(TMath::Abs(mass) < 1.8 )continue;
-       }
-       if(fptc==3){
-         if(TMath::Abs(mass) > 5.0)continue;
-         if(TMath::Abs(mass) < 1.8 )continue;
-       }
-       fhMassTOF->Fill(mass);
+      if(!hasTOF)continue;
+      //      if (hasTOF) {
+      beta = length / (2.99792457999999984e-02 * tof);
+      gamma = 1/TMath::Sqrt(1 - beta*beta);
+      mass = ptot/TMath::Sqrt(gamma*gamma - 1); // using inner TPC mom. as approx.
+  
+      //   cout<<expbeta<<" "<<beta<<" "<<(beta - expbeta)/(0.008*expbeta)<<endl;
+      pullTOF  = (beta - expbeta)/(0.007*expbeta);
+
+      if(TMath::Abs(ptot)< 2)
+       if(TMath::Abs(pullTPC) > 3)continue;
+      
+      // if(TMath::Abs(ptot)< 4)
+      //       if( beta>1)continue;
+      if(TMath::Abs(pullTOF) > 3)continue;
+
+      if(TMath::Sqrt(esdtrack->GetTOFsignalDz()*esdtrack->GetTOFsignalDz() + esdtrack->GetTOFsignalDx()*esdtrack->GetTOFsignalDx()) > 5.)continue; 
+
+      if(fptc==1){
+       if(TMath::Abs(mass) > 2.65)continue;
+       if(TMath::Abs(mass) < 1.05)continue;
       }
+      if(fptc==2){
+       if(TMath::Abs(mass) > 5.0)continue;
+       if(TMath::Abs(mass) < 1.8 )continue;
+      }
+      if(fptc==3){
+       if(TMath::Abs(mass) > 5.0)continue;
+       if(TMath::Abs(mass) < 1.8)continue;
+      }
+      fhMassTOF->Fill(mass);
+      //}
+      fhTOF->Fill(ptot*esdtrack->GetSign(),beta);
+      fhBBDeu->Fill(ptot*esdtrack->GetSign(),TPCSignal);
+      
    
       // Event Plane
       // Remove AutoCorrelation