D+ cuts updated (Renu)

diff --git a/PWG3/vertexingHF/ConfigVertexingHF.C b/PWG3/vertexingHF/ConfigVertexingHF.C
index a29b33c597723d61f8e70e886cf5b285f710c0b3..3ccd950cd428a5e2dd16dd098c8589f108db84dd 100644 (file)
--- a/PWG3/vertexingHF/ConfigVertexingHF.C
+++ b/PWG3/vertexingHF/ConfigVertexingHF.C
@@ -52,8 +52,8 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   cutsJpsitoee->AddTrackCuts(esdTrackCuts);
   vHF->SetCutsJpsitoee(cutsJpsitoee);
   AliRDHFCutsDplustoKpipi *cutsDplustoKpipi = new AliRDHFCutsDplustoKpipi("CutsDplustoKpipi");
-  Float_t cutsArrayDplustoKpipi[12]={0.2,0.3,0.3,0.,0.,0.01,0.06,0.02,0.,0.7,0.,10000000000.};
-  cutsDplustoKpipi->SetCuts(12,cutsArrayDplustoKpipi);
+  Float_t cutsArrayDplustoKpipi[14]={0.2,0.3,0.3,0.,0.,0.01,0.06,0.02,0.,0.7,0.,10000000000.,0.,0.};
+  cutsDplustoKpipi->SetCuts(14,cutsArrayDplustoKpipi);
   cutsDplustoKpipi->AddTrackCuts(esdTrackCuts);
   vHF->SetCutsDplustoKpipi(cutsDplustoKpipi);
   AliRDHFCutsDstoKKpi *cutsDstoKKpi = new AliRDHFCutsDstoKKpi("CutsDstoKKpi");

diff --git a/PWG3/vertexingHF/ConfigVertexingHF_highmult.C b/PWG3/vertexingHF/ConfigVertexingHF_highmult.C
index e191b5598ce910a13c74d2db4e7ceb5916565f32..7582b44db623eaff0a3d33fa277d8feefccc661e 100644 (file)
--- a/PWG3/vertexingHF/ConfigVertexingHF_highmult.C
+++ b/PWG3/vertexingHF/ConfigVertexingHF_highmult.C
@@ -65,9 +65,9 @@ AliAnalysisVertexingHF* ConfigVertexingHF() {
   AliRDHFCutsDplustoKpipi *cutsDplustoKpipi = new AliRDHFCutsDplustoKpipi("CutsDplustoKpipi");
   cutsDplustoKpipi->SetStandardCutsPbPb2010();
   cutsDplustoKpipi->SetUsePID(kFALSE);
-  Float_t cutsArrayDplustoKpipi[12]={0.2,0.8,0.8,0.,0.,0.01,0.06,0.05,1.,0.85,0.,10000000000.};
+  Float_t cutsArrayDplustoKpipi[14]={0.2,0.8,0.8,0.,0.,0.01,0.06,0.05,1.,0.85,0.,10000000000.,0.,0.};
   cutsDplustoKpipi->SetPtBins(nptbins,ptlimits);
-  cutsDplustoKpipi->SetCuts(12,cutsArrayDplustoKpipi);
+  cutsDplustoKpipi->SetCuts(14,cutsArrayDplustoKpipi);
   cutsDplustoKpipi->AddTrackCuts(esdTrackCuts);
   cutsDplustoKpipi->SetMinPtCandidate(3.);
   vHF->SetCutsDplustoKpipi(cutsDplustoKpipi);

diff --git a/PWG3/vertexingHF/macros/makeTFile4CutsDplustoKpipi.C b/PWG3/vertexingHF/macros/makeTFile4CutsDplustoKpipi.C
index 8a55bf8aaef48ea629892b5d20598fb5c403cb23..c9f6eef39c0aa0e36ba2e7febc0fdc11df28b5dd 100644 (file)
--- a/PWG3/vertexingHF/macros/makeTFile4CutsDplustoKpipi.C
+++ b/PWG3/vertexingHF/macros/makeTFile4CutsDplustoKpipi.C
@@ -20,7 +20,7 @@
 //        Giacomo Ortona, ortona@to.infn.it
 //        Renu Bala [Dplus Analysis and CF]
 
-void makeInputAliAnalysisTaskSEDplus(){
+void makeInputAliAnalysisTaskSEDplusPP(){
 
  //  gSystem->SetIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT -I$ALICE_ROOT/include -I$ALICE_ROOT/ITS -I$ALICE_ROOT/TPC -I$ALICE_ROOT/CONTAINERS -I$ALICE_ROOT/STEER -I$ALICE_ROOT/TRD -I$ALICE_ROOT/macros -I$ALICE_ROOT/ANALYSIS -I$ALICE_ROOT/PWG3 -I$ALICE_ROOT/PWG3/vertexingHF -I$ALICE_ROOT/PWG3/vertexingH/macros -g"); 
 
@@ -40,7 +40,7 @@ void makeInputAliAnalysisTaskSEDplus(){
     
     
     
-    const Int_t nptbins=10;
+    const Int_t nptbins=14;
     Float_t* ptbins;
     ptbins=new Float_t[nptbins+1];
     ptbins[0]=0.;
@@ -50,11 +50,15 @@ void makeInputAliAnalysisTaskSEDplus(){
     ptbins[4]=4.;
     ptbins[5]=5.;
     ptbins[6]=6.;
-    ptbins[7]=8.;
-    ptbins[8]=12.;
-    ptbins[9]=16.;
-    ptbins[10]=24.;
-    const Int_t nvars=12;
+    ptbins[7]=7.;
+    ptbins[8]=8.;
+    ptbins[9]=9.;
+    ptbins[10]=10.;
+    ptbins[11]=12.;
+    ptbins[12]=14.;
+    ptbins[13]=16.;
+    ptbins[14]=24.;
+    const Int_t nvars=14;
     
     
     Float_t** prodcutsval;
@@ -73,6 +77,8 @@ void makeInputAliAnalysisTaskSEDplus(){
       prodcutsval[9][ipt]=0.85;
       prodcutsval[10][ipt]=0.;
       prodcutsval[11][ipt]=10000000.0;
+      prodcutsval[12][ipt]=0.0;
+      prodcutsval[13][ipt]=0.0;
       
     }
     
@@ -86,80 +92,76 @@ void makeInputAliAnalysisTaskSEDplus(){
     for(Int_t ipt=0;ipt<nptbins;ipt++){
       anacutsval[ic][ipt]=0.2;
     }
+    ic=1;
+    for(Int_t ipt=3;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.4;
+    }
+    anacutsval[1][0]=0.3;
+    anacutsval[1][1]=0.3;
+    anacutsval[1][2]=0.4;
+
+    ic=2;
+    for(Int_t ipt=3;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.4;
+    }
+    
+    anacutsval[2][0]=0.3;
+    anacutsval[2][1]=0.3;
+    anacutsval[2][2]=0.4;
     
     ic=3;
     for(Int_t ipt=0;ipt<nptbins;ipt++){
       anacutsval[ic][ipt]=0.;
     }
+    
     ic=4;
     for(Int_t ipt=0;ipt<nptbins;ipt++){
       anacutsval[ic][ipt]=0.;
     }
+    
     ic=5;
     for(Int_t ipt=0;ipt<nptbins;ipt++){
       anacutsval[ic][ipt]=0.01;
     }
-    //ic=6;
-    //for(Int_t ipt=0;ipt<nptbins;ipt++){
-    // anacutsval[ic][ipt]=0.06;
-    
-    ic=11;
-    for(Int_t ipt=0;ipt<nptbins;ipt++){
-      anacutsval[ic][ipt]=10000000000.;
-    }
-    
-    
-    anacutsval[1][0]=0.3;
-    anacutsval[1][1]=0.3;
-    anacutsval[1][2]=0.4;
-    
-    anacutsval[2][0]=0.3;
-    anacutsval[2][1]=0.3;
-    anacutsval[2][2]=0.4;
     
     
-    ic=1;
-    for(Int_t ipt=3;ipt<nptbins;ipt++){
-      anacutsval[ic][ipt]=0.4;
+    ic=6;
+    for(Int_t ipt=5;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.023333;
     }
-    ic=2;
-    for(Int_t ipt=3;ipt<nptbins;ipt++){
-      anacutsval[ic][ipt]=0.4;
-    }
-    
     
     anacutsval[6][0]=0.022100;
     anacutsval[6][1]=0.022100;
     anacutsval[6][2]=0.034;
     anacutsval[6][3]=0.020667;
     anacutsval[6][4]=0.020667;
-    anacutsval[6][5]=0.023333;
-    anacutsval[6][6]=0.023333;
-    anacutsval[6][7]=0.023333;
-    anacutsval[6][8]=0.023333;
-    anacutsval[6][9]=0.023333;
-     
+    
+    
+    ic=7;
+    for(Int_t ipt=5;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.115;
+    }
+    
     anacutsval[7][0]=0.08;
     anacutsval[7][1]=0.08;
     anacutsval[7][2]=0.09;  
     anacutsval[7][3]=0.095;
     anacutsval[7][4]=0.095;
-    anacutsval[7][5]=0.115;
-    anacutsval[7][6]=0.115;
-    anacutsval[7][7]=0.115;
-    anacutsval[7][8]=0.115;
-    anacutsval[7][9]=0.115;
+    
+    ic=8;
+    for(Int_t ipt=3;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.5;
+    }
      
     anacutsval[8][0]=0.5;
     anacutsval[8][1]=0.5;
     anacutsval[8][2]=1.0;
-    anacutsval[8][3]=0.5;
-    anacutsval[8][4]=0.5;
-    anacutsval[8][5]=0.5;
-    anacutsval[8][6]=0.5;
-    anacutsval[8][7]=0.5;
-    anacutsval[8][8]=0.5;
-    anacutsval[8][9]=0.5;
+  
+   
+    ic=9;
+    for(Int_t ipt=9;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.90;
+    }
      
     anacutsval[9][0]=0.95;
     anacutsval[9][1]=0.95;
@@ -170,19 +172,36 @@ void makeInputAliAnalysisTaskSEDplus(){
     anacutsval[9][6]=0.92;
     anacutsval[9][7]=0.92;
     anacutsval[9][8]=0.92;
-    anacutsval[9][9]=0.92;
+  
+
+    ic=10;
+    for(Int_t ipt=3;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.000883;
+    }
  
     anacutsval[10][0]=0.0055;
     anacutsval[10][1]=0.0055;
     anacutsval[10][2]= 0.0028;
-    anacutsval[10][3]=0.000883;
-    anacutsval[10][4]=0.000883;
-    anacutsval[10][5]=0.000883;
-    anacutsval[10][6]=0.000883;
-    anacutsval[10][7]=0.000883;
-    anacutsval[10][8]=0.000883;
-    anacutsval[10][9]=0.000883;
+  
+  
+
+
+    ic=11;
+    for(Int_t ipt=0;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=10000000000.;
+    }
+
+    ic=12;
+    for(Int_t ipt=0;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.;
+    }
+
+    ic=13;
+    for(Int_t ipt=0;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.;
+    }
     
+          
     
     AliRDHFCutsDplustoKpipi *prodcuts = new AliRDHFCutsDplustoKpipi();
     prodcuts->SetName("ProdCuts");
@@ -197,6 +216,7 @@ void makeInputAliAnalysisTaskSEDplus(){
     analysiscuts->SetCuts(nvars,nptbins,anacutsval);
     analysiscuts->AddTrackCuts(esdTrackCuts);
     analysiscuts->SetUsePID(kTRUE);
+    analysiscuts->SetUseImpParProdCorrCut(kFALSE);
     cout<<"This is the odject I'm going to save:"<<nptbins<<endl;
     
     analysiscuts->PrintAll();
@@ -210,6 +230,267 @@ void makeInputAliAnalysisTaskSEDplus(){
 }
 
 
+
+
+
+
+void makeInputAliAnalysisTaskSEDplusPbPb(){
+
+ //  gSystem->SetIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT -I$ALICE_ROOT/include -I$ALICE_ROOT/ITS -I$ALICE_ROOT/TPC -I$ALICE_ROOT/CONTAINERS -I$ALICE_ROOT/STEER -I$ALICE_ROOT/TRD -I$ALICE_ROOT/macros -I$ALICE_ROOT/ANALYSIS -I$ALICE_ROOT/PWG3 -I$ALICE_ROOT/PWG3/vertexingHF -I$ALICE_ROOT/PWG3/vertexingH/macros -g"); 
+
+ AliESDtrackCuts* esdTrackCuts=new AliESDtrackCuts();
+  esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
+  //default
+    esdTrackCuts->SetRequireTPCRefit(kTRUE);
+    esdTrackCuts->SetRequireITSRefit(kTRUE);
+    //esdTrackCuts->SetMinNClustersITS(4); // default is 5
+    esdTrackCuts->SetMinNClustersTPC(70);
+    esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
+                                          AliESDtrackCuts::kAny); 
+    // default is kBoth, otherwise kAny
+    esdTrackCuts->SetMinDCAToVertexXY(0.);
+    esdTrackCuts->SetPtRange(0.8,1.e10);
+    
+    TString cent="";
+    //centrality selection (Pb-Pb)                                              
+    Float_t minc=0.,maxc=20.;
+  const Int_t nptbins=12;
+    Float_t* ptbins;
+    ptbins=new Float_t[nptbins+1];
+    ptbins[0]=0.;
+    ptbins[1]=1.;
+    ptbins[2]=2.;
+    ptbins[3]=3.;
+    ptbins[4]=4.;
+    ptbins[5]=5.;
+    ptbins[6]=6.;
+    ptbins[7]=8.;
+    ptbins[8]=9.;
+    ptbins[9]=10.;
+    ptbins[10]=12.;
+    ptbins[11]=16.;
+    ptbins[12]=24.;
+    const Int_t nvars=14;
+        
+        
+    
+    Float_t** prodcutsval;
+    prodcutsval=new Float_t*[nvars];
+    for(Int_t ic=0;ic<nvars;ic++){prodcutsval[ic]=new Float_t[nptbins];}  
+    for(Int_t ipt=0;ipt<nptbins;ipt++){
+      prodcutsval[0][ipt]=0.2;
+      prodcutsval[1][ipt]=0.8;
+      prodcutsval[2][ipt]=0.8;
+      prodcutsval[3][ipt]=0.;
+      prodcutsval[4][ipt]=0.;
+      prodcutsval[5][ipt]=0.01;
+      prodcutsval[6][ipt]=0.06;
+      prodcutsval[7][ipt]=0.02;
+      prodcutsval[8][ipt]=1.;
+      prodcutsval[9][ipt]=0.90;
+      prodcutsval[10][ipt]=0.;
+      prodcutsval[11][ipt]=10000000.0;
+      prodcutsval[12][ipt]=0.0;
+      prodcutsval[13][ipt]=0.0;
+      
+    }
+    
+    
+   
+
+  Float_t** anacutsval;
+    anacutsval=new Float_t*[nvars];
+    for(Int_t ic=0;ic<nvars;ic++){anacutsval[ic]=new Float_t[nptbins];}
+    
+        
+     Int_t ic=0;
+    for(Int_t ipt=0;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.2;
+    }
+    
+    
+    ic=1;
+
+    for(Int_t ipt=7;ipt<nptbins;ipt++){
+      anacutsval[1][ipt]=1.0;
+    }
+    
+    anacutsval[1][0]=0.8;
+    anacutsval[1][1]=0.8;
+    anacutsval[1][2]=0.8;
+    anacutsval[1][3]=0.8;
+    anacutsval[1][4]=0.9;
+    anacutsval[1][5]=0.9;
+    anacutsval[1][6]=0.8;
+ 
+
+
+
+    ic=2;    
+    for(Int_t ipt=0;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.8;
+    }
+    
+
+
+
+    ic=3;
+    for(Int_t ipt=0;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.;
+    }
+    ic=4;
+    for(Int_t ipt=0;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.;
+    }
+    ic=5;
+    for(Int_t ipt=0;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.01;
+    }
+
+
+    ic=6;
+    for(Int_t ipt=5;ipt<nptbins;ipt++){
+     anacutsval[ic][ipt]=0.023333;
+    }
+     anacutsval[6][0]=0.022100;
+     anacutsval[6][1]=0.022100;
+     anacutsval[6][2]=0.034;
+     anacutsval[6][3]=0.020667;
+     anacutsval[6][4]=0.020667;
+    
+    
+     ic=7;
+    for(Int_t ipt=7;ipt<nptbins;ipt++){
+     anacutsval[ic][ipt]=0.19;
+    }
+    
+    anacutsval[7][0]=0.08;
+    anacutsval[7][1]=0.08;
+    anacutsval[7][2]=0.17;  
+    anacutsval[7][3]=0.14;
+    anacutsval[7][4]=0.14;
+    anacutsval[7][5]=0.19;
+    anacutsval[7][6]=0.14;
+  
+
+    
+
+    ic=8;
+    for(Int_t ipt=5;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=2.0;
+    }
+
+    anacutsval[8][0]=0.8;
+    anacutsval[8][1]=0.8;
+    anacutsval[8][2]=1.1;
+    anacutsval[8][3]=0.5;
+    anacutsval[8][4]=0.5;
+   
+    
+
+    ic=9;
+    for(Int_t ipt=7;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.997;
+    }
+
+    anacutsval[9][0]=0.995;
+    anacutsval[9][1]=0.995;
+    anacutsval[9][2]=0.997; 
+    anacutsval[9][3]=0.998; 
+    anacutsval[9][4]= 0.998;
+    anacutsval[9][5]=0.995;
+    anacutsval[9][6]=0.995;
+   
+
+
+
+
+
+    ic=10;
+    for(Int_t ipt=3;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.000883;
+    }   
+
+    anacutsval[10][0]=0.0055;
+    anacutsval[10][1]=0.0055;
+    anacutsval[10][2]= 0.0028;
+       
+
+
+    ic=11;
+    for(Int_t ipt=0;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=10000000000.;
+    }
+   
+         
+     ic=12;
+    for(Int_t ipt=7;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.;
+    }
+
+    anacutsval[12][0]=0.;
+    anacutsval[12][1]=0.;
+
+    anacutsval[12][2]=0.;
+    anacutsval[12][3]=12.;
+    anacutsval[12][4]=12.;
+    anacutsval[12][5]=12.0;
+    anacutsval[12][6]=10.;
+  
+
+    
+  ic=13;
+    for(Int_t ipt=7;ipt<nptbins;ipt++){
+      anacutsval[ic][ipt]=0.;
+    }
+    
+    anacutsval[13][0]=0.;
+    anacutsval[13][1]=0.;
+
+    anacutsval[13][2]=0.;
+    anacutsval[13][3]=0.998571;
+    anacutsval[13][4]=0.998571;
+    anacutsval[13][5]=0.998571;
+    anacutsval[13][6]=0.997143;
+   
+
+    
+    AliRDHFCutsDplustoKpipi *prodcuts = new AliRDHFCutsDplustoKpipi();
+    prodcuts->SetName("ProdCuts");
+    prodcuts->SetPtBins(nptbins+1,ptbins);
+    prodcuts->SetCuts(nvars,nptbins,prodcutsval);
+    
+    
+    AliRDHFCutsDplustoKpipi* analysiscuts=new AliRDHFCutsDplustoKpipi();
+    analysiscuts->SetName("AnalysisCuts");
+    analysiscuts->SetTitle("Cuts for Dplus Analysis and CF");
+    analysiscuts->SetPtBins(nptbins+1,ptbins);
+    analysiscuts->SetCuts(nvars,nptbins,anacutsval);
+    analysiscuts->AddTrackCuts(esdTrackCuts);
+    
+    analysiscuts->SetUsePID(kTRUE);
+    analysiscuts->SetUseImpParProdCorrCut(kFALSE);
+    analysiscuts->SetOptPileup(kFALSE);
+    analysiscuts->SetUseAOD049(kTRUE);
+    analysiscuts->SetMinCentrality(minc);
+    analysiscuts->SetMaxCentrality(maxc);
+    cent=Form("%.0f%.0f",minc,maxc);
+    analysiscuts->SetUseCentrality(AliRDHFCuts::kCentV0M); //kCentOff,kCentV0M,kCentTRK,kCentTKL,kCentCL1,kCentInvalid              
+    analysiscuts->SetMinPtCandidate(3.);
+    analysiscuts->SetMaxPtCandidate(10000.);
+    cout<<"This is the odject I'm going to save:"<<nptbins<<endl;
+    
+    analysiscuts->PrintAll();
+    TFile* fout=new TFile("DplustoKpipiCuts.root","recreate");   
+    fout->cd();
+    prodcuts->Write();
+    analysiscuts->Write();
+    fout->Close();
+    
+}
+
+
+
 void makeInputAliAnalysisTaskSESignificanceMaximization(){
   
   AliRDHFCutsDplustoKpipi* RDHFDplustoKpipi=new AliRDHFCutsDplustoKpipi();