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

#ifndef ALIGENCORRHF_H
#define ALIGENCORRHF_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id$ */

//
// Class to generate correlated Heavy Flavor hadron pairs using paramtrized
// kinematics of quark pairs from some generator and quark fragmentation functions.
// Is a generalisation of AliGenParam class for correlated pairs of hadrons.
// Author: S. Grigoryan, LPC Clermont-Fd & YerPhI, Smbat.Grigoryan@cern.ch
//

#include "AliGenMC.h"

class AliPythia;
class TParticle;
class TRandom;
class TFile;
class TString;

//-------------------------------------------------------------
class AliGenCorrHF : public AliGenMC
{
 public:
    AliGenCorrHF();
    AliGenCorrHF(Int_t npart, Int_t param);
    AliGenCorrHF(char* tname, Int_t npart, Int_t param);
     
    virtual ~AliGenCorrHF();
    virtual void Generate();
    virtual void Init();
    // force decay type
    virtual void SetDecayer(AliDecayer* decayer) {fDecayer = decayer;}

    // Particle type parametrisation functions, needed by GetHadronPair
    static Int_t IpCharm(TRandom* ran);
    static Int_t IpBeauty(TRandom* ran);

    // Computation of cumulative sums of cell wght-s, needed by GetQuarkPair
    static Double_t ComputeIntegral(TFile* fG);

    // fG - input file w QQbar kinematical grid (TTree) and fragm. functions (24 TH2-s)
    static void GetQuarkPair(TFile* fG, Double_t* fInt, Double_t &y1, Double_t &y2, Double_t &pt1, Double_t &pt2, Double_t &dphi);              
    static void GetHadronPair(TFile* fG, Int_t idq, Double_t y1, Double_t y2, Double_t pt1, Double_t pt2, Int_t &id3, Int_t &id4, Double_t &pz3, Double_t &pz4, Double_t &pt3, Double_t &pt4); 

 protected:
    TString     fFileName;    // Name of the input file
    TFile*      fFile;        //! Pointer to input file
    Int_t       fQuark;       // Quark type 
    Float_t     fBias;        // Biasing factor
    Int_t       fTrials;      // Number of trials
    AliDecayer* fDecayer;     //! Pointer to pythia object for decays

 private:
    AliGenCorrHF(const AliGenCorrHF &CorrHF);
    AliGenCorrHF & operator=(const AliGenCorrHF & rhs);

    static Double_t* fgIntegral; //! Pointer to array of cumulative sums of wght-s
    static Int_t  fgnptbins;             // =12 Number of bins for the fragm. 
                                         //   function dependence on quark pt
    // Number of the grid bins in deltaphi, y and pt:  18, 30 and 32
    static Double_t fgdph[19];           // deltaphi bin coordinates  
    static Double_t fgy[31];             // y bin coordinates
    static Double_t fgpt[33];            // pt bin coordinates
    static Double_t fgptbmin[12];        // min & max coordinates of pt bins for
    static Double_t fgptbmax[12];        // the fragm. function
    
    ClassDef(AliGenCorrHF,1)  // Generator using parameterized QQbar & fragm. functions
};
#endif
Commit	Line	Data
2c890605	1	#ifndef ALIGENCORRHF_H
	2	#define ALIGENCORRHF_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/* $Id$ */
	7
	8	//
	9	// Class to generate correlated Heavy Flavor hadron pairs using paramtrized
	10	// kinematics of quark pairs from some generator and quark fragmentation functions.
	11	// Is a generalisation of AliGenParam class for correlated pairs of hadrons.
	12	// Author: S. Grigoryan, LPC Clermont-Fd & YerPhI, Smbat.Grigoryan@cern.ch
	13	//
	14
	15	#include "AliGenMC.h"
	16
	17	class AliPythia;
	18	class TParticle;
	19	class TRandom;
	20	class TFile;
	21	class TString;
	22
	23	//-------------------------------------------------------------
	24	class AliGenCorrHF : public AliGenMC
	25	{
	26	public:
	27	AliGenCorrHF();
	28	AliGenCorrHF(Int_t npart, Int_t param);
	29	AliGenCorrHF(char* tname, Int_t npart, Int_t param);
2c890605	30
	31	virtual ~AliGenCorrHF();
	32	virtual void Generate();
	33	virtual void Init();
	34	// force decay type
	35	virtual void SetDecayer(AliDecayer* decayer) {fDecayer = decayer;}
	36
2c890605	37	// Particle type parametrisation functions, needed by GetHadronPair
	38	static Int_t IpCharm(TRandom* ran);
	39	static Int_t IpBeauty(TRandom* ran);
	40
	41	// Computation of cumulative sums of cell wght-s, needed by GetQuarkPair
	42	static Double_t ComputeIntegral(TFile* fG);
	43
	44	// fG - input file w QQbar kinematical grid (TTree) and fragm. functions (24 TH2-s)
	45	static void GetQuarkPair(TFile* fG, Double_t* fInt, Double_t &y1, Double_t &y2, Double_t &pt1, Double_t &pt2, Double_t &dphi);
	46	static void GetHadronPair(TFile* fG, Int_t idq, Double_t y1, Double_t y2, Double_t pt1, Double_t pt2, Int_t &id3, Int_t &id4, Double_t &pz3, Double_t &pz4, Double_t &pt3, Double_t &pt4);
	47
	48	protected:
	49	TString fFileName; // Name of the input file
	50	TFile* fFile; //! Pointer to input file
	51	Int_t fQuark; // Quark type
	52	Float_t fBias; // Biasing factor
	53	Int_t fTrials; // Number of trials
	54	AliDecayer* fDecayer; //! Pointer to pythia object for decays
	55
	56	private:
93a2041b	57	AliGenCorrHF(const AliGenCorrHF &CorrHF);
	58	AliGenCorrHF & operator=(const AliGenCorrHF & rhs);
	59
2c890605	60	static Double_t* fgIntegral; //! Pointer to array of cumulative sums of wght-s
	61	static Int_t fgnptbins; // =12 Number of bins for the fragm.
	62	// function dependence on quark pt
	63	// Number of the grid bins in deltaphi, y and pt: 18, 30 and 32
	64	static Double_t fgdph[19]; // deltaphi bin coordinates
	65	static Double_t fgy[31]; // y bin coordinates
	66	static Double_t fgpt[33]; // pt bin coordinates
	67	static Double_t fgptbmin[12]; // min & max coordinates of pt bins for
	68	static Double_t fgptbmax[12]; // the fragm. function
	69
	70	ClassDef(AliGenCorrHF,1) // Generator using parameterized QQbar & fragm. functions
	71	};
	72	#endif
	73
	74
	75
	76
	77
	78
	79
	80
	81
	82