]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pPb.C
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
- changed cut selections for all calo related tasks, changed definition of cuts for...
[u/mrichter/AliRoot.git] / PWGGA / GammaConv / macros / AddTask_GammaConvCalo_pPb.C
diff --git a/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pPb.C b/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pPb.C
index a3391cbde6277d6cc9ecfab4ac30130e38c44bd0..4a68a2e958308b3d04829c7e1eff538d997dd860 100644 (file)
--- a/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pPb.C
+++ b/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pPb.C
@@ -114,8 +114,8 @@ void AddTask_GammaConvCalo_pPb(  Int_t trainConfig = 1,  //change different set
        task->SetIsMC(isMC);
        // Cut Numbers to use in Analysis
        Int_t numberOfCuts = 2;
-       if (trainConfig==3 || trainConfig==5  ){ numberOfCuts =6;}
-       if (trainConfig==6 || trainConfig==7 || trainConfig==8 || trainConfig==9 ){ numberOfCuts =5;}
+       if (trainConfig==2 || trainConfig==4  ){ numberOfCuts =6;}
+       if (trainConfig==5 || trainConfig==6 || trainConfig==7 || trainConfig==8 || trainConfig==32 || trainConfig==33){ numberOfCuts =5;}
        
        TString *eventCutArray = new TString[numberOfCuts];
        TString *photonCutArray = new TString[numberOfCuts];
@@ -126,55 +126,68 @@ void AddTask_GammaConvCalo_pPb(  Int_t trainConfig = 1,  //change different set
        // 0 "ClusterType",  1 "EtaMin", 2 "EtaMax", 3 "PhiMin", 4 "PhiMax", 5 "DistanceToBadChannel", 6 "Timing", 7 "TrackMatching", 8 "ExoticCell",
        // 9 "MinEnergy", 10 "MinNCells", 11 "MinM02", 12 "MaxM02", 13 "MinM20", 14 "MaxM20", 15 "MaximumDispersion", 16 "NLM"
        
-       if (trainConfig == 1){ 
-               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000041022030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kINT7 // EMCAL clusters
-               eventCutArray[ 1] = "8005201"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000041022030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMC7 // EMCAL clusters
-       } else if (trainConfig == 2) {  
-               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000031022000000"; mesonCutArray[0] = "01525065000000"; //standart cut, kINT7 // PHOS clusters
-               eventCutArray[ 1] = "8006201"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000031022000000"; mesonCutArray[1] = "01525065000000"; //standard cut, kPHI7       // PHOS clusters        
-       } else if (trainConfig == 3){ 
-               eventCutArray[ 0] = "8008101"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000041022030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kEMCEGA based on INT7 // EMCAL clusters
-               eventCutArray[ 1] = "8008301"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000041022030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMCEG1 based on INT7 // EMCAL clusters
-               eventCutArray[ 2] = "8008501"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000041022030000"; mesonCutArray[2] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 // EMCAL clusters
-               eventCutArray[ 3] = "8009101"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000041022030000"; mesonCutArray[3] = "01525065000000"; //standard cut, kEMCEJE based on INT7 // EMCAL clusters
-               eventCutArray[ 4] = "8009301"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000041022030000"; mesonCutArray[4] = "01525065000000"; //standard cut, kEMCEJ1 based on INT7 // EMCAL clusters
-               eventCutArray[ 5] = "8009501"; photonCutArray[ 5] = "002092970028250400000"; clusterCutArray[5] = "10000041022030000"; mesonCutArray[5] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 // EMCAL clusters
-       } else if (trainConfig == 4){ 
-               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000041052030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kINT7 // EMCAL clusters
-               eventCutArray[ 1] = "8005201"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000041052030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMC7 // EMCAL clusters     
-       } else if (trainConfig == 5){ 
-               eventCutArray[ 0] = "8008101"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000041052030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kEMCEGA based on INT7 // EMCAL clusters
-               eventCutArray[ 1] = "8008301"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000041052030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMCEG1 based on INT7 // EMCAL clusters
-               eventCutArray[ 2] = "8008501"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000041052030000"; mesonCutArray[2] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 // EMCAL clusters
-               eventCutArray[ 3] = "8009101"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000041052030000"; mesonCutArray[3] = "01525065000000"; //standard cut, kEMCEJE based on INT7 // EMCAL clusters
-               eventCutArray[ 4] = "8009301"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000041052030000"; mesonCutArray[4] = "01525065000000"; //standard cut, kEMCEJ1 based on INT7 // EMCAL clusters
-               eventCutArray[ 5] = "8009501"; photonCutArray[ 5] = "002092970028250400000"; clusterCutArray[5] = "10000041052030000"; mesonCutArray[5] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 // EMCAL clusters             
-       } else if (trainConfig == 6){ //EMCAL
-               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000041052030000"; mesonCutArray[0] = "01525065000000"; 
-               eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000042052030000"; mesonCutArray[1] = "01525065000000"; 
-               eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000043052030000"; mesonCutArray[2] = "01525065000000"; 
-               eventCutArray[ 3] = "8000001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000044052030000"; mesonCutArray[3] = "01525065000000"; 
-               eventCutArray[ 4] = "8000001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000045052030000"; mesonCutArray[4] = "01525065000000"; 
-       } else if (trainConfig == 7) { //PHOS
-               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000041053200000"; mesonCutArray[0] = "01525065000000"; 
-               eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000042053200000"; mesonCutArray[1] = "01525065000000"; 
-               eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000043053200000"; mesonCutArray[2] = "01525065000000"; 
-               eventCutArray[ 3] = "8000001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000044053200000"; mesonCutArray[3] = "01525065000000"; 
-               eventCutArray[ 4] = "8000001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000045053200000"; mesonCutArray[4] = "01525065000000"; 
-       } else if (trainConfig == 8){ //EMCAL
-               eventCutArray[ 0] = "8000002"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000041052030000"; mesonCutArray[0] = "01525065000000"; 
-               eventCutArray[ 1] = "8000002"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000042052030000"; mesonCutArray[1] = "01525065000000"; 
-               eventCutArray[ 2] = "8000002"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000043052030000"; mesonCutArray[2] = "01525065000000"; 
-               eventCutArray[ 3] = "8000002"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000044052030000"; mesonCutArray[3] = "01525065000000"; 
-               eventCutArray[ 4] = "8000002"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000045052030000"; mesonCutArray[4] = "01525065000000"; 
-       } else if (trainConfig == 9) { //PHOS
-               eventCutArray[ 0] = "8000002"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000041053200000"; mesonCutArray[0] = "01525065000000"; 
-               eventCutArray[ 1] = "8000002"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000042053200000"; mesonCutArray[1] = "01525065000000"; 
-               eventCutArray[ 2] = "8000002"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000043053200000"; mesonCutArray[2] = "01525065000000"; 
-               eventCutArray[ 3] = "8000002"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000044053200000"; mesonCutArray[3] = "01525065000000"; 
-               eventCutArray[ 4] = "8000002"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000045053200000"; mesonCutArray[4] = "01525065000000"; 
-
+       //************************************************ EMCAL clusters **********************************************************
+       if (trainConfig == 1){ // min energy = 0.3 GeV/c
+               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000042022030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kINT7 
+               eventCutArray[ 1] = "8005201"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000042022030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMC7 
+       } else if (trainConfig == 2){  // min energy = 0.3 GeV/c
+               eventCutArray[ 0] = "8008101"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000042022030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kEMCEGA based on INT7 
+               eventCutArray[ 1] = "8008301"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000042022030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMCEG1 based on INT7 
+               eventCutArray[ 2] = "8008501"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000042022030000"; mesonCutArray[2] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 
+               eventCutArray[ 3] = "8009101"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000042022030000"; mesonCutArray[3] = "01525065000000"; //standard cut, kEMCEJE based on INT7 
+               eventCutArray[ 4] = "8009301"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000042022030000"; mesonCutArray[4] = "01525065000000"; //standard cut, kEMCEJ1 based on INT7 
+               eventCutArray[ 5] = "8009501"; photonCutArray[ 5] = "002092970028250400000"; clusterCutArray[5] = "10000042022030000"; mesonCutArray[5] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 
+       } else if (trainConfig == 3){ // min energy = 0.4 GeV/c
+               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kINT7 
+               eventCutArray[ 1] = "8005201"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMC7       
+       } else if (trainConfig == 4){ // min energy = 0.4 GeV/c
+               eventCutArray[ 0] = "8008101"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kEMCEGA based on INT7 
+               eventCutArray[ 1] = "8008301"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMCEG1 based on INT7 
+               eventCutArray[ 2] = "8008501"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000042032030000"; mesonCutArray[2] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 
+               eventCutArray[ 3] = "8009101"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000042032030000"; mesonCutArray[3] = "01525065000000"; //standard cut, kEMCEJE based on INT7 
+               eventCutArray[ 4] = "8009301"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000042032030000"; mesonCutArray[4] = "01525065000000"; //standard cut, kEMCEJ1 based on INT7 
+               eventCutArray[ 5] = "8009501"; photonCutArray[ 5] = "002092970028250400000"; clusterCutArray[5] = "10000042032030000"; mesonCutArray[5] = "01525065000000"; //standard cut, kEMCEG2 based on INT7
+       } else if (trainConfig == 5){ //EMCAL variation of track matching
+               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000041032030000"; mesonCutArray[0] = "01525065000000"; // 
+               eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01525065000000"; 
+               eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01525065000000"; 
+               eventCutArray[ 3] = "8000001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000044032030000"; mesonCutArray[3] = "01525065000000"; 
+               eventCutArray[ 4] = "8000001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000045032030000"; mesonCutArray[4] = "01525065000000"; 
+       } else if (trainConfig == 6){ //EMCAL added signal
+               eventCutArray[ 0] = "8000002"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000041032030000"; mesonCutArray[0] = "01525065000000"; 
+               eventCutArray[ 1] = "8000002"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01525065000000"; 
+               eventCutArray[ 2] = "8000002"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01525065000000"; 
+               eventCutArray[ 3] = "8000002"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000044032030000"; mesonCutArray[3] = "01525065000000"; 
+               eventCutArray[ 4] = "8000002"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000045032030000"; mesonCutArray[4] = "01525065000000"; 
+       } else if (trainConfig == 7){ //EMCAL minEnergy variation
+               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000042012030000"; mesonCutArray[0] = "01525065000000"; //0.2 GeV/c
+               eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000042022030000"; mesonCutArray[1] = "01525065000000"; //0.3 GeV/c
+               eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000042032030000"; mesonCutArray[2] = "01525065000000"; //0.4 GeV/c default
+               eventCutArray[ 3] = "8000001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000042042030000"; mesonCutArray[3] = "01525065000000"; //0.5 GeV/c
+               eventCutArray[ 4] = "8000001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000042052030000"; mesonCutArray[4] = "01525065000000"; //0.6 GeV/c
+       } else if (trainConfig == 8){ //EMCAL minNCells variation
+               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000042031030000"; mesonCutArray[0] = "01525065000000"; //n cells >= 1
+               eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000042033030000"; mesonCutArray[1] = "01525065000000"; //n cells >= 3
+               eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000042032000000"; mesonCutArray[2] = "01525065000000"; //no M02 cut
+               eventCutArray[ 3] = "8000001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10021042032030000"; mesonCutArray[3] = "01525065000000"; //only modules with TRD infront
+               eventCutArray[ 4] = "8000001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10013042032030000"; mesonCutArray[4] = "01525065000000"; //no modules with TRD infront               
                
+       //************************************************ PHOS clusters **********************************************************     
+       } else if (trainConfig == 31) { // min energy = 0.3 GeV/c
+               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000042033200000"; mesonCutArray[0] = "01525065000000"; //standart cut, kINT7 
+               eventCutArray[ 1] = "8006201"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[1] = "01525065000000"; //standard cut, kPHI7               
+       } else if (trainConfig == 32) { //PHOS
+               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000041033200000"; mesonCutArray[0] = "01525065000000"; 
+               eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[1] = "01525065000000"; 
+               eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[2] = "01525065000000"; 
+               eventCutArray[ 3] = "8000001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000044033200000"; mesonCutArray[3] = "01525065000000"; 
+               eventCutArray[ 4] = "8000001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000045033200000"; mesonCutArray[4] = "01525065000000"; 
+       } else if (trainConfig == 33) { //PHOS
+               eventCutArray[ 0] = "8000002"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000041033200000"; mesonCutArray[0] = "01525065000000"; 
+               eventCutArray[ 1] = "8000002"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[1] = "01525065000000"; 
+               eventCutArray[ 2] = "8000002"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[2] = "01525065000000"; 
+               eventCutArray[ 3] = "8000002"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000044033200000"; mesonCutArray[3] = "01525065000000"; 
+               eventCutArray[ 4] = "8000002"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000045033200000"; mesonCutArray[4] = "01525065000000"; 
        } else {
                Error(Form("GammaConvCalo_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
                return;