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);
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(), );
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);
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){
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;
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;
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));
// 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;
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];
// 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";
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];
// 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;
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;
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];
// 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;
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];
// 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;