Skip to content
Snippets Groups Projects
Commit f7e1144e authored by Zhenyu Ye's avatar Zhenyu Ye
Browse files

add TOF Material Plots

parent cdf2252d
Branches
No related tags found
1 merge request!100Resolve "Add BTOF benchmarks"
This commit is part of merge request !100. Comments created here will be created in the context of that merge request.
...@@ -57,11 +57,12 @@ void materialScanPhi( ...@@ -57,11 +57,12 @@ void materialScanPhi(
} }
// get material scans // get material scans
/*
size_t total{0}; size_t total{0};
std::vector<dd4hep::rec::MaterialVec> scan; std::vector<dd4hep::rec::MaterialVec> scan;
double x0{0}, y0{0}, z0{0}; double x0{0}, y0{0}, z0{0};
for (double phi = phimin; phi <= phimax + 0.5*phistep; phi += phistep) { for (double phi = phimin; phi <= phimax + 0.5*phistep; phi += phistep) {
std::cout << phi << std::endl;
double theta = 2.0 * (atan(1) - atan(exp(-eta))); double theta = 2.0 * (atan(1) - atan(exp(-eta)));
double r = min((theta > 0? zpmax: -znmax) / sin(theta), min(rmax, rhomax / cos(theta))); double r = min((theta > 0? zpmax: -znmax) / sin(theta), min(rmax, rhomax / cos(theta)));
double x = r * cos(theta) * cos(phi); double x = r * cos(theta) * cos(phi);
...@@ -69,7 +70,7 @@ void materialScanPhi( ...@@ -69,7 +70,7 @@ void materialScanPhi(
double z = r * sin(theta); double z = r * sin(theta);
scan.emplace_back(gMaterialScan->scan(x0,y0,z0,x,y,z)); scan.emplace_back(gMaterialScan->scan(x0,y0,z0,x,y,z));
total += scan.back().size(); total += scan.back().size();
}*/ }
// start creating histograms for stacking: // start creating histograms for stacking:
// - start with first material layer at central eta bin // - start with first material layer at central eta bin
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment