Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
//************************************************************************//
// //
// Copyright 2013 Bertram Kopf (bertram@ep1.rub.de) //
// Julian Pychy (julian@ep1.rub.de) //
// - Ruhr-Universität Bochum //
// //
// This file is part of Pawian. //
// //
// Pawian is free software: you can redistribute it and/or modify //
// it under the terms of the GNU General Public License as published by //
// the Free Software Foundation, either version 3 of the License, or //
// (at your option) any later version. //
// //
// Pawian is distributed in the hope that it will be useful, //
// but WITHOUT ANY WARRANTY; without even the implied warranty of //
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the //
// GNU General Public License for more details. //
// //
// You should have received a copy of the GNU General Public License //
// along with Pawian. If not, see <http://www.gnu.org/licenses/>. //
// //
//************************************************************************//
#include "Event/MassRangeCut.hh"
#include "Particle/Particle.hh"
MassRangeCut::MassRangeCut(std::string& rangeAndParticleNames, std::vector<Particle*>& finalStateParticles) :
_massMin(0.),
_massMax(0.)
{
// required order for the string: massMin massMax particleName1 particleName2 ...
std::vector<Particle*>::iterator itParticle;
int counter=0;
std::stringstream stringStrMassRange;
stringStrMassRange << rangeAndParticleNames;
std::string tmpNameMassRange;
while(stringStrMassRange >> tmpNameMassRange){
Info <<"\nMassRangeCont:\t" << tmpNameMassRange << endmsg;
if (counter==0) _massMin=atof(tmpNameMassRange.c_str());
else if(counter==1) _massMax=atof(tmpNameMassRange.c_str());
else{
unsigned int index=0;
for (itParticle=finalStateParticles.begin(); itParticle != finalStateParticles.end(); ++itParticle){
// Info << "tmpNameMassRange: " << tmpNameMassRange << "\t(*itParticle)->name(): " << (*itParticle)->name() <<endmsg;
if( tmpNameMassRange==(*itParticle)->name()){
// Info << "particle " << tmpNameMassRange << " added to list of mass range cuts " <<endmsg;
_particlesMassRange.push_back(*itParticle);
_iDparticlesMassRange.push_back(index);
break;
}
++index;
}
}
++counter;
}
if(_massMin<0. || _massMin > _massMax){
Alert << "massMin<0. or massMin>massMax !!!!!" << endmsg;
exit(0);
}
if(_particlesMassRange.size()<2){
Alert << "It is required to cut on an invariant mass of at least 2 final state particles!"
<< "\n here particlesMassRange.size()= " << _particlesMassRange.size() << endmsg;
exit(0);
}
Info << "mass cut on inv mass of the final state particles:" <<endmsg;
for (itParticle=_particlesMassRange.begin(); itParticle!=_particlesMassRange.end(); ++itParticle){
Info << (*itParticle)->name() << endmsg;
}
Info << "\tmassMin = " << _massMin
<< "\tmassmax = " << _massMax
<< endmsg;
}
MassRangeCut::~MassRangeCut()
{
}
bool MassRangeCut::isMassRangeParticle(Particle* compParticle){
bool result=false;
std::vector<Particle*>::iterator it;
for(it=_particlesMassRange.begin(); it!=_particlesMassRange.end(); ++it){
if ((*it)->sameName( *compParticle )){
result=true;
break;
}
}
return result;
}