//
// plot DP dimuon data
//
//  03/20/2018  MXL
//
//  1 = "mu+" and 2 ="mu-"
//
//  TriggerBits:
//   FPGA1 = 
///////////////////////
{
  TCut XX1 = "chisq1<20 && chisq2<20 && nHits1>0 && nHits2>0 && abs(x1_st1)<100 && abs(x2_st1)<100";  // single muon cuts

  TCut XX2 = "abs((1.6-y1)/ty1)<1000 & abs((1.6-y2)/ty2)<1000 && tx1 < 100 && tx2 < 100 ";  // dimuons cuts


  
   TCanvas * c1 = new TCanvas("c1", "DP: Single Tracks", 10, 20, 600, 800);

   TCanvas * c2 = new TCanvas("c2", "DP: Dimuons", 500, 100, 600, 800);

   save->SetAlias("vz1","(1.6-y1)/ty1");
   save->SetAlias("vz2","(1.6-y2)/ty2");
   
   save->SetAlias("X1","x1_st1+tx1_st1*797");  // (X,Y) at DP Sta-1, z= 797cm
   save->SetAlias("Y1","y1+ty1*797");  // (X,Y) at DP Sta-1, z= 797cm

   save->SetAlias("X2","x2_st1+tx_st12*797");  // (X,Y) at DP Sta-1, z= 797cm
   save->SetAlias("Y2","y2+ty2*797");  // (X,Y) at DP Sta-1, z= 797cm
   
   
   //--- single muons
   
   c1->cd();

   c1->Divide(3,4);

   c1->cd(1);
   save->Draw("runID",XX1);

   c1->cd(2);
   save->Draw("spillID",XX1);

   c1->cd(3);
   save->Draw("eventID",XX1);

   c1->cd(4);
   save->Draw("nTracks",XX1);

   c1->cd(5);
   save->Draw("nHits1",XX1);

   c1->cd(6);
   save->Draw("nHits2",XX1);

   c1->cd(7);
   save->Draw("chisq1",XX1);

   c1->cd(8);
   save->Draw("chisq2",XX1);

   c1->cd(9);  
   save->Draw("pz1",XX1);
   save->Draw("pz2",XX1,"same");

   c1->cd(10)->SetLogy();
   save->Draw("triggerBits",XX1);

   c1->cd(11);
   save->Draw("x1_st1",XX1);
   save->Draw("x2_st1",XX1,"same");
   

   c1->cd(12);
   save->Draw("tx1_st1",XX1);
   save->Draw("tx2_st1",XX1,"same");

   
   //--- c2 dimuons
   c2->cd();
   c2->Divide(3,4);
   
   c2->cd(1)->SetLogy();
   save->Draw("y1",XX1+XX2);

   c2->cd(2);
   save->Draw("ty1",XX1+XX2);

   c2->cd(3);
   save->Draw("vz1",XX1+XX2);

   c2->cd(4)->SetLogy();
   save->Draw("y2",XX1+XX2);
   
   c2->cd(5);
   save->Draw("ty2",XX1+XX2);

   c2->cd(6);
   save->Draw("vz2",XX1+XX2);

	      
   c2->cd(7);
   save->Draw("vz1 - vz2",XX1+XX2);  // delta V_Z
	      
   c2->cd(8);
   save->Draw("vz1:vz2", XX1+XX2,"colz");
	      
   c2->cd(9);
   save->Draw("y2:x2_st1", XX1+XX2,"colz");

   c2->cd(10);
   save->Draw("pz1:pz2", XX1+XX2,"colz");

   c2->cd(11);
   save->Draw("Y1:X1", XX1+XX2,"colz");
   
   c2->cd(12);
   save->Draw("Y2:X2", XX1+XX2,"colz");
     
}