Centrality bins for 2011 PbPb run (Davide)

diff --git a/PWGHF/vertexingHF/macros/AddTaskCFVertexingHF3Prong.C b/PWGHF/vertexingHF/macros/AddTaskCFVertexingHF3Prong.C
index 34d6b47ef50da936d1f61affbc1cc30c0241c6fe..c7e082f343aee1761e6a52199e2693a30020db95 100644 (file)
--- a/PWGHF/vertexingHF/macros/AddTaskCFVertexingHF3Prong.C
+++ b/PWGHF/vertexingHF/macros/AddTaskCFVertexingHF3Prong.C
@@ -21,8 +21,10 @@ const Int_t    minITSClusters = 5;
 
 const Float_t centmin_0_10 = 0.;
 const Float_t centmax_0_10 = 10.;
 
 const Float_t centmin_0_10 = 0.;
 const Float_t centmax_0_10 = 10.;

-const Float_t centmin_10_100 = 10.;
-const Float_t centmax_10_100 = 100.;
+const Float_t centmin_10_60 = 10.;
+const Float_t centmax_10_60 = 60.;
+const Float_t centmin_60_100 = 60.;
+const Float_t centmax_60_100 = 100.;

 const Float_t centmax = 100.;
 const Float_t fakemin = -0.5;
 const Float_t fakemax = 2.5.;
 const Float_t centmax = 100.;
 const Float_t fakemin = -0.5;
 const Float_t fakemax = 2.5.;


@@ -152,9 +154,10 @@ AliCFTaskVertexingHF *AddTaskCFVertexingHF3Prong(const char* cutFile = "./Dplust
        const Int_t nbinpTpi2_4_8  = 4 ; //bins in ptpi2 from 4 to 8 GeV
        const Int_t nbinpTpi2_8_10  = 1 ; //bins in ptpi2 from 8 to 10 GeV
        const Int_t nbinzvtx  = 30 ; //bins in z vertex
        const Int_t nbinpTpi2_4_8  = 4 ; //bins in ptpi2 from 4 to 8 GeV
        const Int_t nbinpTpi2_8_10  = 1 ; //bins in ptpi2 from 8 to 10 GeV
        const Int_t nbinzvtx  = 30 ; //bins in z vertex

-       const Int_t nbincent = 11; //bins in centrality
-       const Int_t nbincent_0_10 = 2;  //bins in centrality between 0 and 10
-       const Int_t nbincent_10_100 = 9;  //bins in centrality between 10 and 100
+       const Int_t nbincent = 18; //bins in centrality
+       const Int_t nbincent_0_10 = 4;  //bins in centrality between 0 and 10
+       const Int_t nbincent_10_60 = 10;  //bins in centrality between 10 and 60
+       const Int_t nbincent_60_100 = 4;  //bins in centrality between 60 and 100

        const Int_t nbinfake = 3;  //bins in fake
        const Int_t nbinpointingXY = 50;  //bins in cosPointingAngleXY
        const Int_t nbinnormDecayLXY = 20;  //bins in NormDecayLengthXY
        const Int_t nbinfake = 3;  //bins in fake
        const Int_t nbinpointingXY = 50;  //bins in cosPointingAngleXY
        const Int_t nbinnormDecayLXY = 20;  //bins in NormDecayLengthXY


@@ -291,11 +294,16 @@ AliCFTaskVertexingHF *AddTaskCFVertexingHF3Prong(const char* cutFile = "./Dplust
        }
        
        // centrality
        }
        
        // centrality

+       

        for(Int_t i=0; i<=nbincent_0_10; i++) binLimcent[i]=(Double_t)centmin_0_10 + (centmax_0_10-centmin_0_10)/nbincent_0_10*(Double_t)i ; 
        for(Int_t i=0; i<=nbincent_0_10; i++) binLimcent[i]=(Double_t)centmin_0_10 + (centmax_0_10-centmin_0_10)/nbincent_0_10*(Double_t)i ; 

-       if (binLimcent[nbincent_0_10] != centmin_10_100)  {
-               Error("AliCFHeavyFlavourTaskMultiVarMultiStep","Calculated bin lim for cent - 1st range - differs from expected!\n");
+       if (binLimcent[nbincent_0_10] != centmin_10_60)  {
+         Error("AliCFHeavyFlavourTaskMultiVarMultiStep","Calculated bin lim for cent - 1st range - differs from expected!\n");

        }
        }

-       for(Int_t i=0; i<=nbincent_10_100; i++) binLimcent[i+nbincent_0_10]=(Double_t)centmin_10_100 + (centmax_10_100-centmin_10_100)/nbincent_10_100*(Double_t)i ; 
+       for(Int_t i=0; i<=nbincent_10_60; i++) binLimcent[i+nbincent_0_10]=(Double_t)centmin_10_60 + (centmax_10_60-centmin_10_60)/nbincent_10_60*(Double_t)i ;
+       if (binLimcent[nbincent_0_10+nbincent_10_60] != centmin_60_100)  {
+         Error("AliCFHeavyFlavourTaskMultiVarMultiStep","Calculated bin lim for cent - 2st range - differs from expected!\n");
+       }
+       for(Int_t i=0; i<=nbincent_60_100; i++) binLimcent[i+nbincent_10_60]=(Double_t)centmin_60_100 + (centmax_60_100-centmin_60_100)/nbincent_60_100*(Double_t)i ;

        
        // fake
        for(Int_t i=0; i<=nbinfake; i++) {
        
        // fake
        for(Int_t i=0; i<=nbinfake; i++) {

diff --git a/PWGHF/vertexingHF/macros/AddTaskCFVertexingHF3ProngDs.C b/PWGHF/vertexingHF/macros/AddTaskCFVertexingHF3ProngDs.C
index 1fda58fefe91f4146d83814c6b0a28a3ac565d52..3fd6ed1c54b8f0aefe663676d44e5b11b1ec08b9 100644 (file)
--- a/PWGHF/vertexingHF/macros/AddTaskCFVertexingHF3ProngDs.C
+++ b/PWGHF/vertexingHF/macros/AddTaskCFVertexingHF3ProngDs.C
@@ -21,8 +21,10 @@ const Int_t    minITSClusters = 5;
 
 const Float_t centmin_0_10 = 0.;
 const Float_t centmax_0_10 = 10.;
 
 const Float_t centmin_0_10 = 0.;
 const Float_t centmax_0_10 = 10.;

-const Float_t centmin_10_100 = 10.;
-const Float_t centmax_10_100 = 100.;
+const Float_t centmin_10_60 = 10.;
+const Float_t centmax_10_60 = 60.;
+const Float_t centmin_60_100 = 60.;
+const Float_t centmax_60_100 = 100.;

 const Float_t centmax = 100.;
 const Float_t fakemin = -0.5;
 const Float_t fakemax = 2.5.;
 const Float_t centmax = 100.;
 const Float_t fakemin = -0.5;
 const Float_t fakemax = 2.5.;


@@ -152,9 +154,10 @@ AliCFTaskVertexingHF *AddTaskCFVertexingHF3ProngDs(Int_t decayOption=AliCFVertex
        const Int_t nbinpT3_4_8  = 4 ; //bins in pt3 from 4 to 8 GeV
        const Int_t nbinpT3_8_10  = 1 ; //bins in pt3 from 8 to 10 GeV
        const Int_t nbinzvtx  = 30 ; //bins in z vertex
        const Int_t nbinpT3_4_8  = 4 ; //bins in pt3 from 4 to 8 GeV
        const Int_t nbinpT3_8_10  = 1 ; //bins in pt3 from 8 to 10 GeV
        const Int_t nbinzvtx  = 30 ; //bins in z vertex

-       const Int_t nbincent = 11; //bins in centrality
-       const Int_t nbincent_0_10 = 2;  //bins in centrality between 0 and 10
-       const Int_t nbincent_10_100 = 9;  //bins in centrality between 10 and 100
+       const Int_t nbincent = 18;  //bins in centrality
+       const Int_t nbincent_0_10 = 4;  //bins in centrality between 0 and 10
+       const Int_t nbincent_10_60 = 10;  //bins in centrality between 10 and 60
+       const Int_t nbincent_60_100 = 4;  //bins in centrality between 60 and 100

        const Int_t nbinfake = 3;  //bins in fake
        const Int_t nbinpointingXY = 50;  //bins in cosPointingAngleXY
        const Int_t nbinnormDecayLXY = 20;  //bins in NormDecayLengthXY
        const Int_t nbinfake = 3;  //bins in fake
        const Int_t nbinpointingXY = 50;  //bins in cosPointingAngleXY
        const Int_t nbinnormDecayLXY = 20;  //bins in NormDecayLengthXY


@@ -290,12 +293,16 @@ AliCFTaskVertexingHF *AddTaskCFVertexingHF3ProngDs(Int_t decayOption=AliCFVertex
                binLimzvtx[i]=(Double_t)zvtxmin  + (zvtxmax-zvtxmin)  /nbinzvtx*(Double_t)i ;
        }
        
                binLimzvtx[i]=(Double_t)zvtxmin  + (zvtxmax-zvtxmin)  /nbinzvtx*(Double_t)i ;
        }
        

-       // centrality
+       // centrality 

        for(Int_t i=0; i<=nbincent_0_10; i++) binLimcent[i]=(Double_t)centmin_0_10 + (centmax_0_10-centmin_0_10)/nbincent_0_10*(Double_t)i ; 
        for(Int_t i=0; i<=nbincent_0_10; i++) binLimcent[i]=(Double_t)centmin_0_10 + (centmax_0_10-centmin_0_10)/nbincent_0_10*(Double_t)i ; 

-       if (binLimcent[nbincent_0_10] != centmin_10_100)  {
-               Error("AliCFHeavyFlavourTaskMultiVarMultiStep","Calculated bin lim for cent - 1st range - differs from expected!\n");
+       if (binLimcent[nbincent_0_10] != centmin_10_60)  {
+         Error("AliCFHeavyFlavourTaskMultiVarMultiStep","Calculated bin lim for cent - 1st range - differs from expected!\n");

        }
        }

-       for(Int_t i=0; i<=nbincent_10_100; i++) binLimcent[i+nbincent_0_10]=(Double_t)centmin_10_100 + (centmax_10_100-centmin_10_100)/nbincent_10_100*(Double_t)i ; 
+       for(Int_t i=0; i<=nbincent_10_60; i++) binLimcent[i+nbincent_0_10]=(Double_t)centmin_10_60 + (centmax_10_60-centmin_10_60)/nbincent_10_60*(Double_t)i ;
+       if (binLimcent[nbincent_0_10+nbincent_10_60] != centmin_60_100)  {
+               Error("AliCFHeavyFlavourTaskMultiVarMultiStep","Calculated bin lim for cent - 2st range - differs from expected!\n");
+       }
+       for(Int_t i=0; i<=nbincent_60_100; i++) binLimcent[i+nbincent_10_60]=(Double_t)centmin_60_100 + (centmax_60_100-centmin_60_100)/nbincent_60_100*(Double_t)i ;

        
        // fake
        for(Int_t i=0; i<=nbinfake; i++) {
        
        // fake
        for(Int_t i=0; i<=nbinfake; i++) {

diff --git a/PWGHF/vertexingHF/macros/AddTaskCFVertexingHF3ProngLc.C b/PWGHF/vertexingHF/macros/AddTaskCFVertexingHF3ProngLc.C
index 7ec0325b8f52d0a901d1c08486d8d6e6e71e3f0c..5130c1bc89cd6a19927a491394d5dd6d5a68c2f5 100644 (file)
--- a/PWGHF/vertexingHF/macros/AddTaskCFVertexingHF3ProngLc.C
+++ b/PWGHF/vertexingHF/macros/AddTaskCFVertexingHF3ProngLc.C
@@ -21,8 +21,10 @@ const Int_t    minITSClusters = 5;
 
 const Float_t centmin_0_10 = 0.;
 const Float_t centmax_0_10 = 10.;
 
 const Float_t centmin_0_10 = 0.;
 const Float_t centmax_0_10 = 10.;

-const Float_t centmin_10_100 = 10.;
-const Float_t centmax_10_100 = 100.;
+const Float_t centmin_10_60 = 10.;
+const Float_t centmax_10_60 = 60.;
+const Float_t centmin_60_100 = 60.;
+const Float_t centmax_60_100 = 100.;

 const Float_t centmax = 100.;
 const Float_t fakemin = -0.5;
 const Float_t fakemax = 2.5.;
 const Float_t centmax = 100.;
 const Float_t fakemin = -0.5;
 const Float_t fakemax = 2.5.;


@@ -159,9 +161,10 @@ AliCFContainer *AddTaskCFVertexingHF3ProngLc(const char* cutFile = "./cuts4Lctop
        const Int_t nbinpTpi2_4_8  = 4 ; //bins in ptpi2 from 4 to 8 GeV
        const Int_t nbinpTpi2_8_10  = 1 ; //bins in ptpi2 from 8 to 10 GeV
        const Int_t nbinzvtx  = 30 ; //bins in z vertex
        const Int_t nbinpTpi2_4_8  = 4 ; //bins in ptpi2 from 4 to 8 GeV
        const Int_t nbinpTpi2_8_10  = 1 ; //bins in ptpi2 from 8 to 10 GeV
        const Int_t nbinzvtx  = 30 ; //bins in z vertex

-       const Int_t nbincent = 11; //bins in centrality
-       const Int_t nbincent_0_10 = 2;  //bins in centrality between 0 and 10
-       const Int_t nbincent_10_100 = 9;  //bins in centrality between 10 and 100
+       const Int_t nbincent = 18;  //bins in centrality
+       const Int_t nbincent_0_10 = 4;  //bins in centrality between 0 and 10
+       const Int_t nbincent_10_60 = 10;  //bins in centrality between 10 and 60
+       const Int_t nbincent_60_100 = 4;  //bins in centrality between 60 and 100

        const Int_t nbinfake = 3;  //bins in fake
        const Int_t nbindist12 = 10; //bins dist12
        const Int_t nbindist23 = 10; //bins dist23
        const Int_t nbinfake = 3;  //bins in fake
        const Int_t nbindist12 = 10; //bins dist12
        const Int_t nbindist23 = 10; //bins dist23


@@ -315,10 +318,14 @@ AliCFContainer *AddTaskCFVertexingHF3ProngLc(const char* cutFile = "./cuts4Lctop
        
        // centrality
        for(Int_t i=0; i<=nbincent_0_10; i++) binLimcent[i]=(Double_t)centmin_0_10 + (centmax_0_10-centmin_0_10)/nbincent_0_10*(Double_t)i ; 
        
        // centrality
        for(Int_t i=0; i<=nbincent_0_10; i++) binLimcent[i]=(Double_t)centmin_0_10 + (centmax_0_10-centmin_0_10)/nbincent_0_10*(Double_t)i ; 

-       if (binLimcent[nbincent_0_10] != centmin_10_100)  {
+       if (binLimcent[nbincent_0_10] != centmin_10_60)  {

                Error("AliCFHeavyFlavourTaskMultiVarMultiStep","Calculated bin lim for cent - 1st range - differs from expected!\n");
        }
                Error("AliCFHeavyFlavourTaskMultiVarMultiStep","Calculated bin lim for cent - 1st range - differs from expected!\n");
        }

-       for(Int_t i=0; i<=nbincent_10_100; i++) binLimcent[i+nbincent_0_10]=(Double_t)centmin_10_100 + (centmax_10_100-centmin_10_100)/nbincent_10_100*(Double_t)i ; 
+       for(Int_t i=0; i<=nbincent_10_60; i++) binLimcent[i+nbincent_0_10]=(Double_t)centmin_10_60 + (centmax_10_60-centmin_10_60)/nbincent_10_60*(Double_t)i ;
+       if (binLimcent[nbincent_0_10+nbincent_10_60] != centmin_60_100)  {
+               Error("AliCFHeavyFlavourTaskMultiVarMultiStep","Calculated bin lim for cent - 2st range - differs from expected!\n");
+       }
+       for(Int_t i=0; i<=nbincent_60_100; i++) binLimcent[i+nbincent_10_60]=(Double_t)centmin_60_100 + (centmax_60_100-centmin_60_100)/nbincent_60_100*(Double_t)i ;

        
        // fake
        for(Int_t i=0; i<=nbinfake; i++) {
        
        // fake
        for(Int_t i=0; i<=nbinfake; i++) {

diff --git a/PWGHF/vertexingHF/macros/AddTaskCFVertexingHFCascade.C b/PWGHF/vertexingHF/macros/AddTaskCFVertexingHFCascade.C
index 8d9e7698e70b8565be2ee8611316516bcaed93fe..be96316deeabb344db5d8727d7270780456e4b8b 100644 (file)
--- a/PWGHF/vertexingHF/macros/AddTaskCFVertexingHFCascade.C
+++ b/PWGHF/vertexingHF/macros/AddTaskCFVertexingHFCascade.C
@@ -27,8 +27,10 @@ const Int_t    minITSClusters = 5;
 
 const Float_t centmin_0_10 = 0.;
 const Float_t centmax_0_10 = 10.;
 
 const Float_t centmin_0_10 = 0.;
 const Float_t centmax_0_10 = 10.;

-const Float_t centmin_10_100 = 10.;
-const Float_t centmax_10_100 = 100.;
+const Float_t centmin_10_60 = 10.;
+const Float_t centmax_10_60 = 60.;
+const Float_t centmin_60_100 = 60.;
+const Float_t centmax_60_100 = 100.;

 const Float_t centmax = 100.;
 const Float_t fakemin = -0.5;
 const Float_t fakemax = 2.5;
 const Float_t centmax = 100.;
 const Float_t fakemin = -0.5;
 const Float_t fakemax = 2.5;


@@ -118,9 +120,10 @@ AliCFTaskVertexingHF *AddTaskCFVertexingHFCascade(const char* cutFile = "./DStar
        const Int_t nbinpointing  = 50 ; //bins in cosPointingAngle
        const Int_t nbinphi  = 18 ; //bins in Phi
        const Int_t nbinzvtx  = 30 ; //bins in z vertex
        const Int_t nbinpointing  = 50 ; //bins in cosPointingAngle
        const Int_t nbinphi  = 18 ; //bins in Phi
        const Int_t nbinzvtx  = 30 ; //bins in z vertex

-       const Int_t nbincent = 11;  //bins in centrality
-       const Int_t nbincent_0_10 = 2;  //bins in centrality between 0 and 10
-       const Int_t nbincent_10_100 = 9;  //bins in centrality between 10 and 100
+       const Int_t nbincent = 18; //bins in centrality
+       const Int_t nbincent_0_10 = 4;  //bins in centrality between 0 and 10
+       const Int_t nbincent_10_60 = 10;  //bins in centrality between 10 and 60
+       const Int_t nbincent_60_100 = 4;  //bins in centrality between 60 and 100

        const Int_t nbinfake = 3;  //bins in fake
        const Int_t nbinpointingXY = 50;  //bins in cosPointingAngleXY
        const Int_t nbinnormDecayLXY = 20;  //bins in NormDecayLengthXY
        const Int_t nbinfake = 3;  //bins in fake
        const Int_t nbinpointingXY = 50;  //bins in cosPointingAngleXY
        const Int_t nbinnormDecayLXY = 20;  //bins in NormDecayLengthXY


@@ -303,11 +306,15 @@ AliCFTaskVertexingHF *AddTaskCFVertexingHFCascade(const char* cutFile = "./DStar
 
        // centrality
        for(Int_t i=0; i<=nbincent_0_10; i++) binLimcent[i]=(Double_t)centmin_0_10 + (centmax_0_10-centmin_0_10)/nbincent_0_10*(Double_t)i ; 
 
        // centrality
        for(Int_t i=0; i<=nbincent_0_10; i++) binLimcent[i]=(Double_t)centmin_0_10 + (centmax_0_10-centmin_0_10)/nbincent_0_10*(Double_t)i ; 

-       if (binLimcent[nbincent_0_10] != centmin_10_100)  {
-               Error("AliCFHeavyFlavourTaskMultiVarMultiStep","Calculated bin lim for cent - 1st range - differs from expected!\n");
+       if (binLimcent[nbincent_0_10] != centmin_10_60)  {
+         Error("AliCFHeavyFlavourTaskMultiVarMultiStep","Calculated bin lim for cent - 1st range - differs from expected!\n");

        }
        }

-       for(Int_t i=0; i<=nbincent_10_100; i++) binLimcent[i+nbincent_0_10]=(Double_t)centmin_10_100 + (centmax_10_100-centmin_10_100)/nbincent_10_100*(Double_t)i ; 
-
+       for(Int_t i=0; i<=nbincent_10_60; i++) binLimcent[i+nbincent_0_10]=(Double_t)centmin_10_60 + (centmax_10_60-centmin_10_60)/nbincent_10_60*(Double_t)i ;
+       if (binLimcent[nbincent_0_10+nbincent_10_60] != centmin_60_100)  {
+         Error("AliCFHeavyFlavourTaskMultiVarMultiStep","Calculated bin lim for cent - 2st range - differs from expected!\n");
+       }
+       for(Int_t i=0; i<=nbincent_60_100; i++) binLimcent[i+nbincent_10_60]=(Double_t)centmin_60_100 + (centmax_60_100-centmin_60_100)/nbincent_60_100*(Double_t)i ;
+        

        // fake
        for(Int_t i=0; i<=nbinfake; i++) {
          binLimfake[i]=(Double_t)fakemin  + (fakemax-fakemin)/nbinfake * (Double_t)i;
        // fake
        for(Int_t i=0; i<=nbinfake; i++) {
          binLimfake[i]=(Double_t)fakemin  + (fakemax-fakemin)/nbinfake * (Double_t)i;