drawRecoP.C
| 1 |
{
|
|---|---|
| 2 |
TFile* f = new TFile("sim_digi.root"); |
| 3 |
TTree* t = (TTree*)f->Get("er");
|
| 4 |
t->AddFriend("er","reco.root"); |
| 5 |
|
| 6 |
TCanvas *c1 = new TCanvas("c1", "Telescope Coordinate reconstruction", 1400, 1000); |
| 7 |
c1->Divide(1,3); |
| 8 |
|
| 9 |
c1->cd(1);
|
| 10 |
gPad->SetLogy(); |
| 11 |
t->Draw("ERQTelescopeParticle_T1_DoubleSi_SD2_XY_0_1000020030.fLVTelescope.P()-(ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPx**2+ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPy**2+ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPz**2)**0.5");
|
| 12 |
|
| 13 |
c1->cd(2);
|
| 14 |
gPad->SetLogy(); |
| 15 |
//t->Draw("ERQTelescopeParticle_T1_DoubleSi_SD2_XY_0_1000020030.fLVTelescope.Px()-ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPx","ERQTelescopeSiDigi_T1_DoubleSi_SD2_XY_0_Y@.GetEntries() == 1 ");
|
| 16 |
t->Draw("ERQTelescopeParticle_T1_DoubleSi_SD2_XY_0_1000020030.fLVTelescope.P()-(ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPx**2+ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPy**2+ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPz**2)**0.5","ERQTelescopeSiDigi_T1_DoubleSi_SD2_XY_0_Y@.GetEntries() == 1 && ERQTelescopeSiDigi_T1_SingleSi_1_X_0@.GetEntries() == 1 && ERQTelescopeSiDigi_T1_SingleSi_1_X_0.fEdep > 0.001 && ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPDG==1000020030"); |
| 17 |
|
| 18 |
c1->cd(3);
|
| 19 |
gPad->SetLogy(); |
| 20 |
//t->Draw("ERQTelescopeParticle_T1_DoubleSi_SD2_XY_0_1000020030.fLVTelescope.Px()-ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPx","ERQTelescopeSiDigi_T1_DoubleSi_SD2_XY_0_Y@.GetEntries() == 1 ");
|
| 21 |
t->Draw("ERQTelescopeParticle_T1_DoubleSi_SD2_XY_0_1000020030.fLVTelescope.P()-(ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPx**2+ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPy**2+ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPz**2)**0.5","ERQTelescopeSiDigi_T1_DoubleSi_SD2_XY_0_Y@.GetEntries() == 1 && ERQTelescopeSiDigi_T1_SingleSi_1_X_0@.GetEntries() == 1 && ERQTelescopeSiDigi_T1_SingleSi_1_X_0.fEdep > 0.001 && Sum$(ERQTelescopeCsIPoint_T1_CsI_1_0.fPx_out)==0 && ERQTelescopeSiPoint_T1_DoubleSi_SD2_XY_0_Y.fPDG==1000020030"); |
| 22 |
|
| 23 |
} |