FairRoot/PandaRoot
Functions | Variables
plot_acceptance_tpc.C File Reference

Go to the source code of this file.

Functions

gStyle SetOptStat (kFALSE)
 
gStyle SetPalette (1)
 
c4 Divide (1, 3)
 
c4 cd (1)
 
nt Draw ("mc_p>>hMcP(100,0,5)")
 
nt Draw ("mc_p>>hMcAccTpcP(100,0,5)","mc_tpc>0","E1")
 
hMcAccTpcP Divide (hMcP)
 
hMcAccTpcP SetTitle ("Geometrical Acceptance")
 
hMcAccTpcP GetXaxis () -> SetTitle("Momentum (GeV)")
 
hMcAccTpcP GetYaxis () -> SetTitle("Acceptance")
 
c4 cd (2)
 
nt Draw ("mc_theta>>hMcTheta(50,0,180)")
 
nt Draw ("mc_theta>>hMcAccTpcTheta(50,0,180)","mc_tpc>0","E1")
 
hMcAccTpcTheta Divide (hMcTheta)
 
c4 cd (3)
 
nt Draw ("mc_phi>>hMcPhi(50,-200,200)")
 
nt Draw ("mc_phi>>hMcAccTpcPhi(50,-200,200)","mc_tpc>0","E1")
 
hMcAccTpcPhi Divide (hMcPhi)
 
c5 cd ()
 
ntEvt Draw ("acc_tpc>>f1")
 
f1 SetTitle ("Geometrical Global Acceptance")
 
f1 Draw ()
 
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")
 
ntEvt Draw ("cand>>f1")
 
f1 SetTitle ("Number of candidates")
 
ntEvt Draw ("eff>>f1")
 
f1 SetTitle ("Number of reconstructed primary tracks")
 
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 * 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 * c7 =new TCanvas("c7","",800,800)
 
TCanvas * c8 =new TCanvas("c8","",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

c4 cd ( 1  )
c4 cd ( 2  )
c4 cd ( 3  )
c5 cd ( )
c9 cd ( 4  )
c9 cd ( 5  )
c4 Divide ( 1  ,
3   
)
hMcAccTpcP Divide ( hMcP  )
hMcAccTpcTheta Divide ( hMcTheta  )
hMcAccTpcPhi Divide ( hMcPhi  )
c9 Divide ( 3  ,
2   
)
nt Draw ( "mc_p>>hMcP(100,0,5)"  )
nt Draw ( "mc_p>>hMcAccTpcP(100,0,5)"  ,
"  mc_tpc,
0"  ,
"E1"   
)
nt Draw ( "mc_theta>>hMcTheta(50,0,180)"  )
nt Draw ( "mc_theta>>hMcAccTpcTheta(50,0,180)"  ,
"  mc_tpc,
0"  ,
"E1"   
)
nt Draw ( "mc_phi>>hMcPhi(50,-200,200)"  )
nt Draw ( "mc_phi>>hMcAccTpcPhi(50,-200,200)"  ,
"  mc_tpc,
0"  ,
"E1"   
)
ntEvt Draw ( "acc_tpc>>f1"  )
f1 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();
69  ExtractData(mvdHitArray, gHitPoints);
70  ExtractData(mvdProjectedArray, gProjections);
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 }
ExtractData
void ExtractData(TClonesArray *array, TGraph *g)
Definition: NHitsPerEvent.C:78
gHitPoints
TGraph * gHitPoints
Definition: NHitsPerEvent.C:15
mvdHitArray
TClonesArray * mvdHitArray
Definition: NHitsPerEvent.C:12
mvdProjectedArray
TClonesArray * mvdProjectedArray
Definition: NHitsPerEvent.C:13
gProjections
TGraph * gProjections
Definition: NHitsPerEvent.C:16
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"   
)
ntEvt Draw ( "cand>>f1"  )
ntEvt Draw ( "eff>>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"   
)
hMcAccTpcP GetXaxis ( ) -> SetTitle("Momentum (GeV)")

Definition at line 653 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by Pic_FEP_DEGASComparison(), Pic_FEP_Eff_lEE_XiManyTargets(), Pic_OmegaQuadrupoleComparison(), PndFtsHoughSpace::setParametersForHsOption(), PndFtsHoughSpace::WriteHistoOfAllPaths(), and PndFtsHoughSpace::WriteHistoOfAllPathsForEachMcTruthTrack().

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

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

131 {reco} - #theta_{truth}, degree");
degree
degree
Definition: macro/detectors/emc/dedicated/reco_analys.C:130
reco
int reco()
Definition: detectors/disc/reco.C:1
f1 Scale ( 1/  numero)
gStyle SetOptStat ( kFALSE  )
gStyle SetPalette ( 1  )
hist3 SetStats ( kTRUE  )
hMcAccTpcP SetTitle ( "Geometrical Acceptance"  )
f1 SetTitle ( "Geometrical Global Acceptance"  )
hRecoEffP SetTitle ( "Pattern Recognition Efficiency"  )
f1 SetTitle ( "Number of candidates"  )
f1 SetTitle ( "Number of reconstructed primary tracks"  )
hresp_p SetTitle ( "Resolution"  )
c9 Update ( )

Variable Documentation

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

Definition at line 139 of file plot_acceptance_tpc.C.

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

Definition at line 145 of file plot_acceptance_tpc.C.

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

Definition at line 151 of file plot_acceptance_tpc.C.

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

Definition at line 11 of file plot_acceptance_tpc.C.

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

Definition at line 36 of file plot_acceptance_tpc.C.

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

Definition at line 49 of file plot_acceptance_tpc.C.

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

Definition at line 76 of file plot_acceptance_tpc.C.

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

Definition at line 88 of file plot_acceptance_tpc.C.

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

Definition at line 100 of file plot_acceptance_tpc.C.

TFile* f = TFile::Open(inFile)

Definition at line 5 of file plot_acceptance_tpc.C.

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

Definition at line 39 of file plot_acceptance_tpc.C.

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

Definition at line 157 of file plot_acceptance_tpc.C.

TH1F* hist =new TH1F("hist","Total Momentum distribution",100,0,10)

Definition at line 141 of file plot_acceptance_tpc.C.

TH1F* hist2 =new TH1F("hist2","#theta Angle Distribution",100,0,180)

Definition at line 147 of file plot_acceptance_tpc.C.

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

Definition at line 153 of file plot_acceptance_tpc.C.

float numero =ntEvt->GetEntries()

Definition at line 37 of file plot_acceptance_tpc.C.

Generated on Wed Apr 3 2019 10:03:14 for FairRoot/PandaRoot by   doxygen 1.8.5