Acceptance

• Created unified script to read reconstructed and thrown information
• Added acceptance calcuation

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,7 +157,8 @@ 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];
+/*
+	momentum2 = px[nmcparticle]*px[nmcparticle]+py[nmcparticle]*py[nmcparticle]+pz[nmcparticle]*pz[nmcparticle];
         momentum  = TMath::Sqrt(momentum2);
 	mass2     = mass[nmcparticle]*mass[nmcparticle];
         total_thr_e = sqrt(momentum2 + mass2)/1.e+3;
@@ -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,9 +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