Eine Selektion starten

Normalerweise möchte man nicht alle Ereignisse eines ROOT tuples anzeigen, sondern nur bestimmte Ereignisse mit speziellen Eigenschaften anzeigen. Folgende Programmzeilen füllen einen TLorentzVector mit Spurinformationen von Myonen.

• Fügen Sie folgende Zeilen am Anfang von MyClass.C ein:
  #include <TLorentzVector.h>
  #include <TMath.h>
  

• Definieren Sie mehrere Histogramme, TLorentzVector und temporäre Variablen mit folgenden Zeilen vor der for-Schleife: TH1D *h_pt_mu1 = new TH1D("ptMuon1","pt of Muon 1", 50,0.,100.*1000.);
  TH1D *h_pt_mu2 = new TH1D("ptMuon2","pt of Muon 2",50,0.,100.*1000.);
  TLorentzVector m1(0,0,0,0), m2(0,0,0,0), m(0,0,0,0);
  

• Fügen Sie folgende Zeilen in der for-Schleife ein:

  m1.SetPxPyPzE(0,0,0,0);
  m2.SetPxPyPzE(0,0,0,0);
  m.SetPxPyPzE(0,0,0,0);
  

  for (int i=0; i<Mu_N;i++) {
  m.SetPxPyPzE((*Mu_px)[i],(*Mu_py)[i],(*Mu_pz)[i],(*Mu_E)[i]);
  

  if (m.Pt()>m1.Pt()) {
  m2 = m1;
  m1 = m;
  }
  else if((m.Pt()>m2.Pt()) && (m.Pt() < m1.Pt())) {
  m2 = m;
  }
  }
  

  if ((m1.Pt()>5.*1000.) &&
  (m2.Pt()>5.*1000.)) {
  h_pt_mu1->Fill(m1.Pt());
  h_pt_mu2->Fill(m2.Pt());
  }
  

  Diese Zeilen füllen die Histogramme mit den Myon Transversalimpulsen pT und selektieren nur Ereignisse mit zwei Myonen und pT>5 GeV.

• Ein Canvas kann auf in mehrer TPad unterteilt werden, so daß mehrere Histogramme gleichzeitig angezeigt werden können: TStyle *style = gROOT->GetStyle("Plain");
  gROOT->SetStyle("Plain");
  TCanvas *can1=new TCanvas("can1","Muons",0,0,800,800);
  can1->Divide(2,2);
  can1->cd(1);
  h_n_mu->Draw();
  can1->cd(2);
  h_pt_mu1->Draw();
  can1->cd(3);
  h_pt_mu2->Draw();