Extracts dE/dx information for calibration testing. More...

#include <HitLevelInfoWriter.h>

Inheritance diagram for HitLevelInfoWriterModule:

Public Types
enum	EModulePropFlags { c_Input = 1 , c_Output = 2 , c_ParallelProcessingCertified = 4 , c_HistogramManager = 8 , c_InternalSerializer = 16 , c_TerminateInAllProcesses = 32 , c_DontCollectStatistics = 64 }
	Each module can be tagged with property flags, which indicate certain features of the module. More...

typedef ModuleCondition::EAfterConditionPath	EAfterConditionPath
	Forward the EAfterConditionPath definition from the ModuleCondition.

Public Member Functions
	HitLevelInfoWriterModule ()
	Default constructor.

virtual	~HitLevelInfoWriterModule ()
	Destructor.

virtual void	initialize () override
	Initialize the module.

virtual void	event () override
	This method is called for each event.

virtual void	terminate () override
	End of the event processing.

void	bookOutput (std::string filename)
	Create the output TFiles and TTrees.

void	recalculateDedx (CDCDedxTrack *dedxTrack, std::map< int, std::vector< double > > &l_var, double(&cdcChi)[Const::ChargedStable::c_SetSize])
	Function to recalculate the dedx with latest constants.

double	GetCorrection (int &adc, int layer, int wireID, double doca, double enta, double costheta, double ring, double time) const
	Function to get the correction factor.

void	HadronCorrection (double costheta, double &dedx) const
	Function to apply the hadron correction.

double	D2I (const double cosTheta, const double D) const
	hadron saturation parameterization part 2

double	I2D (const double cosTheta, const double I) const
	hadron saturation parameterization part 1

void	calculateMeans (double mean, double truncatedMean, double *truncatedMeanErr, const std::vector< double > &dedx) const
	Save arithmetic and truncated mean for the 'dedx' values.

void	saveChiValue (double(&chi)[Const::ChargedStable::c_SetSize], CDCDedxTrack *dedxTrack, double dedx) const
	for all particles, save chi values into 'chi' chi array of chi values to be modified

virtual std::vector< std::string >	getFileNames (bool outputFiles)
	Return a list of output filenames for this modules.

virtual void	beginRun ()
	Called when entering a new run.

virtual void	endRun ()
	This method is called if the current run ends.

const std::string &	getName () const
	Returns the name of the module.

const std::string &	getType () const
	Returns the type of the module (i.e.

const std::string &	getPackage () const
	Returns the package this module is in.

const std::string &	getDescription () const
	Returns the description of the module.

void	setName (const std::string &name)
	Set the name of the module.

void	setPropertyFlags (unsigned int propertyFlags)
	Sets the flags for the module properties.

LogConfig &	getLogConfig ()
	Returns the log system configuration.

void	setLogConfig (const LogConfig &logConfig)
	Set the log system configuration.

void	setLogLevel (int logLevel)
	Configure the log level.

void	setDebugLevel (int debugLevel)
	Configure the debug messaging level.

void	setAbortLevel (int abortLevel)
	Configure the abort log level.

void	setLogInfo (int logLevel, unsigned int logInfo)
	Configure the printed log information for the given level.

void	if_value (const std::string &expression, const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
	Add a condition to the module.

void	if_false (const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
	A simplified version to add a condition to the module.

void	if_true (const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
	A simplified version to set the condition of the module.

bool	hasCondition () const
	Returns true if at least one condition was set for the module.

const ModuleCondition *	getCondition () const
	Return a pointer to the first condition (or nullptr, if none was set)

const std::vector< ModuleCondition > &	getAllConditions () const
	Return all set conditions for this module.

bool	evalCondition () const
	If at least one condition was set, it is evaluated and true returned if at least one condition returns true.

std::shared_ptr< Path >	getConditionPath () const
	Returns the path of the last true condition (if there is at least one, else reaturn a null pointer).

Module::EAfterConditionPath	getAfterConditionPath () const
	What to do after the conditional path is finished.

std::vector< std::shared_ptr< Path > >	getAllConditionPaths () const
	Return all condition paths currently set (no matter if the condition is true or not).

bool	hasProperties (unsigned int propertyFlags) const
	Returns true if all specified property flags are available in this module.

bool	hasUnsetForcedParams () const
	Returns true and prints error message if the module has unset parameters which the user has to set in the steering file.

const ModuleParamList &	getParamList () const
	Return module param list.

template<typename T>
ModuleParam< T > &	getParam (const std::string &name) const
	Returns a reference to a parameter.

bool	hasReturnValue () const
	Return true if this module has a valid return value set.

int	getReturnValue () const
	Return the return value set by this module.

std::shared_ptr< PathElement >	clone () const override
	Create an independent copy of this module.

std::shared_ptr< boost::python::list >	getParamInfoListPython () const
	Returns a python list of all parameters.

Static Public Member Functions
static void	exposePythonAPI ()
	Exposes methods of the Module class to Python.

Protected Member Functions
virtual void	def_initialize ()
	Wrappers to make the methods without "def_" prefix callable from Python.

virtual void	def_beginRun ()
	Wrapper method for the virtual function beginRun() that has the implementation to be used in a call from Python.

virtual void	def_event ()
	Wrapper method for the virtual function event() that has the implementation to be used in a call from Python.

virtual void	def_endRun ()
	This method can receive that the current run ends as a call from the Python side.

virtual void	def_terminate ()
	Wrapper method for the virtual function terminate() that has the implementation to be used in a call from Python.

void	setDescription (const std::string &description)
	Sets the description of the module.

void	setType (const std::string &type)
	Set the module type.

template<typename T>
void	addParam (const std::string &name, T &paramVariable, const std::string &description, const T &defaultValue)
	Adds a new parameter to the module.

template<typename T>
void	addParam (const std::string &name, T &paramVariable, const std::string &description)
	Adds a new enforced parameter to the module.

void	setReturnValue (int value)
	Sets the return value for this module as integer.

void	setReturnValue (bool value)
	Sets the return value for this module as bool.

void	setParamList (const ModuleParamList &params)
	Replace existing parameter list.

Private Member Functions
void	fillTrack (const TrackFitResult *fitResult)
	Fill the TTree with the information from the track fit.

void	fillDedx (CDCDedxTrack *dedxTrack)
	Fill the TTree with the information from a CDCDedxTrack object.

void	clearEntries ()
	Clear the arrays before filling an event.

std::list< ModulePtr >	getModules () const override
	no submodules, return empty list

std::string	getPathString () const override
	return the module name.

void	setParamPython (const std::string &name, const boost::python::object &pyObj)
	Implements a method for setting boost::python objects.

void	setParamPythonDict (const boost::python::dict &dictionary)
	Implements a method for reading the parameter values from a boost::python dictionary.

Private Attributes
std::string	m_strOutputBaseName
	Base name for the output ROOT files.

std::vector< std::string >	m_strParticleList
	Vector of ParticleLists to write out.

std::vector< std::string >	m_filename
	full names of the output ROOT files

std::vector< TFile * >	m_file
	output ROOT files

std::vector< TTree * >	m_tree
	output ROOT trees

StoreArray< CDCDedxTrack >	m_dedxTracks
	Required array of CDCDedxTracks.

StoreArray< Track >	m_tracks
	Required array of input Tracks.

StoreArray< TrackFitResult >	m_trackFitResults
	Required array of input TrackFitResults.

StoreArray< ECLCluster >	m_eclClusters
	Required array of input ECLClusters.

StoreArray< KLMCluster >	m_klmClusters
	Required array of input KLMClusters.

StoreObjPtr< EventLevelTriggerTimeInfo >	m_TTDInfo
	Store Object Ptr: EventLevelTriggerTimeInfo.

int	m_expID { -1}
	experiment in which this Track was found

int	m_runID { -1}
	run in which this Track was found

int	m_eventID { -1}
	event in which this Track was found

double	m_injring { -1.0}
	HER injection status.

double	m_injtime { -1.0}
	time since last injection in micro seconds

double	m_d0 {0.}
	Signed distance to the POCA in the r-phi plane.

double	m_z0 {0.}
	z coordinate of the POCA

double	m_dz { -1.}
	vertex or POCA in case of tracks z in respect to IPs

double	m_dr { -1.}
	track d0 relative to IP

double	m_dphi { -1.}
	POCA in degrees in respect to IP.

double	m_vx0 {0.}
	X coordinate of track POCA to origin.

double	m_vy0 {0.}
	Y coordinate of track POCA to origin.

double	m_vz0 {0.}
	Z coordinate of track POCA to origin.

double	m_tanlambda { -1.}
	Slope of the track in the r-z plane.

double	m_phi0 { -1.}
	Angle of the transverse momentum in the r-phi plane.

double	m_chi2 { -1.}
	chi^2 from track fit

double	m_nCDChits { -1.}
	Number of CDC hits associated to the track.

int	m_inCDC { -1}
	frack is CDC acceptance or not

int	m_trackID { -1}
	ID number of the Track.

double	m_length { -1.}
	total path length of the Track

int	m_charge {0}
	the charge for this Track

double	m_cosTheta { -2.}
	cos(theta) for the track

double	m_pCDC { -1.}
	momentum valid in CDC

double	m_p { -1.}
	momentum from tracking

double	m_pt { -1.}
	transverse momentum from tracking

double	m_phi { -1.}
	phi for the track

double	m_ioasym { -1.}
	asymmetry in increasing vs decreasing layer numbers per track

double	m_theta { -2.}
	cos(theta) for the track

double	m_PDG { -1.}
	MC PID.

double	m_mean { -1.}
	dE/dx averaged

double	m_trunc { -1.}
	dE/dx averaged, truncated mean, with corrections

double	m_truncNoSat { -1.}
	dE/dx averaged, truncated mean, with corrections (not hadron)

double	m_error { -1.}
	standard deviation of the truncated mean

double	m_eop { -1.}
	energy over momentum in the calorimeter

double	m_e { -1.}
	energy in the calorimeter

double	m_e1_9 { -1.}
	ratio of energies of the central 1 crystal vs 3x3 crystals

double	m_e9_21 { -1.}
	ratio of energies of the central 3x3 crystal vs 5x5 crystals

double	m_klmLayers { -1.}
	number of klm layers with hits

double	m_eclsnHits { -1.}

double	m_scale { -1.}
	calibration scale factor

double	m_cosCor { -1.}
	calibration cosine correction

double	m_cosEdgeCor { -1.}
	calibration cosine edge correction

double	m_runGain { -1.}
	calibration run gain

double	m_timeGain { -1.}
	calibration injection time gain

double	m_timeReso { -1.}
	calibration injection time reso

double	m_chieOld { -1.}
	chi value for electron hypothesis

double	m_chimuOld { -1.}
	chi value for muon hypothesis

double	m_chipiOld { -1.}
	chi value for pion hypothesis

double	m_chikOld { -1.}
	chi value for kaon hypothesis

double	m_chipOld { -1.}
	chi value for proton hypothesis

double	m_chidOld { -1.}
	chi value for deuteron hypothesis

double	m_chie { -1.}
	modified chi value for electron hypothesis

double	m_chimu { -1.}
	modified chi value for muon hypothesis

double	m_chipi { -1.}
	modified chi value for pion hypothesis

double	m_chik { -1.}
	modified chi value for kaon hypothesis

double	m_chip { -1.}
	modified chi value for proton hypothesis

double	m_chid { -1.}
	modified chi value for deuteron hypothesis

double	m_prese { -1.}
	pred reso for electron hypothesis

double	m_presmu { -1.}
	pred reso for muon hypothesis

double	m_prespi { -1.}
	pred reso for pion hypothesis

double	m_presk { -1.}
	pred reso for kaon hypothesis

double	m_presp { -1.}
	pred reso for proton hypothesis

double	m_presd { -1.}
	pred reso for deuteron hypothesis

double	m_pmeane { -1.}
	pred mean for electron hypothesis

double	m_pmeanmu { -1.}
	pred mean for muon hypothesis

double	m_pmeanpi { -1.}
	pred mean for pion hypothesis

double	m_pmeank { -1.}
	pred mean for kaon hypothesis

double	m_pmeanp { -1.}
	pred mean for proton hypothesis

double	m_pmeand { -1.}
	pred mean for deuteron hypothesis

int	l_nhits { -1}
	the total number of layer hits for this Track

int	l_nhitsused { -1}
	the total number of layer hits used for this Track

int	l_nhitscombined [kMaxHits] = {}
	the number of hits combined this layer

int	l_wirelongesthit [kMaxHits] = {}
	the wire number of longest hit in this layer

int	l_layer [kMaxHits] = {}
	layer number

double	l_path [kMaxHits] = {}
	distance travelled in this layer

double	l_dedx [kMaxHits] = {}
	dE/dx for this layer

int	h_nhits { -1}
	the number of good hits for this Track

int	h_lwire [kMaxHits] = {}
	sense wire within layer

int	h_wire [kMaxHits] = {}
	sense wire ID

int	h_layer [kMaxHits] = {}
	layer number

double	h_path [kMaxHits] = {}
	path length in cell

double	h_dedx [kMaxHits] = {}
	charge per path length

double	h_adcraw [kMaxHits] = {}
	charge per hit

double	h_adccorr [kMaxHits] = {}
	charge per hit corr by nonlinear ADC

double	h_doca [kMaxHits] = {}
	distance of closest approach

double	h_ndoca [kMaxHits] = {}
	normalized distance of closest approach

double	h_ndocaRS [kMaxHits] = {}
	normalized +RS distance of closest approach

double	h_enta [kMaxHits] = {}
	entrance angle

double	h_entaRS [kMaxHits] = {}
	normalized + RS distance of entrance angle

double	h_driftT [kMaxHits] = {}
	drift time

double	h_driftD [kMaxHits] = {}
	drift distance

double	h_facnladc [kMaxHits] = {}
	calibration hit gain

double	h_wireGain [kMaxHits] = {}
	calibration hit gain

double	h_twodCor [kMaxHits] = {}
	calibration 2D correction

double	h_onedCor [kMaxHits] = {}
	calibration 1D cleanup correction

double	h_WeightPionHypo [kMaxHits] = {}
	weight for pion hypothesis from KalmanFitterInfo

double	h_WeightKaonHypo [kMaxHits] = {}
	weight for kaon hypothesis from KalmanFitterInfo

double	h_WeightProtonHypo [kMaxHits] = {}
	weight for proton hypothesis from KalmanFitterInfo

int	h_foundByTrackFinder [kMaxHits] = {}
	the 'found by track finder' flag for the given hit

DBObjPtr< CDCDedxScaleFactor >	m_DBScaleFactor
	Scale factor to make electrons ~1.

DBObjPtr< CDCDedxWireGain >	m_DBWireGains
	Wire gain DB object.

DBObjPtr< CDCDedxRunGain >	m_DBRunGain
	Run gain DB object.

DBObjPtr< CDCDedxInjectionTime >	m_DBInjectTime
	time gain/reso DB object

DBObjPtr< CDCDedxCosineCor >	m_DBCosineCor
	Electron saturation correction DB object.

DBObjPtr< CDCDedx2DCell >	m_DB2DCell
	2D correction DB object

DBObjPtr< CDCDedx1DCell >	m_DB1DCell
	1D correction DB object

DBObjPtr< CDCDedxADCNonLinearity >	m_DBNonlADC
	hadron saturation non linearity

DBObjPtr< CDCDedxCosineEdge >	m_DBCosEdgeCor
	cosine edge calibration

DBObjPtr< CDCDedxHadronCor >	m_DBHadronCor
	hadron saturation parameters

std::vector< double >	m_hadronpars
	hadron saturation parameters

bool	m_isDeadwire
	write only active wires

bool	m_isHitLevel
	Flag to switch on/off hit level info.

bool	m_isExtraVar
	Flag to switch on/off extra level info and some available w/ release/5 only.

bool	m_isRelative
	Flag to switch on/off relative constants.

bool	m_isCorrection
	Flag to switch on/off corrections.

std::string	m_name
	The name of the module, saved as a string (user-modifiable)

std::string	m_type
	The type of the module, saved as a string.

std::string	m_package
	Package this module is found in (may be empty).

std::string	m_description
	The description of the module.

unsigned int	m_propertyFlags
	The properties of the module as bitwise or (with \|) of EModulePropFlags.

LogConfig	m_logConfig
	The log system configuration of the module.

ModuleParamList	m_moduleParamList
	List storing and managing all parameter of the module.

bool	m_hasReturnValue
	True, if the return value is set.

int	m_returnValue
	The return value.

std::vector< ModuleCondition >	m_conditions
	Module condition, only non-null if set.

Static Private Attributes
static const int	kMaxHits = 200
	default hit level index

Detailed Description

Extracts dE/dx information for calibration testing.

Writes a ROOT file.

Definition at line 52 of file HitLevelInfoWriter.h.

Member Typedef Documentation

◆ EAfterConditionPath

typedef ModuleCondition::EAfterConditionPath EAfterConditionPath

inherited

Forward the EAfterConditionPath definition from the ModuleCondition.

Definition at line 88 of file Module.h.

Member Enumeration Documentation

◆ EModulePropFlags

enum EModulePropFlags

inherited

Each module can be tagged with property flags, which indicate certain features of the module.

Enumerator
c_Input	This module is an input module (reads data).
c_Output	This module is an output module (writes data).
c_ParallelProcessingCertified	This module can be run in parallel processing mode safely (All I/O must be done through the data store, in particular, the module must not write any files.)
c_HistogramManager	This module is used to manage histograms accumulated by other modules.
c_InternalSerializer	This module is an internal serializer/deserializer for parallel processing.
c_TerminateInAllProcesses	When using parallel processing, call this module's terminate() function in all processes(). This will also ensure that there is exactly one process (single-core if no parallel modules found) or at least one input, one main and one output process.
c_DontCollectStatistics	No statistics is collected for this module.

Definition at line 77 of file Module.h.

                          {
      c_Input                       = 1,  
      c_Output                      = 2,  
      c_ParallelProcessingCertified = 4,  
      c_HistogramManager            = 8, 
      c_InternalSerializer          = 16,  
      c_TerminateInAllProcesses     = 32,  
      c_DontCollectStatistics       = 64,  
    };

Constructor & Destructor Documentation

◆ HitLevelInfoWriterModule()

HitLevelInfoWriterModule ( )

Default constructor.

Definition at line 25 of file HitLevelInfoWriter.cc.

                                                   : Module()
{
 
  setDescription("Extract dE/dx information for calibration development.");
 
  addParam("outputBaseName", m_strOutputBaseName, "Suffix for output file name", string("HLInfo.root"));
  addParam("particleLists", m_strParticleList, "Vector of ParticleLists to save", vector<string>());
  addParam("enableHitLevel", m_isHitLevel, "True or False for Hit level variables", false);
  addParam("m_isExtraVar", m_isExtraVar, "True or False for extra track/hit level variables", false);
  addParam("nodeadwire", m_isDeadwire, "True or False for deadwire hit variables", false);
  addParam("relativeCorrections", m_isRelative, "If true, apply corrections relative to those used in production", false);
  addParam("corrections", m_isCorrection, "If true, apply corrections", true);
 
}

◆ ~HitLevelInfoWriterModule()

~HitLevelInfoWriterModule ( )

virtual

Destructor.

Definition at line 40 of file HitLevelInfoWriter.cc.

40{ }

Member Function Documentation

◆ beginRun()

virtual void beginRun ( void )

inlinevirtualinherited

Called when entering a new run.

Called at the beginning of each run, the method gives you the chance to change run dependent constants like alignment parameters, etc.

This method can be implemented by subclasses.

Definition at line 146 of file Module.h.

146{};

◆ bookOutput()

void bookOutput ( std::string filename )

Create the output TFiles and TTrees.

Definition at line 727 of file HitLevelInfoWriter.cc.

{
  // register output root file
  m_file.push_back(new TFile(filename.c_str(), "RECREATE"));
  m_tree.push_back(new TTree("track", "dE/dx information"));
 
  int i = m_tree.size() - 1;
  m_tree[i]->SetDirectory(0);
 
  // event level information (from emd)
  m_tree[i]->Branch("exp", &m_expID, "exp/I");
  m_tree[i]->Branch("run", &m_runID, "run/I");
  m_tree[i]->Branch("event", &m_eventID, "event/I");
 
  // track level information (from tfr)
  m_tree[i]->Branch("dz", &m_dz, "dz/D");
  m_tree[i]->Branch("dr", &m_dr, "dr/D");
  m_tree[i]->Branch("chi2", &m_chi2, "chi2/D");
  m_tree[i]->Branch("injtime", &m_injtime, "injtime/D");
  m_tree[i]->Branch("isher", &m_injring, "isher/D");
 
  if (m_isExtraVar) {
    m_tree[i]->Branch("d0", &m_d0, "d0/D");
    m_tree[i]->Branch("z0", &m_z0, "z0/D");
    m_tree[i]->Branch("dphi", &m_dphi, "dphi/D");
    m_tree[i]->Branch("vx0", &m_vx0, "vx0/D");
    m_tree[i]->Branch("vy0", &m_vy0, "vy0/D");
    m_tree[i]->Branch("vz0", &m_vz0, "vz0/D");
    m_tree[i]->Branch("tanlambda", &m_tanlambda, "tanlambda/D");
    m_tree[i]->Branch("phi0", &m_phi0, "phi0/D");
    m_tree[i]->Branch("e1_9", &m_e1_9, "e1_9/D");
    m_tree[i]->Branch("e9_21", &m_e9_21, "e9_21/D");
    m_tree[i]->Branch("eclsnhits", &m_eclsnHits, "eclsnhits/D");
    m_tree[i]->Branch("klmLayers", &m_klmLayers, "klmLayers/I");
  }
 
  // track level information (from cdt)
  m_tree[i]->Branch("nCDChits", &m_nCDChits, "nCDChits/D");
  m_tree[i]->Branch("inCDC", &m_inCDC, "inCDC/I");
  m_tree[i]->Branch("track", &m_trackID, "track/I");
  m_tree[i]->Branch("length", &m_length, "length/D");
  m_tree[i]->Branch("charge", &m_charge, "charge/I");
  m_tree[i]->Branch("costh", &m_cosTheta, "costh/D");
  m_tree[i]->Branch("pF", &m_pCDC, "pF/D");
  m_tree[i]->Branch("p", &m_p, "p/D");
  m_tree[i]->Branch("pt", &m_pt, "pt/D");
  m_tree[i]->Branch("ioasym", &m_ioasym, "ioasym/D");
  m_tree[i]->Branch("phi", &m_phi, "phi/D");
  m_tree[i]->Branch("pdg", &m_PDG, "pdg/D");
  // track level dE/dx measurements
  m_tree[i]->Branch("mean", &m_mean, "mean/D");
  m_tree[i]->Branch("dedx", &m_trunc, "dedx/D");
  m_tree[i]->Branch("dedxnosat", &m_truncNoSat, "dedxnosat/D");
  m_tree[i]->Branch("dedxerr", &m_error, "dedxerr/D");
 
  // other track level information
  m_tree[i]->Branch("eop", &m_eop, "eop/D");
  m_tree[i]->Branch("e", &m_e, "e/D");
 
  // calibration constants
  m_tree[i]->Branch("scale", &m_scale, "scale/D");
  m_tree[i]->Branch("coscor", &m_cosCor, "coscor/D");
  m_tree[i]->Branch("cosedgecor", &m_cosEdgeCor, "cosedgecor/D");
  m_tree[i]->Branch("rungain", &m_runGain, "rungain/D");
  m_tree[i]->Branch("timegain", &m_timeGain, "timegain/D");
  m_tree[i]->Branch("timereso", &m_timeReso, "timereso/D");
 
  // PID values
  m_tree[i]->Branch("chiE", &m_chie, "chiE/D");
  m_tree[i]->Branch("chiMu", &m_chimu, "chiMu/D");
  m_tree[i]->Branch("chiPi", &m_chipi, "chiPi/D");
  m_tree[i]->Branch("chiK", &m_chik, "chiK/D");
  m_tree[i]->Branch("chiP", &m_chip, "chiP/D");
  m_tree[i]->Branch("chiD", &m_chid, "chiD/D");
 
  m_tree[i]->Branch("chiEOld", &m_chieOld, "chiEOld/D");
  m_tree[i]->Branch("chiMuOld", &m_chimuOld, "chiMuOld/D");
  m_tree[i]->Branch("chiPiOld", &m_chipiOld, "chiPiOld/D");
  m_tree[i]->Branch("chiKOld", &m_chikOld, "chiKOld/D");
  m_tree[i]->Branch("chiPOld", &m_chipOld, "chiPOld/D");
  m_tree[i]->Branch("chiDOld", &m_chidOld, "chiDOld/D");
 
  m_tree[i]->Branch("pmeanE", &m_pmeane, "pmeanE/D");
  m_tree[i]->Branch("pmeanMu", &m_pmeanmu, "pmeanMu/D");
  m_tree[i]->Branch("pmeanPi", &m_pmeanpi, "pmeanPi/D");
  m_tree[i]->Branch("pmeanK", &m_pmeank, "pmeanK/D");
  m_tree[i]->Branch("pmeanP", &m_pmeanp, "pmeanP/D");
  m_tree[i]->Branch("pmeanD", &m_pmeand, "pmeanD/D");
 
  m_tree[i]->Branch("presE", &m_prese, "presE/D");
  m_tree[i]->Branch("presMu", &m_presmu, "presMu/D");
  m_tree[i]->Branch("presPi", &m_prespi, "presPi/D");
  m_tree[i]->Branch("presK", &m_presk, "presK/D");
  m_tree[i]->Branch("presP", &m_presp, "presP/D");
  m_tree[i]->Branch("presD", &m_presd, "presD/D");
 
  // layer level information
  m_tree[i]->Branch("lNHits", &l_nhits, "lNHits/I");
  m_tree[i]->Branch("lNHitsUsed", &l_nhitsused, "lNHitsUsed/I");
  m_tree[i]->Branch("lNHitsCombined", l_nhitscombined, "lNHitsCombined[lNHits]/I");
  m_tree[i]->Branch("lWireLongestHit", l_wirelongesthit, "lWireLongestHit[lNHits]/I");
  m_tree[i]->Branch("lLayer", l_layer, "lLayer[lNHits]/I");
  m_tree[i]->Branch("lPath", l_path, "lPath[lNHits]/D");
  m_tree[i]->Branch("lDedx", l_dedx, "lDedx[lNHits]/D");
 
  // hit level information
  if (m_isHitLevel) {
    m_tree[i]->Branch("hNHits", &h_nhits, "hNHits/I");
    m_tree[i]->Branch("hLWire", h_lwire, "hLWire[hNHits]/I");
    m_tree[i]->Branch("hWire", h_wire, "hWire[hNHits]/I");
    m_tree[i]->Branch("hLayer", h_layer, "hLayer[hNHits]/I");
    m_tree[i]->Branch("hPath", h_path, "hPath[hNHits]/D");
    m_tree[i]->Branch("hDedx", h_dedx, "hDedx[hNHits]/D");
    m_tree[i]->Branch("hADCRaw", h_adcraw, "hADCRaw[hNHits]/D");
    m_tree[i]->Branch("hADCCorr", h_adccorr, "hADCCorr[hNHits]/D");
    m_tree[i]->Branch("hDoca", h_doca, "hDoca[hNHits]/D");
    m_tree[i]->Branch("hNDoca", h_ndoca, "hNDoca[hNHits]/D");
    m_tree[i]->Branch("hNDocaRS", h_ndocaRS, "hNDocaRS[hNHits]/D");
    m_tree[i]->Branch("hEnta", h_enta, "hEnta[hNHits]/D");
    m_tree[i]->Branch("hEntaRS", h_entaRS, "hEntaRS[hNHits]/D");
    m_tree[i]->Branch("hDriftT", h_driftT, "hDriftT[hNHits]/D");
    m_tree[i]->Branch("hDriftD", h_driftD, "hDriftD[hNHits]/D");
    m_tree[i]->Branch("hFacnlADC", h_facnladc, "hFacnlADC[hNHits]/D");
    m_tree[i]->Branch("hWireGain", h_wireGain, "hWireGain[hNHits]/D");
    m_tree[i]->Branch("hTwodcor", h_twodCor, "hTwodcor[hNHits]/D");
    m_tree[i]->Branch("hOnedcor", h_onedCor, "hOnedcor[hNHits]/D");
 
    if (m_isExtraVar) {
      m_tree[i]->Branch("hWeightPionHypo", h_WeightPionHypo, "hWeightPionHypo[hNHits]/D");
      m_tree[i]->Branch("hWeightKaonHypo", h_WeightKaonHypo, "hWeightKaonHypo[hNHits]/D");
      m_tree[i]->Branch("hWeightProtonHypo", h_WeightProtonHypo, "hWeightProtonHypo[hNHits]/D");
      m_tree[i]->Branch("hFoundByTrackFinder", h_foundByTrackFinder, "hFoundByTrackFinder[hNHits]/I");
    }
  }
}

◆ calculateMeans()

void calculateMeans	(	double *	mean,
		double *	truncatedMean,
		double *	truncatedMeanErr,
		const std::vector< double > &	dedx ) const

Save arithmetic and truncated mean for the 'dedx' values.

Parameters

mean	calculated arithmetic mean
truncatedMean	calculated truncated mean
truncatedMeanErr	error for truncatedMean
dedx	input values

Definition at line 630 of file HitLevelInfoWriter.cc.

{
  // Calculate the truncated average by skipping the lowest & highest
  // events in the array of dE/dx values
  vector<double> sortedDedx = dedx;
  sort(sortedDedx.begin(), sortedDedx.end());
  sortedDedx.erase(remove(sortedDedx.begin(), sortedDedx.end(), 0), sortedDedx.end());
  sortedDedx.shrink_to_fit();
 
  double truncMeanTmp = 0.0, meanTmp = 0.0, sumOfSqs = 0.0;
  int nValTrunc = 0;
  const int numDedx = sortedDedx.size();
 
  // add a factor of 0.51 here to make sure we are rounding appropriately...
  const int lEdgeTrunc = int(numDedx * 0.05 + 0.51);
  const int hEdgeTrunc = int(numDedx * (1 - 0.25) + 0.51);
  for (int i = 0; i < numDedx; i++) {
    meanTmp += sortedDedx[i];
    if (i >= lEdgeTrunc and i < hEdgeTrunc) {
      truncMeanTmp += sortedDedx[i];
      sumOfSqs += sortedDedx[i] * sortedDedx[i];
      nValTrunc++;
    }
  }
 
  if (numDedx != 0) meanTmp /= numDedx;
 
  if (nValTrunc != 0) truncMeanTmp /= nValTrunc;
  else truncMeanTmp = meanTmp;
 
  *mean = meanTmp;
  *truncMean = truncMeanTmp;
 
  if (nValTrunc > 1)
    *truncMeanErr = sqrt(sumOfSqs / double(nValTrunc) - truncMeanTmp * truncMeanTmp) / double(nValTrunc - 1);
  else *truncMeanErr = 0;
 
}

◆ clearEntries()

void clearEntries ( )

private

Clear the arrays before filling an event.

Definition at line 688 of file HitLevelInfoWriter.cc.

{
  for (int il = 0; il < 200; ++il) {
    l_nhitscombined[il] = 0;
    l_wirelongesthit[il] = 0;
    l_layer[il] = 0;
    l_path[il] = 0;
    l_dedx[il] = 0;
  }
 
  if (m_isHitLevel) {
    for (int ihit = 0; ihit < 200; ++ihit) {
      h_lwire[ihit] = 0;
      h_wire[ihit] = 0;
      h_layer[ihit] = 0;
      h_path[ihit] = 0;
      h_dedx[ihit] = 0;
      h_adcraw[ihit] = 0;
      h_adccorr[ihit] = 0;
      h_doca[ihit] = 0;
      h_ndoca[ihit] = 0;
      h_ndocaRS[ihit] = 0;
      h_enta[ihit] = 0;
      h_entaRS[ihit] = 0;
      h_driftT[ihit] = 0;
      // h_driftD[ihit] = 0;
      h_facnladc[ihit] = 0;
      h_wireGain[ihit] = 0;
      h_twodCor[ihit] = 0;
      h_onedCor[ihit] = 0;
      h_WeightPionHypo[ihit] = 0;
      h_WeightKaonHypo[ihit] = 0;
      h_WeightProtonHypo[ihit] = 0;
      h_foundByTrackFinder[ihit] = 0;
    }
  }
}

◆ clone()

std::shared_ptr< PathElement > clone ( ) const

overridevirtualinherited

Create an independent copy of this module.

Note that parameters are shared, so changing them on a cloned module will also affect the original module.

Implements PathElement.

Definition at line 179 of file Module.cc.

{
  ModulePtr newModule = ModuleManager::Instance().registerModule(getType());
  newModule->m_moduleParamList.setParameters(getParamList());
  newModule->setName(getName());
  newModule->m_package = m_package;
  newModule->m_propertyFlags = m_propertyFlags;
  newModule->m_logConfig = m_logConfig;
  newModule->m_conditions = m_conditions;
 
  return newModule;
}

◆ D2I()

double D2I	(	const double	cosTheta,
		const double	D ) const

hadron saturation parameterization part 2

Definition at line 568 of file HitLevelInfoWriter.cc.

{
  double absCosTheta   = fabs(cosTheta);
  double projection    = pow(absCosTheta, m_hadronpars[3]) + m_hadronpars[2];
  if (projection == 0) {
    B2WARNING("Something wrong with dE/dx hadron constants!");
    return D;
  }
 
  double chargeDensity = D / projection;
  double numerator     = 1 + m_hadronpars[0] * chargeDensity;
  double denominator   = 1 + m_hadronpars[1] * chargeDensity;
 
  if (denominator == 0) {
    B2WARNING("Something wrong with dE/dx hadron constants!");
    return D;
  }
 
  double I = D * m_hadronpars[4] * numerator / denominator;
  return I;
}

◆ def_beginRun()

virtual void def_beginRun ( )

inlineprotectedvirtualinherited

Wrapper method for the virtual function beginRun() that has the implementation to be used in a call from Python.

Reimplemented in PyModule.

Definition at line 425 of file Module.h.

425{ beginRun(); }

◆ def_endRun()

virtual void def_endRun ( )

inlineprotectedvirtualinherited

This method can receive that the current run ends as a call from the Python side.

For regular C++-Modules that forwards the call to the regular endRun() method.

Reimplemented in PyModule.

Definition at line 438 of file Module.h.

438{ endRun(); }

◆ def_event()

virtual void def_event ( )

inlineprotectedvirtualinherited

Wrapper method for the virtual function event() that has the implementation to be used in a call from Python.

Reimplemented in PyModule.

Definition at line 431 of file Module.h.

431{ event(); }

◆ def_initialize()

virtual void def_initialize ( )

inlineprotectedvirtualinherited

Wrappers to make the methods without "def_" prefix callable from Python.

Overridden in PyModule. Wrapper method for the virtual function initialize() that has the implementation to be used in a call from Python.

Reimplemented in PyModule.

Definition at line 419 of file Module.h.

419{ initialize(); }

◆ def_terminate()

virtual void def_terminate ( )

inlineprotectedvirtualinherited

Wrapper method for the virtual function terminate() that has the implementation to be used in a call from Python.

Reimplemented in PyModule.

Definition at line 444 of file Module.h.

444{ terminate(); }

◆ endRun()

virtual void endRun ( void )

inlinevirtualinherited

This method is called if the current run ends.

Use this method to store information, which should be aggregated over one run.

This method can be implemented by subclasses.

Definition at line 165 of file Module.h.

165{};

◆ evalCondition()

bool evalCondition ( ) const

inherited

If at least one condition was set, it is evaluated and true returned if at least one condition returns true.

If no condition or result value was defined, the method returns false. Otherwise, the condition is evaluated and true returned, if at least one condition returns true. To speed up the evaluation, the condition strings were already parsed in the method if_value().

Returns: True if at least one condition and return value exists and at least one condition expression was evaluated to true.

Definition at line 96 of file Module.cc.

{
  if (m_conditions.empty()) return false;
 
  //okay, a condition was set for this Module...
  if (!m_hasReturnValue) {
    B2FATAL("A condition was set for '" << getName() << "', but the module did not set a return value!");
  }
 
  for (const auto& condition : m_conditions) {
    if (condition.evaluate(m_returnValue)) {
      return true;
    }
  }
  return false;
}

◆ event()

void event ( void )

overridevirtual

This method is called for each event.

All processing of the event takes place in this method.

Reimplemented from Module.

Definition at line 69 of file HitLevelInfoWriter.cc.

{
 
  int nParticleList = m_strParticleList.size();
 
  // **************************************************
  //  LOOP OVER dE/dx measurements for all tracks if
  //    no particle list is specified.
  // **************************************************
  if (nParticleList == 0) {
    for (int idedx = 0; idedx < m_dedxTracks.getEntries(); idedx++) {
 
      CDCDedxTrack* dedxTrack = m_dedxTracks[idedx];
      if (!dedxTrack) {
        B2WARNING("No dedx related track...");
        continue;
      }
 
      const Track* track = dedxTrack->getRelatedFrom<Track>();
      if (!track) {
        B2WARNING("No related track...");
        continue;
      }
 
      const TrackFitResult* fitResult = track->getTrackFitResultWithClosestMass(Const::pion);
      if (!fitResult) {
        B2WARNING("No related fit for this track...");
        continue;
      }
 
      if (dedxTrack->size() == 0 || dedxTrack->size() > 200) continue;
      if (dedxTrack->getCosTheta() < TMath::Cos(150.0 * TMath::DegToRad()))continue; //-0.866
      if (dedxTrack->getCosTheta() > TMath::Cos(17.0 * TMath::DegToRad())) continue; //0.95
 
      // fill the event meta data
      StoreObjPtr<EventMetaData> evtMetaData;
      m_expID = evtMetaData->getExperiment();
      m_runID = evtMetaData->getRun();
      m_eventID = evtMetaData->getEvent();
 
      m_injring = dedxTrack->getInjectionRing();
      m_injtime = dedxTrack->getInjectionTime();
 
      //--------REMOVEBAL--------
      //when CDST are reproduced with injection time
      if (m_injtime == -1 || m_injring == -1) {
        if (m_TTDInfo.isValid() && m_TTDInfo->hasInjection()) {
          m_injring = m_TTDInfo->isHER();
          m_injtime =  m_TTDInfo->getTimeSinceLastInjectionInMicroSeconds();
        }
      }
      //--------REMOVEBAL--------
 
      // fill the E/P
      const ECLCluster* eclCluster = track->getRelated<ECLCluster>();
      if (eclCluster and eclCluster->hasHypothesis(ECLCluster::EHypothesisBit::c_nPhotons)) {
        m_eop = (eclCluster->getEnergy(ECLCluster::EHypothesisBit::c_nPhotons)) / (fitResult->getMomentum().R());
        m_e = eclCluster->getEnergy(ECLCluster::EHypothesisBit::c_nPhotons);
        if (m_isExtraVar) {
          m_e1_9 = eclCluster->getE1oE9();
          m_e9_21 = eclCluster->getE9oE21();
          m_eclsnHits = eclCluster->getNumberOfCrystals();
        }
        // fill the muon depth
        const KLMCluster* klmCluster = eclCluster->getRelated<KLMCluster>();
        if (klmCluster) m_klmLayers = klmCluster->getLayers();
      }
 
      // fill the TTree with the Track information
      fillTrack(fitResult);
 
      // fill the TTree with the CDCDedxTrack information
      fillDedx(dedxTrack);
 
      // fill the TTree
      m_tree[0]->Fill();
    }
  }
 
  // **************************************************
  //
  //  LOOP OVER particles in the given particle lists
  //
  // **************************************************
  for (int iList = 0; iList < nParticleList; iList++) {
    // make sure the list exists and is not empty
    StoreObjPtr<ParticleList> particlelist(m_strParticleList[iList]);
    if (!particlelist or particlelist->getListSize(true) == 0) {
      //B2WARNING("ParticleList " << m_strParticleList[iList] << " not found or empty, skipping");
      continue;
    }
 
    // loop over the particles in the list and follow the links to the
    // dE/dx information (Particle -> PIDLikelihood -> Track -> CDCDedxTrack)
    for (unsigned int iPart = 0; iPart < particlelist->getListSize(true); iPart++) {
      Particle* part = particlelist->getParticle(iPart, true);
      if (!part) {
        B2WARNING("No particles...");
        continue;
      }
      PIDLikelihood* pid = part->getRelatedTo<PIDLikelihood>();
      if (!pid) {
        B2WARNING("No related PID likelihood...");
        continue;
      }
      Track* track = pid->getRelatedFrom<Track>();
      if (!track) {
        B2WARNING("No related track...");
        continue;
      }
      CDCDedxTrack* dedxTrack = track->getRelatedTo<CDCDedxTrack>();
      if (!dedxTrack) {
        B2WARNING("No related CDCDedxTrack...");
        continue;
      }
      string ptype = m_strParticleList[iList].substr(0, m_strParticleList[iList].find(":"));
      const TrackFitResult* fitResult = track->getTrackFitResult(Const::pion);
      if (ptype != "pi+") {
        if (ptype == "K+") fitResult = track->getTrackFitResultWithClosestMass(Const::kaon);
        else if (ptype == "p+") fitResult = track->getTrackFitResultWithClosestMass(Const::proton);
        else if (ptype == "deuteron") fitResult = track->getTrackFitResultWithClosestMass(Const::deuteron);
        else if (ptype == "mu+") fitResult = track->getTrackFitResultWithClosestMass(Const::muon);
        else if (ptype == "e+") fitResult = track->getTrackFitResultWithClosestMass(Const::electron);
      }
      if (!fitResult) {
        B2WARNING("No related fit for this track...");
        continue;
      }
 
      if (dedxTrack->size() == 0 || dedxTrack->size() > 200) continue;
      //if out CDC (dont add as we dont correct via correctionmodules)
      if (dedxTrack->getCosTheta() < TMath::Cos(150.0 * TMath::DegToRad()))continue; //-0.866
      if (dedxTrack->getCosTheta() > TMath::Cos(17.0 * TMath::DegToRad())) continue; //0.95
 
      // fill the event meta data
      StoreObjPtr<EventMetaData> evtMetaData;
      m_expID = evtMetaData->getExperiment();
      m_runID = evtMetaData->getRun();
      m_eventID = evtMetaData->getEvent();
 
      m_injring = dedxTrack->getInjectionRing();
      m_injtime = dedxTrack->getInjectionTime();
 
      //--------REMOVEBAL--------
      //when CDST are reproduced with injection time
      if (m_injtime == -1 || m_injring == -1) {
        if (m_TTDInfo.isValid() && m_TTDInfo->hasInjection()) {
          m_injring = m_TTDInfo->isHER();
          m_injtime =  m_TTDInfo->getTimeSinceLastInjectionInMicroSeconds();
        }
      }
      //--------REMOVEBAL--------
 
      // fill the E/P
      const ECLCluster* eclCluster = track->getRelated<ECLCluster>();
      if (eclCluster and eclCluster->hasHypothesis(ECLCluster::EHypothesisBit::c_nPhotons)) {
        m_eop = (eclCluster->getEnergy(ECLCluster::EHypothesisBit::c_nPhotons)) / (fitResult->getMomentum().R());
        m_e = eclCluster->getEnergy(ECLCluster::EHypothesisBit::c_nPhotons);
        if (m_isExtraVar) {
          m_e1_9 = eclCluster->getE1oE9();
          m_e9_21 = eclCluster->getE9oE21();
          m_eclsnHits = eclCluster->getNumberOfCrystals();
        }
        // fill the muon depth
        const KLMCluster* klmCluster = eclCluster->getRelated<KLMCluster>();
        if (klmCluster) m_klmLayers = klmCluster->getLayers();
      }
 
      // fill the TTree with the Track information
      fillTrack(fitResult);
 
      // fill the TTree with the CDCDedxTrack information
      fillDedx(dedxTrack);
 
      // fill the TTree
      m_tree[iList]->Fill();
    }
  }
}

◆ exposePythonAPI()

void exposePythonAPI ( )

staticinherited

Exposes methods of the Module class to Python.

Definition at line 325 of file Module.cc.

{
  // to avoid confusion between std::arg and boost::python::arg we want a shorthand namespace as well
  namespace bp = boost::python;
 
  docstring_options options(true, true, false); //userdef, py sigs, c++ sigs
 
  void (Module::*setReturnValueInt)(int) = &Module::setReturnValue;
 
  enum_<Module::EAfterConditionPath>("AfterConditionPath",
                                     R"(Determines execution behaviour after a conditional path has been executed:
 
.. attribute:: END
 
  End processing of this path after the conditional path. (this is the default for if_value() etc.)
 
.. attribute:: CONTINUE
 
  After the conditional path, resume execution after this module.)")
  .value("END", Module::EAfterConditionPath::c_End)
  .value("CONTINUE", Module::EAfterConditionPath::c_Continue)
  ;
 
  /* Do not change the names of >, <, ... we use them to serialize conditional paths */
  enum_<Belle2::ModuleCondition::EConditionOperators>("ConditionOperator")
  .value(">", Belle2::ModuleCondition::EConditionOperators::c_GT)
  .value("<", Belle2::ModuleCondition::EConditionOperators::c_ST)
  .value(">=", Belle2::ModuleCondition::EConditionOperators::c_GE)
  .value("<=", Belle2::ModuleCondition::EConditionOperators::c_SE)
  .value("==", Belle2::ModuleCondition::EConditionOperators::c_EQ)
  .value("!=", Belle2::ModuleCondition::EConditionOperators::c_NE)
  ;
 
  enum_<Module::EModulePropFlags>("ModulePropFlags",
                                  R"(Flags to indicate certain low-level features of modules, see :func:`Module.set_property_flags()`, :func:`Module.has_properties()`. Most useful flags are:
 
.. attribute:: PARALLELPROCESSINGCERTIFIED
 
    This module can be run in parallel processing mode safely (All I/O must be done through the data store, in particular, the module must not write any files.)
 
.. attribute:: HISTOGRAMMANAGER
 
  This module is used to manage histograms accumulated by other modules
 
.. attribute:: TERMINATEINALLPROCESSES
 
  When using parallel processing, call this module's terminate() function in all processes. This will also ensure that there is exactly one process (single-core if no parallel modules found) or at least one input, one main and one output process.
)")
  .value("INPUT", Module::EModulePropFlags::c_Input)
  .value("OUTPUT", Module::EModulePropFlags::c_Output)
  .value("PARALLELPROCESSINGCERTIFIED", Module::EModulePropFlags::c_ParallelProcessingCertified)
  .value("HISTOGRAMMANAGER", Module::EModulePropFlags::c_HistogramManager)
  .value("INTERNALSERIALIZER", Module::EModulePropFlags::c_InternalSerializer)
  .value("TERMINATEINALLPROCESSES", Module::EModulePropFlags::c_TerminateInAllProcesses)
  ;
 
  //Python class definition
  class_<Module, PyModule> module("Module", R"(
Base class for Modules.
 
A module is the smallest building block of the framework.
A typical event processing chain consists of a Path containing
modules. By inheriting from this base class, various types of
modules can be created. To use a module, please refer to
:func:`Path.add_module()`.  A list of modules is available by running
``basf2 -m`` or ``basf2 -m package``, detailed information on parameters is
given by e.g. ``basf2 -m RootInput``.
 
The 'Module Development' section in the manual provides detailed information
on how to create modules, setting parameters, or using return values/conditions:
https://xwiki.desy.de/xwiki/rest/p/f4fa4/#HModuleDevelopment
 
)");
  module
  .def("__str__", &Module::getPathString)
  .def("name", &Module::getName, return_value_policy<copy_const_reference>(),
       "Returns the name of the module. Can be changed via :func:`set_name() <Module.set_name()>`, use :func:`type() <Module.type()>` for identifying a particular module class.")
  .def("type", &Module::getType, return_value_policy<copy_const_reference>(),
       "Returns the type of the module (i.e. class name minus 'Module')")
  .def("set_name", &Module::setName, args("name"), R"(
Set custom name, e.g. to distinguish multiple modules of the same type.
 
>>> path.add_module('EventInfoSetter')
>>> ro = path.add_module('RootOutput', branchNames=['EventMetaData'])
>>> ro.set_name('RootOutput_metadata_only')
>>> print(path)
[EventInfoSetter -> RootOutput_metadata_only]
 
)")
  .def("description", &Module::getDescription, return_value_policy<copy_const_reference>(),
       "Returns the description of this module.")
  .def("package", &Module::getPackage, return_value_policy<copy_const_reference>(),
       "Returns the package this module belongs to.")
  .def("available_params", &_getParamInfoListPython,
       "Return list of all module parameters as `ModuleParamInfo` instances")
  .def("has_properties", &Module::hasProperties, (bp::arg("properties")),
       R"DOCSTRING(Allows to check if the module has the given properties out of `ModulePropFlags` set.
 
>>> if module.has_properties(ModulePropFlags.PARALLELPROCESSINGCERTIFIED):
>>>     ...
 
Parameters:
  properties (int): bitmask of `ModulePropFlags` to check for.
)DOCSTRING")
  .def("set_property_flags", &Module::setPropertyFlags, args("property_mask"),
       "Set module properties in the form of an OR combination of `ModulePropFlags`.");
  {
    // python signature is too crowded, make ourselves
    docstring_options subOptions(true, false, false); //userdef, py sigs, c++ sigs
    module
    .def("if_value", &Module::if_value,
         (bp::arg("expression"), bp::arg("condition_path"), bp::arg("after_condition_path")= Module::EAfterConditionPath::c_End),
         R"DOCSTRING(if_value(expression, condition_path, after_condition_path=AfterConditionPath.END)
 
Sets a conditional sub path which will be executed after this
module if the return value set in the module passes the given ``expression``.
 
Modules can define a return value (int or bool) using ``setReturnValue()``,
which can be used in the steering file to split the Path based on this value, for example
 
>>> module_with_condition.if_value("<1", another_path)
 
In case the return value of the ``module_with_condition`` for a given event is
less than 1, the execution will be diverted into ``another_path`` for this event.
 
You could for example set a special return value if an error occurs, and divert
the execution into a path containing :b2:mod:`RootOutput` if it is found;
saving only the data producing/produced by the error.
 
After a conditional path has executed, basf2 will by default stop processing
the path for this event. This behaviour can be changed by setting the
``after_condition_path`` argument.
 
Parameters:
  expression (str): Expression to determine if the conditional path should be executed.
      This should be one of the comparison operators ``<``, ``>``, ``<=``,
      ``>=``, ``==``, or ``!=`` followed by a numerical value for the return value
  condition_path (Path): path to execute in case the expression is fulfilled
  after_condition_path (AfterConditionPath): What to do once the ``condition_path`` has been executed.
)DOCSTRING")
    .def("if_false", &Module::if_false,
         (bp::arg("condition_path"), bp::arg("after_condition_path")= Module::EAfterConditionPath::c_End),
         R"DOC(if_false(condition_path, after_condition_path=AfterConditionPath.END)
 
Sets a conditional sub path which will be executed after this module if
the return value of the module evaluates to False. This is equivalent to
calling `if_value` with ``expression=\"<1\"``)DOC")
    .def("if_true", &Module::if_true,
         (bp::arg("condition_path"), bp::arg("after_condition_path")= Module::EAfterConditionPath::c_End),
         R"DOC(if_true(condition_path, after_condition_path=AfterConditionPath.END)
 
Sets a conditional sub path which will be executed after this module if
the return value of the module evaluates to True. It is equivalent to
calling `if_value` with ``expression=\">=1\"``)DOC");
  }
  module
  .def("has_condition", &Module::hasCondition,
       "Return true if a conditional path has been set for this module "
       "using `if_value`, `if_true` or `if_false`")
  .def("get_all_condition_paths", &_getAllConditionPathsPython,
       "Return a list of all conditional paths set for this module using "
       "`if_value`, `if_true` or `if_false`")
  .def("get_all_conditions", &_getAllConditionsPython,
       "Return a list of all conditional path expressions set for this module using "
       "`if_value`, `if_true` or `if_false`")
  .add_property("logging", make_function(&Module::getLogConfig, return_value_policy<reference_existing_object>()),
                &Module::setLogConfig)

◆ fillDedx()

void fillDedx ( CDCDedxTrack * dedxTrack )

private

Fill the TTree with the information from a CDCDedxTrack object.

Definition at line 303 of file HitLevelInfoWriter.cc.

{
  // clear the containers first
  clearEntries();
 
  m_trackID = dedxTrack->trackID();
  m_length = dedxTrack->getLength();
  m_charge = dedxTrack->getCharge();
  m_cosTheta = dedxTrack->getCosTheta();
  m_PDG = dedxTrack->getPDG();
 
  m_pCDC = dedxTrack->getMomentum();
  if (m_charge < 0) {
    m_pCDC *= -1;
  }
 
  h_nhits = dedxTrack->size();
 
  // Get the calibration constants
  m_scale = m_DBScaleFactor->getScaleFactor();
  m_runGain = m_DBRunGain->getRunGain();
  m_cosCor = m_DBCosineCor->getMean(m_cosTheta);
  m_timeGain = m_DBInjectTime->getCorrection("mean", m_injring, m_injtime);
  m_timeReso = m_DBInjectTime->getCorrection("reso", m_injring, m_injtime);
 
  if (m_cosTheta <= -0.850 || m_cosTheta >= 0.950) {
    m_cosEdgeCor = m_DBCosEdgeCor->getMean(m_cosTheta);
  } else {
    m_cosEdgeCor = 1.0;
  }
  m_hadronpars = m_DBHadronCor->getHadronPars();
 
  //variable to save layer variables
  map<int, vector<double>>l_var;
  double cdcChi[Const::ChargedStable::c_SetSize];
 
  //Modify the hit level dedx
  if (m_isCorrection) {
    recalculateDedx(dedxTrack, l_var, cdcChi);
    m_chie = cdcChi[0];
    m_chimu = cdcChi[1];
    m_chipi = cdcChi[2];
    m_chik = cdcChi[3];
    m_chip = cdcChi[4];
    m_chid = cdcChi[5];
  } else {
    m_mean = dedxTrack->getDedxMean();
    m_trunc = dedxTrack->getDedx();
    m_truncNoSat = dedxTrack->getDedxNoSat();
    m_error = dedxTrack->getDedxError();
    m_chie = dedxTrack->getChi(0);
    m_chimu = dedxTrack->getChi(1);
    m_chipi = dedxTrack->getChi(2);
    m_chik = dedxTrack->getChi(3);
    m_chip = dedxTrack->getChi(4);
    m_chid = dedxTrack->getChi(5);
  }
 
  m_pmeane = dedxTrack->getPmean(0);
  m_pmeanmu = dedxTrack->getPmean(1);
  m_pmeanpi = dedxTrack->getPmean(2);
  m_pmeank = dedxTrack->getPmean(3);
  m_pmeanp = dedxTrack->getPmean(4);
  m_pmeand = dedxTrack->getPmean(5);
 
  m_prese = dedxTrack->getPreso(0);
  m_presmu = dedxTrack->getPreso(1);
  m_prespi = dedxTrack->getPreso(2);
  m_presk = dedxTrack->getPreso(3);
  m_presp = dedxTrack->getPreso(4);
  m_presd = dedxTrack->getPreso(5);
 
  // Get the vector of dE/dx values for all layers
  double lout = 0, lin = 0, increment = 0;
  int lastlayer = 0;
 
  if (m_isCorrection) {
    l_nhits = l_var[0].size();
    const int lEdgeTrunc = int(l_nhits * 0.05 + 0.51);
    const int hEdgeTrunc = int(l_nhits * (1 - 0.25) + 0.51);
    l_nhitsused = hEdgeTrunc - lEdgeTrunc;
  } else {
    l_nhits = dedxTrack->getNLayerHits();
    l_nhitsused = dedxTrack->getNLayerHitsUsed();
  }
 
  for (int il = 0; il < l_nhits; ++il) {
    if (m_isCorrection) {
      l_nhitscombined[il] = l_var[0][il];
      l_wirelongesthit[il] = l_var[1][il];
      l_layer[il] = l_var[2][il];
      l_path[il] = l_var[3][il];
      l_dedx[il] = l_var[4][il];
    } else {
      l_nhitscombined[il] = dedxTrack->getNHitsCombined(il);
      l_wirelongesthit[il] = dedxTrack->getWireLongestHit(il);
      l_layer[il] = dedxTrack->getLayer(il);
      l_path[il] = dedxTrack->getLayerPath(il);
      l_dedx[il] = dedxTrack->getLayerDedx(il);
    }
    if (l_layer[il] > lastlayer) lout++;
    else if (l_layer[il] < lastlayer) lin++;
    else continue;
 
    lastlayer = l_layer[il];
    increment++;
  }
  m_ioasym = (lout - lin) / increment;
 
  // Get the vector of dE/dx values for all hits
  if (m_isHitLevel) {
    for (int ihit = 0; ihit < h_nhits; ++ihit) {
 
      if (m_isDeadwire) continue;
 
      h_lwire[ihit] = dedxTrack->getWireInLayer(ihit);
      h_wire[ihit] = dedxTrack->getWire(ihit);
      h_layer[ihit] = dedxTrack->getHitLayer(ihit);
      h_path[ihit] = dedxTrack->getPath(ihit);
      h_dedx[ihit] = dedxTrack->getDedx(ihit);
      h_adcraw[ihit] = dedxTrack->getADCBaseCount(ihit);
      h_adccorr[ihit] = dedxTrack->getADCCount(ihit);
      h_doca[ihit] = dedxTrack->getDoca(ihit);
      h_ndoca[ihit] = h_doca[ihit] / dedxTrack->getCellHalfWidth(ihit);
      h_ndocaRS[ihit] = dedxTrack->getDocaRS(ihit) / dedxTrack->getCellHalfWidth(ihit);
      h_enta[ihit] = dedxTrack->getEnta(ihit);
      h_entaRS[ihit] = dedxTrack->getEntaRS(ihit);
      h_driftT[ihit] = dedxTrack->getDriftT(ihit);
      h_driftD[ihit] = dedxTrack->getDriftD(ihit);
 
      // Get extra variable from tracking
      if (m_isExtraVar) {
        h_WeightPionHypo[ihit] = dedxTrack->getWeightPionHypo(ihit);
        h_WeightKaonHypo[ihit] = dedxTrack->getWeightKaonHypo(ihit);
        h_WeightProtonHypo[ihit] = dedxTrack->getWeightProtonHypo(ihit);
        h_foundByTrackFinder[ihit] = dedxTrack->getFoundByTrackFinder(ihit);
      }
      // Get the calibration constants
      h_facnladc[ihit] = dedxTrack->getNonLADCCorrection(ihit);
      h_wireGain[ihit] = m_DBWireGains->getWireGain(h_wire[ihit]);
      h_twodCor[ihit] = m_DB2DCell->getMean(h_layer[ihit], h_ndocaRS[ihit], h_entaRS[ihit]);
      h_onedCor[ihit] = m_DB1DCell->getMean(h_layer[ihit], h_entaRS[ihit]);
    }
  }
 
}

◆ fillTrack()

void fillTrack ( const TrackFitResult * fitResult )

private

Fill the TTree with the information from the track fit.

Definition at line 264 of file HitLevelInfoWriter.cc.

{
  ROOT::Math::XYZVector trackMom = fitResult->getMomentum();
  m_p = trackMom.R();
  m_pt = trackMom.Rho();
  m_phi = trackMom.Phi();
 
  m_theta = trackMom.Theta() * 180. / TMath::Pi(); //in degree
  if (m_theta > 17. && m_theta < 150.)m_inCDC =  1;
  else m_inCDC = 0;
 
  if (fitResult->getChargeSign() < 0) {
    m_p *= -1;
    m_pt *= -1;
  }
 
  ROOT::Math::XYZVector trackPos = fitResult->getPosition();
  m_vx0 = trackPos.X();
  m_vy0 = trackPos.Y();
  m_vz0 = trackPos.Z();
 
  m_d0 = fitResult->getD0();
  m_z0 = fitResult->getZ0();
  m_chi2 = fitResult->getPValue();
  m_tanlambda = fitResult->getTanLambda();
  m_phi0 = fitResult->getPhi0();
  m_nCDChits =  fitResult->getHitPatternCDC().getNHits();
 
  static DBObjPtr<BeamSpot> beamSpotDB;
  const auto& frame = ReferenceFrame::GetCurrent();
  UncertainHelix helix = fitResult->getUncertainHelix();
  helix.passiveMoveBy(ROOT::Math::XYZVector(beamSpotDB->getIPPosition()));
  m_dr = frame.getVertex(ROOT::Math::XYZVector(helix.getPerigee())).Rho();
  m_dphi = frame.getVertex(ROOT::Math::XYZVector(helix.getPerigee())).Phi();
  m_dz = frame.getVertex(ROOT::Math::XYZVector(helix.getPerigee())).Z();
 
}

◆ getAfterConditionPath()

Module::EAfterConditionPath getAfterConditionPath ( ) const

inherited

What to do after the conditional path is finished.

(defaults to c_End if no condition is set)

Definition at line 133 of file Module.cc.

{
  if (m_conditions.empty()) return EAfterConditionPath::c_End;
 
  //okay, a condition was set for this Module...
  if (!m_hasReturnValue) {
    B2FATAL("A condition was set for '" << getName() << "', but the module did not set a return value!");
  }
 
  for (const auto& condition : m_conditions) {
    if (condition.evaluate(m_returnValue)) {
      return condition.getAfterConditionPath();
    }
  }
 
  return EAfterConditionPath::c_End;
}

◆ getAllConditionPaths()

std::vector< std::shared_ptr< Path > > getAllConditionPaths ( ) const

inherited

Return all condition paths currently set (no matter if the condition is true or not).

Definition at line 150 of file Module.cc.

{
  std::vector<std::shared_ptr<Path>> allConditionPaths;
  for (const auto& condition : m_conditions) {
    allConditionPaths.push_back(condition.getPath());
  }
 
  return allConditionPaths;
}

◆ getAllConditions()

const std::vector< ModuleCondition > & getAllConditions ( ) const

inlineinherited

Return all set conditions for this module.

Definition at line 323 of file Module.h.

    {
      return m_conditions;
    }

◆ getCondition()

const ModuleCondition * getCondition ( ) const

inlineinherited

Return a pointer to the first condition (or nullptr, if none was set)

Definition at line 313 of file Module.h.

    {
      if (m_conditions.empty()) {
        return nullptr;
      } else {
        return &m_conditions.front();
      }
    }

◆ getConditionPath()

std::shared_ptr< Path > getConditionPath ( ) const

inherited

Returns the path of the last true condition (if there is at least one, else reaturn a null pointer).

Definition at line 113 of file Module.cc.

{
  PathPtr p;
  if (m_conditions.empty()) return p;
 
  //okay, a condition was set for this Module...
  if (!m_hasReturnValue) {
    B2FATAL("A condition was set for '" << getName() << "', but the module did not set a return value!");
  }
 
  for (const auto& condition : m_conditions) {
    if (condition.evaluate(m_returnValue)) {
      return condition.getPath();
    }
  }
 
  // if none of the conditions were true, return a null pointer.
  return p;
}

◆ GetCorrection()

double GetCorrection	(	int &	adc,
		int	layer,
		int	wireID,
		double	doca,
		double	enta,
		double	costheta,
		double	ring,
		double	time ) const

Function to get the correction factor.

Definition at line 542 of file HitLevelInfoWriter.cc.

{
  double correction = 1.0;
  correction *= m_DBScaleFactor->getScaleFactor();
  correction *= m_DBRunGain->getRunGain();
  correction *= m_DB2DCell->getMean(layer, doca, enta);
  correction *= m_DB1DCell->getMean(layer, enta);
  correction *= m_DBCosineCor->getMean(costheta);
  correction *= m_DBInjectTime->getCorrection("mean", ring, time);
  if (costheta <= -0.850 || costheta >= 0.950) correction *= m_DBCosEdgeCor->getMean(costheta);
  if (m_DBWireGains->getWireGain(wireID) > 0)  correction *= m_DBWireGains->getWireGain(wireID);
 
  //last is only for abs constant
  if (!m_isRelative) adc = m_DBNonlADC->getCorrectedADC(adc, layer);
 
  return correction;
}

◆ getDescription()

const std::string & getDescription ( ) const

inlineinherited

Returns the description of the module.

Definition at line 201 of file Module.h.

201{return m_description;}

◆ getFileNames()

virtual std::vector< std::string > getFileNames ( bool outputFiles )

inlinevirtualinherited

Return a list of output filenames for this modules.

This will be called when basf2 is run with "--dry-run" if the module has set either the c_Input or c_Output properties.

If the parameter outputFiles is false (for modules with c_Input) the list of input filenames should be returned (if any). If outputFiles is true (for modules with c_Output) the list of output files should be returned (if any).

If a module has sat both properties this member is called twice, once for each property.

The module should return the actual list of requested input or produced output filenames (including handling of input/output overrides) so that the grid system can handle input/output files correctly.

This function should return the same value when called multiple times. This is especially important when taking the input/output overrides from Environment as they get consumed when obtained so the finalized list of output files should be stored for subsequent calls.

Reimplemented in RootInputModule, RootOutputModule, and StorageRootOutputModule.

Definition at line 133 of file Module.h.

    {
      return std::vector<std::string>();
    }

◆ getLogConfig()

LogConfig & getLogConfig ( )

inlineinherited

Returns the log system configuration.

Definition at line 224 of file Module.h.

224{return m_logConfig;}

◆ getModules()

std::list< ModulePtr > getModules ( ) const

inlineoverrideprivatevirtualinherited

no submodules, return empty list

Implements PathElement.

Definition at line 505 of file Module.h.

505{ return std::list<ModulePtr>(); }

◆ getName()

const std::string & getName ( ) const

inlineinherited

Returns the name of the module.

This can be changed via e.g. set_name() in the steering file to give more useful names if there is more than one module of the same type.

For identifying the type of a module, using getType() (or type() in Python) is recommended.

Definition at line 186 of file Module.h.

186{return m_name;}

◆ getPackage()

const std::string & getPackage ( ) const

inlineinherited

Returns the package this module is in.

Definition at line 196 of file Module.h.

196{return m_package;}

◆ getParamInfoListPython()

std::shared_ptr< boost::python::list > getParamInfoListPython ( ) const

inherited

Returns a python list of all parameters.

Each item in the list consists of the name of the parameter, a string describing its type, a python list of all default values and the description of the parameter.

Returns: A python list containing the parameters of this parameter list.

Definition at line 279 of file Module.cc.

{
  return m_moduleParamList.getParamInfoListPython();
}

◆ getParamList()

const ModuleParamList & getParamList ( ) const

inlineinherited

Return module param list.

Definition at line 362 of file Module.h.

362{ return m_moduleParamList; }

◆ getPathString()

std::string getPathString ( ) const

overrideprivatevirtualinherited

return the module name.

Implements PathElement.

Definition at line 192 of file Module.cc.

{
 
  std::string output = getName();
 
  for (const auto& condition : m_conditions) {
    output += condition.getString();
  }
 
  return output;
}

◆ getReturnValue()

int getReturnValue ( ) const

inlineinherited

Return the return value set by this module.

This value is only meaningful if hasReturnValue() is true

Definition at line 380 of file Module.h.

380{ return m_returnValue; }

◆ getType()

const std::string & getType ( ) const

inherited

Returns the type of the module (i.e.

class name minus 'Module')

Definition at line 41 of file Module.cc.

{
  if (m_type.empty())
    B2FATAL("Module type not set for " << getName());
  return m_type;
}

◆ HadronCorrection()

void HadronCorrection	(	double	costheta,
		double &	dedx ) const

Function to apply the hadron correction.

Definition at line 562 of file HitLevelInfoWriter.cc.

{
  dedx = D2I(costheta, I2D(costheta, 1.00) / 1.00 * dedx);
}

◆ hasCondition()

bool hasCondition ( ) const

inlineinherited

Returns true if at least one condition was set for the module.

Definition at line 310 of file Module.h.

310{ return not m_conditions.empty(); };

◆ hasProperties()

bool hasProperties ( unsigned int propertyFlags ) const

inherited

Returns true if all specified property flags are available in this module.

Parameters

propertyFlags Ored EModulePropFlags which should be compared with the module flags.

Definition at line 160 of file Module.cc.

{
  return (propertyFlags & m_propertyFlags) == propertyFlags;
}

◆ hasReturnValue()

bool hasReturnValue ( ) const

inlineinherited

Return true if this module has a valid return value set.

Definition at line 377 of file Module.h.

377{ return m_hasReturnValue; }

◆ hasUnsetForcedParams()

bool hasUnsetForcedParams ( ) const

inherited

Returns true and prints error message if the module has unset parameters which the user has to set in the steering file.

Definition at line 166 of file Module.cc.

{
  auto missing = m_moduleParamList.getUnsetForcedParams();
  std::string allMissing = "";
  for (const auto& s : missing)
    allMissing += s + " ";
  if (!missing.empty())
    B2ERROR("The following required parameters of Module '" << getName() << "' were not specified: " << allMissing <<
            "\nPlease add them to your steering file.");
  return !missing.empty();
}

◆ I2D()

double I2D	(	const double	cosTheta,
		const double	I ) const

hadron saturation parameterization part 1

Definition at line 591 of file HitLevelInfoWriter.cc.

{
  double absCosTheta = fabs(cosTheta);
  double projection  = pow(absCosTheta, m_hadronpars[3]) + m_hadronpars[2];
 
  if (projection == 0 || m_hadronpars[4] == 0) {
    B2WARNING("Something wrong with dE/dx hadron constants!");
    return I;
  }
 
  double a =  m_hadronpars[0] / projection;
  double b =  1 - m_hadronpars[1] / projection * (I / m_hadronpars[4]);
  double c = -1.0 * I / m_hadronpars[4];
 
  if (b == 0 && a == 0) {
    B2WARNING("both a and b coefficiants for hadron correction are 0");
    return I;
  }
 
  double discr = b * b - 4.0 * a * c;
  if (discr < 0) {
    B2WARNING("negative discriminant; return uncorrectecd value");
    return I;
  }
 
  double D = (a != 0) ? (-b + sqrt(discr)) / (2.0 * a) : -c / b;
  if (D < 0) {
    B2WARNING("D is less 0! will try another solution");
    D = (a != 0) ? (-b - sqrt(discr)) / (2.0 * a) : -c / b;
    if (D < 0) {
      B2WARNING("D is still less 0! just return uncorrectecd value");
      return I;
    }
  }
 
  return D;
}

◆ if_false()

void if_false	(	const std::shared_ptr< Path > &	path,
		EAfterConditionPath	afterConditionPath = EAfterConditionPath::c_End )

inherited

A simplified version to add a condition to the module.

Please note that successive calls of this function will add more than one condition to the module. If more than one condition results in true, only the last of them will be used.

Please be careful: Avoid creating cyclic paths, e.g. by linking a condition to a path which is processed before the path where this module is located in.

It is equivalent to the if_value() method, using the expression "<1". This method is meant to be used together with the setReturnValue(bool value) method.

Parameters

path	Shared pointer to the Path which will be executed if the return value is false.
afterConditionPath	What to do after executing 'path'.

Definition at line 85 of file Module.cc.

{
  if_value("<1", path, afterConditionPath);
}

◆ if_true()

void if_true	(	const std::shared_ptr< Path > &	path,
		EAfterConditionPath	afterConditionPath = EAfterConditionPath::c_End )

inherited

A simplified version to set the condition of the module.

Please note that successive calls of this function will add more than one condition to the module. If more than one condition results in true, only the last of them will be used.

Please be careful: Avoid creating cyclic paths, e.g. by linking a condition to a path which is processed before the path where this module is located in.

It is equivalent to the if_value() method, using the expression ">=1". This method is meant to be used together with the setReturnValue(bool value) method.

Parameters

path	Shared pointer to the Path which will be executed if the return value is true.
afterConditionPath	What to do after executing 'path'.

Definition at line 90 of file Module.cc.

{
  if_value(">=1", path, afterConditionPath);
}

◆ if_value()

void if_value	(	const std::string &	expression,
		const std::shared_ptr< Path > &	path,
		EAfterConditionPath	afterConditionPath = EAfterConditionPath::c_End )

inherited

Add a condition to the module.

Please note that successive calls of this function will add more than one condition to the module. If more than one condition results in true, only the last of them will be used.

See https://xwiki.desy.de/xwiki/rest/p/a94f2 or ModuleCondition for a description of the syntax.

Please be careful: Avoid creating cyclic paths, e.g. by linking a condition to a path which is processed before the path where this module is located in.

Parameters

expression	The expression of the condition.
path	Shared pointer to the Path which will be executed if the condition is evaluated to true.
afterConditionPath	What to do after executing 'path'.

Definition at line 79 of file Module.cc.

{
  m_conditions.emplace_back(expression, path, afterConditionPath);
}

◆ initialize()

void initialize ( void )

overridevirtual

Initialize the module.

Reimplemented from Module.

Definition at line 42 of file HitLevelInfoWriter.cc.

{
 
  B2INFO("Creating a ROOT file for the hit level information...");
 
  // required inputs
  m_dedxTracks.isRequired();
  m_tracks.isRequired();
  m_trackFitResults.isRequired();
  m_eclClusters.isOptional();
  m_klmClusters.isOptional();
 
  // build a map to relate input strings to the right particle type
  map<string, string> pdgMap = {{"pi+", "Const::pion.getPDGCode()"}, {"K+", "Const::kaon.getPDGCode()"}, {"mu+", "Const::muon.getPDGCode()"}, {"e+", "Const::electron.getPDGCode()"}, {"p+", "Const::proton.getPDGCode()"}, {"deuteron", "Const::deuteron.getPDGCode()"}};
 
  // if no particle lists are given, write out all tracks
  if (m_strParticleList.size() == 0) bookOutput(m_strOutputBaseName);
 
  // create a new output file for each particle list specified
  for (unsigned int i = 0; i < m_strParticleList.size(); i++) {
    // strip the name of the particle lists to make this work
    string pdg = pdgMap[m_strParticleList[i].substr(0, m_strParticleList[i].find(":"))];
    string filename = string(m_strOutputBaseName + "_PID" + pdg + ".root");
    bookOutput(filename);
  }
}

◆ recalculateDedx()

void recalculateDedx	(	CDCDedxTrack *	dedxTrack,
		std::map< int, std::vector< double > > &	l_var,
		double(&)	cdcChi[Const::ChargedStable::c_SetSize] )

Function to recalculate the dedx with latest constants.

Definition at line 451 of file HitLevelInfoWriter.cc.

{
  vector<double> newLayerHits;
  double newLayerDe = 0, newLayerDx = 0;
  int nhitscombined = 0; // number of hits combined per layer
  int wirelongesthit = 0; // wire number of longest hit
  double longesthit = 0; // path length of longest hit
 
  for (int ihit = 0; ihit < h_nhits; ++ihit) {
    int jadcbase = dedxTrack->getADCBaseCount(ihit);
    int jLayer = dedxTrack->getHitLayer(ihit);
    double jWire = dedxTrack->getWire(ihit);
    double jNDocaRS = dedxTrack->getDocaRS(ihit) / dedxTrack->getCellHalfWidth(ihit);
    double jEntaRS = dedxTrack->getEntaRS(ihit);
    double jPath = dedxTrack->getPath(ihit);
 
    double correction = dedxTrack->getScaleFactor() * dedxTrack->getRunGain() * dedxTrack->getTimeMean() *
                        dedxTrack->getCosineCorrection() * dedxTrack->getCosEdgeCorrection() *
                        dedxTrack->getTwoDCorrection(ihit) * dedxTrack->getOneDCorrection(ihit) *
                        dedxTrack->getNonLADCCorrection(ihit);
    if (dedxTrack->getWireGain(ihit) > 0) correction *= dedxTrack->getWireGain(ihit); //also keep dead wire
 
    if (m_isRelative) {
      //get same base adc + rel correction factor
      correction *= GetCorrection(jadcbase, jLayer, jWire, jNDocaRS, jEntaRS, m_cosTheta, m_injring, m_injtime);
      if (!m_DBWireGains && dedxTrack->getWireGain(ihit) == 0) correction = 0;
    } else {
      //get modified adc + abs correction factor
      correction = GetCorrection(jadcbase, jLayer, jWire, jNDocaRS, jEntaRS, m_cosTheta, m_injring, m_injtime);
    }
 
    double newhitdedx = 1.0;
 
    newhitdedx *= jadcbase * sqrt(1 - m_cosTheta * m_cosTheta) / jPath;
 
    if (correction != 0) {
      newhitdedx /= correction;
      if (m_DBWireGains->getWireGain(jWire) != 0) {
        newLayerDe += jadcbase / correction;
        newLayerDx += jPath;
        if (jPath > longesthit) {
          longesthit = jPath;
          wirelongesthit = jWire;
        }
        nhitscombined++;
      }
    } else newhitdedx = 0;
 
    dedxTrack->setDedx(ihit, newhitdedx);
 
    if (ihit + 1 < h_nhits && dedxTrack->getHitLayer(ihit + 1) == jLayer) {
      continue;
    } else {
      if (newLayerDx != 0) {
        double totalDistance = newLayerDx / sqrt(1 - m_cosTheta * m_cosTheta);
        double newLayerDedx = newLayerDe / totalDistance ;
        newLayerHits.push_back(newLayerDedx);
        l_var[0].push_back(nhitscombined);
        l_var[1].push_back(wirelongesthit);
        l_var[2].push_back(jLayer);
        l_var[3].push_back(totalDistance);
        l_var[4].push_back(newLayerDedx);
      }
 
      newLayerDe = 0;
      newLayerDx = 0;
      nhitscombined = 0;
      wirelongesthit = 0;
      longesthit = 0;
    }
  }
 
  // recalculate the truncated means
  double dedxmean, dedxtrunc, dedxtruncNoSat, dedxerror;
 
  calculateMeans(&dedxmean, &dedxtruncNoSat, &dedxerror, newLayerHits);
 
  dedxtrunc = dedxtruncNoSat;
 
  HadronCorrection(m_cosTheta, dedxtrunc);
  m_mean = dedxmean;
  m_trunc = dedxtrunc;
  m_truncNoSat = dedxtruncNoSat;
  m_error = dedxerror;
 
  // save the PID information
  saveChiValue(cdcChi, dedxTrack, m_trunc);
}

◆ saveChiValue()

void saveChiValue	(	double(&)	chi[Const::ChargedStable::c_SetSize],
		CDCDedxTrack *	dedxTrack,
		double	dedx ) const

for all particles, save chi values into 'chi' chi array of chi values to be modified

Definition at line 671 of file HitLevelInfoWriter.cc.

{
  // determine a chi value for each particle type
  Const::ParticleSet set = Const::chargedStableSet;
  for (const Const::ChargedStable pdgIter : set) {
 
    // determine the predicted mean and resolution
    double mean = dedxTrack->getPmean(pdgIter.getIndex());
    double sigma = dedxTrack->getPreso(pdgIter.getIndex());
 
    // fill the chi value for this particle type
    if (sigma != 0) chi[pdgIter.getIndex()] = ((dedx - mean) / (sigma));
  }
}

◆ setAbortLevel()

void setAbortLevel ( int abortLevel )

inherited

Configure the abort log level.

Definition at line 67 of file Module.cc.

{
  m_logConfig.setAbortLevel(static_cast<LogConfig::ELogLevel>(abortLevel));
}

◆ setDebugLevel()

void setDebugLevel ( int debugLevel )

inherited

Configure the debug messaging level.

Definition at line 61 of file Module.cc.

{
  m_logConfig.setDebugLevel(debugLevel);
}

◆ setDescription()

void setDescription ( const std::string & description )

protectedinherited

Sets the description of the module.

Parameters

description A description of the module.

Definition at line 214 of file Module.cc.

{
  m_description = description;
}

◆ setLogConfig()

void setLogConfig ( const LogConfig & logConfig )

inlineinherited

Set the log system configuration.

Definition at line 229 of file Module.h.

229{m_logConfig = logConfig;}

◆ setLogInfo()

void setLogInfo	(	int	logLevel,
		unsigned int	logInfo )

inherited

Configure the printed log information for the given level.

Parameters

logLevel	The log level (one of LogConfig::ELogLevel)
logInfo	What kind of info should be printed? ORed combination of LogConfig::ELogInfo flags.

Definition at line 73 of file Module.cc.

{
  m_logConfig.setLogInfo(static_cast<LogConfig::ELogLevel>(logLevel), logInfo);
}

◆ setLogLevel()

void setLogLevel ( int logLevel )

inherited

Configure the log level.

Definition at line 55 of file Module.cc.

{
  m_logConfig.setLogLevel(static_cast<LogConfig::ELogLevel>(logLevel));
}

◆ setName()

void setName ( const std::string & name )

inlineinherited

Set the name of the module.

Note: The module name is set when using the REG_MODULE macro, but the module can be renamed before calling process() using the set_name() function in your steering file.

Parameters

name	The name of the module

Definition at line 213 of file Module.h.

213{ m_name = name; };

◆ setParamList()

void setParamList ( const ModuleParamList & params )

inlineprotectedinherited

Replace existing parameter list.

Definition at line 500 of file Module.h.

500{ m_moduleParamList = params; }

◆ setParamPython()

void setParamPython	(	const std::string &	name,
		const boost::python::object &	pyObj )

privateinherited

Implements a method for setting boost::python objects.

The method supports the following types: list, dict, int, double, string, bool The conversion of the python object to the C++ type and the final storage of the parameter value is done in the ModuleParam class.

Parameters

name	The unique name of the parameter.
pyObj	The object which should be converted and stored as the parameter value.

Definition at line 234 of file Module.cc.

{
  LogSystem& logSystem = LogSystem::Instance();
  logSystem.updateModule(&(getLogConfig()), getName());
  try {
    m_moduleParamList.setParamPython(name, pyObj);
  } catch (std::runtime_error& e) {
    throw std::runtime_error("Cannot set parameter '" + name + "' for module '"
                             + m_name + "': " + e.what());
  }
 
  logSystem.updateModule(nullptr);
}

◆ setParamPythonDict()

void setParamPythonDict ( const boost::python::dict & dictionary )

privateinherited

Implements a method for reading the parameter values from a boost::python dictionary.

The key of the dictionary has to be the name of the parameter and the value has to be of one of the supported parameter types.

Parameters

dictionary The python dictionary from which the parameter values are read.

Definition at line 249 of file Module.cc.

{
 
  LogSystem& logSystem = LogSystem::Instance();
  logSystem.updateModule(&(getLogConfig()), getName());
 
  boost::python::list dictKeys = dictionary.keys();
  int nKey = boost::python::len(dictKeys);
 
  //Loop over all keys in the dictionary
  for (int iKey = 0; iKey < nKey; ++iKey) {
    boost::python::object currKey = dictKeys[iKey];
    boost::python::extract<std::string> keyProxy(currKey);
 
    if (keyProxy.check()) {
      const boost::python::object& currValue = dictionary[currKey];
      setParamPython(keyProxy, currValue);
    } else {
      B2ERROR("Setting the module parameters from a python dictionary: invalid key in dictionary!");
    }
  }
 
  logSystem.updateModule(nullptr);
}

◆ setPropertyFlags()

void setPropertyFlags ( unsigned int propertyFlags )

inherited

Sets the flags for the module properties.

Parameters

propertyFlags bitwise OR of EModulePropFlags

Definition at line 208 of file Module.cc.

{
  m_propertyFlags = propertyFlags;
}

◆ setReturnValue() [1/2]

void setReturnValue ( bool value )

protectedinherited

Sets the return value for this module as bool.

The bool value is saved as an integer with the convention 1 meaning true and 0 meaning false. The value can be used in the steering file to divide the analysis chain into several paths.

Parameters

value The value of the return value.

Definition at line 227 of file Module.cc.

{
  m_hasReturnValue = true;
  m_returnValue = value;
}

◆ setReturnValue() [2/2]

void setReturnValue ( int value )

protectedinherited

Sets the return value for this module as integer.

The value can be used in the steering file to divide the analysis chain into several paths.

Parameters

value The value of the return value.

Definition at line 220 of file Module.cc.

{
  m_hasReturnValue = true;
  m_returnValue = value;
}

◆ setType()

void setType ( const std::string & type )

protectedinherited

Set the module type.

Only for use by internal modules (which don't use the normal REG_MODULE mechanism).

Definition at line 48 of file Module.cc.

{
  if (!m_type.empty())
    B2FATAL("Trying to change module type from " << m_type << " is not allowed, the value is assumed to be fixed.");
  m_type = type;
}

◆ terminate()

void terminate ( void )

overridevirtual

End of the event processing.

Reimplemented from Module.

Definition at line 250 of file HitLevelInfoWriter.cc.

{
 
  for (unsigned int i = 0; i < m_file.size(); i++) {
    B2INFO("Done writing out the hit level information...\t" << m_tree[i]->GetEntries() << " tracks");
 
    // write the ttree to a root file
    m_file[i]->cd();
    m_tree[i]->Write();
    m_file[i]->Close();
  }
}

Member Data Documentation

◆ h_adccorr

double h_adccorr[kMaxHits] = {}

private

charge per hit corr by nonlinear ADC

Definition at line 249 of file HitLevelInfoWriter.h.

249{};

◆ h_adcraw

double h_adcraw[kMaxHits] = {}

private

charge per hit

Definition at line 248 of file HitLevelInfoWriter.h.

248{};

◆ h_dedx

double h_dedx[kMaxHits] = {}

private

charge per path length

Definition at line 247 of file HitLevelInfoWriter.h.

247{};

◆ h_doca

double h_doca[kMaxHits] = {}

private

distance of closest approach

Definition at line 250 of file HitLevelInfoWriter.h.

250{};

◆ h_driftD

double h_driftD[kMaxHits] = {}

private

drift distance

Definition at line 256 of file HitLevelInfoWriter.h.

256{};

◆ h_driftT

double h_driftT[kMaxHits] = {}

private

drift time

Definition at line 255 of file HitLevelInfoWriter.h.

255{};

◆ h_enta

double h_enta[kMaxHits] = {}

private

entrance angle

Definition at line 253 of file HitLevelInfoWriter.h.

253{};

◆ h_entaRS

double h_entaRS[kMaxHits] = {}

private

normalized + RS distance of entrance angle

Definition at line 254 of file HitLevelInfoWriter.h.

254{};

◆ h_facnladc

double h_facnladc[kMaxHits] = {}

private

calibration hit gain

Definition at line 257 of file HitLevelInfoWriter.h.

257{};

◆ h_foundByTrackFinder

int h_foundByTrackFinder[kMaxHits] = {}

private

the 'found by track finder' flag for the given hit

Definition at line 266 of file HitLevelInfoWriter.h.

266{};

◆ h_layer

int h_layer[kMaxHits] = {}

private

layer number

Definition at line 244 of file HitLevelInfoWriter.h.

244{};

◆ h_lwire

int h_lwire[kMaxHits] = {}

private

sense wire within layer

Definition at line 242 of file HitLevelInfoWriter.h.

242{};

◆ h_ndoca

double h_ndoca[kMaxHits] = {}

private

normalized distance of closest approach

Definition at line 251 of file HitLevelInfoWriter.h.

251{};

◆ h_ndocaRS

double h_ndocaRS[kMaxHits] = {}

private

normalized +RS distance of closest approach

Definition at line 252 of file HitLevelInfoWriter.h.

252{};

◆ h_nhits

int h_nhits { -1}

private

the number of good hits for this Track

Definition at line 241 of file HitLevelInfoWriter.h.

241{ -1};

◆ h_onedCor

double h_onedCor[kMaxHits] = {}

private

calibration 1D cleanup correction

Definition at line 260 of file HitLevelInfoWriter.h.

260{};

◆ h_path

double h_path[kMaxHits] = {}

private

path length in cell

Definition at line 246 of file HitLevelInfoWriter.h.

246{};

◆ h_twodCor

double h_twodCor[kMaxHits] = {}

private

calibration 2D correction

Definition at line 259 of file HitLevelInfoWriter.h.

259{};

◆ h_WeightKaonHypo

double h_WeightKaonHypo[kMaxHits] = {}

private

weight for kaon hypothesis from KalmanFitterInfo

Definition at line 264 of file HitLevelInfoWriter.h.

264{};

◆ h_WeightPionHypo

double h_WeightPionHypo[kMaxHits] = {}

private

weight for pion hypothesis from KalmanFitterInfo

Definition at line 263 of file HitLevelInfoWriter.h.

263{};

◆ h_WeightProtonHypo

double h_WeightProtonHypo[kMaxHits] = {}

private

weight for proton hypothesis from KalmanFitterInfo

Definition at line 265 of file HitLevelInfoWriter.h.

265{};

◆ h_wire

int h_wire[kMaxHits] = {}

private

sense wire ID

Definition at line 243 of file HitLevelInfoWriter.h.

243{};

◆ h_wireGain

double h_wireGain[kMaxHits] = {}

private

calibration hit gain

Definition at line 258 of file HitLevelInfoWriter.h.

258{};

◆ kMaxHits

const int kMaxHits = 200

staticprivate

default hit level index

Definition at line 230 of file HitLevelInfoWriter.h.

◆ l_dedx

double l_dedx[kMaxHits] = {}

private

dE/dx for this layer

Definition at line 238 of file HitLevelInfoWriter.h.

238{};

◆ l_layer

int l_layer[kMaxHits] = {}

private

layer number

Definition at line 236 of file HitLevelInfoWriter.h.

236{};

◆ l_nhits

int l_nhits { -1}

private

the total number of layer hits for this Track

Definition at line 232 of file HitLevelInfoWriter.h.

232{ -1};

◆ l_nhitscombined

int l_nhitscombined[kMaxHits] = {}

private

the number of hits combined this layer

Definition at line 234 of file HitLevelInfoWriter.h.

234{};

◆ l_nhitsused

int l_nhitsused { -1}

private

the total number of layer hits used for this Track

Definition at line 233 of file HitLevelInfoWriter.h.

233{ -1};

◆ l_path

double l_path[kMaxHits] = {}

private

distance travelled in this layer

Definition at line 237 of file HitLevelInfoWriter.h.

237{};

◆ l_wirelongesthit

int l_wirelongesthit[kMaxHits] = {}

private

the wire number of longest hit in this layer

Definition at line 235 of file HitLevelInfoWriter.h.

235{};

◆ m_charge

int m_charge {0}

private

the charge for this Track

Definition at line 164 of file HitLevelInfoWriter.h.

164{0};

◆ m_chi2

double m_chi2 { -1.}

private

chi^2 from track fit

Definition at line 158 of file HitLevelInfoWriter.h.

158{ -1.};

◆ m_chid

double m_chid { -1.}

private

modified chi value for deuteron hypothesis

Definition at line 214 of file HitLevelInfoWriter.h.

214{ -1.};

◆ m_chidOld

double m_chidOld { -1.}

private

chi value for deuteron hypothesis

Definition at line 207 of file HitLevelInfoWriter.h.

207{ -1.};

◆ m_chie

double m_chie { -1.}

private

modified chi value for electron hypothesis

Definition at line 209 of file HitLevelInfoWriter.h.

209{ -1.};

◆ m_chieOld

double m_chieOld { -1.}

private

chi value for electron hypothesis

Definition at line 202 of file HitLevelInfoWriter.h.

202{ -1.};

◆ m_chik

double m_chik { -1.}

private

modified chi value for kaon hypothesis

Definition at line 212 of file HitLevelInfoWriter.h.

212{ -1.};

◆ m_chikOld

double m_chikOld { -1.}

private

chi value for kaon hypothesis

Definition at line 205 of file HitLevelInfoWriter.h.

205{ -1.};

◆ m_chimu

double m_chimu { -1.}

private

modified chi value for muon hypothesis

Definition at line 210 of file HitLevelInfoWriter.h.

210{ -1.};

◆ m_chimuOld

double m_chimuOld { -1.}

private

chi value for muon hypothesis

Definition at line 203 of file HitLevelInfoWriter.h.

203{ -1.};

◆ m_chip

double m_chip { -1.}

private

modified chi value for proton hypothesis

Definition at line 213 of file HitLevelInfoWriter.h.

213{ -1.};

◆ m_chipi

double m_chipi { -1.}

private

modified chi value for pion hypothesis

Definition at line 211 of file HitLevelInfoWriter.h.

211{ -1.};

◆ m_chipiOld

double m_chipiOld { -1.}

private

chi value for pion hypothesis

Definition at line 204 of file HitLevelInfoWriter.h.

204{ -1.};

◆ m_chipOld

double m_chipOld { -1.}

private

chi value for proton hypothesis

Definition at line 206 of file HitLevelInfoWriter.h.

206{ -1.};

◆ m_conditions

std::vector<ModuleCondition> m_conditions

privateinherited

Module condition, only non-null if set.

Definition at line 520 of file Module.h.

◆ m_cosCor

double m_cosCor { -1.}

private

calibration cosine correction

Definition at line 195 of file HitLevelInfoWriter.h.

195{ -1.};

◆ m_cosEdgeCor

double m_cosEdgeCor { -1.}

private

calibration cosine edge correction

Definition at line 196 of file HitLevelInfoWriter.h.

196{ -1.};

◆ m_cosTheta

double m_cosTheta { -2.}

private

cos(theta) for the track

Definition at line 165 of file HitLevelInfoWriter.h.

165{ -2.};

◆ m_d0

double m_d0 {0.}

private

Signed distance to the POCA in the r-phi plane.

Definition at line 148 of file HitLevelInfoWriter.h.

148{0.};

◆ m_DB1DCell

DBObjPtr<CDCDedx1DCell> m_DB1DCell

private

1D correction DB object

Definition at line 275 of file HitLevelInfoWriter.h.

◆ m_DB2DCell

DBObjPtr<CDCDedx2DCell> m_DB2DCell

private

2D correction DB object

Definition at line 274 of file HitLevelInfoWriter.h.

◆ m_DBCosEdgeCor

DBObjPtr<CDCDedxCosineEdge> m_DBCosEdgeCor

private

cosine edge calibration

Definition at line 277 of file HitLevelInfoWriter.h.

◆ m_DBCosineCor

DBObjPtr<CDCDedxCosineCor> m_DBCosineCor

private

Electron saturation correction DB object.

Definition at line 273 of file HitLevelInfoWriter.h.

◆ m_DBHadronCor

DBObjPtr<CDCDedxHadronCor> m_DBHadronCor

private

hadron saturation parameters

Definition at line 278 of file HitLevelInfoWriter.h.

◆ m_DBInjectTime

DBObjPtr<CDCDedxInjectionTime> m_DBInjectTime

private

time gain/reso DB object

Definition at line 272 of file HitLevelInfoWriter.h.

◆ m_DBNonlADC

DBObjPtr<CDCDedxADCNonLinearity> m_DBNonlADC

private

hadron saturation non linearity

Definition at line 276 of file HitLevelInfoWriter.h.

◆ m_DBRunGain

DBObjPtr<CDCDedxRunGain> m_DBRunGain

private

Run gain DB object.

Definition at line 271 of file HitLevelInfoWriter.h.

◆ m_DBScaleFactor

DBObjPtr<CDCDedxScaleFactor> m_DBScaleFactor

private

Scale factor to make electrons ~1.

Definition at line 269 of file HitLevelInfoWriter.h.

◆ m_DBWireGains

DBObjPtr<CDCDedxWireGain> m_DBWireGains

private

Wire gain DB object.

Definition at line 270 of file HitLevelInfoWriter.h.

◆ m_dedxTracks

StoreArray<CDCDedxTrack> m_dedxTracks

private

Required array of CDCDedxTracks.

Definition at line 123 of file HitLevelInfoWriter.h.

◆ m_description

std::string m_description

privateinherited

The description of the module.

Definition at line 510 of file Module.h.

◆ m_dphi

double m_dphi { -1.}

private

POCA in degrees in respect to IP.

Definition at line 152 of file HitLevelInfoWriter.h.

152{ -1.};

◆ m_dr

double m_dr { -1.}

private

track d0 relative to IP

Definition at line 151 of file HitLevelInfoWriter.h.

151{ -1.};

◆ m_dz

double m_dz { -1.}

private

vertex or POCA in case of tracks z in respect to IPs

Definition at line 150 of file HitLevelInfoWriter.h.

150{ -1.};

◆ m_e

double m_e { -1.}

private

energy in the calorimeter

Definition at line 187 of file HitLevelInfoWriter.h.

187{ -1.};

◆ m_e1_9

double m_e1_9 { -1.}

private

ratio of energies of the central 1 crystal vs 3x3 crystals

Definition at line 188 of file HitLevelInfoWriter.h.

188{ -1.};

◆ m_e9_21

double m_e9_21 { -1.}

private

ratio of energies of the central 3x3 crystal vs 5x5 crystals

Definition at line 189 of file HitLevelInfoWriter.h.

189{ -1.};

◆ m_eclClusters

StoreArray<ECLCluster> m_eclClusters

private

Required array of input ECLClusters.

Definition at line 126 of file HitLevelInfoWriter.h.

◆ m_eclsnHits

double m_eclsnHits { -1.}

private

of clusters

Definition at line 191 of file HitLevelInfoWriter.h.

191{ -1.};

◆ m_eop

double m_eop { -1.}

private

energy over momentum in the calorimeter

Definition at line 186 of file HitLevelInfoWriter.h.

186{ -1.};

◆ m_error

double m_error { -1.}

private

standard deviation of the truncated mean

Definition at line 183 of file HitLevelInfoWriter.h.

183{ -1.};

◆ m_eventID

int m_eventID { -1}

private

event in which this Track was found

Definition at line 142 of file HitLevelInfoWriter.h.

142{ -1};

◆ m_expID

int m_expID { -1}

private

experiment in which this Track was found

Definition at line 140 of file HitLevelInfoWriter.h.

140{ -1};

◆ m_file

std::vector<TFile* > m_file

private

output ROOT files

Definition at line 120 of file HitLevelInfoWriter.h.

◆ m_filename

std::vector<std::string> m_filename

private

full names of the output ROOT files

Definition at line 119 of file HitLevelInfoWriter.h.

◆ m_hadronpars

std::vector<double> m_hadronpars

private

hadron saturation parameters

Definition at line 280 of file HitLevelInfoWriter.h.

◆ m_hasReturnValue

bool m_hasReturnValue

privateinherited

True, if the return value is set.

Definition at line 517 of file Module.h.

◆ m_inCDC

int m_inCDC { -1}

private

frack is CDC acceptance or not

Definition at line 161 of file HitLevelInfoWriter.h.

161{ -1};

◆ m_injring

double m_injring { -1.0}

private

HER injection status.

Definition at line 144 of file HitLevelInfoWriter.h.

144{ -1.0};

◆ m_injtime

double m_injtime { -1.0}

private

time since last injection in micro seconds

Definition at line 145 of file HitLevelInfoWriter.h.

145{ -1.0};

◆ m_ioasym

double m_ioasym { -1.}

private

asymmetry in increasing vs decreasing layer numbers per track

Definition at line 170 of file HitLevelInfoWriter.h.

170{ -1.};

◆ m_isCorrection

bool m_isCorrection

private

Flag to switch on/off corrections.

Definition at line 287 of file HitLevelInfoWriter.h.

◆ m_isDeadwire

bool m_isDeadwire

private

write only active wires

Definition at line 282 of file HitLevelInfoWriter.h.

◆ m_isExtraVar

bool m_isExtraVar

private

Flag to switch on/off extra level info and some available w/ release/5 only.

Definition at line 285 of file HitLevelInfoWriter.h.

◆ m_isHitLevel

bool m_isHitLevel

private

Flag to switch on/off hit level info.

Definition at line 284 of file HitLevelInfoWriter.h.

◆ m_isRelative

bool m_isRelative

private

Flag to switch on/off relative constants.

Definition at line 286 of file HitLevelInfoWriter.h.

◆ m_klmClusters

StoreArray<KLMCluster> m_klmClusters

private

Required array of input KLMClusters.

Definition at line 127 of file HitLevelInfoWriter.h.

◆ m_klmLayers

double m_klmLayers { -1.}

private

number of klm layers with hits

Definition at line 190 of file HitLevelInfoWriter.h.

190{ -1.};

◆ m_length

double m_length { -1.}

private

total path length of the Track

Definition at line 163 of file HitLevelInfoWriter.h.

163{ -1.};

◆ m_logConfig

LogConfig m_logConfig

privateinherited

The log system configuration of the module.

Definition at line 513 of file Module.h.

◆ m_mean

double m_mean { -1.}

private

dE/dx averaged

Definition at line 180 of file HitLevelInfoWriter.h.

180{ -1.};

◆ m_moduleParamList

ModuleParamList m_moduleParamList

privateinherited

List storing and managing all parameter of the module.

Definition at line 515 of file Module.h.

◆ m_name

std::string m_name

privateinherited

The name of the module, saved as a string (user-modifiable)

Definition at line 507 of file Module.h.

◆ m_nCDChits

double m_nCDChits { -1.}

private

Number of CDC hits associated to the track.

Definition at line 160 of file HitLevelInfoWriter.h.

160{ -1.};

◆ m_p

double m_p { -1.}

private

momentum from tracking

Definition at line 167 of file HitLevelInfoWriter.h.

167{ -1.};

◆ m_package

std::string m_package

privateinherited

Package this module is found in (may be empty).

Definition at line 509 of file Module.h.

◆ m_pCDC

double m_pCDC { -1.}

private

momentum valid in CDC

Definition at line 166 of file HitLevelInfoWriter.h.

166{ -1.};

◆ m_PDG

double m_PDG { -1.}

private

MC PID.

Definition at line 174 of file HitLevelInfoWriter.h.

174{ -1.};

◆ m_phi

double m_phi { -1.}

private

phi for the track

Definition at line 169 of file HitLevelInfoWriter.h.

169{ -1.};

◆ m_phi0

double m_phi0 { -1.}

private

Angle of the transverse momentum in the r-phi plane.

Definition at line 157 of file HitLevelInfoWriter.h.

157{ -1.};

◆ m_pmeand

double m_pmeand { -1.}

private

pred mean for deuteron hypothesis

Definition at line 228 of file HitLevelInfoWriter.h.

228{ -1.};

◆ m_pmeane

double m_pmeane { -1.}

private

pred mean for electron hypothesis

Definition at line 223 of file HitLevelInfoWriter.h.

223{ -1.};

◆ m_pmeank

double m_pmeank { -1.}

private

pred mean for kaon hypothesis

Definition at line 226 of file HitLevelInfoWriter.h.

226{ -1.};

◆ m_pmeanmu

double m_pmeanmu { -1.}

private

pred mean for muon hypothesis

Definition at line 224 of file HitLevelInfoWriter.h.

224{ -1.};

◆ m_pmeanp

double m_pmeanp { -1.}

private

pred mean for proton hypothesis

Definition at line 227 of file HitLevelInfoWriter.h.

227{ -1.};

◆ m_pmeanpi

double m_pmeanpi { -1.}

private

pred mean for pion hypothesis

Definition at line 225 of file HitLevelInfoWriter.h.

225{ -1.};

◆ m_presd

double m_presd { -1.}

private

pred reso for deuteron hypothesis

Definition at line 221 of file HitLevelInfoWriter.h.

221{ -1.};

◆ m_prese

double m_prese { -1.}

private

pred reso for electron hypothesis

Definition at line 216 of file HitLevelInfoWriter.h.

216{ -1.};

◆ m_presk

double m_presk { -1.}

private

pred reso for kaon hypothesis

Definition at line 219 of file HitLevelInfoWriter.h.

219{ -1.};

◆ m_presmu

double m_presmu { -1.}

private

pred reso for muon hypothesis

Definition at line 217 of file HitLevelInfoWriter.h.

217{ -1.};

◆ m_presp

double m_presp { -1.}

private

pred reso for proton hypothesis

Definition at line 220 of file HitLevelInfoWriter.h.

220{ -1.};

◆ m_prespi

double m_prespi { -1.}

private

pred reso for pion hypothesis

Definition at line 218 of file HitLevelInfoWriter.h.

218{ -1.};

◆ m_propertyFlags

unsigned int m_propertyFlags

privateinherited

The properties of the module as bitwise or (with |) of EModulePropFlags.

Definition at line 511 of file Module.h.

◆ m_pt

double m_pt { -1.}

private

transverse momentum from tracking

Definition at line 168 of file HitLevelInfoWriter.h.

168{ -1.};

◆ m_returnValue

int m_returnValue

privateinherited

The return value.

Definition at line 518 of file Module.h.

◆ m_runGain

double m_runGain { -1.}

private

calibration run gain

Definition at line 197 of file HitLevelInfoWriter.h.

197{ -1.};

◆ m_runID

int m_runID { -1}

private

run in which this Track was found

Definition at line 141 of file HitLevelInfoWriter.h.

141{ -1};

◆ m_scale

double m_scale { -1.}

private

calibration scale factor

Definition at line 194 of file HitLevelInfoWriter.h.

194{ -1.};

◆ m_strOutputBaseName

std::string m_strOutputBaseName

private

Base name for the output ROOT files.

Definition at line 117 of file HitLevelInfoWriter.h.

◆ m_strParticleList

std::vector<std::string> m_strParticleList

private

Vector of ParticleLists to write out.

Definition at line 118 of file HitLevelInfoWriter.h.

◆ m_tanlambda

double m_tanlambda { -1.}

private

Slope of the track in the r-z plane.

Definition at line 156 of file HitLevelInfoWriter.h.

156{ -1.};

◆ m_theta

double m_theta { -2.}

private

cos(theta) for the track

Definition at line 171 of file HitLevelInfoWriter.h.

171{ -2.};

◆ m_timeGain

double m_timeGain { -1.}

private

calibration injection time gain

Definition at line 198 of file HitLevelInfoWriter.h.

198{ -1.};

◆ m_timeReso

double m_timeReso { -1.}

private

calibration injection time reso

Definition at line 199 of file HitLevelInfoWriter.h.

199{ -1.};

◆ m_trackFitResults

StoreArray<TrackFitResult> m_trackFitResults

private

Required array of input TrackFitResults.

Definition at line 125 of file HitLevelInfoWriter.h.

◆ m_trackID

int m_trackID { -1}

private

ID number of the Track.

Definition at line 162 of file HitLevelInfoWriter.h.

162{ -1};

◆ m_tracks

StoreArray<Track> m_tracks

private

Required array of input Tracks.

Definition at line 124 of file HitLevelInfoWriter.h.

◆ m_tree

std::vector<TTree* > m_tree

private

output ROOT trees

Definition at line 121 of file HitLevelInfoWriter.h.

◆ m_trunc

double m_trunc { -1.}

private

dE/dx averaged, truncated mean, with corrections

Definition at line 181 of file HitLevelInfoWriter.h.

181{ -1.};

◆ m_truncNoSat

double m_truncNoSat { -1.}

private

dE/dx averaged, truncated mean, with corrections (not hadron)

Definition at line 182 of file HitLevelInfoWriter.h.

182{ -1.};

◆ m_TTDInfo

StoreObjPtr<EventLevelTriggerTimeInfo> m_TTDInfo

private

Store Object Ptr: EventLevelTriggerTimeInfo.

Definition at line 128 of file HitLevelInfoWriter.h.

◆ m_type

std::string m_type

privateinherited

The type of the module, saved as a string.

Definition at line 508 of file Module.h.

◆ m_vx0

double m_vx0 {0.}

private

X coordinate of track POCA to origin.

Definition at line 153 of file HitLevelInfoWriter.h.

153{0.};

◆ m_vy0

double m_vy0 {0.}

private

Y coordinate of track POCA to origin.

Definition at line 154 of file HitLevelInfoWriter.h.

154{0.};

◆ m_vz0

double m_vz0 {0.}

private

Z coordinate of track POCA to origin.

Definition at line 155 of file HitLevelInfoWriter.h.

155{0.};

◆ m_z0

double m_z0 {0.}

private

z coordinate of the POCA

Definition at line 149 of file HitLevelInfoWriter.h.

149{0.};

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

cdc/modules/HitLevelInfoWriter/include/HitLevelInfoWriter.h
cdc/modules/HitLevelInfoWriter/src/HitLevelInfoWriter.cc

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ EAfterConditionPath

Member Enumeration Documentation

◆ EModulePropFlags

Constructor & Destructor Documentation

◆ HitLevelInfoWriterModule()

◆ ~HitLevelInfoWriterModule()

Member Function Documentation

◆ beginRun()

◆ bookOutput()

◆ calculateMeans()

◆ clearEntries()

◆ clone()

◆ D2I()

◆ def_beginRun()

◆ def_endRun()

◆ def_event()

◆ def_initialize()

◆ def_terminate()

◆ endRun()

◆ evalCondition()

◆ event()

◆ exposePythonAPI()

◆ fillDedx()

◆ fillTrack()

◆ getAfterConditionPath()

◆ getAllConditionPaths()

◆ getAllConditions()

◆ getCondition()

◆ getConditionPath()

◆ GetCorrection()

◆ getDescription()

◆ getFileNames()

◆ getLogConfig()

◆ getModules()

◆ getName()

◆ getPackage()

◆ getParamInfoListPython()

◆ getParamList()

◆ getPathString()

◆ getReturnValue()

◆ getType()

◆ HadronCorrection()

◆ hasCondition()

◆ hasProperties()

◆ hasReturnValue()

◆ hasUnsetForcedParams()

◆ I2D()

◆ if_false()

◆ if_true()

◆ if_value()

◆ initialize()

◆ recalculateDedx()

◆ saveChiValue()

◆ setAbortLevel()

◆ setDebugLevel()

◆ setDescription()

◆ setLogConfig()

◆ setLogInfo()

◆ setLogLevel()

◆ setName()

◆ setParamList()

◆ setParamPython()

◆ setParamPythonDict()

◆ setPropertyFlags()

◆ setReturnValue() [1/2]

◆ setReturnValue() [2/2]

◆ setType()

◆ terminate()

Member Data Documentation

◆ h_adccorr

◆ h_adcraw