submitted AddTasks for Dalitz
authorfbock <fbock@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 13 Sep 2013 14:51:18 +0000 (14:51 +0000)
committerfbock <fbock@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 13 Sep 2013 14:51:18 +0000 (14:51 +0000)
PWGGA/GammaConv/AliAnalysisTaskGammaConvV1.cxx
PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb.C
PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb2.C [new file with mode: 0644]
PWGGA/GammaConv/macros/AddTask_gamconv_GammaConvDalitzV1_PbPb.C [new file with mode: 0644]

index 10b809f..dfba6d6 100644 (file)
@@ -1580,9 +1580,10 @@ void AliAnalysisTaskGammaConvV1::CalculatePi0Candidates(){
 //                   cout << "gamma 0: " << gamma0->GetV0Index()<< "\t" << gamma0->GetPx() << "\t" << gamma0->GetPy() << "\t" <<  gamma0->GetPz() << "\t" << endl; 
 //                   cout << "gamma 1: " << gamma1->GetV0Index()<< "\t"<< gamma1->GetPx() << "\t" << gamma1->GetPy() << "\t" <<  gamma1->GetPz() << "\t" << endl; 
 //                    cout << "pi0: "<<fInvMass << "\t" << fPt <<"\t" << fDCAzGammaMin << "\t" << fDCAzGammaMax << "\t" << (Int_t)iFlag << "\t" << (Int_t)iMesonMCInfo <<endl;
-                  if (fIsHeavyIon == 1 && fPt > 0.399 && fPt < 12. ) {
-                     if ( fInvMass > 0.08 && fInvMass < 0.2) tESDMesonsInvMassPtDcazMinDcazMaxFlag[fiCut]->Fill();
-                  } else if (fPt > 0.299 && fPt < 12. )  {
+                  if (fIsHeavyIon == 1 && fPt > 0.399 && fPt < 20. ) {
+                     if (fInvMass > 0.08 && fInvMass < 0.2) tESDMesonsInvMassPtDcazMinDcazMaxFlag[fiCut]->Fill();
+                     if ((fInvMass > 0.45 && fInvMass < 0.6) &&  (fPt > 0.999 && fPt < 20.) )tESDMesonsInvMassPtDcazMinDcazMaxFlag[fiCut]->Fill();
+                  } else if (fPt > 0.299 && fPt < 20. )  {
                      if ( (fInvMass > 0.08 && fInvMass < 0.2) || (fInvMass > 0.45 && fInvMass < 0.6)) tESDMesonsInvMassPtDcazMinDcazMaxFlag[fiCut]->Fill();
                   }   
                }
index 4214082..adfb93f 100644 (file)
@@ -96,457 +96,307 @@ void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set o
    TString *cutarray = new TString[numberOfCuts];
    TString *mesonCutArray = new TString[numberOfCuts];
 
-   if (trainConfig == 1){ // dE/dx E varied to \sigma < +-5
-      cutarray[ 0] = "601000104209197002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000104209197002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000104209197002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000104209197002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 2) { // dE/dx E varied to \sigma < +-5
-      cutarray[ 0] = "524000104209197002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000104209197002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000104209197002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000104209137002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 3) { // dE/dx E varied to \sigma < +-5
-      cutarray[ 0] = "502000104209137002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000104209137002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000104209137002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000104209137002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 4) { // dE/dx E varied to \sigma < +-5
-      cutarray[ 0] = "601000204209197002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000204209197002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000204209197002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000204209197002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 5) { // dE/dx E varied to \sigma < +-5
-      cutarray[ 0] = "524000204209197002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000204209197002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000204209197002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000204209137002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 6) { // dE/dx E varied to \sigma < +-5
-      cutarray[ 0] = "502000204209137002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000204209137002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000204209137002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000204209137002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 7){ // dE/dx E varied to \sigma < -4,5
-      cutarray[ 0] = "601000104209397002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000104209397002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000104209397002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000104209397002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 8) { // dE/dx E varied to \sigma < -4,5
-      cutarray[ 0] = "524000104209397002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000104209397002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000104209397002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000104209337002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 9) { // dE/dx E varied to \sigma < -4,5
-      cutarray[ 0] = "502000104209337002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000104209337002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000104209337002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000104209337002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 10) { // dE/dx E varied to \sigma < -4,5
-      cutarray[ 0] = "601000204209397002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000204209397002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000204209397002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000204209397002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 11) { // dE/dx E varied to \sigma < -4,5
-      cutarray[ 0] = "524000204209397002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000204209397002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000204209397002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000204209337002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 12) { // dE/dx E varied to \sigma < -4,5
-      cutarray[ 0] = "502000204209337002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000204209337002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000204209337002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000204209337002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 13) { // dE/dx E varied to \sigma < +4,-2.5
-      cutarray[ 0] = "601000104209697002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000104209697002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000104209697002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000104209697002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 14) { // dE/dx E varied to \sigma < +4,-2.5
-      cutarray[ 0] = "524000104209697002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000104209697002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000104209697002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000104209637002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 15) {// dE/dx E varied to \sigma < +4,-2.5
-      cutarray[ 0] = "502000104209637002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000104209637002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000104209637002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000104209637002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 16) { // dE/dx E varied to \sigma < +4,-2.5
-      cutarray[ 0] = "601000204209697002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000204209697002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000204209697002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000204209697002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 17) { // dE/dx E varied to \sigma < +4,-2.5
-      cutarray[ 0] = "524000204209697002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000204209697002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000204209697002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000204209637002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 18) {// dE/dx E varied to \sigma < +4,-2.5
-      cutarray[ 0] = "502000204209637002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000204209637002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000204209637002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000204209637002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 19) { // cluster 0.35
-      cutarray[ 0] = "601000104208297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000104208297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000104208297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000104208297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 20) { // cluster 0.35
-      cutarray[ 0] = "524000104208297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000104208297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000104208297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000104208237002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 21) { // cluster 0.35
-      cutarray[ 0] = "502000104208237002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000104208237002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000104208237002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000104208237002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 22) { // cluster 0.35
-      cutarray[ 0] = "601000204208297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000204208297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000204208297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000204208297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 23) { // cluster 0.35
-      cutarray[ 0] = "524000204208297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000204208297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000204208297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000204208237002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 24) { // cluster 0.35
-      cutarray[ 0] = "502000204208237002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000204208237002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000204208237002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000204208237002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 25) { // cluster 0.7
-      cutarray[ 0] = "601000104206297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000104206297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000104206297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000104206297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 26) { // cluster 0.7
-      cutarray[ 0] = "524000104206297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000104206297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000104206297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000104206237002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 27) {  // cluster 0.7
-      cutarray[ 0] = "502000104206237002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000104206237002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000104206237002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000104206237002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 28) {  // cluster 0.7
-      cutarray[ 0] = "601000204206297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000204206297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000204206297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000204206297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 29) {  // cluster 0.7
-      cutarray[ 0] = "524000204206297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000204206297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000204206297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000204206237002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 30) {  // cluster 0.7
-      cutarray[ 0] = "502000204206237002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000204206237002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000204206237002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000204206237002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 31) {  // qt < 0.03
-      cutarray[ 0] = "601000104209297002422000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000104209297002422000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000104209297002422000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000104209297002422000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 32) { // qt < 0.03
-      cutarray[ 0] = "524000104209297002422000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000104209297002422000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000104209297002422000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000104209237002470000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 33) { // qt < 0.03
-      cutarray[ 0] = "502000104209237002470000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000104209237002470000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000104209237002470000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000104209237002470000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 34) { // qt < 0.03
-      cutarray[ 0] = "601000204209297002422000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000204209297002422000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000204209297002422000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000204209297002422000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 35) { // qt < 0.03
-      cutarray[ 0] = "524000204209297002422000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000204209297002422000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000204209297002422000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000204209237002470000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 36) { // qt < 0.03
-      cutarray[ 0] = "502000204209237002470000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000204209237002470000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000204209237002470000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000204209237002470000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%   
-   } else if (trainConfig == 37) {  // qt < 0.07
-      cutarray[ 0] = "601000104209297002222000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000104209297002222000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000104209297002222000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000104209297002222000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 38) { // qt < 0.07
-      cutarray[ 0] = "524000104209297002222000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000104209297002222000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000104209297002222000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000104209237002270000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 39) { // qt < 0.07
-      cutarray[ 0] = "502000104209237002270000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000104209237002270000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000104209237002270000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000104209237002270000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 40) { // qt < 0.07
-      cutarray[ 0] = "601000204209297002222000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000204209297002222000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000204209297002222000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000204209297002222000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 41) { // qt < 0.07
-      cutarray[ 0] = "524000204209297002222000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000204209297002222000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000204209297002222000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000204209237002270000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 42) { // qt < 0.07
-      cutarray[ 0] = "502000204209237002270000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000204209237002270000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000204209237002270000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000204209237002270000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%      
-   } else if (trainConfig == 43) {  // min pt 75
-      cutarray[ 0] = "601000104249297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000104249297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000104249297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000104249297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 44) { // min pt 75
-      cutarray[ 0] = "524000104249297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000104249297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000104249297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000104249237002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 45) { // min pt 75
-      cutarray[ 0] = "502000104249237002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000104249237002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000104249237002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000104249237002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 46) { // min pt 75
-      cutarray[ 0] = "601000204249297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000204249297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000204249297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000204249297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 47) { // min pt 75
-      cutarray[ 0] = "524000204249297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000204249297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000204249297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000204249237002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 48) { // min pt 75
-      cutarray[ 0] = "502000204249237002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000204249237002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000204249237002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000204249237002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%   
-   } else if (trainConfig == 49) {  // min pt 100
-      cutarray[ 0] = "601000104219297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000104219297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000104219297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000104219297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 50) {  // min pt 100
-      cutarray[ 0] = "524000104219297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000104219297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000104219297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000104219237002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 51) {  // min pt 100
-      cutarray[ 0] = "502000104219237002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000104219237002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000104219237002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000104219237002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 52) {  // min pt 100
-      cutarray[ 0] = "601000204219297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000204219297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000204219297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000204219297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 53) {  // min pt 100
-      cutarray[ 0] = "524000204219297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000204219297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000204219297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000204219237002370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 54) {  // min pt 100
-      cutarray[ 0] = "502000204219237002370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000204219237002370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000204219237002370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000204219237002370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%   
-   } else if (trainConfig == 55) {  // TOF \sigma = -3 5
-      cutarray[ 0] = "601000104209297003322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000104209297003322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000104209297003322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000104209297003322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 56) { // TOF \sigma = -3 5
-      cutarray[ 0] = "524000104209297003322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000104209297003322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000104209297003322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000104209237003370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 57) { // TOF \sigma = -3 5
-      cutarray[ 0] = "502000104209237003370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000104209237003370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000104209237003370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000104209237003370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 58) { // TOF \sigma = -3 5
-      cutarray[ 0] = "601000204209297003322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000204209297003322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000204209297003322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000204209297003322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 59) { // TOF \sigma = -3 5
-      cutarray[ 0] = "524000204209297003322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000204209297003322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000204209297003322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000204209237003370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 60) { // TOF \sigma = -3 5
-      cutarray[ 0] = "502000204209237003370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000204209237003370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000204209237003370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000204209237003370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%      
-   } else if (trainConfig == 61) {  // TOF \sigma = -2 3
-      cutarray[ 0] = "601000104209297004322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000104209297004322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000104209297004322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000104209297004322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 62) {  // TOF \sigma = -2 3
-      cutarray[ 0] = "524000104209297004322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000104209297004322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000104209297004322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000104209237004370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 63) {  // TOF \sigma = -2 3
-      cutarray[ 0] = "502000104209237004370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000104209237004370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000104209237004370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000104209237004370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 64) {  // TOF \sigma = -2 3
-      cutarray[ 0] = "601000204209297004322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000204209297004322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000204209297004322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000204209297004322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 65) {  // TOF \sigma = -2 3
-      cutarray[ 0] = "524000204209297004322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000204209297004322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000204209297004322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000204209237004370000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 66) {  // TOF \sigma = -2 3
-      cutarray[ 0] = "502000204209237004370000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000204209237004370000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000204209237004370000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000204209237004370000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%      
-   } else if (trainConfig == 67) {  // alpha central 0.6 (8), peripher 0.75 (5) 
-      cutarray[ 0] = "601000104209297002322000000"; mesonCutArray[ 0] = "01522085009000"; // 0-5%
-      cutarray[ 1] = "612000104209297002322000000"; mesonCutArray[ 1] = "01522085009000"; // 5-10%
-      cutarray[ 2] = "501000104209297002322000000"; mesonCutArray[ 2] = "01522085009000"; // 0-10%
-      cutarray[ 3] = "512000104209297002322000000"; mesonCutArray[ 3] = "01522085009000"; // 10-20%
-   } else if (trainConfig == 68) { // alpha central 0.6 (8), peripher 0.75 (5) 
-      cutarray[ 0] = "524000104209297002322000000"; mesonCutArray[ 0] = "01522085009000"; // 20-40%
-      cutarray[ 1] = "546000104209297002322000000"; mesonCutArray[ 1] = "01522055009000"; // 40-60%
-      cutarray[ 2] = "568000104209297002322000000"; mesonCutArray[ 2] = "01522055009000"; // 60-80%
-      cutarray[ 3] = "524000104209237002370000000"; mesonCutArray[ 3] = "01522085009000"; // 20-40% //direct photon
-   } else if (trainConfig == 69) { // alpha central 0.6 (8), peripher 0.75 (5) 
-      cutarray[ 0] = "502000104209237002370000000"; mesonCutArray[ 0] = "01522085009000"; // 0-20%
-      cutarray[ 1] = "501000104209237002370000000"; mesonCutArray[ 1] = "01522085009000"; // 0-10%
-      cutarray[ 2] = "512000104209237002370000000"; mesonCutArray[ 2] = "01522085009000"; // 10-20%
-      cutarray[ 3] = "548000104209237002370000000"; mesonCutArray[ 3] = "01522055009000"; // 40-80%
-   } else if (trainConfig == 70) { // alpha central 0.6 (8), peripher 0.75 (5) 
-      cutarray[ 0] = "601000204209297002322000000"; mesonCutArray[ 0] = "01522085009000"; // 0-5%
-      cutarray[ 1] = "612000204209297002322000000"; mesonCutArray[ 1] = "01522085009000"; // 5-10%
-      cutarray[ 2] = "501000204209297002322000000"; mesonCutArray[ 2] = "01522085009000"; // 0-10%
-      cutarray[ 3] = "512000204209297002322000000"; mesonCutArray[ 3] = "01522085009000"; // 10-20%
-   } else if (trainConfig == 71) { // alpha central 0.6 (8), peripher 0.75 (5) 
-      cutarray[ 0] = "524000204209297002322000000"; mesonCutArray[ 0] = "01522085009000"; // 20-40%
-      cutarray[ 1] = "546000204209297002322000000"; mesonCutArray[ 1] = "01522055009000"; // 40-60%
-      cutarray[ 2] = "568000204209297002322000000"; mesonCutArray[ 2] = "01522055009000"; // 60-80%
-      cutarray[ 3] = "524000204209237002370000000"; mesonCutArray[ 3] = "01522085009000"; // 20-40% //direct photon
-   } else if (trainConfig == 72) { // alpha central 0.6 (8), peripher 0.75 (5) 
-      cutarray[ 0] = "502000204209237002370000000"; mesonCutArray[ 0] = "01522085009000"; // 0-20%
-      cutarray[ 1] = "501000204209237002370000000"; mesonCutArray[ 1] = "01522085009000"; // 0-10%
-      cutarray[ 2] = "512000204209237002370000000"; mesonCutArray[ 2] = "01522085009000"; // 10-20%
-      cutarray[ 3] = "548000204209237002370000000"; mesonCutArray[ 3] = "01522055009000"; // 40-80%   
-   } else if (trainConfig == 73) {  // alpha central 0.7 (0), peripher 0.85 (7) 
-      cutarray[ 0] = "601000104209297002322000000"; mesonCutArray[ 0] = "01522005009000"; // 0-5%
-      cutarray[ 1] = "612000104209297002322000000"; mesonCutArray[ 1] = "01522005009000"; // 5-10%
-      cutarray[ 2] = "501000104209297002322000000"; mesonCutArray[ 2] = "01522005009000"; // 0-10%
-      cutarray[ 3] = "512000104209297002322000000"; mesonCutArray[ 3] = "01522005009000"; // 10-20%
-   } else if (trainConfig == 74) { // alpha central 0.7 (0), peripher 0.85 (7) 
-      cutarray[ 0] = "524000104209297002322000000"; mesonCutArray[ 0] = "01522005009000"; // 20-40%
-      cutarray[ 1] = "546000104209297002322000000"; mesonCutArray[ 1] = "01522075009000"; // 40-60%
-      cutarray[ 2] = "568000104209297002322000000"; mesonCutArray[ 2] = "01522075009000"; // 60-80%
-      cutarray[ 3] = "524000104209237002370000000"; mesonCutArray[ 3] = "01522005009000"; // 20-40% //direct photon
-   } else if (trainConfig == 75) { // alpha central 0.7 (0), peripher 0.85 (7) 
-      cutarray[ 0] = "502000104209237002370000000"; mesonCutArray[ 0] = "01522005009000"; // 0-20%
-      cutarray[ 1] = "501000104209237002370000000"; mesonCutArray[ 1] = "01522005009000"; // 0-10%
-      cutarray[ 2] = "512000104209237002370000000"; mesonCutArray[ 2] = "01522005009000"; // 10-20%
-      cutarray[ 3] = "548000104209237002370000000"; mesonCutArray[ 3] = "01522075009000"; // 40-80%
-   } else if (trainConfig == 76) { // alpha central 0.7 (0), peripher 0.85 (7) 
-      cutarray[ 0] = "601000204209297002322000000"; mesonCutArray[ 0] = "01522005009000"; // 0-5%
-      cutarray[ 1] = "612000204209297002322000000"; mesonCutArray[ 1] = "01522005009000"; // 5-10%
-      cutarray[ 2] = "501000204209297002322000000"; mesonCutArray[ 2] = "01522005009000"; // 0-10%
-      cutarray[ 3] = "512000204209297002322000000"; mesonCutArray[ 3] = "01522005009000"; // 10-20%
-   } else if (trainConfig == 77) { // alpha central 0.7 (0), peripher 0.85 (7) 
-      cutarray[ 0] = "524000204209297002322000000"; mesonCutArray[ 0] = "01522005009000"; // 20-40%
-      cutarray[ 1] = "546000204209297002322000000"; mesonCutArray[ 1] = "01522075009000"; // 40-60%
-      cutarray[ 2] = "568000204209297002322000000"; mesonCutArray[ 2] = "01522075009000"; // 60-80%
-      cutarray[ 3] = "524000204209237002370000000"; mesonCutArray[ 3] = "01522005009000"; // 20-40% //direct photon
-   } else if (trainConfig == 78) { // alpha central 0.7 (0), peripher 0.85 (7) 
-      cutarray[ 0] = "502000204209237002370000000"; mesonCutArray[ 0] = "01522005009000"; // 0-20%
-      cutarray[ 1] = "501000204209237002370000000"; mesonCutArray[ 1] = "01522005009000"; // 0-10%
-      cutarray[ 2] = "512000204209237002370000000"; mesonCutArray[ 2] = "01522005009000"; // 10-20%
-      cutarray[ 3] = "548000204209237002370000000"; mesonCutArray[ 3] = "01522075009000"; // 40-80%   
-   } else if (trainConfig == 79) {  // chi2 < 20 
-      cutarray[ 0] = "601000104209297002382000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000104209297002382000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000104209297002382000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000104209297002382000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 80) { // chi2 < 20 
-      cutarray[ 0] = "524000104209297002382000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000104209297002382000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000104209297002382000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000104209237002380000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 81) { // chi2 < 20 
-      cutarray[ 0] = "502000104209237002380000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000104209237002380000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000104209237002380000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000104209237002380000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 82) { // chi2 < 20 
-      cutarray[ 0] = "601000204209297002382000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000204209297002382000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000204209297002382000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000204209297002382000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 83) { // chi2 < 20 
-      cutarray[ 0] = "524000204209297002382000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000204209297002382000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000204209297002382000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000204209237002380000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 84) { // chi2 < 20 
-      cutarray[ 0] = "502000204209237002380000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000204209237002380000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000204209237002380000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000204209237002380000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%   
-   } else if (trainConfig == 85) {  // meson: chi2 < 50 (1), gamma: chi2 < 5 (6)
-      cutarray[ 0] = "601000104209297002312000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000104209297002312000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000104209297002312000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000104209297002312000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 86) { // meson: chi2 < 50 (1), gamma: chi2 < 5 (6)
-      cutarray[ 0] = "524000104209297002312000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000104209297002312000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000104209297002312000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000104209237002360000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 87) { // meson: chi2 < 50 (1), gamma: chi2 < 5 (6)
-      cutarray[ 0] = "502000104209237002360000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000104209237002360000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000104209237002360000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000104209237002360000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
-   } else if (trainConfig == 88) { // meson: chi2 < 50 (1), gamma: chi2 < 5 (6)
-      cutarray[ 0] = "601000204209297002312000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
-      cutarray[ 1] = "612000204209297002312000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
-      cutarray[ 2] = "501000204209297002312000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
-      cutarray[ 3] = "512000204209297002312000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
-   } else if (trainConfig == 89) { // meson: chi2 < 50 (1), gamma: chi2 < 5 (6)
-      cutarray[ 0] = "524000204209297002312000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
-      cutarray[ 1] = "546000204209297002312000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
-      cutarray[ 2] = "568000204209297002312000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
-      cutarray[ 3] = "524000204209237002360000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
-   } else if (trainConfig == 90) { // meson: chi2 < 50 (1), gamma: chi2 < 5 (6)
-      cutarray[ 0] = "502000204209237002360000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
-      cutarray[ 1] = "501000204209237002360000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
-      cutarray[ 2] = "512000204209237002360000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
-      cutarray[ 3] = "548000204209237002360000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%   
-
+   if (trainConfig == 1){ 
+      cutarray[ 0] = "601000104209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000104209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000104209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000104209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 2) { 
+      cutarray[ 0] = "524000104209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000104209297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000104209297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000104209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 3) { 
+      cutarray[ 0] = "502000104209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000104209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000104209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000104209297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
+   } else if (trainConfig == 4) { 
+      cutarray[ 0] = "601000204209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000204209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000204209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000204209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 5) { 
+      cutarray[ 0] = "524000204209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000204209297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000204209297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000204209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 6) { 
+      cutarray[ 0] = "502000204209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000204209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000204209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000204209297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
+   } else if (trainConfig == 7){ // Psi pair 10000
+      cutarray[ 0] = "601000104209297002320000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000104209297002320000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000104209297002320000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000104209297002320000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 8) { // Psi pair 10000
+      cutarray[ 0] = "524000104209297002320000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000104209297002320000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000104209297002320000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000104209297002320000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 9) { // Psi pair 10000
+      cutarray[ 0] = "502000104209297002320000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000104209297002320000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000104209297002320000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000104209297002320000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
+   } else if (trainConfig == 10) { // Psi pair 10000
+      cutarray[ 0] = "601000204209297002320000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000204209297002320000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000204209297002320000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000204209297002320000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 11) { // Psi pair 10000
+      cutarray[ 0] = "524000204209297002320000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000204209297002320000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000204209297002320000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000204209297002320000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 12) { // Psi pair 10000
+      cutarray[ 0] = "502000204209297002320000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000204209297002320000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000204209297002320000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000204209297002320000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%   
+   } else if (trainConfig == 13){ // Psi pair 0.1
+      cutarray[ 0] = "601000104209297002321000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000104209297002321000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000104209297002321000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000104209297002321000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 14) { // Psi pair 0.1
+      cutarray[ 0] = "524000104209297002321000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000104209297002321000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000104209297002321000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000104209297002321000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 15) { // Psi pair 0.1
+      cutarray[ 0] = "502000104209297002321000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000104209297002321000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000104209297002321000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000104209297002321000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
+   } else if (trainConfig == 16) { // Psi pair 0.1
+      cutarray[ 0] = "601000204209297002321000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000204209297002321000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000204209297002321000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000204209297002321000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 17) { // Psi pair 0.1
+      cutarray[ 0] = "524000204209297002321000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000204209297002321000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000204209297002321000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000204209297002321000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 18) { // Psi pair 0.1
+      cutarray[ 0] = "502000204209297002321000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000204209297002321000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000204209297002321000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000204209297002321000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%      
+   } else if (trainConfig == 19){ // Psi pair 0.035
+      cutarray[ 0] = "601000104209297002323000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000104209297002323000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000104209297002323000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000104209297002323000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 20) { // Psi pair 0.035
+      cutarray[ 0] = "524000104209297002323000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000104209297002323000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000104209297002323000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000104209297002323000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 21) { // Psi pair 0.035
+      cutarray[ 0] = "502000104209297002323000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000104209297002323000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000104209297002323000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000104209297002323000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
+   } else if (trainConfig == 22) { // Psi pair 0.035
+      cutarray[ 0] = "601000204209297002323000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000204209297002323000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000204209297002323000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000204209297002323000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 23) { // Psi pair 0.035
+      cutarray[ 0] = "524000204209297002323000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000204209297002323000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000204209297002323000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000204209297002323000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 24) { // Psi pair 0.035
+      cutarray[ 0] = "502000204209297002323000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000204209297002323000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000204209297002323000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000204209297002323000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%         
+   } else if (trainConfig == 25){ // Eta 0.75
+      cutarray[ 0] = "601000103209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000103209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000103209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000103209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 26) { // Eta 0.75
+      cutarray[ 0] = "524000103209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000103209297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000103209297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000103209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 27) { // Eta 0.75
+      cutarray[ 0] = "502000103209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000103209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000103209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000103209297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
+   } else if (trainConfig == 28) { // Eta 0.75
+      cutarray[ 0] = "601000203209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000203209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000203209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000203209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 29) { // Eta 0.75
+      cutarray[ 0] = "524000203209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000203209297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000203209297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000203209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 30) { // Eta 0.75
+      cutarray[ 0] = "502000203209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000203209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000203209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000203209297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%         
+   } else if (trainConfig == 31){ // Eta 0.8
+      cutarray[ 0] = "601000100209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000100209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000100209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000100209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 32) { // Eta 0.8
+      cutarray[ 0] = "524000100209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000100209297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000100209297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000100209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 33) { // Eta 0.8
+      cutarray[ 0] = "502000100209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000100209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000100209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000100209297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
+   } else if (trainConfig == 34) { // Eta 0.8
+      cutarray[ 0] = "601000200209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000200209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000200209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000200209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 35) { // Eta 0.8
+      cutarray[ 0] = "524000200209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000200209297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000200209297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000200209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 36) { // Eta 0.8
+      cutarray[ 0] = "502000200209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000200209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000200209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000200209297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%         
+   } else if (trainConfig == 37){ // Eta 0.6
+      cutarray[ 0] = "601000101209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000101209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000101209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000101209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 38) { // Eta 0.6
+      cutarray[ 0] = "524000101209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000101209297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000101209297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000101209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 39) { // Eta 0.6
+      cutarray[ 0] = "502000101209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000101209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000101209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000101209297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
+   } else if (trainConfig == 40) { // Eta 0.6
+      cutarray[ 0] = "601000201209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000201209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000201209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000201209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 41) { // Eta 0.6
+      cutarray[ 0] = "524000201209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000201209297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000201209297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000201209297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 42) { // Eta 0.6
+      cutarray[ 0] = "502000201209297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000201209297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000201209297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000201209297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
+   } else if (trainConfig == 43){ // R 10-180
+      cutarray[ 0] = "601000104509297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000104509297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000104509297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000104509297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 44) {  // R 10-180
+      cutarray[ 0] = "524000104509297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000104509297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000104509297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000104509297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 45) {  // R 10-180
+      cutarray[ 0] = "502000104509297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000104509297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000104509297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000104509297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
+   } else if (trainConfig == 46) {  // R 10-180
+      cutarray[ 0] = "601000204509297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000204509297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000204509297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000204509297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 47) {  // R 10-180
+      cutarray[ 0] = "524000204509297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000204509297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000204509297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000204509297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 48) {  // R 10-180
+      cutarray[ 0] = "502000204509297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000204509297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000204509297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000204509297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%  
+    } else if (trainConfig == 49){ // R 2.8-180
+      cutarray[ 0] = "601000104109297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000104109297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000104109297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000104109297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 50) {  // R 2.8-180
+      cutarray[ 0] = "524000104109297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000104109297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000104109297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000104109297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 51) {  // R 2.8-180
+      cutarray[ 0] = "502000104109297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000104109297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000104109297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000104109297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
+   } else if (trainConfig == 52) {  // R 2.8-180
+      cutarray[ 0] = "601000204109297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000204109297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000204109297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000204109297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 53) {  // R 2.8-180
+      cutarray[ 0] = "524000204109297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000204109297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000204109297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000204109297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 54) {  // R 2.8-180
+      cutarray[ 0] = "502000204109297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000204109297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000204109297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000204109297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%  
+    } else if (trainConfig == 55){ // R 2.8-180
+      cutarray[ 0] = "601000104609297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000104609297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000104609297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000104609297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 56) {  // R 2.8-180
+      cutarray[ 0] = "524000104609297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000104609297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000104609297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000104609297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 57) {  // R 2.8-180
+      cutarray[ 0] = "502000104609297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000104609297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000104609297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000104609297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
+   } else if (trainConfig == 58) {  // R 2.8-180
+      cutarray[ 0] = "601000204609297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
+      cutarray[ 1] = "612000204609297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
+      cutarray[ 2] = "501000204609297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
+      cutarray[ 3] = "512000204609297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
+   } else if (trainConfig == 59) {  // R 2.8-180
+      cutarray[ 0] = "524000204609297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
+      cutarray[ 1] = "546000204609297002322000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
+      cutarray[ 2] = "568000204609297002322000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
+      cutarray[ 3] = "524000204609297002322000000"; mesonCutArray[ 3] = "01522045009000"; // 20-40% //direct photon
+   } else if (trainConfig == 60) {  // R 2.8-180
+      cutarray[ 0] = "502000204609297002322000000"; mesonCutArray[ 0] = "01522045009000"; // 0-20%
+      cutarray[ 1] = "501000204609297002322000000"; mesonCutArray[ 1] = "01522045009000"; // 0-10%
+      cutarray[ 2] = "512000204609297002322000000"; mesonCutArray[ 2] = "01522045009000"; // 10-20%
+      cutarray[ 3] = "548000204609297002322000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%  
+     
    } else {
       Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
       return;
diff --git a/PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb2.C b/PWGGA/GammaConv/macros/AddTask_GammaConvV1_PbPb2.C
new file mode 100644 (file)
index 0000000..8385ed8
--- /dev/null
@@ -0,0 +1,179 @@
+void AddTask_GammaConvV1_PbPb2(  Int_t trainConfig = 1,  //change different set of cuts
+                              Bool_t isMC   = kFALSE, //run MC 
+                              Bool_t enableQAMesonTask = kFALSE, //enable QA in AliAnalysisTaskGammaConvV1
+                              Bool_t enableQAPhotonTask = kFALSE, // enable additional QA task
+                              TString fileNameInputForWeighting = "MCSpectraInput.root", // path to file for weigting input
+                              Bool_t doWeighting = kFALSE,  //enable Weighting
+                              TString cutnumberAODBranch = "1000000060084000001500000" 
+                           ) {
+
+   // ================= Load Librariers =================================
+   gSystem->Load("libCore.so");  
+   gSystem->Load("libTree.so");
+   gSystem->Load("libGeom.so");
+   gSystem->Load("libVMC.so");
+   gSystem->Load("libPhysics.so");
+   gSystem->Load("libMinuit");
+   gSystem->Load("libSTEERBase");
+   gSystem->Load("libESD");
+   gSystem->Load("libAOD");
+   gSystem->Load("libANALYSIS");
+   gSystem->Load("libANALYSISalice");  
+   gSystem->Load("libPWGGAGammaConv.so");
+   gSystem->Load("libCDB.so");
+   gSystem->Load("libSTEER.so");
+   gSystem->Load("libSTEERBase.so");
+   gSystem->Load("libTENDER.so");
+   gSystem->Load("libTENDERSupplies.so");
+      
+   // ================== GetAnalysisManager ===============================
+   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+   if (!mgr) {
+      Error(Form("AddTask_GammaConvV1_%i",trainConfig), "No analysis manager found.");
+      return ;
+   }
+
+   // ================== GetInputEventHandler =============================
+   AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
+   
+   //========= Add PID Reponse to ANALYSIS manager ====
+   if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){
+      gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
+      AddTaskPIDResponse(isMC);
+   }
+   
+   //=========  Set Cutnumber for V0Reader ================================
+   TString cutnumber = "100000000008400100150000000"; 
+   AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
+   //========= Add V0 Reader to  ANALYSIS manager if not yet existent =====
+   if( !(AliV0ReaderV1*)mgr->GetTask("V0ReaderV1") ){
+      AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1");
+      
+      fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE);
+      fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output
+      fV0ReaderV1->SetUseAODConversionPhoton(kTRUE);
+
+      if (!mgr) {
+         Error("AddTask_V0ReaderV1", "No analysis manager found.");
+         return;
+      }
+
+      // Set AnalysisCut Number
+      AliConversionCuts *fCuts=NULL;
+      if(cutnumber!=""){
+         fCuts= new AliConversionCuts(cutnumber.Data(),cutnumber.Data());
+         fCuts->SetPreSelectionCutFlag(kTRUE);
+         if(fCuts->InitializeCutsFromCutString(cutnumber.Data())){
+            fV0ReaderV1->SetConversionCuts(fCuts);
+            fCuts->SetFillCutHistograms("",kTRUE);
+         }
+      }
+      
+      if(inputHandler->IsA()==AliAODInputHandler::Class()){
+      // AOD mode
+         fV0ReaderV1->SetDeltaAODBranchName(Form("GammaConv_%s_gamma",cutnumberAODBranch.Data()));
+      }
+      fV0ReaderV1->Init();
+
+      AliLog::SetGlobalLogLevel(AliLog::kInfo);
+
+      //connect input V0Reader
+      mgr->AddTask(fV0ReaderV1);
+      mgr->ConnectInput(fV0ReaderV1,0,cinput);
+
+   }
+
+   //================================================
+   //========= Add task to the ANALYSIS manager =====
+   //================================================
+   AliAnalysisTaskGammaConvV1 *task=NULL;
+   task= new AliAnalysisTaskGammaConvV1(Form("GammaConvV1_%i",trainConfig));
+   task->SetIsHeavyIon(1);
+   task->SetIsMC(isMC);
+   // Cut Numbers to use in Analysis
+   Int_t numberOfCuts = 1;
+
+   TString *cutarray = new TString[numberOfCuts];
+   TString *mesonCutArray = new TString[numberOfCuts];
+
+   if (trainConfig == 1){ 
+      cutarray[ 0] = "601000104209297002322000000"; mesonCutArray[ 0] = "01522045009000"; 
+   } else if (trainConfig == 2) { 
+      cutarray[ 0] = "612000104209297002322000000"; mesonCutArray[ 0] = "01522045009000"; 
+   } else if (trainConfig == 3) { 
+      cutarray[ 0] = "501000104209297002322000000"; mesonCutArray[ 0] = "01522045009000"; 
+   } else if (trainConfig == 4) { 
+      cutarray[ 0] = "502000104209297002322000000"; mesonCutArray[ 0] = "01522045009000";    
+   } else if (trainConfig == 5) { 
+      cutarray[ 0] = "512000104209297002322000000"; mesonCutArray[ 0] = "01522045009000";    
+   } else if (trainConfig == 6) { 
+      cutarray[ 0] = "524000104209297002322000000"; mesonCutArray[ 0] = "01522045009000";       
+   } else if (trainConfig == 7) {    
+      cutarray[ 0] = "546000104209297002322000000"; mesonCutArray[ 0] = "01522065009000"; 
+   } else if (trainConfig == 8) {    
+      cutarray[ 0] = "548000104209297002322000000"; mesonCutArray[ 0] = "01522065009000";    
+   } else if (trainConfig == 9) {    
+      cutarray[ 0] = "545000104209297002322000000"; mesonCutArray[ 0] = "01522065009000"; 
+   } else if (trainConfig == 10) { 
+      cutarray[ 0] = "556000104209297002322000000"; mesonCutArray[ 0] = "01522065009000";
+   } else if (trainConfig == 11) { 
+      cutarray[ 0] = "568000104209297002322000000"; mesonCutArray[ 0] = "01522065009000";    
+   } else if (trainConfig == 12) { 
+      cutarray[ 0] = "567000104209297002322000000"; mesonCutArray[ 0] = "01522065009000"; 
+   } else if (trainConfig == 13) { 
+      cutarray[ 0] = "578000104209297002322000000"; mesonCutArray[ 0] = "01522065009000"; 
+   } else if (trainConfig == 14) { 
+      cutarray[ 0] = "469000104209297002322000000"; mesonCutArray[ 0] = "01522065009000";
+   } else if (trainConfig == 15) { 
+      cutarray[ 0] = "589000104209297002322000000"; mesonCutArray[ 0] = "01522065009000";    
+   } else {
+      Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
+      return;
+   }
+
+   TList *ConvCutList = new TList();
+   TList *MesonCutList = new TList();
+
+   TList *HeaderList = new TList();
+   TObjString *Header1 = new TObjString("pi0_1");
+   HeaderList->Add(Header1);
+//    TObjString *Header3 = new TObjString("eta_2");
+//    HeaderList->Add(Header3);
+   
+   ConvCutList->SetOwner(kTRUE);
+   AliConversionCuts **analysisCuts = new AliConversionCuts*[numberOfCuts];
+   MesonCutList->SetOwner(kTRUE);
+   AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts];
+
+   for(Int_t i = 0; i<numberOfCuts; i++){
+      analysisCuts[i] = new AliConversionCuts();
+      analysisCuts[i]->InitializeCutsFromCutString(cutarray[i].Data());
+      ConvCutList->Add(analysisCuts[i]);
+
+      analysisCuts[i]->SetFillCutHistograms("",kFALSE);
+      analysisMesonCuts[i] = new AliConversionMesonCuts();
+      analysisMesonCuts[i]->InitializeCutsFromCutString(mesonCutArray[i].Data());
+      MesonCutList->Add(analysisMesonCuts[i]);
+      analysisMesonCuts[i]->SetFillCutHistograms("");
+      analysisCuts[i]->SetAcceptedHeader(HeaderList);
+   }
+
+   task->SetConversionCutList(numberOfCuts,ConvCutList);
+   task->SetMesonCutList(numberOfCuts,MesonCutList);
+   task->SetMoveParticleAccordingToVertex(kTRUE);
+   task->SetDoMesonAnalysis(kTRUE);
+   if (enableQAMesonTask) task->SetDoMesonQA(kTRUE); //Attention new switch for Pi0 QA
+   if (enableQAPhotonTask) task->SetDoPhotonQA(kTRUE);  //Attention new switch small for Photon QA
+
+   //connect containers
+   AliAnalysisDataContainer *coutput =
+      mgr->CreateContainer(Form("GammaConvV1_%i",trainConfig), TList::Class(),
+                           AliAnalysisManager::kOutputContainer,Form("GammaConvV1_%i.root",trainConfig));
+
+   mgr->AddTask(task);
+   mgr->ConnectInput(task,0,cinput);
+   mgr->ConnectOutput(task,1,coutput);
+
+   return;
+
+}
diff --git a/PWGGA/GammaConv/macros/AddTask_gamconv_GammaConvDalitzV1_PbPb.C b/PWGGA/GammaConv/macros/AddTask_gamconv_GammaConvDalitzV1_PbPb.C
new file mode 100644 (file)
index 0000000..960d350
--- /dev/null
@@ -0,0 +1,267 @@
+void AddTask_gamconv_GammaConvDalitzV1_PbPb(   Bool_t isMC   = kFALSE, //run MC 
+                                       Bool_t enableQAMesonTask = kTRUE, //enable QA in AliAnalysisTaskGammaConvDalitzV1
+                                       Bool_t enableDoMesonChic = kFALSE, // enable additional Chic analysis
+                                       TString cutnumberAODBranch = "0000000060084001001500000"
+                                  ) {
+
+
+   
+   cout<<"Entro -1"<<endl;
+
+  // ================= Load Librariers =================================
+   gSystem->Load("libCore.so");  
+   gSystem->Load("libTree.so");
+   gSystem->Load("libGeom.so");
+   gSystem->Load("libVMC.so");
+   gSystem->Load("libPhysics.so");
+   gSystem->Load("libMinuit");
+   gSystem->Load("libSTEERBase");
+   gSystem->Load("libESD");
+   gSystem->Load("libAOD");
+   gSystem->Load("libANALYSIS");
+   gSystem->Load("libANALYSISalice");  
+   gSystem->Load("libPWGGAGammaConv.so");
+   gSystem->Load("libCDB.so");
+   gSystem->Load("libSTEER.so");
+   gSystem->Load("libSTEERBase.so");
+   gSystem->Load("libTENDER.so");
+   gSystem->Load("libTENDERSupplies.so");
+
+
+   cout<<"Entro 0"<<endl;
+
+   // ================== GetAnalysisManager ===============================
+   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+   if (!mgr) {
+      Error("AddTask_GammaConvDalitzV1_PbPb", "No analysis manager found.");
+      return ;
+   }
+
+   // ================== GetInputEventHandler =============================
+   AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
+   
+   //========= Add PID Reponse to ANALYSIS manager ====
+   if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){
+      gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
+      AddTaskPIDResponse(isMC);
+   }
+   
+   //=========  Set Cutnumber for V0Reader ================================
+   TString ConvCutnumber = "108000000008400100150000000";   //Online  V0 finder
+   TString ElecCuts      = "900054000000020000";            //Electron Cuts
+
+
+
+   AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
+   
+   //========= Add V0 Reader to  ANALYSIS manager if not yet existent =====
+   if( !(AliV0ReaderV1*)mgr->GetTask("V0ReaderV1") ){
+      AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1");
+      
+      fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE);
+      fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output
+      fV0ReaderV1->SetUseAODConversionPhoton(kTRUE);
+      
+      if (!mgr) {
+         Error("AddTask_V0ReaderV1", "No analysis manager found.");
+         return;
+      }
+
+      // Set AnalysisCut Number
+      AliConversionCuts *fCuts=NULL;
+      if( ConvCutnumber !=""){
+         fCuts= new AliConversionCuts(ConvCutnumber.Data(),ConvCutnumber.Data());
+         fCuts->SetPreSelectionCutFlag(kTRUE);
+         if(fCuts->InitializeCutsFromCutString(ConvCutnumber.Data())){
+            fV0ReaderV1->SetConversionCuts(fCuts);
+            fCuts->SetFillCutHistograms("",kTRUE);
+         }
+      }
+      if(inputHandler->IsA()==AliAODInputHandler::Class()){
+      // AOD mode
+         fV0ReaderV1->SetDeltaAODBranchName(Form("GammaConv_%s_gamma",cutnumberAODBranch.Data()));
+      }
+      fV0ReaderV1->Init();
+
+      AliLog::SetGlobalLogLevel(AliLog::kInfo);
+
+      //connect input V0Reader
+      mgr->AddTask(fV0ReaderV1);
+      mgr->ConnectInput(fV0ReaderV1,0,cinput);
+
+   }
+
+   //================================================
+   //========= Add Electron Selector ================
+
+
+   if( !(AliDalitzElectronSelector*)mgr->GetTask("ElectronSelector") ){
+
+   AliDalitzElectronSelector *fElectronSelector = new AliDalitzElectronSelector("ElectronSelector");
+
+   // Set AnalysisCut Number
+
+   AliDalitzElectronCuts *fElecCuts=0;
+
+   TString ElecCuts = "900054000000020000";
+
+
+
+    if( ElecCuts!=""){
+
+       fElecCuts= new AliDalitzElectronCuts(ElecCuts.Data(),ElecCuts.Data());
+
+            if(fElecCuts->InitializeCutsFromCutString(ElecCuts.Data())){
+
+                fElectronSelector->SetDalitzElectronCuts(fElecCuts);
+
+                fElecCuts->SetFillCutHistograms("",kTRUE);
+
+            }
+
+    }
+
+    fElectronSelector->Init();
+    mgr->AddTask(fElectronSelector);
+
+    }
+
+
+
+    cout<<"Entro"<<endl;
+   //================================================
+   //========= Add task to the ANALYSIS manager =====
+   //================================================
+   //            find input container
+   
+  
+   AliAnalysisTaskGammaConvDalitzV1 *task=NULL;
+
+   task= new AliAnalysisTaskGammaConvDalitzV1("GammaConvDalitzV1");
+
+   task->SetIsHeavyIon(2);
+   task->SetIsMC(isMC);
+
+
+
+   // Cut Numbers to use in Analysis
+   Int_t numberOfCuts = 5;
+
+   TString *ConvCutarray    = new TString[numberOfCuts];
+
+   TString *ElecCutarray    = new TString[numberOfCuts];
+
+   TString *MesonCutarray   = new TString[numberOfCuts];
+
+   Bool_t doEtaShiftIndCuts = kFALSE;
+
+
+
+
+   TString *cutarray = new TString[numberOfCuts];
+   TString *mesonCutArray = new TString[numberOfCuts];
+   Bool_t doEtaShiftIndCuts = kFALSE;
+   TString stringShift = "";
+
+   // Shifting in pPb direction
+
+   doEtaShiftIndCuts = kTRUE;
+   stringShift = "pPb";
+
+   ConvCutarray[0]  = "312000104209297100322000000"; MesonCutarray[0] = "01522045009000"; ElecCutarray[0]  = "905162001552021217"; //PbPb  5-10%
+   ConvCutarray[1]  = "112000104209297100322000000"; MesonCutarray[1] = "01522045009000"; ElecCutarray[1]  = "905162001552021217"; //PbPb 10-20%
+   ConvCutarray[2]  = "124000104209297100322000000"; MesonCutarray[2] = "01522045009000"; ElecCutarray[2]  = "905162001552021217"; //PbPb 20-40%
+   ConvCutarray[3]  = "146000104209297100322000000"; MesonCutarray[3] = "01522045009000"; ElecCutarray[3]  = "905162001552021217"; //PbPb 40-60%
+   ConvCutarray[4]  = "168000104209297100322000000"; MesonCutarray[4] = "01522045009000"; ElecCutarray[4]  = "905162001552021217"; //PbPb 60-80%
+
+
+   TList *ConvCutList  = new TList();
+   TList *MesonCutList = new TList();
+   TList *ElecCutList  = new TList();
+
+   TList *HeaderList = new TList();
+   TObjString *Header1 = new TObjString("pi0_1");
+   HeaderList->Add(Header1);
+   //TObjString *Header3 = new TObjString("eta_2");
+   //HeaderList->Add(Header3);
+   
+   ConvCutList->SetOwner(kTRUE);
+   AliConversionCuts **analysisCuts             = new AliConversionCuts*[numberOfCuts];
+   MesonCutList->SetOwner(kTRUE);
+   AliConversionMesonCuts **analysisMesonCuts   = new AliConversionMesonCuts*[numberOfCuts];
+   ElecCutList->SetOwner(kTRUE);
+   AliDalitzElectronCuts **analysisElecCuts     = new AliDalitzElectronCuts*[numberOfCuts];
+
+
+
+   for(Int_t i = 0; i<numberOfCuts; i++){
+
+
+      analysisCuts[i] = new AliConversionCuts();
+      if( ! analysisCuts[i]->InitializeCutsFromCutString(ConvCutarray[i].Data()) ) {
+            cout<<"ERROR: analysisCuts [" <<i<<"]"<<endl;
+            return 0;
+      }
+      else {
+      
+      ConvCutList->Add(analysisCuts[i]);
+      analysisCuts[i]->SetFillCutHistograms("",kFALSE);
+      }
+
+
+
+      analysisMesonCuts[i] = new AliConversionMesonCuts();
+    
+      if( ! analysisMesonCuts[i]->InitializeCutsFromCutString(MesonCutarray[i].Data()) ) {
+            cout<<"ERROR: analysisMesonCuts [ " <<i<<" ] "<<endl;
+            return 0;
+      }
+      else {
+        MesonCutList->Add(analysisMesonCuts[i]);
+        analysisMesonCuts[i]->SetFillCutHistograms("");
+      }
+
+
+       TString cutName( Form("%s_%s_%s",ConvCutarray[i].Data(),ElecCutarray[i].Data(),MesonCutarray[i].Data() ) );
+
+
+       analysisElecCuts[i] = new AliDalitzElectronCuts();
+       if( !analysisElecCuts[i]->InitializeCutsFromCutString(ElecCutarray[i].Data())) {
+
+            cout<< "ERROR:  analysisElecCuts [ " <<i<<" ] "<<endl;
+            return 0;
+       }
+       else { 
+        ElecCutList->Add(analysisElecCuts[i]);
+        analysisElecCuts[i]->SetFillCutHistograms("",kFALSE,cutName); 
+       }
+
+
+        analysisCuts[i]->SetAcceptedHeader(HeaderList);
+
+   }
+
+
+   task->SetConversionCutList(numberOfCuts,ConvCutList);
+   task->SetMesonCutList(MesonCutList);
+   task->SetElectronCutList(ElecCutList);
+
+   task->SetMoveParticleAccordingToVertex(kTRUE);
+
+
+   if(enableQAMesonTask) task->SetDoMesonQA(kTRUE);
+   if(enableDoMesonChic) task->SetDoMesonAnalysis(kTRUE);
+
+   //connect containers
+   AliAnalysisDataContainer *coutput =
+   mgr->CreateContainer("GammaConvDalitzV1", TList::Class(),
+                           AliAnalysisManager::kOutputContainer,"GammaConvDalitzV1.root");
+
+   mgr->AddTask(task);
+   mgr->ConnectInput(task,0,cinput);
+   mgr->ConnectOutput(task,1,coutput);
+
+   return;
+
+}