- changes in AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx to run on MC and modific...

diff --git a/PWGGA/GammaConv/AliAODConversionMother.cxx b/PWGGA/GammaConv/AliAODConversionMother.cxx
index 4f2a15075b39a9f6107a9516f010b1f4fece8d50..3ed962be17088144d6df46f0139b8decd53cb9cf 100644 (file)
--- a/PWGGA/GammaConv/AliAODConversionMother.cxx
+++ b/PWGGA/GammaConv/AliAODConversionMother.cxx
@@ -19,7 +19,8 @@ AliAODConversionParticle(),
     fdcaBetweenPhotons(1),
     fdcaZPrimVtx(100),
     fdcaRPrimVtx(100),
-    fQuality(0)
+    fQuality(0),
+    fTrueMeson(0)
    
 {
        fLabel[0] = -1;
@@ -42,7 +43,8 @@ fWeight(1),
 fdcaBetweenPhotons(100),
 fdcaZPrimVtx(100),
 fdcaRPrimVtx(100),
-fQuality(0)
+fQuality(0),
+fTrueMeson(0)
 {
     // Set 4momentu
     SetPxPyPzE(kf->GetPx(),kf->GetPy(),kf->GetPz(),kf->GetE());
@@ -67,7 +69,8 @@ fWeight(1),
 fdcaBetweenPhotons(1),
 fdcaZPrimVtx(100),
 fdcaRPrimVtx(100),
-fQuality(0)
+fQuality(0),
+fTrueMeson(0)
 {
     // Set 4momentum
     SetPxPyPzE(y1->Px()+y2->Px(),y1->Py()+y2->Py(),y1->Pz()+y2->Pz(),y1->E()+y2->E());
@@ -80,7 +83,7 @@ fQuality(0)
     DetermineMesonQuality(y1,y2);
     // Calculate Alpha
     if((y1->E()+y2->E()) != 0){
-       fAlpha=TMath::Abs((y1->E()-y2->E())/(y1->E()+y2->E()));
+               fAlpha=TMath::Abs((y1->E()-y2->E())/(y1->E()+y2->E()));
     }
 
     // Set Chi2 to the mean chi2 of gammas
@@ -102,7 +105,8 @@ fWeight(1),
 fdcaBetweenPhotons(1),
 fdcaZPrimVtx(100),
 fdcaRPrimVtx(100),
-fQuality(0)
+fQuality(0),
+fTrueMeson(0)
 {
     // Set 4momentum
     SetPxPyPzE(meson->Px()+gamma->Px(),meson->Py()+gamma->Py(),meson->Pz()+gamma->Pz(),meson->E()+gamma->E());
@@ -122,7 +126,7 @@ fQuality(0)
     }
 
     // Set Chi2 to the mean chi2 of gammas
- //   fChi2=0.5*(y1->GetChi2perNDF()+y2->GetChi2perNDF());
+       // fChi2=0.5*(y1->GetChi2perNDF()+y2->GetChi2perNDF());
 
     //Set Decay Photon Labels
     fLabel[0]=-1;

diff --git a/PWGGA/GammaConv/AliAODConversionMother.h b/PWGGA/GammaConv/AliAODConversionMother.h
index de2a99d11f79592ee0d0cf4a751958583598e751..a5634876435756d1b70379f926a3119f4983e2c6 100644 (file)
--- a/PWGGA/GammaConv/AliAODConversionMother.h
+++ b/PWGGA/GammaConv/AliAODConversionMother.h
@@ -89,7 +89,11 @@ class AliAODConversionMother : public AliAODConversionParticle{
                Float_t CalculateDistanceBetweenPhotons(AliAODConversionPhoton* y1, AliAODConversionPhoton* y2 , Double_t prodPoint[3]);
                void CalculateDistanceOfClossetApproachToPrimVtx(const AliVVertex* primVertex);
                void DetermineMesonQuality(AliAODConversionPhoton* y1, AliAODConversionPhoton* y2);
-                       
+               
+               void SetTrueMesonValue(Int_t trueMeson) {fTrueMeson = trueMeson;}
+               Int_t GetTrueMesonValue()const {return fTrueMeson;}
+               
+               
        private:
                Int_t fLabel[3];                                                // Labels of the decay photons
                Int_t fMCLabel;                                                 // MC Label
@@ -110,8 +114,14 @@ class AliAODConversionMother : public AliAODConversionParticle{
                                                                                                        // 5: 1 photon quality 2, 1 photon quality 3
                                                                                                        // 6: both photons quality 3
                
+               Int_t fTrueMeson;                                               // is true meson
+                                                                                                       // 0 : no
+                                                                                                       // 1 : pi0
+                                                                                                       // 2 : eta
+                                                                                                       // 3 : eta'
+                                                                                                       // 4 : omega
                
-    ClassDef(AliAODConversionMother,4)
+    ClassDef(AliAODConversionMother,5)
 };
 
 #endif

diff --git a/PWGGA/GammaConv/AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx b/PWGGA/GammaConv/AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx
index 197955d73060e60d72fec665d4dace4b280639c3..075171b0122516458be931a9a761d1e420dcbb27 100644 (file)
--- a/PWGGA/GammaConv/AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx
+++ b/PWGGA/GammaConv/AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero.cxx
@@ -951,6 +951,7 @@ void AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueNeutralPionCandidat
                Bool_t gamma1DalitzCand = kFALSE;
                Int_t gamma0MCLabel = TrueGammaCandidate0->GetMCParticleLabel(MCStack);
                Int_t gamma0MotherLabel = -1;
+               Int_t motherRealLabel = -1;
                if(gamma0MCLabel != -1){ // Gamma is Combinatorial; MC Particles don't belong to the same Mother
                        // Daughters Gamma 0
                        TParticle * negativeMC = (TParticle*)TrueGammaCandidate0->GetNegativeMCDaughter(MCStack);
@@ -960,6 +961,7 @@ void AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueNeutralPionCandidat
                                if(negativeMC->GetUniqueID() == 5 && positiveMC->GetUniqueID() ==5){ // ... From Conversion ...
                                        if(gammaMC0->GetPdgCode() == 22){ // ... with Gamma Mother
                                                gamma0MotherLabel=gammaMC0->GetFirstMother();
+                                               motherRealLabel=gammaMC0->GetFirstMother();
                                        }
                                }
                                if(gammaMC0->GetPdgCode() ==111){ // Dalitz candidate
@@ -1006,6 +1008,8 @@ void AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueNeutralPionCandidat
                        
                        
                        if(isTruePi0 || isTruePi0Dalitz){// True Pion 
+                               Pi0Candidate->SetTrueMesonValue(1);
+                               Pi0Candidate->SetMCLabel(motherRealLabel);
                                fHistoTrueMotherGammaGammaInvMassPt[fiCut]->Fill(Pi0Candidate->M(),Pi0Candidate->Pt()); 
                        }
                }       
@@ -1021,7 +1025,8 @@ void AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueNeutralPionCandidat
        Bool_t isTruePi0Dalitz = kFALSE;
        Bool_t gamma0DalitzCand = kFALSE;
        Bool_t gamma1DalitzCand = kFALSE;
-   
+       Int_t motherRealLabel = -1;
+               
        if (AODMCTrackArray!=NULL && TrueGammaCandidate0 != NULL){
                AliAODMCParticle *positiveMC = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(TrueGammaCandidate0->GetMCLabelPositive()));
                AliAODMCParticle *negativeMC = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(TrueGammaCandidate0->GetMCLabelNegative()));
@@ -1041,7 +1046,8 @@ void AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueNeutralPionCandidat
                        if(abs(negativeMC->GetPdgCode())==11 && abs(positiveMC->GetPdgCode())==11){  // Electrons ...
                                if(((positiveMC->GetMCProcessCode())) == 5 && ((negativeMC->GetMCProcessCode())) == 5){ // ... From Conversion ...     
                                        if(gammaMC0->GetPdgCode() == 22){ // ... with Gamma Mother
-                                       gamma0MotherLabel=gammaMC0->GetMother();
+                                               gamma0MotherLabel=gammaMC0->GetMother();
+                                               motherRealLabel=gammaMC0->GetMother();
                                        }
                                }
                                if(gammaMC0->GetPdgCode() ==111){ // Dalitz candidate
@@ -1093,6 +1099,8 @@ void AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueNeutralPionCandidat
                }
                                        
                if(isTruePi0 || isTruePi0Dalitz){// True Pion 
+                       Pi0Candidate->SetTrueMesonValue(1);
+                       Pi0Candidate->SetMCLabel(motherRealLabel);
                        fHistoTrueMotherGammaGammaInvMassPt[fiCut]->Fill(Pi0Candidate->M(),Pi0Candidate->Pt());
                }       
        }
@@ -1500,77 +1508,64 @@ void AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::CalculateBackground(){
 //______________________________________________________________________
 void AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueMesonCandidates(AliAODConversionMother *mesoncand, AliAODConversionMother *TrueNeutralPionCandidate, AliAODConversionPhoton *TrueVirtualParticleCandidate){
 
-//     // Process True Mesons
-// 
-//     AliStack *MCStack = fMCEvent->Stack();
-//     
-//     Bool_t isTrueEta = kFALSE;
-//     Bool_t isTrueOmega = kFALSE;
-//     Int_t gammaMCLabel = TrueNeutralPionCandidate->GetMCParticleLabel(MCStack);
-//     Int_t gammaMotherLabel = -1;
-//     
-//     if(gammaMCLabel != -1){ // Gamma is Combinatorial; MC Particles don't belong to the same Mother
-//             // Daughters Gamma 0
-//             TParticle * negativeMC = (TParticle*)TrueNeutralPionCandidate->GetNegativeMCDaughter(MCStack);
-//             TParticle * positiveMC = (TParticle*)TrueNeutralPionCandidate->GetPositiveMCDaughter(MCStack);
-//             TParticle * gammaMC = (TParticle*)MCStack->Particle(gammaMCLabel);
-// 
-//             if(TMath::Abs(negativeMC->GetPdgCode())==11 && TMath::Abs(positiveMC->GetPdgCode())==11){  // Electrons ...
-//                     if(negativeMC->GetUniqueID() == 5 && positiveMC->GetUniqueID() ==5){ // ... From Conversion ...
-//                             if(gammaMC->GetPdgCode() == 22){ // ... with Gamma Mother
-//                                     gammaMotherLabel=gammaMC->GetFirstMother();
+       // Process True Mesons
+       AliStack *MCStack = fMCEvent->Stack();
+       
+       Bool_t isTrueEta = kFALSE;
+       Bool_t isTrueOmega = kFALSE;
+       Int_t trueMesonFlag = TrueNeutralPionCandidate->GetTrueMesonValue();
+       Int_t pi0MCLabel= TrueNeutralPionCandidate->GetMCLabel();
+
+       
+       if ( !(trueMesonFlag == 1 && pi0MCLabel != -1)) return;
+//     cout << trueMesonFlag << "\t" << pi0MCLabel << endl;
+
+ 
+       Int_t virtualParticleMCLabel = TrueVirtualParticleCandidate->GetMCParticleLabel(MCStack);
+       Int_t virtualParticleMotherLabel = -1;
+
+       Bool_t isPiPiDecay = kFALSE;
+       
+       if (fDoMesonQA){
+               TParticle * negativeMC = (TParticle*)TrueVirtualParticleCandidate->GetNegativeMCDaughter(MCStack);
+               TParticle * positiveMC = (TParticle*)TrueVirtualParticleCandidate->GetPositiveMCDaughter(MCStack);
+               if(TMath::Abs(negativeMC->GetPdgCode())==211 && TMath::Abs(positiveMC->GetPdgCode())==211){  // Pions ...
+                       fHistoTruePionPionInvMassPt[fiCut]->Fill(TrueVirtualParticleCandidate->GetMass(),TrueVirtualParticleCandidate->Pt());
+               }
+       }
+       
+       if(virtualParticleMCLabel != -1){ // if virtualParticleMCLabel==-1 particles don't have same mother 
+               TParticle * negativeMC = (TParticle*)TrueVirtualParticleCandidate->GetNegativeMCDaughter(MCStack);
+               TParticle * positiveMC = (TParticle*)TrueVirtualParticleCandidate->GetPositiveMCDaughter(MCStack);
+//                     TParticle * virtualParticleMotherMC = (TParticle*)MCStack->Particle(virtualParticleMCLabel);
+//                     cout << "pdg code same mother - " << virtualParticleMotherMC->GetPdgCode() << endl;
+               
+               if(TMath::Abs(negativeMC->GetPdgCode())==211 && TMath::Abs(positiveMC->GetPdgCode())==211){  // Pions ...
+                       virtualParticleMotherLabel=virtualParticleMCLabel;
+                       isPiPiDecay=kTRUE;
+//                     } else if(TMath::Abs(negativeMC->GetPdgCode())==11 && TMath::Abs(positiveMC->GetPdgCode())==11){  // Electrons ...
+//                             if( virtualParticleMotherMC->GetPdgCode() != 22 ){
+//                                     virtualParticleMotherLabel=virtualParticleMCLabel;
+//                                     isDalitz = kTRUE;
+//                             } else if(negativeMC->GetUniqueID() == 5 && positiveMC->GetUniqueID() ==5){ // ... From Conversion ...
+//                                     virtualParticleMotherLabel=virtualParticleMotherMC->GetFirstMother();
+//                                     isRealGamma = kTRUE; //no virtual gamma
 //                             }
-//                     }
-//             }
-//             
-//             
-//     }
-// 
-//     Int_t virtualParticleMCLabel = TrueVirtualParticleCandidate->GetMCParticleLabel(MCStack);
-//     Int_t virtualParticleMotherLabel = -1;
-// 
-//     Bool_t isPiPiDecay = kFALSE;
-//     
-//     if (fDoMesonQA){
-//             TParticle * negativeMC = (TParticle*)TrueVirtualParticleCandidate->GetNegativeMCDaughter(MCStack);
-//             TParticle * positiveMC = (TParticle*)TrueVirtualParticleCandidate->GetPositiveMCDaughter(MCStack);
-//             if(TMath::Abs(negativeMC->GetPdgCode())==211 && TMath::Abs(positiveMC->GetPdgCode())==211){  // Pions ...
-//                     fHistoTruePionPionInvMassPt[fiCut]->Fill(TrueVirtualParticleCandidate->GetMass(),TrueVirtualParticleCandidate->Pt());
-//             }
-//     }
-//     
-//     if(virtualParticleMCLabel != -1){ // if virtualParticleMCLabel==-1 particles don't have same mother 
-//             TParticle * negativeMC = (TParticle*)TrueVirtualParticleCandidate->GetNegativeMCDaughter(MCStack);
-//             TParticle * positiveMC = (TParticle*)TrueVirtualParticleCandidate->GetPositiveMCDaughter(MCStack);
-// //                  TParticle * virtualParticleMotherMC = (TParticle*)MCStack->Particle(virtualParticleMCLabel);
-// //                  cout << "pdg code same mother - " << virtualParticleMotherMC->GetPdgCode() << endl;
-//             
-//             if(TMath::Abs(negativeMC->GetPdgCode())==211 && TMath::Abs(positiveMC->GetPdgCode())==211){  // Pions ...
-//                     virtualParticleMotherLabel=virtualParticleMCLabel;
-//                     isPiPiDecay=kTRUE;
-// //                  } else if(TMath::Abs(negativeMC->GetPdgCode())==11 && TMath::Abs(positiveMC->GetPdgCode())==11){  // Electrons ...
-// //                          if( virtualParticleMotherMC->GetPdgCode() != 22 ){
-// //                                  virtualParticleMotherLabel=virtualParticleMCLabel;
-// //                                  isDalitz = kTRUE;
-// //                          } else if(negativeMC->GetUniqueID() == 5 && positiveMC->GetUniqueID() ==5){ // ... From Conversion ...
-// //                                  virtualParticleMotherLabel=virtualParticleMotherMC->GetFirstMother();
-// //                                  isRealGamma = kTRUE; //no virtual gamma
-// //                          }
-//             }       
-//     }
-// 
-//     if(gammaMotherLabel >= 0 && ( gammaMotherLabel == virtualParticleMotherLabel) ){
-//             if(((TParticle*)MCStack->Particle(virtualParticleMotherLabel))->GetPdgCode() == 221){
-//                     isTrueEta=kTRUE;
-//             }
-//             if(((TParticle*)MCStack->Particle(virtualParticleMotherLabel))->GetPdgCode() == 223){
-//                     isTrueOmega=kTRUE;
-//             }
-//     }
-// 
-//     if( isTrueEta || isTrueOmega ){ // True Eta or Omega
-//             if ( isPiPiDecay) { //real eta -> Pi+ Pi- Pi0
-//                     Float_t weighted= 1;
+               }       
+       }
+ 
+       if( pi0MCLabel == virtualParticleMotherLabel ){
+               if(((TParticle*)MCStack->Particle(virtualParticleMotherLabel))->GetPdgCode() == 221){
+                       isTrueEta=kTRUE;
+               }
+               if(((TParticle*)MCStack->Particle(virtualParticleMotherLabel))->GetPdgCode() == 223){
+                       isTrueOmega=kTRUE;
+               }
+       }
+
+       if( isTrueEta || isTrueOmega ){ // True Eta or Omega
+               if ( isPiPiDecay) { //real eta -> Pi+ Pi- Pi0
+                       Float_t weighted= 1;
 //                     if( ((AliPrimaryPionCuts*) fPionCutArray->At(fiCut))->DoWeights() ) { 
 //                             if(((AliConvEventCuts*)fEventCutArray->At(fiCut))->IsParticleFromBGEvent(gammaMotherLabel, fMCStack,fInputEvent)){
 //                                     if (((TParticle*)MCStack->Particle(gammaMotherLabel))->Pt()>0.005){
@@ -1578,10 +1573,10 @@ void AliAnalysisTaskNeutralMesonToPiPlPiMiPiZero::ProcessTrueMesonCandidates(Ali
 //                                     }
 //                             }
 //                     }
-//                     fHistoTruePionPionFromNeutralMesonInvMassPt[fiCut]->Fill(TrueVirtualParticleCandidate->GetMass(),TrueVirtualParticleCandidate->Pt());
-//                     fHistoTrueMotherPiPlPiMiPiZeroInvMassPt[fiCut]->Fill(EtaCandidate->M(),EtaCandidate->Pt(),weighted);
-//             }       
-//     }
+                       fHistoTruePionPionFromNeutralMesonInvMassPt[fiCut]->Fill(TrueVirtualParticleCandidate->GetMass(),TrueVirtualParticleCandidate->Pt());
+                       fHistoTrueMotherPiPlPiMiPiZeroInvMassPt[fiCut]->Fill(mesoncand->M(),mesoncand->Pt(),weighted);
+               }       
+       }
 
 }
 

diff --git a/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pPb.C b/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pPb.C
index 47b7ed10f4f030af413888155f174e6423eb462f..efeaa40309523964fd19087b9d205c5f7fbd344e 100644 (file)
--- a/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pPb.C
+++ b/PWGGA/GammaConv/macros/AddTask_GammaConvCalo_pPb.C
@@ -111,7 +111,7 @@ void AddTask_GammaConvCalo_pPb(  Int_t trainConfig = 1,  //change different set
        task->SetIsMC(isMC);
        // Cut Numbers to use in Analysis
        Int_t numberOfCuts = 2;
-       if (trainConfig==3){ numberOfCuts =6;}
+       if (trainConfig==3 || trainConfig==5  ){ numberOfCuts =6;}
 
        TString *eventCutArray = new TString[numberOfCuts];
        TString *photonCutArray = new TString[numberOfCuts];
@@ -135,6 +135,16 @@ void AddTask_GammaConvCalo_pPb(  Int_t trainConfig = 1,  //change different set
                eventCutArray[ 3] = "8009101"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000040022030000"; mesonCutArray[3] = "01525065000000"; //standard cut, kEMCEJE based on INT7 // EMCAL clusters
                eventCutArray[ 4] = "8009301"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000040022030000"; mesonCutArray[4] = "01525065000000"; //standard cut, kEMCEJ1 based on INT7 // EMCAL clusters
                eventCutArray[ 5] = "8009501"; photonCutArray[ 5] = "002092970028250400000"; clusterCutArray[5] = "10000040022030000"; mesonCutArray[5] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 // EMCAL clusters
+       } else if (trainConfig == 4){ 
+               eventCutArray[ 0] = "8000001"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000040052030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kINT7 // EMCAL clusters
+               eventCutArray[ 1] = "8005201"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000040052030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMC7 // EMCAL clusters     
+       } else if (trainConfig == 5){ 
+               eventCutArray[ 0] = "8008101"; photonCutArray[ 0] = "002092970028250400000"; clusterCutArray[0] = "10000040052030000"; mesonCutArray[0] = "01525065000000"; //standart cut, kEMCEGA based on INT7 // EMCAL clusters
+               eventCutArray[ 1] = "8008301"; photonCutArray[ 1] = "002092970028250400000"; clusterCutArray[1] = "10000040052030000"; mesonCutArray[1] = "01525065000000"; //standard cut, kEMCEG1 based on INT7 // EMCAL clusters
+               eventCutArray[ 2] = "8008501"; photonCutArray[ 2] = "002092970028250400000"; clusterCutArray[2] = "10000040052030000"; mesonCutArray[2] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 // EMCAL clusters
+               eventCutArray[ 3] = "8009101"; photonCutArray[ 3] = "002092970028250400000"; clusterCutArray[3] = "10000040052030000"; mesonCutArray[3] = "01525065000000"; //standard cut, kEMCEJE based on INT7 // EMCAL clusters
+               eventCutArray[ 4] = "8009301"; photonCutArray[ 4] = "002092970028250400000"; clusterCutArray[4] = "10000040052030000"; mesonCutArray[4] = "01525065000000"; //standard cut, kEMCEJ1 based on INT7 // EMCAL clusters
+               eventCutArray[ 5] = "8009501"; photonCutArray[ 5] = "002092970028250400000"; clusterCutArray[5] = "10000040052030000"; mesonCutArray[5] = "01525065000000"; //standard cut, kEMCEG2 based on INT7 // EMCAL clusters             
        } else {
                Error(Form("GammaConvCalo_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
                return;

diff --git a/PWGGA/GammaConv/macros/AddTask_GammaConvV1_pPb.C b/PWGGA/GammaConv/macros/AddTask_GammaConvV1_pPb.C
index e39e326884f54f11bd386c046f8ba5bb1e343fe3..ee540eba23ed532a35cc83363fb0372be6582913 100644 (file)
--- a/PWGGA/GammaConv/macros/AddTask_GammaConvV1_pPb.C
+++ b/PWGGA/GammaConv/macros/AddTask_GammaConvV1_pPb.C
@@ -1019,6 +1019,11 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
                eventCutArray[ 1] = "8000012"; photonCutArray[ 1] = "005092170008260400000"; mesonCutArray[ 1] = "01621035009000"; //new standard RCut=10cm
                eventCutArray[ 2] = "8000012"; photonCutArray[ 2] = "008092170008260400000"; mesonCutArray[ 2] = "01621035009000"; //new standard RCut=12.5cm
                eventCutArray[ 3] = "8000012"; photonCutArray[ 3] = "006092170008260400000"; mesonCutArray[ 3] = "01621035009000"; //new standard RCut=20cm
+       } else if (trainConfig == 181) {
+               eventCutArray[ 0] = "8000011"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; 
+               eventCutArray[ 1] = "8000011"; photonCutArray[ 1] = "002092970028260400000"; mesonCutArray[ 1] = "01525065000000"; 
+               eventCutArray[ 2] = "8000011"; photonCutArray[ 2] = "002092970028270400000"; mesonCutArray[ 2] = "01525065000000"; 
+               eventCutArray[ 3] = "8000011"; photonCutArray[ 3] = "002092970028250000000"; mesonCutArray[ 3] = "01525065000000";                       
        } else {
                Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
                return;
@@ -1060,7 +1065,7 @@ void AddTask_GammaConvV1_pPb(  Int_t trainConfig = 1,  //change different set of
        for(Int_t i = 0; i<numberOfCuts; i++){
                
                analysisEventCuts[i] = new AliConvEventCuts();
-               if ( trainConfig == 1 || trainConfig == 3 || trainConfig == 5 || trainConfig == 7 || trainConfig == 9 || trainConfig == 11 || trainConfig == 13 || trainConfig == 15|| trainConfig == 17|| trainConfig == 19 || trainConfig == 21 || trainConfig == 133 || trainConfig == 135 || trainConfig == 137 || trainConfig == 139 || trainConfig == 141 || trainConfig == 143 || trainConfig == 145 || trainConfig == 147 || trainConfig == 149 || trainConfig == 151 || trainConfig == 173 || trainConfig == 175 || trainConfig == 177 || trainConfig == 179 ){
+               if ( trainConfig == 1 || trainConfig == 3 || trainConfig == 5 || trainConfig == 7 || trainConfig == 9 || trainConfig == 11 || trainConfig == 13 || trainConfig == 15|| trainConfig == 17|| trainConfig == 19 || trainConfig == 21 || trainConfig == 133 || trainConfig == 135 || trainConfig == 137 || trainConfig == 139 || trainConfig == 141 || trainConfig == 143 || trainConfig == 145 || trainConfig == 147 || trainConfig == 149 || trainConfig == 151 || trainConfig == 173 || trainConfig == 175 || trainConfig == 177 || trainConfig == 179 || trainConfig == 181  ){
                        if (doWeighting){
                                if (generatorName.CompareTo("DPMJET")==0){
                                        analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_DPMJET_LHC13b2_efix_pPb_5023GeV_MBV0A", "Eta_DPMJET_LHC13b2_efix_pPb_5023GeV_MBV0A", "","Pi0_Fit_Data_pPb_5023GeV_MBV0A","Eta_Fit_Data_pPb_5023GeV_MBV0A");