FixedTargetGenerator Class Reference

#include <FixedTargetGenerator.h>

Inheritance diagram for FixedTargetGenerator:

Collaboration diagram for FixedTargetGenerator:

Public Member Functions
	FixedTargetGenerator ()
	~FixedTargetGenerator () override
Bool_t	ReadEvent (FairPrimaryGenerator *) override
void	SetParameters (char *)
void	Print ()
Bool_t	Init (const char *inFile) override
Bool_t	Init (const char *inFile, int startEvent) override
Bool_t	Init () override
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	SetTarget (TString s, Double_t x, Double_t y)
void	SetZoffset (Double_t z)
void	SetXoffset (Double_t x)
void	SetYoffset (Double_t y)
void	SetSmearBeam (Double_t sb)
void	SetPaintRadius (Double_t r)
void	SetTargetCoordinates (Double_t start_z, Double_t end_z)
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 ()
Public Member Functions inherited from SHiP::Generator
	Generator ()
virtual	~Generator ()
virtual Bool_t	Init (const char *, int)=0
virtual Bool_t	Init (const char *)=0
virtual Bool_t	Init (const std::vector< std::string > &inFiles, int startNumber)
virtual Bool_t	Init (const std::vector< std::string > &inFiles)
virtual void	UseExternalFile (std::string x, Int_t i)
virtual void	UseExternalFile (std::vector< std::string > &inFiles, Int_t i)

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
FairLogger *	fLogger
Pythia8::Pythia *	fPythiaN
	don't make it persistent, magic ROOT command
Pythia8::Pythia *	fPythiaP
Pythia8::EvtGenDecays *	evtgenN
Pythia8::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	zOff
Double_t	fsmearBeam
Double_t	fPaintBeam
Double_t	start [3]
Double_t	end [3]
Double_t	bparam
Double_t	mparam [10]
Double_t	startZ
Double_t	endZ
Bool_t	targetFromGeometry
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
Protected Attributes inherited from SHiP::Generator
std::optional< std::string >	fextFile
Int_t	firstEvent = 0

Private Attributes
std::shared_ptr< Pythia8::RndmEngine >	fRandomEngine

Detailed Description

Definition at line 22 of file FixedTargetGenerator.h.

Constructor & Destructor Documentation

FixedTargetGenerator()

FixedTargetGenerator::FixedTargetGenerator

(

)

default constructor

Definition at line 33 of file FixedTargetGenerator.cxx.

                                           {
  fUseRandom1 = kFALSE;
  fUseRandom3 = kTRUE;
  fSeed = 0;
  fMom = 400;  // proton
  fLogger = FairLogger::GetLogger();
  targetName = "cave_1/target_vacuum_box_1/TargetArea_1/HeVolume_1";
  xOff = 0;
  yOff = 0;
  zOff = 0;
  fsmearBeam = 8 * mm;  // default value for smearing beam (8 mm)
  fPaintBeam = 5 * cm;  // default value for painting beam (5 cm)
  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;
  targetFromGeometry = kFALSE;
  Option = "Primary";
  wspill = 1.;  // event weight == 1 for primary events
  heartbeat = 1000;
}

~FixedTargetGenerator()

FixedTargetGenerator::~FixedTargetGenerator ( )

override

destructor

Definition at line 325 of file FixedTargetGenerator.cxx.

{}

Member Function Documentation

GetPotForCharm()

Double_t FixedTargetGenerator::GetPotForCharm ( )

inline

Definition at line 104 of file FixedTargetGenerator.h.

{ return nrpotspill / wspill; }

GetPythia()

Pythia8::Pythia * FixedTargetGenerator::GetPythia ( )

inline

Definition at line 105 of file FixedTargetGenerator.h.

{ return fPythiaP; }

GetPythiaN()

Pythia8::Pythia * FixedTargetGenerator::GetPythiaN ( )

inline

Definition at line 106 of file FixedTargetGenerator.h.

{ return fPythiaN; }

Init() [1/3]

Bool_t FixedTargetGenerator::Init ( )

override

Definition at line 127 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(error) << "Option not known " << Option.Data() << ", abort";
  }
  if (fUseRandom1) fRandomEngine = std::make_shared<PyTr1Rng>();
  if (fUseRandom3) fRandomEngine = std::make_shared<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);  // codespell:ignore parm
      fPythia->settings.parm("Beams:eB", 0.);    // codespell:ignore parm
    }
    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);
      std::shared_ptr<Pythia8::ParticleDataEntry> p =
          fPythia->particleData.particleDataEntryPtr(n);
      if (p->tau0() > 1) {
        std::string particle = std::to_string(n) + ":mayDecay = false";
        fPythia->readString(particle);
        LOG(info) << "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) << "Rescale BRs of dimuon decays in Pythia: " << fBoost;
      for (unsigned int i = 0; i < r.size(); ++i) {
        std::shared_ptr<Pythia8::ParticleDataEntry> V =
            fPythia->particleData.particleDataEntryPtr(r[i]);
        Pythia8::DecayChannel ch = V->channel(c[i]);
        if (TMath::Abs(ch.product(0)) != 13 ||
            TMath::Abs(ch.product(1)) != 13) {
          LOG(info) << "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) {
    // Use EVTGENDATA environment variable to find EvtGen data files
    const char* evtgendata = getenv("EVTGENDATA");
    if (!evtgendata) {
      LOG(fatal) << "EVTGENDATA environment variable not set";
    }
    TString DecayFile = TString(evtgendata) + "/DECAY.DEC";
 
    TString ParticleFile = TString(evtgendata) + "/evt.pdl";
    std::cout << "Using $EVTGENDATA " << evtgendata << std::endl;
    EvtExternalGenList* extPtr = new EvtExternalGenList();
    std::list<EvtDecayBase*> models = extPtr->getListOfModels();
    TString UdecayFile = getenv("FAIRSHIP");
    UdecayFile += "/gconfig/USERDECAY.DEC";
    evtgenP = new Pythia8::EvtGenDecays(fPythiaP, DecayFile.Data(),
                                        ParticleFile.Data(), extPtr);
    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 Pythia8::EvtGenDecays(fPythiaN, DecayFile.Data(), "", extPtr);
    }
  }
  if (targetFromGeometry) {
    // Use geometry-based coordinates passed from run_simScript.py
    fMaterialInvestigator = new GenieGenerator();
    start[0] = xOff;
    start[1] = yOff;
    start[2] = startZ + zOff;
    end[0] = xOff;
    end[1] = yOff;
    end[2] = endZ;
    LOG(info) << "FixedTargetGenerator: Using geometry-based target "
                 "coordinates startZ="
              << startZ << " endZ=" << endZ;
    // find maximum interaction length
    bparam = shipgen::MeanMaterialBudget(start, end, mparam);
    maxCrossSection = mparam[9];
  } else if (targetName != "") {
    // Fallback to fragile TGeo navigation for backward compatibility
    fMaterialInvestigator = new GenieGenerator();
    TGeoNavigator* nav = gGeoManager->GetCurrentNavigator();
    if (nav->CheckPath(targetName)) {
      nav->cd(targetName);
    } else {
      LOG(fatal) << "Invalid target volume specified";
    }
    TGeoNode* target = nav->GetCurrentNode();
    TObjArray* nodes = target->GetVolume()->GetNodes();
    // Get the first and last node of the target to calculate the material seen
    TGeoNode* first = static_cast<TGeoNode*>(nodes->At(0));
    TGeoNode* last = static_cast<TGeoNode*>(nodes->At(nodes->GetSize() - 1));
    nav->cd(targetName + "/" + first->GetName());
    TGeoBBox* sha = static_cast<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 = static_cast<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 + zOff;
    end[0] = xOff;
    end[1] = yOff;
    end[2] = endZ;
    // find maximum interaction length
    bparam = shipgen::MeanMaterialBudget(start, end, mparam);
    maxCrossSection = mparam[9];
  } else {
    LOG(fatal) << "No target set.";
  }
 
  return kTRUE;
}

Init() [2/3]

Bool_t FixedTargetGenerator::Init ( const char *  inFile )

inlineoverridevirtual

Implements SHiP::Generator.

Definition at line 34 of file FixedTargetGenerator.h.

{ return Init(inFile, 0); };

Init() [3/3]

Bool_t FixedTargetGenerator::Init ( const char *  inFile, int  startEvent  )

inlineoverridevirtual

Implements SHiP::Generator.

Definition at line 36 of file FixedTargetGenerator.h.

                                                           {
    LOG(warning) << "Init with files not implemented for FixedTargetGenerator. "
                    "Using default Init() instead";
    return Init();
  };

InitForCharmOrBeauty()

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

Definition at line 66 of file FixedTargetGenerator.cxx.

                                                                               {
  Option = "charm";
  nEntry = nStart;
  // open input file with charm or beauty
  fin = TFile::Open(fInName);
  if (!fin) {
    LOG(fatal) << "Could not open input file: " << fInName.Data();
    return kFALSE;
  }
  nTree = dynamic_cast<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) << "+++has branch+++";
    nTree->SetBranchAddress("k", &ck);
  }
  // check if we deal with charm or beauty:
  nTree->GetEvent(0);
  if (!setByHand && n_M > 5) {
    chicc = chibb;
    LOG(info) << "automatic detection of beauty, configured for beauty";
    LOG(info) << "bb cross section / mbias " << chicc;
  } else {
    LOG(info) << "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.
  auto* potHist = dynamic_cast<TH1F*>(fin->Get("2"));
  if (!potHist) {
    LOG(fatal) << "Histogram '2' not found in input file";
    return kFALSE;
  }
  Int_t nrcpot =
      potHist->GetBinContent(1) / 2.;  // number of primary interactions
  wspill = nrpotspill * chicc / nrcpot * nEvents / nev;
  LOG(info) << "Input file: " << fInName.Data() << " with " << nEvents
            << " entries, corresponding to nr-pot=" << (nrcpot / chicc);
  LOG(info) << "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 )

override

public method ReadEvent

Definition at line 329 of file FixedTargetGenerator.cxx.

                                                                {
  // Calculate beam smearing and painting
  auto [dx, dy] = CalculateBeamOffset(fsmearBeam, fPaintBeam);
 
  Double_t zinter = 0;
  Double_t ZoverA = 1.;
  if (!targetName.IsNull()) {
    // 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;
    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 = shipgen::MeanMaterialBudget(start, point, mparam);
        Double_t interLength = mparam[8];
        TGeoNode* node = gGeoManager->FindNode(point[0], point[1], point[2]);
        TGeoMaterial* mat = nullptr;
        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.);
      }
      zinterStart = zinter;
      ck -= 1;
    }
    zinter = zinter * cm;
  }
  Pythia8::Pythia* fPythia;
  if (G4only) {
    cpg->AddTrack(2212, 0., 0., fMom, (xOff + dx) / cm, (yOff + dy) / 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) << "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 = static_cast<int>(TMath::Abs(n_id));
    if (idabs == 431) {
      nDsprim += 1;
    }
    fPythiaP->event.reset();
    Int_t nID1 = n_id;
    fPythiaP->event.append(static_cast<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(static_cast<int>(n_mid), n_mpx, n_mpy, n_mpz,
                  (xOff + dx) * cm, (yOff + dy) * 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) << "same sign charm: " << nEntry << ", " << nID1 << ", "
                << n_id;
    }
    nEntry += 1;
    fPythiaP->event.append(static_cast<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() * mm + zinter * cm;
    Double_t x = fPythia->event[ii].xProd() * mm + xOff * cm + dx * cm;
    Double_t y = fPythia->event[ii].yProd() * mm + yOff * cm + dy * cm;
    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 71 of file FixedTargetGenerator.h.

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

SetChibb()

void FixedTargetGenerator::SetChibb ( Double_t  x )

inline

Definition at line 88 of file FixedTargetGenerator.h.

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

SetChicc()

void FixedTargetGenerator::SetChicc ( Double_t  x )

inline

Definition at line 91 of file FixedTargetGenerator.h.

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

SetDebug()

void FixedTargetGenerator::SetDebug ( Bool_t  x )

inline

Definition at line 99 of file FixedTargetGenerator.h.

{ Debug = x; }

SetDrellYan()

void FixedTargetGenerator::SetDrellYan ( )

inline

Definition at line 80 of file FixedTargetGenerator.h.

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

SetEnergyCut()

void FixedTargetGenerator::SetEnergyCut ( Float_t  emax )

inline

Definition at line 96 of file FixedTargetGenerator.h.

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

SetG4only()

void FixedTargetGenerator::SetG4only ( )

inline

Definition at line 74 of file FixedTargetGenerator.h.

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

SetHeartBeat()

void FixedTargetGenerator::SetHeartBeat ( Int_t  x )

inline

Definition at line 95 of file FixedTargetGenerator.h.

{ heartbeat = x; }

SetJpsiMainly()

void FixedTargetGenerator::SetJpsiMainly ( )

inline

Definition at line 78 of file FixedTargetGenerator.h.

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

SetMom()

void FixedTargetGenerator::SetMom ( Double_t  mom )

inline

Definition at line 47 of file FixedTargetGenerator.h.

{ fMom = mom; };

SetOnlyMuons()

void FixedTargetGenerator::SetOnlyMuons ( )

inline

Definition at line 79 of file FixedTargetGenerator.h.

{ OnlyMuons = true; }    // only transport muons

SetOpt4DP()

void FixedTargetGenerator::SetOpt4DP ( TNtuple *  t )

inline

Definition at line 100 of file FixedTargetGenerator.h.

                                    {
    withNtuple = kTRUE;
    fNtuple = t;
  }

SetPaintRadius()

void FixedTargetGenerator::SetPaintRadius ( Double_t  r )

inline

Definition at line 65 of file FixedTargetGenerator.h.

{ fPaintBeam = r; };

SetParameters()

void FixedTargetGenerator::SetParameters

(

char *

)

SetPhotonCollision()

void FixedTargetGenerator::SetPhotonCollision ( )

inline

Definition at line 81 of file FixedTargetGenerator.h.

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

SetSeed()

void FixedTargetGenerator::SetSeed ( Double_t  seed )

inline

Definition at line 94 of file FixedTargetGenerator.h.

{ fSeed = seed; }

SetSmearBeam()

void FixedTargetGenerator::SetSmearBeam ( Double_t  sb )

inline

Definition at line 64 of file FixedTargetGenerator.h.

{ fsmearBeam = sb; };

SetTarget()

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

inline

Definition at line 56 of file FixedTargetGenerator.h.

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

SetTargetCoordinates()

void FixedTargetGenerator::SetTargetCoordinates ( Double_t  start_z, Double_t  end_z  )

inline

Definition at line 66 of file FixedTargetGenerator.h.

                                                              {
    startZ = start_z;
    endZ = end_z;
    targetFromGeometry = true;
  };

SetTauOnly()

void FixedTargetGenerator::SetTauOnly ( )

inline

Definition at line 77 of file FixedTargetGenerator.h.

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

SetXoffset()

void FixedTargetGenerator::SetXoffset ( Double_t  x )

inline

Definition at line 62 of file FixedTargetGenerator.h.

{ xOff = x; };

SetYoffset()

void FixedTargetGenerator::SetYoffset ( Double_t  y )

inline

Definition at line 63 of file FixedTargetGenerator.h.

{ yOff = y; };

SetZoffset()

void FixedTargetGenerator::SetZoffset ( Double_t  z )

inline

Definition at line 61 of file FixedTargetGenerator.h.

{ zOff = z; };

UseRandom1()

void FixedTargetGenerator::UseRandom1 ( )

inline

Definition at line 48 of file FixedTargetGenerator.h.

                    {
    fUseRandom1 = kTRUE;
    fUseRandom3 = kFALSE;
  };

UseRandom3()

void FixedTargetGenerator::UseRandom3 ( )

inline

Definition at line 52 of file FixedTargetGenerator.h.

                    {
    fUseRandom1 = kFALSE;
    fUseRandom3 = kTRUE;
  };

WithEvtGen()

void FixedTargetGenerator::WithEvtGen ( )

inline

Definition at line 84 of file FixedTargetGenerator.h.

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

Member Data Documentation

bparam

Double_t FixedTargetGenerator::bparam

protected

Definition at line 135 of file FixedTargetGenerator.h.

chibb

Double_t FixedTargetGenerator::chibb

protected

Definition at line 115 of file FixedTargetGenerator.h.

chicc

Double_t FixedTargetGenerator::chicc

protected

Definition at line 115 of file FixedTargetGenerator.h.

ck

Float_t FixedTargetGenerator::ck

protected

Definition at line 144 of file FixedTargetGenerator.h.

Debug

Bool_t FixedTargetGenerator::Debug

protected

Definition at line 118 of file FixedTargetGenerator.h.

DrellYan

Bool_t FixedTargetGenerator::DrellYan

protected

Definition at line 117 of file FixedTargetGenerator.h.

EMax

Double_t FixedTargetGenerator::EMax

protected

Definition at line 115 of file FixedTargetGenerator.h.

end

Double_t FixedTargetGenerator::end[3]

protected

Definition at line 134 of file FixedTargetGenerator.h.

endZ

Double_t FixedTargetGenerator::endZ

protected

Definition at line 138 of file FixedTargetGenerator.h.

evtgenN

Pythia8::EvtGenDecays* FixedTargetGenerator::evtgenN

protected

Definition at line 122 of file FixedTargetGenerator.h.

evtgenP

Pythia8::EvtGenDecays* FixedTargetGenerator::evtgenP

protected

Definition at line 123 of file FixedTargetGenerator.h.

fBoost

Double_t FixedTargetGenerator::fBoost

protected

Definition at line 115 of file FixedTargetGenerator.h.

fin

TFile* FixedTargetGenerator::fin

protected

Definition at line 141 of file FixedTargetGenerator.h.

fLogger

FairLogger* FixedTargetGenerator::fLogger

protected

Definition at line 119 of file FixedTargetGenerator.h.

fMaterialInvestigator

GenieGenerator* FixedTargetGenerator::fMaterialInvestigator

protected

Definition at line 124 of file FixedTargetGenerator.h.

fMom

Double_t FixedTargetGenerator::fMom

protected

Definition at line 112 of file FixedTargetGenerator.h.

fNtuple

TNtuple* FixedTargetGenerator::fNtuple

protected

special option for Dark Photon physics studies

Definition at line 126 of file FixedTargetGenerator.h.

fPaintBeam

Double_t FixedTargetGenerator::fPaintBeam

protected

Definition at line 132 of file FixedTargetGenerator.h.

fPythiaN

Pythia8::Pythia* FixedTargetGenerator::fPythiaN

protected

don't make it persistent, magic ROOT command

Definition at line 120 of file FixedTargetGenerator.h.

fPythiaP

Pythia8::Pythia* FixedTargetGenerator::fPythiaP

protected

Definition at line 121 of file FixedTargetGenerator.h.

fRandomEngine

std::shared_ptr<Pythia8::RndmEngine> FixedTargetGenerator::fRandomEngine

private

Definition at line 109 of file FixedTargetGenerator.h.

fSeed

Double_t FixedTargetGenerator::fSeed

protected

Definition at line 115 of file FixedTargetGenerator.h.

fsmearBeam

Double_t FixedTargetGenerator::fsmearBeam

protected

Definition at line 131 of file FixedTargetGenerator.h.

fUseRandom1

Bool_t FixedTargetGenerator::fUseRandom1

protected

Definition at line 113 of file FixedTargetGenerator.h.

fUseRandom3

Bool_t FixedTargetGenerator::fUseRandom3

protected

Definition at line 114 of file FixedTargetGenerator.h.

G4only

Bool_t FixedTargetGenerator::G4only

protected

Definition at line 117 of file FixedTargetGenerator.h.

heartbeat

Int_t FixedTargetGenerator::heartbeat

protected

Definition at line 145 of file FixedTargetGenerator.h.

JpsiMainly

Bool_t FixedTargetGenerator::JpsiMainly

protected

Definition at line 117 of file FixedTargetGenerator.h.

maxCrossSection

Double_t FixedTargetGenerator::maxCrossSection

protected

Definition at line 140 of file FixedTargetGenerator.h.

mparam

Double_t FixedTargetGenerator::mparam[10]

protected

Definition at line 136 of file FixedTargetGenerator.h.

n_E

Float_t FixedTargetGenerator::n_E

protected

Definition at line 143 of file FixedTargetGenerator.h.

n_id

Float_t FixedTargetGenerator::n_id

protected

Definition at line 143 of file FixedTargetGenerator.h.

n_M

Float_t FixedTargetGenerator::n_M

protected

Definition at line 143 of file FixedTargetGenerator.h.

n_mE

Float_t FixedTargetGenerator::n_mE

protected

Definition at line 143 of file FixedTargetGenerator.h.

n_mid

Float_t FixedTargetGenerator::n_mid

protected

Definition at line 143 of file FixedTargetGenerator.h.

n_mpx

Float_t FixedTargetGenerator::n_mpx

protected

Definition at line 143 of file FixedTargetGenerator.h.

n_mpy

Float_t FixedTargetGenerator::n_mpy

protected

Definition at line 143 of file FixedTargetGenerator.h.

n_mpz

Float_t FixedTargetGenerator::n_mpz

protected

Definition at line 143 of file FixedTargetGenerator.h.

n_px

Float_t FixedTargetGenerator::n_px

protected

Definition at line 143 of file FixedTargetGenerator.h.

n_py

Float_t FixedTargetGenerator::n_py

protected

Definition at line 143 of file FixedTargetGenerator.h.

n_pz

Float_t FixedTargetGenerator::n_pz

protected

Definition at line 143 of file FixedTargetGenerator.h.

nDsprim

Int_t FixedTargetGenerator::nDsprim

protected

Definition at line 116 of file FixedTargetGenerator.h.

nEntry

Int_t FixedTargetGenerator::nEntry

protected

Definition at line 116 of file FixedTargetGenerator.h.

nEvents

Int_t FixedTargetGenerator::nEvents

protected

Definition at line 116 of file FixedTargetGenerator.h.

nrpotspill

Double_t FixedTargetGenerator::nrpotspill

protected

Definition at line 115 of file FixedTargetGenerator.h.

ntotprim

Int_t FixedTargetGenerator::ntotprim

protected

Definition at line 116 of file FixedTargetGenerator.h.

nTree

TNtuple* FixedTargetGenerator::nTree

protected

Definition at line 142 of file FixedTargetGenerator.h.

OnlyMuons

Bool_t FixedTargetGenerator::OnlyMuons

protected

Definition at line 118 of file FixedTargetGenerator.h.

Option

TString FixedTargetGenerator::Option

protected

Definition at line 127 of file FixedTargetGenerator.h.

PhotonCollision

Bool_t FixedTargetGenerator::PhotonCollision

protected

Definition at line 117 of file FixedTargetGenerator.h.

pot

Int_t FixedTargetGenerator::pot

protected

Definition at line 116 of file FixedTargetGenerator.h.

setByHand

Bool_t FixedTargetGenerator::setByHand

protected

Definition at line 117 of file FixedTargetGenerator.h.

start

Double_t FixedTargetGenerator::start[3]

protected

Definition at line 133 of file FixedTargetGenerator.h.

startZ

Double_t FixedTargetGenerator::startZ

protected

Definition at line 137 of file FixedTargetGenerator.h.

targetFromGeometry

Bool_t FixedTargetGenerator::targetFromGeometry

protected

Definition at line 139 of file FixedTargetGenerator.h.

targetName

TString FixedTargetGenerator::targetName

protected

Definition at line 127 of file FixedTargetGenerator.h.

tauOnly

Bool_t FixedTargetGenerator::tauOnly

protected

Definition at line 117 of file FixedTargetGenerator.h.

withEvtGen

Bool_t FixedTargetGenerator::withEvtGen

protected

Definition at line 118 of file FixedTargetGenerator.h.

withNtuple

Bool_t FixedTargetGenerator::withNtuple

protected

Definition at line 125 of file FixedTargetGenerator.h.

wspill

Double_t FixedTargetGenerator::wspill

protected

Definition at line 115 of file FixedTargetGenerator.h.

xOff

Double_t FixedTargetGenerator::xOff

protected

Definition at line 128 of file FixedTargetGenerator.h.

yOff

Double_t FixedTargetGenerator::yOff

protected

Definition at line 129 of file FixedTargetGenerator.h.

zOff

Double_t FixedTargetGenerator::zOff

protected

Definition at line 130 of file FixedTargetGenerator.h.

---

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

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