From 22bb217700c336821c0acb2617ef513ff07f1feb Mon Sep 17 00:00:00 2001 From: "jihee.kim" <jihee.kim@anl.gov> Date: Wed, 21 Oct 2020 16:08:30 -0500 Subject: [PATCH] Dropped mcparticles2 and CrystalEcalHits2 branches to read --- ecal/scripts/rec_emcal_electrons_reader.C | 11 ++++++++++- 1 file changed, 10 insertions(+), 1 deletion(-) diff --git a/ecal/scripts/rec_emcal_electrons_reader.C b/ecal/scripts/rec_emcal_electrons_reader.C index 63490be4..38d59ba3 100644 --- a/ecal/scripts/rec_emcal_electrons_reader.C +++ b/ecal/scripts/rec_emcal_electrons_reader.C @@ -22,18 +22,21 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c // Set Branch status and addressed t->SetMakeClass(1); t->SetBranchStatus("*", 0); +/* Int_t mcparticles2_; t->SetBranchStatus("mcparticles2", 1); t->SetBranchAddress("mcparticles2", &mcparticles2_); Int_t CrystalEcalHits2_; t->SetBranchStatus("CrystalEcalHits2", 1); t->SetBranchAddress("CrystalEcalHits2", &CrystalEcalHits2_); +*/ Int_t RecoEcalHits_; t->SetBranchStatus("RecoEcalHits", 1); t->SetBranchAddress("RecoEcalHits", &RecoEcalHits_); Int_t EcalClusters_; t->SetBranchStatus("EcalClusters", 1); t->SetBranchAddress("EcalClusters", &EcalClusters_); +/* const Int_t kMaxmcparticles2 = 100000; Double_t px[kMaxmcparticles2]; Double_t py[kMaxmcparticles2]; @@ -54,6 +57,7 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c t->SetBranchStatus("CrystalEcalHits2.energyDeposit",1); t->SetBranchAddress("CrystalEcalHits2.truth.deposit",truth_deposit); t->SetBranchAddress("CrystalEcalHits2.energyDeposit",energyDeposit); +*/ const Int_t kMaxRecoEcalHits = 100000; Double_t rec_x_pos[kMaxRecoEcalHits]; Double_t rec_y_pos[kMaxRecoEcalHits]; @@ -141,8 +145,10 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c // Read event by event t->GetEntry(ievent); // Read number of hits/clusters +/* Int_t nmcparticle = 2; Int_t nCrystalEcalHits = CrystalEcalHits2_; +*/ Int_t nreconhits = RecoEcalHits_; Int_t ncluster = EcalClusters_; // Initialize total energy variables @@ -151,6 +157,7 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c total_sim_e = 0.0; total_cluster_e = 0.0; // Thrown energy, momentum, and mass +/* momentum2 = px[nmcparticle]*px[nmcparticle]+py[nmcparticle]*py[nmcparticle]+pz[nmcparticle]*pz[nmcparticle]; momentum = TMath::Sqrt(momentum2); mass2 = mass[nmcparticle]*mass[nmcparticle]; @@ -163,7 +170,7 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c total_truth_sim_e += truth_deposit[isimhit]; total_sim_e += energyDeposit[isimhit]/1.e+3; } - +*/ // Loop over reconstructed hit by reconstructed hit for(int ireconhit=0; ireconhit < nreconhits; ireconhit++) h9->Fill(rec_x_pos[ireconhit],rec_y_pos[ireconhit], 1.0); @@ -203,8 +210,10 @@ int rec_emcal_electrons_reader(double e_start = 1.0, double e_end = 1.0, const c h7->Fill(eta, 1.0); h11->Fill(eng_res, 1.0); } +/* if(total_cluster_e > 0.9*total_thr_e) h13->Fill(momentum, 1.0); +*/ } } // Drawing and Saving figures -- GitLab