TENDER becomes Tender

[u/mrichter/AliRoot.git] / PWGGA / GammaConv / macros / AddTask_GammaConvV1_PbPb.C

diff --git a/PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb.C b/PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb.C
index b92843be5a33499496acca582a70f974c28caee6..69a2d4cbe3848b4e0779a937a108faafbec91dae 100644 (file)
--- a/PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb.C
+++ b/PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb.C
@@ -24,8 +24,8 @@ void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set o
        gSystem->Load("libCDB.so");
        gSystem->Load("libSTEER.so");
        gSystem->Load("libSTEERBase.so");
-       gSystem->Load("libTENDER.so");
-       gSystem->Load("libTENDERSupplies.so");
+       gSystem->Load("libTender.so");
+       gSystem->Load("libTenderSupplies.so");
        gSystem->Load("libPWGflowBase.so");
        gSystem->Load("libPWGflowTasks.so");
        gSystem->Load("libPWGGAGammaConv.so");
@@ -110,7 +110,8 @@ void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set o
        task->SetIsHeavyIon(isHeavyIon);
        task->SetIsMC(isMC);
        // Cut Numbers to use in Analysis
-       Int_t numberOfCuts = 5;
+       if (trainConfig == 135 || trainConfig == 136 || trainConfig == 137 ) Int_t numberOfCuts = 7;
+        else Int_t numberOfCuts = 5; 
 
        TString *eventCutArray = new TString[numberOfCuts];
        TString *photonCutArray = new TString[numberOfCuts];
@@ -680,7 +681,7 @@ void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set o
                eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002000092970022250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
                eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002000092970022250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
                eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002000092970022250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
-       } else if ( trainConfig         == 95){ // qt 0.07 no2D  added signal
+       } else if ( trainConfig == 95){ // qt 0.07 no2D  added signal
                eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002000092970022250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
                eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002000092970022250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
                eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002000092970022250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
@@ -759,11 +760,11 @@ void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set o
                eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002000092970028250400000"; mesonCutArray[ 3]= "01525075000000"; // 20-40%
                eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002000092970028250400000"; mesonCutArray[ 4]= "01525075000000"; // 20-50%
        } else if ( trainConfig == 108){ // psi pair 0.2
-               eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002000092970028280400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
-               eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002000092970028280400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
-               eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002000092970028280400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
-               eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002000092970028280400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
-               eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002000092970028280400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
+               eventCutArray[ 0] = "6013301"; photonCutArray[ 0] = "002111092970028280400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
+               eventCutArray[ 1] = "6123301"; photonCutArray[ 1] = "002111092970028280400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
+               eventCutArray[ 2] = "5013301"; photonCutArray[ 2] = "002111092970028280400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
+               eventCutArray[ 3] = "5243601"; photonCutArray[ 3] = "002111092970028280400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
+               eventCutArray[ 4] = "5253601"; photonCutArray[ 4] = "002111092970028280400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
        } else if ( trainConfig == 109){ // psi pair 0.2  added signal
                eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002000092970028280400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
                eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002000092970028280400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
@@ -879,12 +880,12 @@ void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set o
                eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002000092370028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
                eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002000092370028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
        } else if ( trainConfig == 128){ // cleaner cuts, pion line at 2.0 sigma
-               eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002000092570028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
-               eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002000092570028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
-               eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002000092570028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
-               eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002000092570028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
-               eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002000092570028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
-       } else if ( trainConfig == 129){ // cleaner cuts, pion line at 2.0 sigma added signals
+               eventCutArray[ 0] = "6013301"; photonCutArray[ 0] = "002111092573028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               eventCutArray[ 1] = "6123301"; photonCutArray[ 1] = "002111092573028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               eventCutArray[ 2] = "5013301"; photonCutArray[ 2] = "002111092573028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
+               eventCutArray[ 3] = "5243601"; photonCutArray[ 3] = "002111092573028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
+               eventCutArray[ 4] = "5253601"; photonCutArray[ 4] = "002111092573028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
+       } else if ( trainConfig == 129){ // cleaner cuts, pion line at 2.0 sigma added signals
                eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002000092570028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
                eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002000092570028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
                eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002000092570028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
@@ -908,12 +909,48 @@ void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set o
                eventCutArray[ 2] = "5340001"; photonCutArray[ 2] = "002000092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
                eventCutArray[ 3] = "5450001"; photonCutArray[ 3] = "002000092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
                eventCutArray[ 4] = "5560001"; photonCutArray[ 4] = "002000092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%  
-       } else if ( trainConfig == 133){ // cleaner cuts
-               eventCutArray[ 0] = "5080001"; photonCutArray[ 0] = "002000092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // no phi cut
-               eventCutArray[ 1] = "5080001"; photonCutArray[ 1] = "002111092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // oroc phi cut
-               eventCutArray[ 2] = "5080001"; photonCutArray[ 2] = "002022092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // EMCal phi cut
-               eventCutArray[ 3] = "5080001"; photonCutArray[ 3] = "002033092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // PHOS phi cut
-               eventCutArray[ 4] = "5080001"; photonCutArray[ 4] = "002223092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // inclusion phi
+       } else if ( trainConfig == 133){ // cleaner cuts with oroc phi cut & specific centrality selection
+               eventCutArray[ 0] = "6013301"; photonCutArray[ 0] = "002111092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
+               eventCutArray[ 1] = "6123301"; photonCutArray[ 1] = "002111092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
+               eventCutArray[ 2] = "5013301"; photonCutArray[ 2] = "002111092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% central
+               eventCutArray[ 3] = "5013601"; photonCutArray[ 3] = "002111092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 0-10% semicentral
+               eventCutArray[ 4] = "5123601"; photonCutArray[ 4] = "002111092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 10-20%
+       } else if ( trainConfig == 134){ // cleaner cuts
+               eventCutArray[ 0] = "5243601"; photonCutArray[ 0] = "002111092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 20-40%
+               eventCutArray[ 1] = "5253601"; photonCutArray[ 1] = "002111092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 20-50%
+               eventCutArray[ 2] = "5233601"; photonCutArray[ 2] = "002111092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 20-30%
+               eventCutArray[ 3] = "5343601"; photonCutArray[ 3] = "002111092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 30-40%
+               eventCutArray[ 4] = "5453601"; photonCutArray[ 4] = "002111092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 40-50%
+       } else if ( trainConfig == 135){ // flow cuts with eta = 0.9, y = 0.85
+               eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002000092970028250400000"; mesonCutArray[ 0] = "01525065000000";
+               eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002000092970028250400000"; mesonCutArray[ 1] = "01525065000000";
+               eventCutArray[ 2] = "5120001"; photonCutArray[ 2] = "002000092970028250400000"; mesonCutArray[ 2] = "01525065000000";
+               eventCutArray[ 3] = "5230001"; photonCutArray[ 3] = "002000092970028250400000"; mesonCutArray[ 3] = "01525065000000";
+               eventCutArray[ 4] = "5340001"; photonCutArray[ 4] = "002000092970028250400000"; mesonCutArray[ 4] = "01525065000000";
+               eventCutArray[ 5] = "5460001"; photonCutArray[ 5] = "002000092970028250400000"; mesonCutArray[ 5] = "01525065000000";
+               eventCutArray[ 6] = "5680001"; photonCutArray[ 6] = "002000092970028250400000"; mesonCutArray[ 6] = "01525065000000";
+       } else if ( trainConfig == 136){ // flow cuts with eta = 0.65, y = 0.6
+               eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "032000092970028250400000"; mesonCutArray[ 0] = "01523065000000";
+               eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "032000092970028250400000"; mesonCutArray[ 1] = "01523065000000";
+               eventCutArray[ 2] = "5120001"; photonCutArray[ 2] = "032000092970028250400000"; mesonCutArray[ 2] = "01523065000000";
+               eventCutArray[ 3] = "5230001"; photonCutArray[ 3] = "032000092970028250400000"; mesonCutArray[ 3] = "01523065000000";
+               eventCutArray[ 4] = "5340001"; photonCutArray[ 4] = "032000092970028250400000"; mesonCutArray[ 4] = "01523065000000";
+               eventCutArray[ 5] = "5460001"; photonCutArray[ 5] = "032000092970028250400000"; mesonCutArray[ 5] = "01523065000000";
+               eventCutArray[ 6] = "5680001"; photonCutArray[ 6] = "032000092970028250400000"; mesonCutArray[ 6] = "01523065000000";
+       } else if ( trainConfig == 137){ // flow cuts with eta = 0.6, y = 0.5
+               eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "012000092970028250400000"; mesonCutArray[ 0] = "01524065000000";
+               eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "012000092970028250400000"; mesonCutArray[ 1] = "01524065000000";
+               eventCutArray[ 2] = "5120001"; photonCutArray[ 2] = "012000092970028250400000"; mesonCutArray[ 2] = "01524065000000";
+               eventCutArray[ 3] = "5230001"; photonCutArray[ 3] = "012000092970028250400000"; mesonCutArray[ 3] = "01524065000000";
+               eventCutArray[ 4] = "5340001"; photonCutArray[ 4] = "012000092970028250400000"; mesonCutArray[ 4] = "01524065000000";
+               eventCutArray[ 5] = "5460001"; photonCutArray[ 5] = "012000092970028250400000"; mesonCutArray[ 5] = "01524065000000";
+               eventCutArray[ 6] = "5680001"; photonCutArray[ 6] = "012000092970028250400000"; mesonCutArray[ 6] = "01524065000000";  
+       } else if ( trainConfig == 138){
+               eventCutArray[ 0] = "6013301"; photonCutArray[ 0] = "002111092573028280400000"; mesonCutArray[ 0] = "01525065000000"; 
+               eventCutArray[ 1] = "6123301"; photonCutArray[ 1] = "002111092573028280400000"; mesonCutArray[ 1] = "01525065000000"; 
+               eventCutArray[ 2] = "5013301"; photonCutArray[ 2] = "002111092573028280400000"; mesonCutArray[ 2] = "01525065000000"; 
+               eventCutArray[ 3] = "5243601"; photonCutArray[ 3] = "002111092573028280400000"; mesonCutArray[ 3] = "01525065000000"; 
+               eventCutArray[ 4] = "5243601"; photonCutArray[ 4] = "002111092573028280400000"; mesonCutArray[ 4] = "01525065000000"; 
        } else {
                Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
                return;