Slow nucleon generator and model.

[u/mrichter/AliRoot.git] / EVGEN / AliGenSlowNucleons.h

diff --git a/EVGEN/AliGenSlowNucleons.h b/EVGEN/AliGenSlowNucleons.h
new file mode 100644 (file)
index 0000000..aa48c3e
--- /dev/null
+++ b/EVGEN/AliGenSlowNucleons.h
@@ -0,0 +1,67 @@
+#ifndef ALIGENSLOWNUCLEONS_H
+#define ALIGENSLOWNUCLEONS_H
+/* Copyright(c) 198-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/* $Id$ */
+
+#include "AliGenerator.h"
+class AliSlowNucleonModel;
+class TH2F;
+
+
+class AliGenSlowNucleons : public AliGenerator
+{
+public:
+    AliGenSlowNucleons();
+    AliGenSlowNucleons(Int_t npart);
+    virtual ~AliGenSlowNucleons();
+    virtual void Init();
+    virtual void FinishRun();
+    virtual void Generate();
+    virtual void SetPmax(Float_t pmax = 10.) {fPmax = pmax;}
+    virtual void SetNominalCmsEnergy(Float_t energy = 14000.) {fCMS = energy;}
+    virtual void SetTarget(Float_t a=208, Float_t z=82) {fATarget = a; fZTarget = z;}
+    virtual void SetCharge(Int_t c = 1) {fCharge = c;}
+    virtual void SetTemperature(Double_t t1 = 0.04, Double_t t2 = 0.004)
+       {fTemperatureG = t1; fTemperatureB = t2;}
+    virtual void SetBetaSource(Double_t b1 = 0.05, Double_t b2 = 0.)
+       {fBetaSourceG = b1; fBetaSourceB = b2;}
+    //
+    virtual void SetSlowNucleonModel(AliSlowNucleonModel* model) 
+       {fSlowNucleonModel = model;}
+    virtual Bool_t NeedsCollisionGeometry() const {return kTRUE;}
+    virtual void   SetCollisionGeometry(AliCollisionGeometry* geom)
+       {fCollisionGeometry = geom;}
+    virtual void   SetDebug(Int_t flag = 0) {fDebug = flag;}
+           
+ protected:
+    void     GenerateSlow(Int_t charge, Double_t T, Double_t beta, Float_t* q);
+    Double_t Maxwell(Double_t m, Double_t p, Double_t t);
+    void     Lorentz(Double_t m, Double_t beta, Float_t* q);
+ protected:
+    Float_t  fCMS;          // Center of mass energy
+    Float_t  fMomentum;     // Target nucleus momentum
+    Float_t  fBeta;         // Target nucleus beta
+    Float_t  fPmax;         // Maximum slow nucleon momentum
+    Float_t  fATarget;      // Target nucleus mass number
+    Float_t  fZTarget;      // Target nucleus charge number
+    Int_t    fCharge;       // Slow nucleon charge
+    Float_t  fTemperatureG; // Source Temperature for gray nucleons
+    Float_t  fBetaSourceG;  // Source beta for gray nucleons
+    Float_t  fTemperatureB; // Source Temperature for black nucleons
+    Float_t  fBetaSourceB;  // Source beta for black nucleons
+    Int_t    fDebug;        // Debug flag
+    TH2F*    fDebugHist;    // Histogram for debugging
+    
+    //
+    AliSlowNucleonModel* fSlowNucleonModel; // The slow nucleon model
+    ClassDef(AliGenSlowNucleons,1) // Slow Nucleon Generator
+};
+#endif
+
+
+
+
+
+