rtang/ANASEN_analysis
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

modified: .gitignore

Browse Source 
modified:   Analyzer.C
This commit is contained in:
vs19g 2025-01-17 11:09:46 -05:00
parent f5a4b750fc
commit 20bf766dc5
2 changed files with 1 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
.gitignore vendored
View File

@ -10,3 +10,4 @@ AnasenMS
data/ data/
data_proton/ data_proton/
root_data/

4
Analyzer.C
View File

@ -302,18 +302,14 @@ Bool_t Analyzer::Process(Long64_t entry)
pwinstance.ConstructGeo(); pwinstance.ConstructGeo();
Coord Crossover[24][24]; Coord Crossover[24][24];
std::vector<std::tuple<int, TVector3, TVector3>> AnGeo;
std::vector<std::tuple<int, TVector3, TVector3>> CatGeo;
TVector3 a, c, diff; TVector3 a, c, diff;
double a2, ac, c2, adiff, cdiff, denom, alpha, beta; double a2, ac, c2, adiff, cdiff, denom, alpha, beta;
int index = 0; int index = 0;
for (int i = 0; i < pwinstance.An.size(); i++) for (int i = 0; i < pwinstance.An.size(); i++)
{ {
AnGeo.push_back(std::make_tuple(i, pwinstance.An[i].first, pwinstance.An[i].second));
a = pwinstance.An[i].first - pwinstance.An[i].second; a = pwinstance.An[i].first - pwinstance.An[i].second;
for (int j = 0; j < pwinstance.Ca.size(); j++) for (int j = 0; j < pwinstance.Ca.size(); j++)
{ {
CatGeo.push_back(std::make_tuple(j, pwinstance.Ca[j].first, pwinstance.Ca[j].second));
c = pwinstance.Ca[j].first - pwinstance.Ca[j].second; c = pwinstance.Ca[j].first - pwinstance.Ca[j].second;
diff = pwinstance.An[i].first - pwinstance.Ca[j].first; diff = pwinstance.An[i].first - pwinstance.Ca[j].first;
a2 = a.Dot(a); a2 = a.Dot(a);