Skip to content
Snippets Groups Projects
Commit 42516b6e authored by lintao@ihep.ac.cn's avatar lintao@ihep.ac.cn
Browse files

WIP: add the debug info.

parent 27c290e5
No related branches found
No related tags found
No related merge requests found
......@@ -59,6 +59,14 @@ from Configurables import HepMCRdr
from Configurables import GenPrinter
# gun = GtGunTool("GtGunTool")
# gun.Particles = ["pi+"]
# gun.Energies = [100.] # GeV
# gun.ThetaMins = [] # rad; 45deg
# gun.ThetaMaxs = [] # rad; 45deg
# gun.PhiMins = [] # rad; 0deg
# gun.PhiMaxs = [] # rad; 360deg
stdheprdr = StdHepRdr("StdHepRdr")
stdheprdr.Input = "/cefs/data/stdhep/CEPC250/2fermions/E250.Pbhabha.e0.p0.whizard195/bhabha.e0.p0.00001.stdhep"
......
......@@ -87,6 +87,7 @@ GtGunTool::mutate(MyHepMC::GenEvent& event) {
mcp.setSimulatorStatus(1);
// mcp.setCharge();
mcp.setTime(0.0);
mcp.setMass(mass);
// mcp.setVertex();
// mcp.setEndpoint();
......@@ -123,6 +124,11 @@ GtGunTool::mutate(MyHepMC::GenEvent& event) {
info() << "phi is fixed: " << phimin << endmsg;
}
}
debug() << "Direction: "
<< " cos(theta): " << costheta
<< " phi: " << phi
<< endmsg;
double px = p*sintheta*cos(phi);
double py = p*sintheta*sin(phi);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment