TTree *tTmp;
void SetAlias()
{
  tTmp = (TTree*) gROOT->FindObject("PhysicsTree");
  //  tTmp->SetAlias("mp","1.00866491582");
  // tTmp->SetAlias("m14o","14.008596706");
  // tTmp->SetAlias("Elab14o","71.");//MeV
  // tTmp->SetAlias("Erel14o","1.00866491582*14.008596706/(1.00866491582+14.008596706)/(1.00866491582+14.008596706)*71."); 
  tTmp->SetAlias("EPCM","ELab+4.4485492-2*sqrt(ELab*4.4485492)*cos(ThetaLab/57.4)");
  tTmp->SetAlias("ThetaPCM","57.4*(3.14-acos((EPCM+4.4485492-ELab)/2./sqrt(EPCM*4.4485492)))");


  //tTmp->SetAlias("EPCM","ELab+Erel14o-2*sqrt(ELab*Erel14o)*cos(ThetaLab)");
  //tTmp->SetAlias("ThetaPCM","57.4*(3.14-atan((EPCM+Erel14o-ELab)/2./sqrt(EPCM*Erel14o)))");


}