- Make sure that all cuts are taken into account independent of process flags.
[u/mrichter/AliRoot.git] / TDPMjet / AliGenDPMjet.cxx
index d20bfcfa8ee467f43835f149feb0eb09fe1bd548..5902af61040985c5c867ea3df5f15ec589b1684d 100644 (file)
@@ -83,6 +83,10 @@ AliGenDPMjet::AliGenDPMjet(Int_t npart)
     AliDpmJetRndm::SetDpmJetRandom(GetRandom());
 }
 
     AliDpmJetRndm::SetDpmJetRandom(GetRandom());
 }
 
+AliGenDPMjet::AliGenDPMjet(const AliGenDPMjet &/*Dpmjet*/)
+:AliGenMC()
+{
+}
 
 //______________________________________________________________________________
 AliGenDPMjet::~AliGenDPMjet()
 
 //______________________________________________________________________________
 AliGenDPMjet::~AliGenDPMjet()
@@ -100,7 +104,7 @@ void AliGenDPMjet::Init()
     SetMC(new TDPMjet(fAProjectile, fZProjectile, fATarget, fZTarget, 
                      fBeamEn,fEnergyCMS));
 
     SetMC(new TDPMjet(fAProjectile, fZProjectile, fATarget, fZTarget, 
                      fBeamEn,fEnergyCMS));
 
-    fDPMjet=(TDPMjet*) fgMCEvGen;
+    fDPMjet=(TDPMjet*) fMCEvGen;
     //
     // **** Flag to force central production
     // fICentr=1. central production forced 
     //
     // **** Flag to force central production
     // fICentr=1. central production forced 
@@ -108,6 +112,8 @@ void AliGenDPMjet::Init()
     // fICentr<-99 -> fraction of x-sec. = XSFRAC                
     // fICentr=-1. -> evaporation/fzc suppressed                 
     // fICentr<-1. -> evaporation/fzc suppressed                 
     // fICentr<-99 -> fraction of x-sec. = XSFRAC                
     // fICentr=-1. -> evaporation/fzc suppressed                 
     // fICentr<-1. -> evaporation/fzc suppressed                 
+    if (fAProjectile == 1 && fZProjectile == 1) fDPMjet->SetfIdp(1);
+    
     fDPMjet->SetfFCentr(fICentr);  
     fDPMjet->SetbRange(fMinImpactParam,fMaxImpactParam); 
     
     fDPMjet->SetfFCentr(fICentr);  
     fDPMjet->SetbRange(fMinImpactParam,fMaxImpactParam); 
     
@@ -267,7 +273,7 @@ void AliGenDPMjet::Generate()
       for (i = 0; i<np; i++) {
          TParticle *  iparticle = (TParticle *) fParticles->At(i);
          Bool_t  hasMother   = (iparticle->GetFirstMother()>=0);
       for (i = 0; i<np; i++) {
          TParticle *  iparticle = (TParticle *) fParticles->At(i);
          Bool_t  hasMother   = (iparticle->GetFirstMother()>=0);
-         Bool_t  hasDaughter = (iparticle->GetFirstDaughter()>=0);
+//       Bool_t  hasDaughter = (iparticle->GetFirstDaughter()>=0);
 
          if (pSelected[i]) {
              
 
          if (pSelected[i]) {
              
@@ -544,7 +550,7 @@ void AliGenDPMjet::MakeHeader()
 
 
 //______________________________________________________________________________
 
 
 //______________________________________________________________________________
-AliGenDPMjet& AliGenDPMjet::operator=(const  AliGenDPMjet& rhs)
+AliGenDPMjet& AliGenDPMjet::operator=(const  AliGenDPMjet& /*rhs*/)
 {
 // Assignment operator
     return *this;
 {
 // Assignment operator
     return *this;