Update (Chiara)

diff --git a/PWG3/vertexingHF/macros/makeTFile4CutsD0toKpi.C b/PWG3/vertexingHF/macros/makeTFile4CutsD0toKpi.C
index 4df7536c54532c49da285af612a94a6203c66605..2c5a7e306a44e701b0c89ef4d52d43fd12b20578 100644 (file)
--- a/PWG3/vertexingHF/macros/makeTFile4CutsD0toKpi.C
+++ b/PWG3/vertexingHF/macros/makeTFile4CutsD0toKpi.C
@@ -28,8 +28,9 @@ void makeInputAliAnalysisTaskSED0Mass(){
   //default
   esdTrackCuts->SetRequireTPCRefit(kTRUE);
   esdTrackCuts->SetRequireITSRefit(kTRUE);
+  esdTrackCuts->SetEtaRange(-0.8,0.8);
   //esdTrackCuts->SetMinNClustersITS(4); // default is 5
-  //esdTrackCuts->SetMinNClustersTPC(70);
+  //esdTrackCuts->SetMinNClustersTPC(120);
   esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
                                         AliESDtrackCuts::kAny); 
  // default is kBoth, otherwise kAny
@@ -252,17 +253,18 @@ void makeInputAliAnalysisTaskSESignificanceMaximization(){
   esdTrackCuts->SetRequireTPCRefit(kTRUE);
   esdTrackCuts->SetRequireITSRefit(kTRUE);
   //esdTrackCuts->SetMinNClustersITS(4);
-  
+  //esdTrackCuts->SetMinNClustersTPC(120);
+
   esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
   esdTrackCuts->SetMinDCAToVertexXY(0.);
-  esdTrackCuts->SetEtaRange(-0.9,0.9);
+  esdTrackCuts->SetEtaRange(-0.8,0.8);
   esdTrackCuts->SetPtRange(0.1,1.e10);
   
   RDHFD0toKpi->AddTrackCuts(esdTrackCuts);
 
   const Int_t nvars=9;
 
-  const Int_t nptbins=10; //change this when adding pt bins!
+  const Int_t nptbins=14; //change this when adding pt bins!
   Float_t ptbins[nptbins+1];
   ptbins[0]=0.;
   ptbins[1]=0.5;
@@ -273,10 +275,12 @@ void makeInputAliAnalysisTaskSESignificanceMaximization(){
   ptbins[6]=5.;
   ptbins[7]=6.;
   ptbins[8]=8.;
-  ptbins[9]=12.;
-  //ptbins[10]=16.;
-  ptbins[10]=9999.;
-
+  ptbins[9]=10.;
+  ptbins[10]=12.;
+  ptbins[11]=14.;
+  ptbins[12]=16.;
+  ptbins[13]=24.;
+  ptbins[14]=9999.;
   
   RDHFD0toKpi->SetPtBins(nptbins+1,ptbins);
   
@@ -299,139 +303,28 @@ void makeInputAliAnalysisTaskSESignificanceMaximization(){
     //     printf("    d0d0  [cm^2] < %f\n",fD0toKpiCuts[7]);
     //     printf("    cosThetaPoint    > %f\n",fD0toKpiCuts[8]);
 
-  /*
-  //setting PPR cut values
-  rdcutsvalPPR[0][0]=0.7;
-  rdcutsvalPPR[1][0]=0.04;
-  rdcutsvalPPR[2][0]=0.8;
-  rdcutsvalPPR[3][0]=0.5;
-  rdcutsvalPPR[4][0]=0.5;
-  rdcutsvalPPR[5][0]=0.05;
-  rdcutsvalPPR[6][0]=0.05;
-  rdcutsvalPPR[7][0]=-0.0002;
-  rdcutsvalPPR[8][0]=0.5;
-
-  rdcutsvalPPR[0][1]=rdcutsvalPPR[0][2]=0.7;
-  rdcutsvalPPR[1][1]=rdcutsvalPPR[1][2]=0.02;
-  rdcutsvalPPR[2][1]=rdcutsvalPPR[2][2]=0.8;
-  rdcutsvalPPR[3][1]=rdcutsvalPPR[3][2]=0.7;
-  rdcutsvalPPR[4][1]=rdcutsvalPPR[4][2]=0.7;
-  rdcutsvalPPR[5][1]=rdcutsvalPPR[5][2]=0.05;
-  rdcutsvalPPR[6][1]=rdcutsvalPPR[6][2]=0.05;
-  rdcutsvalPPR[7][1]=rdcutsvalPPR[7][2]=-0.0002;
-  rdcutsvalPPR[8][1]=rdcutsvalPPR[8][2]=0.6;
-
-  rdcutsvalPPR[0][3]=0.7;
-  rdcutsvalPPR[1][3]=0.02;
-  rdcutsvalPPR[2][3]=0.8;
-  rdcutsvalPPR[3][3]=0.7;
-  rdcutsvalPPR[4][3]=0.7;
-  rdcutsvalPPR[5][3]=0.05;
-  rdcutsvalPPR[6][3]=0.05;
-  rdcutsvalPPR[7][3]=-0.0001;
-  rdcutsvalPPR[8][3]=0.8;
-
-  rdcutsvalPPR[0][4]=0.7;
-  rdcutsvalPPR[1][4]=0.02;
-  rdcutsvalPPR[2][4]=0.8;
-  rdcutsvalPPR[3][4]=0.7;
-  rdcutsvalPPR[4][4]=0.7;
-  rdcutsvalPPR[5][4]=0.05;
-  rdcutsvalPPR[6][4]=0.05;
-  rdcutsvalPPR[7][4]=-0.00005;
-  rdcutsvalPPR[8][4]=0.8;
-  */
-  /*
-  //setting my cut values
 
-  Float_t passcut[9]={0.7,0.04,0.8,0.7,0.7,0.1,0.1,-0.00015,0.5};
-
-  //setting my cut values
-  //0-1
-  rdcutsvalmine[0][0]=passcut[0];//0.7;
-  rdcutsvalmine[1][0]=passcut[1];//0.04;
-  rdcutsvalmine[2][0]=passcut[2];//0.8;
-  rdcutsvalmine[3][0]=passcut[3];//0.5;
-  rdcutsvalmine[4][0]=passcut[4];//0.5;
-  rdcutsvalmine[5][0]=passcut[5];//0.05;
-  rdcutsvalmine[6][0]=passcut[6];//0.05;
-  rdcutsvalmine[7][0]=passcut[7];//-0.00025;
-  rdcutsvalmine[8][0]=passcut[8];//0.7;
-
-  //1-2;2-3
-  passcut[1]=0.02;
-  passcut[8]=0.7;
-  rdcutsvalmine[0][1]=rdcutsvalmine[0][2]=passcut[0];//0.7;
-  rdcutsvalmine[1][1]=rdcutsvalmine[1][2]=passcut[1];//0.02;
-  rdcutsvalmine[2][1]=rdcutsvalmine[2][2]=passcut[2];//0.8;
-  rdcutsvalmine[3][1]=rdcutsvalmine[3][2]=passcut[3];//0.7;
-  rdcutsvalmine[4][1]=rdcutsvalmine[4][2]=passcut[4];//0.7;
-  rdcutsvalmine[5][1]=rdcutsvalmine[5][2]=passcut[5];//1.;
-  rdcutsvalmine[6][1]=rdcutsvalmine[6][2]=passcut[6];//1.;
-  rdcutsvalmine[7][1]=rdcutsvalmine[7][2]=passcut[7];//-0.00025;
-  rdcutsvalmine[8][1]=rdcutsvalmine[8][2]=passcut[8];//0.8;
-
-  //3-5
-  passcut[7]=-0.00005;
-  passcut[8]=0.6;
-  rdcutsvalmine[0][3]=passcut[0];//0.7;
-  rdcutsvalmine[1][3]=passcut[1];//0.02;
-  rdcutsvalmine[2][3]=passcut[2];//0.8;
-  rdcutsvalmine[3][3]=passcut[3];//0.7;
-  rdcutsvalmine[4][3]=passcut[4];//0.7;
-  rdcutsvalmine[5][3]=passcut[5];//0.05;
-  rdcutsvalmine[6][3]=passcut[6];//0.05;
-  rdcutsvalmine[7][3]=passcut[7];//-0.00015;
-  rdcutsvalmine[8][3]=passcut[8];//0.8;
-
-  //5-8
-  passcut[7]=-0.00001;
-  passcut[8]=0.6;
-  rdcutsvalmine[0][4]=passcut[0];//0.7;
-  rdcutsvalmine[1][4]=passcut[1];//0.02;
-  rdcutsvalmine[2][4]=passcut[2];//0.8;
-  rdcutsvalmine[3][4]=passcut[3];//0.7;
-  rdcutsvalmine[4][4]=passcut[4];//0.7;
-  rdcutsvalmine[5][4]=passcut[5];//0.05;
-  rdcutsvalmine[6][4]=passcut[6];//0.05;
-  rdcutsvalmine[7][4]=passcut[7];//-0.00015;
-  rdcutsvalmine[8][4]=passcut[8];//0.9;
-
-  //8-12
-  passcut[7]=-0.00000;
-  rdcutsvalmine[0][5]=passcut[0];
-  rdcutsvalmine[1][5]=passcut[1];
-  rdcutsvalmine[2][5]=passcut[2];
-  rdcutsvalmine[3][5]=passcut[3];
-  rdcutsvalmine[4][5]=passcut[4];
-  rdcutsvalmine[5][5]=passcut[5];
-  rdcutsvalmine[6][5]=passcut[6];
-  rdcutsvalmine[7][5]=passcut[7];
-  rdcutsvalmine[8][5]=passcut[8];
-
-  //>12
-  passcut[7]=-0.00000;
-  rdcutsvalmine[0][6]=passcut[0];
-  rdcutsvalmine[1][6]=passcut[1];
-  rdcutsvalmine[2][6]=passcut[2];
-  rdcutsvalmine[3][6]=passcut[3];
-  rdcutsvalmine[4][6]=passcut[4];
-  rdcutsvalmine[5][6]=passcut[5];
-  rdcutsvalmine[6][6]=passcut[6];
-  rdcutsvalmine[7][6]=passcut[7];
-  rdcutsvalmine[8][6]=passcut[8];
-  */
-
-  Float_t cutsMatrixD0toKpiStand[nptbins][nvars]={{0.2,400.*1E-4,0.7,0.3,0.3,1000.*1E-4,1000.*1E-4,-0.00005,0.85},/* pt<0.5*/
-                                                 {0.2,400.*1E-4,0.7,0.3,0.3,1000.*1E-4,1000.*1E-4,-0.00005,0.85},/* 0.5<pt<1*/ 
-                                                 {0.2,400.*1E-4,0.7,0.4,0.4,1000.*1E-4,1000.*1E-4,-0.00005,0.85},/* 1<pt<2 */
-                                                 {0.2,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,-0.00005,0.7},/* 2<pt<3 */
-                                                 {0.2,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,-0.00005,0.7},/* 3<pt<4 */
-                                                 {0.2,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,-0.00005,0.7},/* 4<pt<5 */
-                                                 {0.2,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,-0.00005,0.7},/* 5<pt<6 */
-                                                 {0.2,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,-0.00005,0.7},/* 6<pt<8 */
-                                                 {0.2,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.0002,0.7},/* 8<pt<12 */ 
-                                                 {0.2,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.0002,0.7}};/*pt>12 */ 
+  
+  Float_t cutsMatrixD0toKpiStand[nptbins][nvars]={{0.2,400.*1E-4,0.7,0.3,0.3,1000.*1E-4,1000.*1E-4,-0.00005,0.8},/* pt<0.5*/
+                                                 {0.2,400.*1E-4,0.7,0.3,0.3,1000.*1E-4,1000.*1E-4,-0.00005,0.8},/* 0.5<pt<1*/ 
+                                                 {0.2,400.*1E-4,0.7,0.4,0.4,1000.*1E-4,1000.*1E-4,-0.00005,0.8},/* 1<pt<2 */
+                                                 {0.2,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,-0.00005,0.7},/* 2<pt<3 */
+                                                 {0.2,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,-0.00005,0.7},/* 3<pt<4 */
+                                                 {0.2,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,-0.00005,0.7},/* 4<pt<5 */
+                                                 {0.2,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,-0.00005,0.7},/* 5<pt<6 */
+                                                 {0.3,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.00005,0.6},/* 6<pt<8 */
+                                                 {0.3,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.00005,0.6},/* 8<pt<10 */ 
+                                                 {0.3,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.00005,0.6}, /*10<pt<12*/
+                                                 {0.3,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.00005,0.6}, /*12<pt<14*/
+                                                 {0.3,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.00005,0.6}, /*14<pt<16*/
+                                                 {0.3,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.00005,0.6}, /*16<pt<24*/
+                                                 {0.3,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.00005,0.6}}; /*pt>24*/
+
+  // Float_t cutsMatrixD0toKpiStand[nptbins][nvars]={{0.3,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.0002,0.6},/* 6<pt<8 */
+  //                                             {0.3,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.0002,0.6},/* 8<pt<12 */ 
+  //                                             {0.3,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.0002,0.6}, /*12<pt<16*/
+  //                                             {0.3,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.0002,0.6}, /*16<pt<24*/
+  //                                             {0.3,400.*1E-4,0.7,0.7,0.7,1000.*1E-4,1000.*1E-4,0.0002,0.6}};/*pt>24 */ 
 
   
   //CREATE TRANSPOSE MATRIX...REVERSE INDICES as required by AliRDHFCuts
@@ -446,7 +339,7 @@ void makeInputAliAnalysisTaskSESignificanceMaximization(){
 
 
   Int_t nvarsforopt=RDHFD0toKpi->GetNVarsForOpt();
-  Int_t dim=3; //set this!!
+  Int_t dim=2; //set this!!
   Bool_t *boolforopt;
   boolforopt=new Bool_t[nvars];
   if(dim>nvarsforopt){
@@ -480,64 +373,85 @@ void makeInputAliAnalysisTaskSESignificanceMaximization(){
 
 
   Float_t tighterval[dim][nptbins];
-  //dca  <- this
+  //dca  
   //costhetastar
   //d0d0 <-this 
   //costhetapoint <-this 
 
   
-  //number of steps for each variable is 4 now
+  //number of steps for each variable is set in the AddTask arguments (default=8)
   //set this!!
   // tighterval[0][0]=0.01;
-  tighterval[0][0]=100e-4;
-  tighterval[1][0]=-0.0006;
-  tighterval[2][0]=1.;
+  // tighterval[0][0]=100e-4;
+  // tighterval[1][0]=-0.0006;
+  // tighterval[2][0]=1.;
+  tighterval[0][0]=-0.00065;
+  tighterval[1][0]=1.;
 
   // tighterval[0][1]=0.01;
-  tighterval[0][1]=100e-4;
-  tighterval[1][1]=-0.0007; //try with tighter in this bin
-  tighterval[2][1]=1.;
- 
+  // tighterval[0][1]=100e-4;
+  // tighterval[1][1]=-0.0007; //try with tighter in this bin
+  // tighterval[2][1]=1.;
+  tighterval[0][1]=-0.00065;
+  tighterval[1][1]=1.;
+
  // tighterval[0][2]=0.01;
-  tighterval[0][2]=100e-4;
-  tighterval[1][2]=-0.0006;
-  tighterval[2][2]=1.;
+  // tighterval[0][2]=100e-4;
+  // tighterval[1][2]=-0.0006;
+  // tighterval[2][2]=1.;
+  tighterval[0][2]=-0.00065;
+  tighterval[1][2]=1.;
  
   // tighterval[0][3]=0.01;
-  tighterval[0][3]=100e-4;
-  tighterval[1][3]=-0.0006;
-  tighterval[2][3]=0.95;
+  // tighterval[0][3]=100e-4;
+  // tighterval[1][3]=-0.0006;
+  // tighterval[2][3]=0.95;
+  tighterval[0][3]=-0.00065;
+  tighterval[1][3]=1.;
 
   // tighterval[0][4]=0.01;
-  tighterval[0][4]=100e-4;
-  tighterval[1][4]=-0.0006;
-  tighterval[2][4]=0.95;
-
+  // tighterval[0][4]=100e-4;
+  // tighterval[1][4]=-0.0006;
+  // tighterval[2][4]=0.95;
+  tighterval[0][4]=-0.00065;
+  tighterval[1][4]=1.;
+ 
   // tighterval[0][5]=0.01;
-  tighterval[0][5]=100e-4;
-  tighterval[1][5]=-0.0006;
-  tighterval[2][5]=0.95;
-
+  //  tighterval[0][5]=100e-4;
+  tighterval[0][5]=-0.00065;
+  tighterval[1][5]=1.;
+ 
   // tighterval[0][6]=0.01;
-  tighterval[0][6]=100e-4;
-  tighterval[1][6]=-0.0006;
-  tighterval[2][6]=0.95;
+  //  tighterval[0][6]=100e-4;
+  tighterval[0][6]=-0.00065;
+  tighterval[1][6]=1.;
 
   // tighterval[0][6]=0.01;
-  tighterval[0][7]=100e-4;
-  tighterval[1][7]=-0.0006;
-  tighterval[2][7]=0.95;
+  //tighterval[0][7]=100e-4;
+  tighterval[0][7]=-0.00027;
+  tighterval[1][7]=1.;
 
   // tighterval[0][6]=0.01;
-  tighterval[0][8]=100e-4;
-  tighterval[1][8]=-0.0006;
-  tighterval[2][8]=0.95;
+  //tighterval[0][8]=100e-4;
+  tighterval[0][8]=-0.00027;
+  tighterval[1][8]=1.;
 
   // tighterval[0][6]=0.01;
-  tighterval[0][9]=100e-4;
-  tighterval[1][9]=-0.0006;
-  tighterval[2][9]=0.95;
+  //tighterval[0][9]=100e-4;
+  tighterval[0][9]=-0.00027;
+  tighterval[1][9]=1.;
+ 
+  tighterval[0][10]=-0.00027;
+  tighterval[1][10]=1.;
 
+  tighterval[0][11]=-0.00027;
+  tighterval[1][11]=1.;
+
+  tighterval[0][12]=-0.00027;
+  tighterval[1][12]=1.;
+
+  tighterval[0][13]=-0.00027;
+  tighterval[1][13]=1.;
 
   TString name=""; 
   Int_t arrdim=dim*nptbins;
@@ -580,7 +494,8 @@ void makeInputAliAnalysisTaskSESignificanceMaximization(){
   RDHFD0toKpi->SetOptPileup(AliRDHFCuts::kRejectPileupEvent);
 
 
-  TFile* fout=new TFile("cuts4SignifMaxim10binspileup.root","recreate");   //set this!! 
+  TFile* fout=new TFile("cuts4SignifMaxim.root","recreate");   //set this!! 
+
   fout->cd();
   RDHFD0toKpi->Write();
   max.Write();