check_aligned.C File Reference

Functions
void	a ()
void	check_aligned ()
void	check_diff ()
void	check_shrinkage ()

Function Documentation

a()

void a

(

)

{
  cshr_1->cd();   couples->Draw("(s2.eY-s1.eY)/(s2.eZ-s1.eZ):(s2.eX-s1.eX)/(s2.eZ-s1.eZ)");
  cshr_3->cd();   couples->Draw("eCHI2");
}

check_aligned()

void check_aligned

(

)

{
  check_diff();
  //  check_shrinkage();
}

check_diff()

void check_diff

(

)

{
  TCanvas *diff = new TCanvas("diff","diff");
  diff->Clear();
  diff->Divide(2,2);
  gStyle->SetOptFit(0001);
 
  diff_1->cd();   
  couples->Draw("s1.eTX-(s2.eX-s1.eX)/(s2.eZ-s1.eZ):(s2.eX-s1.eX)/(s2.eZ-s1.eZ)","pid2>-1&&eCHI2<1.5&&abs(s1.eTY)<.1&&abs(s2.eTY)<.1");
  diff_2->cd();   
  couples->Draw("s2.eTX-(s2.eX-s1.eX)/(s2.eZ-s1.eZ):(s2.eX-s1.eX)/(s2.eZ-s1.eZ)","pid2>-1&&eCHI2<1.5&&abs(s1.eTY)<.1&&abs(s2.eTY)<.1");
 
  diff_3->cd();   
  couples->Draw("eCHI2","pid2>-1&&abs(s1.eTY)<.1&&abs(s2.eTY)<.1");
  diff_4->cd();   
  couples->Draw("eCHI2P","pid2>-1&&abs(s1.eTY)<.1&&abs(s2.eTY)<.1");
}

check_shrinkage()

void check_shrinkage

(

)

{
  // this function check shrinkage and/or distance between linked planes 
  // Note: do not use s.* (linked segment parameters), because them could 
  // be different from the "base angle" calculated here directly
 
  TCanvas *cshr = new TCanvas("cshr","shrinkage");
  cshr->Clear();
  cshr->Divide(2,2);
  gStyle->SetOptFit(0001);
 
  cshr_1->cd();   
  couples->Draw("s1.eTX-(s2.eX-s1.eX)/(s2.eZ-s1.eZ):(s2.eX-s1.eX)/(s2.eZ-s1.eZ)","pid2>-1&&eCHI2<1.5","prof");
  htemp->Fit("pol1","w","",-.4,.4);
  float p0 = htemp->GetFunction("pol1")->GetParameter(0);
  float p1 = htemp->GetFunction("pol1")->GetParameter(1);
  char str[160]="";
  sprintf(str,"s1.eTX*(1-(%f))-(%f)-(s2.eX-s1.eX)/(s2.eZ-s1.eZ):(s2.eX-s1.eX)/(s2.eZ-s1.eZ)",p1,p0);
  cshr_3->cd();   couples->Draw(str,"pid2>-1&&eCHI2<1.5","prof");
  htemp->Fit("pol1","w","",-.4,.4);
 
  cshr_2->cd();   
  couples->Draw("s2.eTX-(s2.eX-s1.eX)/(s2.eZ-s1.eZ):(s2.eX-s1.eX)/(s2.eZ-s1.eZ)","pid2>-1&&eCHI2<1.5","prof");
  htemp->Fit("pol1","w","",-.4,.4);
  p0 = htemp->GetFunction("pol1")->GetParameter(0);
  p1 = htemp->GetFunction("pol1")->GetParameter(1);
  sprintf(str,"s2.eTX*(1-(%f))-(%f)-(s2.eX-s1.eX)/(s2.eZ-s1.eZ):(s2.eX-s1.eX)/(s2.eZ-s1.eZ)",p1,p0);
  cshr_4->cd();   couples->Draw(str,"pid2>-1&&eCHI2<1.5","prof");
  htemp->Fit("pol1","w","",-.4,.4);
}