char cond[1000];
TTree *t1;
void SetAlias()
{
  gStyle->SetPalette(1);

  t1 = (TTree*) gROOT->FindObject("AD")
;
  t1->SetAlias("fExist_DC_Y","DC0_Y>-1500&&DC1_Y>-1500&&DC2_Y>-1500&&DC3_Y>-1500"); 
  t1->SetAlias("fExist_IC","IC[0]>0&&IC[1]>0&&IC[2]>0&&IC[3]>0"); 

  t1->SetAlias("mT","(T_FPMW_CATS2_C*-1 + 305.3 + -4*(MWTVN==8)-1*(MWTVN==9)-1*(MWTVN==10)-2.2*(MWTVN==11)-0.7*(MWTVN==12)-0.5*(MWTVN==15)-0.6*(MWTVN==16)-1.1*(MWTVN==17)+0.6*(MWTVN==18)-1*(MWTVN==6)-0.5*(MWTVN==5)-0.5*(MWTVN==4)-1.5*(MWTVN==3))*(T_FPMW_CATS1_C>10)");
  t1->SetAlias("mD","D");
  t1->SetAlias("mV","(mD)/mT");
  t1->SetAlias("mBeta","mV/29.9792");
  t1->SetAlias("mGamma","1./sqrt(1.0-mBeta*mBeta)");
  
  t1->SetAlias("mM_Q","Brho/3.107/mBeta/mGamma");
  t1->SetAlias("mM","(mE)/931.5016/(mGamma-1.)");
  t1->SetAlias("mMr","mM_Q*((13.*qmq13+14.*qmq14+15.*qmq15+16.*qmq16))"); 
  t1->SetAlias("mQ","mM/mM_Q");
 
  //t1->SetAlias("mT1","(T_FPMW_HF_C*-1. + 559. + 3*(MWTVN==8))*(T_FPMW_HF_C>10)");
  t1->SetAlias("mT1","1.5+((23.7593+1.03259*MWTV)*(MWTVN==0)+(-9.53168+1.0711*MWTV)*(MWTVN==1)+(-30.486+1.11554*MWTV)*(MWTVN==2)+(-31.1261+1.11624*MWTV)*(MWTVN==3)+(-62.99829+1.1569*MWTV)*(MWTVN==4)+(-47.9021+1.13413*MWTV)*(MWTVN==5)+(-14.3957+1.08086*MWTV)*(MWTVN==6)+(-24.6655+1.1058*MWTV)*(MWTVN==7)+(60.4375+0.981809*MWTV)*(MWTVN==8)+(4.44472+1.06168*MWTV)*(MWTVN==9)+(11.408+1.04492*MWTV)*(MWTVN==10)+(56.8282+0.990287*MWTV)*(MWTVN==11)+(59.6833+0.995356*MWTV)*(MWTVN==12)+(27.6344+1.02846*MWTV)*(MWTVN==13)+(6.02794+1.06554*MWTV)*(MWTVN==14)+(20.1561+1.0342*MWTV)*(MWTVN==15)+(51.0718+1.01008*MWTV)*(MWTVN==16)+(36.5605+1.0391*MWTV)*(MWTVN==17)+(-6.01158+1.09581*MWTV)*(MWTVN==18)+(37.6143+1.00939*MWTV)*(MWTVN==19)+(29.3289+1.01291*MWTV)*(MWTVN==20) + 3*83.4*(GATCONF_MASTER==8)-0*(GATCONF_MASTER==1)-2*83.4*(GATCONF_MASTER==4))");
  t1->SetAlias("mT1","540-2*T_FPMW_CATS2_C");
  t1->SetAlias("mE1","(IC[0]+IC[1]+IC[2]+IC[3]+IC[4])*(IC[0]>0.1&&IC[1]>0.1&&IC[2]>0.1&&IC[3]>0.1)");
  //  t1->SetAlias("mE1","(IC[1]+IC[2]+IC[3])*(IC[1]>1&&IC[2]>1&&IC[3]>1)");
  t1->SetAlias("mdE0","IC[0]*(IC[0]>1)");
  t1->SetAlias("mdE1","IC[0]*(IC[0]>1)+IC[1]*(IC[1]>1)");
  t1->SetAlias("mD1","Path");
  t1->SetAlias("mV1","mD1/(mT1)");
  t1->SetAlias("mBeta1","mV1/29.9792");
  t1->SetAlias("mGamma1","1./sqrt(1.0-mBeta1*mBeta1)");
  
  t1->SetAlias("mM_Q1","Brho/3.107/mBeta1/mGamma1");
  t1->SetAlias("mM1","1.3*(mE1)/931.5016/(mGamma1-1.)");
  t1->SetAlias("mMr","mM_Q*((13.*qmq13+14.*qmq14+15.*qmq15+16.*qmq16))"); 
  t1->SetAlias("mQ1","mM1/mM_Q1");
  t1->SetAlias("ProjD","-174.*1.");
  t1->SetAlias("Xt","TMW1_X+ProjD*tan(TP_Theta/1000.)");
  t1->SetAlias("Yt","TMW1_Y+ProjD*tan(TP_Phi/1000.)");

 
  t1->SetAlias("mT2","540.5-2.*T_FPMW_CATS2_C+2.7*(MWTVN==16)+2.7*(MWTVN==15)+2.9*(MWTVN==14)+2.9*(MWTVN==13)+2.4*(MWTVN==12)+1.3*(MWTVN==11)+1.5*(MWTVN==10)+1.6*(MWTVN==9)-0.6*(MWTVN==8)+2.5*(MWTVN==7)+2.*(MWTVN==6)+1.6*(MWTVN==5)+1.1*(MWTVN==4)-0.6*(MWTVN==3)-1.2*(MWTVN==2)-4.0*(MWTVN==1)");

 //t1->SetAlias("mT2","mT1-TMWT1_MWT2_C+31");
  //  t1->SetAlias("mT2","mT1-4.83");

  t1->SetAlias("mdE2","(IC[0]+IC[1]*(IC[0]>0.1)");

 t1->SetAlias("mE2","0.85*(pow(Xf+450*tan(Tf/1000.),2)*0.0000015*IC[0] + 2.2*IC[0]+1.1*IC[1]*(IC[0]>0.1)+1.34*1.95*IC[2]*(IC[0]>0.1)*(IC[1]>0.1)+3.1*IC[3]*(IC[0]>0.1)*(IC[1]>0.1)*(IC[2]>0.1)+1.05*IC[4]*(IC[0]>0.1)*(IC[1]>0.1)*(IC[2]>0.1)*(IC[3]>0.1)+IC[5]*(IC[0]>0.1)*(IC[1]>0.1)*(IC[2]>0.1)*(IC[3]>0.1)*(IC[4]>0.1))");

 // t1->SetAlias("mE2","0.825*(2.4*IC[0]+IC[1]*(IC[0]>0.1)+1.5*1.95*IC[2]*(IC[0]>0.1)*(IC[1]>0.1)+3.16*IC[3]*(IC[0]>0.1)*(IC[1]>0.1)*(IC[2]>0.1)+1.1*IC[4]*(IC[0]>0.1)*(IC[1]>0.1)*(IC[2]>0.1)*(IC[3]>0.1)+IC[5]*(IC[0]>0.1)*(IC[1]>0.1)*(IC[2]>0.1)*(IC[3]>0.1)*(IC[4]>0.1))");

//t1->SetAlias("mE2","IC[0]+IC[1]*(IC[0]>0.1)+2*IC[2]*(IC[0]>0.1)*(IC[1]>0.1)+2*IC[3]*(IC[0]>0.1)*(IC[1]>0.1)*(IC[2]>0.1)+IC[4]*(IC[0]>0.1)*(IC[1]>0.1)*(IC[2]>0.1)*(IC[3]>0.1)+IC[5]*(IC[0]>0.1)*(IC[1]>0.1)*(IC[2]>0.1)*(IC[3]>0.1)*(IC[4]>0.1)");
   t1->SetAlias("mdE21","(IC[0]+IC[1])*(IC[0]>1&&IC[1]>0.1)");
 t1->SetAlias("mE21","(IC[0]+IC[1]+2*IC[2]+2*IC[3]+IC[4]+IC[5])*(IC[0]>0.1)*(IC[1]>0.1)*(IC[2]>0.1)*(IC[3]>0.1)");
 t1->SetAlias("mD2","Path+5");
  t1->SetAlias("mV2","mD2/(mT2)");
  t1->SetAlias("mBeta2","mV2/29.9792");
  t1->SetAlias("mGamma2","1./sqrt(1.0-mBeta2*mBeta2)");
  
  t1->SetAlias("mM_Q2","Brho/3.105/mBeta2/mGamma2");
  t1->SetAlias("mM2","(mE2)/931.5016/(mGamma2-1.)");
  t1->SetAlias("mQ2","mM2/mM_Q2");
  t1->SetAlias("mMrec","int(mQ2+0.5)*mM_Q2");
  
  t1->SetAlias("Gamma_from_E","mE2/14./931.5016+1");
  t1->SetAlias("Beta_from_E","sqrt(1.-1./pow(Gamma_from_E,2))");
  t1->SetAlias("Brho_from_E","14./8*3.105*Beta_from_E*Gamma_from_E");

  t1->SetAlias("mEM","(3.7*IC[0]+IC[1]+IC[2]+IC[3]*(IC[2]>0.1)+18./55*IC[4]*(IC[3]>0.1&&IC[2]>0.1) + 7.*IC[5]*(IC[4]>0.1&&IC[3]>0.1&&IC[2]>0.1) )*(IC[0]>1&&IC[1]>0.1&&IC[2]>1&&IC[3]>1)");

}