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","535.1-2*T_FPMW_CATS2_C-1.2*(MWTVN==11)-0.3*(MWTVN==12)-0.5*(MWTVN==10)-0.7*(MWTVN==9)-2.8*(MWTVN==8)+0.2*(MWTVN==7)-0.4*(MWTVN==5)-1.0*(MWTVN==4)-3.0*(MWTVN==3)-3.0*(MWTVN==2)-8.0*(MWTVN==1)");
 //t1->SetAlias("mT2","mT1-TMWT1_MWT2_C+31");
  //  t1->SetAlias("mT2","mT1-4.83");

  //t1->SetAlias("mD2","D-10.5/cos(TP_Theta/1000.)/cos(TP_Phi/1000.)");
  t1->SetAlias("mE2","1.42*(IC[0]+IC[1]+IC[2]+IC[3]+0.3*IC[4])*(IC[0]>0.1&&IC[1]>0.1&&IC[2]>0.1&&IC[3]>0.1)");
  t1->SetAlias("mdE2","(IC[0]+IC[1])*(IC[0]>0.1&&IC[1]>0.1)");
  //  t1->SetAlias("mE1","(IC[1]+IC[2]+IC[3])*(IC[1]>1&&IC[2]>1&&IC[3]>1)");
  t1->SetAlias("mD2","Path");
  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.107/mV2*29.9792");
  t1->SetAlias("mM2","(mE2)/931.5016/(mGamma2-1.)");
  t1->SetAlias("mQ2","mM2/mM_Q2");
  t1->SetAlias("mMrec","int(mQ2+0.5)*mM_Q2");
  
}