Possibility to cut in the total momentum. (B. Vulpescu)
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 5 Oct 2007 20:48:08 +0000 (20:48 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 5 Oct 2007 20:48:08 +0000 (20:48 +0000)
EVGEN/AliGenMUONCocktailpp.cxx
EVGEN/AliGenMUONCocktailpp.h

index 596fcb6..2a8ecde 100644 (file)
@@ -24,6 +24,7 @@
 //      production cross-sections in pp collisions 
 // July 07:added heavy quark production from AliGenCorrHF and heavy quark 
 //         production switched off in Pythia 
+// Aug. 07: added trigger cut on total momentum
 
 #include <TObjArray.h>
 #include <TParticle.h>
@@ -55,8 +56,9 @@ AliGenMUONCocktailpp::AliGenMUONCocktailpp()
      fTotalRate(0),
      fMuonMultiplicity(0),
      fMuonPtCut(0.),
+     fMuonPCut(0.),
      fMuonThetaMinCut(0.),
-     fMuonThetaMaxCut(0.),
+     fMuonThetaMaxCut(180.),
      fMuonOriginCut(-999.),
      fNSucceded(0),
      fNGenerated(0)
@@ -327,7 +329,8 @@ void AliGenMUONCocktailpp::Generate()
            if((part->Vz() > fMuonOriginCut) && //take only the muons that decayed before the abs + 1 int. length in C abs
               (part->Theta()*180./TMath::Pi()>fMuonThetaMinCut) &&
               (part->Theta()*180./TMath::Pi()<fMuonThetaMaxCut) &&
-              (part->Pt()>fMuonPtCut)) {
+              (part->Pt()>fMuonPtCut) &&
+              (part->P()>fMuonPCut)) {
              numberOfMuons++;
            }
          }
index a46290e..524d42b 100644 (file)
@@ -20,6 +20,7 @@ class AliGenMUONCocktailpp : public AliGenCocktail
     Int_t   GetNGenerated()        const {return fNGenerated;}
     Int_t   GetMuonMultiplicity()  const {return fMuonMultiplicity;}
     Float_t GetMuonPtCut()         const {return fMuonPtCut;}
+    Float_t GetMuonPCut()          const {return fMuonPCut;}    
     Float_t GetMuonThetaMin()      const {return fMuonThetaMinCut;}
     Float_t GetMuonThetaMax()      const {return fMuonThetaMaxCut;}        
     Float_t GetMuonOriginCut()     const {return fMuonOriginCut;}          
@@ -28,6 +29,7 @@ class AliGenMUONCocktailpp : public AliGenCocktail
     
     void    SetMuonMultiplicity(Int_t MuonMultiplicity) { fMuonMultiplicity = MuonMultiplicity;}
     void    SetMuonPtCut(Float_t PtCut) { fMuonPtCut = PtCut;}
+    void    SetMuonPCut(Float_t PCut) { fMuonPCut = PCut;}    
     void    SetMuonOriginCut(Float_t originCut) { fMuonOriginCut = originCut;}
     void    SetMuonThetaRange(Float_t ThetaMin, Float_t ThetaMax){
        fMuonThetaMinCut=ThetaMin;
@@ -47,6 +49,7 @@ class AliGenMUONCocktailpp : public AliGenCocktail
     Float_t fTotalRate;// Total rate of the full cocktail processes
     Int_t   fMuonMultiplicity; // Muon multiplicity for the primordial trigger
     Float_t fMuonPtCut;// Transverse momentum cut for muons
+    Float_t fMuonPCut;// Momentum cut for muons    
     Float_t fMuonThetaMinCut;// Minimum theta cut for muons
     Float_t fMuonThetaMaxCut; // Maximum theta cut for muons
     Float_t fMuonOriginCut; //use only muons whose "part->Vz()" value is larger than fMuonOrigin
@@ -54,7 +57,7 @@ class AliGenMUONCocktailpp : public AliGenCocktail
     Int_t   fNGenerated;// Number of generated cocktails
     
 
-    ClassDef(AliGenMUONCocktailpp,1)  //  cocktail for physics in the Alice
+    ClassDef(AliGenMUONCocktailpp,2)  //  cocktail for physics in the Alice
 };
 
 #endif