PndDpmDirect Class Reference

#include <PndDpmDirect.h>

Inheritance diagram for PndDpmDirect:

Public Member Functions
	PndDpmDirect ()
	PndDpmDirect (Double_t Mom, Int_t Mode, Long_t Seed=-1)
	PndDpmDirect (Double_t Mom, Int_t Mode, Long_t Seed, Double_t ThtMin)
	PndDpmDirect (Double_t Mom, Int_t Mode, Double_t Rsigma, TF1 *DensityFunction, Long_t Seed=-1, Double_t ThtMin=0.001)
virtual	~PndDpmDirect ()
virtual Bool_t	ReadEvent (FairPrimaryGenerator *primGen)
void	SetUnstable (int pdg)
void	SetStable (int pdg)

Private Member Functions
	ClassDef (PndDpmDirect, 1)

Private Attributes
double	fMom
double	fMode
double	fSeed
int	fGasmode
double	fRsigma
double	fThtMin
TF1 *	fDensityFunction

Detailed Description

Definition at line 25 of file PndDpmDirect.h.

Constructor & Destructor Documentation

PndDpmDirect::PndDpmDirect

(

)

Default constructor (should not be used)

Definition at line 45 of file PndDpmDirect.cxx.

                           {

}

PndDpmDirect::PndDpmDirect	(	Double_t	Mom,
		Int_t	Mode,
		Long_t	Seed = -1
	)

Standard constructor

Parameters

Mom	in GeV/C
Mode	= 0. - No elastic scattering, only inelastic
Mode	= 1. - Elastic and inelastic interactions
Mode	= 2. - Only elastic scattering, no inelastic one

Definition at line 52 of file PndDpmDirect.cxx.

References Double_t, and CAMath::Log().

                                                                {

//
// Calculate ThtMin first. For this we make a cut-off on the value of -t of 1e-2 GeV^2 (~100 MeV/c momentum)
// This estimated from a parametrization found in thesis of Thomas Wuerschig (figure 6.4, page 121): 
// Roughly: 0.4 deg at 15 GeV/c and 4 deg at 1.5 GeV/c, lineair interpolation in double log-scale.
//

  Double_t logangle = TMath::Log(0.4)+(TMath::Log(15.)-TMath::Log(Mom))*(TMath::Log(4)-TMath::Log(0.4))/(TMath::Log(15)-TMath::Log(1.5));
  Double_t ThtMin = TMath::Exp(logangle);

  PndDpmDirect(Mom, Mode, Seed, ThtMin) ;
}

PndDpmDirect::PndDpmDirect	(	Double_t	Mom,
		Int_t	Mode,
		Long_t	Seed,
		Double_t	ThtMin
	)

Definition at line 66 of file PndDpmDirect.cxx.

References a, fDensityFunction, fSeed, and init1_().

                                                                                 {
   fMom    = Mom;
   fMode   = Mode;
   fSeed   = Seed;
   fThtMin = ThtMin;
      
   if (fSeed < 0)
   {
     Long_t iSeed = gRandom->GetSeed();
     int a = iSeed/100000;   
     fSeed=iSeed - a*100000 + a/100000.;
   }
   
   fGasmode = 0;
   fRsigma = 0.;
   fDensityFunction = new TF1();

   cout << "<I> PndDpmDirect initialization" << endl;
   cout << "<I> Momentum  = " << fMom << endl;
   cout << "<I> Seed      = " << fSeed << endl;
   cout << "<I> Mode      = " << fMode << endl; 
   cout << "<I> Theta min = " << fThtMin<<endl;

   init1_(&fMom,&fSeed,&fMode, &fThtMin);     // init the DPM generator  
}

PndDpmDirect::PndDpmDirect	(	Double_t	Mom,
		Int_t	Mode,
		Double_t	Rsigma,
		TF1 *	DensityFunction,
		Long_t	Seed = -1,
		Double_t	ThtMin = 0.001
	)

Definition at line 94 of file PndDpmDirect.cxx.

References a, fDensityFunction, fSeed, and init1_().

                                                                                                                        {
   fMom    = Mom;
   fMode   = Mode;
   fSeed   = Seed;
   fThtMin = ThtMin;
   
   if (fSeed < 0)
   {
         Long_t iSeed = gRandom->GetSeed();
     int a = iSeed/100000;   
     fSeed=iSeed - a*100000 + a/100000.;
   }
   
   fGasmode = 1;
   fRsigma = Rsigma;
   fDensityFunction = DensityFunction;

   cout << "<I> PndDpmDirect initialization" << endl;
   cout << "<I> Momentum  = " << fMom << endl;
   cout << "<I> Seed      = " << fSeed << endl;
   cout << "<I> Mode      = " << fMode << endl;
   cout << "<I> Gasmode   = " << fGasmode << endl;
   cout << "<I> Theta min = " << fThtMin<<endl;

   init1_(&fMom,&fSeed,&fMode, &fThtMin);     // init the DPM generator     
}

PndDpmDirect::~PndDpmDirect ( )

virtual

Destructor

Definition at line 123 of file PndDpmDirect.cxx.

                            {

}

Member Function Documentation

PndDpmDirect::ClassDef ( PndDpmDirect  , 1    )

private

Bool_t PndDpmDirect::ReadEvent ( FairPrimaryGenerator *  primGen )

virtual

Generate one event using DPM

Parameters

primGen

pointer to the FairPrimaryGenerator

Definition at line 143 of file PndDpmDirect.cxx.

References Double_t, dpm_gen_(), fDensityFunction, fSeed, fX, fY, fZ, and i.

                                                            {

  int  npart, i;        
        Double_t fX, fY, fZ, radius;
  double Px[1000],Py[1000],Pz[1000];  //,E[1000],Pm[1000],Wh[1000];  //[R.K. 01/2017] unused variables
  int Id[1000];
  
  double Generator=0.; //  Format in which events are produced (0=pythia, 1=pluto)

  //Double_t weight = 1.0; //[R.K. 01/2017] unused variable?
  //Int_t activeCnt=0; //[R.K. 01/2017] unused variable?

  // run generator
        dpm_gen_(&Generator, &fSeed);
  
        // Loop over all produced particles 
   npart = lujets_.n;
      
   for (i= 0; i< npart; ++i) {        

                Id[i]=lujets_.k[i+1000];
                Px[i]=lujets_.p[i];
    Py[i]=lujets_.p[i+1000];
    Pz[i]=lujets_.p[i+2000];
    //Pm[i]=lujets_.p[i+4000];
    //E[i]=lujets_.p[i+3000];
    //Wh[i]=1.0;

    /* Check if fGasmode is set */
    fX = 0.;
    fY = 0.;
    fZ = 0.;
                if (fGasmode == 1) {
     
        // define position of track start
        // Random 2D point in a circle of radius r (simple beamprofile)         
        radius = gRandom->Gaus(0,fRsigma);
                        gRandom->Circle(fX, fY, radius);
        
                        // calculate fZ according to some (probability) density function of the gas
        fZ=fDensityFunction->GetRandom();
                        
    }
                
                // add track
                //printf("- I -: new particle at: %f, %f, %f ...\n", fX, fY, fZ);
                primGen->AddTrack(Id[i], Px[i], Py[i], Pz[i], fX, fY, fZ);

   }

  return kTRUE;

}

void PndDpmDirect::SetStable

(

int

pdg

)

Definition at line 127 of file PndDpmDirect.cxx.

References chstatus_(), and status.

Referenced by prod_fsim().

{
        int status=1;
        chstatus_(&pdg, &status);
}

void PndDpmDirect::SetUnstable

(

int

pdg

)

Definition at line 134 of file PndDpmDirect.cxx.

References chstatus_(), and status.

Referenced by prod_fsim(), quickfsimana(), simfast(), simfast_opt(), and tut_fastsim().

{
        int status=0;
        chstatus_(&pdg, &status);
}

Member Data Documentation

TF1* PndDpmDirect::fDensityFunction

private

Definition at line 73 of file PndDpmDirect.h.

int PndDpmDirect::fGasmode

private

Definition at line 69 of file PndDpmDirect.h.

double PndDpmDirect::fMode

private

0. - No elastic scattering, only inelastic

1. - Elastic and inelastic interactions
2. - Only elastic scattering, no inelastic one

Definition at line 67 of file PndDpmDirect.h.

double PndDpmDirect::fMom

private

P_lab(GeV/c)

Definition at line 61 of file PndDpmDirect.h.

double PndDpmDirect::fRsigma

private

Definition at line 70 of file PndDpmDirect.h.

double PndDpmDirect::fSeed

private

Definition at line 68 of file PndDpmDirect.h.

double PndDpmDirect::fThtMin

private

Definition at line 71 of file PndDpmDirect.h.

---

The documentation for this class was generated from the following files:

• PndDpmDirect.h
• PndDpmDirect.cxx