diff --git a/test/full/fitcentroid/cylinder/cylinder.geo b/test/full/fitcentroid/cylinder/cylinder.geo index 4cb29979..24e907ed 100644 --- a/test/full/fitcentroid/cylinder/cylinder.geo +++ b/test/full/fitcentroid/cylinder/cylinder.geo @@ -12,7 +12,7 @@ invisible: 1, { name: "GEO", -index: "water", +index: "detector", valid_begin: [0, 0], valid_end: [0, 0], mother: "world", @@ -27,10 +27,10 @@ color: [0.4, 0.4, 0.6, 0.1], // use an inner volume to speed up simulation. { name: "GEO", -index: "water_inner", +index: "detector_inner", valid_begin: [0, 0], valid_end: [0, 0], -mother: "water", +mother: "detector", type: "tube", r_max: 2500.0, size_z: 2500.0, @@ -45,7 +45,7 @@ name: "GEO", index: "pmts", valid_begin: [0, 0], valid_end: [0, 0], -mother: "water", +mother: "detector", type: "pmtarray", pmt_model: "r1408", pmt_detector_type: "idpmt", diff --git a/test/full/fitcentroid/fitcentroid.C b/test/full/fitcentroid/fitcentroid.C index 7521c445..a34d0a73 100644 --- a/test/full/fitcentroid/fitcentroid.C +++ b/test/full/fitcentroid/fitcentroid.C @@ -1,29 +1,3 @@ -// void make_plots(TFile *event_file, TTree *T, TFile *out_file) -// { -// TH1F *hCenterRes = new TH1F("hCenterRes", "20 keVee electrons at the center;Reconstructed X (mm);Events per bin", -// 40, -800, 800); TProfile *hReconVsRadius = new TProfile("hReconVsRadius", "20 keVee electrons (error bars are -// RMS);(Radius/437 mm)**3;Reconstructed X - True X (mm)", 5, 0, 1, "s"); TProfile *hReconVsEnergy = new -// TProfile("hReconVsEnergy", "Electrons at center (error bars are RMS);Kinetic energy (keVee);Reconstructed X (mm)", -// 5, 20, 100,"s"); -// -// T->Draw("ev.centroid.pos.fX>>hCenterRes","mc.particle.ke < 0.0205 && mc.particle.pos.fX < 0.05","goff"); -// T->Draw("ev.centroid.pos.fX-mc.particle.pos.fX:(mc.particle.pos.Mag()/437)>>hReconVsRadius","mc.particle.ke < -// 0.0205 && mc.particle.pos.fX > 0.05","goff prof"); -// T->Draw("ev.centroid.pos.fX:mc.particle.ke*1000>>hReconVsEnergy","mc.particle.ke > 0.0205 && mc.particle.pos.fX < -// 0.05","goff prof"); -// -// hCenterRes->Fit("gaus"); -// hReconVsRadius->SetMinimum(-200); -// hReconVsRadius->SetMaximum(200); -// hReconVsEnergy->SetMinimum(-200); -// hReconVsEnergy->SetMaximum(200); -// -// out_file->cd(); -// hCenterRes->Write(); -// hReconVsRadius->Write(); -// hReconVsEnergy->Write(); -// } - void fitcentroid(std::string event_filename, std::string out_filename) { TFile *event_file = new TFile(event_filename.c_str(), "READ"); TTree *T = (TTree *)event_file->Get("T");