This module generates time vs crystal 2D histograms to later (in eclBhabhaTAlgorithm) find time crystal/crate offsets from bhabha events. More...

#include <ECLBhabhaTCollectorModule.h>

Inheritance diagram for ECLBhabhaTCollectorModule:

[legend]

Collaboration diagram for ECLBhabhaTCollectorModule:

[legend]

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
	ECLBhabhaTCollectorModule ()
	Module constructor.

virtual	~ECLBhabhaTCollectorModule ()
	Module destructor.

void	inDefineHisto ()
	Replacement for defineHisto() in CalibrationCollector modules.

void	prepare ()
	Define histograms and read payloads from DB.

void	collect ()
	Select events and crystals and accumulate histograms. More...

void	initialize () final
	Set up a default RunRange object in datastore and call prepare()

void	event () final
	Check current experiment and run and update if needed, fill into RunRange and collect()

void	beginRun () final
	Reset the m_runCollectOnRun flag, if necessary, to begin collection again. More...

void	endRun () final
	Write the current collector objects to a file and clear their memory.

void	terminate () final
	Write the final objects to the file.

void	defineHisto () final
	Runs due to HistoManager, allows us to discover the correct file.

template<class T >
void	registerObject (std::string name, T *obj)
	Register object with a name, takes ownership, do not access the pointer beyond prepare()

template<class T >
T *	getObjectPtr (std::string name)
	Calls the CalibObjManager to get the requested stored collector data.

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

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

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

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. More...

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

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. More...

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

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

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

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. More...

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). More...

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

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. More...

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. More...

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. More...

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

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

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

Protected Member Functions
virtual void	startRun ()
	Replacement for beginRun(). Do anything you would normally do in beginRun here.

virtual void	closeRun ()
	Replacement for endRun(). Do anything you would normally do in endRun here.

virtual void	finish ()
	Replacement for terminate(). Do anything you would normally do in terminate here.

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

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. More...

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. More...

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

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. More...

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

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

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

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

Protected Attributes
TDirectory *	m_dir
	The top TDirectory that collector objects for this collector will be stored beneath.

CalibObjManager	m_manager
	Controls the creation, collection and access to calibration objects.

RunRange *	m_runRange
	Overall list of runs processed.

Calibration::ExpRun	m_expRun
	Current ExpRun for object retrieval (becomes -1,-1 for granularity=all)

StoreObjPtr< EventMetaData >	m_emd
	Current EventMetaData.

Private Member Functions
bool	getPreScaleChoice ()
	I'm a little worried about floating point precision when comparing to 0.0 and 1.0 as special values. More...

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. More...

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

Private Attributes
bool	m_saveTree
	If true, save TTree with more detailed event info.

StoreArray< Track >	tracks
	StoreArray for tracks.

StoreObjPtr< EventT0 >	m_eventT0
	StoreObjPtr for T0. More...

DBObjPtr< ECLCrystalCalib >	m_ElectronicsDB
	electronics amplitude calibration from database Scale amplitudefor each crystal and for dead pre-amps More...

std::vector< float >	m_Electronics
	vector obtained from DB object

DBObjPtr< ECLCrystalCalib >	m_ElectronicsTimeDB
	Time offset from electronics calibration from database. More...

std::vector< float >	m_ElectronicsTime
	vector obtained from DB object

DBObjPtr< ECLCrystalCalib >	m_FlightTimeDB
	Time offset from flight time b/w IP and crystal from database. More...

std::vector< float >	m_FlightTime
	vector obtained from DB object

DBObjPtr< ECLCrystalCalib >	m_PreviousCrystalTimeDB
	Time offset from previous crystal time calibration (this calibration) from database. More...

std::vector< float >	m_PreviousCrystalTime
	vector obtained from DB object

std::vector< float >	m_PreviousCrystalTimeUnc
	vector obtained from DB object

DBObjPtr< ECLCrystalCalib >	m_CrateTimeDB
	Time offset from crate time calibration (also this calibration) from database. More...

std::vector< float >	m_CrateTime
	vector obtained from DB object

std::vector< float >	m_CrateTimeUnc
	uncertainty vector obtained from DB object

TTree *	m_dbgTree_electrons
	Output tree with detailed event data.

TTree *	m_dbgTree_tracks

TTree *	m_dbgTree_crystals

TTree *	m_dbgTree_event

TTree *	m_dbgTree_allCuts

TTree *	m_dbgTree_evt_allCuts

TTree *	m_dbgTree_crys_allCuts

int	m_tree_evtNum
	Event number for debug TTree output.

int	m_tree_cid
	ECL Cell ID (1..8736) for debug TTree output.

double	m_tree_phi
	phi position for debug TTree output

double	m_tree_theta
	theta position for debug TTree output

int	m_tree_amp
	Fitting amplitude from ECL for debug TTree output.

double	m_tree_en
	Energy of crystal with maximum energy within ECL cluster, GeV for debug TTree output.

double	m_tree_E1Etot
	Energy of crystal with maximum energy within ECL cluster divided by total cluster energy, unitless for debug TTree output.

double	m_tree_E1E2
	Energy of crystal with maximum energy within ECL cluster divided by second most energetic crystal in the cluster, unitless for debug TTree output.

double	m_tree_E1p
	Energy of crystal with maximum energy within ECL cluster divided by total cluster energy divided by the track momentum, unitless for debug TTree output.

int	m_tree_quality
	ECL fit quality for debug TTree output.

double	m_tree_timeF
	ECL fitting time for debug TTree output.

double	m_tree_time
	Time for Ts distribution for debug TTree output.

double	m_tree_timetsPreviousTimeCalibs
	Time for Ts distribution after application of previous time calibrations for debug TTree output.

double	m_tree_t0
	EventT0 (not from ECL) for debug TTree output.

double	m_tree_t0_unc
	EventT0 uncertainty for debug TTree output.

double	m_tree_t0_ECLclosestCDC
	EventT0 (from ECL) closest to CDC for debug TTree output.

double	m_tree_t0_ECL_minChi2
	EventT0 (from ECL) min chi2 for debug TTree output.

double	m_tree_d0

double	m_tree_z0
	Track d0 for debug TTree output.

double	m_tree_p
	Track z0 for debug TTree output.

double	m_tree_nCDChits
	Track momentum for debug TTree output.

double	m_tree_clustCrysE_DIV_maxEcrys
	Number of CDC hits along the track for debug TTree output.

double	m_tree_clustCrysE
	ratio of crystal energy to energy of the crystal that has the maximum energy, only for the crystals that meet all the selection criteria for debug TTree output

double	m_tree_enPlus
	crystal energy, only for the crystals that meet all the selection criteria for debug TTree output More...

double	m_tree_enNeg
	Energy of cluster associated to negatively charged track, GeV for debug TTree output.

double	m_tree_tClustPos
	Cluster time of cluster associated to positively charged track, ns for debug TTree output.

double	m_tree_tClustNeg
	Cluster time of cluster associated to negatively charged track, ns for debug TTree output.

double	m_tree_maxEcrystPosClust
	Time of the highest energy crystal in the cluster associated to positively charged track, ns for debug TTree output.

double	m_tree_maxEcrystNegClust
	Time of the highest energy crystal in the cluster associated to negatively charged track, ns for debug TTree output.

double	m_tree_tClust
	Cluster time of a cluster, ns for debug TTree output.

double	m_tree_ECLCalDigitTime
	Time of an ECLCalDigit within a cluster, ns for debug TTree output.

double	m_tree_ECLCalDigitE
	Energy of an ECLCalDigit within a cluster, GeV for debug TTree output.

double	m_tree_ECLDigitAmplitude
	Amplitude (used to calculate energy) of an ECLDigit within a cluster, for debug TTree output.

int	m_charge
	particle charge, for debug TTree output

double	m_E_DIV_p
	Energy divided by momentum, for debug TTree output.

double	m_massInvTracks
	invariant mass of the two tracks, for debug TTree output

StoreArray< ECLDigit >	m_eclDigitArray
	Required input array of ECLDigits.

StoreArray< ECLCalDigit >	m_eclCalDigitArray
	Required input array of ECLCalDigits.

StoreArray< ECLCluster >	m_eclClusterArray
	Required input array of ECLClusters.

std::vector< float >	m_EperCrys
	ECL digit energy for each crystal.

std::vector< int >	m_eclCalDigitID
	ECL cal digit id sorter.

std::vector< int >	m_eclDigitID
	ECL digit id sorter.

short	m_timeAbsMax
	Events with abs(time) > m_timeAbsMax are excluded, mostly for histogram x-range purposes.

int	m_minCrystal
	First CellId to handle.

int	m_maxCrystal
	Last CellId to handle.

double	m_looseTrkZ0

double	m_tightTrkZ0

double	m_looseTrkD0

double	m_tightTrkD0

int	m_crystalCrate
	Crate id for the crystal.

int	m_runNum
	run number

bool	m_storeCalib = true
	Boolean for whether or not to store the previous calibration calibration constants.

std::unique_ptr< Belle2::ECL::ECLTimingUtilities >	m_ECLTimeUtil
	ECL timing tools. More...

double	m_energyDependenceTimeOffsetFitParam_p1 = 0
	p1 in "energy dependence equation"

double	m_energyDependenceTimeOffsetFitParam_p2 = 88449.
	p2 in "energy dependence equation"

double	m_energyDependenceTimeOffsetFitParam_p3 = 0.20867E+06
	p3 in "energy dependence equation"

double	m_energyDependenceTimeOffsetFitParam_p4 = 3.1482
	p4 in "energy dependence equation"

double	m_energyDependenceTimeOffsetFitParam_p5 = 7.4747
	p5 in "energy dependence equation"

double	m_energyDependenceTimeOffsetFitParam_p6 = 1279.3
	p6 in "energy dependence equation"

std::string	m_granularity
	Granularity of data collection = run\|all(= no granularity, exp,run=-1,-1)

int	m_maxEventsPerRun
	Maximum number of events to be collected at the start of each run (-1 = no maximum)

float	m_preScale
	Prescale module parameter, this fraction of events will have collect() run on them [0.0 -> 1.0].

StoreObjPtr< EventMetaData >	m_evtMetaData
	Required input for EventMetaData.

bool	m_runCollectOnRun = true
	Whether or not we will run the collect() at all this run, basically skips the event() function if false.

std::map< Calibration::ExpRun, int >	m_expRunEvents
	How many events processed for each ExpRun so far, stops counting up once max is hit Only used/incremented if m_maxEventsPerRun > -1.

int *	m_eventsCollectedInRun
	Will point at correct value in m_expRunEvents.

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.

Detailed Description

This module generates time vs crystal 2D histograms to later (in eclBhabhaTAlgorithm) find time crystal/crate offsets from bhabha events.

Definition at line 45 of file ECLBhabhaTCollectorModule.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 79 of file Module.h.

Member Function Documentation

◆ beginRun()

void beginRun ( )

finalvirtualinherited

Reset the m_runCollectOnRun flag, if necessary, to begin collection again.

It seems that the beginRun() function is called in each basf2 subprocess when the run changes in each process. This is nice because it allows us to write the new (exp,run) object creation in the beginRun function as though the other processes don't exist.

Reimplemented from HistoModule.

Definition at line 70 of file CalibrationCollectorModule.cc.

◆ 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 181 of file Module.cc.

◆ collect()

void collect ( )

virtual

Select events and crystals and accumulate histograms.

< vector derived from DB object

Record the input database constants for the first call

< number of loose tracks

< number of tight tracks

Reimplemented from CalibrationCollectorModule.

Definition at line 324 of file ECLBhabhaTCollectorModule.cc.

 {
   int cutIndexPassed = 0;
   getObjectPtr<TH1F>("cutflow")->Fill(cutIndexPassed);
   B2DEBUG(10, "Cutflow: no cuts: index = " << cutIndexPassed);
  
  
   /* Use ECLChannelMapper to get other detector indices for the crystals */
   /* For conversion from CellID to crate, shaper, and channel ids. */
  
   // Use smart pointer to avoid memory leak when the ECLChannelMapper object needs destroying at the end of the event.
   shared_ptr< ECL::ECLChannelMapper > crystalMapper(new ECL::ECLChannelMapper());
   crystalMapper->initFromDB();
  
   //== Get expected energies and calibration constants from DB. Need to call
   //   hasChanged() for later comparison
   if (m_ElectronicsDB.hasChanged()) {
     m_Electronics = m_ElectronicsDB->getCalibVector();
   }
   if (m_ElectronicsTimeDB.hasChanged()) {
     m_ElectronicsTime = m_ElectronicsTimeDB->getCalibVector();
   }
   if (m_FlightTimeDB.hasChanged()) {
     m_FlightTime = m_FlightTimeDB->getCalibVector();
   }
  
   // Get the previous crystal time offset (the same thing that this calibration is meant to calculate).
   // This can be used for testing purposes, and for the crate time offset.
   if (m_PreviousCrystalTimeDB.hasChanged()) {
     m_PreviousCrystalTime = m_PreviousCrystalTimeDB->getCalibVector();
     m_PreviousCrystalTimeUnc = m_PreviousCrystalTimeDB->getCalibUncVector();
   }
  
   B2DEBUG(35, "Finished checking if previous crystal time payload has changed");
   if (m_CrateTimeDB.hasChanged()) {
     m_CrateTime = m_CrateTimeDB->getCalibVector();
     m_CrateTimeUnc = m_CrateTimeDB->getCalibUncVector();
   }
   B2DEBUG(35, "Finished checking if previous crate time payload has changed");
   B2DEBUG(35, "m_CrateTime size = " << m_CrateTime.size());
   B2DEBUG(29, "Crate time +- uncertainty [0]= " << m_CrateTime[0] << " +- " << m_CrateTimeUnc[0]);
   B2DEBUG(29, "Crate time +- uncertainty [8735]= " << m_CrateTime[8735] << " +- " << m_CrateTimeUnc[8735]);
  
  
   // Conversion coefficient from ADC ticks to nanoseconds
   // TICKS_TO_NS ~ 0.4931 ns/clock tick
   // 1/(4fRF) = 0.4913 ns/clock tick, where fRF is the accelerator RF frequency, fRF=508.889 MHz.
   const double TICKS_TO_NS = 1.0 / (4.0 * EclConfiguration::m_rf) * 1e3;
  
  
   vector<float> Crate_time_ns(52, 0.0); 
   // Make a crate time offset vector with an entry per crate (instead of per crystal) and convert from ADC counts to ns.
   for (int crysID = 1; crysID <= 8736; crysID++) {
     int crateID_temp = crystalMapper->getCrateID(crysID);
     Crate_time_ns[crateID_temp - 1] = m_CrateTime[crysID] * TICKS_TO_NS;
   }
  
  
   if (m_storeCalib) {
     for (int crysID = 0; crysID < 8736; crysID++) {
       getObjectPtr<TH1F>("TsDatabase")->Fill(crysID + 0.001, m_PreviousCrystalTime[crysID]);
       getObjectPtr<TH1F>("TsDatabaseUnc")->Fill(crysID + 0.001, m_PreviousCrystalTimeUnc[crysID]);
       getObjectPtr<TH1F>("TcrateDatabase")->Fill(crysID + 0.001, m_CrateTime[crysID]);
       getObjectPtr<TH1F>("TcrateUncDatabase")->Fill(crysID + 0.001, m_CrateTimeUnc[crysID]);
  
       B2INFO("cid = " << crysID + 1 << ", Ts previous = " << m_PreviousCrystalTime[crysID]);
     }
  
     for (int crateID_temp = 1; crateID_temp <= 52; crateID_temp++) {
       getObjectPtr<TH1F>("tcrateDatabase_ns")->Fill(crateID_temp + 0.001, Crate_time_ns[crateID_temp - 1]);
     }
  
     // Use a histogram with only one bin as a counter to know the number of times the database histograms were filled.
     //    This is mostly useful for the talg when running over multiple runs and trying to read ts values.
     getObjectPtr<TH1F>("databaseCounter")->Fill(1.001, 1);
  
     m_storeCalib = false;
   }
  
  
   /* Getting the event t0 using the full event t0 rather than from the CDC specifically */
   double evt_t0 = -1;
   double evt_t0_unc = -1;
   double evt_t0_ECL_closestCDC = -1;
   double evt_t0_ECL_minChi2 = -1;
  
   // Determine if there is an event t0 to use and then extract the information about it
   if (!m_eventT0.isValid()) {
     //cout << "event t0 not valid\n";
     return;
   } else if (!m_eventT0->hasTemporaryEventT0(Const::EDetector::CDC)) {
     //cout << "no event t0\n";
     return;
   } else {
     // Event has a t0 from CDC
     cutIndexPassed++;
     getObjectPtr<TH1F>("cutflow")->Fill(cutIndexPassed);
     B2DEBUG(10, "Cutflow: Event t0 exists: index = " << cutIndexPassed);
  
  
     // Get event t0 from CDC.  We don't want event t0 from ECL as we are calibrating the ECL wrt the more accurately measured time measurements of the time.  Start with the CDC since it has an event t0 but in the future we may switch to the TOP detector.
     // Based on the information from Thomas Hauth <Thomas.Hauth@kit.edu> (leaving physics) we should take the last event t0 in the list of event t0's from the CDC as the later event t0 measurements are calculated in slower but more accurate ways.
     vector<EventT0::EventT0Component> evt_t0_list = m_eventT0->getTemporaryEventT0s(Const::EDetector::CDC);
     evt_t0 = evt_t0_list.back().eventT0;   // time value
     evt_t0_unc = evt_t0_list.back().eventT0Uncertainty;   // uncertainty on event t0
  
  
     // Get the ECL event t0 for comparison - validations
     if (m_eventT0->hasTemporaryEventT0(Const::EDetector::ECL)) {
       vector<EventT0::EventT0Component> evt_t0_list_ECL = m_eventT0->getTemporaryEventT0s(Const::EDetector::ECL);
  
  
       double smallest_CDC_ECL_t0_diff = fabs(evt_t0_list_ECL[0].eventT0 - evt_t0);
       int smallest_CDC_ECL_t0_diff_idx = 0;
       for (long unsigned int ECLi = 0; ECLi < evt_t0_list_ECL.size(); ECLi++) {
         double tempt_ECL_t0 = evt_t0_list_ECL[ECLi].eventT0;
         if (fabs(tempt_ECL_t0 - evt_t0) < smallest_CDC_ECL_t0_diff) {
           smallest_CDC_ECL_t0_diff = fabs(tempt_ECL_t0 - evt_t0);
           smallest_CDC_ECL_t0_diff_idx = ECLi;
         }
       }
  
       evt_t0_ECL_closestCDC = evt_t0_list_ECL[smallest_CDC_ECL_t0_diff_idx].eventT0;   // time value
       B2DEBUG(30, "evt_t0_ECL_closestCDC = " << evt_t0_ECL_closestCDC);
  
  
  
       double smallest_ECL_t0_minChi2 = evt_t0_list_ECL[0].quality;
       int smallest_ECL_t0_minChi2_idx = 0;
  
       B2DEBUG(30, "evt_t0_list_ECL[0].quality = " << evt_t0_list_ECL[0].quality
               << ", with ECL event t0 = " << evt_t0_list_ECL[0].eventT0);
  
       for (long unsigned int ECLi = 0; ECLi < evt_t0_list_ECL.size(); ECLi++) {
         B2DEBUG(30, "evt_t0_list_ECL[" << ECLi << "].quality = " << evt_t0_list_ECL[ECLi].quality
                 << ", with ECL event t0 = " <<
                 evt_t0_list_ECL[ECLi].eventT0);
         if (evt_t0_list_ECL[ECLi].quality < smallest_ECL_t0_minChi2) {
           smallest_ECL_t0_minChi2 = evt_t0_list_ECL[ECLi].quality;
           smallest_ECL_t0_minChi2_idx = ECLi;
         }
       }
  
       evt_t0_ECL_minChi2 = evt_t0_list_ECL[smallest_ECL_t0_minChi2_idx].eventT0;   // time value
  
       B2DEBUG(30, "evt_t0_ECL_minChi2 = " << evt_t0_ECL_minChi2);
       B2DEBUG(30, "smallest_ECL_t0_minChi2_idx = " << smallest_ECL_t0_minChi2_idx);
     }
   }
  
  
  
   /* Determine the energies for each of the crystals since this isn't naturally connected to the cluster.
      Also determine the indexing of the ecl cal digits and the ecl digits
      Taken from Chris's ec/modules/eclGammaGammaECollector   */
  
   // Resize vectors
   m_EperCrys.resize(8736);
   m_eclCalDigitID.resize(8736);
   m_eclDigitID.resize(8736);
  
  
   int idx = 0;
   for (auto& eclCalDigit : m_eclCalDigitArray) {
     int tempCrysID = eclCalDigit.getCellId() - 1;
     m_EperCrys[tempCrysID] = eclCalDigit.getEnergy();
     m_eclCalDigitID[tempCrysID] = idx;
     idx++;
   }
  
   idx = 0;
   for (auto& eclDigit : m_eclDigitArray) {
     int tempCrysID = eclDigit.getCellId() - 1;
     m_eclDigitID[tempCrysID] = idx;
     idx++;
   }
  
  
  
  
   //---------------------------------------------------------------------
   //..Some utilities
   ClusterUtils cUtil;
   const TVector3 clustervertex = cUtil.GetIPPosition();
   PCmsLabTransform boostrotate;
  
   //---------------------------------------------------------------------
   //..Track properties, including 2 maxp tracks. Use pion (211) mass hypothesis,
   //     which is the only particle hypothesis currently available???
   double maxp[2] = {0., 0.};
   int maxiTrk[2] = { -1, -1};
   int nTrkAll = tracks.getEntries();
  
   int nTrkLoose = 0; 
   int nTrkTight = 0; 
   /* Loop over all the tracks to define the tight and loose selection tracks.
      We will select events with only 2 tight tracks and no additional loose tracks.
      Tight tracks are a subset of looses tracks. */
   for (int iTrk = 0; iTrk < nTrkAll; iTrk++) {
     // Get track and assume it is a pion for now ... because it is the only particle we can assume?
     const TrackFitResult* tempTrackFit = tracks[iTrk]->getTrackFitResult(Const::ChargedStable(211));
     if (not tempTrackFit) {continue;}
  
     // Collect track info to be used for categorizing
     short charge = tempTrackFit->getChargeSign();
     double z0 = tempTrackFit->getZ0();
     double d0 = tempTrackFit->getD0();
     int nCDChits = tempTrackFit->getHitPatternCDC().getNHits();
     double p = tempTrackFit->getMomentum().Mag();
  
     // >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     //== Save debug TTree with detailed information if necessary.
     m_tree_d0 = d0;
     m_tree_z0 = z0;
     m_tree_p  = p;
     m_charge = charge;
     m_tree_nCDChits = nCDChits;
  
     if (m_saveTree) {
       m_dbgTree_tracks->Fill();
     }
     //<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
  
  
     /* Test if loose track  */
  
     // d0 and z0 cuts
     if (abs(d0) > m_looseTrkD0) {
       continue;
     }
     if (abs(z0) > m_looseTrkZ0) {
       continue;
     }
     // Number of hits in the CDC
     if (nCDChits < 1) {
       continue;
     }
     nTrkLoose++;
  
  
  
     /* Test if the loose track is also a tight track */
  
     // Number of hits in the CDC
     if (nCDChits < 20) {
       continue;
     }
  
  
     // d0 and z0 cuts
     if (abs(d0) > m_tightTrkD0) {
       continue;
     }
     if (abs(z0) > m_tightTrkZ0) {
       continue;
     }
     nTrkTight++;
  
     // Sorting of tight tracks.  Not really required as we only want two tight tracks (at the moment) but okay.
     //..Find the maximum p negative [0] and positive [1] tracks
     int icharge = 0;
     if (charge > 0) {icharge = 1;}
     if (p > maxp[icharge]) {
       maxp[icharge] = p;
       maxiTrk[icharge] = iTrk;
     }
  
   }
   /* After that last section the numbers of loose and tight tracks are known as well as the
      index of the loose tracks that have the highest p negatively charged and highest p positively
      charged tracks as measured in the centre of mass frame */
  
  
   if (nTrkTight != 2) {
     return;
   }
   // There are exactly two tight tracks
   cutIndexPassed++;
   getObjectPtr<TH1F>("cutflow")->Fill(cutIndexPassed);
   B2DEBUG(10, "Cutflow: Two tight tracks: index = " << cutIndexPassed);
  
  
   if (nTrkLoose != 2) {
     return;
   }
   // There are exactly two loose tracks as well, i.e. no additional loose tracks
   cutIndexPassed++;
   getObjectPtr<TH1F>("cutflow")->Fill(cutIndexPassed);
   B2DEBUG(10, "Cutflow: No additional loose tracks: index = " << cutIndexPassed);
  
   /* Determine if the two tracks have the opposite electric charge.
      We know this because the track indices stores the max pt track in [0] for negatively charged track
      and [1] fo the positively charged track.  If both are filled then both a negatively charged
      and positively charged track were found.   */
   bool oppositelyChargedTracksPassed = maxiTrk[0] != -1  &&  maxiTrk[1] != -1;
   if (!oppositelyChargedTracksPassed) {
     return;
   }
   // The two tracks have the opposite electric charges.
   cutIndexPassed++;
   getObjectPtr<TH1F>("cutflow")->Fill(cutIndexPassed);
   B2DEBUG(10, "Cutflow: Oppositely charged tracks: index = " << cutIndexPassed);
  
  
  
  
   //---------------------------------------------------------------------
   /* Determine associated energy clusters to each of the two tracks.  Sum the energies of the
      multiple clusters to each track and find the crystal with the maximum energy within all
      the sets of clusters associated to the tracks*/
   double trkEClustLab[2] = {0., 0.};
   double trkEClustCOM[2] = {0., 0.};
   double trkpLab[2];
   double trkpCOM[2];
   //double trkThetaLab[2];
   TLorentzVector trkp4Lab[2];
   TLorentzVector trkp4COM[2];
  
   // Index of the cluster and the crystal that has the highest energy crystal for the two tracks
   int crysIDMax[2] = { -1, -1 };
   double crysEMax[2] = { -1, -1 };
   double crysE2Max[2] = { -1, -1 };
   int numClustersPerTrack[2] = { 0, 0 };
  
   double clusterTime[2] = {0, 0};
  
   double E_DIV_p[2];
  
   vector<double> time_ECLCaldigits_bothClusters;
   vector<int> cid_ECLCaldigits_bothClusters;
   vector<double> E_ECLCaldigits_bothClusters;
   vector<double> amp_ECLDigits_bothClusters;
   vector<int> chargeID_ECLCaldigits_bothClusters;
  
   for (int icharge = 0; icharge < 2; icharge++) {
     if (maxiTrk[icharge] > -1) {
       B2DEBUG(10, "looping over the 2 max pt tracks");
  
       const TrackFitResult* tempTrackFit = tracks[maxiTrk[icharge]]->getTrackFitResult(Const::ChargedStable(211));
       trkp4Lab[icharge] = tempTrackFit->get4Momentum();
       trkp4COM[icharge] = boostrotate.rotateLabToCms() * trkp4Lab[icharge];
       trkpLab[icharge] = trkp4Lab[icharge].Rho();
       trkpCOM[icharge] = trkp4COM[icharge].Rho();
  
  
       /* For each cluster associated to the current track, sum up the energies to get the total
          energy of all clusters associated to the track and find which crystal has the highest
          energy from all those clusters*/
       auto eclClusterRelationsFromTracks = tracks[maxiTrk[icharge]]->getRelationsTo<ECLCluster>();
       for (unsigned int clusterIdx = 0; clusterIdx < eclClusterRelationsFromTracks.size(); clusterIdx++) {
  
         B2DEBUG(10, "Looking at clusters.  index = " << clusterIdx);
         auto cluster = eclClusterRelationsFromTracks[clusterIdx];
         bool goodClusterType = false;
  
         if (cluster->hasHypothesis(Belle2::ECLCluster::EHypothesisBit::c_nPhotons)) {
           trkEClustLab[icharge] += cluster->getEnergy(Belle2::ECLCluster::EHypothesisBit::c_nPhotons);
           goodClusterType = true;
           numClustersPerTrack[icharge]++;
         }
  
         if (goodClusterType) {
  
           clusterTime[icharge] = cluster->getTime();
  
           auto eclClusterRelations = cluster->getRelationsTo<ECLCalDigit>("ECLCalDigits");
  
           // Find the crystal that has the largest energy
           for (unsigned int ir = 0; ir < eclClusterRelations.size(); ir++) {
             const auto calDigit = eclClusterRelations.object(ir);
             int tempCrysID = calDigit->getCellId() - 1;
             double tempE = m_EperCrys[tempCrysID];
  
             int eclDigitIndex = m_eclDigitID[tempCrysID];
             ECLDigit*    ecl_dig = m_eclDigitArray[eclDigitIndex];
  
             // for the max E crystal
             if (tempE > crysEMax[icharge]) {
               // Set 2nd highest E crystal to the info from the highest E crystal
               crysE2Max[icharge] = crysEMax[icharge];
               // Set the highest E crystal to the current crystal
               crysEMax[icharge] = tempE;
               crysIDMax[icharge] = tempCrysID;
             }
             // for the 2nd highest E crystal
             if (tempE > crysE2Max[icharge] && tempCrysID != crysIDMax[icharge]) {
               crysE2Max[icharge] = tempE;
             }
  
  
             B2DEBUG(30,  "calDigit(ir" << ir << ") time = " << calDigit->getTime() << "ns , with E = " << tempE << " GeV");
             time_ECLCaldigits_bothClusters.push_back(calDigit->getTime());
             cid_ECLCaldigits_bothClusters.push_back(tempCrysID);
             E_ECLCaldigits_bothClusters.push_back(tempE);
             amp_ECLDigits_bothClusters.push_back(ecl_dig->getAmp());
             chargeID_ECLCaldigits_bothClusters.push_back(icharge);
  
           }
         }
       }
       trkEClustCOM[icharge] = trkEClustLab[icharge] * trkpCOM[icharge] / trkpLab[icharge];
  
       // Send to the histogram the number of photon clusters associated to the track
       getObjectPtr<TH1F>("numPhotonClustersPerTrack")->Fill(numClustersPerTrack[icharge]);
  
  
       // Check both electrons to see if their cluster energy / track momentum is good.
       // The Belle II physics book shows that this is the main way of separating electrons from other particles
       // Done in the centre of mass reference frame although I believe E/p is invariant under a boost.
       E_DIV_p[icharge] = trkEClustCOM[icharge] / trkpCOM[icharge];
  
     }
   }
   /* At the end of this section the 3-momenta magnitudes and the cluster energies are known
      for the two saved track indices for both the lab and COM frames.
      The crystal with the maximum energy, one associated to each track, is recorded*/
  
  
  
   //=== Check each crystal in the processed event and fill histogram.
  
   int cid;
   double time;
   int numCrystalsPassingCuts = 0;
  
   int crystalIDs[2] = { -1, -1};
   int crateIDs[2] = { -1, -1};
   double ts_prevCalib[2] = { -1, -1};
   double tcrate_prevCalib[2] = { -1, -1};
   double times_noPreviousCalibrations[2] = { -1, -1};
   bool crystalCutsPassed[2] = {false, false};
   double crystalEnergies[2] = { -1, -1};
   double crystalEnergies2[2] = { -1, -1};
  
   for (int iCharge = 0; iCharge < 2; iCharge++) {
     int crystal_idx = crysIDMax[iCharge];
     int eclCalDigitIndex = m_eclCalDigitID[crystal_idx];
     int eclDigitIndex = m_eclDigitID[crystal_idx];
  
     ECLDigit*    ecl_dig = m_eclDigitArray[eclDigitIndex];
     ECLCalDigit* ecl_cal = m_eclCalDigitArray[eclCalDigitIndex];
  
     //== Check whether specific ECLDigits should be excluded.
  
     auto en = ecl_cal->getEnergy();
     auto amplitude = ecl_dig->getAmp();
     crystalEnergies[iCharge] = en;
  
     cid   = ecl_dig->getCellId();
     time  = ecl_dig->getTimeFit() * TICKS_TO_NS - evt_t0;
  
     // Offset time by electronics calibration and flight time calibration.
     time -= m_ElectronicsTime[cid - 1] * TICKS_TO_NS;
     time -= m_FlightTime[cid - 1];
  
  
     // Apply the time walk correction: time shift as a function of the amplitude corrected by the electronics calibration.
     //    The electronics calibration also accounts for crystals that have a dead pre-amp and thus half the normal amplitude.
     double energyTimeShift = m_ECLTimeUtil->energyDependentTimeOffsetElectronic(amplitude * m_Electronics[cid - 1]) * TICKS_TO_NS;
  
     B2DEBUG(35, "cellid = " << cid << ", amplitude = " << amplitude << ", time before t(E) shift = " << time <<
             ", t(E) shift = " << energyTimeShift << " ns");
     time -= energyTimeShift;
  
  
     // Cell ID should be within specified range.
     if (cid < m_minCrystal || cid > m_maxCrystal) continue;
  
     // Absolute time should be in specified range condition.
     if (fabs(time) > m_timeAbsMax) continue;
  
     // Fit quality flag -- choose only events with best fit quality
     if (ecl_dig->getQuality() != 0) continue;
  
     //== Save time and crystal information.  Fill plot after both electrons are tested
     crystalIDs[iCharge] = cid;
     crateIDs[iCharge] = crystalMapper->getCrateID(ecl_cal->getCellId());
  
  
     ts_prevCalib[iCharge] = m_PreviousCrystalTime[cid - 1] * TICKS_TO_NS;
     tcrate_prevCalib[iCharge] = m_CrateTime[cid - 1] * TICKS_TO_NS;
     times_noPreviousCalibrations[iCharge] = time;
  
  
     B2DEBUG(30, "iCharge = " << iCharge);
     B2DEBUG(30, "crateIDs[iCharge] = " << crateIDs[iCharge]);
     B2DEBUG(30, "times_noPreviousCalibrations[iCharge] = " << times_noPreviousCalibrations[iCharge]);
     B2DEBUG(30, "tcrate_prevCalib[iCharge] = " << tcrate_prevCalib[iCharge]);
     B2DEBUG(30, "ts_prevCalib[iCharge] = " << ts_prevCalib[iCharge]);
  
  
     crystalCutsPassed[iCharge] = true;
  
  
     // For second most energetic energy crystal
     crystalEnergies2[iCharge] = crysE2Max[iCharge];
  
  
 //  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Tree saving
     //== Save debug TTree with detailed information if necessary.
     StoreObjPtr<EventMetaData> evtMetaData;
  
     m_tree_cid      = ecl_dig->getCellId();
     m_tree_amp      = ecl_dig->getAmp();
     m_tree_en       = en;
     m_tree_E1Etot   = en / trkEClustLab[iCharge];
     m_tree_E1E2     = en / crystalEnergies2[iCharge];
     m_tree_E1p      = en / trkpLab[iCharge];
     m_tree_timetsPreviousTimeCalibs = time - ts_prevCalib[iCharge] - tcrate_prevCalib[iCharge];
     m_tree_timeF    = ecl_dig->getTimeFit() * TICKS_TO_NS;
     m_tree_time     = time;
     m_tree_quality  = ecl_dig->getQuality();
     m_tree_t0       = evt_t0;
     m_tree_t0_unc   = evt_t0_unc;
     m_E_DIV_p       = E_DIV_p[iCharge];
     m_tree_evtNum  = evtMetaData->getEvent();
     m_crystalCrate  = crystalMapper->getCrateID(ecl_cal->getCellId());
     m_runNum        = evtMetaData->getRun();
  
     if (m_saveTree) {
       m_dbgTree_electrons->Fill();
     }
 //  <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Tree saving
  
     // Fill histogram with information about maximum energy crystal energy fraction
     getObjectPtr<TH1F>("maxEcrsytalEnergyFraction")->Fill(en / trkEClustLab[iCharge]);
  
  
   }
  
  
  
   // Check both electrons to see if their cluster energy / track momentum is good.
   // The Belle II physics book shows that this is the main way of separating electrons from other particles
   // Done in the centre of mass reference frame although I believe E/p is invariant under a boost.
   bool E_DIV_p_instance_passed[2] = {false, false};
   double E_DIV_p_CUT = 0.7;
   for (int icharge = 0; icharge < 2; icharge++) {
     E_DIV_p_instance_passed[icharge] = E_DIV_p[icharge] > E_DIV_p_CUT;
   }
   if (!E_DIV_p_instance_passed[0] || !E_DIV_p_instance_passed[1]) {
     return;
   }
   // E/p sufficiently large
   cutIndexPassed++;
   getObjectPtr<TH1F>("cutflow")->Fill(cutIndexPassed);
   B2DEBUG(10, "Cutflow: E_i/p_i > " << E_DIV_p_CUT  << ": index = " << cutIndexPassed);
  
  
  
  
  
  
   // Now find energy clusters independently of the tracks
   //------------------------------------------------------------------------
  
   vector<double> goodClusters_E;
  
  
  
  
   double clusterE_minCut = 0.06;  // GeV
   int nclust = m_eclClusterArray.getEntries();
   int nGoodClusts = 0;
   vector<int> goodPhotonClusterIdxs;
   vector<int> goodECLClusterIds;
   for (int ic = 0; ic < nclust; ic++) {
     if (m_eclClusterArray[ic]->hasHypothesis(Belle2::ECLCluster::EHypothesisBit::c_nPhotons)) {
       double eClust = m_eclClusterArray[ic]->getEnergy(Belle2::ECLCluster::EHypothesisBit::c_nPhotons);
       if (eClust > clusterE_minCut) {
         goodPhotonClusterIdxs.push_back(ic);
         goodECLClusterIds.push_back(m_eclClusterArray[ic]->getClusterId());
         nGoodClusts++;
       }
     }
   }
  
   int numTrackless = 0;
   for (int clustId = 0; clustId < 2; clustId++) {
     B2DEBUG(30, "m_eclClusterArray[goodPhotonClusterIdxs[clustId]]->isTrack() = " <<
             m_eclClusterArray[goodPhotonClusterIdxs[clustId]]->isTrack());
     if (!m_eclClusterArray[goodPhotonClusterIdxs[clustId]]->isTrack()) {
       numTrackless++;
     }
   }
   if (numTrackless != 0) {
     B2DEBUG(20, "Number of trackless ECL clusters != 0");
   }
  
   // There are exactly three energy clusters
   cutIndexPassed++;
   getObjectPtr<TH1F>("cutflow")->Fill(cutIndexPassed);
   B2DEBUG(10, "Cutflow: ECL cluster associated to photon does not have a track associated to it: index = " << cutIndexPassed);
   B2DEBUG(10, "Cutflow: NEW CUT TO BHABHA CODE FROM eegamma code");
  
  
  
  
  
   // Start of cuts on both the combined system of tracks and energy clusters
  
   double invMassTrk = (trkp4Lab[0] + trkp4Lab[1]).M();
   double invMass_CUT = 0.9;
   m_massInvTracks = invMassTrk;   // invariant mass of the two tracks
  
 //  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Tree saving
   if (m_saveTree) {
     m_dbgTree_event->Fill();
   }
 //  <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Tree saving
  
   bool invMassCutsPassed = invMassTrk > (invMass_CUT * boostrotate.getCMSEnergy());
   if (!invMassCutsPassed) {
     return;
   }
   // Invariable mass of the two tracks are above the minimum
   cutIndexPassed++;
   getObjectPtr<TH1F>("cutflow")->Fill(cutIndexPassed);
   B2DEBUG(10, "Cutflow: m(track 1+2) > " << invMass_CUT << "*E_COM = " << invMass_CUT << " * " << boostrotate.getCMSEnergy() <<
           " : index = " << cutIndexPassed);
  
  
  
   //== Fill output histogram.
   for (int iCharge = 0; iCharge < 2; iCharge++) {
     if (crystalCutsPassed[iCharge]) {
       getObjectPtr<TH2F>("TimevsCrysPrevCrateCalibPrevCrystCalib")->Fill((crystalIDs[iCharge] - 1) + 0.001,
           times_noPreviousCalibrations[iCharge] - ts_prevCalib[iCharge] - tcrate_prevCalib[iCharge] , 1);
       getObjectPtr<TH2F>("TimevsCratePrevCrateCalibPrevCrystCalib")->Fill((crateIDs[iCharge]) + 0.001,
           times_noPreviousCalibrations[iCharge]  - ts_prevCalib[iCharge] - tcrate_prevCalib[iCharge], 1);
       getObjectPtr<TH2F>("TimevsCrysNoCalibrations")->Fill((crystalIDs[iCharge] - 1) + 0.001, times_noPreviousCalibrations[iCharge], 1);
       getObjectPtr<TH2F>("TimevsCrateNoCalibrations")->Fill((crateIDs[iCharge]) + 0.001, times_noPreviousCalibrations[iCharge], 1);
       getObjectPtr<TH2F>("TimevsCrysPrevCrateCalibNoCrystCalib")->Fill((crystalIDs[iCharge] - 1) + 0.001,
           times_noPreviousCalibrations[iCharge] - tcrate_prevCalib[iCharge], 1);
       getObjectPtr<TH2F>("TimevsCrateNoCrateCalibPrevCrystCalib")->Fill((crateIDs[iCharge]) + 0.001,
           times_noPreviousCalibrations[iCharge]  - ts_prevCalib[iCharge] , 1);
  
       getObjectPtr<TH1F>("numEntriesPerCrystal")->Fill((crystalIDs[iCharge] - 1) + 0.001);
  
       // Record number of crystals used from the event.  Should be exactly two.
       numCrystalsPassingCuts++;
  
  
       /* Store information about one random specific crystal and one random specific crate for
          studies of how they vary over time.*/
       if (crystalIDs[iCharge] == 1338) {
         getObjectPtr<TH2F>("TimevsRunPrevCrateCalibPrevCrystCalibCID1338")->Fill(m_runNum + 0.001,
             times_noPreviousCalibrations[iCharge] - ts_prevCalib[iCharge] - tcrate_prevCalib[iCharge] , 1);
       }
       if (crateIDs[iCharge] == 8) {
         getObjectPtr<TH2F>("TimevsRunPrevCrateCalibPrevCrystCalibCrate8")->Fill(m_runNum + 0.001,
             times_noPreviousCalibrations[iCharge] - ts_prevCalib[iCharge] - tcrate_prevCalib[iCharge] , 1);
       }
  
     }
   }
  
  
   // Change cutflow method for this bit ... don't call return because we used to call the hadron cluster stuff afterwards
   //
   if (crystalCutsPassed[0] || crystalCutsPassed[1]) {
     // At least one ECL crystal time and quality cuts passed
     cutIndexPassed++;
     getObjectPtr<TH1F>("cutflow")->Fill(cutIndexPassed);
     B2DEBUG(10, "Cutflow: At least one crystal time and quality cuts passed: index = " << cutIndexPassed);
  
     getObjectPtr<TH1F>("numCrystalEntriesPerEvent")->Fill(numCrystalsPassingCuts);
   }
  
  
   // Save final information to the tree after all cuts are applied
   for (int iCharge = 0; iCharge < 2; iCharge++) {
     if (crystalCutsPassed[iCharge]) {
       //  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Tree saving
       StoreObjPtr<EventMetaData> evtMetaData;
       m_tree_evtNum  = evtMetaData->getEvent();
       m_tree_cid      = crystalIDs[iCharge];
       //m_tree_time     = times[iCharge];
       m_tree_time     = times_noPreviousCalibrations[iCharge];
       m_crystalCrate  = crateIDs[iCharge];
       m_runNum        = evtMetaData->getRun();
       m_tree_en       = crystalEnergies[iCharge];  // for studies of ts as a function of energy
       m_tree_E1Etot   = crystalEnergies[iCharge] / trkEClustLab[iCharge];
       m_tree_E1E2     = crystalEnergies[iCharge] / crystalEnergies2[iCharge];
       m_tree_E1p      = crystalEnergies[iCharge] / trkpLab[iCharge];
       m_tree_timetsPreviousTimeCalibs = times_noPreviousCalibrations[iCharge] - ts_prevCalib[iCharge] - tcrate_prevCalib[iCharge];
       m_tree_t0       = evt_t0;
       m_tree_t0_ECLclosestCDC   = evt_t0_ECL_closestCDC;
       m_tree_t0_ECL_minChi2   = evt_t0_ECL_minChi2;
       m_tree_tClust = clusterTime[iCharge];
  
       m_massInvTracks = invMassTrk;   // This is probably already set but I'll set it again anyways just so that it is clear
  
       if (m_saveTree) {
         m_dbgTree_allCuts->Fill();
       }
       //  <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Tree saving
     }
   }
  
  
  
  
   if (crystalCutsPassed[0]       &&  crystalCutsPassed[1]        &&
       numClustersPerTrack[0] == 1  &&  numClustersPerTrack[1] == 1) {
     m_tree_enNeg = trkEClustLab[0];
     m_tree_enPlus = trkEClustLab[1];
     m_tree_tClustNeg = clusterTime[0];
     m_tree_tClustPos = clusterTime[1];
     m_tree_maxEcrystPosClust = times_noPreviousCalibrations[0] - ts_prevCalib[0] - tcrate_prevCalib[0];
     m_tree_maxEcrystNegClust = times_noPreviousCalibrations[1] - ts_prevCalib[1] - tcrate_prevCalib[1];
     m_tree_t0       = evt_t0;
     m_tree_t0_ECLclosestCDC   = evt_t0_ECL_closestCDC;
     m_tree_t0_ECL_minChi2   = evt_t0_ECL_minChi2;
  
     //  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Tree saving
     if (m_saveTree) {
       m_dbgTree_evt_allCuts->Fill();
     }
     //  <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Tree saving
   }
  
  
   B2DEBUG(30, "m_tree_maxEcrystPosClust + evt_t0 = " << m_tree_maxEcrystPosClust + evt_t0);
   B2DEBUG(30, "m_tree_maxEcrystNegClust + evt_t0 = " << m_tree_maxEcrystNegClust + evt_t0);
   B2DEBUG(30, "CDC evt_t0 = " << evt_t0);
   B2DEBUG(30, "ECL min chi2 even t0, m_tree_t0_ECL_minChi2 = " << m_tree_t0_ECL_minChi2);
  
  
  
   for (long unsigned int digit_i = 0; digit_i < time_ECLCaldigits_bothClusters.size(); digit_i++) {
     StoreObjPtr<EventMetaData> evtMetaData;
     m_runNum        = evtMetaData->getRun();
     m_tree_evtNum  = evtMetaData->getEvent();
     m_tree_ECLCalDigitTime  = time_ECLCaldigits_bothClusters[digit_i];
     m_tree_ECLCalDigitE = E_ECLCaldigits_bothClusters[digit_i];
     m_tree_ECLDigitAmplitude = amp_ECLDigits_bothClusters[digit_i];
     m_tree_t0       = evt_t0;
     m_tree_t0_ECLclosestCDC   = evt_t0_ECL_closestCDC;
     m_tree_t0_ECL_minChi2   = evt_t0_ECL_minChi2;
     m_tree_timetsPreviousTimeCalibs = times_noPreviousCalibrations[chargeID_ECLCaldigits_bothClusters[digit_i]] -
                                       ts_prevCalib[chargeID_ECLCaldigits_bothClusters[digit_i]] -
                                       tcrate_prevCalib[chargeID_ECLCaldigits_bothClusters[digit_i]];
     m_tree_cid = cid_ECLCaldigits_bothClusters[digit_i];
  
     //  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Tree saving
     if (m_saveTree) {
       m_dbgTree_crys_allCuts->Fill();
     }
     //  <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Tree saving
  
   }
  
  
   B2DEBUG(30, "This was for event number = " << m_tree_evtNum);
  
 }

◆ 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 441 of file Module.h.

◆ 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 422 of file Module.h.

◆ 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 98 of file Module.cc.

◆ 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 135 of file Module.cc.

◆ 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 115 of file Module.cc.

◆ getFileNames()

virtual std::vector<std::string> getFileNames ( __attribute__((unused)) 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.

Definition at line 136 of file Module.h.

◆ 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 189 of file Module.h.

◆ 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 281 of file Module.cc.

◆ getPreScaleChoice()

bool getPreScaleChoice ( )

inlineprivateinherited

I'm a little worried about floating point precision when comparing to 0.0 and 1.0 as special values.

But since a user will have set them (or left them as default) as exactly equal to 0.0 or 1.0 rather than calculating them in almost every case, I think we can assume that the equalities hold.

Definition at line 134 of file CalibrationCollectorModule.h.

◆ 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 383 of file Module.h.

◆ getType()

const std::string & getType ( ) const

inherited

Returns the type of the module (i.e.

class name minus 'Module')

Definition at line 43 of file Module.cc.

◆ 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 162 of file Module.cc.

◆ 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 87 of file Module.cc.

◆ 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 92 of file Module.cc.

◆ 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://confluence.desy.de/display/BI/Software+ModCondTut 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 81 of file Module.cc.

◆ setDescription()

void setDescription ( const std::string & description )

protectedinherited

Sets the description of the module.

Parameters

description A description of the module.

Definition at line 216 of file Module.cc.

◆ 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 75 of file Module.cc.

◆ 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 216 of file Module.h.

◆ 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 236 of file Module.cc.

◆ 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 251 of file Module.cc.

◆ setPropertyFlags()

void setPropertyFlags ( unsigned int propertyFlags )

inherited

Sets the flags for the module properties.

Parameters

propertyFlags bitwise OR of EModulePropFlags

Definition at line 210 of file Module.cc.

◆ 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 229 of file Module.cc.

◆ 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 222 of file Module.cc.

◆ 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 50 of file Module.cc.

Member Data Documentation

◆ m_CrateTimeDB

DBObjPtr<ECLCrystalCalib> m_CrateTimeDB

private

Time offset from crate time calibration (also this calibration) from database.

database object

Definition at line 99 of file ECLBhabhaTCollectorModule.h.

◆ m_ECLTimeUtil

std::unique_ptr< Belle2::ECL::ECLTimingUtilities > m_ECLTimeUtil

private

Initial value:

=

std::make_unique<Belle2::ECL::ECLTimingUtilities>()

ECL timing tools.

Definition at line 209 of file ECLBhabhaTCollectorModule.h.

◆ m_ElectronicsDB

DBObjPtr<ECLCrystalCalib> m_ElectronicsDB

private

electronics amplitude calibration from database Scale amplitudefor each crystal and for dead pre-amps

database object

Definition at line 82 of file ECLBhabhaTCollectorModule.h.

◆ m_ElectronicsTimeDB

DBObjPtr<ECLCrystalCalib> m_ElectronicsTimeDB

private

Time offset from electronics calibration from database.

database object

Definition at line 86 of file ECLBhabhaTCollectorModule.h.

◆ m_eventT0

StoreObjPtr<EventT0> m_eventT0

private

StoreObjPtr for T0.

The event t0 class has an overall event t0

Definition at line 78 of file ECLBhabhaTCollectorModule.h.

◆ m_FlightTimeDB

DBObjPtr<ECLCrystalCalib> m_FlightTimeDB

private

Time offset from flight time b/w IP and crystal from database.

database object

Definition at line 90 of file ECLBhabhaTCollectorModule.h.

◆ m_PreviousCrystalTimeDB

DBObjPtr<ECLCrystalCalib> m_PreviousCrystalTimeDB

private

Time offset from previous crystal time calibration (this calibration) from database.

database object

Definition at line 94 of file ECLBhabhaTCollectorModule.h.

◆ m_tree_enPlus

double m_tree_enPlus

private

crystal energy, only for the crystals that meet all the selection criteria for debug TTree output

Energy of cluster associated to positively charged track, GeV for debug TTree output

Definition at line 152 of file ECLBhabhaTCollectorModule.h.

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

ecl/modules/eclBhabhaTCollector/include/ECLBhabhaTCollectorModule.h
ecl/modules/eclBhabhaTCollector/src/ECLBhabhaTCollectorModule.cc

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ EModulePropFlags

Member Function Documentation

◆ beginRun()

◆ clone()

◆ collect()

◆ def_endRun()

◆ def_initialize()

◆ evalCondition()

◆ getAfterConditionPath()

◆ getConditionPath()

◆ getFileNames()

◆ getName()

◆ getParamInfoListPython()

◆ getPreScaleChoice()

◆ getReturnValue()

◆ getType()

◆ hasProperties()

◆ if_false()

◆ if_true()

◆ if_value()

◆ setDescription()

◆ setLogInfo()

◆ setName()

◆ setParamPython()

◆ setParamPythonDict()

◆ setPropertyFlags()

◆ setReturnValue() [1/2]

◆ setReturnValue() [2/2]

◆ setType()

Member Data Documentation

◆ m_CrateTimeDB

◆ m_ECLTimeUtil

◆ m_ElectronicsDB

◆ m_ElectronicsTimeDB

◆ m_eventT0

◆ m_FlightTimeDB

◆ m_PreviousCrystalTimeDB

◆ m_tree_enPlus