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");
  
}

void chkef( ULong64_t low, ULong64_t high,   ULong64_t n ){
  cout<<endl;
  cout<< "=========================================================\n" ;
  cout<<"i... HOUR: "<<n<<endl; //"     "<<low<<" "<<high<<endl;
  string sep = "----------------------------\n" ;



  TTree *t1 = (TTree*) gROOT->FindObject("AD");

  TH1F * hc = new TH1F("hc","hc",2000,0.1,50);
  TH1F * ht = new TH1F("ht","ht",2000,0.1,50);


  TString cond=  Form("(((RTS<<32)+(VTS<<16)+LTS)>=%lld)&&(((RTS<<32)+(VTS<<16)+LTS)<%lld)",low,high);
  TString condad=Form("&&( (((RTS<<32)+(VTS<<16)+LTS)>=%lld)&&(((RTS<<32)+(VTS<<16)+LTS)<%lld) )",low,high);


  // cout<<cond<<endl<<condad<<endl;



      cout << "Forward (condition by before)" << endl; 
  for(int i =0; i<1;i++)
    {
      t1->Draw("T_FPMW_HF_C:T_FPMW_CATS2_C",cond.Data(),"goff"); 
      long long a=AD->GetSelectedRows();

      t1->Draw("T_FPMW_HF_C:T_FPMW_CATS2_C",Form("T_FPMW_CATS2_C<1 && T_FPMW_HF_C>2 && %s",cond.Data() ),"goff");
      long long b=AD->GetSelectedRows();

	cout << "CATS lost efficiency with resect to VAMOS" << 100.*b/a  << "%" << endl;

    }

  cout << sep ;




      cout << "Forward (condition by before)" << endl; 
  for(int i =0; i<4;i++)
    {
      t1->Draw(Form("IC[%d]>>ht",i+1),cond.Data() ,"goff");
      t1->Draw(Form("IC[%d]>>hc",i+1),Form("IC[%d]>0.1 %s",i,condad.Data() ),"goff");
      cout << "IC " << i << " " << hc->Integral()/ht->Integral()*100. << "%"<< endl;

    }

  cout << sep ;


  cout << "Backward  (condition by after)" << endl; 

  for(int i =0; i<4;i++)
    {
      t1->Draw(Form("IC[%d]>>ht",i),cond.Data() ,"goff");
      t1->Draw(Form("IC[%d]>>hc",i),Form("IC[%d]>0.1 %s",i+1,  condad.Data() ),"goff");
      cout << "IC " << i+1 << " " << hc->Integral()/ht->Integral()*100. << "%"<< endl;

    }
cout << sep ;

 cout << "DC  (condition by IC[0])" << endl; 

  for(int i =0; i<4;i++)
    {
      t1->Draw(Form("IC[0]>>ht"),cond.Data(),"goff");
      t1->Draw(Form("IC[0]>>hc"),Form("EWIRE_%d>100 %s",i, condad.Data() ),"goff");
      cout << "DC WIRE " << i << " " << hc->Integral()/ht->Integral()*100. << "%"<< endl;

    }
cout << sep ;

cout << "DC X  (condition by IC[0])" << endl; 

  for(int i =0; i<4;i++)
    {
      t1->Draw(Form("IC[0]>>ht"),cond.Data(),"goff");
      t1->Draw(Form("IC[0]>>hc"),Form("DC%d_X>-1500  %s",i, condad.Data() ),"goff");
      cout << "DCX" << i << " " << hc->Integral()/ht->Integral()*100. << "%"<< endl;

    }

cout << sep ;


cout << "Xf  (condition by IC[0])" << endl; 

      t1->Draw(Form("IC[0]>>ht"),cond.Data() ,"goff");
      t1->Draw(Form("IC[0]>>hc"),Form("Xf>-1500 %s", condad.Data() ),"goff");
      cout << "Xf " << hc->Integral()/ht->Integral()*100. << "%"<< endl;


cout << sep ;


 delete hc;
 delete ht;
}
//==============================================================end





/*
ORIGINAL NAME but  now the original code for VAMOS eff is above
 */
void EffCheck_fhjm()
{

  ULong64_t hour;
  ULong64_t LTSmin=-1;
  ULong64_t LTSmax=-1;
  ULong64_t LTS,VTS,RTS, tmin,tmax,tmp;
 
  cout<< "============================================================\n" ;
  cout<<"          "<<gFile->GetName()<<endl;
  cout<< "============================================================\n" ;


  cout<<"Looking for min and max time stamp... very slow"<<endl;
  TTree *tree= (TTree*)gDirectory->Get("AD");
  tree->SetBranchAddress("LTS",&LTS);
  tree->SetBranchAddress("VTS",&VTS);
  tree->SetBranchAddress("RTS",&RTS);
  
  
 for (int i=1;i<tree->GetEntries();i++){
   tree->GetEntry(i);
   tmp=(RTS<<32)+(VTS<<16)+LTS;
   if ((LTS>0)&&(tmp>1e+14)){  //================== kill some low values?
     if ((LTSmin>tmp) || (LTSmin==-1)){
       LTSmin=tmp;
     }
     if ((LTSmax<tmp) || (LTSmax==-1)){
       LTSmax=tmp;
     }
   }
 }


 //tree->Draw("(RTS<<32)+(VTS<<16)+LTS>>h(10000)","","");
 //  Long64_t LTSmin=AD->GetV1()[0];
 //  AD->Draw("RTS<<32+VTS<<16+LTS","LTS>0","goff", 1, AD->GetEntries()-1 );
 //  Long64_t LTSmax=AD->GetV1()[0];
 cout<<"min    "<<LTSmin<<" "<<endl;
 cout<<"max    "<<LTSmax<<" "<< endl;
 cout<<"delta  "<<LTSmax-LTSmin<<endl;
 cout<<"delta  "<<(LTSmax-LTSmin)/100/1e+6/3600/282162*4.3<<" hours  EMPIRICAL RUN LENGHT"<<endl;
 hour=143567947269965277/6.0775;
 cout<<"1hour=="<<hour<<endl<<endl<<endl;
 //==========================================================================================
//   MAIN LOOP OVER TIME STEPS
 ULong64_t n=0;
 int hourstep=2;
 for (ULong64_t i=LTSmin; i<LTSmax; i+=hour*hourstep){
   n=n+1;
   chkef( i, i+hour , n*hourstep);
 }

  return NULL;




}