54 |
|
|
55 |
|
pAngle::pAngle(SimInfo* info, const std::string& filename, |
56 |
|
const std::string& sele1, int nthetabins) |
57 |
< |
: StaticAnalyser(info, filename), selectionScript1_(sele1), |
58 |
< |
evaluator1_(info), evaluator2_(info), seleMan1_(info), seleMan2_(info), |
59 |
< |
nThetaBins_(nthetabins), |
60 |
< |
doVect_(true), doOffset_(false) { |
57 |
> |
: StaticAnalyser(info, filename), doVect_(true), doOffset_(false), |
58 |
> |
selectionScript1_(sele1), seleMan1_(info), seleMan2_(info), |
59 |
> |
evaluator1_(info), evaluator2_(info), |
60 |
> |
nThetaBins_(nthetabins) { |
61 |
|
|
62 |
|
setOutputName(getPrefix(filename) + ".pAngle"); |
63 |
|
|
73 |
|
pAngle::pAngle(SimInfo* info, const std::string& filename, |
74 |
|
const std::string& sele1, const std::string& sele2, |
75 |
|
int nthetabins) |
76 |
< |
: StaticAnalyser(info, filename), selectionScript1_(sele1), |
77 |
< |
selectionScript2_(sele2), evaluator1_(info), evaluator2_(info), |
78 |
< |
seleMan1_(info), seleMan2_(info), nThetaBins_(nthetabins), |
79 |
< |
doVect_(false), doOffset_(false) { |
76 |
> |
: StaticAnalyser(info, filename), doVect_(false), doOffset_(false), |
77 |
> |
selectionScript1_(sele1), selectionScript2_(sele2), |
78 |
> |
seleMan1_(info), seleMan2_(info), evaluator1_(info), evaluator2_(info), |
79 |
> |
nThetaBins_(nthetabins) { |
80 |
|
|
81 |
|
setOutputName(getPrefix(filename) + ".pAngle"); |
82 |
|
|
96 |
|
|
97 |
|
pAngle::pAngle(SimInfo* info, const std::string& filename, |
98 |
|
const std::string& sele1, int seleOffset, int nthetabins) |
99 |
< |
: StaticAnalyser(info, filename), selectionScript1_(sele1), |
100 |
< |
evaluator1_(info), evaluator2_(info), seleMan1_(info), seleMan2_(info), |
101 |
< |
nThetaBins_(nthetabins), seleOffset_(seleOffset), |
102 |
< |
doVect_(false), doOffset_(true), doOffset2_(false) { |
99 |
> |
: StaticAnalyser(info, filename), doVect_(false), doOffset_(true), |
100 |
> |
doOffset2_(false), selectionScript1_(sele1), |
101 |
> |
seleMan1_(info), seleMan2_(info), evaluator1_(info), evaluator2_(info), |
102 |
> |
seleOffset_(seleOffset), nThetaBins_(nthetabins) { |
103 |
|
|
104 |
|
setOutputName(getPrefix(filename) + ".pAngle"); |
105 |
|
|
115 |
|
pAngle::pAngle(SimInfo* info, const std::string& filename, |
116 |
|
const std::string& sele1, int seleOffset, int seleOffset2, |
117 |
|
int nthetabins) |
118 |
< |
: StaticAnalyser(info, filename), selectionScript1_(sele1), |
119 |
< |
evaluator1_(info), evaluator2_(info), seleMan1_(info), seleMan2_(info), |
120 |
< |
nThetaBins_(nthetabins), seleOffset_(seleOffset), |
121 |
< |
seleOffset2_(seleOffset2), |
122 |
< |
doVect_(false), doOffset_(true), doOffset2_(true) { |
118 |
> |
: StaticAnalyser(info, filename), doVect_(false), doOffset_(true), |
119 |
> |
doOffset2_(true), selectionScript1_(sele1), |
120 |
> |
seleMan1_(info), seleMan2_(info), evaluator1_(info), evaluator2_(info), |
121 |
> |
seleOffset_(seleOffset), seleOffset2_(seleOffset2), |
122 |
> |
nThetaBins_(nthetabins) { |
123 |
|
|
124 |
|
setOutputName(getPrefix(filename) + ".pAngle"); |
125 |
|
|
182 |
|
// only do this if the stunt double actually has a vector associated |
183 |
|
// with it |
184 |
|
if (sd1->isDirectional()) { |
185 |
< |
Vector3d vec = sd1->getA().getColumn(2); |
186 |
< |
RealType distance = r1.length(); |
187 |
< |
|
185 |
> |
Vector3d vec = sd1->getA().getColumn(2); |
186 |
|
vec.normalize(); |
187 |
|
r1.normalize(); |
188 |
|
RealType cosangle = dot(r1, vec); |