FixedTargetGenerator Class Reference

#include <FixedTargetGenerator.h>

Inheritance diagram for FixedTargetGenerator:

Collaboration diagram for FixedTargetGenerator:

Public Member Functions
	FixedTargetGenerator ()
virtual	~FixedTargetGenerator ()
Bool_t	ReadEvent (FairPrimaryGenerator *)
void	SetParameters (char *)
void	Print ()
virtual Bool_t	Init ()
Bool_t	InitForCharmOrBeauty (TString fInName, Int_t nev, Double_t npots=5E13, Int_t nStart=0)
void	SetMom (Double_t mom)
void	UseRandom1 ()
void	UseRandom3 ()
void	SetCharmTarget (Bool_t charmtarget=true)
void	SetTarget (TString s, Double_t x, Double_t y)
void	SetBoost (Double_t f)
void	SetG4only ()
void	SetTauOnly ()
void	SetJpsiMainly ()
void	SetOnlyMuons ()
void	SetDrellYan ()
void	SetPhotonCollision ()
void	WithEvtGen ()
void	SetChibb (Double_t x)
void	SetChicc (Double_t x)
void	SetSeed (Double_t seed)
void	SetHeartBeat (Int_t x)
void	SetEnergyCut (Float_t emax)
void	SetDebug (Bool_t x)
void	SetOpt4DP (TNtuple *t)
Double_t	GetPotForCharm ()
Pythia8::Pythia *	GetPythia ()
Pythia8::Pythia *	GetPythiaN ()

Protected Member Functions
	ClassDef (FixedTargetGenerator, 2)

Protected Attributes
Double_t	fMom
Bool_t	fUseRandom1
Bool_t	fUseRandom3
Double_t	fSeed
Double_t	EMax
Double_t	fBoost
Double_t	chicc
Double_t	chibb
Double_t	wspill
Double_t	nrpotspill
Int_t	nEvents
Int_t	nEntry
Int_t	pot
Int_t	nDsprim
Int_t	ntotprim
Bool_t	tauOnly
Bool_t	JpsiMainly
Bool_t	DrellYan
Bool_t	PhotonCollision
Bool_t	G4only
Bool_t	setByHand
Bool_t	Debug
Bool_t	withEvtGen
Bool_t	OnlyMuons
Bool_t	fcharmtarget
FairLogger *	fLogger
Pythia8::Pythia *	fPythiaN
	don't make it persistent, magic ROOT command
Pythia8::Pythia *	fPythiaP
EvtGenDecays *	evtgenN
EvtGenDecays *	evtgenP
GenieGenerator *	fMaterialInvestigator
Bool_t	withNtuple
TNtuple *	fNtuple
	special option for Dark Photon physics studies
TString	targetName
TString	Option
Double_t	xOff
Double_t	yOff
Double_t	start [3]
Double_t	end [3]
Double_t	bparam
Double_t	mparam [10]
Double_t	startZ
Double_t	endZ
Double_t	maxCrossSection
TFile *	fin
TNtuple *	nTree
Float_t	n_id
Float_t	n_px
Float_t	n_py
Float_t	n_pz
Float_t	n_M
Float_t	n_E
Float_t	n_mpx
Float_t	n_mpy
Float_t	n_mpz
Float_t	n_mE
Float_t	n_mid
Float_t	ck
Int_t	heartbeat

Private Attributes
Pythia8::RndmEngine *	fRandomEngine

Detailed Description

Definition at line 15 of file FixedTargetGenerator.h.

Constructor & Destructor Documentation

FixedTargetGenerator()

FixedTargetGenerator::FixedTargetGenerator

(

)

default constructor

Definition at line 23 of file FixedTargetGenerator.cxx.

{
  fUseRandom1 = kFALSE; 
  fUseRandom3 = kTRUE;
  fSeed = 0;
  fMom        = 400;  // proton
  targetName = ""; 
  fcharmtarget = false;
  xOff=0; yOff=0;
  tauOnly = false;
  JpsiMainly = false;
  DrellYan = false;
  PhotonCollision = false;
  G4only = false;
  OnlyMuons = false;
  maxCrossSection = 0.;
  EMax = 0;
  fBoost = 1.;
  withEvtGen = kFALSE;
  withNtuple = kFALSE;
  chicc=1.7e-3;     //prob to produce primary ccbar pair/pot
  chibb=1.6e-7;     //prob to produce primary bbbar pair/pot
  nrpotspill = 5.E13; // nr of protons / spill
  setByHand = kFALSE;
  Debug = kFALSE;
  Option = "Primary";
  wspill = 1.; // event weight == 1 for primary events
  heartbeat = 1000;
}

~FixedTargetGenerator()

FixedTargetGenerator::~FixedTargetGenerator ( )

virtual

destructor

Definition at line 296 of file FixedTargetGenerator.cxx.

{
}

Member Function Documentation

ClassDef()

FixedTargetGenerator::ClassDef ( FixedTargetGenerator  , 2    )

protected

GetPotForCharm()

Double_t FixedTargetGenerator::GetPotForCharm ( )

inline

Definition at line 53 of file FixedTargetGenerator.h.

{return nrpotspill/wspill;}

GetPythia()

Pythia8::Pythia * FixedTargetGenerator::GetPythia ( )

inline

Definition at line 54 of file FixedTargetGenerator.h.

{return fPythiaP;}

GetPythiaN()

Pythia8::Pythia * FixedTargetGenerator::GetPythiaN ( )

inline

Definition at line 55 of file FixedTargetGenerator.h.

{return fPythiaN;}

Init()

Bool_t FixedTargetGenerator::Init ( )

virtual

Definition at line 99 of file FixedTargetGenerator.cxx.

{
  fPythiaP =  new Pythia8::Pythia(); // pythia needed also for G4only, to copy 2mu BRs
  if (Option == "Primary" && !G4only){
   fPythiaN =  new Pythia8::Pythia();
  }else if (Option != "charm" && Option != "beauty" && !G4only) {
   LOG(INFO) << "FixedTargetGenerator::Init Option not known "<< Option.Data() <<", abort";
  }
  if (fUseRandom1) fRandomEngine = new PyTr1Rng();
  if (fUseRandom3) fRandomEngine = new PyTr3Rng();
  std::vector<int> r = { 221, 221, 223, 223,   113, 331, 333};
  std::vector<int> c = {  6 ,  7,   5 ,  7,     5,   6,   9}; // decay channel mumu mumuX
 
  if (Option == "Primary" && !G4only){
   fPythiaP->settings.mode("Beams:idB",  2212);
   fPythiaN->settings.mode("Beams:idB",  2112);
  }else{
   fPythiaP->readString("ProcessLevel:all = off");
  }
  std::array<Pythia8::Pythia*,2> plist = {fPythiaP,fPythiaN};
  Int_t pcount = 0;
  for(const auto& fPythia : plist) {
   if (pcount > 0 && (Option != "Primary" || G4only)){continue;}
   pcount+=1; 
   fPythia->setRndmEnginePtr(fRandomEngine);
   fPythia->settings.mode("Random:seed",fSeed);
   fPythia->settings.mode("Next:numberCount",heartbeat);
   if (Option == "Primary"){
    fPythia->settings.mode("Beams:idA",  2212);
    fPythia->settings.mode("Beams:frameType",  2);
    fPythia->settings.parm("Beams:eA",fMom);
    fPythia->settings.parm("Beams:eB",0.);
   }
   if (JpsiMainly){
// use this for all onia productions
     fPythia->readString("Onia:all(3S1) = on");
     fPythia->readString("443:new  J/psi  J/psi  3   0   0    3.09692    0.00009    3.09602    3.09782  0.   1   1   0   1   0");
     fPythia->readString("443:addChannel = 1   1.    0      -13       13");
     fPythia->readString("553:new  Upsilon  Upsilon  3   0   0    9.46030    0.00005    9.45980    9.46080  0.00000e+00   0   1   0   1   0");
     fPythia->readString("553:addChannel = 1   1.    0      -13       13");
   }
   if (DrellYan){
     fPythia->readString("WeakSingleBoson:ffbar2gmZ = on");
//   Z0/gamma  -> mu+ mu-
     fPythia->readString("23:onMode = off");
     fPythia->readString("23:onIfAll = 13 13");
     fPythia->readString("23:mMin = 0.5");
     fPythia->readString("PhaseSpace:mHatMin = 0.5");
     fPythia->readString("PartonLevel:FSR = on");
     fPythia->readString("PDF:pSet = 4");
   }
   if (PhotonCollision){
     fPythia->readString("PhotonCollision:gmgm2mumu = on");
   }
   if (!DrellYan && !PhotonCollision){
     fPythia->readString("SoftQCD:inelastic = on");
     fPythia->readString("PhotonCollision:gmgm2mumu = on");
     fPythia->readString("PromptPhoton:all = on");
     fPythia->readString("WeakBosonExchange:all = on");
     fPythia->readString("WeakSingleBoson:all = on");
    }
   if (tauOnly){
    fPythia->readString("431:new  D_s+  D_s-  1   3   0    1.96849    0.00000    0.00000    0.00000  1.49900e-01   0   1   0   1   0");
    fPythia->readString("431:addChannel = 1   0.0640000    0      -15       16");
   }
 
// find all long lived particles in pythia
   Int_t n = 1;
   while(n!=0){
    n = fPythia->particleData.nextId(n);
#ifndef PYTHIA8_V
    ParticleDataEntry* p = fPythia->particleData.particleDataEntryPtr(n);
#else
#if PYTHIA8_V<8309
    ParticleDataEntry* p = fPythia->particleData.particleDataEntryPtr(n);
#else
    Pythia8::ParticleDataEntryPtr p = fPythia->particleData.particleDataEntryPtr(n);
#endif
#endif
    if (p->tau0()>1){
     std::string particle = std::to_string(n)+":mayDecay = false";
     fPythia->readString(particle);
     LOG(INFO) << "FixedTargetGenerator::Init Made "<< p->name().c_str() <<" stable for Pythia, should decay in Geant4";
    }
   } 
// boost branching fraction of rare di-muon decays
//                       eta  omega rho0  eta' phi   
   if (fBoost != 1.){
    LOG(INFO) << "FixedTargetGenerator::Init Rescale BRs of dimuon decays in Pythia: "<<fBoost;
    for (unsigned int i=0; i<r.size(); ++i) {  
#ifndef PYTHIA8_V
     Pythia8::ParticleDataEntry* V = fPythia->particleData.particleDataEntryPtr(r[i]);
#else
#if PYTHIA8_V<8309
     Pythia8::ParticleDataEntry* V = fPythia->particleData.particleDataEntryPtr(r[i]);
#else
     Pythia8::ParticleDataEntryPtr V = fPythia->particleData.particleDataEntryPtr(r[i]);
#endif
#endif     
     Pythia8::DecayChannel ch = V->channel(c[i]);
     if (TMath::Abs(ch.product(0))!=13 || TMath::Abs(ch.product(1))!=13){
      LOG(INFO) << "FixedTargetGenerator::Init this is not the right decay channel:"<<r[i] << " "<<c[i];
     }else{     
     TString tmp="";
     tmp+=r[i];tmp+=":";tmp+= c[i];
     tmp+=":bRatio =";
     tmp+=fBoost*ch.bRatio();
     fPythia->readString(tmp.Data());
     }
    }
   }
   fPythia->init();
  }
  // Initialize EvtGen.
  if (withEvtGen){
   TString SIMPATH=getenv("SIMPATH");
   TString DecayFile =SIMPATH+"/share/EvtGen/DECAY.DEC";
   TString ParticleFile =SIMPATH + "/share/EvtGen/evt.pdl";
   if(SIMPATH == "")
     {
       std::cout << "Using $EVTGENDATA "<< getenv("EVTGENDATA") << std::endl;
       DecayFile =TString(getenv("EVTGENDATA"))+"/DECAY.DEC";
       ParticleFile =TString(getenv("EVTGENDATA"))+ "/evt.pdl";
     }
   EvtAbsRadCorr *fsrPtrIn = 0;
   EvtExternalGenList *extPtr = new EvtExternalGenList();
   std::list<EvtDecayBase*> models = extPtr->getListOfModels();
 // Define the random number generator
   EvtRandomEngine* eng = new EvtSimpleRandomEngine();
   EvtRandom::setRandomEngine(eng);
   EvtGen *myEvtGenPtr = new EvtGen(DecayFile.Data(), ParticleFile.Data(),eng, fsrPtrIn, &models, 1, false);
   TString UdecayFile    = getenv("FAIRSHIP");UdecayFile +="/gconfig/USERDECAY.DEC";
   evtgenP = new EvtGenDecays(fPythiaP, DecayFile.Data(), ParticleFile.Data(),myEvtGenPtr);
   evtgenP->readDecayFile(UdecayFile.Data()); // will make update of EvtGen with user decay file
   // use one instance of EvtGen, requires patch to Pythia8Plugins/EvtGen.h
   if (Option == "Primary"){
    evtgenN = new EvtGenDecays(fPythiaN, DecayFile.Data(), ParticleFile.Data(),myEvtGenPtr);
   }
  }
  if (targetName!=""){
   fMaterialInvestigator = new GenieGenerator();
   if (!fcharmtarget){   
   TGeoNavigator* nav = gGeoManager->GetCurrentNavigator();
   nav->cd(targetName);
   TGeoNode* target = nav->GetCurrentNode(); 
   TObjArray* nodes =  target->GetVolume()->GetNodes();
   TGeoNode* first = (TGeoNode*)nodes->At(0);
   Int_t ilast = nodes->GetSize()-5; // assumes that the last 5 nodes are for the shielding around the target.
   TGeoNode* last  = (TGeoNode*)nodes->At(ilast);
   nav->cd(targetName+"/"+first->GetName()); 
   TGeoBBox* sha = (TGeoBBox*)first->GetVolume()->GetShape();
   Double_t dz   = sha->GetDZ();
   Double_t origin[3] = {0,0,-dz};
   Double_t master[3] = {0,0,0};
   nav->LocalToMaster(origin,master);
   startZ =  master[2];
   nav->cd(targetName+"/"+last->GetName()); 
   sha = (TGeoBBox*)first->GetVolume()->GetShape();
   dz   = sha->GetDZ();
   origin[2] = +dz;
   nav->LocalToMaster(origin,master);
   endZ =  master[2];
   start[0]=xOff;
   start[1]=yOff;
   start[2]=startZ;
   end[0]=xOff;
   end[1]=yOff;
   end[2]=endZ;
   }
   else{  //charm geometry uses a different target 
   TGeoVolume* top = gGeoManager->GetTopVolume();
   TGeoNode* target = top->FindNode(targetName);
   if (!target){
       LOG(ERROR) << "FixedTargetGenerator::Init target not found, "<< targetName.Data() <<", program will crash";
   }
   Double_t z_middle = target->GetMatrix()->GetTranslation()[2];
   TGeoBBox* sha = (TGeoBBox*)target->GetVolume()->GetShape();
   startZ =  z_middle - sha->GetDZ();
   endZ   =  z_middle + sha->GetDZ();
   start[0]=xOff;
   start[1]=yOff;
   start[2]=startZ;
   end[0]=xOff;
   end[1]=yOff;
   end[2]=endZ;
   }
//find maximum interaction length
   bparam = fMaterialInvestigator->MeanMaterialBudget(start, end, mparam);
   maxCrossSection =  mparam[9];
  }
 
  return kTRUE;
}

InitForCharmOrBeauty()

Bool_t FixedTargetGenerator::InitForCharmOrBeauty	(	TString	fInName,
		Int_t	nev,
		Double_t	npots = 5E13,
		Int_t	nStart = 0
	)

Definition at line 52 of file FixedTargetGenerator.cxx.

{
  Option = "charm";
  nEntry = nStart;
  // open input file with charm or beauty
  fin   = TFile::Open(fInName);
  nTree = (TNtuple*)fin->FindObjectAny("pythia6"); // old format, simple ntuple
  nEvents = nTree->GetEntries();
  nTree->SetBranchAddress("id",&n_id);
  nTree->SetBranchAddress("px",&n_px);
  nTree->SetBranchAddress("py",&n_py);
  nTree->SetBranchAddress("pz",&n_pz);
  nTree->SetBranchAddress("M",&n_M);
  nTree->SetBranchAddress("E",&n_E);
  nTree->SetBranchAddress("mid",&n_mid);
  nTree->SetBranchAddress("mpx",&n_mpx);
  nTree->SetBranchAddress("mpy",&n_mpy);
  nTree->SetBranchAddress("mpz",&n_mpz);
  nTree->SetBranchAddress("mE",&n_mE);
  if (nTree->GetBranch("k")){
   LOG(INFO) << "FixedTargetGenerator::InitForCharmOrBeauty: +++has branch+++ ";
   nTree->SetBranchAddress("k",&ck);}
// check if we deal with charm or beauty:
  nTree->GetEvent(0);
  if (!setByHand and n_M>5){ 
    chicc = chibb;
    LOG(INFO) << "FixedTargetGenerator::InitForCharmOrBeauty: automatic detection of beauty, configured for beauty";
    LOG(INFO) << "FixedTargetGenerator::InitForCharmOrBeauty: bb cross section / mbias"<<chicc;
  }else{
    LOG(INFO) << "FixedTargetGenerator::InitForCharmOrBeauty: cc cross section / mbias"<<chicc;
  }
// convert pot to weight corresponding to one spill of 5e13 pot
 // get histogram with number of pot to normalise
 // pot are counted double, i.e. for each signal, i.e. pot/2.
  Int_t nrcpot=((TH1F*)fin->Get("2"))->GetBinContent(1)/2.; // number of primary interactions
  wspill = nrpotspill*chicc/nrcpot*nEvents/nev;
  LOG(INFO) << "FixedTargetGenerator::InitForCharmOrBeauty: Input file:"<<fInName.Data()<<" with "<<nEvents<< " entries, corresponding to nr-pot="<<nrcpot/chicc;
  LOG(INFO) << "FixedTargetGenerator::InitForCharmOrBeauty: weight "<< wspill <<", corresponding to "<<nrpotspill<<" p.o.t. per spill for "<< nev <<" events to process";
 
  pot=0.;
  //Determine fDs on this file for primaries
  nDsprim=0;
  ntotprim=0;
 
  return kTRUE;
}

Print()

void FixedTargetGenerator::Print

(

)

ReadEvent()

Bool_t FixedTargetGenerator::ReadEvent

(

FairPrimaryGenerator *

cpg

)

public method ReadEvent

Definition at line 302 of file FixedTargetGenerator.cxx.

{
  Double_t zinter=0;
  Double_t ZoverA = 1.;
  if (targetName.Data() !=""){
// calculate primary proton interaction point:
// loop over trajectory between start and end to pick an interaction point, copied from GenieGenerator and adapted to hadrons
   Double_t prob2int = -1.;
   Double_t rndm = 0.;
   Double_t sigma;
   Int_t count=0;
   Double_t zinterStart = start[2];
   if (Option == "charm" || Option == "beauty"){ 
// simulate more downstream interaction points for interactions down in the cascade
    if (!(nTree->GetBranch("k"))){ck=1;}
   }
   else {ck=1;} 
   while (ck>0.5){
    while (prob2int<rndm) {
 //place x,y,z uniform along path
      zinter = gRandom->Uniform(zinterStart,end[2]);
      Double_t point[3]={xOff,yOff,zinter};
      bparam = fMaterialInvestigator->MeanMaterialBudget(start, point, mparam);
      Double_t interLength = mparam[8]; 
      TGeoNode *node = gGeoManager->FindNode(point[0],point[1],point[2]);
      TGeoMaterial *mat = 0;
      if (node && !gGeoManager->IsOutside()) {
         mat = node->GetVolume()->GetMaterial();
         Double_t n = mat->GetDensity()/mat->GetA();
         ZoverA =  mat->GetZ() / mat->GetA();
         sigma = 1./(n*mat->GetIntLen())/mbarn;
         prob2int = TMath::Exp(-interLength)*sigma/maxCrossSection;
      }else{
         prob2int=0.;
      }
      rndm = gRandom->Uniform(0.,1.); 
      count+=1;
    }
    zinterStart = zinter;
    ck-=1;
  } 
  zinter = zinter*cm;
  }
  Pythia8::Pythia* fPythia;
  if (G4only){
   cpg->AddTrack(2212,0.,0.,fMom,xOff/cm,yOff/cm,start[2],-1,kTRUE,-1.,0.,1.);
   return kTRUE;
  }else if (Option == "Primary"){
   if (gRandom->Uniform(0.,1.) < ZoverA ){
    fPythiaP->next();
    if (withEvtGen){evtgenP->decay();}
    fPythia = fPythiaP;
   }else{
    fPythiaN->next();
    if (withEvtGen){evtgenN->decay();}
    fPythia = fPythiaN;
   }
  }else{
    if (nEntry==nEvents){
      LOG(INFO) << "FixedTargetGenerator::readEvent Rewind input file: "<<nEntry;
      nEntry=0;}
    nTree->GetEvent(nEntry);
    nEntry+=1;
    // sanity check, count number of p.o.t. on input file.
    Double_t pt=TMath::Sqrt( (n_mpx*n_mpx)+(n_mpy*n_mpy));
    // every event appears twice, i.e.
    if (pt<1.e-5 && n_mid==2212){ 
     pot+=+0.5;
     ntotprim+=1;}
    Int_t idabs=int(TMath::Abs(n_id));
    if (idabs==431){ nDsprim+=1;}
    fPythiaP->event.reset();
    Int_t nID1 = n_id;
    fPythiaP->event.append(int(n_id),1,0,0,n_px,n_py,n_pz,n_E,n_M,0.,9.);
    TMCProcess procID  = kPTransportation;
    if (n_mid==2212 && (n_mpx*n_mpx+n_mpy*n_mpy)<1E-5) {procID = kPPrimary;} // probably primary and not from cascade
    cpg->AddTrack(int(n_mid),n_mpx,n_mpy,n_mpz, xOff/cm,yOff/cm,zinter/cm,-1,kFALSE,n_mE,0.,wspill,procID);
// second charm hadron in the event
    nTree->GetEvent(nEntry);
    if (nID1 * n_id > 0){LOG(INFO) << "FixedTargetGenerator::readEvent same sign charm: "<<nEntry<<" "<< nID1<<" "<< n_id;}
    nEntry+=1;
    fPythiaP->event.append(int(n_id),1,0,0,n_px,n_py,n_pz,n_E,n_M,0.,9.);
    fPythiaP->next();
    fPythia = fPythiaP;
  }
  if (withEvtGen){
   fPythia->moreDecays();} // let the very short lived resonances decay via Pythia8
  if(Debug && !G4only){fPythia->event.list();}
  TMCProcess procID; 
  for(Int_t ii=1; ii<fPythia->event.size(); ii++){
     Double_t  e = fPythia->event[ii].e();
     Double_t  m = fPythia->event[ii].m();
     Double_t pz = fPythia->event[ii].pz();
     Int_t id = fPythia->event[ii].id(); 
     Bool_t wanttracking=kTRUE;
     if (e-m<EMax || !fPythia->event[ii].isFinal() || pz<0) {wanttracking=kFALSE;}
     if (DrellYan || PhotonCollision || OnlyMuons){
// don't track underlying event
      if (fabs(id)!=13){wanttracking=kFALSE;}
     }
     Double_t z  = fPythia->event[ii].zProd()+zinter;
     Double_t x  = fPythia->event[ii].xProd()+xOff;
     Double_t y  = fPythia->event[ii].yProd()+yOff;
     Double_t tof = fPythia->event[ii].tProd() / (10*c_light) ; // to go from mm to s
     Double_t px = fPythia->event[ii].px();
     Double_t py = fPythia->event[ii].py();
     Int_t im = fPythia->event[ii].mother1()-1;
     procID = kPPrimary;
     if (Option != "Primary"){
        procID = kPDecay;
        im+=1;
        if (ii==1) {procID  = kPHadronic;}
     }else{
      if (ii<3){im=-1;}
     }
     cpg->AddTrack(id,px,py,pz,x/cm,y/cm,z/cm,im,wanttracking,e,tof,wspill,procID);
     if(withNtuple){
          int idabs = TMath::Abs(id);
          if (idabs<10)  {continue;}
          if (idabs<18 || idabs==22 || idabs==111 || idabs==221 || idabs==223 || idabs==331 
                 || idabs==211 || idabs==113 || idabs==333  || idabs==321   || idabs==2212 ){
          Int_t moID = fPythia->event[im+1].id();
          fNtuple->Fill(0,id,px,py,pz,e,x/cm,y/cm,z/cm,moID);
         }
     }
    }
  return kTRUE;
}

SetBoost()

void FixedTargetGenerator::SetBoost ( Double_t  f )

inline

Definition at line 38 of file FixedTargetGenerator.h.

{ fBoost  = f; }  // boost factor for rare di-muon decays

SetCharmTarget()

void FixedTargetGenerator::SetCharmTarget ( Bool_t  charmtarget = true )

inline

Definition at line 36 of file FixedTargetGenerator.h.

{fcharmtarget = charmtarget;}; //charm geometry uses a different target, default one is usual

SetChibb()

void FixedTargetGenerator::SetChibb ( Double_t  x )

inline

Definition at line 46 of file FixedTargetGenerator.h.

{ chibb = x; }  // chibb = bbbar over mbias cross section   

SetChicc()

void FixedTargetGenerator::SetChicc ( Double_t  x )

inline

Definition at line 47 of file FixedTargetGenerator.h.

{ chicc = x; }  // chicc = ccbar over mbias cross section   

SetDebug()

void FixedTargetGenerator::SetDebug ( Bool_t  x )

inline

Definition at line 51 of file FixedTargetGenerator.h.

{Debug=x;}

SetDrellYan()

void FixedTargetGenerator::SetDrellYan ( )

inline

Definition at line 43 of file FixedTargetGenerator.h.

{ DrellYan  = true; }  // only generate prompt Z0* processes

SetEnergyCut()

void FixedTargetGenerator::SetEnergyCut ( Float_t  emax )

inline

Definition at line 50 of file FixedTargetGenerator.h.

{EMax=emax;}// min energy to be copied to Geant4

SetG4only()

void FixedTargetGenerator::SetG4only ( )

inline

Definition at line 39 of file FixedTargetGenerator.h.

{ G4only  = true; }  // only run Geant4, no pythia primary interaction

SetHeartBeat()

void FixedTargetGenerator::SetHeartBeat ( Int_t  x )

inline

Definition at line 49 of file FixedTargetGenerator.h.

{heartbeat=x;}

SetJpsiMainly()

void FixedTargetGenerator::SetJpsiMainly ( )

inline

Definition at line 41 of file FixedTargetGenerator.h.

{ JpsiMainly  = true; }  // let all Jpsi decay to mumu

SetMom()

void FixedTargetGenerator::SetMom ( Double_t  mom )

inline

Definition at line 33 of file FixedTargetGenerator.h.

{ fMom = mom; };

SetOnlyMuons()

void FixedTargetGenerator::SetOnlyMuons ( )

inline

Definition at line 42 of file FixedTargetGenerator.h.

{ OnlyMuons = true; }  // only transport muons

SetOpt4DP()

void FixedTargetGenerator::SetOpt4DP ( TNtuple *  t )

inline

Definition at line 52 of file FixedTargetGenerator.h.

{withNtuple=kTRUE;  fNtuple = t ; }

SetParameters()

void FixedTargetGenerator::SetParameters

(

char *

)

SetPhotonCollision()

void FixedTargetGenerator::SetPhotonCollision ( )

inline

Definition at line 44 of file FixedTargetGenerator.h.

{ PhotonCollision   = true; }  // only generate prompt photon processes

SetSeed()

void FixedTargetGenerator::SetSeed ( Double_t  seed )

inline

Definition at line 48 of file FixedTargetGenerator.h.

{fSeed=seed;}

SetTarget()

void FixedTargetGenerator::SetTarget ( TString  s, Double_t  x, Double_t  y  )

inline

Definition at line 37 of file FixedTargetGenerator.h.

{ targetName = s; xOff=x; yOff=y; };

SetTauOnly()

void FixedTargetGenerator::SetTauOnly ( )

inline

Definition at line 40 of file FixedTargetGenerator.h.

{ tauOnly  = true; }  // only have Ds decay to tau 

UseRandom1()

void FixedTargetGenerator::UseRandom1 ( )

inline

Definition at line 34 of file FixedTargetGenerator.h.

{ fUseRandom1 = kTRUE; fUseRandom3 = kFALSE; };

UseRandom3()

void FixedTargetGenerator::UseRandom3 ( )

inline

Definition at line 35 of file FixedTargetGenerator.h.

{ fUseRandom1 = kFALSE; fUseRandom3 = kTRUE; };

WithEvtGen()

void FixedTargetGenerator::WithEvtGen ( )

inline

Definition at line 45 of file FixedTargetGenerator.h.

{ withEvtGen = true;} // use EvtGen as external decayer to Pythia, experimental phase, only works for one Pythia instance

Member Data Documentation

bparam

Double_t FixedTargetGenerator::bparam

protected

Definition at line 82 of file FixedTargetGenerator.h.

chibb

Double_t FixedTargetGenerator::chibb

protected

Definition at line 65 of file FixedTargetGenerator.h.

chicc

Double_t FixedTargetGenerator::chicc

protected

Definition at line 65 of file FixedTargetGenerator.h.

ck

Float_t FixedTargetGenerator::ck

protected

Definition at line 89 of file FixedTargetGenerator.h.

Debug

Bool_t FixedTargetGenerator::Debug

protected

Definition at line 67 of file FixedTargetGenerator.h.

DrellYan

Bool_t FixedTargetGenerator::DrellYan

protected

Definition at line 67 of file FixedTargetGenerator.h.

EMax

Double_t FixedTargetGenerator::EMax

protected

Definition at line 65 of file FixedTargetGenerator.h.

end

Double_t FixedTargetGenerator::end[3]

protected

Definition at line 81 of file FixedTargetGenerator.h.

endZ

Double_t FixedTargetGenerator::endZ

protected

Definition at line 85 of file FixedTargetGenerator.h.

evtgenN

EvtGenDecays* FixedTargetGenerator::evtgenN

protected

Definition at line 72 of file FixedTargetGenerator.h.

evtgenP

EvtGenDecays* FixedTargetGenerator::evtgenP

protected

Definition at line 73 of file FixedTargetGenerator.h.

fBoost

Double_t FixedTargetGenerator::fBoost

protected

Definition at line 65 of file FixedTargetGenerator.h.

fcharmtarget

Bool_t FixedTargetGenerator::fcharmtarget

protected

Definition at line 68 of file FixedTargetGenerator.h.

fin

TFile* FixedTargetGenerator::fin

protected

Definition at line 87 of file FixedTargetGenerator.h.

fLogger

FairLogger* FixedTargetGenerator::fLogger

protected

Definition at line 69 of file FixedTargetGenerator.h.

fMaterialInvestigator

GenieGenerator* FixedTargetGenerator::fMaterialInvestigator

protected

Definition at line 74 of file FixedTargetGenerator.h.

fMom

Double_t FixedTargetGenerator::fMom

protected

Definition at line 62 of file FixedTargetGenerator.h.

fNtuple

TNtuple* FixedTargetGenerator::fNtuple

protected

special option for Dark Photon physics studies

Definition at line 76 of file FixedTargetGenerator.h.

fPythiaN

Pythia8::Pythia* FixedTargetGenerator::fPythiaN

protected

don't make it persistent, magic ROOT command

Definition at line 70 of file FixedTargetGenerator.h.

fPythiaP

Pythia8::Pythia* FixedTargetGenerator::fPythiaP

protected

Definition at line 71 of file FixedTargetGenerator.h.

fRandomEngine

Pythia8::RndmEngine* FixedTargetGenerator::fRandomEngine

private

Definition at line 58 of file FixedTargetGenerator.h.

fSeed

Double_t FixedTargetGenerator::fSeed

protected

Definition at line 65 of file FixedTargetGenerator.h.

fUseRandom1

Bool_t FixedTargetGenerator::fUseRandom1

protected

Definition at line 63 of file FixedTargetGenerator.h.

fUseRandom3

Bool_t FixedTargetGenerator::fUseRandom3

protected

Definition at line 64 of file FixedTargetGenerator.h.

G4only

Bool_t FixedTargetGenerator::G4only

protected

Definition at line 67 of file FixedTargetGenerator.h.

heartbeat

Int_t FixedTargetGenerator::heartbeat

protected

Definition at line 90 of file FixedTargetGenerator.h.

JpsiMainly

Bool_t FixedTargetGenerator::JpsiMainly

protected

Definition at line 67 of file FixedTargetGenerator.h.

maxCrossSection

Double_t FixedTargetGenerator::maxCrossSection

protected

Definition at line 86 of file FixedTargetGenerator.h.

mparam

Double_t FixedTargetGenerator::mparam[10]

protected

Definition at line 83 of file FixedTargetGenerator.h.

n_E

Float_t FixedTargetGenerator::n_E

protected

Definition at line 89 of file FixedTargetGenerator.h.

n_id

Float_t FixedTargetGenerator::n_id

protected

Definition at line 89 of file FixedTargetGenerator.h.

n_M

Float_t FixedTargetGenerator::n_M

protected

Definition at line 89 of file FixedTargetGenerator.h.

n_mE

Float_t FixedTargetGenerator::n_mE

protected

Definition at line 89 of file FixedTargetGenerator.h.

n_mid

Float_t FixedTargetGenerator::n_mid

protected

Definition at line 89 of file FixedTargetGenerator.h.

n_mpx

Float_t FixedTargetGenerator::n_mpx

protected

Definition at line 89 of file FixedTargetGenerator.h.

n_mpy

Float_t FixedTargetGenerator::n_mpy

protected

Definition at line 89 of file FixedTargetGenerator.h.

n_mpz

Float_t FixedTargetGenerator::n_mpz

protected

Definition at line 89 of file FixedTargetGenerator.h.

n_px

Float_t FixedTargetGenerator::n_px

protected

Definition at line 89 of file FixedTargetGenerator.h.

n_py

Float_t FixedTargetGenerator::n_py

protected

Definition at line 89 of file FixedTargetGenerator.h.

n_pz

Float_t FixedTargetGenerator::n_pz

protected

Definition at line 89 of file FixedTargetGenerator.h.

nDsprim

Int_t FixedTargetGenerator::nDsprim

protected

Definition at line 66 of file FixedTargetGenerator.h.

nEntry

Int_t FixedTargetGenerator::nEntry

protected

Definition at line 66 of file FixedTargetGenerator.h.

nEvents

Int_t FixedTargetGenerator::nEvents

protected

Definition at line 66 of file FixedTargetGenerator.h.

nrpotspill

Double_t FixedTargetGenerator::nrpotspill

protected

Definition at line 65 of file FixedTargetGenerator.h.

ntotprim

Int_t FixedTargetGenerator::ntotprim

protected

Definition at line 66 of file FixedTargetGenerator.h.

nTree

TNtuple* FixedTargetGenerator::nTree

protected

Definition at line 88 of file FixedTargetGenerator.h.

OnlyMuons

Bool_t FixedTargetGenerator::OnlyMuons

protected

Definition at line 67 of file FixedTargetGenerator.h.

Option

TString FixedTargetGenerator::Option

protected

Definition at line 77 of file FixedTargetGenerator.h.

PhotonCollision

Bool_t FixedTargetGenerator::PhotonCollision

protected

Definition at line 67 of file FixedTargetGenerator.h.

pot

Int_t FixedTargetGenerator::pot

protected

Definition at line 66 of file FixedTargetGenerator.h.

setByHand

Bool_t FixedTargetGenerator::setByHand

protected

Definition at line 67 of file FixedTargetGenerator.h.

start

Double_t FixedTargetGenerator::start[3]

protected

Definition at line 80 of file FixedTargetGenerator.h.

startZ

Double_t FixedTargetGenerator::startZ

protected

Definition at line 84 of file FixedTargetGenerator.h.

targetName

TString FixedTargetGenerator::targetName

protected

Definition at line 77 of file FixedTargetGenerator.h.

tauOnly

Bool_t FixedTargetGenerator::tauOnly

protected

Definition at line 67 of file FixedTargetGenerator.h.

withEvtGen

Bool_t FixedTargetGenerator::withEvtGen

protected

Definition at line 67 of file FixedTargetGenerator.h.

withNtuple

Bool_t FixedTargetGenerator::withNtuple

protected

Definition at line 75 of file FixedTargetGenerator.h.

wspill

Double_t FixedTargetGenerator::wspill

protected

Definition at line 65 of file FixedTargetGenerator.h.

xOff

Double_t FixedTargetGenerator::xOff

protected

Definition at line 78 of file FixedTargetGenerator.h.

yOff

Double_t FixedTargetGenerator::yOff

protected

Definition at line 79 of file FixedTargetGenerator.h.

---

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

• shipgen/FixedTargetGenerator.h
• shipgen/FixedTargetGenerator.cxx