Add Boost() method to boost all particles to LHC lab frame. Needed for asymmetric...

[u/mrichter/AliRoot.git] / THbtp / THBTprocessor.h

#ifndef THBTPROCESSOR_H
#define THBTPROCESSOR_H
/*******************************************************/
//Author Piotr Krzysztof Skowronski e-mial: Piotr.Skowronski@cern.ch
// C++ Wrapper Class for fortran made HBT Processor written by Lanny Ray
// 
// Comunication is done via COMMON BLOCKS declared in HBTprocCOMMON.h
// using cfortran.h routines
// User should use class AliGenHBTprocessor and all their interface
// see there for more description
//

#include <TGenerator.h>
#include "HBTprocCOMMON.h"


class THBTprocessor: public TGenerator
  {
    public:
      THBTprocessor(); 
      virtual ~THBTprocessor() {};
      
      virtual void  Initialize();
      virtual void  GenerateEvent();
      virtual void  PrintEvent();
      virtual Int_t ImportParticles(TClonesArray *particles, Option_t *option="");

      //Set/Gets
      //comprehensive description off all these methods 
      //can be found in proper methods of AliGenHBTprocessor
      
      virtual void SetTrackRejectionFactor(Float_t trf = 1.0) {PARAMETERS.trk_accep = trf;}
      virtual void SetRefControl(Int_t rc =2) {PARAMETERS.ref_control = rc;}
      virtual void SetPIDs(Int_t pid1 = 8,Int_t pid2 = 9) {PARAMETERS.pid[0]=pid1; PARAMETERS.pid[1]=pid2;}
      virtual void SetNPIDtypes(Int_t npidt = 2){PARAMETERS.n_pid_types = npidt;}
      virtual void SetDeltap(Float_t deltp = 0.1) {PARAMETERS.deltap=deltp; }
      virtual void SetMaxIterations(Int_t maxiter = 50) {PARAMETERS.maxit = maxiter;}
      virtual void SetDelChi(Float_t dc = 0.1){PARAMETERS.delchi = dc;}
      virtual void SetIRand(Int_t irnd = 76564) {PARAMETERS.irand = irnd;}
      virtual void SetLambda(Float_t lam = 0.6) { PARAMETERS.lambda = lam;}
      virtual void SetR1d(Float_t r = 7.0) {PARAMETERS.R_1d=r;}
      virtual void SetRSide(Float_t rs = 6.0) {PARAMETERS.Rside=rs;}
      virtual void SetROut(Float_t ro = 7.0) {PARAMETERS.Rout=ro;}
      virtual void SetRLong(Float_t rl = 4.0) {PARAMETERS.Rlong=rl;}
      virtual void SetRPerp(Float_t rp = 6.0) {PARAMETERS.Rperp=rp;}
      virtual void SetRParallel(Float_t rprl = 4.0) {PARAMETERS.Rparallel=rprl;}
      virtual void SetR0(Float_t r0 = 4.0) {PARAMETERS.R0=r0;}
      virtual void SetQ0(Float_t q0 = 9.0) {PARAMETERS.Q0=q0;}
      virtual void SetSwitch1D(Int_t s1d = 3) {PARAMETERS.switch_1d = s1d;}
      virtual void SetSwitch3D(Int_t s3d = 0) {PARAMETERS.switch_3d = s3d;}
      virtual void SetSwitchType(Int_t st = 3) {PARAMETERS.switch_type = st;}
      virtual void SetSwitchCoherence(Int_t sc = 0) {PARAMETERS.switch_coherence = sc;}
      virtual void SetSwitchCoulomb(Int_t scol = 2) {PARAMETERS.switch_coulomb = scol;}
      virtual void SetSwitchFermiBose(Int_t sfb = 1) {PARAMETERS.switch_fermi_bose = sfb;}
      
      
      virtual void SetPtRange(Float_t ptmin = 0.1, Float_t ptmax = 0.98) //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) 
                       { MESH.px_min=pxmin;MESH.px_max=pxmax;}

      virtual void SetPyRange(Float_t pymin = -1.0, Float_t pymax = 1.0) 
                       { MESH.py_min=pymin;MESH.py_max=pymax;}

      virtual void SetPzRange(Float_t pzmin = -3.6, Float_t pzmax = 3.6) 
                       { MESH.pz_min=pzmin;MESH.pz_max=pzmax;}


      virtual void SetPhiRange(Float_t phimin = 0.0, Float_t phimax = 360.0) //Angle in degrees
                       { MESH.phi_min=phimin;MESH.phi_max=phimax;}

      virtual void SetEtaRange(Float_t etamin = -1.5, Float_t etamax = 1.5) //Pseudorapidity  !!!!!!!!!
                       { MESH.eta_min=etamin;MESH.eta_max=etamax;}
		      
      virtual void SetNPtBins(Int_t nptbin = 50){MESH.n_pt_bins=nptbin;}
      virtual void SetNPhiBins(Int_t nphibin = 50){MESH.n_phi_bins=nphibin;}
      virtual void SetNEtaBins(Int_t netabin = 50){MESH.n_eta_bins=netabin;}
      
      virtual void SetNPxBins(Int_t npxbin = 20){MESH.n_px_bins=npxbin;}
      virtual void SetNPyBins(Int_t npybin = 20){MESH.n_py_bins=npybin;}
      virtual void SetNPzBins(Int_t npzbin = 70){MESH.n_pz_bins=npzbin;}

      virtual void SetNBins1DFineMesh(Int_t n = 10){ MESH.n_1d_fine=n;}
      virtual void SetBinSize1DFineMesh(Float_t x=0.01){MESH.binsize_1d_fine=x;}
      
      virtual void SetNBins1DCoarseMesh(Int_t n =2 ){MESH.n_1d_coarse =n;}
      virtual void SetBinSize1DCoarseMesh(Float_t x=0.05){MESH.binsize_1d_coarse=x;}
      
      virtual void SetNBins3DFineMesh(Int_t n = 8){MESH.n_3d_fine =n;}
      virtual void SetBinSize3DFineMesh(Float_t x=0.01){MESH.binsize_3d_fine=x;}
      
      virtual void SetNBins3DCoarseMesh(Int_t n = 2){ MESH.n_3d_coarse=n;}
      virtual void SetBinSize3DCoarseMesh(Float_t x=0.08){MESH.binsize_3d_coarse=x;}      
      
      virtual void SetNBins3DFineProjectMesh(Int_t n =3 ){ MESH.n_3d_fine_project=n;}
    
    
    protected:

    private:
    ClassDef(THBTprocessor,0)
  };

#endif

/**************************************************************
COMMENTS APROPOS HBTPROCESSOR FORTRAN CODE:

SUBROUTINE INITIALIZE HAS TO BE OMITTED, IT RESETS (ZEROES) COMMON BLOCKS

SUBROUTINE KINEMATICS HAS TO BE RENAMED EITHER IN MEVSIM OR IN PROCESSOR 
AS THEY HAVE IDENTICAL DECLARATIONS AND SOMTIMES IS USED THE BAD ONE IN CASE BOTH LIBRARIES ARE LOADED

****************************************************************/