From: fbock Date: Tue, 2 Dec 2014 11:31:44 +0000 (+0100) Subject: - changed cut selections for all calo related tasks, changed definition of cuts for... X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=0e006b2c90fca0b8e1fd223df4fd7fa6acbf08d2;p=u%2Fmrichter%2FAliRoot.git - changed cut selections for all calo related tasks, changed definition of cuts for min energy, and definition of phi cuts to follow ALICE convention reaching from 0-2Pi --- diff --git a/PWGGA/GammaConv/AliCaloPhotonCuts.cxx b/PWGGA/GammaConv/AliCaloPhotonCuts.cxx index 4eade9df878..327fb92b61c 100644 --- a/PWGGA/GammaConv/AliCaloPhotonCuts.cxx +++ b/PWGGA/GammaConv/AliCaloPhotonCuts.cxx @@ -320,11 +320,11 @@ void AliCaloPhotonCuts::InitCutHistograms(TString name){ fHistograms->Add(fHistClusterIdentificationCuts); // Acceptance related histogramms - fHistClusterEtavsPhiBeforeAcc=new TH2F(Form("EtaPhi_beforeAcceptance %s",GetCutNumber().Data()),"EtaPhi_beforeAcceptance",430,-TMath::Pi(),TMath::Pi(),96,-0.7,0.7); // 462->430, 110->96 + fHistClusterEtavsPhiBeforeAcc=new TH2F(Form("EtaPhi_beforeAcceptance %s",GetCutNumber().Data()),"EtaPhi_beforeAcceptance",430,0,2*TMath::Pi(),96,-0.7,0.7); // 462->430, 110->96 fHistograms->Add(fHistClusterEtavsPhiBeforeAcc); - fHistClusterEtavsPhiAfterAcc=new TH2F(Form("EtaPhi_afterAcceptance %s",GetCutNumber().Data()),"EtaPhi_afterAcceptance",430,-TMath::Pi(),TMath::Pi(),96,-0.7,0.7); // 462->430, 110->96 + fHistClusterEtavsPhiAfterAcc=new TH2F(Form("EtaPhi_afterAcceptance %s",GetCutNumber().Data()),"EtaPhi_afterAcceptance",430,0,2*TMath::Pi(),96,-0.7,0.7); // 462->430, 110->96 fHistograms->Add(fHistClusterEtavsPhiAfterAcc); - fHistClusterEtavsPhiAfterQA=new TH2F(Form("EtaPhi_afterClusterQA %s",GetCutNumber().Data()),"EtaPhi_afterClusterQA",430,-TMath::Pi(),TMath::Pi(),96,-0.7,0.7); // 462->430, 110->96 + fHistClusterEtavsPhiAfterQA=new TH2F(Form("EtaPhi_afterClusterQA %s",GetCutNumber().Data()),"EtaPhi_afterClusterQA",430,0,2*TMath::Pi(),96,-0.7,0.7); // 462->430, 110->96 fHistograms->Add(fHistClusterEtavsPhiAfterQA); //fHistDistanceToBadChannelBeforeAcc = new TH1F(Form("DistanceToBadChannel_beforeAcceptance %s",GetCutNumber().Data()),"DistanceToBadChannel_beforeAcceptance",200,0,40); //fHistograms->Add(fHistDistanceToBadChannelBeforeAcc); @@ -649,7 +649,8 @@ Bool_t AliCaloPhotonCuts::ClusterQualityCuts(AliVCluster* cluster, AliVEvent *ev TVector3 clusterVector(clusPos[0],clusPos[1],clusPos[2]); Double_t etaCluster = clusterVector.Eta(); Double_t phiCluster = clusterVector.Phi(); - + if (phiCluster < 0) phiCluster= phiCluster+ 2*TMath::Pi(); + if(fHistClusterEtavsPhiAfterQA) fHistClusterEtavsPhiAfterQA->Fill(phiCluster,etaCluster); if(fHistClusterTimevsEAfterQA) fHistClusterTimevsEAfterQA->Fill(cluster->GetTOF(), cluster->E()); // if(fHistExoticCellAfterQA) fHistExoticCellAfterQA->Fill(cluster->E(), ); @@ -684,6 +685,7 @@ Bool_t AliCaloPhotonCuts::ClusterIsSelected(AliVCluster *cluster, AliVEvent * ev TVector3 clusterVector(clusPos[0],clusPos[1],clusPos[2]); Double_t etaCluster = clusterVector.Eta(); Double_t phiCluster = clusterVector.Phi(); + if (phiCluster < 0) phiCluster= phiCluster+ 2*TMath::Pi(); // Histos before cuts if(fHistClusterEtavsPhiBeforeAcc) fHistClusterEtavsPhiBeforeAcc->Fill(phiCluster,etaCluster); @@ -741,6 +743,7 @@ Bool_t AliCaloPhotonCuts::AcceptanceCuts(AliVCluster *cluster, AliVEvent* event) TVector3 clusterVector(clusPos[0],clusPos[1],clusPos[2]); Double_t etaCluster = clusterVector.Eta(); Double_t phiCluster = clusterVector.Phi(); + if (phiCluster < 0) phiCluster= phiCluster+ 2*TMath::Pi(); // check eta range if (fUseEtaCut){ @@ -1226,10 +1229,20 @@ Bool_t AliCaloPhotonCuts::SetMinPhiCut(Int_t minPhi) if (!fUsePhiCut) fUsePhiCut=0; fMinPhiCut=-10000; break; - case 1: + case 1: // min EMCAL if (!fUsePhiCut) fUsePhiCut=1; fMinPhiCut=1.39626; break; + case 2: // min EMCAL with TRD 2012 + if (!fUsePhiCut) fUsePhiCut=1; + fMinPhiCut=2.10; + break; + case 3: // min EMCAL with TRD 2011 + if (!fUsePhiCut) fUsePhiCut=1; + fMinPhiCut=2.45; + break; + + default: AliError(Form("MinPhi Cut not defined %d",minPhi)); return kFALSE; @@ -1245,10 +1258,20 @@ Bool_t AliCaloPhotonCuts::SetMaxPhiCut(Int_t maxPhi) if (!fUsePhiCut) fUsePhiCut=0; fMaxPhiCut=10000; break; - case 1: + case 1: // max EMCAL + if (!fUsePhiCut) fUsePhiCut=1; + fMaxPhiCut=3.15; + break; + case 2: // max EMCAL with TRD 2011 if (!fUsePhiCut) fUsePhiCut=1; - fMaxPhiCut=3.125; + fMaxPhiCut=2.45; break; + case 3: // max EMCAL with TRD 2012 + if (!fUsePhiCut) fUsePhiCut=1; + fMaxPhiCut=2.10; + break; + + default: AliError(Form("Max Phi Cut not defined %d",maxPhi)); return kFALSE; @@ -1408,43 +1431,43 @@ Bool_t AliCaloPhotonCuts::SetMinEnergyCut(Int_t minEnergy) switch(minEnergy){ case 0: if (!fUseMinEnergy) fUseMinEnergy=0; - fMinEnergy=0; + fMinEnergy=0.1; break; case 1: if (!fUseMinEnergy) fUseMinEnergy=1; - fMinEnergy=0.05; + fMinEnergy=0.2; break; case 2: if (!fUseMinEnergy) fUseMinEnergy=1; - fMinEnergy=0.1; + fMinEnergy=0.3; break; case 3: if (!fUseMinEnergy) fUseMinEnergy=1; - fMinEnergy=0.15; + fMinEnergy=0.4; break; case 4: if (!fUseMinEnergy) fUseMinEnergy=1; - fMinEnergy=0.2; + fMinEnergy=0.5; break; case 5: if (!fUseMinEnergy) fUseMinEnergy=1; - fMinEnergy=0.3; + fMinEnergy=0.6; break; case 6: if (!fUseMinEnergy) fUseMinEnergy=1; - fMinEnergy=0.4; + fMinEnergy=4.5; break; case 7: if (!fUseMinEnergy) fUseMinEnergy=1; - fMinEnergy=0.5; + fMinEnergy=5.0; break; case 8: if (!fUseMinEnergy) fUseMinEnergy=1; - fMinEnergy=0.75; + fMinEnergy=5.5; break; case 9: if (!fUseMinEnergy) fUseMinEnergy=1; - fMinEnergy=1.; + fMinEnergy=6.0; break; default: AliError(Form("Minimum Energy Cut not defined %d",minEnergy)); diff --git a/PWGGA/GammaConv/macros/AddTask_GammaCalo_PbPb.C b/PWGGA/GammaConv/macros/AddTask_GammaCalo_PbPb.C index 40eb2476904..8d61b3983ec 100644 --- a/PWGGA/GammaConv/macros/AddTask_GammaCalo_PbPb.C +++ b/PWGGA/GammaConv/macros/AddTask_GammaCalo_PbPb.C @@ -121,41 +121,41 @@ void AddTask_GammaCalo_PbPb( Int_t trainConfig = 1, //change different set // meson type (Dalitz or not), BG scheme, pool depth, rotation degrees, rapidity cut, radius cut, alpha, chi2, shared electrons, reject to close v0, MC smearing, dca, dca, dca if (trainConfig == 1){ // EMCAL clusters - eventCutArray[ 0] = "6010001"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% - eventCutArray[ 1] = "6120001"; clusterCutArray[1] = "10000040022030000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% - eventCutArray[ 2] = "5010001"; clusterCutArray[2] = "10000040022030000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% - eventCutArray[ 3] = "5240001"; clusterCutArray[3] = "10000040022030000"; mesonCutArray[ 3] = "01525065000000"; // 20-40% - eventCutArray[ 4] = "5250001"; clusterCutArray[4] = "10000040022030000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% + eventCutArray[ 0] = "6010001"; clusterCutArray[0] = "10000040032030000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% + eventCutArray[ 1] = "6120001"; clusterCutArray[1] = "10000040032030000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% + eventCutArray[ 2] = "5010001"; clusterCutArray[2] = "10000040032030000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% + eventCutArray[ 3] = "5240001"; clusterCutArray[3] = "10000040032030000"; mesonCutArray[ 3] = "01525065000000"; // 20-40% + eventCutArray[ 4] = "5250001"; clusterCutArray[4] = "10000040032030000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% } else if (trainConfig == 2){ // EMCAL clusters - eventCutArray[ 0] = "6010001"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% - eventCutArray[ 1] = "6120001"; clusterCutArray[1] = "10000040022030000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% - eventCutArray[ 2] = "5010001"; clusterCutArray[2] = "10000040022030000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% - eventCutArray[ 3] = "5120001"; clusterCutArray[3] = "10000040022030000"; mesonCutArray[ 3] = "01525065000000"; // 10-20% - eventCutArray[ 4] = "5240001"; clusterCutArray[4] = "10000040022030000"; mesonCutArray[ 4] = "01525065000000"; // 20-40% + eventCutArray[ 0] = "6010001"; clusterCutArray[0] = "10000040032030000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% + eventCutArray[ 1] = "6120001"; clusterCutArray[1] = "10000040032030000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% + eventCutArray[ 2] = "5010001"; clusterCutArray[2] = "10000040032030000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% + eventCutArray[ 3] = "5120001"; clusterCutArray[3] = "10000040032030000"; mesonCutArray[ 3] = "01525065000000"; // 10-20% + eventCutArray[ 4] = "5240001"; clusterCutArray[4] = "10000040032030000"; mesonCutArray[ 4] = "01525065000000"; // 20-40% } else if (trainConfig == 3){ // EMCAL clusters - eventCutArray[ 0] = "5460001"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[ 0] = "01525065000000"; // 40-60% - eventCutArray[ 1] = "5680001"; clusterCutArray[1] = "10000040022030000"; mesonCutArray[ 1] = "01525065000000"; // 60-80% - eventCutArray[ 2] = "5260001"; clusterCutArray[2] = "10000040022030000"; mesonCutArray[ 2] = "01525065000000"; // 20-60% - eventCutArray[ 3] = "5480001"; clusterCutArray[3] = "10000040022030000"; mesonCutArray[ 3] = "01525065000000"; // 40-80% - eventCutArray[ 4] = "5250001"; clusterCutArray[4] = "10000040022030000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% + eventCutArray[ 0] = "5460001"; clusterCutArray[0] = "10000040032030000"; mesonCutArray[ 0] = "01525065000000"; // 40-60% + eventCutArray[ 1] = "5680001"; clusterCutArray[1] = "10000040032030000"; mesonCutArray[ 1] = "01525065000000"; // 60-80% + eventCutArray[ 2] = "5260001"; clusterCutArray[2] = "10000040032030000"; mesonCutArray[ 2] = "01525065000000"; // 20-60% + eventCutArray[ 3] = "5480001"; clusterCutArray[3] = "10000040032030000"; mesonCutArray[ 3] = "01525065000000"; // 40-80% + eventCutArray[ 4] = "5250001"; clusterCutArray[4] = "10000040032030000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% } else if (trainConfig == 31){ // PHOS clusters - eventCutArray[ 0] = "6010001"; clusterCutArray[0] = "20000040053200000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% - eventCutArray[ 1] = "6120001"; clusterCutArray[1] = "20000040053200000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% - eventCutArray[ 2] = "5010001"; clusterCutArray[2] = "20000040053200000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% - eventCutArray[ 3] = "5240001"; clusterCutArray[3] = "20000040053200000"; mesonCutArray[ 3] = "01525065000000"; // 20-40% - eventCutArray[ 4] = "5250001"; clusterCutArray[4] = "20000040053200000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% + eventCutArray[ 0] = "6010001"; clusterCutArray[0] = "20000040033200000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% + eventCutArray[ 1] = "6120001"; clusterCutArray[1] = "20000040033200000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% + eventCutArray[ 2] = "5010001"; clusterCutArray[2] = "20000040033200000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% + eventCutArray[ 3] = "5240001"; clusterCutArray[3] = "20000040033200000"; mesonCutArray[ 3] = "01525065000000"; // 20-40% + eventCutArray[ 4] = "5250001"; clusterCutArray[4] = "20000040033200000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% } else if (trainConfig == 32){ // PHOS clusters - eventCutArray[ 0] = "6010001"; clusterCutArray[0] = "20000040053200000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% - eventCutArray[ 1] = "6120001"; clusterCutArray[1] = "20000040053200000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% - eventCutArray[ 2] = "5010001"; clusterCutArray[2] = "20000040053200000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% - eventCutArray[ 3] = "5120001"; clusterCutArray[3] = "20000040053200000"; mesonCutArray[ 3] = "01525065000000"; // 10-20% - eventCutArray[ 4] = "5240001"; clusterCutArray[4] = "20000040053200000"; mesonCutArray[ 4] = "01525065000000"; // 20-40% + eventCutArray[ 0] = "6010001"; clusterCutArray[0] = "20000040033200000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% + eventCutArray[ 1] = "6120001"; clusterCutArray[1] = "20000040033200000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% + eventCutArray[ 2] = "5010001"; clusterCutArray[2] = "20000040033200000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% + eventCutArray[ 3] = "5120001"; clusterCutArray[3] = "20000040033200000"; mesonCutArray[ 3] = "01525065000000"; // 10-20% + eventCutArray[ 4] = "5240001"; clusterCutArray[4] = "20000040033200000"; mesonCutArray[ 4] = "01525065000000"; // 20-40% } else if (trainConfig == 33){ // PHOS clusters - eventCutArray[ 0] = "5460001"; clusterCutArray[0] = "20000040053200000"; mesonCutArray[ 0] = "01525065000000"; // 40-60% - eventCutArray[ 1] = "5680001"; clusterCutArray[1] = "20000040053200000"; mesonCutArray[ 1] = "01525065000000"; // 60-80% - eventCutArray[ 2] = "5260001"; clusterCutArray[2] = "20000040053200000"; mesonCutArray[ 2] = "01525065000000"; // 20-60% - eventCutArray[ 3] = "5480001"; clusterCutArray[3] = "20000040053200000"; mesonCutArray[ 3] = "01525065000000"; // 40-80% - eventCutArray[ 4] = "5250001"; clusterCutArray[4] = "20000040053200000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% + eventCutArray[ 0] = "5460001"; clusterCutArray[0] = "20000040033200000"; mesonCutArray[ 0] = "01525065000000"; // 40-60% + eventCutArray[ 1] = "5680001"; clusterCutArray[1] = "20000040033200000"; mesonCutArray[ 1] = "01525065000000"; // 60-80% + eventCutArray[ 2] = "5260001"; clusterCutArray[2] = "20000040033200000"; mesonCutArray[ 2] = "01525065000000"; // 20-60% + eventCutArray[ 3] = "5480001"; clusterCutArray[3] = "20000040033200000"; mesonCutArray[ 3] = "01525065000000"; // 40-80% + eventCutArray[ 4] = "5250001"; clusterCutArray[4] = "20000040033200000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% } else { Error(Form("GammaConvCalo_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration"); return; diff --git a/PWGGA/GammaConv/macros/AddTask_GammaCalo_pPb.C b/PWGGA/GammaConv/macros/AddTask_GammaCalo_pPb.C index 4af9b4c6464..d12497a3ab1 100644 --- a/PWGGA/GammaConv/macros/AddTask_GammaCalo_pPb.C +++ b/PWGGA/GammaConv/macros/AddTask_GammaCalo_pPb.C @@ -111,8 +111,9 @@ void AddTask_GammaCalo_pPb( 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 = 1;} + if (trainConfig == 2 || trainConfig == 4 ){ numberOfCuts = 6;} + if (trainConfig == 5 || trainConfig == 6 ){ numberOfCuts = 5;} + if (trainConfig == 7 || trainConfig == 8 || trainConfig == 32 || trainConfig == 33){ numberOfCuts = 1;} TString *eventCutArray = new TString[numberOfCuts]; TString *clusterCutArray = new TString[numberOfCuts]; @@ -122,36 +123,51 @@ void AddTask_GammaCalo_pPb( // 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){ + //************************************************ EMCAL clusters ************************************************* + if (trainConfig == 1){ // min energy = 0.3 GeV/c eventCutArray[ 0] = "8000001"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kINT7 // EMCAL clusters eventCutArray[ 1] = "8005201"; clusterCutArray[1] = "10000040022030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMC7 // EMCAL clusters - } else if (trainConfig == 2) { - eventCutArray[ 0] = "8000001"; clusterCutArray[0] = "20000040053200000"; mesonCutArray[0] = "01525065000000"; //standart cut, kINT7 // PHOS clusters - eventCutArray[ 1] = "8006201"; clusterCutArray[1] = "20000040053200000"; mesonCutArray[1] = "01525065000000"; //standard cut, kPHI7 // PHOS clusters - } else if (trainConfig == 3){ + } else if (trainConfig == 2){ // min energy = 0.3 GeV/c eventCutArray[ 0] = "8008101"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kEMCEGA based on INT7 // EMCAL clusters eventCutArray[ 1] = "8008301"; clusterCutArray[1] = "10000040022030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMCEG1 based on INT7 // EMCAL clusters eventCutArray[ 2] = "8008501"; clusterCutArray[2] = "10000040022030000"; mesonCutArray[2] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 // EMCAL clusters eventCutArray[ 3] = "8009101"; clusterCutArray[3] = "10000040022030000"; mesonCutArray[3] = "01525065000000"; //standard cut, kEMCEJE based on INT7 // EMCAL clusters eventCutArray[ 4] = "8009301"; clusterCutArray[4] = "10000040022030000"; mesonCutArray[4] = "01525065000000"; //standard cut, kEMCEJ1 based on INT7 // EMCAL clusters eventCutArray[ 5] = "8009501"; clusterCutArray[5] = "10000040022030000"; mesonCutArray[5] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 // EMCAL clusters - } else if (trainConfig == 4){ - eventCutArray[ 0] = "8000001"; clusterCutArray[0] = "10000040052030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kINT7 // EMCAL clusters - eventCutArray[ 1] = "8005201"; clusterCutArray[1] = "10000040052030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMC7 // EMCAL clusters - } else if (trainConfig == 5){ - eventCutArray[ 0] = "8008101"; clusterCutArray[0] = "10000040052030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kEMCEGA based on INT7 // EMCAL clusters - eventCutArray[ 1] = "8008301"; clusterCutArray[1] = "10000040052030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMCEG1 based on INT7 // EMCAL clusters - eventCutArray[ 2] = "8008501"; clusterCutArray[2] = "10000040052030000"; mesonCutArray[2] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 // EMCAL clusters - eventCutArray[ 3] = "8009101"; clusterCutArray[3] = "10000040052030000"; mesonCutArray[3] = "01525065000000"; //standard cut, kEMCEJE based on INT7 // EMCAL clusters - eventCutArray[ 4] = "8009301"; clusterCutArray[4] = "10000040052030000"; mesonCutArray[4] = "01525065000000"; //standard cut, kEMCEJ1 based on INT7 // EMCAL clusters - eventCutArray[ 5] = "8009501"; clusterCutArray[5] = "10000040052030000"; mesonCutArray[5] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 // EMCAL clusters - } else if (trainConfig == 6){ // Validation EMCAL - eventCutArray[ 0] = "8000001"; clusterCutArray[0] = "10000040062000000"; mesonCutArray[0] = "01630031009000"; - } else if (trainConfig == 7){ // Validation PHOS - eventCutArray[ 0] = "8000001"; clusterCutArray[0] = "20000040053200000"; mesonCutArray[0] = "01630031009000"; + } else if (trainConfig == 3){ // min energy = 0.4 GeV/c + eventCutArray[ 0] = "8000001"; clusterCutArray[0] = "10000040032030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kINT7 // EMCAL clusters + eventCutArray[ 1] = "8005201"; clusterCutArray[1] = "10000040032030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMC7 // EMCAL clusters + } else if (trainConfig == 4){ // min energy = 0.4 GeV/c + eventCutArray[ 0] = "8008101"; clusterCutArray[0] = "10000040032030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kEMCEGA based on INT7 // EMCAL clusters + eventCutArray[ 1] = "8008301"; clusterCutArray[1] = "10000040032030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMCEG1 based on INT7 // EMCAL clusters + eventCutArray[ 2] = "8008501"; clusterCutArray[2] = "10000040032030000"; mesonCutArray[2] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 // EMCAL clusters + eventCutArray[ 3] = "8009101"; clusterCutArray[3] = "10000040032030000"; mesonCutArray[3] = "01525065000000"; //standard cut, kEMCEJE based on INT7 // EMCAL clusters + eventCutArray[ 4] = "8009301"; clusterCutArray[4] = "10000040032030000"; mesonCutArray[4] = "01525065000000"; //standard cut, kEMCEJ1 based on INT7 // EMCAL clusters + eventCutArray[ 5] = "8009501"; clusterCutArray[5] = "10000040032030000"; mesonCutArray[5] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 // EMCAL clusters + } else if (trainConfig == 5){ //EMCAL minEnergy variation + eventCutArray[ 0] = "8000001"; clusterCutArray[0] = "10000040012030000"; mesonCutArray[0] = "01525065000000"; //0.2 GeV/c + eventCutArray[ 1] = "8000001"; clusterCutArray[1] = "10000040022030000"; mesonCutArray[1] = "01525065000000"; //0.3 GeV/c + eventCutArray[ 2] = "8000001"; clusterCutArray[2] = "10000040032030000"; mesonCutArray[2] = "01525065000000"; //0.4 GeV/c default + eventCutArray[ 3] = "8000001"; clusterCutArray[3] = "10000040042030000"; mesonCutArray[3] = "01525065000000"; //0.5 GeV/c + eventCutArray[ 4] = "8000001"; clusterCutArray[4] = "10000040052030000"; mesonCutArray[4] = "01525065000000"; //0.6 GeV/c + } else if (trainConfig == 6){ //EMCAL minNCells variation + eventCutArray[ 0] = "8000001"; clusterCutArray[0] = "10000040031030000"; mesonCutArray[0] = "01525065000000"; //n cells >= 1 + eventCutArray[ 1] = "8000001"; clusterCutArray[1] = "10000040033030000"; mesonCutArray[1] = "01525065000000"; //n cells >= 3 + eventCutArray[ 2] = "8000001"; clusterCutArray[2] = "10000040032000000"; mesonCutArray[2] = "01525065000000"; //no M02 cut + eventCutArray[ 3] = "8000001"; clusterCutArray[3] = "10021040032030000"; mesonCutArray[3] = "01525065000000"; //only modules with TRD infront + eventCutArray[ 4] = "8000001"; clusterCutArray[4] = "10013040032030000"; mesonCutArray[4] = "01525065000000"; //no modules with TRD infront + } else if (trainConfig == 7){ // Validation EMCAL + eventCutArray[ 0] = "8000001"; clusterCutArray[0] = "10000040032000000"; mesonCutArray[0] = "01630031009000"; } else if (trainConfig == 8){ // Validation EMCAL, only added signals - eventCutArray[ 0] = "8000002"; clusterCutArray[0] = "10000040062000000"; mesonCutArray[0] = "01630031009000"; - } else if (trainConfig == 9){ // Validation PHOS, only added signals + eventCutArray[ 0] = "8000002"; clusterCutArray[0] = "10000040032000000"; mesonCutArray[0] = "01630031009000"; + + //************************************************ PHOS clusters ************************************************* + } else if (trainConfig == 31) { // min energy = 0.3 GeV/c + eventCutArray[ 0] = "8000001"; clusterCutArray[0] = "20000040033200000"; mesonCutArray[0] = "01525065000000"; //standart cut, kINT7 // PHOS clusters + eventCutArray[ 1] = "8006201"; clusterCutArray[1] = "20000040033200000"; mesonCutArray[1] = "01525065000000"; //standard cut, kPHI7 // PHOS clusters + } else if (trainConfig == 32){ // Validation PHOS + eventCutArray[ 0] = "8000001"; clusterCutArray[0] = "20000040053200000"; mesonCutArray[0] = "01630031009000"; + } else if (trainConfig == 33){ // Validation PHOS, only added signals eventCutArray[ 0] = "8000002"; clusterCutArray[0] = "20000040053200000"; mesonCutArray[0] = "01630031009000"; diff --git a/PWGGA/GammaConv/macros/AddTask_GammaCalo_pp.C b/PWGGA/GammaConv/macros/AddTask_GammaCalo_pp.C index ad66807678c..54b07c3a7b1 100644 --- a/PWGGA/GammaConv/macros/AddTask_GammaCalo_pp.C +++ b/PWGGA/GammaConv/macros/AddTask_GammaCalo_pp.C @@ -107,9 +107,11 @@ void AddTask_GammaCalo_pp( Int_t trainConfig = 1, //change different set of cu task->SetIsHeavyIon(isHeavyIon); task->SetIsMC(isMC); // Cut Numbers to use in Analysis - Int_t numberOfCuts = 4; - if (trainConfig == 9) numberOfCuts = 6; - if (trainConfig == 10) numberOfCuts = 1; + Int_t numberOfCuts = 2; + if (trainConfig == 2 || trainConfig == 3) numberOfCuts = 5; + if (trainConfig == 5) numberOfCuts = 5; + if (trainConfig == 31) numberOfCuts = 4; + if (trainConfig == 32) numberOfCuts = 1; TString *eventCutArray = new TString[numberOfCuts]; TString *clusterCutArray = new TString[numberOfCuts]; @@ -118,62 +120,41 @@ void AddTask_GammaCalo_pp( Int_t trainConfig = 1, //change different set of cu // cluster cuts // 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){ // EMCAL clusters 2.76 TeV LHC11a, with SDD (0,1), kEMC1 (2,3) - eventCutArray[ 0] = "0000311"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[0] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 1] = "0000311"; clusterCutArray[1] = "10000040052030000"; mesonCutArray[1] = "01631031009000"; // 300 MeV cluster min energy - eventCutArray[ 2] = "0005111"; clusterCutArray[2] = "10000040022030000"; mesonCutArray[2] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 3] = "0005111"; clusterCutArray[3] = "10000040052030000"; mesonCutArray[3] = "01631031009000"; // 300 MeV cluster min energy - } else if (trainConfig == 2){ // EMCAL clusters, EMCEGA triggers - eventCutArray[ 0] = "0008311"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[0] = "01631031009000"; // EMCEG1, 100 MeV cluster min energy - eventCutArray[ 1] = "0008311"; clusterCutArray[1] = "10000040052030000"; mesonCutArray[1] = "01631031009000"; // EMCEG1, 300 MeV cluster min energy - eventCutArray[ 2] = "0008511"; clusterCutArray[2] = "10000040022030000"; mesonCutArray[2] = "01631031009000"; // EMCEG2, 100 MeV cluster min energy - eventCutArray[ 3] = "0008511"; clusterCutArray[3] = "10000040052030000"; mesonCutArray[3] = "01631031009000"; // EMCEG2, 300 MeV cluster min energy - } else if (trainConfig == 3){ // EMCAL clusters, EMCEJE triggers - eventCutArray[ 0] = "0009311"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[0] = "01631031009000"; // EMCEJ1, 100 MeV cluster min energy - eventCutArray[ 1] = "0009311"; clusterCutArray[1] = "10000040052030000"; mesonCutArray[1] = "01631031009000"; // EMCEJ1, 300 MeV cluster min energy - eventCutArray[ 2] = "0009511"; clusterCutArray[2] = "10000040022030000"; mesonCutArray[2] = "01631031009000"; // EMCEJ2, 100 MeV cluster min energy - eventCutArray[ 3] = "0009511"; clusterCutArray[3] = "10000040052030000"; mesonCutArray[3] = "01631031009000"; // EMCEJ2, 300 MeV cluster min energy - } else if (trainConfig == 4){ // EMCAL clusters 2.76 TeV LHC13g, kINT7 (0,1), kEMC7 (2,3), track - eventCutArray[ 0] = "0000011"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[0] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 1] = "0000011"; clusterCutArray[1] = "10000040052030000"; mesonCutArray[1] = "01631031009000"; // 300 MeV cluster min energy - eventCutArray[ 2] = "0005211"; clusterCutArray[2] = "10000040022030000"; mesonCutArray[2] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 3] = "0005211"; clusterCutArray[3] = "10000040052030000"; mesonCutArray[3] = "01631031009000"; // 300 MeV cluster min energy - } else if (trainConfig == 5){ // EMCAL clusters 2.76 TeV LHC11a, with SDD (0,1), kEMC1 (2,3) - eventCutArray[ 0] = "0000311"; clusterCutArray[0] = "10000040062030000"; mesonCutArray[0] = "01631031009000"; // 400 MeV cluster min energy - eventCutArray[ 1] = "0000311"; clusterCutArray[1] = "10000040062031000"; mesonCutArray[1] = "01631031009000"; // 400 MeV cluster min energy, min M20 > 0.02 - eventCutArray[ 2] = "0005111"; clusterCutArray[2] = "10000040062030000"; mesonCutArray[2] = "01631031009000"; // 500 MeV cluster min energy - eventCutArray[ 3] = "0005111"; clusterCutArray[3] = "10000040062031000"; mesonCutArray[3] = "01631031009000"; // 500 MeV cluster min energy, min M20 > 0.02 - } else if (trainConfig == 6){ // EMCAL clusters, EMCEGA triggers - eventCutArray[ 0] = "0008311"; clusterCutArray[0] = "10000040062030000"; mesonCutArray[0] = "01631031009000"; // EMCEG1, 400 MeV cluster min energy - eventCutArray[ 1] = "0008311"; clusterCutArray[1] = "10000040062031000"; mesonCutArray[1] = "01631031009000"; // EMCEG1, 400 MeV cluster min energy, min M20 > 0.02 - eventCutArray[ 2] = "0008511"; clusterCutArray[2] = "10000040062030000"; mesonCutArray[2] = "01631031009000"; // EMCEG2, 400 MeV cluster min energy - eventCutArray[ 3] = "0008511"; clusterCutArray[3] = "10000040062031000"; mesonCutArray[3] = "01631031009000"; // EMCEG2, 400 MeV cluster min energy, min M20 > 0.02 - } else if (trainConfig == 7){ // EMCAL clusters, EMCEJE triggers - eventCutArray[ 0] = "0009311"; clusterCutArray[0] = "10000040062030000"; mesonCutArray[0] = "01631031009000"; // EMCEJ1, 400 MeV cluster min energy - eventCutArray[ 1] = "0009311"; clusterCutArray[1] = "10000040062031000"; mesonCutArray[1] = "01631031009000"; // EMCEJ1, 400 MeV cluster min energy, min M20 > 0.02 - eventCutArray[ 2] = "0009511"; clusterCutArray[2] = "10000040062030000"; mesonCutArray[2] = "01631031009000"; // EMCEJ2, 400 MeV cluster min energy - eventCutArray[ 3] = "0009511"; clusterCutArray[3] = "10000040062031000"; mesonCutArray[3] = "01631031009000"; // EMCEJ2, 400 MeV cluster min energy, min M20 > 0.02 - } else if (trainConfig == 8){ // EMCAL clusters 2.76 TeV LHC13g, kINT7 (0,1), kEMC7 (2,3), track - eventCutArray[ 0] = "0000011"; clusterCutArray[0] = "10000040062030000"; mesonCutArray[0] = "01631031009000"; // 400 MeV cluster min energy - eventCutArray[ 1] = "0000011"; clusterCutArray[1] = "10000040062031000"; mesonCutArray[1] = "01631031009000"; // 400 MeV cluster min energy, min M20 > 0.02 - eventCutArray[ 2] = "0005211"; clusterCutArray[2] = "10000040062030000"; mesonCutArray[2] = "01631031009000"; // 400 MeV cluster min energy - eventCutArray[ 3] = "0005211"; clusterCutArray[3] = "10000040062031000"; mesonCutArray[3] = "01631031009000"; // 400 MeV cluster min energy, min M20 > 0.02 - } else if (trainConfig == 9){ // Trigger test - eventCutArray[ 0] = "0000311"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[0] = "01631031009000"; - eventCutArray[ 1] = "0005011"; clusterCutArray[1] = "10000040022030000"; mesonCutArray[1] = "01631031009000"; - eventCutArray[ 2] = "0008311"; clusterCutArray[2] = "10000040022030000"; mesonCutArray[2] = "01631031009000"; - eventCutArray[ 3] = "0009311"; clusterCutArray[3] = "10000040022030000"; mesonCutArray[3] = "01631031009000"; - eventCutArray[ 4] = "0008511"; clusterCutArray[4] = "10000040022030000"; mesonCutArray[4] = "01631031009000"; - eventCutArray[ 5] = "0009511"; clusterCutArray[5] = "10000040022030000"; mesonCutArray[5] = "01631031009000"; - } else if (trainConfig == 10){ // Validation EMCAL - eventCutArray[ 0] = "0000311"; clusterCutArray[0] = "10000040062000000"; mesonCutArray[0] = "01630031009000"; + + // ************************************* EMCAL cuts **************************************************** + // LHC11a + if (trainConfig == 1){ // EMCAL clusters 2.76 TeV LHC11a, with SDD (0), kEMC1 (1) + eventCutArray[ 0] = "0000311"; clusterCutArray[0] = "10000040032030000"; mesonCutArray[0] = "01631031009000"; // 400 MeV cluster min energy + eventCutArray[ 1] = "0005111"; clusterCutArray[1] = "10000040032030000"; mesonCutArray[1] = "01631031009000"; // 400 MeV cluster min energy + } else if (trainConfig == 2){ //EMCAL minEnergy variation + eventCutArray[ 0] = "0000311"; clusterCutArray[0] = "10000040012030000"; mesonCutArray[0] = "01525065000000"; //0.2 GeV/c + eventCutArray[ 1] = "0000311"; clusterCutArray[1] = "10000040022030000"; mesonCutArray[1] = "01525065000000"; //0.3 GeV/c + eventCutArray[ 2] = "0000311"; clusterCutArray[2] = "10000040032030000"; mesonCutArray[2] = "01525065000000"; //0.4 GeV/c default + eventCutArray[ 3] = "0000311"; clusterCutArray[3] = "10000040042030000"; mesonCutArray[3] = "01525065000000"; //0.5 GeV/c + eventCutArray[ 4] = "0000311"; clusterCutArray[4] = "10000040052030000"; mesonCutArray[4] = "01525065000000"; //0.6 GeV/c + } else if (trainConfig == 3){ //EMCAL minNCells variation + eventCutArray[ 0] = "0000311"; clusterCutArray[0] = "10000040031030000"; mesonCutArray[0] = "01525065000000"; //n cells >= 1 + eventCutArray[ 1] = "0000311"; clusterCutArray[1] = "10000040033030000"; mesonCutArray[1] = "01525065000000"; //n cells >= 3 + eventCutArray[ 2] = "0000311"; clusterCutArray[2] = "10000040032000000"; mesonCutArray[2] = "01525065000000"; //no M02 cut + eventCutArray[ 3] = "0000311"; clusterCutArray[3] = "10031040032030000"; mesonCutArray[3] = "01525065000000"; //only modules with TRD infront + eventCutArray[ 4] = "0000311"; clusterCutArray[4] = "10012040032030000"; mesonCutArray[4] = "01525065000000"; //no modules with TRD infront + // LHC13g + } else if (trainConfig == 5){ // EMCAL clusters, EMCEGA triggers + eventCutArray[ 0] = "0008311"; clusterCutArray[0] = "10000040032030000"; mesonCutArray[0] = "01631031009000"; // EMCEG1, + eventCutArray[ 1] = "0008511"; clusterCutArray[1] = "10000040032030000"; mesonCutArray[1] = "01631031009000"; // EMCEG2, + eventCutArray[ 2] = "0009311"; clusterCutArray[2] = "10000040032030000"; mesonCutArray[2] = "01631031009000"; // EMCEJ1, + eventCutArray[ 3] = "0009511"; clusterCutArray[3] = "10000040032030000"; mesonCutArray[3] = "01631031009000"; // EMCEJ2, + eventCutArray[ 4] = "0000011"; clusterCutArray[4] = "10000040032030000"; mesonCutArray[4] = "01631031009000"; // INT7 + eventCutArray[ 5] = "0005211"; clusterCutArray[5] = "10000040032030000"; mesonCutArray[5] = "01631031009000"; // EMC7 + + // ************************************* PHOS cuts **************************************************** } else if (trainConfig == 31) { //PHOS clusters - eventCutArray[ 0] = "0000311"; clusterCutArray[0] = "20000040053200000"; mesonCutArray[0] = "01631031009000"; //pp LHC11a with SDD, PHOS - eventCutArray[ 1] = "0000011"; clusterCutArray[1] = "20000040053200000"; mesonCutArray[1] = "01631031009000"; //pp LHC13g default MB - eventCutArray[ 2] = "0006111"; clusterCutArray[2] = "20000040053200000"; mesonCutArray[2] = "01631031009000"; //pp LHC11a PHI1 - eventCutArray[ 3] = "0006211"; clusterCutArray[3] = "20000040053200000"; mesonCutArray[3] = "01631031009000"; //pp LHC11a PHI7 + eventCutArray[ 0] = "0000311"; clusterCutArray[0] = "20000040033200000"; mesonCutArray[0] = "01631031009000"; //pp LHC11a with SDD, PHOS + eventCutArray[ 1] = "0000011"; clusterCutArray[1] = "20000040033200000"; mesonCutArray[1] = "01631031009000"; //pp LHC13g default MB + eventCutArray[ 2] = "0006111"; clusterCutArray[2] = "20000040033200000"; mesonCutArray[2] = "01631031009000"; //pp LHC11a PHI1 + eventCutArray[ 3] = "0006211"; clusterCutArray[3] = "20000040033200000"; mesonCutArray[3] = "01631031009000"; //pp LHC11a PHI7 } else if (trainConfig == 32){ // Validation PHOS - eventCutArray[ 0] = "0000311"; clusterCutArray[0] = "20000040053200000"; mesonCutArray[0] = "01630031009000"; + eventCutArray[ 0] = "0000311"; clusterCutArray[0] = "20000040033200000"; mesonCutArray[0] = "01630031009000"; } else { Error(Form("GammaCalo_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration"); return; diff --git a/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_PbPb.C b/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_PbPb.C index f654ecad748..e8cd1b4a31e 100644 --- a/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_PbPb.C +++ b/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_PbPb.C @@ -118,82 +118,79 @@ void AddTask_GammaConvCalo_PbPb( Int_t trainConfig = 1, //change different s TString *photonCutArray = new TString[numberOfCuts]; TString *clusterCutArray = new TString[numberOfCuts]; TString *mesonCutArray = new TString[numberOfCuts]; - - // meson cuts - // meson type (Dalitz or not), BG scheme, pool depth, rotation degrees, rapidity cut, radius cut, alpha, chi2, shared electrons, reject to close v0, MC smearing, dca, dca, dca if (trainConfig == 1){ // EMCAL clusters - eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% - eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000040022030000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% - eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000040022030000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% - eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000040022030000"; mesonCutArray[ 3] = "01525065000000"; // 20-40% - eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000040022030000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% + eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000040032030000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% + eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000040032030000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% + eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000040032030000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% + eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000040032030000"; mesonCutArray[ 3] = "01525065000000"; // 20-40% + eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000040032030000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% } else if (trainConfig == 2){ // EMCAL clusters - eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% - eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000040022030000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% - eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000040022030000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% - eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000040022030000"; mesonCutArray[ 3] = "01525065000000"; // 10-20% - eventCutArray[ 4] = "5240001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000040022030000"; mesonCutArray[ 4] = "01525065000000"; // 20-40% + eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000040032030000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% + eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000040032030000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% + eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000040032030000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% + eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000040032030000"; mesonCutArray[ 3] = "01525065000000"; // 10-20% + eventCutArray[ 4] = "5240001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000040032030000"; mesonCutArray[ 4] = "01525065000000"; // 20-40% } else if (trainConfig == 3){ // EMCAL clusters - eventCutArray[ 0] = "5460001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[ 0] = "01525065000000"; // 40-60% - eventCutArray[ 1] = "5680001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000040022030000"; mesonCutArray[ 1] = "01525065000000"; // 60-80% - eventCutArray[ 2] = "5260001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000040022030000"; mesonCutArray[ 2] = "01525065000000"; // 20-60% - eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000040022030000"; mesonCutArray[ 3] = "01525065000000"; // 40-80% - eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000040022030000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% + eventCutArray[ 0] = "5460001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000040032030000"; mesonCutArray[ 0] = "01525065000000"; // 40-60% + eventCutArray[ 1] = "5680001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000040032030000"; mesonCutArray[ 1] = "01525065000000"; // 60-80% + eventCutArray[ 2] = "5260001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000040032030000"; mesonCutArray[ 2] = "01525065000000"; // 20-60% + eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000040032030000"; mesonCutArray[ 3] = "01525065000000"; // 40-80% + eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000040032030000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% } else if (trainConfig == 4){ // EMCAL clusters - eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000043022030000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% - eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000043022030000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% - eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000043022030000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% - eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000043022030000"; mesonCutArray[ 3] = "01525065000000"; // 20-40% - eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000043022030000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% + eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000043032030000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% + eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000043032030000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% + eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% + eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000043032030000"; mesonCutArray[ 3] = "01525065000000"; // 20-40% + eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000043032030000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% } else if (trainConfig == 5){ // EMCAL clusters - eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000043022030000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% - eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000043022030000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% - eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000043022030000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% - eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000043022030000"; mesonCutArray[ 3] = "01525065000000"; // 10-20% - eventCutArray[ 4] = "5240001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000043022030000"; mesonCutArray[ 4] = "01525065000000"; // 20-40% + eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000043032030000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% + eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000043032030000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% + eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% + eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000043032030000"; mesonCutArray[ 3] = "01525065000000"; // 10-20% + eventCutArray[ 4] = "5240001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000043032030000"; mesonCutArray[ 4] = "01525065000000"; // 20-40% } else if (trainConfig == 6){ // EMCAL clusters - eventCutArray[ 0] = "5460001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000043022030000"; mesonCutArray[ 0] = "01525065000000"; // 40-60% - eventCutArray[ 1] = "5680001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000043022030000"; mesonCutArray[ 1] = "01525065000000"; // 60-80% - eventCutArray[ 2] = "5260001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000043022030000"; mesonCutArray[ 2] = "01525065000000"; // 20-60% - eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000043022030000"; mesonCutArray[ 3] = "01525065000000"; // 40-80% - eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000043022030000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% + eventCutArray[ 0] = "5460001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000043032030000"; mesonCutArray[ 0] = "01525065000000"; // 40-60% + eventCutArray[ 1] = "5680001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000043032030000"; mesonCutArray[ 1] = "01525065000000"; // 60-80% + eventCutArray[ 2] = "5260001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[ 2] = "01525065000000"; // 20-60% + eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000043032030000"; mesonCutArray[ 3] = "01525065000000"; // 40-80% + eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000043032030000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% } else if (trainConfig == 31){ // PHOS clusters - eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000042053200000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% - eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000042053200000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% - eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000042053200000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% - eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000042053200000"; mesonCutArray[ 3] = "01525065000000"; // 20-40% - eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000042053200000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% + eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000042033200000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% + eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% + eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000042033200000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% + eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000042033200000"; mesonCutArray[ 3] = "01525065000000"; // 20-40% + eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000042033200000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% } else if (trainConfig == 32){ // PHOS clusters - eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000042053200000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% - eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000042053200000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% - eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000042053200000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% - eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000042053200000"; mesonCutArray[ 3] = "01525065000000"; // 10-20% - eventCutArray[ 4] = "5240001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000042053200000"; mesonCutArray[ 4] = "01525065000000"; // 20-40% + eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000042033200000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% + eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% + eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000042033200000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% + eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000042033200000"; mesonCutArray[ 3] = "01525065000000"; // 10-20% + eventCutArray[ 4] = "5240001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000042033200000"; mesonCutArray[ 4] = "01525065000000"; // 20-40% } else if (trainConfig == 33){ // PHOS clusters - eventCutArray[ 0] = "5460001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000042053200000"; mesonCutArray[ 0] = "01525065000000"; // 40-60% - eventCutArray[ 1] = "5680001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000042053200000"; mesonCutArray[ 1] = "01525065000000"; // 60-80% - eventCutArray[ 2] = "5260001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000042053200000"; mesonCutArray[ 2] = "01525065000000"; // 20-60% - eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000042053200000"; mesonCutArray[ 3] = "01525065000000"; // 40-80% - eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000042053200000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% + eventCutArray[ 0] = "5460001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000042033200000"; mesonCutArray[ 0] = "01525065000000"; // 40-60% + eventCutArray[ 1] = "5680001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[ 1] = "01525065000000"; // 60-80% + eventCutArray[ 2] = "5260001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000042033200000"; mesonCutArray[ 2] = "01525065000000"; // 20-60% + eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000042033200000"; mesonCutArray[ 3] = "01525065000000"; // 40-80% + eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000042033200000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% } else if (trainConfig == 34){ // PHOS clusters - eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000043053200000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% - eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000043053200000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% - eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000043053200000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% - eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000043053200000"; mesonCutArray[ 3] = "01525065000000"; // 20-40% - eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000043053200000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% + eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000043033200000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% + eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000043033200000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% + eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% + eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000043033200000"; mesonCutArray[ 3] = "01525065000000"; // 20-40% + eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000043033200000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% } else if (trainConfig == 35){ // PHOS clusters - eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000043053200000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% - eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000043053200000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% - eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000043053200000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% - eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000043053200000"; mesonCutArray[ 3] = "01525065000000"; // 10-20% - eventCutArray[ 4] = "5240001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000043053200000"; mesonCutArray[ 4] = "01525065000000"; // 20-40% + eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000043033200000"; mesonCutArray[ 0] = "01525065000000"; // 0-5% + eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000043033200000"; mesonCutArray[ 1] = "01525065000000"; // 5-10% + eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[ 2] = "01525065000000"; // 0-10% + eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000043033200000"; mesonCutArray[ 3] = "01525065000000"; // 10-20% + eventCutArray[ 4] = "5240001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000043033200000"; mesonCutArray[ 4] = "01525065000000"; // 20-40% } else if (trainConfig == 36){ // PHOS clusters - eventCutArray[ 0] = "5460001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000043053200000"; mesonCutArray[ 0] = "01525065000000"; // 40-60% - eventCutArray[ 1] = "5680001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000043053200000"; mesonCutArray[ 1] = "01525065000000"; // 60-80% - eventCutArray[ 2] = "5260001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000043053200000"; mesonCutArray[ 2] = "01525065000000"; // 20-60% - eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000043053200000"; mesonCutArray[ 3] = "01525065000000"; // 40-80% - eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000043053200000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% + eventCutArray[ 0] = "5460001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "20000043033200000"; mesonCutArray[ 0] = "01525065000000"; // 40-60% + eventCutArray[ 1] = "5680001"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "20000043033200000"; mesonCutArray[ 1] = "01525065000000"; // 60-80% + eventCutArray[ 2] = "5260001"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[ 2] = "01525065000000"; // 20-60% + eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "20000043033200000"; mesonCutArray[ 3] = "01525065000000"; // 40-80% + eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "20000043033200000"; mesonCutArray[ 4] = "01525065000000"; // 20-50% } else { Error(Form("GammaConvCalo_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration"); return; diff --git a/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pPb.C b/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pPb.C index a3391cbde62..4a68a2e9583 100644 --- 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; diff --git a/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pp.C b/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pp.C index d7bffd2f366..f47f44c4aae 100644 --- a/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pp.C +++ b/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pp.C @@ -111,8 +111,10 @@ void AddTask_GammaConvCalo_pp( Int_t trainConfig = 1, //change different set o task->SetIsHeavyIon(isHeavyIon); task->SetIsMC(isMC); // Cut Numbers to use in Analysis - Int_t numberOfCuts = 4; - + Int_t numberOfCuts = 2; + if (trainConfig==2 || trainConfig==3 || trainConfig==4 ){ numberOfCuts =5;} + if (trainConfig==5 || trainConfig==31 || trainConfig==32 ){ numberOfCuts =6;} + TString *eventCutArray = new TString[numberOfCuts]; TString *photonCutArray = new TString[numberOfCuts]; TString *clusterCutArray = new TString[numberOfCuts]; @@ -121,82 +123,56 @@ void AddTask_GammaConvCalo_pp( Int_t trainConfig = 1, //change different set o // cluster cuts // 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){ // EMCAL clusters 2.76 TeV LHC11a, with SDD (0,1), kEMC1 (2,3) - eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000041022030000"; mesonCutArray[0] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000041052030000"; mesonCutArray[1] = "01631031009000"; // 300 MeV cluster min energy - eventCutArray[ 2] = "0005111"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000041022030000"; mesonCutArray[2] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 3] = "0005111"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000041052030000"; mesonCutArray[3] = "01631031009000"; // 300 MeV cluster min energy - } else if (trainConfig == 2){ // EMCAL clusters, EMCEGA triggers - eventCutArray[ 0] = "0008311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000041022030000"; mesonCutArray[0] = "01631031009000"; // EMCEG1, 100 MeV cluster min energy - eventCutArray[ 1] = "0008311"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000041052030000"; mesonCutArray[1] = "01631031009000"; // EMCEG1, 300 MeV cluster min energy - eventCutArray[ 2] = "0008511"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000041022030000"; mesonCutArray[2] = "01631031009000"; // EMCEG2, 100 MeV cluster min energy - eventCutArray[ 3] = "0008511"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000041052030000"; mesonCutArray[3] = "01631031009000"; // EMCEG2, 300 MeV cluster min energy - } else if (trainConfig == 3){ // EMCAL clusters, EMCEJE triggers - eventCutArray[ 0] = "0009311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000041022030000"; mesonCutArray[0] = "01631031009000"; // EMCEJ1, 100 MeV cluster min energy - eventCutArray[ 1] = "0009311"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000041052030000"; mesonCutArray[1] = "01631031009000"; // EMCEJ1, 300 MeV cluster min energy - eventCutArray[ 2] = "0009511"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000041022030000"; mesonCutArray[2] = "01631031009000"; // EMCEJ2, 100 MeV cluster min energy - eventCutArray[ 3] = "0009511"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000041052030000"; mesonCutArray[3] = "01631031009000"; // EMCEJ2, 300 MeV cluster min energy - } else if (trainConfig == 4){ // EMCAL clusters 2.76 TeV LHC11a, with SDD (0,1), kEMC1 (2,3), track matching 0.035 - eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000042022030000"; mesonCutArray[0] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000042052030000"; mesonCutArray[1] = "01631031009000"; // 300 MeV cluster min energy - eventCutArray[ 2] = "0005111"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000042022030000"; mesonCutArray[2] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 3] = "0005111"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000042052030000"; mesonCutArray[3] = "01631031009000"; // 300 MeV cluster min energy + + // ************************************* EMCAL cuts **************************************************** + // LHC11a + if (trainConfig == 1){ // EMCAL clusters 2.76 TeV LHC11a, with SDD (0), kEMC1 (1) + eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031009000"; // 400 MeV cluster min energy + eventCutArray[ 1] = "0005111"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031009000"; // 400 MeV cluster min energy + } else if (trainConfig == 2){ //EMCAL minEnergy variation + eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000042012030000"; mesonCutArray[0] = "01525065000000"; //0.2 GeV/c + eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000042022030000"; mesonCutArray[1] = "01525065000000"; //0.3 GeV/c + eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000042032030000"; mesonCutArray[2] = "01525065000000"; //0.4 GeV/c default + eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000042042030000"; mesonCutArray[3] = "01525065000000"; //0.5 GeV/c + eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002093663003800000000"; clusterCutArray[4] = "10000042052030000"; mesonCutArray[4] = "01525065000000"; //0.6 GeV/c + } else if (trainConfig == 3){ //EMCAL minNCells variation + eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000042031030000"; mesonCutArray[0] = "01525065000000"; //n cells >= 1 + eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000042033030000"; mesonCutArray[1] = "01525065000000"; //n cells >= 3 + eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000042032000000"; mesonCutArray[2] = "01525065000000"; //no M02 cut + eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10031042032030000"; mesonCutArray[3] = "01525065000000"; //only modules with TRD infront + eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002093663003800000000"; clusterCutArray[4] = "10012042032030000"; mesonCutArray[4] = "01525065000000"; //no modules with TRD infront + } else if (trainConfig == 4){ // EMCAL track matching variations + eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000041032030000"; mesonCutArray[0] = "01631031009000"; // + eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031009000"; // + eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000043032030000"; mesonCutArray[2] = "01631031009000"; // + eventCutArray[ 3] = "0000311"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000044032030000"; mesonCutArray[3] = "01631031009000"; // + eventCutArray[ 4] = "0000311"; photonCutArray[ 4] = "002093663003800000000"; clusterCutArray[4] = "10000045032030000"; mesonCutArray[4] = "01631031009000"; // + // LHC13g } else if (trainConfig == 5){ // EMCAL clusters, EMCEGA triggers, track matching 0.035 - eventCutArray[ 0] = "0008311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000042022030000"; mesonCutArray[0] = "01631031009000"; // EMCEG1, 100 MeV cluster min energy - eventCutArray[ 1] = "0008311"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000042052030000"; mesonCutArray[1] = "01631031009000"; // EMCEG1, 300 MeV cluster min energy - eventCutArray[ 2] = "0008511"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000042022030000"; mesonCutArray[2] = "01631031009000"; // EMCEG2, 100 MeV cluster min energy - eventCutArray[ 3] = "0008511"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000042052030000"; mesonCutArray[3] = "01631031009000"; // EMCEG2, 300 MeV cluster min energy - } else if (trainConfig == 6){ // EMCAL clusters, EMCEJE triggers, track matching 0.035 - eventCutArray[ 0] = "0009311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000042022030000"; mesonCutArray[0] = "01631031009000"; // EMCEJ1, 100 MeV cluster min energy - eventCutArray[ 1] = "0009311"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000042052030000"; mesonCutArray[1] = "01631031009000"; // EMCEJ1, 300 MeV cluster min energy - eventCutArray[ 2] = "0009511"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000042022030000"; mesonCutArray[2] = "01631031009000"; // EMCEJ2, 100 MeV cluster min energy - eventCutArray[ 3] = "0009511"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000042052030000"; mesonCutArray[3] = "01631031009000"; // EMCEJ2, 300 MeV cluster min energy - } else if (trainConfig == 7){ // EMCAL clusters 2.76 TeV LHC11a, with SDD (0,1), kEMC1 (2,3), track matching 0.04 - eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000043022030000"; mesonCutArray[0] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000043052030000"; mesonCutArray[1] = "01631031009000"; // 300 MeV cluster min energy - eventCutArray[ 2] = "0005111"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000043022030000"; mesonCutArray[2] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 3] = "0005111"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000043052030000"; mesonCutArray[3] = "01631031009000"; // 300 MeV cluster min energy - } else if (trainConfig == 8){ // EMCAL clusters, EMCEGA triggers, track matching 0.04 - eventCutArray[ 0] = "0008311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000043022030000"; mesonCutArray[0] = "01631031009000"; // EMCEG1, 100 MeV cluster min energy - eventCutArray[ 1] = "0008311"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000043052030000"; mesonCutArray[1] = "01631031009000"; // EMCEG1, 300 MeV cluster min energy - eventCutArray[ 2] = "0008511"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000043022030000"; mesonCutArray[2] = "01631031009000"; // EMCEG2, 100 MeV cluster min energy - eventCutArray[ 3] = "0008511"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000043052030000"; mesonCutArray[3] = "01631031009000"; // EMCEG2, 300 MeV cluster min energy - } else if (trainConfig == 9){ // EMCAL clusters, EMCEJE triggers, track matching 0.04 - eventCutArray[ 0] = "0009311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000043022030000"; mesonCutArray[0] = "01631031009000"; // EMCEJ1, 100 MeV cluster min energy - eventCutArray[ 1] = "0009311"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000043052030000"; mesonCutArray[1] = "01631031009000"; // EMCEJ1, 300 MeV cluster min energy - eventCutArray[ 2] = "0009511"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000043022030000"; mesonCutArray[2] = "01631031009000"; // EMCEJ2, 100 MeV cluster min energy - eventCutArray[ 3] = "0009511"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000043052030000"; mesonCutArray[3] = "01631031009000"; // EMCEJ2, 300 MeV cluster min energy - } else if (trainConfig == 10){ // EMCAL clusters 2.76 TeV LHC11g, INT7 (0,1), kEMC7 (2,3) - eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000040022030000"; mesonCutArray[0] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000040052030000"; mesonCutArray[1] = "01631031009000"; // 300 MeV cluster min energy - eventCutArray[ 2] = "0005211"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000040022030000"; mesonCutArray[2] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 3] = "0005211"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000040052030000"; mesonCutArray[3] = "01631031009000"; // 300 MeV cluster min energy - } else if (trainConfig == 11){ // EMCAL clusters 2.76 TeV LHC11g, INT7 (0,1), kEMC7 (2,3), track matching 0.035 - eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000042022030000"; mesonCutArray[0] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000042052030000"; mesonCutArray[1] = "01631031009000"; // 300 MeV cluster min energy - eventCutArray[ 2] = "0005211"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000042022030000"; mesonCutArray[2] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 3] = "0005211"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000042052030000"; mesonCutArray[3] = "01631031009000"; // 300 MeV cluster min energy - } else if (trainConfig == 12){ // EMCAL clusters 2.76 TeV LHC11g, INT7 (0,1), kEMC7 (2,3), track matching 0.04 - eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000043022030000"; mesonCutArray[0] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000043052030000"; mesonCutArray[1] = "01631031009000"; // 300 MeV cluster min energy - eventCutArray[ 2] = "0005211"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000043022030000"; mesonCutArray[2] = "01631031009000"; // 100 MeV cluster min energy - eventCutArray[ 3] = "0005211"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000043052030000"; mesonCutArray[3] = "01631031009000"; // 300 MeV cluster min energy - } else if (trainConfig == 31) { //PHOS clusters, track matching 0.03 - eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "20000041053200000"; mesonCutArray[0] = "01631031009000"; //pp LHC11a with SDD, PHOS - eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "20000041053200000"; mesonCutArray[1] = "01631031009000"; //pp LHC13g default MB - eventCutArray[ 2] = "0006111"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "20000041053200000"; mesonCutArray[2] = "01631031009000"; //pp LHC11a PHI1 - eventCutArray[ 3] = "0006211"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "20000041053200000"; mesonCutArray[3] = "01631031009000"; //pp LHC11a PHI7 - } else if (trainConfig == 32) { //PHOS clusters, track matching 0.035 - eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "20000042053200000"; mesonCutArray[0] = "01631031009000"; //pp LHC11a with SDD, PHOS - eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "20000042053200000"; mesonCutArray[1] = "01631031009000"; //pp LHC13g default MB - eventCutArray[ 2] = "0006111"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "20000042053200000"; mesonCutArray[2] = "01631031009000"; //pp LHC11a PHI1 - eventCutArray[ 3] = "0006211"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "20000042053200000"; mesonCutArray[3] = "01631031009000"; //pp LHC11a PHI7 - } else if (trainConfig == 33) { //PHOS clusters, track matching 0.04 - eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "20000043053200000"; mesonCutArray[0] = "01631031009000"; //pp LHC11a with SDD, PHOS - eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "20000043053200000"; mesonCutArray[1] = "01631031009000"; //pp LHC13g default MB - eventCutArray[ 2] = "0006111"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "20000043053200000"; mesonCutArray[2] = "01631031009000"; //pp LHC11a PHI1 - eventCutArray[ 3] = "0006211"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "20000043053200000"; mesonCutArray[3] = "01631031009000"; //pp LHC11a PHI7 + eventCutArray[ 0] = "0008311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "10000042032030000"; mesonCutArray[0] = "01631031009000"; // EMCEG1, + eventCutArray[ 1] = "0008511"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "10000042032030000"; mesonCutArray[1] = "01631031009000"; // EMCEG2, + eventCutArray[ 2] = "0009311"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "10000042032030000"; mesonCutArray[2] = "01631031009000"; // EMCEJ1, + eventCutArray[ 3] = "0009511"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "10000042032030000"; mesonCutArray[3] = "01631031009000"; // EMCEJ2, + eventCutArray[ 4] = "0000011"; photonCutArray[ 4] = "002093663003800000000"; clusterCutArray[4] = "10000042032030000"; mesonCutArray[4] = "01631031009000"; // INT7 + eventCutArray[ 5] = "0005211"; photonCutArray[ 5] = "002093663003800000000"; clusterCutArray[5] = "10000042032030000"; mesonCutArray[5] = "01631031009000"; // EMC7 + + // ************************************* PHOS cuts **************************************************** + // LHC11a + } else if (trainConfig == 31) { //PHOS clusters + eventCutArray[ 0] = "0000311"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "20000041033200000"; mesonCutArray[0] = "01631031009000"; + eventCutArray[ 1] = "0000311"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[1] = "01631031009000"; + eventCutArray[ 2] = "0000311"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[2] = "01631031009000"; + eventCutArray[ 3] = "0006111"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "20000041033200000"; mesonCutArray[3] = "01631031009000"; + eventCutArray[ 4] = "0006111"; photonCutArray[ 4] = "002093663003800000000"; clusterCutArray[4] = "20000042033200000"; mesonCutArray[4] = "01631031009000"; + eventCutArray[ 5] = "0006111"; photonCutArray[ 5] = "002093663003800000000"; clusterCutArray[5] = "20000043033200000"; mesonCutArray[5] = "01631031009000"; + // LHC13g & LHC12x + } else if (trainConfig == 32) { //PHOS clusters + eventCutArray[ 0] = "0000011"; photonCutArray[ 0] = "002093663003800000000"; clusterCutArray[0] = "20000041033200000"; mesonCutArray[0] = "01631031009000"; + eventCutArray[ 1] = "0000011"; photonCutArray[ 1] = "002093663003800000000"; clusterCutArray[1] = "20000042033200000"; mesonCutArray[1] = "01631031009000"; + eventCutArray[ 2] = "0000011"; photonCutArray[ 2] = "002093663003800000000"; clusterCutArray[2] = "20000043033200000"; mesonCutArray[2] = "01631031009000"; + eventCutArray[ 3] = "0006211"; photonCutArray[ 3] = "002093663003800000000"; clusterCutArray[3] = "20000041033200000"; mesonCutArray[3] = "01631031009000"; + eventCutArray[ 4] = "0006211"; photonCutArray[ 4] = "002093663003800000000"; clusterCutArray[4] = "20000042033200000"; mesonCutArray[4] = "01631031009000"; + eventCutArray[ 5] = "0006211"; photonCutArray[ 5] = "002093663003800000000"; clusterCutArray[5] = "20000043033200000"; mesonCutArray[5] = "01631031009000"; } else { Error(Form("GammaConvCalo_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration"); return;