Alignment/OfflineValidation/scripts/Chi2_mu.C - master - githubmirror/cmssw

Alignment/OfflineValidation/scripts/Chi2_mu.C
146 строк4 KB
Giulio Eulisse
Snapshot of CMSSW_6_2_0_pre8.
26 июн 2013, 18:12
26 июн 2013, 18:12214ab73

class Plots{

  public:
  Plots();

  void Legend(TString histoname1,TString histoname2,TString histoname3, TH1F *histo1, TH1F *histo2, TH1F *histo3);  void Save();
  
  float convert(float num);

};

void Plots::Plots()
{

gROOT->Reset();
gROOT->Clear();

gStyle->SetNdivisions(10);
gStyle->SetCanvasBorderMode(0); 
gStyle->SetPadBorderMode(1);
gStyle->SetOptTitle(1);
gStyle->SetStatFont(42);
gStyle->SetCanvasColor(10);
gStyle->SetPadColor(0);
gStyle->SetTitleFont(62,"xy");
gStyle->SetLabelFont(62,"xy");
gStyle->SetTitleFontSize(0.07);
gStyle->SetTitleSize(0.046,"xy");
gStyle->SetLabelSize(0.052,"xy");
// gStyle->SetTitleFillColor(0);
gStyle->SetHistFillStyle(1001);
gStyle->SetHistFillColor(0);
gStyle->SetHistLineStyle(1);
gStyle->SetHistLineWidth(2);
gStyle->SetHistLineColor(2);
gStyle->SetTitleXOffset(1.1);
gStyle->SetTitleYOffset(1.15);
gStyle->SetOptStat(1110);
gStyle->SetOptStat(kFALSE);
gStyle->SetOptFit(0111);
gStyle->SetStatH(0.1); 

TCanvas *c1 = new TCanvas("c1","c1",129,17,926,703);
c1->SetBorderSize(2);
c1->SetFrameFillColor(0);
c1->SetLogy(0);
c1->cd(); 

TFile *f[4];
TTree *MyTree[4];


f[0]= new TFile("../../singlemu_310607/Misalignment_scenarioIdeal_singlemu131.root");
MyTree[0]=Tracks;
 
f[1]=new TFile("Misalignment_SurveyLASOnlyScenario_refitter_zmumu_singlemuSurveyLASCosmics.root");
MyTree[1]=Tracks;
 
f[2]=new TFile("../../singlemu_310607/Misalignment100.root");
MyTree[2]=Tracks;
 
f[3]=new TFile("../../singlemu_310607/Misalignment1000.root");
MyTree[3]=Tracks;

////&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

//&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
/// EFFICIENCIES VS ETA ALIGNED
//&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
char histoname[128];
char name[128];
TH1F *Chi2[4];

for(int i=0; i<4; i++){
  sprintf(name,"Chi2[%d]",i);
  Chi2[i] = new TH1F(name,name,100,0,15);
  sprintf(histoname,"Chi2[%d]",i);
  MyTree[i]->Project(histoname,"CHISQ","eff==1 && TrackID==13");

  cout << "Entries " << Chi2[i]->GetEntries() <<endl;
  Chi2[i]->Scale(1/Chi2[i]->GetEntries());
  Chi2[i]->SetTitle("Track Fit #chi^{2}/ndf");
  Chi2[i]->SetXTitle("#chi^{2}/ndf");
  Chi2[i]->SetYTitle("arb. units");
  Chi2[i]->SetLineColor(i+2);
  Chi2[i]->SetLineStyle(i+1);
  Chi2[i]->SetLineWidth(i+2);
  //if (i==0) Chi2[i]->Draw();
  //else Chi2[i]->Draw("same");
  c1->Update();
}

Chi2[0]->Draw("");
//Chi2[0]->Draw("same");
Chi2[1]->Draw("same");
Chi2[2]->Draw("same");
Chi2[3]->Draw("same");


Legend("Chi2[0]","Chi2[1]","Chi2[2]","Chi2[3]",Chi2[0],Chi2[1],Chi2[2],Chi2[3]);

c1->SaveAs("Chi2_mu.eps");
c1->SaveAs("Chi2_mu.gif");
gROOT->Reset();
gROOT->Clear();

delete c1;

}

void Plots::Legend(TString histoname1,TString histoname2,TString histoname3, TString histoname4, TH1F *histo1, TH1F *histo2, TH1F *histo3, TH1F *histo4)
{

TLegend *leg = new TLegend(0.25,0.7,0.85,0.85.); 
leg->SetTextAlign(32);
leg->SetTextColor(1);
leg->SetTextSize(0.025);

char  label[128];
sprintf(label,"perfect alignment;      mean = %1.3f, RMS = %1.3f",convert(histo1->GetMean()),convert(histo1->GetRMS()));
leg->AddEntry(histoname1, label, "l");
sprintf(label,"SuveyLASCosmics alignment; mean = %1.3f, RMS = %1.3f",convert(histo2->GetMean()),convert(histo2->GetRMS()));
leg->AddEntry(histoname2, label, "l");
sprintf(label,"100 pb-1 alignment;  mean = %1.3f, RMS = %1.3f",convert(histo3->GetMean()),convert(histo3->GetRMS()));
leg->AddEntry(histoname3, label, "l");
sprintf(label,"1000 pb-1 alignment;  mean = %1.3f, RMS = %1.3f",convert(histo4->GetMean()),convert(histo4->GetRMS()));
leg->AddEntry(histoname4, label, "l");
leg->Draw();

}

float Plots::convert(float num){
  int mean1 = num;
  float res = num - mean1;
  int res2 = res*1000;
  float res3 = res2*0.001;
  float mean2 = mean1 + res3;
                                                                                                                  
  float res4 = res - res3;
  int res5 = res4*10000;
                                                                                                                  
  if(res5>5)
    mean2 = mean2 + 0.001;
                                                                                                                  
  return mean2;
}