QwPromptSummary Class Reference

#include <QwPromptSummary.h>

Inheritance diagram for QwPromptSummary:

Collaboration diagram for QwPromptSummary:

Public Member Functions
	QwPromptSummary ()
	QwPromptSummary (Int_t run_number, Int_t runlet_number)
	QwPromptSummary (Int_t run_number, Int_t runlet_number, const std::string &parameter_file)
virtual	~QwPromptSummary ()
void	SetRunNumber (const Int_t in)
Int_t	GetRunNumber ()
void	SetRunletNumber (const Int_t in)
Int_t	GetRunletNumber ()
void	SetPatternSize (const Int_t in)
Int_t	GetPatternSize ()
void	AddElement (PromptSummaryElement *in)
PromptSummaryElement *	GetElementByName (TString name)
void	FillDataInElement (TString name, Double_t yield, Double_t yield_err, Double_t yield_width, TString yield_unit, Double_t asym_diff, Double_t asym_diff_err, Double_t asym_diff_width, TString asym_diff_unit)
void	FillYieldToElement (TString name, Double_t yield, Double_t yield_error, Double_t yield_width, TString yield_unit)
void	FillAsymDiffToElement (TString name, Double_t asym_diff, Double_t asym_diff_err, Double_t asym_diff_width, TString asym_diff_unit)
void	FillDoubleDifference (TString type, TString name1, TString name2)
Int_t	GetSize () const
Int_t	Size () const
Int_t	HowManyElements () const
void	PrintCSV (Int_t nEvents, TString start_time, TString end_time)
void	PrintTextSummary ()

Data Fields
std::map< TString, PromptSummaryElement * >	fElementList
PromptSummaryElement *	fReferenceElement {nullptr}

Private Member Functions
TString	PrintTextSummaryHeader ()
TString	PrintTextSummaryTailer ()
TString	PrintCSVHeader (Int_t nEvents, TString start_time, TString end_time)
void	SetupElementList ()
void	LoadElementsFromParameterFile (const std::string &parameter_file)
void	LoadElementsFromParameterFile (QwParameterFile &parameterfile)
	ClassDefOverride (QwPromptSummary, 0)

Private Attributes
Int_t	fPatternSize
Int_t	fRunNumber
Int_t	fRunletNumber
Bool_t	fLocalDebug

Detailed Description

Definition at line 132 of file QwPromptSummary.h.

Constructor & Destructor Documentation

QwPromptSummary() [1/3]

QwPromptSummary::QwPromptSummary

(

)

Definition at line 179 of file QwPromptSummary.cc.

{
  fRunNumber    = 0;
  fRunletNumber = 0;
 
  fLocalDebug = kTRUE;
 
  this->SetupElementList();
 
};

References fLocalDebug, fRunletNumber, fRunNumber, and SetupElementList().

Referenced by ClassDefOverride().

Here is the call graph for this function:

Here is the caller graph for this function:

QwPromptSummary() [2/3]

QwPromptSummary::QwPromptSummary	(	Int_t	run_number,
		Int_t	runlet_number )

Definition at line 191 of file QwPromptSummary.cc.

{
  fRunNumber    = run_number;
  fRunletNumber = runlet_number;
 
  fLocalDebug = kFALSE;
 
  this->SetupElementList();
 
};

References fLocalDebug, fRunletNumber, fRunNumber, and SetupElementList().

Here is the call graph for this function:

QwPromptSummary() [3/3]

QwPromptSummary::QwPromptSummary	(	Int_t	run_number,
		Int_t	runlet_number,
		const std::string &	parameter_file )

Definition at line 203 of file QwPromptSummary.cc.

{
  fRunNumber    = run_number;
  fRunletNumber = runlet_number;
 
  fLocalDebug = kFALSE;
 
  this->LoadElementsFromParameterFile(parameter_file);
 
};

References fLocalDebug, fRunletNumber, fRunNumber, and LoadElementsFromParameterFile().

Here is the call graph for this function:

~QwPromptSummary()

QwPromptSummary::~QwPromptSummary ( )

virtual

Definition at line 215 of file QwPromptSummary.cc.

{
 
  for (auto i=fElementList.begin() ; i!=fElementList.end();i++){
    delete i->second;
  }
  fElementList.clear();
 
};

References fElementList.

Member Function Documentation

AddElement()

void QwPromptSummary::AddElement

(

PromptSummaryElement *

in

)

Definition at line 330 of file QwPromptSummary.cc.

{
  TString name = in->GetElementName(); name.ToLower();
  fElementList[name] = in;
  if(fLocalDebug) {
    printf("AddElement %s at pos %lu\n", in->GetElementName().Data(), fElementList.size()-1);
  }
  if (fElementList.size() == 1) {
    QwMessage << "First element: " << in->GetElementName() << " will be reference" << QwLog::endl;
    fReferenceElement = in;
  }
};

References QwLog::endl(), fElementList, fLocalDebug, fReferenceElement, PromptSummaryElement::GetElementName(), and QwMessage.

Referenced by LoadElementsFromParameterFile(), QwBeamLine::WritePromptSummary(), VQwDataHandler::WritePromptSummary(), and VQwDetectorArray::WritePromptSummary().

Here is the call graph for this function:

Here is the caller graph for this function:

ClassDefOverride()

QwPromptSummary::ClassDefOverride ( QwPromptSummary , 0  )

private

References QwPromptSummary().

Here is the call graph for this function:

FillAsymDiffToElement()

void QwPromptSummary::FillAsymDiffToElement	(	TString	name,
		Double_t	asym_diff,
		Double_t	asym_diff_err,
		Double_t	asym_diff_width,
		TString	asym_diff_unit )

Definition at line 480 of file QwPromptSummary.cc.

{
  PromptSummaryElement* an_element = NULL;
  an_element = this->GetElementByName(name);
  if(an_element) {
    an_element->SetAsymmetry(asym_diff);
    an_element->SetAsymmetryError(asym_diff_err);
    an_element->SetAsymmetryWidth(asym_diff_width);
  }
  else {
    if(fLocalDebug) {
      std::cout
    << "QwPromptSummary::FillYieldToElement : No Element with the name "
    << name
    <<  std::endl;
    }
  }
 
  return;
};

References fLocalDebug, GetElementByName(), PromptSummaryElement::SetAsymmetry(), PromptSummaryElement::SetAsymmetryError(), and PromptSummaryElement::SetAsymmetryWidth().

Here is the call graph for this function:

FillDataInElement()

void QwPromptSummary::FillDataInElement	(	TString	name,
		Double_t	yield,
		Double_t	yield_err,
		Double_t	yield_width,
		TString	yield_unit,
		Double_t	asym_diff,
		Double_t	asym_diff_err,
		Double_t	asym_diff_width,
		TString	asym_diff_unit )

Definition at line 424 of file QwPromptSummary.cc.

{
 
  PromptSummaryElement* an_element = NULL;
  an_element = this->GetElementByName(name);
  if(an_element) {
    an_element->FillData(yield, yield_err, yield_width, yield_unit, asym_diff, asym_diff_err, asym_diff_width, asym_diff_unit);
  }
  else {
    if(fLocalDebug) {
      std::cout
    << "QwPromptSummary::FillDataInElement : No Element with the name "
    << name
    <<  std::endl;
    }
  }
  return;
};

References PromptSummaryElement::FillData(), fLocalDebug, and GetElementByName().

Here is the call graph for this function:

FillDoubleDifference()

void QwPromptSummary::FillDoubleDifference	(	TString	type,
		TString	name1,
		TString	name2 )

Definition at line 507 of file QwPromptSummary.cc.

{
  PromptSummaryElement* an_element = NULL;
 
 
  PromptSummaryElement* one_element = NULL;
  PromptSummaryElement* two_element = NULL;
 
  one_element = this->GetElementByName(name1);
  two_element = this->GetElementByName(name2);
 
  if(one_element && two_element ) {
 
    an_element = this->GetElementByName(name1+"-"+name2);
 
    if(an_element) {
 
 
      Double_t diff       = 0.0;
      Double_t error_diff = 0.0;
      Double_t width_diff = 0.0;
 
      Double_t a = 0.0;
      Double_t b = 0.0;
      Double_t a_err = 0.0;
      Double_t b_err = 0.0;
      Double_t a_wit = 0.0;
      Double_t b_wit = 0.0;
 
 
      if(type.Contains("yield")) {
    a     = one_element -> GetYield();
    b     = two_element -> GetYield();
    a_err = one_element -> GetYieldError();
    b_err = two_element -> GetYieldError();
    a_wit = one_element -> GetYieldWidth();
    b_wit = two_element -> GetYieldWidth();
 
    diff       = a - b;
    error_diff = TMath::Sqrt(a_err*a_err + b_err*b_err);
    width_diff = a_wit - b_wit;
 
    an_element -> SetYield(diff);
    an_element -> SetYieldError(error_diff);
    an_element -> SetYieldWidth(width_diff);
 
      } else {
    a     = one_element -> GetAsymmetry();
    b     = two_element -> GetAsymmetry();
    a_err = one_element -> GetAsymmetryError();
    b_err = two_element -> GetAsymmetryError();
    a_wit = one_element -> GetAsymmetryWidth();
    b_wit = two_element -> GetAsymmetryWidth();
 
    diff       = a - b;
    error_diff = TMath::Sqrt(a_err*a_err + b_err*b_err);
    width_diff = a_wit - b_wit;
 
    an_element -> SetAsymmetry(diff);
    an_element -> SetAsymmetryError(error_diff);
    an_element -> SetAsymmetryWidth(width_diff);
 
 
      }
 
    }
    else {
      if(fLocalDebug) {
    std::cout
      << "QwPromptSummary::FillYieldToElement : No Element with the name "
      << name1 + "-" + name2
      <<  std::endl;
      }
  }
 
  }
  else {
    if(fLocalDebug) {
      std::cout
    << "QwPromptSummary::FillDoubleDifference: No Elements with the name "
    << name1 << " and " << name2
    <<  std::endl;
    }
  }
}

References fLocalDebug, and GetElementByName().

Here is the call graph for this function:

FillYieldToElement()

void QwPromptSummary::FillYieldToElement	(	TString	name,
		Double_t	yield,
		Double_t	yield_error,
		Double_t	yield_width,
		TString	yield_unit )

Definition at line 454 of file QwPromptSummary.cc.

{
  PromptSummaryElement* an_element = NULL;
  an_element = this->GetElementByName(name);
  if(an_element) {
    an_element->SetYield(yield);
    an_element->SetYieldError(yield_error);
    an_element->SetYieldWidth(yield_width);
  }
  else {
    if(fLocalDebug) {
      std::cout
    << "QwPromptSummary::FillYieldToElement : No Element with the name "
    << name
    <<  std::endl;
    }
  }
  return;
};

References fLocalDebug, GetElementByName(), PromptSummaryElement::SetYield(), PromptSummaryElement::SetYieldError(), and PromptSummaryElement::SetYieldWidth().

Here is the call graph for this function:

GetElementByName()

PromptSummaryElement * QwPromptSummary::GetElementByName

(

TString

name

)

Definition at line 345 of file QwPromptSummary.cc.

{
  name.ToLower();
  if (fElementList.find(name) != fElementList.end()) {
    PromptSummaryElement* an_element = fElementList[name];
    if(fLocalDebug) {
      std::cout << "System " << an_element->GetElementName()
                << " QwPromptSummary::GetElementByName address at " << an_element << std::endl;
    }
    return an_element;
  } else {
      QwDebug << "System " << name
              << " QwPromptSummary::GetElementByName not found" << std::endl;
  }
  return NULL;
};

References fElementList, fLocalDebug, PromptSummaryElement::GetElementName(), and QwDebug.

Referenced by FillAsymDiffToElement(), FillDataInElement(), FillDoubleDifference(), FillYieldToElement(), QwBeamLine::WritePromptSummary(), VQwDataHandler::WritePromptSummary(), and VQwDetectorArray::WritePromptSummary().

Here is the call graph for this function:

Here is the caller graph for this function:

GetPatternSize()

Int_t QwPromptSummary::GetPatternSize ( )

inline

Definition at line 151 of file QwPromptSummary.h.

{ return fPatternSize; };

References fPatternSize.

GetRunletNumber()

Int_t QwPromptSummary::GetRunletNumber ( )

inline

Definition at line 148 of file QwPromptSummary.h.

{return fRunletNumber;};

References fRunletNumber.

GetRunNumber()

Int_t QwPromptSummary::GetRunNumber ( )

inline

Definition at line 145 of file QwPromptSummary.h.

{return fRunNumber;};

References fRunNumber.

GetSize()

Int_t QwPromptSummary::GetSize ( ) const

inline

Definition at line 168 of file QwPromptSummary.h.

{return fElementList.size();};

References fElementList.

HowManyElements()

Int_t QwPromptSummary::HowManyElements ( ) const

inline

Definition at line 170 of file QwPromptSummary.h.

{return fElementList.size();};

References fElementList.

LoadElementsFromParameterFile() [1/2]

void QwPromptSummary::LoadElementsFromParameterFile ( const std::string & parameter_file )

private

Definition at line 250 of file QwPromptSummary.cc.

{
  try {
    QwParameterFile paramfile(parameter_file);
    LoadElementsFromParameterFile(paramfile);
  } catch (const std::exception& e) {
    QwError << "QwPromptSummary::LoadElementsFromParameterFile: Unable to open parameter file: "
            << parameter_file << " - " << e.what() << QwLog::endl;
    QwMessage << "Falling back to default (empty) element list." << QwLog::endl;
  }
};

References QwLog::endl(), LoadElementsFromParameterFile(), QwError, and QwMessage.

Referenced by LoadElementsFromParameterFile(), QwPromptSummary(), and SetupElementList().

Here is the call graph for this function:

Here is the caller graph for this function:

LoadElementsFromParameterFile() [2/2]

void QwPromptSummary::LoadElementsFromParameterFile ( QwParameterFile & parameterfile )

private

Definition at line 264 of file QwPromptSummary.cc.

{
  QwMessage << "QwPromptSummary::LoadElementsFromParameterFile: Loading prompt summary elements" << QwLog::endl;
 
  // Read preamble
  std::unique_ptr<QwParameterFile> preamble = parameterfile.ReadSectionPreamble();
  if (preamble) {
    QwVerbose << "PromptSummary preamble:" << QwLog::endl;
    QwVerbose << *preamble << QwLog::endl;
  }
 
  // Read sections
  std::unique_ptr<QwParameterFile> section;
  std::string section_name;
  while ((section = parameterfile.ReadNextSection(section_name))) {
    QwVerbose << "Processing section: " << section_name << QwLog::endl;
 
    // Check if this is a prompt summary elements section
    if (section_name == "prompt_summary_elements" || section_name == "elements") {
      // Process individual elements in this section
      while (section->ReadNextLine()) {
        section->TrimWhitespace();
        section->TrimComment();
        if (section->LineIsEmpty()) continue;
 
        std::string line = section->GetLine();
 
        // Parse element definitions
        // Format 1: Simple element name
        // Format 2: element_name = type (where type could be single, difference, etc.)
 
        std::string element_name, element_type;
        if (section->HasVariablePair("=", element_name, element_type)) {
          // Format: element_name = type
          // Trim whitespace manually since TrimWhitespace is protected
          element_name.erase(0, element_name.find_first_not_of(" \t\r\n"));
          element_name.erase(element_name.find_last_not_of(" \t\r\n") + 1);
          element_type.erase(0, element_type.find_first_not_of(" \t\r\n"));
          element_type.erase(element_type.find_last_not_of(" \t\r\n") + 1);
 
          if (fLocalDebug) {
            QwMessage << "Adding element: " << element_name << " (type: " << element_type << ")" << QwLog::endl;
          }
        } else {
          // Format: simple element name
          element_name = line;
          element_type = "single"; // default type
 
          if (fLocalDebug) {
            QwMessage << "Adding element: " << element_name << " (default type)" << QwLog::endl;
          }
        }
 
        // Create and add the element
        if (!element_name.empty()) {
          this->AddElement(new PromptSummaryElement(TString(element_name.c_str())));
        }
      }
    }
  }
 
  QwMessage << "QwPromptSummary: Loaded " << fElementList.size() << " elements from parameter file" << QwLog::endl;
};

References AddElement(), QwLog::endl(), fElementList, fLocalDebug, QwMessage, QwVerbose, QwParameterFile::ReadNextSection(), and QwParameterFile::ReadSectionPreamble().

Here is the call graph for this function:

PrintCSV()

void QwPromptSummary::PrintCSV	(	Int_t	nEvents,
		TString	start_time,
		TString	end_time )

Definition at line 594 of file QwPromptSummary.cc.

{
  printf("-----------------------\n");
  TString filename = gQwOptions.GetValue<std::string>("rootfiles");
  filename+=Form("/summary_%d.txt", fRunNumber);
  TString header= this->PrintCSVHeader(nEvents, start_time, end_time);
  std::ofstream output;
  output.open(filename.Data());
  output<< header.Data();
 
  TString secheader= "=========================================================================\n";
  secheader+=Form("%40s \n","Yields");
  secheader+="=========================================================================\n" ;
  output << secheader.Data() ;
 
  for (auto i=fElementList.begin(); i!=fElementList.end(); i++  )
    {
      output << i->second->GetCSVSummary("yield") ;
    }
 
 
  secheader= "=========================================================================\n";
  secheader+=Form("%50s\n","Asymmetries/Differences");
  secheader+="=========================================================================\n";
  output << secheader.Data();
 
  for ( auto j=fElementList.begin(); j!=fElementList.end(); j++ )
    {
      output << j->second->GetCSVSummary("asymmetry");
    }
 
 
 
  secheader= "=========================================================================\n";
  secheader+=Form("%55s\n", "Combined Differences/Averages");
  secheader+="=========================================================================\n";
  output << secheader.Data();
 
  for ( auto j=fElementList.begin(); j!=fElementList.end(); j++ )
    {
      output << j->second->GetCSVSummary("double");
    }
 
  output<< "=========================================================================\n";
  output<< Form("%45s\n"," End of Summary");
  output<< "=========================================================================\n";
  output.close();
 
  return;
};

References fElementList, fRunNumber, gQwOptions, and PrintCSVHeader().

Referenced by main().

Here is the call graph for this function:

Here is the caller graph for this function:

PrintCSVHeader()

TString QwPromptSummary::PrintCSVHeader ( Int_t nEvents, TString start_time, TString end_time )

private

Definition at line 389 of file QwPromptSummary.cc.

{
  TString out = "";
 
  Double_t goodEvents = 0.0;
  TString referenceElementName = "N/A";
 
  // Use the first element in the list to determine good events
  if (fReferenceElement) {
    goodEvents = fReferenceElement->GetNumGoodEvents() * fPatternSize;
    referenceElementName = fReferenceElement->GetElementName();
  } else {
    QwError << "Warning: No reference element found in QwPromptSummary. Setting goodEvents=0." << QwLog::endl;
  }
 
  out += Form("Run: %d \n",fRunNumber);
  out += "Start Time: "+start_time+"\nEnd Time: "+end_time+"\n";
  out += Form("Number of events processed: %i\n",nEvents);
  out += Form("Number of events in good multiplicity patterns: %3.0f\n", goodEvents);
  out += Form("Percentage of good events: %3.1f %%\n", goodEvents/nEvents*100);
  out += Form("Good events reference: %s (first element from parameter file)\n", referenceElementName.Data());
  out += "=========================================================================\n";
  out += "Yield Units: bcm(uA), cavq(uA), bpm(mm), sam(mV/uA)\n";
  out += "Asymmetry/Difference Units: bcm(ppm), cavq(ppm), bpm(um), sam(ppm)\n";
 
  out += "=========================================================================\n";
 
 
 
  return out;
};

References QwLog::endl(), fPatternSize, fReferenceElement, fRunNumber, and QwError.

Referenced by PrintCSV().

Here is the call graph for this function:

Here is the caller graph for this function:

PrintTextSummary()

void QwPromptSummary::PrintTextSummary

(

)

Definition at line 646 of file QwPromptSummary.cc.

{
  return;
};

PrintTextSummaryHeader()

TString QwPromptSummary::PrintTextSummaryHeader ( )

private

Definition at line 364 of file QwPromptSummary.cc.

{
  TString out = "";
  //  TString filename = "";
 
  //  filename = Form("summary_%d_%d.txt", fRunNumber, fRunletNumber);
 
  out = "======= BEGIN ======= \n";
  out += Form("        RUN = %5d RUNLET %2d \n\n\n\n\n\n", fRunNumber, fRunletNumber);
  return out;
};

References fRunletNumber, and fRunNumber.

PrintTextSummaryTailer()

TString QwPromptSummary::PrintTextSummaryTailer ( )

private

Definition at line 379 of file QwPromptSummary.cc.

{
  TString out = "";
 
  out =  "==================================================\n";
  return out;
};

SetPatternSize()

void QwPromptSummary::SetPatternSize ( const Int_t in )

inline

Definition at line 150 of file QwPromptSummary.h.

{ fPatternSize=in; };

References fPatternSize.

Referenced by QwHelicityPattern::WritePromptSummary().

Here is the caller graph for this function:

SetRunletNumber()

void QwPromptSummary::SetRunletNumber ( const Int_t in )

inline

Definition at line 147 of file QwPromptSummary.h.

{fRunletNumber = in;};

References fRunletNumber.

SetRunNumber()

void QwPromptSummary::SetRunNumber ( const Int_t in )

inline

Definition at line 144 of file QwPromptSummary.h.

{fRunNumber = in;};

References fRunNumber.

SetupElementList()

void QwPromptSummary::SetupElementList ( )

private

Definition at line 228 of file QwPromptSummary.cc.

{
  // Try to load from default parameter file first
  std::string default_param_file = "prompt_summary.map";
 
  try {
    QwParameterFile paramfile(default_param_file);
    LoadElementsFromParameterFile(paramfile);
    if (fElementList.size() > 0) {
      QwMessage << "QwPromptSummary: Loaded " << fElementList.size() << " elements from " << default_param_file << QwLog::endl;
      return; // Successfully loaded from file
    }
  } catch (const std::exception& e) {
    if (fLocalDebug) {
      QwMessage << "QwPromptSummary: Could not load from " << default_param_file
                << ", using default elements: " << e.what() << QwLog::endl;
    }
  }
};

References QwLog::endl(), fElementList, fLocalDebug, LoadElementsFromParameterFile(), and QwMessage.

Referenced by QwPromptSummary(), and QwPromptSummary().

Here is the call graph for this function:

Here is the caller graph for this function:

Size()

Int_t QwPromptSummary::Size ( ) const

inline

Definition at line 169 of file QwPromptSummary.h.

{return fElementList.size();};

References fElementList.

Field Documentation

fElementList

std::map<TString, PromptSummaryElement*> QwPromptSummary::fElementList

Definition at line 141 of file QwPromptSummary.h.

Referenced by AddElement(), GetElementByName(), GetSize(), HowManyElements(), LoadElementsFromParameterFile(), PrintCSV(), SetupElementList(), Size(), and ~QwPromptSummary().

fLocalDebug

Bool_t QwPromptSummary::fLocalDebug

private

Definition at line 192 of file QwPromptSummary.h.

Referenced by AddElement(), FillAsymDiffToElement(), FillDataInElement(), FillDoubleDifference(), FillYieldToElement(), GetElementByName(), LoadElementsFromParameterFile(), QwPromptSummary(), QwPromptSummary(), QwPromptSummary(), and SetupElementList().

fPatternSize

Int_t QwPromptSummary::fPatternSize

private

Definition at line 179 of file QwPromptSummary.h.

Referenced by GetPatternSize(), PrintCSVHeader(), and SetPatternSize().

fReferenceElement

PromptSummaryElement* QwPromptSummary::fReferenceElement {nullptr}

Definition at line 142 of file QwPromptSummary.h.

{nullptr};

Referenced by AddElement(), and PrintCSVHeader().

fRunletNumber

Int_t QwPromptSummary::fRunletNumber

private

Definition at line 190 of file QwPromptSummary.h.

Referenced by GetRunletNumber(), PrintTextSummaryHeader(), QwPromptSummary(), QwPromptSummary(), QwPromptSummary(), and SetRunletNumber().

fRunNumber

Int_t QwPromptSummary::fRunNumber

private

Definition at line 189 of file QwPromptSummary.h.

Referenced by GetRunNumber(), PrintCSV(), PrintCSVHeader(), PrintTextSummaryHeader(), QwPromptSummary(), QwPromptSummary(), QwPromptSummary(), and SetRunNumber().

---

The documentation for this class was generated from the following files:

• QwPromptSummary.h
• QwPromptSummary.cc