#ifndef L1AddMaterial_h
#define L1AddMaterial_h

#include "L1Types.h"
#include "L1MaterialInfo.h"
#include "L1TrackPar.h"

#define cnst static const fvec 

inline void L1AddMaterial( L1TrackPar &T, L1MaterialInfo &info, fvec qp0 )
{
  cnst ZERO = 0.0, ONE = 1.;
  cnst mass2 = 0.1395679*0.1395679;

  fvec tx = T.tx;
  fvec ty = T.ty;
  fvec txtx = tx*tx;
  fvec tyty = ty*ty;
  fvec txtx1 = txtx + ONE;
  fvec h = txtx + tyty;
  fvec t = sqrt(txtx1 + tyty);
  fvec h2 = h*h;
  fvec qp0t = qp0*t;
  
  cnst c1=0.0136, c2=c1*0.038, c3=c2*0.5, c4=-c3/2.0, c5=c3/3.0, c6=-c3/4.0;
    
  fvec s0 = (c1+c2*info.logRadThick + c3*h + h2*(c4 + c5*h +c6*h2) )*qp0t;    
  fvec a = ( (ONE+mass2*qp0*qp0t)*info.RadThick*s0*s0 );

  T.C22 += txtx1*a;
  T.C32 += tx*ty*a; T.C33 += (ONE+tyty)*a; 
}

#undef cnst

#endif