Skip to content
Snippets Groups Projects

Acceptance

Merged Acceptance
jihee.kim/reconstruction_benchmarks:acceptance into master
Merged Jihee Kim requested to merge jihee.kim/reconstruction_benchmarks:acceptance into master
Viewing commit a886aad9
Show latest version
1 file
+ 34
− 34
Compare changes
  • Side-by-side
  • Inline
ecal/scripts/rec_emcal_electrons_reader.C
+ 34
− 34
@@ -29,37 +29,37 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c
@@ -29,37 +29,37 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c
t->SetMakeClass(1);
t->SetMakeClass(1);
t->SetBranchStatus("*", 0);
t->SetBranchStatus("*", 0);
Int_t mcparticles2_;
Int_t mcparticles2_;
t->SetBranchStatus("mcparticles2", 1);
//t->SetBranchStatus("mcparticles2", 1);
t->SetBranchAddress("mcparticles2", &mcparticles2_);
//t->SetBranchAddress("mcparticles2", &mcparticles2_);
Int_t CrystalEcalHits2_;
//Int_t CrystalEcalHits2_;
t->SetBranchStatus("CrystalEcalHits2", 1);
//t->SetBranchStatus("CrystalEcalHits2", 1);
t->SetBranchAddress("CrystalEcalHits2", &CrystalEcalHits2_);
//t->SetBranchAddress("CrystalEcalHits2", &CrystalEcalHits2_);
Int_t RecoEcalHits_;
Int_t RecoEcalHits_;
t->SetBranchStatus("RecoEcalHits", 1);
t->SetBranchStatus("RecoEcalHits", 1);
t->SetBranchAddress("RecoEcalHits", &RecoEcalHits_);
t->SetBranchAddress("RecoEcalHits", &RecoEcalHits_);
Int_t EcalClusters_;
Int_t EcalClusters_;
t->SetBranchStatus("EcalClusters", 1);
t->SetBranchStatus("EcalClusters", 1);
t->SetBranchAddress("EcalClusters", &EcalClusters_);
t->SetBranchAddress("EcalClusters", &EcalClusters_);
const Int_t kMaxmcparticles2 = 100000;
//const Int_t kMaxmcparticles2 = 100000;
Double_t px[kMaxmcparticles2];
//Double_t px[kMaxmcparticles2];
Double_t py[kMaxmcparticles2];
//Double_t py[kMaxmcparticles2];
Double_t pz[kMaxmcparticles2];
//Double_t pz[kMaxmcparticles2];
Double_t mass[kMaxmcparticles2];
//Double_t mass[kMaxmcparticles2];
t->SetBranchStatus("mcparticles2.psx",1);
//t->SetBranchStatus("mcparticles2.psx",1);
t->SetBranchStatus("mcparticles2.psy",1);
//t->SetBranchStatus("mcparticles2.psy",1);
t->SetBranchStatus("mcparticles2.psz",1);
//t->SetBranchStatus("mcparticles2.psz",1);
t->SetBranchStatus("mcparticles2.mass",1);
//t->SetBranchStatus("mcparticles2.mass",1);
t->SetBranchAddress("mcparticles2.psx",px);
//t->SetBranchAddress("mcparticles2.psx",px);
t->SetBranchAddress("mcparticles2.psy",py);
//t->SetBranchAddress("mcparticles2.psy",py);
t->SetBranchAddress("mcparticles2.psz",pz);
//t->SetBranchAddress("mcparticles2.psz",pz);
t->SetBranchAddress("mcparticles2.mass",mass);
//t->SetBranchAddress("mcparticles2.mass",mass);
const Int_t kMaxCrystalEcalHits2 = 100000;
//const Int_t kMaxCrystalEcalHits2 = 100000;
Double_t truth_deposit[kMaxCrystalEcalHits2];
//Double_t truth_deposit[kMaxCrystalEcalHits2];
Double_t energyDeposit[kMaxCrystalEcalHits2];
//Double_t energyDeposit[kMaxCrystalEcalHits2];
t->SetBranchStatus("CrystalEcalHits2.truth.deposit",1);
//t->SetBranchStatus("CrystalEcalHits2.truth.deposit",1);
t->SetBranchStatus("CrystalEcalHits2.energyDeposit",1);
//t->SetBranchStatus("CrystalEcalHits2.energyDeposit",1);
t->SetBranchAddress("CrystalEcalHits2.truth.deposit",truth_deposit);
//t->SetBranchAddress("CrystalEcalHits2.truth.deposit",truth_deposit);
t->SetBranchAddress("CrystalEcalHits2.energyDeposit",energyDeposit);
//t->SetBranchAddress("CrystalEcalHits2.energyDeposit",energyDeposit);
const Int_t kMaxRecoEcalHits = 100000;
const Int_t kMaxRecoEcalHits = 100000;
Double_t rec_x_pos[kMaxRecoEcalHits];
Double_t rec_x_pos[kMaxRecoEcalHits];
Double_t rec_y_pos[kMaxRecoEcalHits];
Double_t rec_y_pos[kMaxRecoEcalHits];
@@ -147,8 +147,8 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c
@@ -147,8 +147,8 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c
// Read event by event
// Read event by event
t->GetEntry(ievent);
t->GetEntry(ievent);
// Read number of hits/clusters
// Read number of hits/clusters
Int_t nmcparticle = 2;
//Int_t nmcparticle = 2;
Int_t nCrystalEcalHits = CrystalEcalHits2_;
//Int_t nCrystalEcalHits = CrystalEcalHits2_;
Int_t nreconhits = RecoEcalHits_;
Int_t nreconhits = RecoEcalHits_;
Int_t ncluster = EcalClusters_;
Int_t ncluster = EcalClusters_;
// Initialize total energy variables
// Initialize total energy variables
@@ -157,11 +157,11 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c
@@ -157,11 +157,11 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c
total_sim_e = 0.0;
total_sim_e = 0.0;
total_cluster_e = 0.0;
total_cluster_e = 0.0;
// Thrown energy, momentum, and mass
// Thrown energy, momentum, and mass
momentum2 = px[nmcparticle]*px[nmcparticle]+py[nmcparticle]*py[nmcparticle]+pz[nmcparticle]*pz[nmcparticle];
//momentum2 = px[nmcparticle]*px[nmcparticle]+py[nmcparticle]*py[nmcparticle]+pz[nmcparticle]*pz[nmcparticle];
momentum = TMath::Sqrt(momentum2);
//momentum = TMath::Sqrt(momentum2);
mass2 = mass[nmcparticle]*mass[nmcparticle];
//mass2 = mass[nmcparticle]*mass[nmcparticle];
total_thr_e = sqrt(momentum2 + mass2)/1.e+3;
//total_thr_e = sqrt(momentum2 + mass2)/1.e+3;
h12->Fill(momentum, 1.0);
//h12->Fill(momentum, 1.0);
// Loop over simulated hit by simulated hit
// Loop over simulated hit by simulated hit
for (int isimhit=0; isimhit < nCrystalEcalHits; isimhit++)
for (int isimhit=0; isimhit < nCrystalEcalHits; isimhit++)
@@ -212,8 +212,8 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c
@@ -212,8 +212,8 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c
h11->Fill(eng_res, 1.0);
h11->Fill(eng_res, 1.0);
}
}
if(total_cluster_e > 0.9*total_thr_e)
//if(total_cluster_e > 0.9*total_thr_e)
h13->Fill(momentum, 1.0);
// h13->Fill(momentum, 1.0);
}
}
}
}
// Drawing and Saving figures
// Drawing and Saving figures