| 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); |