FairRoot/PandaRoot
Functions | Variables
plot_acceptance_stt.C File Reference

Go to the source code of this file.

Functions

gStyle SetOptStat (kFALSE)
 
gStyle SetPalette (1)
 
c1 Divide (1, 3)
 
c1 cd (1)
 
hMcSttP SetTitle ("Geometrical coverage for the momentum")
 
hMcSttP GetXaxis () -> SetTitle("Momentum (GeV)")
 
hMcSttP Draw ()
 
hMcSttP SetLineColor (kRed)
 
hMcSttP Draw ("same")
 
hMcMvdP SetLineColor (8)
 
hMcGemP SetLineColor (6)
 
leg SetFillColor (kWhite)
 
leg AddEntry (hMcP,"ALL","l")
 
leg AddEntry (hMcSttP,"STT","l")
 
leg AddEntry (hMcMvdP,"MVD","l")
 
leg AddEntry (hMcGemP,"GEM","l")
 
c1 cd (2)
 
hMcSttTheta SetTitle ("Geometrical coverage for the theta angle")
 
leg AddEntry (hMcTheta,"ALL","l")
 
leg AddEntry (hMcSttTheta,"STT","l")
 
leg AddEntry (hMcMvdTheta,"MVD","l")
 
leg AddEntry (hMcGemTheta,"GEM","l")
 
c1 cd (3)
 
hMcSttPhi SetTitle ("Geometrical coverage for the phi angle")
 
leg AddEntry (hMcPhi,"ALL","l")
 
leg AddEntry (hMcSttPhi,"STT","l")
 
leg AddEntry (hMcMvdPhi,"MVD","l")
 
leg AddEntry (hMcGemPhi,"GEM","l")
 
nt Draw ("mc_p>>hMcP(100,0,5)")
 
nt Draw ("mc_p>>hMcAccSttP(100,0,5)","mc_stt>0","E1")
 
hMcAccSttP Divide (hMcP)
 
hMcAccSttP SetTitle ("Geometrical Acceptance")
 
hMcAccSttP GetYaxis () -> SetTitle("Acceptance")
 
nt Draw ("mc_theta>>hMcTheta(50,0,180)")
 
nt Draw ("mc_theta>>hMcAccSttTheta(50,0,180)","mc_stt>0","E1")
 
hMcAccSttTheta Divide (hMcTheta)
 
nt Draw ("mc_phi>>hMcPhi(50,-200,200)")
 
nt Draw ("mc_phi>>hMcAccSttPhi(50,-200,200)","mc_stt>0","E1")
 
hMcAccSttPhi Divide (hMcPhi)
 
c5 cd ()
 
ntEvt Draw ("acc_stt>>f1")
 
f1 SetTitle ("Geometrical Global Acceptance")
 
f1 Scale (1/numero)
 
nt Draw ("mc_p>>hRecoEffP(100,0,5)","mult>0","E1")
 
hRecoEffP SetTitle ("Pattern Recognition Efficiency")
 
nt Draw ("mc_theta>>hRecoEffTheta(50,0,180)","mult>0","E1")
 
nt Draw ("mc_phi>>hRecoEffPhi(50,-200,200)","mult>0","E1")
 
nt Draw ("mc_p>>hRecoEffCtP(100,0,5)","mult>0&&(stt>0||tpc>0)","E1")
 
hRecoEffCtP SetTitle ("Pattern Recognition Efficiency Only Central Tracker")
 
nt Draw ("mc_theta>>hRecoEffCtTheta(50,0,180)","mult>0&&(stt>0||tpc>0)","E1")
 
nt Draw ("mc_phi>>hRecoEffCtPhi(50,-200,200)","mult>0&&(stt>0||tpc>0)","E1")
 
ntEvt Draw ("cand>>f1")
 
f1 SetTitle ("Number of candidates")
 
ntEvt Draw ("eff>>f1")
 
f1 SetTitle ("Number of reconstructed primary tracks STT+MVD+GEM")
 
ntEvt Draw ("effct>>f1")
 
f1 SetTitle ("Number of reconstructed primary tracks Only Central Tracker")
 
c9 Divide (3, 2)
 
nt Draw ("(mc_p-p):mc_p>>hresp_p(100,0,5,100,-0.5,0.5)","mult>0 && (stt>0||tpc>0)","colz")
 
hresp_p SetTitle ("Resolution")
 
nt Draw ("(mc_theta-theta):mc_theta>>hrestheta_theta(100,0,100,100,-2,2)","mult>0 && (stt>0||tpc>0)","colz")
 
nt Draw ("(mc_p-p):mc_theta>>hresp_theta(100,0,100,100,-0.5,0.5)","mult>0 && (stt>0||tpc>0)","colz")
 
c9 cd (4)
 
nt Draw ("(mc_theta-theta):mc_p>>hrestheta_p(100,0,5,100,-2,2)","mult>0 && (stt>0||tpc>0)","colz")
 
c9 cd (5)
 
nt Draw ("(mc_phi-phi):mc_phi>>hresphi_phi(100,-200,200,100,-2,2)","mult>0 && (stt>0||tpc>0)","colz")
 
c9 Update ()
 
nt Draw ("p>>hist","mult>0 && (stt>0 || tpc>0)")
 
nt Draw ("theta>>hist2","mult>0 && (stt>0 || tpc>0)")
 
hist3 SetStats (kTRUE)
 
nt Draw ("(mc_p-p)/mc_p>>hist3","mult>0 && (stt>0 || tpc>0)")
 
hist3 Fit (f10,"R")
 

Variables

TFile * f = TFile::Open(inFile)
 
TCanvas * c1 =new TCanvas("c1","",800,800)
 
 leg =new TLegend(0.6,0.7,0.89,0.89)
 
TCanvas * c4 =new TCanvas("c4","",800,800)
 
TCanvas * c5 =new TCanvas("c5","",800,800)
 
float numero =ntEvt->GetEntries()
 
TH1F * f1 =new TH1F("f1","",100,0,5)
 
TCanvas * c6 =new TCanvas("c6","",800,800)
 
TCanvas * c14 =new TCanvas("c14","",800,800)
 
TCanvas * c7 =new TCanvas("c7","",800,800)
 
TCanvas * c8 =new TCanvas("c8","",800,800)
 
TCanvas * c11 =new TCanvas("c11","",800,800)
 
TCanvas * c9 =new TCanvas("c9","",800,800)
 
TCanvas * c20 =new TCanvas("c20","",800,800)
 
TH1F * hist =new TH1F("hist","Total Momentum distribution",100,0,10)
 
TCanvas * c21 =new TCanvas("c21","",800,800)
 
TH1F * hist2 =new TH1F("hist2","#theta Angle Distribution",100,0,180)
 
TCanvas * c22 =new TCanvas("c22","",800,800)
 
TH1F * hist3 =new TH1F("hist3","(Reconstructed Momentum - MC Momentum )/MC Momentum",100,-1,1)
 
TF1 * f10 = new TF1("f10","gaus",-0.09,0.09)
 

Function Documentation

leg AddEntry ( hMcP  ,
"ALL"  ,
"l"   
)
leg AddEntry ( hMcSttP  ,
"STT"  ,
"l"   
)
leg AddEntry ( hMcMvdP  ,
"MVD"  ,
"l"   
)
leg AddEntry ( hMcGemP  ,
"GEM"  ,
"l"   
)
leg AddEntry ( hMcTheta  ,
"ALL"  ,
"l"   
)
leg AddEntry ( hMcSttTheta  ,
"STT"  ,
"l"   
)
leg AddEntry ( hMcMvdTheta  ,
"MVD"  ,
"l"   
)
leg AddEntry ( hMcGemTheta  ,
"GEM"  ,
"l"   
)
leg AddEntry ( hMcPhi  ,
"ALL"  ,
"l"   
)
leg AddEntry ( hMcSttPhi  ,
"STT"  ,
"l"   
)
leg AddEntry ( hMcMvdPhi  ,
"MVD"  ,
"l"   
)
leg AddEntry ( hMcGemPhi  ,
"GEM"  ,
"l"   
)
c1 cd ( )
c1 cd ( )
c1 cd ( )
c5 cd ( )
c9 cd ( )
c9 cd ( )
c1 Divide ( ,
 
)
hRecoEffP Divide ( hMcP  )
hRecoEffTheta Divide ( hMcTheta  )
hRecoEffPhi Divide ( hMcPhi  )
c9 Divide ( ,
 
)
hMcSttP Draw ( )

Definition at line 64 of file NHitsPerEvent.C.

64  {
65 
66  std::cout << "mvdHitArray: " << mvdHitArray->GetEntries() << " projectedHitArray " << mvdProjectedArray->GetEntries() << std::endl;
67 
68 // gHitPoints->Clear();
71  //if (gHitPoints->GetN() != x.size()) gHitPoints->Expand(x.size());
72 
73  std::cout << "NHits: " << gHitPoints->GetN() << std::endl;
74  gHitPoints->Draw("apl");
75  gProjections->Draw("pl");
76 }
void ExtractData(TClonesArray *array, TGraph *g)
Definition: NHitsPerEvent.C:78
TGraph * gHitPoints
Definition: NHitsPerEvent.C:15
TClonesArray * mvdHitArray
Definition: NHitsPerEvent.C:12
TClonesArray * mvdProjectedArray
Definition: NHitsPerEvent.C:13
TGraph * gProjections
Definition: NHitsPerEvent.C:16
hMcSttP Draw ( "same"  )
nt Draw ( "mc_p>>hMcP(100,0,5)"  )
nt Draw ( "mc_p>>hMcAccSttP(100,0,5)"  ,
mc_stt,
0"  ,
"E1"   
)
nt Draw ( "mc_theta>>hMcTheta(50,0,180)"  )
nt Draw ( "mc_theta>>hMcAccSttTheta(50,0,180)"  ,
mc_stt,
0"  ,
"E1"   
)
nt Draw ( "mc_phi>>hMcPhi(50,-200,200)"  )
nt Draw ( "mc_phi>>hMcAccSttPhi(50,-200,200)"  ,
mc_stt,
0"  ,
"E1"   
)
ntEvt Draw ( "acc_stt>>f1 )
nt Draw ( "mc_p>>hRecoEffP(100,0,5)"  ,
mult,
0"  ,
"E1"   
)
nt Draw ( "mc_theta>>hRecoEffTheta(50,0,180)"  ,
mult,
0"  ,
"E1"   
)
nt Draw ( "mc_phi>>hRecoEffPhi(50,-200,200)"  ,
mult,
0"  ,
"E1"   
)
nt Draw ( "mc_p>>hRecoEffCtP(100,0,5)"  ,
mult,
0 &&(stt >0||tpc >0)"  ,
"E1"   
)
nt Draw ( "mc_theta>>hRecoEffCtTheta(50,0,180)"  ,
mult,
0 &&(stt >0||tpc >0)"  ,
"E1"   
)
nt Draw ( "mc_phi>>hRecoEffCtPhi(50,-200,200)"  ,
mult,
0 &&(stt >0||tpc >0)"  ,
"E1"   
)
ntEvt Draw ( "cand>>f1 )
ntEvt Draw ( "eff>>f1 )
ntEvt Draw ( "effct>>f1 )
nt Draw ( "(mc_p-p):mc_p>>hresp_p(100,0,5,100,-0.5,0.5)"  ,
mult,
0 &&(stt >0||tpc >0)"  ,
"colz"   
)
nt Draw ( "(mc_theta-theta):mc_theta>>hrestheta_theta(100,0,100,100,-2,2)"  ,
mult,
0 &&(stt >0||tpc >0)"  ,
"colz"   
)
nt Draw ( "(mc_p-p):mc_theta>>hresp_theta(100,0,100,100,-0.5,0.5)"  ,
mult,
0 &&(stt >0||tpc >0)"  ,
"colz"   
)
nt Draw ( "(mc_theta-theta):mc_p>>hrestheta_p(100,0,5,100,-2,2)"  ,
mult,
0 &&(stt >0||tpc >0)"  ,
"colz"   
)
nt Draw ( "(mc_phi-phi):mc_phi>>hresphi_phi(100,-200,200,100,-2,2)"  ,
mult,
0 &&(stt >0||tpc >0)"  ,
"colz"   
)
nt Draw ( "p>>hist ,
mult,
0 &&(stt >0||tpc >0)"   
)
nt Draw ( "theta>>hist2 ,
mult,
0 &&(stt >0||tpc >0)"   
)
nt Draw ( "(mc_p-p)/mc_p>>hist3 ,
mult,
0 &&(stt >0||tpc >0)"   
)
hist3 Fit ( f10  ,
"R"   
)
hMcSttP GetXaxis ( ) -> SetTitle("Momentum (GeV)")

Definition at line 653 of file anal_hit_digi_cluster_fwendcap.C.

653 {Reconstructed}-#theta_{MC} (#circ)");
hMcAccSttP GetYaxis ( ) -> SetTitle("Acceptance")

Definition at line 131 of file macro/detectors/emc/dedicated/reco_analys.C.

131 {reco} - #theta_{truth}, degree");
int reco()
f1 Scale ( 1/  numero)

Referenced by hypGe_DoubleGeo(), and hypGe_TripleGeo().

leg SetFillColor ( kWhite  )
hMcSttP SetLineColor ( kRed  )
hMcMvdPhi SetLineColor ( )
hMcGemP SetLineColor ( )
gStyle SetOptStat ( kFALSE  )
gStyle SetPalette ( )
hist3 SetStats ( kTRUE  )
hMcSttP SetTitle ( "Geometrical coverage for the momentum"  )
hMcSttTheta SetTitle ( "Geometrical coverage for the theta angle )
hMcSttPhi SetTitle ( "Geometrical coverage for the phi angle )
hMcAccTpcPhi SetTitle ( "Geometrical Acceptance"  )
f1 SetTitle ( "Geometrical Global Acceptance"  )
hRecoEffPhi SetTitle ( "Pattern Recognition Efficiency"  )
hRecoEffCtPhi SetTitle ( "Pattern Recognition Efficiency Only Central Tracker"  )
f1 SetTitle ( "Number of candidates"  )
f1 SetTitle ( "Number of reconstructed primary tracks STT+MVD+GEM )
f1 SetTitle ( "Number of reconstructed primary tracks Only Central Tracker"  )
hresphi_phi SetTitle ( "Resolution"  )
c9 Update ( )

Variable Documentation

TCanvas* c1 =new TCanvas("c1","",800,800)

Definition at line 11 of file plot_acceptance_stt.C.

TCanvas* c11 =new TCanvas("c11","",800,800)
TCanvas* c14 =new TCanvas("c14","",800,800)

Definition at line 135 of file plot_acceptance_stt.C.

Referenced by main().

TCanvas* c20 =new TCanvas("c20","",800,800)
TCanvas* c21 =new TCanvas("c21","",800,800)
TCanvas* c22 =new TCanvas("c22","",800,800)
TCanvas* c4 =new TCanvas("c4","",800,800)

Definition at line 69 of file plot_acceptance_stt.C.

TCanvas* c5 =new TCanvas("c5","",800,800)

Definition at line 94 of file plot_acceptance_stt.C.

TCanvas* c6 =new TCanvas("c6","",800,800)

Definition at line 107 of file plot_acceptance_stt.C.

TCanvas* c7 =new TCanvas("c7","",800,800)

Definition at line 164 of file plot_acceptance_stt.C.

TCanvas* c8 =new TCanvas("c8","",800,800)
TCanvas* c9 =new TCanvas("c9","",800,800)

Definition at line 202 of file plot_acceptance_stt.C.

TFile* f = TFile::Open(inFile)

Definition at line 5 of file plot_acceptance_stt.C.

TH1F* f1 =new TH1F("f1","",100,0,5)

Definition at line 97 of file plot_acceptance_stt.C.

TF1* f10 = new TF1("f10","gaus",-0.09,0.09)

Definition at line 259 of file plot_acceptance_stt.C.

Referenced by error_matrix_fit().

TH1F* hist =new TH1F("hist","Total Momentum distribution",100,0,10)
TH1F* hist2 =new TH1F("hist2","#theta Angle Distribution",100,0,180)

Definition at line 249 of file plot_acceptance_stt.C.

TH1F* hist3 =new TH1F("hist3","(Reconstructed Momentum - MC Momentum )/MC Momentum",100,-1,1)

Definition at line 255 of file plot_acceptance_stt.C.

leg =new TLegend(0.6,0.7,0.89,0.89)
float numero =ntEvt->GetEntries()

Definition at line 95 of file plot_acceptance_stt.C.