X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=blobdiff_plain;f=EVGEN%2FAliGenParam.h;h=a9f2ebfc896c99af2ac5fb309149ad2aaa88e34f;hp=5f2fe0a23ee0e22b77b63dae2f4727911fa0bf66;hb=f2db2f8e718f6a8dc32978165907cfa46eb6a0a1;hpb=fe4da5cc22f890b04843f1aebec0f1bf4f9c3fc9

diff --git a/EVGEN/AliGenParam.h b/EVGEN/AliGenParam.h
index 5f2fe0a23ee..a9f2ebfc896 100644
--- a/EVGEN/AliGenParam.h
+++ b/EVGEN/AliGenParam.h
@@ -1,55 +1,72 @@
+#ifndef ALIGENPARAM_H
+#define ALIGENPARAM_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
 
-#ifndef AliGenParam_H
-#define AliGenParam_H
-#include "AliGenerator.h"
-#include "AliPythia.h"
-#include "TNamed.h"
-#include "TF1.h"
-#include "TArrayF.h"
-#include "TArrayI.h"
-#include "TTree.h"
-#include "TMCParticle.h"
+/* $Id$ */
 
-//-------------------------------------------------------------
-// Generators specific to MUON Arm
+//
+// Class to generate particles from using parametrized pT and y distributions.
+// Distributions are obtained from pointer to object of type AliGenLib.
+// (For example AliGenMUONlib)
+//
+// andreas.morsch@cern.ch
+//
+
+#include "AliGenMC.h"
 
-// Generate heavy mesons - J/Psi, Upsilon, Phi
+class AliPythia;
+class TParticle;
+class AliGenLib;
+class TF1;
 
-class AliGenParam : public AliGenerator
+typedef enum { kAnalog, kNonAnalog} Weighting_t;
+//-------------------------------------------------------------
+class AliGenParam : public AliGenMC
 {
-protected:
-    Double_t (*fPtParaFunc)(Double_t*, Double_t*);
-    Double_t (*fYParaFunc )(Double_t*, Double_t*);
-    Int_t    (*fIpParaFunc )();    
-    TF1* fPtPara;
-    TF1* fYPara;
-    Int_t fIpart;
-    Float_t fdNdy0;
-    Float_t fYWgt;
-    Float_t fPtWgt;
-    Float_t fBias;
-    Int_t   fTrials;
-    Decay_t fForceDecay;
-    TArrayI   fChildSelect;
-    AliPythia *fPythia;
- private:
-    // check if particle is selected as child
-    Bool_t ChildSelected(Int_t ip);
-    // all kinematic selection goes here
-    Bool_t KinematicSelection(TMCParticle *particle);
  public:
-  AliGenParam();
-  AliGenParam(Int_t npart, Int_t ipart);
-//		   Double_t (*PtPara)(Double_t*, Double_t*),
-//		   Double_t (*YPara )(Double_t*, Double_t*));
-  virtual ~AliGenParam();
-  virtual void Generate();
-  virtual void Init();
-  // select particle type
-  virtual void SetPart(Int_t part=443) {fIpart=part;}
-  // force decay type
-  virtual void ForceDecay(Decay_t decay=dimuon) {fForceDecay=decay;}
-  ClassDef(AliGenParam,1)
+    AliGenParam();
+    AliGenParam(Int_t npart, AliGenLib * Library, Int_t param,   const char*  tname = 0);
+    AliGenParam(Int_t npart, Int_t param, const char* tname = 0, const char*  name  = 0);
+    AliGenParam(Int_t npart, Int_t param,
+		Double_t (*PtPara)(const Double_t*, const Double_t*),
+		Double_t (*YPara )(const Double_t*, const Double_t*),
+		Int_t    (*IpPara)(TRandom*)           );
+     
+    virtual ~AliGenParam();
+    virtual void Generate();
+    virtual void Init();
+    // select particle type
+    virtual void SetParam(Int_t param) {fParam = param;}
+    // force decay type
+    virtual void SetWeighting(Weighting_t flag = kAnalog) {fAnalog = flag;}	
+    virtual void SetDeltaPt(Float_t delta=0.01) {fDeltaPt = delta;}
+    virtual void SetDecayer(AliDecayer* decayer) {fDecayer = decayer;}
+    virtual void Draw(const char * opt);
+    TF1 *  GetPt() { return fPtPara;}
+    TF1 *  GetY() {return fYPara;}
+    Float_t GetRelativeArea(Float_t ptMin, Float_t ptMax, Float_t yMin, Float_t yMax, Float_t phiMin, Float_t phiMax);
+
+ protected:
+    Double_t (*fPtParaFunc)(const Double_t*, const Double_t*); //! Pointer to Pt parametrisation function
+    Double_t (*fYParaFunc )(const Double_t*, const Double_t*); //! Pointer to Y parametrisation function
+    Int_t    (*fIpParaFunc )(TRandom*);    //! Pointer to particle type parametrisation function
+    TF1* fPtPara;              // Transverse momentum parameterisation
+    TF1* fYPara;               // Rapidity parameterisation
+    Int_t       fParam;        // Parameterisation type 
+    Float_t     fdNdy0;        // central multiplicity per event
+    Float_t     fYWgt;         // Y-weight
+    Float_t     fPtWgt;        // Pt-weight
+    Float_t     fBias;         // Biasing factor
+    Int_t       fTrials;       // Number of trials
+    Float_t     fDeltaPt;      // pT sampling in steps of fDeltaPt
+    AliDecayer  *fDecayer;     // ! Pointer to pythia object for decays
+
+ private:
+    AliGenParam(const AliGenParam &Param);
+    AliGenParam & operator=(const AliGenParam & rhs);
+
+    ClassDef(AliGenParam,1) // Generator using parameterised pt- and y-distribution
 };
 #endif