simple analysis - adding cluster energy calibration

benchmarks/diffractive_vm/analysis/diffractive_vm_simple_analysis.cxx

@@ -103,6 +103,7 @@ int diffractive_vm_simple_analysis(const std::string& config_name)
    	//energy clus
     TH2D* h_emClus_position_REC = new TH2D("h_emClus_position_REC",";x (cm);y (cm)",400,-800,800,400,-800,800);
 	TH2D* h_emHits_position_REC = new TH2D("h_emHits_position_REC",";x (cm);y (cm)",400,-800,800,400,-800,800);
+	TH2D* h_emHits_position_2_REC = new TH2D("h_emHits_position_2_REC",";x (cm);y (cm)",400,-800,800,400,-800,800);
     TH2D* h_energy_res = new TH2D("h_energy_res",";E_{MC} (GeV); E_{MC}-E_{REC}/E_{MC} emcal",100,0,20,1000,-1,1);
     TH1D* h_energy_calibration_REC = new TH1D("h_energy_calibration_REC",";E (GeV)",200,0,2);
 
@@ -193,6 +194,9 @@ int diffractive_vm_simple_analysis(const std::string& config_name)
 		h_energy_REC->Fill(maxEnergy);
 		h_emClus_position_REC->Fill(xpos,ypos);
 		h_emHits_position_REC->Fill(xhitpos,yhitpos);
+		if(fabs(xpos)<100 && fabs(ypos)<100){
+	        h_emHits_position_2_REC->Fill(xhitpos,yhitpos,maxEnergy);
+		}
 
 		//association of rec level scat' e
 		int rec_elect_index=-1;