added a check for TPC gaussianity

diff --git a/PWGPP/pid/doeffAllSpecies.C b/PWGPP/pid/doeffAllSpecies.C
index deb79ef77011b38ad085775abd41326e2c764619..4946626ad3c61830f6b6e1d75840fb6ba9115479 100644 (file)
--- a/PWGPP/pid/doeffAllSpecies.C
+++ b/PWGPP/pid/doeffAllSpecies.C
@@ -33,6 +33,7 @@ performAllPi(){
   rebinsize = 1; // don't change this, not choice here!!!
   kGoodMatch=kFALSE;
   kSigma2vs3 = kFALSE;
+  kSigma2vs3TPC = kFALSE;
   require5sigma = kFALSE;
   bayesVsigma = kFALSE;
   kTOFmatch = kFALSE;
@@ -68,6 +69,12 @@ performAllPi(){
   doeffPi(1,0.1,etamin,etamax);
   doeffPi(0,0.1,etamin,etamax);
   kSigma2vs3=kFALSE;
+
+  // eff 2 sigma / eff 3 sigma TPC
+  kSigma2vs3TPC=kTRUE;
+  doeffPi(1,0.1,etamin,etamax);
+  doeffPi(0,0.1,etamin,etamax);
+  kSigma2vs3TPC=kFALSE;
   
   kOverAll=kTRUE;
   // TPC|TOF overall eff
@@ -130,6 +137,7 @@ performAllPr(){
   rebinsize = 1; // don't change this, not choice here!!!
   kGoodMatch=kFALSE;
   kSigma2vs3 = kFALSE;
+  kSigma2vs3TPC = kFALSE;
   require5sigma = kFALSE;
   bayesVsigma = kFALSE;
   kTOFmatch = kFALSE;
@@ -166,6 +174,12 @@ performAllPr(){
   doeffPr(0,0.1,etamin,etamax);
   kSigma2vs3=kFALSE;
 
+  // eff 2 sigma / eff 3 sigma TPC
+  kSigma2vs3TPC=kTRUE;
+  doeffPr(1,0.1,etamin,etamax);
+  doeffPr(0,0.1,etamin,etamax);
+  kSigma2vs3TPC=kFALSE;
+
   kOverAll=kTRUE;
   // TPC|TOF overall eff
   doeffPr(1,0.2,etamin,etamax);
@@ -227,6 +241,7 @@ performAllKa(){
   rebinsize = 1; // don't change this, not choice here!!!
   kGoodMatch=kFALSE;
   kSigma2vs3 = kFALSE;
+  kSigma2vs3TPC = kFALSE;
   require5sigma = kFALSE;
   bayesVsigma = kFALSE;
   kTOFmatch = kFALSE;
@@ -263,6 +278,12 @@ performAllKa(){
   doeffKa(0,0.1,etamin,etamax);
   kSigma2vs3=kFALSE;
 
+  // eff 2 sigma / eff 3 sigma TPC
+  kSigma2vs3TPC=kTRUE;
+  doeffKa(1,0.1,etamin,etamax);
+  doeffKa(0,0.1,etamin,etamax);
+  kSigma2vs3TPC=kFALSE;
+
   kOverAll=kTRUE;
   // TPC|TOF overall eff
   doeffKa(1,0.2,etamin,etamax);

diff --git a/PWGPP/pid/doeffKa.C b/PWGPP/pid/doeffKa.C
index 7ed6f6b660874a7173fcce6db23d01491e296b32..bea1198b696dceb6bab986bb9cd849fea72849ea 100644 (file)
--- a/PWGPP/pid/doeffKa.C
+++ b/PWGPP/pid/doeffKa.C
@@ -42,7 +42,7 @@ Float_t fitmax = 1.045;
 Int_t cmin = 1;// min 1
 Int_t cmax = 8;// max 10
 
-Float_t weightS = -1.;
+Float_t weightS = -0.9;
 
 Int_t rebinsize = 1;
 
@@ -55,6 +55,7 @@ Bool_t keepTrue = kFALSE; // put it to false to fit only background (only for MC
 Bool_t kGoodMatch = kFALSE; // to check good matching
 
 Bool_t kSigma2vs3 = kFALSE; // to check good matching
+Bool_t kSigma2vs3TPC = kFALSE; // to check good matching
 
 Bool_t require5sigma = kFALSE; // don't touch this flag
 
@@ -71,7 +72,7 @@ TH2F *htracked;
 
 Bool_t kLoaded=kFALSE;
 int LoadLib(){
-  weightS = -1.;
+  weightS = -0.9;
 
   require5sigma = kFALSE;
 
@@ -144,10 +145,16 @@ int LoadLib(){
     kOverAll = 0;
   }
 
+  if(kSigma2vs3TPC){
+    kGoodMatch=kFALSE;
+    kOverAll = 0;
+  }
+
   if(bayesVsigma){
     kOverAll = 0;
     kGoodMatch=kFALSE;
     kSigma2vs3=kFALSE;
+    kSigma2vs3TPC=kFALSE;
     kTOFmatch=kTRUE;
     weightS = -0.7;
   }
@@ -379,8 +386,9 @@ void doeffKa(Int_t pos,Float_t prob,Float_t etaminkp,Float_t etamaxkp){
   if(prob > 0.1999){
     kGoodMatch = kFALSE;
     kSigma2vs3 = kFALSE;
+    kSigma2vs3TPC=kFALSE;
     if(! kOverAll) require5sigma = kTRUE;
-    if(!isMC && !kOverAll) weightS = -0.95;
+    if(!isMC && !kOverAll) weightS = -0.9;
   }
 
   TCanvas *c1 = new TCanvas();
@@ -530,6 +538,9 @@ void doeffKa(Int_t pos,Float_t prob,Float_t etaminkp,Float_t etamaxkp){
   if(kSigma2vs3)
     sprintf(flag2,"Sigma2vs3");
 
+  if(kSigma2vs3TPC)
+    sprintf(flag2,"Sigma2vs3TPC");
+
   if(kOverAll)
     sprintf(flag2,"OverAll");
   if(kOverAllTOFmatch)
@@ -541,7 +552,7 @@ void doeffKa(Int_t pos,Float_t prob,Float_t etaminkp,Float_t etamaxkp){
     if(prob >=0.2) sprintf(name,"kaonPos%sP%iEff%i_%i%s%s.root",etarange,Int_t(prob*100),(cmin-1)*10,cmax*10,flag,flag2);
     else{
       sprintf(name,"kaonPos%sMatchEff%i_%i%s%s.root",etarange,(cmin-1)*10,cmax*10,flag,flag2);
-      if(!(kOverAll || bayesVsigma || kGoodMatch || kSigma2vs3)) kWriteME = kTRUE;
+      if(!(kOverAll || bayesVsigma || kGoodMatch || kSigma2vs3 || kSigma2vs3TPC)) kWriteME = kTRUE;
     }
   }
   else{
@@ -597,6 +608,14 @@ TH2F *GetHistoKap(Float_t pt,Float_t ptM,Float_t pMinkp,Float_t pMinkn,Float_t e
       x2[9] = 2;
   }
 
+  if(kSigma2vs3TPC){
+    x[6] = 0.0001;
+    x2[6] = 0.0002;
+    x2[9] = 3;
+    if(pMinkp > 0)
+      x2[9] = 2;
+  }
+
   if(bayesVsigma){
     if(pMinkp > 0){
       x[4] = 0.2001;
@@ -659,6 +678,14 @@ TH2F *GetHistoKan(Float_t pt,Float_t ptM,Float_t pMinkn,Float_t pMinkp,Float_t e
     if(pMinkn > 0)
       x2[10] = 2;
   }
+
+  if(kSigma2vs3TPC){
+    x[7] = 0.0001;
+    x2[7] = 0.0002;
+    x2[10] = 3;
+    if(pMinkn > 0)
+      x2[10] = 2;
+  }
  
  if(bayesVsigma){
     if(pMinkn > 0){

diff --git a/PWGPP/pid/doeffPi.C b/PWGPP/pid/doeffPi.C
index 780acddf799a1502f6ed131687ee63ee0276690f..a80af449a53133cef3dce02487fc1fc723d6eb68 100644 (file)
--- a/PWGPP/pid/doeffPi.C
+++ b/PWGPP/pid/doeffPi.C
@@ -43,7 +43,7 @@ Float_t fitmax = 0.7;
 Int_t cmin = 1; // min 1
 Int_t cmax = 10;//max 10
 
-Float_t weightS = -1.;
+Float_t weightS = -0.9;
 
 Int_t rebinsize = 1;
 
@@ -56,6 +56,7 @@ Bool_t keepTrue = kFALSE; // put it to false to fit only background (only for MC
 Bool_t kGoodMatch = kFALSE; // to check good matching
 
 Bool_t kSigma2vs3 = kFALSE; // to check good matching
+Bool_t kSigma2vs3TPC = kFALSE; // to check good matching
 
 Bool_t require5sigma = kFALSE; // don't touch this flag
 
@@ -73,7 +74,7 @@ TH2F *htracked;
 Bool_t kLoaded=kFALSE;
 
 Int_t LoadLib(){
-  weightS = -1.;
+  weightS = -0.9;
 
   require5sigma = kFALSE;
 
@@ -137,10 +138,16 @@ Int_t LoadLib(){
     kOverAll = 0;
   }
 
+  if(kSigma2vs3TPC){
+    kGoodMatch=kFALSE;
+    kOverAll = 0;
+  }
+
   if(bayesVsigma){
     kOverAll = 0;
     kGoodMatch=kFALSE;
     kSigma2vs3=kFALSE;
+    kSigma2vs3TPC=kFALSE;
     kTOFmatch=kTRUE;
     weightS = -0.7;
   }
@@ -368,8 +375,9 @@ void doeffPi(Int_t pos,Float_t prob,Float_t etaminkp,Float_t etamaxkp){
   if(prob > 0.1999){
     kGoodMatch = kFALSE;
     kSigma2vs3 = kFALSE;
+    kSigma2vs3TPC = kFALSE;
     if(! kOverAll) require5sigma = kTRUE;
-    if(!isMC) weightS = -0.95;
+    if(!isMC) weightS = -0.9;
   }
 
   TCanvas *c1 = new TCanvas();
@@ -519,6 +527,9 @@ void doeffPi(Int_t pos,Float_t prob,Float_t etaminkp,Float_t etamaxkp){
   if(kSigma2vs3)
     sprintf(flag2,"Sigma2vs3");
 
+  if(kSigma2vs3TPC)
+    sprintf(flag2,"Sigma2vs3TPC");
+
   if(kOverAll)
     sprintf(flag2,"OverAll");
   if(kOverAllTOFmatch)
@@ -530,7 +541,7 @@ void doeffPi(Int_t pos,Float_t prob,Float_t etaminkp,Float_t etamaxkp){
     if(prob >=0.2) sprintf(name,"pionPos%sP%iEff%i_%i%s%s.root",etarange,Int_t(prob*100),(cmin-1)*10,cmax*10,flag,flag2);
     else{
       sprintf(name,"pionPos%sMatchEff%i_%i%s%s.root",etarange,(cmin-1)*10,cmax*10,flag,flag2);
-      if(!(kOverAll || bayesVsigma || kGoodMatch || kSigma2vs3)) kWriteME = kTRUE;
+      if(!(kOverAll || bayesVsigma || kGoodMatch || kSigma2vs3 || kSigma2vs3TPC)) kWriteME = kTRUE;
     }
   }
   else{
@@ -586,6 +597,14 @@ TH2F *GetHistoPip(Float_t pt,Float_t ptM,Float_t pMinkp,Float_t pMinkn,Float_t e
       x2[9] = 2;
   }
 
+  if(kSigma2vs3TPC){
+    x[6] = 0.0001;
+    x2[6] = 0.0002;
+    x2[9] = 3;
+    if(pMinkp > 0)
+      x2[9] = 2;
+  }
+
   if(bayesVsigma){
     if(pMinkp > 0){
       x[4] = 0.2001;
@@ -649,6 +668,14 @@ TH2F *GetHistoPin(Float_t pt,Float_t ptM,Float_t pMinkn,Float_t pMinkp,Float_t e
       x2[10] = 2;
   }
  
+  if(kSigma2vs3TPC){
+    x[7] = 0.0001;
+    x2[7] = 0.0002;
+    x2[10] = 3;
+    if(pMinkn > 0)
+      x2[10] = 2;
+  }
+
  if(bayesVsigma){
     if(pMinkn > 0){
       x[5] = 0.2001;

diff --git a/PWGPP/pid/doeffPr.C b/PWGPP/pid/doeffPr.C
index 74cf5e25bb7da372e720773a3ede2a246afc9bf1..074674951f917be8f9f2d6fc8d8367a80a6518fe 100644 (file)
--- a/PWGPP/pid/doeffPr.C
+++ b/PWGPP/pid/doeffPr.C
@@ -43,7 +43,7 @@ Float_t fitmax = 1.15;
 Int_t cmin = 1; // min 1
 Int_t cmax = 10; // max 10
 
-Float_t weightS = -1.;
+Float_t weightS = -0.9;
 
 Int_t rebinsize = 1;
 
@@ -56,6 +56,7 @@ Bool_t keepTrue = kFALSE; // put it to false to fit only background (only for MC
 Bool_t kGoodMatch = kFALSE; // to check good matching
 
 Bool_t kSigma2vs3 = kFALSE; // to check good matching
+Bool_t kSigma2vs3TPC = kFALSE; // to check good matching
 
 Bool_t require5sigma = kFALSE; // don't touch this flag
 
@@ -72,7 +73,7 @@ TH2F *htracked;
 
 Bool_t kLoaded=kFALSE;
 Int_t LoadLib(){
-  weightS = -1.;
+  weightS = -0.9;
 
   require5sigma = kFALSE;
 
@@ -137,10 +138,16 @@ Int_t LoadLib(){
     kOverAll = 0;
   }
 
+  if(kSigma2vs3TPC){
+    kGoodMatch=kFALSE;
+    kOverAll = 0;
+  }
+
   if(bayesVsigma){
     kOverAll = 0;
     kGoodMatch=kFALSE;
     kSigma2vs3=kFALSE;
+    kSigma2vs3TPC = kFALSE;
     kTOFmatch=kTRUE;
     weightS = -0.7;
   }
@@ -368,8 +375,9 @@ void doeffPr(Int_t pos,Float_t prob,Float_t etaminkp,Float_t etamaxkp){
   if(prob > 0.1999){
     kGoodMatch = kFALSE;
     kSigma2vs3 = kFALSE;
+    kSigma2vs3TPC = kFALSE;
     if(! kOverAll) require5sigma = kTRUE;
-    if(!isMC) weightS = -0.95;
+    if(!isMC) weightS = -0.9;
   }
 
   TCanvas *c1 = new TCanvas();
@@ -519,6 +527,9 @@ void doeffPr(Int_t pos,Float_t prob,Float_t etaminkp,Float_t etamaxkp){
   if(kSigma2vs3)
     sprintf(flag2,"Sigma2vs3");
 
+  if(kSigma2vs3TPC)
+    sprintf(flag2,"Sigma2vs3TPC");
+
   if(kOverAll)
     sprintf(flag2,"OverAll");
   if(kOverAllTOFmatch)
@@ -530,7 +541,7 @@ void doeffPr(Int_t pos,Float_t prob,Float_t etaminkp,Float_t etamaxkp){
     if(prob >=0.2) sprintf(name,"protonPos%sP%iEff%i_%i%s%s.root",etarange,Int_t(prob*100),(cmin-1)*10,cmax*10,flag,flag2);
     else{
       sprintf(name,"protonPos%sMatchEff%i_%i%s%s.root",etarange,(cmin-1)*10,cmax*10,flag,flag2);
-      if(!(kOverAll || bayesVsigma || kGoodMatch || kSigma2vs3)) kWriteME = kTRUE;
+      if(!(kOverAll || bayesVsigma || kGoodMatch || kSigma2vs3 || kSigma2vs3TPC)) kWriteME = kTRUE;
     }
   }
   else{
@@ -586,6 +597,14 @@ TH2F *GetHistoPrp(Float_t pt,Float_t ptM,Float_t pMinkp,Float_t pMinkn,Float_t e
       x2[9] = 2;
   }
 
+  if(kSigma2vs3TPC){
+    x[6] = 0.0001;
+    x2[6] = 0.0002;
+    x2[9] = 3;
+    if(pMinkp > 0)
+      x2[9] = 2;
+  }
+
   if(bayesVsigma){
     if(pMinkp > 0){
       x[4] = 0.2001;
@@ -648,6 +667,14 @@ TH2F *GetHistoPrn(Float_t pt,Float_t ptM,Float_t pMinkn,Float_t pMinkp,Float_t e
     if(pMinkn > 0)
       x2[10] = 2;
   }
+
+  if(kSigma2vs3TPC){
+    x[7] = 0.0001;
+    x2[7] = 0.0002;
+    x2[10] = 3;
+    if(pMinkn > 0)
+      x2[10] = 2;
+  }
  
  if(bayesVsigma){
     if(pMinkn > 0){