Go to the source code of this file.
|
TFile * | f = new TFile("full_emc.root") |
|
TTree * | t =(TTree *) f->Get("pndsim") |
|
TClonesArray * | bump_array =new TClonesArray("PndEmcBump") |
|
TFile * | fsim = new TFile("sim_emc.root") |
|
TTree * | tsim =(TTree *) fsim->Get("pndsim") |
|
TClonesArray * | mctrack_array =new TClonesArray("PndMCTrack") |
|
double | bump_energy =bump->energy() |
|
double | bump_energy_corrected =bump->GetEnergyCorrected() |
|
TH1F * | h1 = new TH1F("h1","Cluster energy",100,0.0,1.10) |
|
TH1F * | h2 = new TH1F("h2","Cluster energy corrected",100,0.0,1.10) |
|
TCanvas * | c1 = new TCanvas("c1", "Cluster Energy", 100, 100, 800, 800) |
|
Definition at line 64 of file NHitsPerEvent.C.
73 std::cout <<
"NHits: " <<
gHitPoints->GetN() << std::endl;
void ExtractData(TClonesArray *array, TGraph *g)
TClonesArray * mvdHitArray
TClonesArray * mvdProjectedArray
Definition at line 1 of file outdated/mpiTools/macros/emc/rootlogon.C.
3 gSystem->AddIncludePath(
"-I$VMCWORKDIR/emc");
4 gSystem->AddIncludePath(
"-I$VMCWORKDIR/emc/EmcDigi");
5 gSystem->AddIncludePath(
"-I$VMCWORKDIR/emc/EmcMC");
6 gSystem->AddIncludePath(
"-I$VMCWORKDIR/emc/EmcTools");
7 gSystem->AddIncludePath(
"-I$VMCWORKDIR/emc/EmcReco");
8 gSystem->AddIncludePath(
"-I$VMCWORKDIR/emc/EmcData");
9 gSystem->AddIncludePath(
"-I$VMCWORKDIR/pnddata");
10 gSystem->AddIncludePath(
"-I$VMCWORKDIR/base");
11 gROOT->LoadMacro(
"$VMCWORKDIR/gconfig/basiclibs.C");
14 gSystem->Load(
"libDpmEvtGen.so");
15 gSystem->Load(
"libGeoBase");
16 gSystem->Load(
"libParBase");
17 gSystem->Load(
"libBase");
18 gSystem->Load(
"libTrkBase");
19 gSystem->Load(
"libPndData");
20 gSystem->Load(
"libField");
21 gSystem->Load(
"libPassive");
22 gSystem->Load(
"libGen");
23 gSystem->Load(
"libPGen");
24 gSystem->Load(
"libEmc");
27 gRandom=
new TRandom3(0);
t SetBranchAddress |
( |
"EmcBump" |
, |
|
|
& |
bump_array |
|
) |
| |
tsim SetBranchAddress |
( |
"MCTrack" |
, |
|
|
& |
mctrack_array |
|
) |
| |
h1 SetTitle |
( |
"Cluster energy of 1 GeV photon" |
| ) |
|
h2 SetTitle |
( |
"Cluster energy of 1 GeV photon (corrected)" |
| ) |
|
TClonesArray* bump_array =new TClonesArray("PndEmcBump") |
bump_energy_corrected =bump->GetEnergyCorrected() |
TCanvas* c1 = new TCanvas("c1", "Cluster Energy", 100, 100, 800, 800) |
TFile* f = new TFile("full_emc.root") |
TFile* fsim = new TFile("sim_emc.root") |
TH1F* h1 = new TH1F("h1","Cluster energy",100,0.0,1.10) |
TH1F* h2 = new TH1F("h2","Cluster energy corrected",100,0.0,1.10) |
TClonesArray* mctrack_array =new TClonesArray("PndMCTrack") |
TTree* t =(TTree *) f->Get("pndsim") |
TTree* tsim =(TTree *) fsim->Get("pndsim") |