Skip to content
Snippets Groups Projects
Commit c13cb14b authored by Michael Leyhe's avatar Michael Leyhe
Browse files

BugFix: Initialized cmsMass in epemParser to jpsi-mass

parent 060e5f12
No related branches found
No related tags found
No related merge requests found
......@@ -67,7 +67,6 @@ Meson f2Hat(1800) \hat{f}_2(1800) 0 4 1 -1 1 0 0 0 0 1.800 0.3000
#
# type name texName 3*Q 2*J P C G I I3 S C mass width dyn.f comment
#
Meson pion \pi 0 0 -1 1 -1 1 0 0 0 0.139 0.0001 relBW [PDG nonett 1^1S_0]
Meson pion0 \pi^0 0 0 -1 1 -1 1 0 0 0 0.13498 0.0001 relBW [PDG nonett 1^1S_0]
Meson pion+ \pi^+ 3 0 -1 0 -1 1 1 0 0 0.13957 0.0001 relBW [PDG nonett 1^1S_0]
Meson pion- \pi^- -3 0 -1 0 -1 1 -1 0 0 0.13957 0.0001 relBW [PDG nonett 1^1S_0]
......
......@@ -70,6 +70,7 @@ void epemHist::initRootStuff(){
double massMin = 0;
// double massMax = psiMass;
double massMax = epemEnv::instance()->cmsMass();
DebugMsg << "epemEnv::instance()->cmsMass():\t" << epemEnv::instance()->cmsMass() << endmsg;
std::vector<std::string> fspNames=tmpMassHistData->_fspNames;
std::vector<Particle*> allFsp = epemEnv::instance()->finalStateParticles();
......@@ -91,21 +92,27 @@ void epemHist::initRootStuff(){
massMax += (massMax - massMin) * 0.02;
massMin -= (massMax - massMin) * 0.02;
TH1F* currentMassDataHist=new TH1F(histName.c_str(), histDescription.c_str(), 100., massMin, massMax);
TH1F* currentMassDataHist=new TH1F(histName.c_str(), histDescription.c_str(), 100, massMin, massMax);
currentMassDataHist->Sumw2();
_massDataHistMap[tmpMassHistData]=currentMassDataHist;
DebugMsg << histName.c_str() << "\t" << massMin << "\t" << massMax << endmsg;
histName="MC"+tmpBaseName;
histDescription = "M("+tmpMassHistData->_name+") (MC)";
TH1F* currentMassMcHist=new TH1F(histName.c_str(), histDescription.c_str(), 100., massMin, massMax);
TH1F* currentMassMcHist=new TH1F(histName.c_str(), histDescription.c_str(), 100, massMin, massMax);
currentMassMcHist->Sumw2();
_massMcHistMap[tmpMassHistData]=currentMassMcHist;
DebugMsg << histName.c_str() << "\t" << massMin << "\t" << massMax << endmsg;
histName="Fit"+tmpBaseName;
histDescription = "M("+tmpMassHistData->_name+") (fit)";
TH1F* currentMassFitHist=new TH1F(histName.c_str(), histDescription.c_str(), 100., massMin, massMax);
TH1F* currentMassFitHist=new TH1F(histName.c_str(), histDescription.c_str(), 100, massMin, massMax);
currentMassFitHist->Sumw2();
_massFitHistMap[tmpMassHistData]=currentMassFitHist;
DebugMsg << histName.c_str() << "\t" << massMin << "\t" << massMax << endmsg;
}
}
......
......@@ -36,6 +36,7 @@ using namespace std;
epemParser::epemParser(int argc,char **argv):
ParserBase(argc,argv)
,_cmsMass(3.096916)
{
po::options_description common("Common Options");
common.add_options()
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment