set cuts for conversion + calo to be restricted to the calorimeters as option
authorfbock <friederike.bock@cern.ch>
Fri, 5 Dec 2014 14:41:57 +0000 (15:41 +0100)
committerfbock <friederike.bock@cern.ch>
Fri, 5 Dec 2014 14:46:40 +0000 (15:46 +0100)
PWGGA/GammaConv/AliConversionPhotonCuts.cxx
PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pPb.C
PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pp.C

index 5bddb7e..4ee0b1c 100644 (file)
@@ -1897,7 +1897,7 @@ Bool_t AliConversionPhotonCuts::SetMinPhiSectorCut(Int_t minPhiCut) {
                if (!fDoShrinkTPCAcceptance) fDoShrinkTPCAcceptance = kTRUE;
                fMinPhiCut = 1.0; //PHOS
                break;
-       case 5:
+       case 4:
                if (!fDoShrinkTPCAcceptance) fDoShrinkTPCAcceptance = kTRUE;
                fMinPhiCut = 3.4; //EMCal tight
                break;
index 039781a..a5a9ca2 100644 (file)
@@ -171,6 +171,9 @@ void AddTask_GammaConvCalo_pPb(  Int_t trainConfig = 1,  //change different set
                eventCutArray[ 2] = "8000001"; photonCutArray[ 2] = "002000093270028250400000"; clusterCutArray[2] = "10000042032000000"; mesonCutArray[2] = "01631031000000"; //no M02 cut
                eventCutArray[ 3] = "8000001"; photonCutArray[ 3] = "002000093270028250400000"; clusterCutArray[3] = "10021042032030000"; mesonCutArray[3] = "01631031000000"; //only modules with TRD infront
                eventCutArray[ 4] = "8000001"; photonCutArray[ 4] = "002000093270028250400000"; clusterCutArray[4] = "10013042032030000"; mesonCutArray[4] = "01631031000000"; //no modules with TRD infront            
+       } else if (trainConfig == 9){ //PCM restriction in acceptance 
+               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002022093270028250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031000000"; // PCM photons pointing to EMCAL loose
+               eventCutArray[ 1] = "8000001"; photonCutArray[ 1] = "002044093270028250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; // PCM photons pointing to EMCAL tight
                
        //************************************************ PHOS clusters **********************************************************     
        } else if (trainConfig == 31) { // min energy = 0.3 GeV/c
index c962d08..feb63e4 100644 (file)
@@ -171,7 +171,9 @@ void AddTask_GammaConvCalo_pp(  Int_t trainConfig = 1,  //change different set o
                eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002000083270008250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; // min TPC cl > 0.35
                eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002000093270008250400000"; clusterCutArray[2] = "10000042032030000"; mesonCutArray[2] = "01631061000000"; // alpha < 0.8
                eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002000093270008250400000"; clusterCutArray[3] = "10000042032030000"; mesonCutArray[3] = "01631051000000"; // alpha < 0.75
-
+       } else if (trainConfig == 9){ // PCM variations
+               eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002022093270008250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031000000"; // restrict acceptance to EMCAL loose
+               eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002044093270008250400000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031000000"; // restrict acceptance to EMCAL tight
                // LHC13g       
        } else if (trainConfig == 10){  // EMCAL clusters, EMCEGA triggers, track matching 0.035
                eventCutArray[ 0] = "0008311"; photonCutArray[ 0] = "002000093270008250400000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031000000"; // EMCEG1,