pionana Namespace Reference

Classes
class	TruthAnalyzer

Functions
std::vector< int >	MakeSlices (double E0, double Ef, double p, const simb::MCParticle *part)
std::map< size_t, double >	GetEDepByTraj (const simb::MCParticle *part, int id, const std::vector< sim::SimEnergyDeposit > &dep_vec, const sim::ParticleList &plist)
std::map< size_t, std::vector< int > >	GetEMDaughterByTraj (const simb::MCParticle *part, const sim::ParticleList &plist)
double	gamma (double KE, const simb::MCParticle *part)
double	beta (double KE, const simb::MCParticle *part)
double	Tmax (double KE, const simb::MCParticle *part)
double	dEdX (double KE, const simb::MCParticle *part)
double	MPV (double KE, double p, const simb::MCParticle *part)

Variables
double	K = .30705
double	A = 39.95
double	Z = 18
double	I = 188.e-6
double	me = .511
double	rho = 1.39

Function Documentation

double pionana::beta	(	double	KE,
		const simb::MCParticle *	part
	)

Definition at line 59 of file TruthAnalyzer_module.cc.

                                                      {
    return sqrt(1. - 1./(gamma(KE, part)*gamma(KE, part)));
  }

double pionana::dEdX	(	double	KE,
		const simb::MCParticle *	part
	)

Definition at line 67 of file TruthAnalyzer_module.cc.

                                                      {
    double dedx = (rho*K*Z/A);
    dedx /= (beta(KE, part)*beta(KE, part));
    double post_factor = .5*log(2*me*(beta(KE, part)*beta(KE, part))*(gamma(KE, part)*gamma(KE, part))*Tmax(KE, part)/(I*I));
    post_factor -= beta(KE, part)*beta(KE, part);
    dedx *= post_factor;
    return dedx;
  }

double pionana::gamma	(	double	KE,
		const simb::MCParticle *	part
	)

Definition at line 55 of file TruthAnalyzer_module.cc.

                                                       {
    return (KE + part->Mass())/part->Mass();
  }

std::map< size_t, double > pionana::GetEDepByTraj	(	const simb::MCParticle *	part,
		int	id,
		const std::vector< sim::SimEnergyDeposit > &	dep_vec,
		const sim::ParticleList &	plist
	)

Definition at line 96 of file TruthAnalyzer_module.cc.

                                   {

  const simb::MCTrajectory & traj = part->Trajectory();
  
  //First build up the trajectory points
  std::map<size_t, double> results;
  for (size_t i = 0; i < traj.size() - 1; ++i) {
    results[i] = 0.;
  }

  std::map<size_t, std::vector<int>> daughters_by_traj
      = GetEMDaughterByTraj(part, plist);
  
  //Iterate over the deposits
  for (const sim::SimEnergyDeposit & dep : dep_vec) {
    if (dep.TrackID() == id) {
      for (size_t i = 1; i < traj.size(); ++i) {
        if (traj.Z(i-1) <= dep.MidPointZ() && traj.Z(i) > dep.MidPointZ()) {
          results[i-1] += dep.Energy();
          break;
        }
      }
    }
    else {
      for (auto it = daughters_by_traj.begin();
           it != daughters_by_traj.end(); ++it) {
        for (size_t i = 0; i < it->second.size(); ++i) {
          if (it->second[i] == dep.TrackID()) {
            results[it->first] += dep.Energy();
          }
        }
      }
    }
  }
  return results;
}

std::map< size_t, std::vector< int > > pionana::GetEMDaughterByTraj	(	const simb::MCParticle *	part,
		const sim::ParticleList &	plist
	)

Definition at line 136 of file TruthAnalyzer_module.cc.

                                   {

  const simb::MCTrajectory & traj = part->Trajectory();
  std::map<size_t, std::vector<int>> results;

  for (int i = 0; i < part->NumberDaughters(); ++i) {
    int id = part->Daughter(i);
    auto daughter = plist[id];   
    std::string process = daughter->Process();

    if ((process.find("Ioni") == std::string::npos) &&
        (process.find("Brem") == std::string::npos) &&
        (process.find("annihil") == std::string::npos)) {
      continue;
    }

    for (size_t j = 1; j < traj.size(); ++j) {
      if (traj.Z(j-1) < daughter->Position().Z() &&
          daughter->Position().Z() < traj.Z(j)) {
        //results[j-1].push_back(id);
        std::deque<int> downstream; 
        downstream.push_back(id);


        while (!downstream.empty()) {
          int d_id = downstream.front();
          auto d_part = plist[d_id];
          results[j-1].push_back(d_id);
          for (int k = 0; k < d_part->NumberDaughters(); ++k) {
            downstream.push_back(d_part->Daughter(k));
          }
          downstream.pop_front();
        }
        break; 
      }
    }
  }

  return results;
}

std::vector< int > pionana::MakeSlices	(	double	E0,
		double	Ef,
		double	p,
		const simb::MCParticle *	part
	)

Definition at line 82 of file TruthAnalyzer_module.cc.

                                                                                              {
  std::vector<int> slices;
  int s = 0;
  while (E0 > Ef) {
    slices.push_back(s);
    ++s;
    //std::cout << E0 - 1.e3*part->Mass() << " " << dEdX(E0 - 1.e3*part->Mass(), part) << " " << MPV(E0 - 1.e3*part->Mass(), part)<< std::endl;
    //E0 -= dEdX(E0 - 1.e3*part->Mass(), part)*p;
    E0 -= MPV(E0 - 1.e3*part->Mass(), p, part);
  }

  return slices;
}

double pionana::MPV	(	double	KE,
		double	p,
		const simb::MCParticle *	part
	)

Definition at line 76 of file TruthAnalyzer_module.cc.

                                                               {
    return (K/2)*(Z/A)*(p*rho/(beta(KE, part)*beta(KE, part)))*(log(2*me*beta(KE, part)*beta(KE, part)*gamma(KE, part)*gamma(KE, part)/I) + log((K/2)*(Z/A)*p*rho/(beta(KE, part)*beta(KE, part)*I)) + .2 - beta(KE, part)*beta(KE, part));
  }

double pionana::Tmax	(	double	KE,
		const simb::MCParticle *	part
	)

Definition at line 63 of file TruthAnalyzer_module.cc.

                                                      {
    return 2*me*(beta(KE, part)*beta(KE, part))*(gamma(KE, part)*gamma(KE, part))/(1 + 2*gamma(KE, part)*me/part->Mass() + (me*me)/(part->Mass()*part->Mass()));
  }

Variable Documentation

double pionana::A = 39.95

Definition at line 38 of file TruthAnalyzer_module.cc.

double pionana::I = 188.e-6

Definition at line 40 of file TruthAnalyzer_module.cc.

double pionana::K = .30705

Definition at line 37 of file TruthAnalyzer_module.cc.

double pionana::me = .511

Definition at line 41 of file TruthAnalyzer_module.cc.

double pionana::rho = 1.39

Definition at line 42 of file TruthAnalyzer_module.cc.

double pionana::Z = 18

Definition at line 39 of file TruthAnalyzer_module.cc.