void LumiLinCheck() {

  TCanvas *c1 = new TCanvas("lc1","LinCheck",0,0,600,600);
  

  P->Draw("asym_bcm3","reg.ok_cut");
  Double_t Aqrms = htemp->GetRMS();
  char linetxt[50];
  sprintf(linetxt,"reg.ok_cut&&abs(asym_bcm3)<%d",4*Aqrms);
  TString ctstr = linetxt;
  P->SetAlias("cut1",ctstr);
  cout << P->GetAlias("cut1") << endl;

  c1->Clear();
  c1->Divide(2,2);


  TString plotstr;
  for(Int_t ilum=0; ilum<4; ilum++) {
    c1->cd(ilum+1);
    P->SetLineColor(1);
    sprintf(linetxt,"reg_asym_n_blumi%1d:asym_bcm3",(ilum+1)*2);
    plotstr = linetxt;
    cout << plotstr << endl;
    P->Draw(plotstr,"cut1","prof");
  }
  c1->cd();

  TCanvas *c2 = new TCanvas("lc2","LinCheck2",0,0,600,600);
  c2->Divide(2,2);
  TString plotstr;
  for(Int_t ilum=0; ilum<4; ilum++) {
    c2->cd(ilum+1);
    P->SetLineColor(1);
    sprintf(linetxt,"reg_asym_n_blumi%1d:asym_bcm3",(ilum+1)*2);
    plotstr = linetxt;
    cout << plotstr << endl;
    P->Draw(plotstr,"cut1","prof");
    P->SetLineColor(2);
    P->Draw(plotstr,"cut1&&evt_pairsynch[0]==0","prof same");
    P->SetLineColor(4);
    P->Draw(plotstr,"cut1&&evt_pairsynch[0]==1","prof same");
  }
  c2->cd();
  c1->Update();
  c2->Update();

  c1->Print("LumLinCheck1.png");
  c2->Print("LumLinCheck2.png");
}