本文整理汇总了C++中AliGenPythia::SetForceDecay方法的典型用法代码示例。如果您正苦于以下问题：C++ AliGenPythia::SetForceDecay方法的具体用法？C++ AliGenPythia::SetForceDecay怎么用？C++ AliGenPythia::SetForceDecay使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类AliGenPythia的用法示例。
在下文中一共展示了AliGenPythia::SetForceDecay方法的4个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: Pythia
AliGenerator* MbPythiaTunePerugia0BtoJpsi2e()
{
comment = comment.Append(" pp: Pythia (Perugia0) BtoJpsi (1 bbbar per event, 1 b-hadron in |y|<2, 1 J/psi in |y|<2");
//
// Pythia
AliGenPythia* pythia = new AliGenPythia(-1);
pythia->SetMomentumRange(0, 999999.);
pythia->SetThetaRange(0., 180.);
pythia->SetYRange(-2.,2.);
pythia->SetPtRange(0,1000.);
pythia->SetProcess(kPyBeautyppMNRwmi);
pythia->SetEnergyCMS(energy);
// Tune
// 320 Perugia 0
pythia->SetTune(320);
pythia->UseNewMultipleInteractionsScenario();
//
// decays
pythia->SetCutOnChild(1);
pythia->SetPdgCodeParticleforAcceptanceCut(443);
pythia->SetChildYRange(-2,2);
pythia->SetChildPtRange(0,10000.);
//
// decays
pythia->SetForceDecay(kBJpsiDiElectron);
return pythia;
}示例2: MbPythia
AliGenerator* MbPythia() {
comment = comment.Append(" pp: Pythia low-pt");
// Pythia
AliGenPythia* pythia = new AliGenPythia(-1);
pythia->SetMomentumRange(0, 999999.);
// pythia->SetThetaRange(0., 180.);
// pythia->SetChildYRange(-12.,0.);
// pythia->SetPtRange(0,1000.);
// pythia->SetCutOnChild(1);
pythia->SetProcess(kPyMb);
pythia->SetEnergyCMS(energy);
pythia->SetForceDecay(kSemiMuonic);
return pythia;
}示例3: switch
AliGenPythia *PythiaHard(PDC07Proc_t proc) {
//*******************************************************************//
// Configuration file for hard QCD processes generation with PYTHIA //
// //
//*******************************************************************//
AliGenPythia * gener = 0x0;
switch(proc) {
case kPyJetJet:
comment = comment.Append(" pp->jet + jet over at 14 TeV, no restriction");
AliGenPythia * gener = new AliGenPythia(nEvts);
gener->SetEnergyCMS(eCMS);// Centre of mass energy
gener->SetProcess(kPyJets);// Process type
gener->SetJetEtaRange(-1.5, 1.5);// Final state kinematic cuts
gener->SetJetPhiRange(0., 360.);
gener->SetJetEtRange(10., 1000.);
gener->SetPtHard(ptHardMin, ptHardMax);// Pt transfer of the hard scattering
gener->SetStrucFunc(kCTEQ4L);
return gener;
case kPyGammaJetPHOS:
comment = comment.Append(" pp->jet + gamma over PHOS");
gener = new AliGenPythia(nEvts);
gener->SetEnergyCMS(eCMS);
gener->SetProcess(kPyDirectGamma);
gener->SetStrucFunc(kCTEQ4L);
gener->SetPtHard(ptHardMin,ptHardMax);
//gener->SetYHard(-1.0,1.0);
gener->SetGammaEtaRange(-0.13,0.13);
gener->SetGammaPhiRange(218.,322.);//Over 5 modules +-2 deg
break;
case kPyJetJetPHOS:
comment = comment.Append(" pp->jet + jet over PHOS");
gener = new AliGenPythia(nEvts);
gener->SetEnergyCMS(eCMS);
gener->SetProcess(kPyJets);
gener->SetStrucFunc(kCTEQ4L);
gener->SetPtHard(ptHardMin,ptHardMax);
//gener->SetYHard(-1.0,1.0);
gener->SetJetEtaRange(-1.,1.);
gener->SetJetPhiRange(200.,340.);
gener->SetPi0InPHOS(kTRUE);
gener->SetFragPhotonOrPi0MinPt(ptGammaPi0Min);
printf("\n \n Event generator: Minimum pT of particle in calorimeter %f \n \n", ptGammaPi0Min);
break;
case kPyGammaBremsPHOS:
comment = comment.Append(" pp->jet + jet+bremsphoton over PHOS at 14 TeV");
gener = new AliGenPythia(nEvts);
gener->SetEnergyCMS(eCMS);
gener->SetProcess(kPyJets);
gener->SetStrucFunc(kCTEQ4L);
gener->SetPtHard(ptHardMin,ptHardMax);
//gener->SetYHard(-1.0,1.0);
gener->SetJetEtaRange(-1.,1.);
gener->SetJetPhiRange(200.,340.);
gener->SetFragPhotonInPHOS(kTRUE);
gener->SetFragPhotonOrPi0MinPt(ptGammaPi0Min);
printf("\n \n Event generator: Minimum pT of particle in calorimeter %f \n \n", ptGammaPi0Min);
break;
case kPyJetJetPHOSv2:
comment = comment.Append(" pp->jet + jet over PHOS version2 ");
gener = new AliGenPythia(nEvts);
gener->SetEnergyCMS(eCMS);
gener->SetProcess(kPyJets);
gener->SetStrucFunc(kCTEQ4L);
gener->SetPtHard(ptHardMin,ptHardMax);
//gener->SetYHard(-1.0,1.0);
gener->SetJetEtaRange(-1.,1.);
gener->SetJetPhiRange(200.,340.);
//gener->SetPi0InPHOS(kTRUE);
gener->SetPhotonInPHOSeta(kTRUE);
gener->SetPhotonMinPt(ptGammaPi0Min);
gener->SetForceDecay(kAll);
break;
case kPyGammaJetEMCAL:
comment = comment.Append(" pp->jet + gamma over EMCAL at 14 TeV");
gener = new AliGenPythia(nEvts);
gener->SetEnergyCMS(eCMS);
gener->SetProcess(kPyDirectGamma);
gener->SetStrucFunc(kCTEQ4L);
gener->SetPtHard(ptHardMin,ptHardMax);
//gener->SetYHard(-1.0,1.0);
gener->SetGammaEtaRange(-0.71,0.71);
gener->SetGammaPhiRange(78.,192.);//Over 6 supermodules +-2 deg
break;
case kPyJetJetEMCAL:
comment = comment.Append(" pp->jet + jet over EMCAL at 14 TeV");
gener = new AliGenPythia(nEvts);
gener->SetEnergyCMS(eCMS);
gener->SetProcess(kPyJets);
gener->SetStrucFunc(kCTEQ4L);
gener->SetPtHard(ptHardMin,ptHardMax);
//gener->SetYHard(-1.0,1.0);
gener->SetJetEtaRange(-1,1);
gener->SetJetPhiRange(60.,210.);
gener->SetPi0InEMCAL(kTRUE);
gener->SetFragPhotonOrPi0MinPt(ptGammaPi0Min);
//.........这里部分代码省略.........
示例4: Config
void Config()
{
// Get settings from environment variables
ProcessEnvironmentVars();
// Random Number seed
if (gSystem->Getenv("CONFIG_SEED")) {
seed = atoi(gSystem->Getenv("CONFIG_SEED"));
}
else if(gSystem->Getenv("DC_EVENT") && gSystem->Getenv("DC_RUN")){
seed = runNumber * 100000 + eventNumber;
}
gRandom->SetSeed(seed);
cout<<"Seed for random number generation= "<<gRandom->GetSeed()<<endl;
// libraries required by geant321
#if defined(__CINT__)
gSystem->Load("liblhapdf"); // Parton density functions
gSystem->Load("libEGPythia6"); // TGenerator interface
gSystem->Load("libpythia6"); // Pythia
gSystem->Load("libAliPythia6"); // ALICE specific implementations
gSystem->Load("libgeant321");
#endif
new TGeant3TGeo("C++ Interface to Geant3");
// Output every 100 tracks
((TGeant3*)gMC)->SetSWIT(4,100);
//=======================================================================
// Run loader
AliRunLoader* rl=0x0;
AliLog::Message(AliLog::kInfo, "Creating Run Loader", "", "", "Config()"," ConfigPPR.C", __LINE__);
rl = AliRunLoader::Open("galice.root",
AliConfig::GetDefaultEventFolderName(),
"recreate");
if (rl == 0x0)
{
gAlice->Fatal("Config.C","Can not instatiate the Run Loader");
return;
}
rl->SetCompressionLevel(2);
rl->SetNumberOfEventsPerFile(100);
gAlice->SetRunLoader(rl);
// Set the trigger configuration
AliSimulation::Instance()->SetTriggerConfig(TrigConfName[trig]);
cout<<"Trigger configuration is set to "<<TrigConfName[trig]<<endl;
//======================//
// Set External decayer //
//======================//
TVirtualMCDecayer* decayer = new AliDecayerPythia();
// DECAYS
//
switch(decHvFl) {
case kNature:
decayer->SetForceDecay(kAll);
break;
case kHadr:
decayer->SetForceDecay(kHadronicD);
break;
case kSemiEl:
decayer->SetForceDecay(kSemiElectronic);
break;
case kSemiMu:
decayer->SetForceDecay(kSemiMuonic);
break;
}
decayer->Init();
gMC->SetExternalDecayer(decayer);
if(proc == kPyJetJetPHOSv2) // in this case we decay the pi0
decayer->SetForceDecay(kNeutralPion);
//
//=======================================================================
// ************* STEERING parameters FOR ALICE SIMULATION **************
// --- Specify event type to be tracked through the ALICE setup
// --- All positions are in cm, angles in degrees, and P and E in GeV
gMC->SetProcess("DCAY",1);
gMC->SetProcess("PAIR",1);
gMC->SetProcess("COMP",1);
gMC->SetProcess("PHOT",1);
gMC->SetProcess("PFIS",0);
gMC->SetProcess("DRAY",0);
gMC->SetProcess("ANNI",1);
gMC->SetProcess("BREM",1);
gMC->SetProcess("MUNU",1);
gMC->SetProcess("CKOV",1);
gMC->SetProcess("HADR",1);
gMC->SetProcess("LOSS",2);
gMC->SetProcess("MULS",1);
gMC->SetProcess("RAYL",1);
Float_t cut = 1.e-3; // 1MeV cut by default
Float_t tofmax = 1.e10;
//.........这里部分代码省略.........