SwimdEdxTable Class Reference

#include <SwimdEdxTable.h>

List of all members.

Public Member Functions
virtual	~SwimdEdxTable ()
void	FilldEdxTable (SwimGeo::SwimMaterial_t material)
Static Public Member Functions
static double	GetdEdx (const SwimParticle &particle, SwimGeo::SwimMaterial_t material)
Private Member Functions
	SwimdEdxTable ()
Static Private Member Functions
static double	CaldEdx (double p, SwimGeo::SwimMaterial_t material, double *dEdxTable)

---

Detailed Description

Definition at line 16 of file SwimdEdxTable.h.

---

Constructor & Destructor Documentation

virtual SwimdEdxTable::~SwimdEdxTable

(

)

[inline, virtual]

Definition at line 19 of file SwimdEdxTable.h.

{ }

SwimdEdxTable::SwimdEdxTable

(

)

[inline, private]

Definition at line 26 of file SwimdEdxTable.h.

{ }

---

Member Function Documentation

double SwimdEdxTable::CaldEdx	(	double	p,
		SwimGeo::SwimMaterial_t	material,
		double *	dEdxTable
	)			[static, private]

Definition at line 331 of file SwimdEdxTable.cxx.

References Munits::cm, Munits::cm2, Munits::g, SwimGeo::GetSwimMaterialDensity(), Munits::GeV, SwimGeo::kIron, SwimGeo::kScint, and Munits::m.

Referenced by GetdEdx().

{
  int    i; // array index
  double p0,p1;
  double slope;
  double eLoss;
  double dEdxTableUnit = Munits::GeV/(Munits::g)*Munits::cm2;
  double density = SwimGeo::GetSwimMaterialDensity(material);

  // momentum (in GeV) are evenly distributed in log10
  // log10(p1) - log10(p0) = 0.0472440944881889757
  i = int((log10(p)+2.0)/0.0472440944881889757);
  if (i>=127)
    i = 126;
  p0 = pow(10, (-2.0 + 0.0472440944881889757*i));
  p1 = pow(10, (-2.0 + 0.0472440944881889757*(i+1)));
  slope = (dEdxTable[i+1]-dEdxTable[i])*dEdxTableUnit/
          ((p1-p0)*Munits::GeV);

  // eLoss unit: GeV/m
  eLoss = (density*(slope*(p-p0*Munits::GeV)
                   + dEdxTable[i]*dEdxTableUnit))*Munits::cm/Munits::m;

  if(material == SwimGeo::kIron){
    if(p>0) eLoss *= (1.005 + 0.005*log10(p));
    if(fabs(p)<0.1*Munits::GeV) eLoss *= 0.995;
  }
    
    if(material == SwimGeo::kScint){
      if(p>0) eLoss *= (1.0158 + 0.00511*log10(p));
    }

  
  return eLoss;
}

void SwimdEdxTable::FilldEdxTable

(

SwimGeo::SwimMaterial_t

material

)

Definition at line 372 of file SwimdEdxTable.cxx.

References MuELoss::a, SwimGeo::kIron, SwimGeo::kScint, Munits::m, and Mphysical::pi.

                                                                {
  double logpMin = log10(0.01);
  double logpMax = log10(10000.);
  double logpTable[128];
  for (int n=0; n<128; ++n)
    logpTable[n] = logpMin + (logpMax-logpMin)*n/127;

  // energyTable units: GeV
  double energyTable[45] =
  {1.0,   2.0,   4.0,   6.0,   8.0,
   10.0,  15.0,  20.0,  30.0,  40.0,
   50.0,  60.0,  70.0,  80.0,  90.0,
   100.,  120.,  140.,  160.,  180.,
   200.,  220.,  240.,  260.,  280.,
   300.,  350.,  400.,  450.,  500.,
   600.,  700.,  800.,  900.,  1000.,
   1500., 2000., 3000., 4000., 5000.,
   6000., 7000., 8000., 9000., 10000.};

  // eLossTable units: GeV*cm^2/g
  double eLossTableScint[45] =
  {0.002027, 0.002147, 0.002265, 0.002330, 0.002375,
   0.002409, 0.002470, 0.002514, 0.002581, 0.002633, 
   0.002677, 0.002718, 0.002756, 0.002792, 0.002826, 
   0.002860, 0.002925, 0.002987, 0.003049, 0.003109, 
   0.003169, 0.003228, 0.003287, 0.003345, 0.003403, 
   0.003461, 0.003605, 0.003749, 0.003892, 0.004035, 
   0.004320, 0.004606, 0.004891, 0.005178, 0.005464,
   0.006899, 0.008341, 0.01124,  0.01417,  0.01710,
   0.02005,  0.02300,  0.02597,  0.02893,  0.03191};

  double eLossTableFe[45] = 
  {0.001561, 0.001682, 0.001802, 0.001870, 0.001919,
   0.001958, 0.002032, 0.002091, 0.002189, 0.002275, 
   0.002357, 0.002436, 0.002513, 0.002590, 0.002666,
   0.002741, 0.002893, 0.003044, 0.003195, 0.003347, 
   0.003499, 0.003652, 0.003805, 0.003958, 0.004112, 
   0.004265, 0.004652, 0.005040, 0.005429, 0.005821,
   0.006607, 0.007395, 0.008189, 0.008985, 0.009782,
   0.01380,  0.01785,  0.02600,  0.03419,  0.04241,
   0.05065,  0.05892,  0.06719,  0.07548,  0.08377};

  double E;
  double p;
  double beta;
  double gamma;
  double X;
  double Em;
  double dEdx  = 0;
  double alpha = 1/137.036;
  double N     = 6.023e23;
  double me    = 0.510998902e-3;
  double mu    = 0.105658357;
  double lamda = 3.8616e-11;
  double pi    = 3.141592654;
  double delta  = 0.0;;

  double ZOverA = 0;
  double I      = 0;
  double C      = 0;
  double X0     = 0;
  double X1     = 0;
  double a      = 0;
  double m      = 0;

  if (material==SwimGeo::kScint) {
    // Scint information
    ZOverA = .5377;
    I      = 68.7e-9;
    C      = -3.300;
    X0     = .165;
    X1     = 2.503;
    a      = .165;
    m      = 3.222;
  }
  else if (material==SwimGeo::kIron) {
    // Fe information
    ZOverA = 26./55.847;
    I      = 286.0e-9;
    C      = -4.291;
    X0     = -.001;
    X1     = 3.153;
    a      = .147;
    m      = 2.963;
  }

  for (int i=0; i<128; ++i) {
    p = pow(10,logpTable[i]);
    E = TMath::Sqrt(p*p+mu*mu);

    if (p<1.0) {
      beta  = p/E;
      gamma = E/mu;
      X     = log10(beta*gamma);
      
      if (X>X0 && X<X1)
        delta = 4.6052*X+a*pow((X1-X),m)+C;
      else if (X>X1)
        delta = 4.6052*X+C;

      Em = 2*me*p*p/(me*me+mu*mu+2*me*E);

      dEdx = alpha*alpha*2*pi*N*lamda*lamda*ZOverA*me/(beta*beta)*
             (log(2*me*beta*beta*gamma*gamma*Em/(I*I))-2*beta*beta+Em*Em/(E*E)/4-delta);
    }
    else {
      int j;
      double slope;

      for (j=0; j<45; ++j)
        if (p<TMath::Sqrt(energyTable[j]*energyTable[j]-mu*mu))
          break;
      if (p>TMath::Sqrt(energyTable[44]*energyTable[44]-mu*mu))
        j = 44;

      double logp1 = log10(TMath::Sqrt(energyTable[j]*energyTable[j]-mu*mu));
      double logp2 = log10(TMath::Sqrt(energyTable[j-1]*energyTable[j-1]-mu*mu));
      if (material==SwimGeo::kScint) {
        slope = (log10(eLossTableScint[j])-log10(eLossTableScint[j-1]))/(logp1-logp2);
        dEdx = pow(10,(slope*(logpTable[i]-logp2)+log10(eLossTableScint[j-1])));
      }
      else if (material==SwimGeo::kIron) {
        slope = (log10(eLossTableFe[j])-log10(eLossTableFe[j-1]))/(logp1-logp2);
        dEdx = pow(10,(slope*(logpTable[i]-logp2)+log10(eLossTableFe[j-1])));
      }
    }
    cout << p << " " << dEdx << endl;
  }
}

double SwimdEdxTable::GetdEdx	(	const SwimParticle &	particle,
		SwimGeo::SwimMaterial_t	material
	)			[static]

Definition at line 19 of file SwimdEdxTable.cxx.

References CaldEdx(), SwimParticle::GetMomentumModulus(), Munits::GeV, SwimGeo::kAluminum, SwimGeo::kCarbon, SwimGeo::kCoil_far, SwimGeo::kCoil_near, SwimGeo::kCopper, SwimGeo::kHydrogen, SwimGeo::kIron, SwimGeo::kNitrogen, SwimGeo::kOxygen, SwimGeo::kScint, and SwimGeo::kWater.

Referenced by SwimdEdxAction::Perform().

{
  double p = particle.GetMomentumModulus();

  if (material==SwimGeo::kScint) {
    static double dEdxTableScint[128] =
      {0.090744, 0.076389, 0.064217, 0.053925, 0.045244, 
       0.037938, 0.031804, 0.026664, 0.022364, 0.018775, 
       0.015784, 0.013296, 0.011229, 0.009517, 0.008099, 
       0.006929, 0.005964, 0.005171, 0.004521, 0.003989, 
       0.003555, 0.003204, 0.002919, 0.002691, 0.002510, 
       0.002366, 0.002248, 0.002156, 0.002086, 0.002032, 
       0.001994, 0.001966, 0.001949, 0.001939, 0.001935, 
       0.001937, 0.001942, 0.001951, 0.001962, 0.001975, 
       0.001990, 0.002006, 0.002023, 0.002040, 0.002058, 
       0.002077, 0.002096, 0.002115, 0.002134, 0.002153, 
       0.002171, 0.002189, 0.002207, 0.002226, 0.002245, 
       0.002264, 0.002281, 0.002298, 0.002316, 0.002333, 
       0.002350, 0.002367, 0.002384, 0.002401, 0.002417, 
       0.002433, 0.002450, 0.002466, 0.002483, 0.002499, 
       0.002516, 0.002534, 0.002552, 0.002570, 0.002589, 
       0.002608, 0.002628, 0.002649, 0.002671, 0.002694, 
       0.002719, 0.002746, 0.002774, 0.002804, 0.002837, 
       0.002873, 0.002912, 0.002954, 0.002999, 0.003050, 
       0.003105, 0.003167, 0.003234, 0.003309, 0.003392, 
       0.003484, 0.003586, 0.003700, 0.003826, 0.003966,
       0.004124, 0.004296, 0.004492, 0.004707, 0.004947, 
       0.005217, 0.005521, 0.005878, 0.006257, 0.006661, 
       0.007120, 0.007650, 0.008219, 0.008889, 0.009629, 
       0.010432, 0.011307, 0.012342, 0.013472, 0.014730, 
       0.016144, 0.017715, 0.019479, 0.021448, 0.023643, 
       0.026103, 0.028839, 0.031910};

    return(SwimdEdxTable::CaldEdx(p, material, dEdxTableScint));
  }
  else if (material==SwimGeo::kIron) {
    static double dEdxTableFe[128] =
      {0.055599, 0.047622, 0.040664, 0.034633, 0.029434, 
       0.024973, 0.021162, 0.017918, 0.015166, 0.012839, 
       0.010878, 0.009229, 0.007846, 0.006690, 0.005726, 
       0.004924, 0.004259, 0.003710, 0.003257, 0.002885, 
       0.002581, 0.002334, 0.002116, 0.001954, 0.001824, 
       0.001722, 0.001641, 0.001580, 0.001533, 0.001499, 
       0.001476, 0.001461, 0.001453, 0.001450, 0.001452, 
       0.001458, 0.001467, 0.001478, 0.001492, 0.001506, 
       0.001522, 0.001539, 0.001556, 0.001574, 0.001593, 
       0.001611, 0.001630, 0.001649, 0.001668, 0.001688, 
       0.001706, 0.001724, 0.001743, 0.001762, 0.001781, 
       0.001801, 0.001819, 0.001837, 0.001855, 0.001874, 
       0.001892, 0.001911, 0.001929, 0.001948, 0.001968, 
       0.001987, 0.002007, 0.002027, 0.002049, 0.002071, 
       0.002094, 0.002120, 0.002146, 0.002173, 0.002202, 
       0.002234, 0.002267, 0.002305, 0.002345, 0.002390, 
       0.002438, 0.002492, 0.002551, 0.002617, 0.002689, 
       0.002771, 0.002861, 0.002962, 0.003073, 0.003197, 
       0.003337, 0.003493, 0.003668, 0.003864, 0.004082, 
       0.004327, 0.004600, 0.004907, 0.005248, 0.005631, 
       0.006063, 0.006539, 0.007076, 0.007674, 0.008343, 
       0.009092, 0.009934, 0.010894, 0.011948, 0.013104, 
       0.014403, 0.015875, 0.017497, 0.019341, 0.021394, 
       0.023665, 0.026183, 0.029039, 0.032207, 0.035744, 
       0.039702, 0.044114, 0.049044, 0.054555, 0.060704, 
       0.067561, 0.075225, 0.083770};

    return(SwimdEdxTable::CaldEdx(p,material, dEdxTableFe));
  }
  else if( material==SwimGeo::kAluminum ) {
    static double dEdxTableAl[128] =
      {0.066606, 0.056594, 0.047979, 0.040601, 0.034306,
       0.028953, 0.024417, 0.020583, 0.017352, 0.014636,
       0.012358, 0.010452, 0.008861, 0.007535, 0.006434,
       0.005520, 0.004765, 0.004142, 0.003630, 0.003210,
       0.002867, 0.002589, 0.002364, 0.002183, 0.002039,
       0.001926, 0.001827, 0.001759, 0.001708, 0.001671,
       0.001645, 0.001628, 0.001618, 0.001615, 0.001616,
       0.001622, 0.001630, 0.001641, 0.001654, 0.001669,
       0.001685, 0.001702, 0.001719, 0.001736, 0.001754,
       0.001773, 0.001792, 0.001811, 0.001830, 0.001848,
       0.001866, 0.001884, 0.001902, 0.001920, 0.001938,
       0.001957, 0.001974, 0.001991, 0.002008, 0.002025,
       0.002042, 0.002058, 0.002075, 0.002092, 0.002108,
       0.002125, 0.002142, 0.002159, 0.002177, 0.002195,
       0.002214, 0.002233, 0.002253, 0.002273, 0.002294,
       0.002318, 0.002342, 0.002368, 0.002394, 0.002425,
       0.002458, 0.002494, 0.002533, 0.002575, 0.002621,
       0.002673, 0.002730, 0.002792, 0.002861, 0.002937,
       0.003023, 0.003117, 0.003224, 0.003341, 0.003472,
       0.003620, 0.003783, 0.003966, 0.004170, 0.004397,
       0.004653, 0.004934, 0.005254, 0.005608, 0.006003,
       0.006444, 0.006944, 0.007520, 0.008143, 0.008819,
       0.009582, 0.010453, 0.011404, 0.012501, 0.013719,
       0.015056, 0.016529, 0.018221, 0.020088, 0.022172,
       0.024514, 0.027120, 0.030033, 0.033292, 0.036926,
       0.040980, 0.045510, 0.050570};
    return(SwimdEdxTable::CaldEdx(p, material, dEdxTableAl));
  }
  else if( material==SwimGeo::kCopper ) {
    static double dEdxTableCu[128] =
      {0.052630, 0.045172, 0.038644, 0.032967, 0.028059,
       0.023838, 0.020225, 0.017143, 0.014525, 0.012307,
       0.010436, 0.008860, 0.007538, 0.006432, 0.005508,
       0.004739, 0.004102, 0.003574, 0.003139, 0.002782,
       0.002489, 0.002252, 0.002048, 0.001890, 0.001765,
       0.001665, 0.001588, 0.001528, 0.001483, 0.001450,
       0.001427, 0.001413, 0.001405, 0.001403, 0.001405,
       0.001411, 0.001420, 0.001431, 0.001444, 0.001459,
       0.001474, 0.001491, 0.001508, 0.001525, 0.001544,
       0.001562, 0.001581, 0.001600, 0.001619, 0.001639,
       0.001657, 0.001675, 0.001694, 0.001712, 0.001731,
       0.001751, 0.001769, 0.001787, 0.001806, 0.001825,
       0.001843, 0.001862, 0.001880, 0.001899, 0.001919,
       0.001940, 0.001960, 0.001981, 0.002003, 0.002026,
       0.002049, 0.002076, 0.002104, 0.002131, 0.002162,
       0.002196, 0.002231, 0.002271, 0.002313, 0.002361,
       0.002411, 0.002469, 0.002532, 0.002602, 0.002680,
       0.002767, 0.002863, 0.002971, 0.003090, 0.003223,
       0.003373, 0.003540, 0.003728, 0.003938, 0.004172,
       0.004436, 0.004728, 0.005058, 0.005426, 0.005837,
       0.006301, 0.006812, 0.007391, 0.008034, 0.008754,
       0.009558, 0.010463, 0.011493, 0.012626, 0.013870,
       0.015268, 0.016852, 0.018600, 0.020581, 0.022790,
       0.025235, 0.027946, 0.031017, 0.034426, 0.038230,
       0.042487, 0.047231, 0.052532, 0.058457, 0.065069,
       0.072449, 0.080686, 0.089870};
    
    return(SwimdEdxTable::CaldEdx(p, material, dEdxTableCu));
  }
  else if( material==SwimGeo::kHydrogen ) {
    static double dEdxTableH[128] =
      {0.206816, 0.172698, 0.144099, 0.120169, 0.100177,
       0.083501, 0.069612, 0.058058, 0.048461, 0.040500,
       0.033904, 0.028447, 0.023937, 0.020216, 0.017149,
       0.014627, 0.012555, 0.010857, 0.009468, 0.008335,
       0.007413, 0.006665, 0.006062, 0.005577, 0.005191,
       0.004885, 0.004645, 0.004461, 0.004321, 0.004218,
       0.004145, 0.004098, 0.004069, 0.004052, 0.004044,
       0.004045, 0.004051, 0.004063, 0.004079, 0.004098,
       0.004119, 0.004142, 0.004167, 0.004191, 0.004219,
       0.004247, 0.004275, 0.004303, 0.004332, 0.004361,
       0.004389, 0.004418, 0.004447, 0.004477, 0.004506,
       0.004536, 0.004564, 0.004593, 0.004622, 0.004651,
       0.004679, 0.004707, 0.004735, 0.004762, 0.004790,
       0.004817, 0.004844, 0.004872, 0.004898, 0.004925,
       0.004951, 0.004978, 0.005004, 0.005031, 0.005058,
       0.005085, 0.005112, 0.005140, 0.005169, 0.005198,
       0.005227, 0.005258, 0.005290, 0.005323, 0.005358,
       0.005394, 0.005433, 0.005473, 0.005516, 0.005563,
       0.005612, 0.005666, 0.005723, 0.005786, 0.005854,
       0.005929, 0.006011, 0.006100, 0.006198, 0.006306,
       0.006427, 0.006557, 0.006704, 0.006864, 0.007042,
       0.007239, 0.007464, 0.007732, 0.008010, 0.008297,
       0.008628, 0.009015, 0.009419, 0.009915, 0.010457,
       0.011029, 0.011639, 0.012392, 0.013194, 0.014091,
       0.015111, 0.016236, 0.017497, 0.018911, 0.020483,
       0.022237, 0.024215, 0.026420};
    
     return(SwimdEdxTable::CaldEdx(p, material, dEdxTableH));
  }
  else if( material==SwimGeo::kOxygen ) {
    static double dEdxTableO[128] =
      {0.078780, 0.066518, 0.056073, 0.047204, 0.039697,
       0.033358, 0.028020, 0.023534, 0.019773, 0.016626,
       0.013998, 0.011807, 0.009985, 0.008472, 0.007218,
       0.006181, 0.005326, 0.004621, 0.004043, 0.003571,
       0.003185, 0.002872, 0.002619, 0.002416, 0.002254,
       0.002126, 0.002027, 0.001951, 0.001894, 0.001853,
       0.001823, 0.001804, 0.001793, 0.001788, 0.001788,
       0.001793, 0.001801, 0.001812, 0.001824, 0.001839,
       0.001854, 0.001871, 0.001887, 0.001904, 0.001921,
       0.001939, 0.001957, 0.001976, 0.001994, 0.002012,
       0.002028, 0.002045, 0.002062, 0.002079, 0.002096,
       0.002113, 0.002129, 0.002144, 0.002160, 0.002176,
       0.002191, 0.002207, 0.002223, 0.002238, 0.002254,
       0.002269, 0.002285, 0.002300, 0.002317, 0.002333,
       0.002350, 0.002368, 0.002387, 0.002405, 0.002424,
       0.002445, 0.002466, 0.002489, 0.002513, 0.002539,
       0.002566, 0.002595, 0.002627, 0.002662, 0.002700,
       0.002742, 0.002787, 0.002838, 0.002892, 0.002952,
       0.003019, 0.003093, 0.003176, 0.003267, 0.003368,
       0.003482, 0.003607, 0.003748, 0.003903, 0.004077,
       0.004273, 0.004487, 0.004730, 0.004999, 0.005299,
       0.005634, 0.006013, 0.006454, 0.006927, 0.007436,
       0.008011, 0.008672, 0.009388, 0.010222, 0.011148,
       0.012157, 0.013263, 0.014547, 0.015956, 0.017531,
       0.019305, 0.021275, 0.023476, 0.025941, 0.028690,
       0.031757, 0.035186, 0.039030};

     return(SwimdEdxTable::CaldEdx(p, material, dEdxTableO));
  }
  else if( material==SwimGeo::kWater ) {
    static double dEdxTableH20[128]=
      {0.091993, 0.077501, 0.065198, 0.054784, 0.045992,
       0.038587, 0.032365, 0.027147, 0.022780, 0.019132,
       0.016090, 0.013558, 0.011455, 0.009711, 0.008267,
       0.007074, 0.006091, 0.005282, 0.004619, 0.004076,
       0.003634, 0.003275, 0.002985, 0.002752, 0.002567,
       0.002420, 0.002306, 0.002218, 0.002145, 0.002091,
       0.002051, 0.002023, 0.002005, 0.001995, 0.001991,
       0.001993, 0.001998, 0.002008, 0.002019, 0.002033,
       0.002049, 0.002065, 0.002083, 0.002101, 0.002120,
       0.002140, 0.002159, 0.002179, 0.002199, 0.002219,
       0.002238, 0.002257, 0.002276, 0.002295, 0.002315,
       0.002335, 0.002353, 0.002371, 0.002390, 0.002408,
       0.002426, 0.002444, 0.002462, 0.002479, 0.002497,
       0.002514, 0.002532, 0.002550, 0.002568, 0.002586,
       0.002604, 0.002624, 0.002643, 0.002662, 0.002683,
       0.002704, 0.002726, 0.002749, 0.002774, 0.002800,
       0.002827, 0.002857, 0.002889, 0.002924, 0.002961,
       0.003003, 0.003048, 0.003096, 0.003149, 0.003208,
       0.003273, 0.003344, 0.003424, 0.003513, 0.003610,
       0.003720, 0.003840, 0.003974, 0.004125, 0.004291,
       0.004478, 0.004682, 0.004915, 0.005171, 0.005457,
       0.005777, 0.006139, 0.006562, 0.007014, 0.007497,
       0.008045, 0.008675, 0.009356, 0.010152, 0.011034,
       0.011992, 0.013039, 0.014266, 0.015609, 0.017108,
       0.018796, 0.020670, 0.022770, 0.025117, 0.027732,
       0.027483, 0.003627, 0.037580};

    return(SwimdEdxTable::CaldEdx(p, material, dEdxTableH20));
  }
  else if( material==SwimGeo::kNitrogen ) {
    static double dEdxTableN[128] =
      {0.081280, 0.068531, 0.057696, 0.048513, 0.040753,
       0.034212, 0.028710, 0.024093, 0.020227, 0.016995,
       0.014299, 0.012053, 0.010187, 0.008639, 0.007356,
       0.006297, 0.005423, 0.004704, 0.004114, 0.003631,
       0.003238, 0.002919, 0.002661, 0.002454, 0.002289,
       0.002159, 0.002057, 0.001980, 0.001922, 0.001879,
       0.001850, 0.001830, 0.001818, 0.001813, 0.001813,
       0.001818, 0.001826, 0.001836, 0.001849, 0.001863,
       0.001879, 0.001895, 0.001912, 0.001929, 0.001946,
       0.001964, 0.001982, 0.002001, 0.002019, 0.002037,
       0.002053, 0.002070, 0.002087, 0.002104, 0.002121,
       0.002138, 0.002154, 0.002169, 0.002185, 0.002201,
       0.002216, 0.002231, 0.002246, 0.002261, 0.002277,
       0.002292, 0.002308, 0.002323, 0.002339, 0.002356,
       0.002372, 0.002390, 0.002407, 0.002425, 0.002444,
       0.002463, 0.002483, 0.002505, 0.002527, 0.002552,
       0.002578, 0.002606, 0.002636, 0.002668, 0.002703,
       0.002742, 0.002784, 0.002831, 0.002881, 0.002936,
       0.002998, 0.002099, 0.003141, 0.003225, 0.003317,
       0.003421, 0.003536, 0.003664, 0.003806, 0.003964,
       0.004143, 0.004338, 0.004559, 0.004803, 0.005076,
       0.005380, 0.005725, 0.006128, 0.006559, 0.007020,
       0.007543, 0.008145, 0.008794, 0.009554, 0.010397,
       0.011313, 0.012315, 0.013485, 0.014766, 0.016197,
       0.017807, 0.019597, 0.021601, 0.023841, 0.026339,
       0.029132, 0.032256, 0.035740};

    return(SwimdEdxTable::CaldEdx(p, material, dEdxTableN));
  }
  else if( material==SwimGeo::kCarbon ) {
    static double dEdxTableC[128] =
      {0.082111, 0.069200, 0.058234, 0.048947, 0.041103,
       0.034494, 0.028938, 0.024278, 0.020376, 0.017116,
       0.014398, 0.012134, 0.010253, 0.008693, 0.007402,
       0.006334, 0.005455, 0.004731, 0.004137, 0.003651,
       0.003256, 0.002934, 0.002429, 0.002244, 0.002096,
       0.001979, 0.001887, 0.001816, 0.001763, 0.001724,
       0.001697, 0.001679, 0.001670, 0.001666, 0.001668,
       0.001675, 0.001684, 0.001696, 0.001711, 0.001727,
       0.001744, 0.001763, 0.001782, 0.001842, 0.001860,
       0.001877, 0.001895, 0.001913, 0.001930, 0.001948,
       0.001965, 0.001982, 0.001999, 0.002017, 0.002034,
       0.002052, 0.002068, 0.002084, 0.002100, 0.002116,
       0.002132, 0.002148, 0.002164, 0.002179, 0.002195,
       0.002210, 0.002226, 0.002241, 0.002257, 0.002273,
       0.002289, 0.002306, 0.002323, 0.002340, 0.002358,
       0.002376, 0.002395, 0.002416, 0.002438, 0.002461,
       0.002485, 0.002511, 0.002539, 0.002569, 0.002601,
       0.002637, 0.002676, 0.002718, 0.002764, 0.002815,
       0.002871, 0.002933, 0.003001, 0.003078, 0.003162,
       0.003256, 0.003360, 0.003476, 0.003604, 0.003746,
       0.003907, 0.004083, 0.004283, 0.004503, 0.004748,
       0.005023, 0.005334, 0.005697, 0.006085, 0.006500,
       0.006969, 0.007511, 0.008095, 0.008779, 0.009536,
       0.010359, 0.011258, 0.012313, 0.013466, 0.009176,
       0.003280, 0.002969, 0.012916, 0.021635, 0.023886,
       0.026397, 0.029207, 0.032350};

    return(SwimdEdxTable::CaldEdx(p, material, dEdxTableC));
  }
  else if ( material==SwimGeo::kCoil_near) {
    // Fraction by weight 
    Float_t f_H20 = 0.085313;
    Float_t f_Al  = 0.914687;
  
    return
      f_Al*SwimdEdxTable::GetdEdx(particle,SwimGeo::kAluminum)+
      f_H20*SwimdEdxTable::GetdEdx(particle,SwimGeo::kWater);
  }
  else if ( material==SwimGeo::kCoil_far) {
    // materials in far detector coil foam from mix_farcoil.F
    Float_t f_H =  0.012814;                  
    Float_t f_C =  0.060215;
    Float_t f_N =  0.000197;
    Float_t f_O =  0.021346;
    Float_t f_Cu = 0.905425; 
    return 
      f_H*SwimdEdxTable::GetdEdx(particle,SwimGeo::kHydrogen)+
      f_C*SwimdEdxTable::GetdEdx(particle,SwimGeo::kCarbon) +
      f_N*SwimdEdxTable::GetdEdx(particle,SwimGeo::kNitrogen)+
      f_O*SwimdEdxTable::GetdEdx(particle,SwimGeo::kOxygen)  +
      f_Cu*SwimdEdxTable::GetdEdx(particle,SwimGeo::kCopper);
  }
  else
    return 0.0*Munits::GeV;
}

---

The documentation for this class was generated from the following files:

• SwimdEdxTable.h
• SwimdEdxTable.cxx

---