Skip to content
Snippets Groups Projects
Commit ea874a40 authored by Bertram Kopf's avatar Bertram Kopf
Browse files

few modifications in BwShape.cc

parent 5bbef7b7
No related branches found
No related tags found
No related merge requests found
...@@ -55,20 +55,20 @@ BwShape::BwShape(double MassRes, double MassWidth, double MassDec1, double MassD ...@@ -55,20 +55,20 @@ BwShape::BwShape(double MassRes, double MassWidth, double MassDec1, double MassD
std::stringstream Lstrstr; std::stringstream Lstrstr;
Lstrstr << i; Lstrstr << i;
std::string histName="BreitWigner_L"+Lstrstr.str(); // std::string histName="BreitWigner_L"+Lstrstr.str();
_histMap[i]= new TH1F(histName.c_str(),histName.c_str(),301, MassRes-deltaMass, MassRes+deltaMass); // _histMap[i]= new TH1F(histName.c_str(),histName.c_str(),301, MassRes-deltaMass, MassRes+deltaMass);
std::string histNameNew="BreitWigner_Lnew"+Lstrstr.str(); std::string histNameNew="BreitWigner_Lnew"+Lstrstr.str();
_histMapNew[i]= new TH1F(histNameNew.c_str(),histNameNew.c_str(),301, MassRes-deltaMass, MassRes+deltaMass); _histMapNew[i]= new TH1F(histNameNew.c_str(),histNameNew.c_str(),301, MassRes-deltaMass, MassRes+deltaMass);
histName="Argand_L"+Lstrstr.str(); // histName="Argand_L"+Lstrstr.str();
TH2F* currentArgandHist=new TH2F(histName.c_str(),histName.c_str(),301, -1., 1., 301, 0., 1.3); // TH2F* currentArgandHist=new TH2F(histName.c_str(),histName.c_str(),301, -1., 1., 301, 0., 1.3);
currentArgandHist->SetXTitle("Re(Bw)"); // currentArgandHist->SetXTitle("Re(Bw)");
currentArgandHist->SetYTitle("Im(Bw)"); // currentArgandHist->SetYTitle("Im(Bw)");
currentArgandHist->SetMarkerStyle(6); // currentArgandHist->SetMarkerStyle(6);
_argandHistMap[i]=currentArgandHist; // _argandHistMap[i]=currentArgandHist;
histNameNew="Argand_Lnew"+Lstrstr.str(); histNameNew="Argand_Lnew"+Lstrstr.str();
TH2F* currentArgandHistNew=new TH2F(histNameNew.c_str(),histNameNew.c_str(),301, -1., 1., 301, 0., 1.3); TH2F* currentArgandHistNew=new TH2F(histNameNew.c_str(),histNameNew.c_str(),301, -1., 1., 301, 0., 1.3);
...@@ -84,9 +84,9 @@ BwShape::BwShape(double MassRes, double MassWidth, double MassDec1, double MassD ...@@ -84,9 +84,9 @@ BwShape::BwShape(double MassRes, double MassWidth, double MassDec1, double MassD
for (unsigned int lIt=0; lIt<=Lmax; ++lIt){ for (unsigned int lIt=0; lIt<=Lmax; ++lIt){
TH1F* currentHist=_histMap[lIt]; // TH1F* currentHist=_histMap[lIt];
TH1F* currentHistNew=_histMapNew[lIt]; TH1F* currentHistNew=_histMapNew[lIt];
TH2F* currentArgandHist=_argandHistMap[lIt]; // TH2F* currentArgandHist=_argandHistMap[lIt];
TH2F* currentArgandHistNew=_argandHistMapNew[lIt]; TH2F* currentArgandHistNew=_argandHistMapNew[lIt];
...@@ -94,13 +94,13 @@ BwShape::BwShape(double MassRes, double MassWidth, double MassDec1, double MassD ...@@ -94,13 +94,13 @@ BwShape::BwShape(double MassRes, double MassWidth, double MassDec1, double MassD
Vector4<double> res4V(massIt, 0., 0., 0.); Vector4<double> res4V(massIt, 0., 0., 0.);
complex<double> currentBW=BreitWignerBlattW(res4V, MassDec1, MassDec2, MassRes, MassWidth, lIt); // complex<double> currentBW=BreitWignerBlattW(res4V, MassDec1, MassDec2, MassRes, MassWidth, lIt);
complex<double> currentBWnew=BreitWignerFunction::BlattWRel(lIt, res4V.Mass(), MassRes, MassWidth, MassDec1, MassDec2); complex<double> currentBWnew=BreitWignerFunction::BlattWRel(lIt, res4V.Mass(), MassRes, MassWidth, MassDec1, MassDec2);
double weight=norm(currentBW); // double weight=norm(currentBW);
currentHist->Fill(massIt,weight); // currentHist->Fill(massIt,weight);
currentHistNew->Fill(massIt,norm(currentBWnew)); currentHistNew->Fill(massIt,norm(currentBWnew));
currentArgandHist->Fill(currentBW.real(),currentBW.imag()); // currentArgandHist->Fill(currentBW.real(),currentBW.imag());
currentArgandHistNew->Fill(currentBWnew.real(),currentBWnew.imag()); currentArgandHistNew->Fill(currentBWnew.real(),currentBWnew.imag());
......
...@@ -60,11 +60,11 @@ protected: ...@@ -60,11 +60,11 @@ protected:
private: private:
TFile* _theTFile; TFile* _theTFile;
std::map <unsigned int, TH1F* > _histMap; // std::map <unsigned int, TH1F* > _histMap;
std::map <unsigned int, TH1F* > _histMapNew; std::map <unsigned int, TH1F* > _histMapNew;
std::map <unsigned int, TH2F* > _argandHistMap; // std::map <unsigned int, TH2F* > _argandHistMap;
std::map <unsigned int, TH2F* > _argandHistMapNew; std::map <unsigned int, TH2F* > _argandHistMapNew;
std::map <unsigned int, TH2F* > _phaseHistMap; // std::map <unsigned int, TH2F* > _phaseHistMap;
}; };
#endif #endif
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