diff --git a/Armory/anasenMS.cpp b/Armory/anasenMS.cpp index 3d0add5..7d225c1 100644 --- a/Armory/anasenMS.cpp +++ b/Armory/anasenMS.cpp @@ -203,8 +203,8 @@ int main(int argc, char **argv){ //TGraph* elossLight = LoadELoss("../ELoss/HeLoss/E_vs_x_" + b + ".dat"); // define vertex position uniform distribution ranges (mm) - double vertexXRange[2] = { -3.5 / 2, 3.5 / 2}; // mm - 5, 5// 2.5 - double vertexYRange[2] = { -3.5 / 2, 3.5 / 2}; // -5, 5 + double vertexXRange[2] = { -1.75, 1.75}; // mm - 5, 5// 2.5 + double vertexYRange[2] = { -1.75, 1.75}; // -5, 5 double vertexZRange[2] = { -174.3, 174.3}; // -174.3, 174.3 (full length of gas volume, centered at 0) const double beamEntranceZ = -280 - 174.3; //vertexZRange[0]; // mm, assumed beam entrance into the gas diff --git a/Armory/testMacro.C b/Armory/testMacro.C index 3aca5e0..5472b1e 100644 --- a/Armory/testMacro.C +++ b/Armory/testMacro.C @@ -36,10 +36,10 @@ new TCanvas("c1", "Canvas", 900, 600); tree1->Draw("Ex_kin:vZ_track>>h2(800,-500,300,110,-1,10)", "sx3ID >=0 && EPC > .02", "colz"); - //tree1->Draw("Ex:vZ", "sx3ID >=0 && EPC > .02", "same"); + tree1->Draw("Ex:vZ", "sx3ID >=0 && EPC > .02", "same"); c1->SaveAs("Ex_kin_vz.png"); //{0, 2.235, 3.498, 3.774, 4.831, 4.809, 5.279, 5.614, 5.951, 6.550}; -/* + new TCanvas("c2", "Canvas", 900, 600); tree1->Draw("theta_recon:thetab", "sx3ID >=0 && Esx3 > 0", ""); c2->SaveAs("theta_recon_vs_thetab.png"); @@ -66,7 +66,7 @@ tree1->Draw("EBeam_kin:beamEnergy", "sx3ID >=0 && !TMath::IsNaN(EBeam_kin)", ""); TF1 *line_kin_kin = new TF1("line_beam_kin", "x", gPad->GetUxmin(), gPad->GetUxmax()); line_kin_kin->SetLineColor(kRed); line_kin_kin->Draw("SAME"); - c6->SaveAs("Ebeam_kin_vs_beamEnergy.png");*/ + c6->SaveAs("Ebeam_kin_vs_beamEnergy.png"); new TCanvas("c7", "Canvas", 900, 600); //c7->SetLogy();