#include <TParticle.h>
#include <TMath.h>
+#include "HBTprocCOMMON.h"
#include <Riostream.h>
#ifndef WIN32
//prints values set in common blocks
ctest();
}
+void THBTprocessor::SetTrackRejectionFactor(Float_t trf) const
+{
+ //setter
+ PARAMETERS.trk_accep = trf;
+}
+void THBTprocessor::SetRefControl(Int_t rc) const
+{
+ //setter
+ PARAMETERS.ref_control = rc;
+}
+void THBTprocessor::SetPIDs(Int_t pid1,Int_t pid2) const
+{
+ //setter
+ PARAMETERS.pid[0]=pid1; PARAMETERS.pid[1]=pid2;
+}
+void THBTprocessor::SetNPIDtypes(Int_t npidt)const
+{
+ //setter
+ PARAMETERS.n_pid_types = npidt;
+}
+void THBTprocessor::SetDeltap(Float_t deltp) const
+{
+ //setter
+ PARAMETERS.deltap=deltp;
+}
+void THBTprocessor::SetMaxIterations(Int_t maxiter) const
+{
+ //setter
+ PARAMETERS.maxit = maxiter;
+}
+void THBTprocessor::SetDelChi(Float_t dc)const
+{
+ //setter
+ PARAMETERS.delchi = dc;
+}
+void THBTprocessor::SetIRand(Int_t irnd) const
+{
+ //setter
+ PARAMETERS.irand = irnd;
+}
+void THBTprocessor::SetLambda(Float_t lam) const
+{
+ //setter
+ PARAMETERS.lambda = lam;
+}
+void THBTprocessor::SetR1d(Float_t r) const
+{
+ //setter
+ PARAMETERS.R_1d=r;
+}
+void THBTprocessor::SetRSide(Float_t rs) const
+{
+ //setter
+ PARAMETERS.Rside=rs;
+}
+void THBTprocessor::SetROut(Float_t ro) const
+{
+ //setter
+ PARAMETERS.Rout=ro;
+}
+void THBTprocessor::SetRLong(Float_t rl) const
+{
+ //setter
+ PARAMETERS.Rlong=rl;
+}
+void THBTprocessor::SetRPerp(Float_t rp) const
+{
+ //setter
+ PARAMETERS.Rperp=rp;
+}
+void THBTprocessor::SetRParallel(Float_t rprl) const
+{
+ //setter
+ PARAMETERS.Rparallel=rprl;
+}
+void THBTprocessor::SetR0(Float_t r0) const
+{
+ //setter
+ PARAMETERS.R0=r0;
+}
+void THBTprocessor::SetQ0(Float_t q0) const
+{
+ //setter
+PARAMETERS.Q0=q0;
+}
+void THBTprocessor::SetSwitch1D(Int_t s1d) const
+{
+ //setter
+ PARAMETERS.switch_1d = s1d;
+}
+void THBTprocessor::SetSwitch3D(Int_t s3d) const
+{
+ //setter
+ PARAMETERS.switch_3d = s3d;
+}
+void THBTprocessor::SetSwitchType(Int_t st) const
+{
+ //setter
+ PARAMETERS.switch_type = st;
+}
+
+void THBTprocessor::SetSwitchCoherence(Int_t sc) const
+{
+ //setter
+ PARAMETERS. switch_coherence = sc;
+}
+void THBTprocessor::SetSwitchCoulomb(Int_t scol) const
+{
+ //setter
+ PARAMETERS. switch_coulomb = scol;
+}
+void THBTprocessor::SetSwitchFermiBose(Int_t sfb) const
+{
+ //setter
+ PARAMETERS.switch_fermi_bose = sfb;
+}
+
+void THBTprocessor::SetPtRange(Float_t ptmin, Float_t ptmax)const //Pt in GeV/c
+{
+ //setter
+ MESH.pt_min=ptmin;MESH.pt_max=ptmax;
+}
+
+void THBTprocessor::SetPxRange(Float_t pxmin, Float_t pxmax) const
+{
+ //setter
+ MESH.px_min=pxmin;MESH.px_max=pxmax;
+}
+
+void THBTprocessor::SetPyRange(Float_t pymin, Float_t pymax)const
+{
+ //setter
+ MESH.py_min=pymin;MESH.py_max=pymax;
+}
+
+void THBTprocessor::SetPzRange(Float_t pzmin, Float_t pzmax) const
+{
+ //setter
+ MESH.pz_min=pzmin;MESH.pz_max=pzmax;
+}
+
+void THBTprocessor::SetPhiRange(Float_t phimin, Float_t phimax)const //Angle in degrees
+{
+ //setter
+ MESH.phi_min=phimin;MESH.phi_max=phimax;
+}
+
+void THBTprocessor::SetEtaRange(Float_t etamin, Float_t etamax)const //Pseudorapidity !!!!!!!!!
+{
+ //setter
+ MESH.eta_min=etamin;MESH.eta_max=etamax;
+}
+
+void THBTprocessor::SetNPtBins(Int_t nptbin)const
+{
+//setter
+MESH.n_pt_bins=nptbin;
+}
+void THBTprocessor::SetNPhiBins(Int_t nphibin)const
+{
+//setter
+MESH.n_phi_bins=nphibin;
+}
+void THBTprocessor::SetNEtaBins(Int_t netabin)const
+{
+//setter
+MESH.n_eta_bins=netabin;
+}
+
+void THBTprocessor::SetNPxBins(Int_t npxbin)const
+{
+//setter
+MESH.n_px_bins=npxbin;
+}
+void THBTprocessor::SetNPyBins(Int_t npybin)const
+{
+//setter
+MESH.n_py_bins=npybin;
+}
+void THBTprocessor::SetNPzBins(Int_t npzbin)const
+{
+//setter
+MESH.n_pz_bins=npzbin;
+}
+
+void THBTprocessor::SetNBins1DFineMesh(Int_t n)const
+{
+//setter
+MESH.n_1d_fine=n;
+}
+void THBTprocessor::SetBinSize1DFineMesh(Float_t x)const
+{
+//setter
+MESH.binsize_1d_fine=x;
+}
+
+void THBTprocessor::SetNBins1DCoarseMesh(Int_t n)const
+{
+//setter
+MESH.n_1d_coarse =n;
+}
+void THBTprocessor::SetBinSize1DCoarseMesh(Float_t x)const
+{
+//setter
+MESH.binsize_1d_coarse=x;
+}
+
+void THBTprocessor::SetNBins3DFineMesh(Int_t n)const
+{
+//setter
+MESH.n_3d_fine =n;
+}
+void THBTprocessor::SetBinSize3DFineMesh(Float_t x)const
+{
+//setter
+MESH.binsize_3d_fine=x;
+}
+
+void THBTprocessor::SetNBins3DCoarseMesh(Int_t n)const
+{
+//setter
+MESH.n_3d_coarse=n;
+}
+void THBTprocessor::SetBinSize3DCoarseMesh(Float_t x)const
+{
+//setter
+MESH.binsize_3d_coarse=x;
+}
+
+void THBTprocessor::SetNBins3DFineProjectMesh(Int_t n )const
+{
+//setter
+MESH.n_3d_fine_project=n;
+}
+
+void THBTprocessor::SetPrintFull(Int_t flag) const
+{
+ //setter
+ PARAMETERS.print_full=flag;
+}
//
#include <TGenerator.h>
-#include "HBTprocCOMMON.h"
-
class THBTprocessor: public TGenerator
{
//comprehensive description off all these methods
//can be found in proper methods of AliGenHBTprocessor
- virtual void SetTrackRejectionFactor(Float_t trf = 1.0) const {PARAMETERS.trk_accep = trf;}
- virtual void SetRefControl(Int_t rc =2) const {PARAMETERS.ref_control = rc;}
- virtual void SetPIDs(Int_t pid1 = 8,Int_t pid2 = 9) const {PARAMETERS.pid[0]=pid1; PARAMETERS.pid[1]=pid2;}
- virtual void SetNPIDtypes(Int_t npidt = 2)const {PARAMETERS.n_pid_types = npidt;}
- virtual void SetDeltap(Float_t deltp = 0.1) const {PARAMETERS.deltap=deltp; }
- virtual void SetMaxIterations(Int_t maxiter = 50) const {PARAMETERS.maxit = maxiter;}
- virtual void SetDelChi(Float_t dc = 0.1)const {PARAMETERS.delchi = dc;}
- virtual void SetIRand(Int_t irnd = 76564) const {PARAMETERS.irand = irnd;}
- virtual void SetLambda(Float_t lam = 0.6) const { PARAMETERS.lambda = lam;}
- virtual void SetR1d(Float_t r = 7.0) const {PARAMETERS.R_1d=r;}
- virtual void SetRSide(Float_t rs = 6.0) const {PARAMETERS.Rside=rs;}
- virtual void SetROut(Float_t ro = 7.0) const {PARAMETERS.Rout=ro;}
- virtual void SetRLong(Float_t rl = 4.0) const {PARAMETERS.Rlong=rl;}
- virtual void SetRPerp(Float_t rp = 6.0) const {PARAMETERS.Rperp=rp;}
- virtual void SetRParallel(Float_t rprl = 4.0) const {PARAMETERS.Rparallel=rprl;}
- virtual void SetR0(Float_t r0 = 4.0) const {PARAMETERS.R0=r0;}
- virtual void SetQ0(Float_t q0 = 9.0) const {PARAMETERS.Q0=q0;}
- virtual void SetSwitch1D(Int_t s1d = 3) const {PARAMETERS.switch_1d = s1d;}
- virtual void SetSwitch3D(Int_t s3d = 0) const {PARAMETERS.switch_3d = s3d;}
- virtual void SetSwitchType(Int_t st = 3) const {PARAMETERS.switch_type = st;}
- virtual void SetSwitchCoherence(Int_t sc = 0) const {PARAMETERS.switch_coherence = sc;}
- virtual void SetSwitchCoulomb(Int_t scol = 2) const {PARAMETERS.switch_coulomb = scol;}
- virtual void SetSwitchFermiBose(Int_t sfb = 1) const {PARAMETERS.switch_fermi_bose = sfb;}
-
+ virtual void SetTrackRejectionFactor(Float_t trf = 1.0) const;
+ virtual void SetRefControl(Int_t rc =2) const;
+ virtual void SetPIDs(Int_t pid1 = 8,Int_t pid2 = 9) const;
+ virtual void SetNPIDtypes(Int_t npidt = 2)const;
+ virtual void SetDeltap(Float_t deltp = 0.1) const;
+ virtual void SetMaxIterations(Int_t maxiter = 50) const;
+ virtual void SetDelChi(Float_t dc = 0.1)const;
+ virtual void SetIRand(Int_t irnd = 76564) const;
+ virtual void SetLambda(Float_t lam = 0.6) const;
+ virtual void SetR1d(Float_t r = 7.0) const;
+ virtual void SetRSide(Float_t rs = 6.0) const;
+ virtual void SetROut(Float_t ro = 7.0) const ;
+ virtual void SetRLong(Float_t rl = 4.0) const;
+ virtual void SetRPerp(Float_t rp = 6.0) const;
+ virtual void SetRParallel(Float_t rprl = 4.0) const;
+ virtual void SetR0(Float_t r0 = 4.0) const;
+ virtual void SetQ0(Float_t q0 = 9.0) const ;
+ virtual void SetSwitch1D(Int_t s1d = 3) const ;
+ virtual void SetSwitch3D(Int_t s3d = 0) const ;
+ virtual void SetSwitchType(Int_t st = 3) const ;
+ virtual void SetSwitchCoherence(Int_t sc = 0) const ;
+ virtual void SetSwitchCoulomb(Int_t scol = 2) const ;
+ virtual void SetSwitchFermiBose(Int_t sfb = 1) const ;
- virtual void SetPtRange(Float_t ptmin = 0.1, Float_t ptmax = 0.98)const //Pt in GeV/c
- { MESH.pt_min=ptmin;MESH.pt_max=ptmax;}
-
- virtual void SetPxRange(Float_t pxmin = -1.0, Float_t pxmax = 1.0) const
- { MESH.px_min=pxmin;MESH.px_max=pxmax;}
-
- virtual void SetPyRange(Float_t pymin = -1.0, Float_t pymax = 1.0)const
- { MESH.py_min=pymin;MESH.py_max=pymax;}
-
- virtual void SetPzRange(Float_t pzmin = -3.6, Float_t pzmax = 3.6) const
- { MESH.pz_min=pzmin;MESH.pz_max=pzmax;}
-
-
- virtual void SetPhiRange(Float_t phimin = 0.0, Float_t phimax = 360.0)const //Angle in degrees
- { MESH.phi_min=phimin;MESH.phi_max=phimax;}
-
- virtual void SetEtaRange(Float_t etamin = -1.5, Float_t etamax = 1.5)const //Pseudorapidity !!!!!!!!!
- { MESH.eta_min=etamin;MESH.eta_max=etamax;}
-
- virtual void SetNPtBins(Int_t nptbin = 50)const{MESH.n_pt_bins=nptbin;}
- virtual void SetNPhiBins(Int_t nphibin = 50)const{MESH.n_phi_bins=nphibin;}
- virtual void SetNEtaBins(Int_t netabin = 50)const{MESH.n_eta_bins=netabin;}
-
- virtual void SetNPxBins(Int_t npxbin = 20)const{MESH.n_px_bins=npxbin;}
- virtual void SetNPyBins(Int_t npybin = 20)const{MESH.n_py_bins=npybin;}
- virtual void SetNPzBins(Int_t npzbin = 70)const{MESH.n_pz_bins=npzbin;}
-
- virtual void SetNBins1DFineMesh(Int_t n = 10)const{ MESH.n_1d_fine=n;}
- virtual void SetBinSize1DFineMesh(Float_t x=0.01)const{MESH.binsize_1d_fine=x;}
-
- virtual void SetNBins1DCoarseMesh(Int_t n =2 )const{MESH.n_1d_coarse =n;}
- virtual void SetBinSize1DCoarseMesh(Float_t x=0.05)const{MESH.binsize_1d_coarse=x;}
-
- virtual void SetNBins3DFineMesh(Int_t n = 8)const{MESH.n_3d_fine =n;}
- virtual void SetBinSize3DFineMesh(Float_t x=0.01)const{MESH.binsize_3d_fine=x;}
-
- virtual void SetNBins3DCoarseMesh(Int_t n = 2)const{ MESH.n_3d_coarse=n;}
- virtual void SetBinSize3DCoarseMesh(Float_t x=0.08)const{MESH.binsize_3d_coarse=x;}
-
- virtual void SetNBins3DFineProjectMesh(Int_t n =3 )const{ MESH.n_3d_fine_project=n;}
-
- virtual void SetPrintFull(Int_t flag = 1)const{PARAMETERS.print_full=flag;}
+ virtual void SetPtRange(Float_t ptmin = 0.1, Float_t ptmax = 0.98)const ; //Pt in GeV/c
+ virtual void SetPxRange(Float_t pxmin = -1.0, Float_t pxmax = 1.0) const;
+ virtual void SetPyRange(Float_t pymin = -1.0, Float_t pymax = 1.0)const ;
+ virtual void SetPzRange(Float_t pzmin = -3.6, Float_t pzmax = 3.6) const;
+ virtual void SetPhiRange(Float_t phimin = 0.0, Float_t phimax = 360.0)const; //Angle in degrees
+ virtual void SetEtaRange(Float_t etamin = -1.5, Float_t etamax = 1.5)const; //Pseudorapidity !!!!!!!!!
+ virtual void SetNPtBins(Int_t nptbin = 50)const;
+ virtual void SetNPhiBins(Int_t nphibin = 50)const;
+ virtual void SetNEtaBins(Int_t netabin = 50)const;
+ virtual void SetNPxBins(Int_t npxbin = 20)const;
+ virtual void SetNPyBins(Int_t npybin = 20)const;
+ virtual void SetNPzBins(Int_t npzbin = 70)const;
+ virtual void SetNBins1DFineMesh(Int_t n = 10)const;
+ virtual void SetBinSize1DFineMesh(Float_t x=0.01)const;
+ virtual void SetNBins1DCoarseMesh(Int_t n =2 )const;
+ virtual void SetBinSize1DCoarseMesh(Float_t x=0.05)const;
+ virtual void SetNBins3DFineMesh(Int_t n = 8)const;
+ virtual void SetBinSize3DFineMesh(Float_t x=0.01)const;
+ virtual void SetNBins3DCoarseMesh(Int_t n = 2)const;
+ virtual void SetBinSize3DCoarseMesh(Float_t x=0.08)const;
+ virtual void SetNBins3DFineProjectMesh(Int_t n =3 )const;
+ virtual void SetPrintFull(Int_t flag = 1)const;
void DumpSettings() const;