本文整理汇总了C++中AliGenPythia::SetVertexSmear方法的典型用法代码示例。如果您正苦于以下问题:C++ AliGenPythia::SetVertexSmear方法的具体用法?C++ AliGenPythia::SetVertexSmear怎么用?C++ AliGenPythia::SetVertexSmear使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类AliGenPythia
的用法示例。
在下文中一共展示了AliGenPythia::SetVertexSmear方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: Config
//.........这里部分代码省略.........
pythia->SetMomentumRange(0,99999999);
// pythia->SetJetEtaRange(-1.5, 1.5);// Final state kinematic cuts
// pythia->SetJetPhiRange(0., 360.);
// pythia->SetThetaRange(45,135);
if(proc == kPyJetJetPHOSv2)
pythia->SetForceDecay(kNeutralPion);
else
pythia->SetForceDecay(kAll);
pythia->SetPycellParameters(2., 274, 432, 0., 4., 5., 1.0);
pythia->SetPtKick(5); // set the intrinsic kt to 5 GeV/c
gener = pythia;
}
// PRIMARY VERTEX
gener->SetOrigin(0., 0., 0.); // vertex position
// Size of the interaction diamond
// Longitudinal
Float_t sigmaz = 7.55 / TMath::Sqrt(2.); // [cm]
// Transverse
Float_t betast = 10; // beta* [m]
Float_t eps = 3.75e-6; // emittance [m]
Float_t gamma = 7000. / 0.938272; // relativistic gamma [1]
Float_t sigmaxy = TMath::Sqrt(eps * betast / gamma) / TMath::Sqrt(2.) * 100.; // [cm]
printf("\n \n Diamond size x-y: %10.3e z: %10.3e\n \n", sigmaxy, sigmaz);
gener->SetSigma(sigmaxy, sigmaxy, sigmaz); // Sigma in (X,Y,Z) (cm) on IP position
gener->SetCutVertexZ(3.); // Truncate at 3 sigma
gener->SetVertexSmear(kPerEvent);
gener->Init();
//Quenching
gener->SetQuench(iquenching);
if(iquenching == 1){
Float_t k = 6e5*(qhat/1.7) ; //qhat=1.7, k = 6e5, default value
AliPythia::Instance()->InitQuenching(0.,0.1,k,0,0.95,6);
}
// FIELD
if (mag == AliMagF::k2kG) {
comment = comment.Append(" | L3 field 0.2 T");
} else if (mag == AliMagF::k5kG) {
comment = comment.Append(" | L3 field 0.5 T");
}
printf("\n \n Comment: %s \n \n", comment.Data());
// to use constant field in the barrel use:
// AliMagF* field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kGUniform);
AliMagF* field = new AliMagF("Maps","Maps", -1., -1., mag);
TGeoGlobalMagField::Instance()->SetField(field);
rl->CdGAFile();
Int_t iABSO = 1;
Int_t iACORDE= 1;
Int_t iDIPO = 1;
Int_t iEMCAL = 1;
Int_t iFMD = 1;
Int_t iFRAME = 1;
Int_t iHALL = 1;
Int_t iITS = 1;