Skip to content
Snippets Groups Projects
Commit a4dd0069 authored by Gianfranco Ingratta's avatar Gianfranco Ingratta
Browse files

creating a tEvent tree beside the anaevt tree (restored since last commit)

parent 55a7a62b
No related branches found
No related tags found
1 merge request!2#1 harmonize output format with agreed one
......@@ -12,7 +12,6 @@ find_package(ROOT REQUIRED COMPONENTS Geom EG)
#find_package(SandReco)
include_directories(/data/thesis/sand-reco/include)
ROOT_GENERATE_DICTIONARY(G__Struct struct.h LINKDEF Linkdef.h)
add_library(Struct SHARED G__Struct.cxx)
......
# how to compile
source compile.sh
# how to build
mkdir build
/cmake -S . -B build
cd build; make; cd ..
# how run
./build/sandSmearGo cc 0 ${EDEPSIM_INPUT} ${FILE_OUTPUT}
......
......@@ -293,6 +293,7 @@ void fill1Par(double recoPx, double recoPy, double recoPz, int trackid, double l
fast_reco::particle par;
int pdg=Event->Trajectories[trackid].PDGCode;
double m=dbpdg->GetParticle(pdg)->Mass()*1000; // MeV
double q=dbpdg->GetParticle(pdg)->Charge()/3;
par.pdg=pdg;
par.recoP4[0]=recoPx;
par.recoP4[1]=recoPy;
......@@ -313,9 +314,9 @@ void fill1Par(double recoPx, double recoPy, double recoPz, int trackid, double l
particle p;
//p.primary =
p.pdg = pdg;
//p.tid = trackid;
p.tid = trackid;
p.parent_tid=findTopParent(trackid);
//p.charge =
p.charge = q;
p.mass = m;
p.pxtrue = Event->Trajectories[trackid].InitialMomentum.X();
p.pytrue = Event->Trajectories[trackid].InitialMomentum.Y();
......
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