Z production and possibility for cuts on secondaries added for

kPyMb and kPyMbNonDiffr. (Zaida Conesa)

diff --git a/PYTHIA6/AliGenPythia.cxx b/PYTHIA6/AliGenPythia.cxx
index c6405e2b906a188d0d1308edb2f3acb903c2c041..f1dc4b74474d3997b8634eb86d9090fedbdb9f90 100644 (file)
--- a/PYTHIA6/AliGenPythia.cxx
+++ b/PYTHIA6/AliGenPythia.cxx
@@ -329,6 +329,7 @@ void AliGenPythia::Init()
     case kPyDirectGamma:
        break;
     case kPyW:
+    case kPyZ:
         break;
     }
 //
@@ -466,7 +467,7 @@ void AliGenPythia::Generate()
        if (fProcess != kPyMb && fProcess != kPyJets && 
            fProcess != kPyDirectGamma &&
            fProcess != kPyMbNonDiffr  &&
-           fProcess != kPyW) {
+           fProcess != kPyW && fProcess != kPyZ ) {
            
            for (i = 0; i < np; i++) {
                TParticle* iparticle = (TParticle *) fParticles->At(i);
@@ -701,9 +702,13 @@ Int_t  AliGenPythia::GenerateMB()
     }
     
 
-    //Introducing child cuts in case kPyW
-    if ( (fProcess == kPyW)  && (fCutOnChild == 1) ) {
-       if ( !CheckKinematicsOnChild() ) return 0;
+    //Introducing child cuts in case kPyW, kPyZ, kPyMb, and kPyMbNonDiff
+    if ( (fProcess == kPyW || fProcess == kPyZ || fProcess == kPyMb || fProcess == kPyMbNonDiffr)  
+        && (fCutOnChild == 1) ) {
+      if ( !CheckKinematicsOnChild() ) {
+       if (pParent)   delete[] pParent;
+       return 0;
+      }
     }
   
 
@@ -958,7 +963,7 @@ Bool_t AliGenPythia::CheckKinematicsOnChild(){
     if (fNumberOfAcceptedParticles != 0) { numberOfAcceptedParticles = fNumberOfAcceptedParticles; }
     Int_t npart = fParticles->GetEntriesFast();
     
-    for (j = 0; j<npart; j++) {  
+    for (j = 0; j<npart; j++) {
        TParticle *  jparticle = (TParticle *) fParticles->At(j);
        kcode = TMath::Abs( CheckPDGCode(jparticle->GetPdgCode()) );
        ks = jparticle->GetStatusCode();
@@ -967,14 +972,13 @@ Bool_t AliGenPythia::CheckKinematicsOnChild(){
        if( (ks == 1)  &&  (kcode == fPdgCodeParticleforAcceptanceCut)  &&  (KinematicSelection(jparticle,1)) ){
            nPartAcc++;
        }
+       if( numberOfAcceptedParticles <= nPartAcc){
+         checking = kTRUE;
+         break;
+       }
     }
-  
-    if( numberOfAcceptedParticles <= nPartAcc){
-       checking = kTRUE;
-    }
-    
+
     return checking;
-    
 }
 
          

diff --git a/PYTHIA6/AliGenPythia.h b/PYTHIA6/AliGenPythia.h
index 92e21b02152477f7b13063eb6af68ab3d28ed93a..3a003ca2c1217dbc506debb9fe33fcc624e39a4b 100644 (file)
--- a/PYTHIA6/AliGenPythia.h
+++ b/PYTHIA6/AliGenPythia.h
@@ -126,7 +126,7 @@ class AliGenPythia : public AliGenMC
     virtual void FinishRun();
     Bool_t CheckTrigger(TParticle* jet1, TParticle* jet2);
 
-    //Used in case kPyW when limiting selected child properties
+    //Used in some processes to selected child properties
     Bool_t CheckKinematicsOnChild();
     
 

diff --git a/PYTHIA6/AliPythia.cxx b/PYTHIA6/AliPythia.cxx
index f5c45413dda07416cc47fb28b70ede96e0c46612..73cf02d493454ae35070fd88293895b4ba039d7b 100644 (file)
--- a/PYTHIA6/AliPythia.cxx
+++ b/PYTHIA6/AliPythia.cxx
@@ -643,6 +643,33 @@ void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfun
       SetMSTP(71,1); //Final QCD & QED showers on
       
       break;  
+
+    case kPyZ:
+
+      //Inclusive production of Z
+      SetMSEL(0);
+      //f fbar -> Z/gamma
+      SetMSUB(1,1);
+      
+      //       // f fbar -> g Z/gamma
+      //       SetMSUB(15,1);
+      //       // f fbar -> gamma Z/gamma
+      //       SetMSUB(19,1);
+      //       // f g -> f Z/gamma
+      //       SetMSUB(30,1);
+      //       // f gamma -> f Z/gamma
+      //       SetMSUB(35,1);
+      
+      //only Z included, not gamma
+      SetMSTP(43,2);
+      
+      // Initial/final parton shower on (Pythia default)
+      // With parton showers on we are generating "Z inclusive process"
+      SetMSTP(61,1); //Initial QCD & QED showers on
+      SetMSTP(71,1); //Final QCD & QED showers on
+      
+      break;  
+
     }
 //
 //  Initialize PYTHIA

diff --git a/PYTHIA6/AliPythia.h b/PYTHIA6/AliPythia.h
index 0f7c2c5d37533a6ad6af36d7ecadfa4941cb70cd..5bf1944587dd2aaf8cc0c640bcc1fcdbc24927bc 100644 (file)
--- a/PYTHIA6/AliPythia.h
+++ b/PYTHIA6/AliPythia.h
@@ -13,7 +13,7 @@ typedef enum
  kPyJpsi, kPyJpsiChi, kPyMb, kPyMbNonDiffr, kPyJets, kPyDirectGamma,
  kPyCharmPbPbMNR, kPyD0PbPbMNR, kPyDPlusPbPbMNR, kPyBeautyPbPbMNR,
  kPyCharmpPbMNR, kPyD0pPbMNR, kPyDPluspPbMNR, kPyBeautypPbMNR,
- kPyCharmppMNR, kPyD0ppMNR, kPyDPlusppMNR, kPyBeautyppMNR, kPyW,
+ kPyCharmppMNR, kPyD0ppMNR, kPyDPlusppMNR, kPyBeautyppMNR, kPyW, kPyZ,
  kPyOldUEQ2ordered, kPyOldUEQ2ordered2, kPyOldPopcorn}
 Process_t;
 /*