PndEvtGenDirect Class Reference

#include <PndEvtGenDirect.h>

Inheritance diagram for PndEvtGenDirect:

Public Member Functions
	PndEvtGenDirect ()
	PndEvtGenDirect (TString particle, TString decfile="", Double_t Mom=0, Long_t Seed=-1, TString defDECAY="", TString defPDL="", Double_t ATarg=0.0)
virtual	~PndEvtGenDirect ()
virtual Bool_t	ReadEvent (FairPrimaryGenerator *primGen)
void	SetStoreTree (Bool_t store=true)
void	SetVerbose (int v=1)

Private Member Functions
	ClassDef (PndEvtGenDirect, 1)
	The mother particle. More...

Private Attributes
Bool_t	fStoreTree
Int_t	verbose
Double_t	fEnergy
Double_t	fMomentum
	Energy of System. More...
EvtGen *	myGenerator
	Momentum of System. More...
EvtStdHep	evtstdhep
	Pointer to the actual EvtGen. More...
EvtId	PART
	The decay tree. More...

Detailed Description

Definition at line 28 of file PndEvtGenDirect.h.

Constructor & Destructor Documentation

PndEvtGenDirect::PndEvtGenDirect

(

)

Default constructor (should not be used)

Definition at line 65 of file PndEvtGenDirect.cxx.

References SetName, and verbose.

                                 {
  SetName("PndEvtGenDirect");
  fStoreTree=false;
  verbose=0;
}

PndEvtGenDirect::PndEvtGenDirect	(	TString	particle,
		TString	decfile = "",
		Double_t	Mom = 0,
		Long_t	Seed = -1,
		TString	defDECAY = "",
		TString	defPDL = "",
		Double_t	ATarg = 0.0
	)

Definition at line 73 of file PndEvtGenDirect.cxx.

References exit(), mp, sqrt(), TString, and val.

                                                                                                                                                 {
  PndEvtGenDirect();

  cout << "<I> PndEvtGenDirect"<<endl;
  cout << "<I> Particle: "<<particle<<endl;
  cout << "<I> decfile: "<<decfile<<endl;
  if(Mom>0) cout << "<I> pbar-Momentum: "<<Mom<<endl;
  if(Mom==0) cout << "<I> Momentum: "<<Mom<<endl;
  if(Mom<0) cout << "<I> CMS energy: "<<Mom<<endl;
  cout << "<I> Rnd Seed: "<<Seed<<endl;
  TString work = getenv("VMCWORKDIR");
  if (defaultDECAY=="") defaultDECAY = work + "/pgenerators/EvtGen/EvtGen/Private/DECAY.DEC";
  if (defaultPDL=="") defaultPDL = work + "/pgenerators/EvtGen/EvtGen/Private/evt.pdl";

  if ( (particle.Contains("pbarASystem") || particle.Contains("pASystem")) && (ATarg<3.||ATarg>238.))
  {
    cerr <<"******  FATAL ERROR: nuclear target mass number MUST be between 3 and 238! ******"<<endl;
    exit(0);
  }

  //Initialize the generator - read in the decay table and particle properties

  EvtRandomEngine* myRandomEngine=0;

  // Make sure that the seed is always set if default value (-1) is used, JGM, August 2011
  if (Seed<0)
  {
    Seed = gRandom->GetSeed();
    cout << "<I> Rnd Seed changed to " << Seed << endl;
  }
  myRandomEngine=new EvtRootRandomEngine(Seed);


  // Set up the default external generator list: Photos, Pythia and/or Tauola ONLY if available
#if EVTGEN_EXTERNAL
  EvtExternalGenList genList;
  EvtAbsRadCorr* radCorrEngine = genList.getPhotosModel();
  std::list<EvtDecayBase*> extraModels = genList.getListOfModels();

  // Create the EvtGen generator object
  myGenerator=new EvtGen(defaultDECAY,defaultPDL,myRandomEngine,
                         radCorrEngine, &extraModels);
#else
  //If you don't want to use external generators, use the following:
  myGenerator=new EvtGen(defaultDECAY,defaultPDL,myRandomEngine);
#endif

  //If I wanted a user decay file, I would read it in now.
  if(decfile!="") myGenerator->readUDecay(decfile.Data());

  PART=EvtPDL::getId(std::string(particle.Data()));

  double val=-3.0969;
  fMomentum = 0.0;
  fEnergy = 0.0;
  double mp=0.938272;
  double md=1.875613;
  double mu=0.931494;
  double mA=ATarg*mu;

  if ( particle.Contains("pbarp") || particle.Contains("pbard") || particle.Contains("pbarA") || particle.Contains("pp") || particle.Contains("pd") || particle.Contains("pA") ) {
    if (Mom==0) {
      cerr <<"\033[5m\033[31m -E  ******  FATAL ERROR: <particle> is '" << particle.Data() << "'; MUST give pbar momentum or cms energy!\033[0m"<<endl;
      exit(0);
    } else {
      val=Mom;
    }
  } else {
    if(PART.getId()==-1) {
      cerr << "Particle \""<<particle<<"\" is unknown!!!"<<endl<<"Check your Macro for spelling mistake."<<endl;
      exit(0);
    }
    val=-EvtPDL::getMass(PART);
  }

  // val is the momentum of the beam
  if (val>0) {
    fMomentum = val;
    if ( particle.Contains("pbarpSystem") || particle.Contains("ppSystem")) fEnergy = mp+sqrt(fMomentum*fMomentum+mp*mp);
    if ( particle.Contains("pbardSystem") || particle.Contains("pdSystem")) fEnergy = md+sqrt(fMomentum*fMomentum+mp*mp);
    if ( particle.Contains("pbarASystem") || particle.Contains("pASystem")) fEnergy = mA+sqrt(fMomentum*fMomentum+mp*mp);
  } else {  //val is -E_cm
    val=-val;
    fEnergy = val*val/(2*mp);
    fMomentum = sqrt(fEnergy*fEnergy-val*val);
  }

  cout <<"\n############# Generating with following conditions:\n\n";
  cout <<"incident 4-mom : ("<<fEnergy<<", 0, 0, "<<fMomentum<<"), m = "<<sqrt(fEnergy*fEnergy-fMomentum*fMomentum)<<endl;
  cout <<"\n######################\n\n"<<endl;

}

PndEvtGenDirect::~PndEvtGenDirect ( )

virtual

Destructor

Definition at line 168 of file PndEvtGenDirect.cxx.

                                  {

}

Member Function Documentation

PndEvtGenDirect::ClassDef ( PndEvtGenDirect  , 1    )

private

The mother particle.

Bool_t PndEvtGenDirect::ReadEvent ( FairPrimaryGenerator *  primGen )

virtual

Generate one event using EVT

Parameters

primGen

pointer to the FairPrimaryGenerator

Definition at line 174 of file PndEvtGenDirect.cxx.

References Bool_t, C(), Double_t, fT, fX, fY, fZ, i, printf(), and verbose.

                                                               {
  static Int_t evtnr=0;
  // Loop to create nEvents, starting from an Upsilon(4S)
  // Set up the parent particle
  EvtParticle *parent;

  EvtVector4R pInit(fEnergy,  0.0000, -0.0000,  fMomentum);
  parent=EvtParticleFactory::particleFactory(PART,pInit);
  parent->setDiagonalSpinDensity();

  // Generate the event
  myGenerator->generateDecay(parent);
  // Write out the results

  // we should better use the EvtParticle properties here... and kick the evtstd structure!!!!
  evtstdhep.init();
  parent->makeStdHep(evtstdhep);

  Bool_t plotflag;
  plotflag=false;
  //print out some status info
  if (verbose>1 ||(verbose==1 && (evtnr<10 || ((evtnr+1)%100)==0))) {
    cout << "PndEvtGenDirect::ReadEvent "<<evtnr <<" "<<fEnergy<<" "<<fMomentum << endl;
    parent->printParticle();
    // Write out the results
    EvtHepMCEvent theEvent;
    theEvent.constructEvent(parent);
    HepMC::GenEvent* genEvent = theEvent.getEvent();
    genEvent->print(std::cout);

    //    report(INFO,"EvtGen") << "event Number\t"<< evtnr << evtstdhep << endl;
    //     cout << evtnr << "\t" << evtstdhep.getNPart();
    //     cout <<evtstdhep<<endl;

    cout <<"==== now compare ==="<<endl;
    plotflag=true;
  }

  // Write the output

  Int_t  npart;
  Double_t fX, fY, fZ, fT;
  double Px,Py,Pz, fE; // ,Pm[1000],Wh[1000];
  int Id;
  npart=evtstdhep.getNPart();
  EvtVector4R vxyz,pxyz;

  for(Int_t i=0; i<npart; i++) {
    Int_t nFD, nLD;
    // add track
    nFD=evtstdhep.getFirstDaughter(i);
    nLD=evtstdhep.getLastDaughter(i);
    if(fStoreTree ||(nFD==-1 && nLD==-1))
    {
      Id=evtstdhep.getStdHepID(i);
      vxyz=evtstdhep.getX4(i);
      pxyz=evtstdhep.getP4(i);
      fT=vxyz.get(0)/(1000*TMath::C());  //mm - > s conversion
      fX=vxyz.get(1)/10.; // mm -> cm conversion
      fY=vxyz.get(2)/10.; // mm -> cm conversion
      fZ=vxyz.get(3)/10.; // mm -> cm conversion
      fE=pxyz.get(0);
      Px=pxyz.get(1);
      Py=pxyz.get(2);
      Pz=pxyz.get(3);
      if(plotflag) printf("- I -: new particle %d at: %f, %f, %f (%f)-> %f %f %f (%f) ID %d ##Daughters %d %d Mothers %d %d\n", i,
                            fX, fY, fZ, fT,Px, Py, Pz, fE, Id, nFD, nLD,evtstdhep.getFirstMother(i),evtstdhep.getLastMother(i));
      if(fStoreTree) {
        primGen->AddTrack(Id, Px, Py, Pz, fX, fY, fZ, evtstdhep.getFirstMother(i),(nFD==-1 && nLD==-1),fE,fT);
      } else {
        primGen->AddTrack(Id, Px, Py, Pz, fX, fY, fZ,-1,true,fE,fT);// default -1, true
      }
    }
  }
  if(plotflag) cout <<"==== compare end ==="<<endl;

  parent->deleteTree();

  evtnr++;

  return kTRUE;

}

void PndEvtGenDirect::SetStoreTree ( Bool_t  store = true )

inline

Definition at line 47 of file PndEvtGenDirect.h.

References fStoreTree, and store().

Referenced by prod_fsim(), prod_sim(), QAmacro_evtgen_1(), QAmacro_evtgen_2(), quickfsimana(), run_sim(), run_sim_complete(), run_sim_stt_evt(), run_sim_stt_evt_fix(), run_sim_sttcombi_dpm(), run_sim_sttcombi_evtgen(), run_sim_tpc_evt(), run_sim_tpccombi(), sim(), sim_complete(), sim_complete_pbarA(), sim_complete_runs(), sim_complete_vis(), sim_evtgen(), sim_filter_ex1(), sim_filter_ex2(), sim_filter_inv_mass(), sim_rich(), sim_scit(), sim_test(), sim_toy(), simfast(), simfast_cmp(), simfast_opt(), tb_sim_complete(), tb_sim_scit(), tut_fastsim(), and PndMasterRunSim::UseEvtGenGenerator().

{ fStoreTree=store;};

void PndEvtGenDirect::SetVerbose ( int  v = 1 )

inline

Definition at line 48 of file PndEvtGenDirect.h.

References v, and verbose.

Referenced by runSim().

{verbose=v;};

Member Data Documentation

EvtStdHep PndEvtGenDirect::evtstdhep

private

Pointer to the actual EvtGen.

Definition at line 61 of file PndEvtGenDirect.h.

Double_t PndEvtGenDirect::fEnergy

private

P_lab(GeV/c)

Definition at line 57 of file PndEvtGenDirect.h.

Double_t PndEvtGenDirect::fMomentum

private

Energy of System.

Definition at line 58 of file PndEvtGenDirect.h.

Bool_t PndEvtGenDirect::fStoreTree

private

Definition at line 48 of file PndEvtGenDirect.h.

Referenced by SetStoreTree().

EvtGen* PndEvtGenDirect::myGenerator

private

Momentum of System.

Definition at line 60 of file PndEvtGenDirect.h.

EvtId PndEvtGenDirect::PART

private

The decay tree.

Definition at line 62 of file PndEvtGenDirect.h.

Int_t PndEvtGenDirect::verbose

private

Definition at line 53 of file PndEvtGenDirect.h.

Referenced by SetVerbose().

---

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

• PndEvtGenDirect.h
• PndEvtGenDirect.cxx