Skip to content

add nCapture time, add hittime, hit charge #11

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
dustyattenborough wants to merge 1 commit into
base: RENE
Choose a base branch
from
Open

add nCapture time, add hittime, hit charge #11

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
7 changes: 7 additions & 0 deletions src/user/include/RAT/UserAnalProc.hh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -44,6 +44,11 @@ class UserAnalProc : public Processor {
double b_Gen_Px, b_Gen_Py, b_Gen_Pz, b_Gen_KE; // [MeV/C],[MeV/C],[MeV/C],[MeV]
int b_Gen_PdgId;

std::vector<double> b_MC_HitTime_PMT0; // photon hit time at PMT one by one
std::vector<double> b_MC_HitTime_PMT1;
std::vector<double> b_MC_HitCharge_PMT0; // photon charge one by one
std::vector<double> b_MC_HitCharge_PMT1;

protected:
// For the neutron capture analysis... to be moved in a separate sub-class
std::vector<double> b_Ncap_GammaKE; // [MeV]
Expand All @@ -52,6 +57,8 @@ class UserAnalProc : public Processor {
std::string b_Ncap_pName;
std::string b_Ncap_Volume; // Volume where neutron captured , ex) "target_GdLS", "target_vessel"
double b_Ncap_X, b_Ncap_Y, b_Ncap_Z; // [mm],[mm],[mm]
double b_Ncap_gTime; // [ns], captured time of neutron to Gd or sth.

};

} // namespace RAT
Expand Down
27 changes: 27 additions & 0 deletions src/user/src/UserAnalProc.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -82,6 +82,12 @@ bool UserAnalProc::OpenFile(const std::string fileName) {

outTree_->Branch("Gen_PdgId", &b_Gen_PdgId, "Gen_PdgId/I");

outTree_->Branch("MC_HitTime_PMT0", &b_MC_HitTime_PMT0);
outTree_->Branch("MC_HitTime_PMT1", &b_MC_HitTime_PMT1);
outTree_->Branch("MC_HitCharge_PMT0", &b_MC_HitCharge_PMT0);
outTree_->Branch("MC_HitCharge_PMT1", &b_MC_HitCharge_PMT1);


// Variables for neutron capture analysis
outTree_->Branch("Ncap_GammaKE", &b_Ncap_GammaKE);
outTree_->Branch("Ncap_SumGammaKE", &b_Ncap_SumGammaKE, "Ncap_SumGammaKE/D");
Expand All @@ -92,6 +98,7 @@ bool UserAnalProc::OpenFile(const std::string fileName) {
outTree_->Branch("Ncap_X", &b_Ncap_X, "Ncap_X/D");
outTree_->Branch("Ncap_Y", &b_Ncap_Y, "Ncap_Y/D");
outTree_->Branch("Ncap_Z", &b_Ncap_Z, "Ncap_Z/D");
outTree_->Branch("Ncap_gTime", &b_Ncap_gTime, "Ncap_gTime/D");

return true;
}
Expand Down Expand Up @@ -122,12 +129,19 @@ Processor::Result UserAnalProc::DSEvent(DS::Root *ds) {
b_Gen_Px = b_Gen_Py = b_Gen_Pz = b_Gen_KE = 0;
b_Gen_PdgId = 0;

b_MC_HitTime_PMT0.clear();
b_MC_HitTime_PMT1.clear();
b_MC_HitCharge_PMT0.clear();
b_MC_HitCharge_PMT1.clear();


b_Ncap_GammaKE.clear();
b_Ncap_SumGammaKE = b_Ncap_Edeposit = 0;
b_Ncap_PdgId = 0;
b_Ncap_pName = "";
b_Ncap_Volume = "";
b_Ncap_X = b_Ncap_Y = b_Ncap_Z = 0;
b_Ncap_gTime = 0;

for (int subev = 0; subev < ds->GetEVCount(); ++subev) {
auto event = ds->GetEV(subev);
Expand All @@ -151,6 +165,7 @@ Processor::Result UserAnalProc::DSEvent(DS::Root *ds) {
const std::string pName_Ncap = track_Ncap->GetParticleName();
auto step_Ncap = track_Ncap->GetMCTrackStep(track_Ncap->GetMCTrackStepCount() - 1);
b_Ncap_Volume = step_Ncap->GetVolume();
b_Ncap_gTime = step_Ncap->GetGlobalTime();

auto endPoint_Ncap = step_Ncap->GetEndpoint();
b_Ncap_X = endPoint_Ncap.X();
Expand Down Expand Up @@ -223,6 +238,18 @@ Processor::Result UserAnalProc::DSEvent(DS::Root *ds) {
}
const int nHit = mcPMT->GetMCPhotonCount();
b_PMT_NPE[pmtID] += nHit;
for ( int iHit=0; iHit<nHit; ++iHit ) {
auto hit = mcPMT->GetMCPhoton(iHit);
if( pmtID==0){
b_MC_HitTime_PMT0.push_back(hit->GetHitTime());
b_MC_HitCharge_PMT0.push_back(hit->GetCharge());
}
if( pmtID==1){
b_MC_HitTime_PMT1.push_back(hit->GetHitTime());
b_MC_HitCharge_PMT1.push_back(hit->GetCharge());
}
}

}

outTree_->Fill();
Expand Down