PndPythia6Direct Class Reference

#include <PndPythia6Direct.h>

Inheritance diagram for PndPythia6Direct:

Public Member Functions
	PndPythia6Direct ()
virtual	~PndPythia6Direct ()
Bool_t	ReadEvent (FairPrimaryGenerator *)
void	Print ()
virtual Bool_t	Init ()
void	SetMom (Double_t mom)
void	SetStoreTree (Bool_t store=true)
void	SetVerbose (int v=1)
void	SetPandaDefaults ()
TPythia6 *	GetTPythia6 ()

Protected Member Functions
	ClassDef (PndPythia6Direct, 1)

Protected Attributes
Double_t	fMom
bool	fStoreTree
int	fVerbose

Private Attributes
TPythia6 *	fPythia
TClonesArray *	fParticleList
TMCParticle *	fParticle

Detailed Description

Definition at line 16 of file PndPythia6Direct.h.

Constructor & Destructor Documentation

PndPythia6Direct::PndPythia6Direct

(

)

default constructor

Definition at line 12 of file PndPythia6Direct.cxx.

References fMom, fPythia, fStoreTree, and SetPandaDefaults().

{
  fMom = 15.;
  fStoreTree=false;
  fPythia=TPythia6::Instance();
  SetPandaDefaults();
}

PndPythia6Direct::~PndPythia6Direct ( )

virtual

destructor

Definition at line 51 of file PndPythia6Direct.cxx.

References fPythia.

{
  if(fPythia!=NULL) delete fPythia;
}

Member Function Documentation

PndPythia6Direct::ClassDef ( PndPythia6Direct  , 1    )

protected

TPythia6* PndPythia6Direct::GetTPythia6 ( )

inline

Definition at line 42 of file PndPythia6Direct.h.

References fPythia.

                            {
      return fPythia;
    };

Bool_t PndPythia6Direct::Init ( )

virtual

Definition at line 38 of file PndPythia6Direct.cxx.

References fMom, fParticleList, and fPythia.

{
  std::cout<<"PndPythia6Direct::Init(): Beam Momentum "<<fMom<<std::endl;

  fPythia->Initialize("FIXT", "pbar", "p", fMom);
  fParticleList = (TClonesArray*)fPythia->GetListOfParticles();

  return kTRUE;
}

void PndPythia6Direct::Print

(

)

Definition at line 132 of file PndPythia6Direct.cxx.

                             {
  //fPythia.settings.listAll();
}

Bool_t PndPythia6Direct::ReadEvent

(

FairPrimaryGenerator *

primGen

)

public method ReadEvent

Definition at line 58 of file PndPythia6Direct.cxx.

References fParticle, fParticleList, fPythia, fStoreTree, fT, fVerbose, fX, fY, fZ, i, printf(), and status.

{
  static int evtn=0;
  if(fVerbose>0) std::cout<<"PndPythia6Direct::ReadEvent() No. "<<evtn<<std::endl;
  fPythia->GenerateEvent();
  if(fVerbose>2) std::cout<<"PndPythia6Direct::ReadEvent() Generation done."<<std::endl;
  double fX, fY, fZ, fT;
  double Px,Py,Pz, fE; // ,Pm[1000],Wh[1000];
  int Id,nFD,nLD,nMO,status;
  bool dotracking;
  //fParticleList=(TClonesArray*)fPythia->ImportParticles();
  int npart=fParticleList->GetEntriesFast();
  if(fVerbose>0) std::cout<<"PndPythia6Direct::ReadEvent() We have "<<npart<<" tracks."<<std::endl;
  for(int i=0; i<npart; ++i)
  {
    if(fVerbose>4) std::cout<<"PndPythia6Direct::ReadEvent() particle No. "<<i<<std::endl;
    fParticle=(TMCParticle*)fParticleList->At(i);
    if(fVerbose>4) std::cout<<"PndPythia6Direct::ReadEvent() particle pointer "<<fParticle<<std::endl;
    if(fParticle==NULL) continue;
    // add track
    nFD=fParticle->GetFirstChild();
    nLD=fParticle->GetLastChild();
    nMO=fParticle->GetParent();
    status=fParticle->GetKS();
      //Id=TDatabasePDG::Instance()->GetParticle(fParticle->GetName())->PdgCode();
    Id=fParticle->GetKF(); // "Flavour code" is pdg code
    if(fVerbose>4) {
      std::cout<<"PndPythia6Direct::ReadEvent() Particle:"
               //<<" pdg="<<TDatabasePDG::Instance()->GetParticle(fParticle->GetName())->PdgCode()
               <<" name=\""<<fParticle->GetName()<<"\""
               <<" nFD= "<<nFD
               <<" nLD= "<<nLD
               <<" nMO= "<<nMO
               <<" pdg= "<<Id
               <<" status= "<<status
               <<"\n                                       "
               <<" Time()*1000="<<fParticle->GetTime()*1000
               <<" Vx()*0.1= "<<fParticle->GetVx()*0.1
               <<" Vy()*0.1= "<<fParticle->GetVy()*0.1
               <<" Vz()*0.1= "<<fParticle->GetVz()*0.1
               <<"\n                                       "
               <<" Energy()= "<<fParticle->GetEnergy()
               <<" Px()= "<<fParticle->GetPx()
               <<" Py()= "<<fParticle->GetPy()
               <<" Pz()= "<<fParticle->GetPz()
               <<std::endl;
    }

    if(status==1 || fStoreTree )         //fStoreTree && status==11 && nFD!=nLD)
    {
      fT=fParticle->GetTime()*1000; //[mm/c]->s  with c in [m/s]           //   /(1000*TMath::C());  //mm - > s conversion
      fX=fParticle->GetVx()*0.1; // mm -> cm conversion
      fY=fParticle->GetVy()*0.1; // mm -> cm conversion
      fZ=fParticle->GetVz()*0.1; // mm -> cm conversion
      fE=fParticle->GetEnergy();
      Px=fParticle->GetPx();
      Py=fParticle->GetPy();
      Pz=fParticle->GetPz();
      if(fVerbose>0&&(evtn<10||(evtn%100)==0)) printf("- I -: new particle %d at: %f, %f, %f (%f)-> %f %f %f (%f) ID %d ##Daughters %d %d Mother %d\n",
            i,fX, fY, fZ, fT,Px, Py, Pz, fE, Id, nFD, nLD, nMO);
      if(fStoreTree) {
        // here we write all particles to the stack, but do processing only for final state particles
        dotracking = (status==1);
        primGen->AddTrack(Id, Px, Py, Pz, fX, fY, fZ, nMO,dotracking,fE,fT);
      } else {
        primGen->AddTrack(Id, Px, Py, Pz, fX, fY, fZ,-1,true,fE,fT);// set mother id to default: -1, do tracking: true
      }
    }
  }
  evtn++;
  return kTRUE;
}

void PndPythia6Direct::SetMom ( Double_t  mom )

inline

Definition at line 32 of file PndPythia6Direct.h.

References fMom, and mom.

                              {
      fMom = mom;
    };

void PndPythia6Direct::SetPandaDefaults

(

)

Definition at line 22 of file PndPythia6Direct.cxx.

References fPythia.

Referenced by PndPythia6Direct().

                                        {
  std::cout<<"PndPythia6Direct::SetPandaDefaults()"<<std::endl;
  // Those are the FORTRAN settings. Let's access them via TPythia6
  //c----------  Choice of the nessesary process --------------------
  fPythia->SetMSEL(0);       // full User Access
  fPythia->SetMSUB(1,1);     // qQ-->A*/Z0*
  //c-------------------  Nessesary limits for PANDA ENERGY region ---------
  fPythia->SetPARP(2,0.5);   // Lowest inv. Mass (D=10.GeV)
  fPythia->SetPARP(111,0.1); // Minimum invar. mass of the remnant hadronic system (D=2.GeV)
  fPythia->SetCKIN(1,1.);    // limits for shat
  //fPythia->SetCKIN(2,-1.);   // no upper limit for S-HAT
  fPythia->SetCKIN(5,0.1);   // low cut limits for PT_-range
  fPythia->SetCKIN(6,0.1);   // low cut limits for intermediate mass
}

void PndPythia6Direct::SetStoreTree ( Bool_t  store = true )

inline

Definition at line 35 of file PndPythia6Direct.h.

References fStoreTree, and store().

                                         {
      fStoreTree=store;
    };

void PndPythia6Direct::SetVerbose ( int  v = 1 )

inline

Definition at line 38 of file PndPythia6Direct.h.

References fVerbose, and v.

                                    {
      fVerbose=v;
    };

Member Data Documentation

Double_t PndPythia6Direct::fMom

protected

Definition at line 53 of file PndPythia6Direct.h.

Referenced by Init(), PndPythia6Direct(), and SetMom().

TMCParticle* PndPythia6Direct::fParticle

private

Definition at line 50 of file PndPythia6Direct.h.

Referenced by ReadEvent().

TClonesArray* PndPythia6Direct::fParticleList

private

Definition at line 49 of file PndPythia6Direct.h.

Referenced by Init(), and ReadEvent().

TPythia6* PndPythia6Direct::fPythia

private

Definition at line 44 of file PndPythia6Direct.h.

Referenced by GetTPythia6(), Init(), PndPythia6Direct(), ReadEvent(), SetPandaDefaults(), and ~PndPythia6Direct().

bool PndPythia6Direct::fStoreTree

protected

Definition at line 54 of file PndPythia6Direct.h.

Referenced by PndPythia6Direct(), ReadEvent(), and SetStoreTree().

int PndPythia6Direct::fVerbose

protected

Definition at line 55 of file PndPythia6Direct.h.

Referenced by ReadEvent(), and SetVerbose().

---

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

• PndPythia6Direct.h
• PndPythia6Direct.cxx