void NumEvLimAlignConsts(){
 TString pathG="/home/karavdin/datastorage/AlignmentLMDpixel/";
 TString pathGb=pathG+"/RECdataAUG2014/BOX/mom_15/DOF110001/";
 // TString pathGb=pathG+"/RECdataAUG2014/BOX/mom_15/DOF111111/";
 const int nAS1=5; // number of misalignment consts samples
 // const int nAS1=10; // number of misalignment consts samples
 const int nAS2=10; // number of misalignment consts samples
 const int nS=100; // number of data samples for each misalignment consts sample
 // const int nParDt=11;
 // double v_Dt[nParDt]={0, 10, 50, 100,200,300,400,500,600,700,800};
 const int nParDt=8;
 double v_Dt[nParDt]={0, 10, 50, 100,200,300,400,500};
 //  const int nParDa=5;
 //  double v_Da[nParDa]={0,0.2,1,3,5};
 //  TString strv_Da[nParDa]={"0","0_2","1","3","5"};
 // TString strv_Da_leg[nParDa]={"0","0.2","1","3","5"};
 
 const int nParDa=4;
 double v_Da[nParDa]={0,1,3,5};
 TString strv_Da[nParDa]={"0","1","3","5"};
 TString strv_Da_leg[nParDa]={"0","1","3","5"};

 // const int nParDa=1;
 // double v_Da[nParDa]={5};
 // TString strv_Da[nParDa]={"5"};
 // TString strv_Da_leg[nParDa]={"5"};

 double NevForKnossos[nParDa][nParDt];
 const double trans_lim=0.5*15.;//[mum] 15 GeV
 const double rot_lim=0.5*1.;//[mrad] 15 GeV

  for(int iDt=0;iDt<nParDt;iDt++){
    double tr_sc = v_Dt[iDt];
    for(int iDa=0;iDa<nParDa;iDa++){
      double rt_sc= v_Da[iDa];
      TString strrt_sc= strv_Da[iDa];
      //Read data
      int nAS = nAS1;
      if(tr_sc>599.) nAS = nAS2;
      else{
	if(tr_sc>499. && rt_sc>1) 
	  nAS = nAS2;
	else{
	  if(rt_sc>3)
	    nAS = nAS2;
	}
      }
      TString inmane= pathGb+"/ResultsSamplesSummary_Av_";
      inmane += nAS;
      inmane += "AlignSamples_";
      inmane += nS;
      inmane += "DataSamples_";
      inmane+=tr_sc;
      inmane+="um_";
      inmane+=strrt_sc;
      inmane+="mrad.root";
      TFile *fin = new TFile(inmane,"READ");
      TGraphErrors *misX_a = (TGraphErrors *)fin->Get("dX_a");
      TGraphErrors *misY_a = (TGraphErrors *)fin->Get("dY_a");
      TGraphErrors *misZ_a = (TGraphErrors *)fin->Get("dZ_a");
      TGraphErrors *misA_a = (TGraphErrors *)fin->Get("dA_a");
      TGraphErrors *misB_a = (TGraphErrors *)fin->Get("dB_a");
      TGraphErrors *misG_a = (TGraphErrors *)fin->Get("dG_a");
      int Npoints = misX_a->GetN();
      Double_t *nval =  misX_a->GetX();
      Double_t *xval =  misX_a->GetY();
      Double_t *yval =  misY_a->GetY();
      Double_t *gval =  misG_a->GetY();
      NevForKnossos[iDa][iDt]=-1e-9;
      for(int i=1;i<Npoints;i++){
	//	cout<<"["<<v_Dt[iDt]<<"mum]["<<v_Da[iDa]<<"]="<<gval[i]<<endl;
	if(xval[i]<trans_lim && yval[i]<trans_lim && gval[i]<trans_lim){
	//	if((xval[i]<trans_lim || fabs(xval[i]-xval[i-1])<0.1) && (yval[i]<trans_lim || fabs(yval[i]-yval[i-1])<0.1) && (gval[i]<trans_lim|| fabs(gval[i]-gval[i-1])<0.01)){
	 NevForKnossos[iDa][iDt]=nval[i];
	 //	 cout<<"NevForKnossos = "<<NevForKnossos[iDa][iDt]<<endl;
	}
	if(NevForKnossos[iDa][iDt]>1) break;
      }

    } // loop over dA  [rotation]
  }// loop over dT [translation]

  //Draw resutls
  const int colors_a[5] = {kViolet-6,kOrange+7,kGreen+3,kRed+1,kBlue+1};
  const int colors[5] = {28,kAzure+2,46,8,kGreen+8};
  TLegend *leg_a = new TLegend(0.73,0.15,0.87,0.37);
  // TLegend *leg_a = new TLegend(0.13,0.65,0.27,0.87);
  leg_a->SetFillColor(0);
  TMultiGraph *mgr = new TMultiGraph();
  TGraphErrors *gr[nParDa];   

 
  for(int iDa=0;iDa<nParDa;iDa++){
    gr[iDa] = new TGraphErrors(nParDt,v_Dt,NevForKnossos[iDa],0,0);
    gr[iDa]->SetMarkerSize(2.5);
    //    gr[iDa]->SetMarkerStyle(29+iDa);
    gr[iDa]->SetMarkerStyle(22+iDa);
    gr[iDa]->SetMarkerColor(colors_a[iDa]);
    mgr->Add(gr[iDa]);

    TString dRName = "#Delta_{r} = ";
    dRName +=strv_Da_leg[iDa];
    leg_a->AddEntry(gr[iDa],dRName,"lep");
  } // loop over dA  [rotation]
  //mgr_G->Draw("AP");
  //Save plots
  TString resname= pathGb+"/ResultsMinTrksKnossos_";
  resname += nAS;
  resname += "AlignSamples_";
  resname += nS;
  resname += "DataSamples";
  TString resname_pdf = resname+".pdf";
  TString resname_pdf_o = resname_pdf+"(";
  TString resname_pdf_c = resname_pdf+")";
 

  TCanvas c1;
  c1.SetLogy();
  mgr->Draw("AP");
  mgr->GetXaxis()->SetTitle("N trks/sector");
  mgr->GetYaxis()->SetTitle("#Delta_{x}, #mum");
  leg_a->Draw();
  mgr->SetName("Nmin");
  mgr->SetMinimum(10);
  mgr->SetMaximum(100000);
  c1.Print(resname_pdf);

  TString out = resname+".root";
  TFile *f = new TFile(out,"RECREATE");
  mgr->Write();
  leg_a->SetName("leg");
  leg_a->Write();
  f->Write();
  f->Close();
}