ROOTコマンド集ページ

    {
      TFile *f = new TFile("ggh4mu.root");
      f->cd("ATLFAST");
      h51->Draw("Nmuo");
    }

http://www.icepp.s.u-tokyo.ac.jp/~sakamoto/work/Pukiwiki/pukiwiki.php?UserAlgorithm
にある。
-----------------------------------------------------
.hの中では以下のようにTH1F*で定義しておくようだ
-----------------------------------------------------
#include "TH1.h"
 private:  
  /** Histograms */
  TH1F* m_h_cluset;
  TH1F* m_h_cluse;
  TH1F* m_h_cluseta;
  TH1F* m_h_clusphi;
----------------------------------------------------
そしてcxx中では
----------------------------------------------------
#include "GaudiKernel/ITHistSvc.h" かなあ？？？
////////////////////////////////////////////////////////////////////////////////////
/// Initialize

  addBranch("ElectronEt",    m_aan_et);
  addBranch("ElectronPhi",   m_aan_phi);
  addBranch("ElectronE",    m_aan_e);
  addBranch("ElecPtRatio",   m_aan_elecetres);

  /// ROOT histograms ---------------------------------------

  m_h_cluset     = new TH1F("EMTopoCluster430_et","et cl",50,0,5.*GeV);
  sc = m_thistSvc->regHist("/AANT/EMTopoCluster430/clus_et",m_h_cluset);
　m_h_electron_eta_vs_phi = new TH2F( "electron_eta_vs_phi", "Electron Eta vs Phi", 100, -5.0, 5.0, 90, 0, 360 );
　m_thistSvc -> regHist( "/MyAnalysis/Electron/electron_eta_vs_phi", m_h_electron_eta_vs_phi );
 }
  m_aan_et->clear();
  m_aan_phi->clear();

/// Execute - on event by event

    m_h_cluset->Fill( (*elecItr)->et(), 1.);
    m_h_cluse->Fill( (*elecItr)->e(), 1.);
    m_h_cluseta->Fill( (*elecItr)->eta(), 1.);
    m_h_clusphi->Fill( (*elecItr)->phi(), 1.);
　　m_h_electron_eta_vs_phi -> Fill( elec -> eta( ), phi, 1.0 );

    /** fill Athena-Aware NTuple */
    m_aan_eta->push_back((*elecItr)->eta());
-----------------------------------------------------------------








root> c1->Print("test.ps")
   Info in : PostScript file test.ps has been created
rppt> c1->Update()     まだ意味不明
root> TPostScript zzz("test-zzz.ps"
)	ページを開く
root [17] h51->Draw("Nmuo")
root [18] c1->Update()     
root [19] h51->Draw("Pxmuo")
root [20] c1->Update()      
root [21] h51->Draw("Pymuo") 
root [22] c1->Update() 
Info in : PostScript file test.ps has been created

画面分割法

root [13] TFile *f = new TFile("ggh4mu.root")  fは変数名、 
root [14] f->cd("ATLFAST")
(Bool_t)1 
root [15] TCanvas *c1 = new TCanvas("sakamoto","c1 title");
  sakamotoはリソース名、c1はプログラム上の変数名、c1 titlはタイトル。
root [16] c1->Divide(3,2)
root [17] c1->Clear();
root [18] c1->Divide(3,2);
root [19] c1->cd(1);
root [20] h51->Draw("Nmuo");
root [21] h51->Draw("Pxmuo");
root [22] c1->cd(2);         
root [23] h51->Draw("Eemuo");
root [24] c1->cd(3);         
root [25] h51->Draw("Kfmuo");



 STL

std::vector<Particle*> &AAParticle型の可変長配列(vextor)のAAをSTLのコンテナvectorで設定する。AA.push_back(a), AA.size(), AA.......などがある。
（蔵重の中級入門＠岡山）


 CLASS

class B : public A { } クラスAを基底として派生クラスBを定義する。Aのpublic memberはすべてBのpublicメンバになる。


 その他

./ZZ < input.script >& loginput.script入力でZZを実行し、すべての出力をlogに出す。
&& || ! .and.   .or.   not.
double xdouble precision float型の変数 x を定義。"float" "double" "long double"がある。
const double pi=3.14定数識別子piの定義
int main(void)void識別子は引数がないことを示す。Cの流れで引数のvoidは省略できる。
void f(int);関数fは整数引数をもつが値を返さない関数なのでvoidと書く。



p.eta(), p.phi(),  p.pt(),  
Particleクラス
 eta(), phi(), pt(), pId().islsoMuon(),  pType().pdgid(),  


---------------------------------------------------
2003.7.22

9:11    mkdir ana
cd ana
cp /atlas/tutorial/root/atlfast/codes.tar.gz .
tar zxvf codes.tar.gz
cd samples/
source atlas/tutorial/root/etc/atlj.csh  
echo $ATLJ_TOP_DIR

cd code01
gmake    
./ZZ
現在６個のコマンドのみ  .q, exit, quitなどで出る

code02がヒストグラム
文字列をptype unitで読み込んでいる

activate allで全ておんになり時間がかかる。

ーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーー
2003.7.22  
ここでZZー＞Hを組むプログラムを入れる。隣の蔵重がたすけてくれた。

ZZ.h ,  ZZ.cxx  を変える。
make clean  が必要みたいだ。
gmake
./ZZ < signal.script > & log
root result.root
root [1] TBrowser tw
で見えるはず。
孫のmuonが重なっていないかcheckしたら39 eventsになるはず。

ーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーー
バックグランドファイルは大きいので(１００ファイル)、興味あるイベント
のみをのこして新しいファイルをつくる。code05

---------------------------------------------------------------
バックグランドを重ねるーーー＞ code04.ansをコピーする。




ーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーー
長野
mkdir tgc
cd tgc
rootlogon.C


{
  if (!TClassTable::GetDict("TgcDataObject")) {
    gSystem->Load("libTgcDataObject.so");
  }

  TFile *f = new TFile("tgc_2.root");
  f->ls();

  TgcDataObject *data = new TgcDataObject();

  TTree *tree = (TTree*)gROOT->FindObject("tree");
  TBranch *branch = tree->GetBranch("data");
  branch->SetAddress(&data);
}


setenv LD_LIBRARY_PATH ${LD_LIBRARY_PATH}:
./ed
TGC H8 ED > next
------------------------------------------
    Event:        26
------------------------------------------
 L1ID:    27  BCID:   114  nSLB:     4
------------------------------------------

/usr/local/Acrobat5/bin/acroread RootTutorial.pdf &
に解説あり。

hits.C を変えて１イベント毎にプリントしてみ







Send your comments on this page to 
kondo@kekvax.kek.jp