trg modules

Classes
class	CDCTriggerNeuroDQMModule
	CDC Trigger DQM Module. More...
class	CDCTrigger2DFitterModule
	Module for the 2D Fitter of the CDC trigger. More...
class	CDCTrigger3DFitterModule
	Module for the 3D Fitter of the CDC trigger. More...
class	CDCTriggerHoughCand
	Hough Candidates class. More...
class	CDCTriggerHoughETFModule
class	CDCTrigger2DFinderModule
class	CDCTriggerMCMatcherModule
	A module to match CDCTriggerTracks to MCParticles. More...
class	CDCTriggerRecoHitMatcherModule
	A module to match CDCTriggerTracks to RecoTracks. More...
class	CDCTriggerRecoMatcherModule
	A module to match CDCTriggerTracks to RecoTracks. More...
class	CDCTriggerNDFinderModule
	CDC Trigger NDFinder Module. More...
class	CDCTriggerNeuroDataModule
class	CDCTriggerNeuroIDHistModule
	Description. More...
class	CDCTriggerNeuroModule
	The neural network module of the CDC trigger. More...
class	CDCTriggerETFModule
	Module for the Event Time Finder of the CDC trigger. More...
class	CDCTriggerTrackCombinerModule
	Module to combine the information from the various track trigger stages. More...
class	CDCTriggerTSFFirmwareModule
	This class is the interface between TSim/basf2 TSF module and the firmware simulation core of XSim/ISim. More...
class	CDCTriggerTSFModule
	Module for the Track Segment Finder of the CDC trigger. More...
class	TRGCDCModule
	A module to simulate the CDC trigger. More...
class	TRGCDCT2DDQMModule
	DQM module of TRGCDCT2D. More...
class	TRGCDCT3DDQMModule
	DQM module for TRGCDCT3D. More...
class	TRGCDCTSFDQMModule
	DQM module of TRGCDCTSF. More...
class	TRGCDCTSStreamModule
	A module to process CDC trigger data. More...
struct	SubTrigger
	enum class SubTriggerType : unsigned char {Merger, TSF, T2D, T3D, Neuro, ETF}; More...
class	CDCTriggerUnpackerModule
	Unpack the trigger data recorded in B2L. More...
struct	Merger
	unpacker for the merger reader (TSF which reads the merger output) More...
struct	Tracker2D
	unpacker for the 2D tracker More...
struct	Neuro
	unpacker for the Neuro More...
class	MCMatcherTRGECLModule
	Class to represent the hit of one cell. More...
class	TRGECLFAMModule
	A module of FAM. More...
class	TRGECLModule
	A module of ETM. More...
class	TRGECLBGTCHitModule
class	TRGECLDQMModule
class	TRGECLEventTimingDQMModule
	This module is for ecl trigger DQM related to L1 event timing. More...
class	TRGECLQAMModule
class	TRGECLRawdataAnalysisModule
	A module of ETM. More...
class	TRGECLTimingCalModule
	Class ECL Trigger Timiing Calibration Module. More...
class	TRGECLUnpackerModule
	A module of TRG ECL Unpacker. More...
class	TRGGDLCosmicRunModule
	Module that returns true if the trigger condition for the 2017 cosmic runs is fulfilled. More...
class	TRGGDLModule
	A module to simulate the Global Decision Logic. More...
class	TRGEFFDQMModule
	Make summary of data quality from reconstruction. More...
class	TRGGDLDQMModule
class	TRGGDLDSTModule
class	TRGGDLSummaryModule
class	TRGGRLMatchModule
	Match between CDC trigger track and ECL trigger cluster. More...
class	TRGGRLModule
	A module to simulate the Global Decision Logic. More...
class	TRGGRLProjectsModule
	A module to simulate the Global Decision Logic. More...
class	TRGGRLDQMModule
class	GRLNeuroModule
	The module for application of the neural networks of the GRL. More...
class	GRLNeuroTrainerModule
	The trainer module for the neural networks of the CDC trigger. More...
class	klmtrg_layer_counter_t
class	klm_trig_linear_fit_t
class	KLMTriggerModule
struct	track_maker_t
struct	Linear_fit_of_Hits_t
class	TRGTOPDQMModule
class	TRGTOPTRD2TTSConverterModule
	TRG TOP Raw Digits to TimeStamps Converter. More...
class	TRGTOPUnpackerModule
	TRG TOP Unpacker. More...
class	TRGTOPUnpackerWaveformModule
	TRG TOP Unpacker for Timestamps. More...
class	TRGTOPWaveformPlotterModule
class	TRGRAWDATAModule

Typedefs
typedef std::pair< unsigned short, ROOT::Math::XYVector >	cdcPair
	Pair of <iSuperLayer, (x, y)>, for hits in conformal space.
typedef std::map< int, cdcPair >	cdcMap
	Map of <counter, cdcPair>, for hits with indices.
typedef std::pair< ROOT::Math::XYVector, ROOT::Math::XYVector >	coord2dPair
	Hough Tuples.
using	NodeList = std::vector< std::vector< int > >
	Node list.
using	MergerBus = std::array< std::bitset< mergerWidth >, nAllMergers >
	Merger data bus.
using	MergerBits = Bitstream< MergerBus >
	Merger data bus Bitstream.

Functions
bool	operator== (const CDCTriggerHoughCand &a, const CDCTriggerHoughCand &b)
	Two cells are identical if they have the same coordinates.
	REG_MODULE (CDCTriggerNeuroData)
	REG_MODULE (CDCTriggerNeuroIDHist)
	REG_MODULE (TRGCDC)
	register the basf2 module TRGCDC
	REG_MODULE (TRGCDCTSStream)
	register the basf2 TRGCDCTSStream module
	REG_MODULE (TRGECLFAM)
	Register module name.
	REG_MODULE (TRGECL)
	Register Module.
	REG_MODULE (TRGECLBGTCHit)
	Register module name.
	REG_MODULE (TRGECLQAM)
	Register module name.
	REG_MODULE (TRGECLRawdataAnalysis)
	Register Module.
	REG_MODULE (TRGGDL)
	REG_MODULE (TRGGRL)
	Register module for TRGGRL.
int	countBits (uint64_t n)
template<typename AXIS_NAME_T , typename CONTAINER_T >
uint64_t	to_bit_mask (const CONTAINER_T &container)
void	bitshift (int64_t &out, int64_t &shift, int64_t in, int64_t deltashift)
template<typename T >
constexpr int64_t	get_index (const T &e, int track_id=0)
	CDCTriggerNeuroDataModule ()
	Constructor, for setting module description and parameters.
virtual void	initialize () override
	Initialize the module.
virtual void	event () override
	Called once for each event.
virtual void	terminate () override
	Do the training for all sectors.
	CDCTriggerNeuroIDHistModule ()
	Constructor, for setting module description and parameters.
virtual void	initialize () override
	Initialize the module.
virtual void	event () override
	Called once for each event.
virtual void	terminate () override
	Do the training for all sectors.
std::string	version (void) const
	returns version of TRGCDCModule.
	TRGCDCModule ()
	Constructor.
virtual	~TRGCDCModule ()
	Destructor.
virtual void	initialize () override
	Initilizes TRGCDCModule.
virtual void	beginRun () override
	Called when new run started.
virtual void	event () override
	Called event by event.
virtual void	endRun () override
	Called when run ended.
virtual void	terminate () override
	Called when processing ended.
std::string	version (void) const
	returns version of TRGCDCTSStreamModule.
	TRGCDCTSStreamModule ()
	Constructor.
virtual	~TRGCDCTSStreamModule ()
	Destructor.
virtual void	initialize () override
	Initilizes TRGCDCTSStreamModule.
virtual void	beginRun () override
	Called when new run started.
virtual void	event () override
	Called event by event.
virtual void	endRun () override
	Called when run ended.
virtual void	terminate () override
	Called when processing ended.
std::string	version (void) const
	returns version of TRGECLFAMModule.
	TRGECLFAMModule ()
	Constructor.
virtual	~TRGECLFAMModule ()
	Destructor.
virtual void	initialize () override
	Initilizes TRGECLFAMModule.
virtual void	beginRun () override
	Called when new run started.
virtual void	event () override
	Called event by event.
virtual void	endRun () override
	Called when run ended.
virtual void	terminate () override
	Called when processing ended.
std::string	version (void) const
	returns version of TRGECLModule.
	TRGECLModule ()
	Constructor.
virtual	~TRGECLModule ()
	Destructor
virtual void	initialize () override
	Initilizes TRGECLModule.
virtual void	beginRun () override
	Called when new run started.
virtual void	event () override
	Called event by event.
virtual void	endRun () override
	Called when run ended.
virtual void	terminate () override
	Called when processing ended.
std::string	version (void) const
	returns version of TRGECLBGTCHitModule.
	TRGECLBGTCHitModule ()
	Constructor.
virtual	~TRGECLBGTCHitModule ()
	Destructor.
virtual void	initialize () override
	Initilizes TRGECLBGTCHitModule.
virtual void	beginRun () override
	Called when new run started.
virtual void	event () override
	Called event by event.
virtual void	endRun () override
	Called when run ended.
virtual void	terminate () override
	Called when processing ended.
void	genSimulationObj ()
	Parameters.
void	genRandomTrgObj ()
std::string	version () const
	version
	TRGECLQAMModule ()
	Costructor.
virtual	~TRGECLQAMModule ()
	Destrunctor.
virtual void	initialize () override
	initialize
virtual void	beginRun () override
	begin Run
virtual void	event () override
	Event.
virtual void	endRun () override
	End Run.
virtual void	terminate () override
	terminate
std::string	version (void) const
	returns version of TRGECLModule.
	TRGECLRawdataAnalysisModule ()
	Constructor.
virtual	~TRGECLRawdataAnalysisModule ()
	Destructor
virtual void	initialize () override
	Initilizes TRGECLModule.
virtual void	beginRun () override
	Called when new run started.
virtual void	event () override
	Called event by event.
virtual void	endRun () override
	Called when run ended.
virtual void	terminate () override
	Called when processing ended.
std::string	version (void) const
	returns version of TRGGDLModule.
	TRGGDLModule ()
	Constructor.
virtual	~TRGGDLModule ()
	Destructor.
virtual void	defineHisto () override
	Define Histogram.
virtual void	initialize () override
	Initilizes TRGGDLModule.
virtual void	beginRun () override
	Called when new run started.
virtual void	event () override
	Called event by event.
virtual void	endRun () override
	Called when run ended.
virtual void	terminate () override
	Called when processing ended.
std::string	version (void) const
	returns version of TRGGRLModule.
	TRGGRLModule ()
	Constructor.
virtual	~TRGGRLModule ()
	Destructor.
virtual void	initialize () override
	Initilizes TRGGRLModule.
virtual void	beginRun () override
	Called when new run started.
virtual void	event () override
	Called event by event.
virtual void	endRun () override
	Called when run ended.
virtual void	terminate () override
	Called when processing ended.
void	clear_layersUsed ()
void	add_layersUsed (int layersUsed)
void	set_NLayerTrigger (int NlayerTrigger)
void	run (const KLM_TRG_definitions::KLM_Digit_compact_ts &hits)
int	get_n_sections_trig (int subdetector)
int	get_triggermask (int subdetector, int section)
int	get_BKLM_back_to_back_flag (int subdetector)
void	clear_geometry ()
void	add_geometry (const KLM_TRG_definitions::KLM_geo_fit_t &geometry)
void	run (const KLM_TRG_definitions::KLM_Digit_compact_ts &hits)
const KLM_TRG_definitions::KLM_trig_linear_fits &	get_result () const
int	get_triggermask (int subdetector, int section)
int	get_triggermask_or (int subdetector, int section)
void	set_y_cutoff (int cutoff)
void	set_intercept_cutoff (int cutoff)

Variables
static constexpr int	mergerWidth = 256
	Merger data width.
static constexpr int	nAllMergers = 146
	Number of Mergers.
static constexpr int	wordWidth = 32
	width of a single word in the raw int buffer
static constexpr int	nFinesse = 48
	Number of FINESSE in a PCIe40.
constexpr Linear_fit_of_Hits_t	Linear_fit_of_Hits
static const int	intNaN = std::numeric_limits<int>::quiet_NaN()
static TRGGDL *	_gdl = 0
	A pointer to a TRGGDL;.
static TRGGRL *	_grl = 0
	A pointer to a TRGGRL;.

Detailed Description

Typedef Documentation

cdcMap

typedef std::map< int, cdcPair > cdcMap

Map of <counter, cdcPair>, for hits with indices.

Definition at line 37 of file CDCTriggerHoughETFModule.h.

cdcPair

typedef std::pair< unsigned short, ROOT::Math::XYVector > cdcPair

Pair of <iSuperLayer, (x, y)>, for hits in conformal space.

Definition at line 35 of file CDCTriggerHoughETFModule.h.

coord2dPair

typedef std::pair< ROOT::Math::XYVector, ROOT::Math::XYVector > coord2dPair

Hough Tuples.

Definition at line 39 of file CDCTriggerHoughETFModule.h.

MergerBits

using MergerBits = Bitstream<MergerBus>

Merger data bus Bitstream.

Definition at line 55 of file CDCTriggerUnpackerModule.h.

MergerBus

using MergerBus = std::array<std::bitset<mergerWidth>, nAllMergers>

Merger data bus.

Definition at line 53 of file CDCTriggerUnpackerModule.h.

NodeList

using NodeList = std::vector<std::vector<int> >

Node list.

Definition at line 42 of file CDCTriggerUnpackerModule.h.

Function Documentation

add_geometry()

void add_geometry

(

const KLM_TRG_definitions::KLM_geo_fit_t &

geometry

)

Definition at line 154 of file klm_trig_linear_fit.cc.

  {
 
    m_KLMgeomap[get_index(geometry)] = geo_KLM_t{
      geometry.slopeX,
      geometry.offsetX,
      geometry.slopeY,
      geometry.offsetY
    };
 
 
  }

add_layersUsed()

void add_layersUsed

(

int

layersUsed

)

Definition at line 44 of file KLM_Trig.cc.

  {
    m_layersUsed.emplace_back(layersUsed);
  }

beginRun() [1/9]

void beginRun ( void  )

overridevirtual

Called when new run started.

Reimplemented from Module.

Definition at line 353 of file TRGCDCModule.cc.

  {
    if (TRGDebug::level())
      cout << "TRGCDCModule ... beginRun called " << endl;
  }

beginRun() [2/9]

void beginRun ( void  )

overridevirtual

Called when new run started.

Reimplemented from Module.

Definition at line 102 of file TRGCDCTSStreamModule.cc.

  {
 
    _cdc = TRGCDC::getTRGCDC();
 
    //...Super layer loop...
    for (unsigned l = 0; l < _cdc->nSegmentLayers(); l++) {
      const Belle2::TRGCDCLayer* lyr = _cdc->segmentLayer(l);
      const unsigned nWires = lyr->nCells();
 
      //...Clear old pointers...
      _wires[l].clear();
 
      //...TS loop...
      for (unsigned i = 0; i < nWires; i++) {
        const TCSegment& s = (TCSegment&) * (* lyr)[i];
        _wires[l].push_back(s.wires()[5]);
      }
    }
 
    if (_out) {
      unsigned val = TRGBSRecord_BeginRun;
      _out->write((char*) & val, 4);
      val = 0;
      _out->write((char*) & val, 4);
    }
 
    if (TRGDebug::level())
      cout << "TRGCDCTSStreamModule ... beginRun called. TRGCDC version="
           << _cdc->version() << endl;
  }

beginRun() [3/9]

void beginRun ( void  )

overridevirtual

Called when new run started.

Reimplemented from Module.

Definition at line 131 of file TRGECLFAMModule.cc.

  {
    if (_ConditionDB == 0) {
      Threshold.resize(576, _threshold);
    } else if (_ConditionDB == 1) { //Use global tag
      Threshold.resize(576, 0);
      for (const auto& para : m_FAMPara) {
        Threshold[para.getTCId() - 1] = (int)((para.getThreshold()) * (para.getConversionFactor()));
      }
    }
 
    B2DEBUG(200, "TRGECLFAMModule ... beginRun called ");
 
  }

beginRun() [4/9]

void beginRun ( void  )

overridevirtual

Called when new run started.

Reimplemented from Module.

Definition at line 255 of file TRGECLModule.cc.

  {
    if (_ConditionDB == 1) {
      for (const auto& para : m_ETMPara) {
        _TotalEnergy = {(double)para.getELow(),
                        (double)para.getELow(),
                        (double)para.getELow()
                       };
        for (int index = 0; index < 14; index ++) {
          _2DBhabhaThresholdFWD.push_back((double)para.get2DBhabhaFWD(index));
          _2DBhabhaThresholdBWD.push_back((double)para.get2DBhabhaBWD(index));
        }
        _3DBhabhaVetoThreshold = {
          (double)para.get3DBhabhaVetoThreshold(0),
          (double)para.get3DBhabhaVetoThreshold(1)
        }; //  /100 MeV
        _3DBhabhaSelectionThreshold = {
          (double)para.get3DBhabhaSelectionThreshold(0),
          (double)para.get3DBhabhaSelectionThreshold(1)
        }; //  /100 MeV
        _3DBhabhaVetoAngle = {
          (double)para.get3DBhabhaVetoAngle(0),
          (double)para.get3DBhabhaVetoAngle(1),
          (double)para.get3DBhabhaVetoAngle(2),
          (double)para.get3DBhabhaVetoAngle(3)
        };
        _3DBhabhaSelectionAngle = {
          (double)para.get3DBhabhaSelectionAngle(0),
          (double)para.get3DBhabhaSelectionAngle(1),
          (double)para.get3DBhabhaSelectionAngle(2),
          (double)para.get3DBhabhaSelectionAngle(3)
        };
        _mumuThreshold =
          (double)(para.getmumuThreshold());
        _mumuAngle = {
          (double)para.getmumuAngle(0),
          (double)para.getmumuAngle(1),
          (double)para.getmumuAngle(2)
        };
        _LowMultiThreshold = {
          (double)para.getLowMultiThreshold(0),
          (double)para.getLowMultiThreshold(1),
          (double)para.getLowMultiThreshold(2),
          (double)para.getLowMultiThreshold(3)
        }; //  /100 MeV
      }
    }
 
    B2INFO("[TRGECLModule] 3DBhabhaVetoInTrackThetaRegion (low : high) = ("
           << m_3DBhabhaVetoInTrackThetaRegion[0] << " : "
           << m_3DBhabhaVetoInTrackThetaRegion[1] << ")");
 
    B2INFO("[TRGECLModule] EventTimingQualityThresholds (low : high)= ("
           << m_EventTimingQualityThresholds[0] << " : "
           << m_EventTimingQualityThresholds[1] << ") (GeV)");
 
    B2INFO("[TRGECLModule] Taub2b two CL Angle cut in CM (degree) (dphi L, H, theta sum L, H ) = ("
           << m_taub2bAngleCut[0] << ","
           << m_taub2bAngleCut[1] << ","
           << m_taub2bAngleCut[2] << ","
           << m_taub2bAngleCut[3] << ")");
    B2INFO("[TRGECLModule] Taub2b Total Energy Cut in lab. (GeV) = "
           << m_taub2bEtotCut);
    B2INFO("[TRGECLModule] Taub2b Cluster Energy Cut in lab. (GeV) : (E(CL1), E(CL2)) = ("
           << m_taub2bClusterECut1 << ", "
           << m_taub2bClusterECut2 << ")");
    B2INFO("[TRGECLModule] Bhabha Add Angle cut in CM (degree) (dphi L, H, theta sum L, H) = ("
           << m_3DBhabhaAddAngleCut[0] << ","
           << m_3DBhabhaAddAngleCut[1] << ","
           << m_3DBhabhaAddAngleCut[2] << ","
           << m_3DBhabhaAddAngleCut[3] << ")");
 
    B2INFO("[TRGECLModule] Taub2b2 two CL Angle cut in CM (degree) (dphi L, H, theta sum L, H) = ("
           << m_taub2b2AngleCut[0] << ","
           << m_taub2b2AngleCut[1] << ","
           << m_taub2b2AngleCut[2] << ","
           << m_taub2b2AngleCut[3] << ")");
    B2INFO("[TRGECLModule] Taub2b2 Total Energy Cut in lab. (GeV) = "
           << m_taub2b2EtotCut);
    B2INFO("[TRGECLModule] Taub2b2 Energy Cut in lab. for cluster in endcap (GeV) = "
           << m_taub2b2CLEEndcapCut);
    B2INFO("[TRGECLModule] Taub2b2 Cluster energy cut in lab. (GeV) = "
           << m_taub2b2CLECut);
 
    B2INFO("[TRGECLModule] Taub2b3 two CL Angle cut in CM (degree) (dphi L, H, theta sum L, H) = ("
           << m_taub2b3AngleCut[0] << ","
           << m_taub2b3AngleCut[1] << ","
           << m_taub2b3AngleCut[2] << ","
           << m_taub2b3AngleCut[3] << ")");
    B2INFO("[TRGECLModule] Taub2b3 Total Energy Cut in lab. (GeV) = "
           << m_taub2b3EtotCut);
    B2INFO("[TRGECLModule] Taub2b3 Cluster energy cut in lab. for one of b2b clusters (GeV) = "
           << m_taub2b3CLEb2bCut);
    B2INFO("[TRGECLModule] Taub2b3 Cluster energy low  cut in lab. for all clusters (GeV) = "
           << m_taub2b3CLELowCut);
    B2INFO("[TRGECLModule] Taub2b3 Cluster energy high cut in lab. for all clusters(GeV) = "
           << m_taub2b3CLEHighCut);
 
    if (TRGDebug::level()) {
      std::cout << "TRGECLModule ... beginRun called " << std::endl;
    }
    //  _ecl = TrgEcl::getTrgEcl();
 
  }

beginRun() [5/9]

void beginRun ( void  )

overridevirtual

Called when new run started.

Reimplemented from Module.

Definition at line 103 of file TRGECLBGTCHitModule.cc.

  {
    B2DEBUG(200, "TRGECLBGTCHitModule ... beginRun called ");
  }

beginRun() [6/9]

void beginRun ( void  )

overridevirtual

begin Run

Reimplemented from Module.

Definition at line 132 of file TRGECLQAMModule.cc.

  {
 
    B2DEBUG(200, "TRGECLQAMModule ... beginRun called ");
 
  }

beginRun() [7/9]

void beginRun ( void  )

overridevirtual

Called when new run started.

Reimplemented from Module.

Definition at line 97 of file TRGECLRawdataAnalysisModule.cc.

  {
 
    if (TRGDebug::level()) {
      std::cout << "TRGECLRawdataAnalysisModule ... beginRun called " << std::endl;
    }
    //  _ecl = TrgEcl::getTrgEcl();
 
  }

beginRun() [8/9]

void beginRun ( void  )

overridevirtual

Called when new run started.

Reimplemented from HistoModule.

Definition at line 147 of file TRGGDLModule.cc.

  {
 
    //...GDL config. name...
    string cfn = _configFilename;
 
    //...GDL...
    if (_gdl == 0) {
      _gdl = TRGGDL::getTRGGDL(cfn,
                               _simulationMode,
                               _fastSimulationMode,
                               _firmwareSimulationMode,
                               _Phase,
                               _algFromDB,
                               _algFilePath,
                               _debugLevel,
                               _timquality_threshold_sfin,
                               _timquality_threshold_fine);
    } else if (cfn != _gdl->configFile()) {
      _gdl = TRGGDL::getTRGGDL(cfn,
                               _simulationMode,
                               _fastSimulationMode,
                               _firmwareSimulationMode,
                               _Phase,
                               _algFromDB,
                               _algFilePath,
                               _debugLevel,
                               _timquality_threshold_sfin,
                               _timquality_threshold_fine);
    }
    if (_debugLevel > 9) printf("TRGGDLModule::beginRun() ends.\n");
 
    _gdl->checkDatabase();
 
    B2DEBUG(100, "TRGGDLModule ... beginRun called  configFile = " << cfn);
  }

beginRun() [9/9]

void beginRun ( void  )

overridevirtual

Called when new run started.

Reimplemented from Module.

Definition at line 86 of file TRGGRLModule.cc.

  {
 
    //...GDL config. name...
    string cfn = _configFilename;
 
    //...GRL...
    if (_grl == 0) {
      _grl = TRGGRL::getTRGGRL(cfn,
                               _simulationMode,
                               _fastSimulationMode,
                               _firmwareSimulationMode);
    } else if (cfn != _grl->configFile()) {
      _grl = TRGGRL::getTRGGRL(cfn,
                               _simulationMode,
                               _fastSimulationMode,
                               _firmwareSimulationMode);
    }
 
    B2DEBUG(300, "TRGGDLModule ... beginRun called configFile = " << cfn);
 
  }

bitshift()

void bitshift	(	int64_t &	out,
		int64_t &	shift,
		int64_t	in,
		int64_t	deltashift
	)

Definition at line 31 of file klm_trig_linear_fit.cc.

  {
    int64_t mask = 0xFFFFFF;
    if (((mask << deltashift) & in) > 0) {
      throw std::runtime_error("bitshift error");
    }
 
    out |= (in << shift);
    shift += deltashift;
  }

CDCTriggerNeuroDataModule()

CDCTriggerNeuroDataModule

(

)

Constructor, for setting module description and parameters.

Definition at line 42 of file CDCTriggerNeuroDataModule.cc.

                                                       : Module()
  {
    setDescription(
      "This module takes 2dtracks, track segments, and targettracks (either recotracks or mcparticles) as input and generates training data for the neurotrigger in a tab separated, gzip compressed file."
    );
    // parameters for saving / loading
    addParam("hitCollectionName", m_hitCollectionName,
             "Name of the input StoreArray of CDCTriggerSegmentHits. The Axials need to have a relation to inputtracks",
             std::string(""));
    addParam("inputCollectionName", m_inputCollectionName,
             "Name of the StoreArray holding the 2D input tracks.",
             std::string("CDCTriggerNNInput2DTracks"));
    addParam("targetCollectionName", m_targetCollectionName,
             "Name of the MCParticle/RecoTrack collection used as target values.",
             std::string("RecoTracks"));
    addParam("trainOnRecoTracks", m_trainOnRecoTracks,
             "If true, use RecoTracks as targets instead of MCParticles.",
             true);
    addParam("EventTimeName", m_EventTimeName,
             "Name of the event time object.",
             std::string("CDCTriggerNeuroETFT0"));
    addParam("NeuroTrackInputMode", m_neuroTrackInputMode,
             "When using real tracks, use neurotracks instead of 2dtracks as input to the neurotrigger. this is important to get the relations right.",
             false);
    addParam("singleUse", m_singleUse,
             "Only use a track for a single expert", true);
    addParam("configFileName", m_configFileName,
             "Name of the configuration file. This File should be created by the CDCTriggerIDHistModule and will be extended in this module",
             std::string(""));
    addParam("writeconfigFileName", m_writeconfigFileName,
             "Name of the configuration file, which will be written. If left blank, the same file as the input configuration file is used and it will be overwritten.",
             std::string(""));
    addParam("gzipFilename", m_filename,
             "Name of the gzip file, where the training samples will be saved.",
             std::string("out.gz"));
 
 
 
 
  }

CDCTriggerNeuroIDHistModule()

CDCTriggerNeuroIDHistModule

(

)

Constructor, for setting module description and parameters.

Definition at line 38 of file CDCTriggerNeuroIDHistModule.cc.

                                                           : Module()
  {
    setDescription(
      "description" //TODO
    );
    // parameters for saving / loading
    addParam("hitCollectionName", m_hitCollectionName,
             "Name of the input StoreArray of CDCTriggerSegmentHits.",
             std::string(""));
    addParam("inputCollectionName", m_inputCollectionName,
             "Name of the StoreArray holding the 2D input tracks.",
             std::string("TRGCDC2DFinderTracks"));
    addParam("trainOnRecoTracks", m_trainOnRecoTracks,
             "If true, use RecoTracks as targets instead of MCParticles.",
             false);
    addParam("targetCollectionName", m_targetCollectionName,
             "Name of the MCParticle/RecoTrack collection used as target values.",
             std::string("MCParticles"));
    addParam("writeconfigfile", m_writeconfigFileName,
             "Name of the config file, where all the parameters and the IDHist configuration is written.",
             std::string(""));
    addParam("configfile", m_configFileName,
             "Name of the config file, where all the parameters and the IDHist configuration is read in from.",
             std::string(""));
    addParam("MaxEvents", m_nPrepare,
             "amount of events used for creating the IDHist. If it is 0, "
             "all Events are used.",
             0);
 
  }

clear_geometry()

void clear_geometry

(

)

Definition at line 148 of file klm_trig_linear_fit.cc.

  {
    m_KLMgeomap.clear();
 
  }

clear_layersUsed()

void clear_layersUsed

(

)

Definition at line 39 of file KLM_Trig.cc.

  {
    m_layersUsed.clear();
  }

countBits()

int countBits ( uint64_t  n )

inline

Definition at line 21 of file bit_operations.h.

  {
    return static_cast<int>(std::bitset<64>(n).count());
  }

defineHisto()

void defineHisto ( )

overridevirtual

Define Histogram.

Reimplemented from HistoModule.

Definition at line 109 of file TRGGDLModule.cc.

  {
 
    oldDir = gDirectory;
    newDir = gDirectory;
    oldDir->mkdir("TRGGDLModule");
    newDir->cd("TRGGDLModule");
 
    h_inp  = new TH1I("hTRGGDL_inp", "input bits from TRGGDLModule", 200, 0, 200);
    h_ftd  = new TH1I("hTRGGDL_ftd", "ftdl  bits from TRGGDLModule", 200, 0, 200);
    h_psn  = new TH1I("hTRGGDL_psn", "psnm  bits from TRGGDLModule", 200, 0, 200);
 
    oldDir->cd();
 
  }

endRun() [1/9]

void endRun ( void  )

overridevirtual

Called when run ended.

Reimplemented from Module.

Definition at line 391 of file TRGCDCModule.cc.

  {
    if (TRGDebug::level())
      cout << "TRGCDCModule ... endRun called " << endl;
  }

endRun() [2/9]

void endRun ( void  )

overridevirtual

Called when run ended.

Reimplemented from Module.

Definition at line 203 of file TRGCDCTSStreamModule.cc.

  {
    if (_out) {
      unsigned val = TRGBSRecord_EndRun;
      _out->write((char*) & val, 4);
      val = 0;
      _out->write((char*) & val, 4);
    }
 
    if (TRGDebug::level())
      cout << "TRGCDCTSStreamModule ... endRun called " << endl;
  }

endRun() [3/9]

void endRun ( void  )

overridevirtual

Called when run ended.

Reimplemented from Module.

Definition at line 233 of file TRGECLFAMModule.cc.

  {
    B2DEBUG(200, "TRGECLFAMModule ... endRun called ");
  }

endRun() [4/9]

void endRun ( void  )

overridevirtual

Called when run ended.

Reimplemented from Module.

Definition at line 428 of file TRGECLModule.cc.

  {
    if (TRGDebug::level()) {
      std::cout << "TRGECLModule ... endRun called " << std::endl;
    }
  }

endRun() [5/9]

void endRun ( void  )

overridevirtual

Called when run ended.

Reimplemented from Module.

Definition at line 125 of file TRGECLBGTCHitModule.cc.

  {
    B2DEBUG(200, "TRGECLBGTCHitModule ... endRun called ");
  }

endRun() [6/9]

void endRun ( void  )

overridevirtual

End Run.

Reimplemented from Module.

Definition at line 201 of file TRGECLQAMModule.cc.

  {
 
    double mean_FWD = 0;
    double mean_BAR = 0;
    double mean_BWD = 0;
 
    //QAM TC Hit Map Check
    for (int TCId = 1; TCId < 577; TCId++) {
      if (TCId < 81) { //Forward Endcap
        mean_FWD += TCID[TCId - 1];
      } else if (TCId < 513) { //Barrel
        mean_BAR += TCID[TCId - 1];
      } else { //Backward Endcap
        mean_BWD += TCID[TCId - 1];
      }
    }
    mean_FWD /= 80;
    mean_BAR /= (512 - 80);
    mean_BWD /= (576 - 512);
 
    for (int TCId = 1; TCId < 577; TCId++) {
      if (TCId < 81) { //Forward Endcap
        if (TCID[TCId - 1] < mean_FWD * 0.1)m_FWD++;
      } else if (TCId < 513) { //Barrel
        if (TCID[TCId - 1] < mean_BAR * 0.1)m_BAR++;
      } else { //Backward Endcap
        if (TCID[TCId - 1] < mean_BWD * 0.1)m_BWD++;
      }
    }
 
    m_ALL = m_FWD + m_BAR + m_BWD;
    const int etot_size = etot.size();
    for (int i = 0; i < etot_size; i++) {
      h_etot->Fill(etot[i]);
    }
    //Caluate E_total peak and width
 
    RooRealVar* E_tot = new RooRealVar("E_tot", "E_tot", 0, 4000);
    RooDataHist* data_E_tot = new RooDataHist("data_E_tot", "dataset with E_tot", *E_tot, h_etot);
    RooRealVar* mean_check = new RooRealVar("mean_check", "Mean for checking", 2000, 1000, 3000);
    RooRealVar* sigma_check = new RooRealVar("sigma_check", "Sigma for checking", 100, 0, 300);
    RooGaussian* gauss_check = new RooGaussian("gauss_check", "Gaussian for checking", *E_tot, *mean_check, *sigma_check);
    gauss_check->fitTo(*data_E_tot);
    double mean_ch = mean_check->getVal();
 
    RooRealVar* mean_E_tot = new RooRealVar("mean_E_tot", "Mean of Gaussian", mean_ch, mean_ch - 200, mean_ch + 200);
    //  RooRealVar sigma_E_tot("sigma_E_tot","Width of Gaussian", 10, 0, 50);
 
    RooRealVar* sigma1_E_tot = new RooRealVar("sigma1_E_tot", "Width of Gaussian", 30, 0, 60);
    RooGaussian* gauss1_E_tot = new RooGaussian("gauss1_E_tot", "gauss(x,mean,sigma)", *E_tot, *mean_E_tot, *sigma1_E_tot);
    RooRealVar* sigma2_E_tot = new RooRealVar("sigma2_E_tot", "Width of Gaussian", 100, 30, 250);
    RooGaussian* gauss2_E_tot = new RooGaussian("gauss2_E_tot", "gauss(x,mean,sigma)", *E_tot, *mean_E_tot, *sigma2_E_tot);
 
    RooRealVar* frac_core_E_tot = new RooRealVar("frac_core", "core fraction", 0.85, 0.0, 1.0);
    RooAddModel* gaussm_E_tot = new RooAddModel("gaussm", "core+tail gauss", RooArgList(*gauss1_E_tot, *gauss2_E_tot),
                                                RooArgList(*frac_core_E_tot));
 
    //  RooGaussian gauss_E_tot("gauss_E_tot","gauss(x,mean,sigma)",E_tot,mean_E_tot,sigma_E_tot);
 
    gaussm_E_tot->fitTo(*data_E_tot, RooFit::Range(mean_ch - 300, mean_ch + 300));
    //RooPlot* xframe_E_tot = E_tot->frame();
    RooPlot* xframe_E_tot = E_tot->frame(RooFit::Title("energy_tot"));
    //data_E_tot->plotOn(xframe_E_tot);
    data_E_tot->plotOn(xframe_E_tot, RooFit::MarkerStyle(7));
    //  gauss_E_tot->plotOn(xframe_E_tot);
    gaussm_E_tot->plotOn(xframe_E_tot, RooFit::Name("Cal_energy"), RooFit::LineColor(2));
    double v_m_E_tot  = mean_E_tot->getVal();
    double e_m_E_tot  = mean_E_tot->getError();
 
    double v1_s1_E_tot = sigma1_E_tot->getVal();
    //    double e1_s1_E_tot = sigma1_E_tot->getError();
    double v2_s1_E_tot = sigma2_E_tot->getVal();
    //    double e2_s1_E_tot = sigma2_E_tot->getError();
 
    double frac_E_tot = frac_core_E_tot->getVal();
    double v_s1_E_tot = v1_s1_E_tot * (frac_E_tot) + v2_s1_E_tot * (1 - frac_E_tot);
 
 
 
    for (long unsigned int ii = 0; ii < caltime.size(); ii++) {
      h_caltime->Fill(caltime[ii]);
    }
 
 
    //Calculate Caltime peak and width
    RooRealVar T_c("T_c", "T_c", -100, 100);
    RooDataHist data_T_c("data_T_c", "dataset with T_c", T_c, h_caltime);
 
    RooRealVar mean_T_c("mean_T_c", "Mean of Gaussian", 0, -100, 100);
    RooRealVar sigma1_T_c("sigma1_T_c", "Width of Gaussian", 10, 0, 50);
    RooGaussian gauss1_T_c("gauss1_T_c", "gauss(x,mean,sigma)", T_c, mean_T_c, sigma1_T_c);
    RooRealVar sigma2_T_c("sigma2_T_c", "Width of Gaussian", 30, 10, 100);
    RooGaussian gauss2_T_c("gauss2_T_c", "gauss(x,mean,sigma)", T_c, mean_T_c, sigma2_T_c);
 
    RooRealVar frac_core_T_c("frac_core", "core fraction", 0.85, 0.0, 1.0);
    RooAddModel gaussm_T_c("gaussm", "core+tail gauss", RooArgList(gauss1_T_c, gauss2_T_c), RooArgList(frac_core_T_c));
 
    gaussm_T_c.fitTo(data_T_c, RooFit::Range(-45, 45));
    RooPlot* xframe_T_c = T_c.frame(RooFit::Title("Caltime"));
    data_T_c.plotOn(xframe_T_c, RooFit::MarkerStyle(7));
    gaussm_T_c.plotOn(xframe_T_c, RooFit::Name("Caltime"), RooFit::LineColor(2));
    gaussm_T_c.plotOn(xframe_T_c, RooFit::Components("gauss1_T_c"), RooFit::LineColor(kGreen));
    gaussm_T_c.plotOn(xframe_T_c, RooFit::Components("gauss2_T_c"), RooFit::LineColor(kMagenta));
 
    xframe_T_c->GetXaxis()->SetTitle("Caltime[ns]");
    xframe_T_c->GetYaxis()->SetTitle("Entries");
 
    double v_m_T_c  = mean_T_c.getVal();
    double e_m_T_c  = mean_T_c.getError();
 
    double v1_s1_T_c = sigma1_T_c.getVal();
    //    double e1_s1_T_c = sigma1_T_c.getError();
    double v2_s1_T_c = sigma2_T_c.getVal();
    //    double e2_s1_T_c = sigma2_T_c.getError();
 
    double frac_T_c = frac_core_T_c.getVal();
    double v_s1_T_c = v1_s1_T_c * (frac_T_c) + v2_s1_T_c * (1 - frac_T_c);
 
 
 
    //Cluster E sum for Bhabha skim
    const int cluster_size = cluster.size();
    for (int ii = 0; ii < cluster_size; ii++) {
      h_clusterE->Fill(cluster[ii]);
    }
 
    RooRealVar x("x", "esum", 0, 4000);
    RooDataHist dh("dh", "ecltrg",   x, h_clusterE);
 
    RooRealVar mean("mean",    "mean", 2000, 1000, 2500);
    RooRealVar sigma("sigma", "sigma",  150,   0,   300);
    RooRealVar tail("tail",    "tail",  0.45,     0,     1);
    RooNovosibirsk novo("novo", "", x, mean, sigma, tail);
    novo.fitTo(dh, RooFit::Extended(0), RooFit::Range(1500, 2200));
 
 
    double clusterE_vm, clusterE_em, clusterE_vs;
 
    clusterE_vm = mean.getVal();
    clusterE_em = mean.getError();
    clusterE_vs = sigma.getVal();
 
 
    TString outputfile = m_outputfile;
    TFile file(outputfile, "RECREATE");
    TTree* tree = new TTree("tree", "tree");
    int nevent;
    double FWD;
    double BAR;
    double BWD;
    double ALL;
 
    nevent = m_nevent;
    FWD = m_FWD;
    BAR = m_BAR;
    BWD = m_BWD;
    ALL = m_ALL;
 
    m_etot_mean = v_m_E_tot;
    m_etot_error = e_m_E_tot;
    m_etot_sigma = v_s1_E_tot;
 
    m_caltime_mean = v_m_T_c;
    m_caltime_error = e_m_T_c;
    m_caltime_sigma = v_s1_T_c;
 
    tree->Branch("m_nRun",  &m_nRun,   "m_nRun/I");
    tree->Branch("m_nevent",  &nevent,   "m_nevent/I");
    tree->Branch("m_FWD",  &FWD,   "m_FWD/D");
    tree->Branch("m_BAR",  &BAR,   "m_BAR/D");
    tree->Branch("m_BWD",  &BWD,   "m_BWD/D");
    tree->Branch("m_ALL",  &ALL,   "m_ALL/D");
    tree->Branch("m_etot_mean",  &m_etot_mean,   "m_etot_mean/D");
    tree->Branch("m_etot_error",  &m_etot_error,   "m_etot_error/D");
    tree->Branch("m_etot_sigma",  &m_etot_sigma,   "m_etot_sigma/D");
    tree->Branch("m_caltime_mean",  &m_caltime_mean,   "m_caltime_mean/D");
    tree->Branch("m_caltime_error",  &m_caltime_error,   "m_caltime_error/D");
    tree->Branch("m_caltime_sigma",  &m_caltime_sigma,   "m_caltime_sigma/D");
    tree->Branch("m_clusterE_mean",  &clusterE_vm,   "m_clusterE_mean/D");
    tree->Branch("m_clusterE_error",  &clusterE_em,   "m_clusterE_error/D");
    tree->Branch("m_clusterE_sigma",  &clusterE_vs,   "m_clusterE_sigma/D");
 
    tree->Fill();
    tree->Write();
    file.Close();
 
 
 
    B2INFO("ECL Trigger QAM result "
           << LogVar("Run Number",  m_nRun)
           << LogVar("Total Event",  nevent)
           << LogVar("The # of Low Hit TC in Forward Endcap",  FWD)
           << LogVar("The # of Low Hit TC in Barrel",  BAR)
           << LogVar("The # of Low Hit TC in Backward Endcap",  BWD)
           << LogVar("Total Energy Peak",  m_etot_mean)
           << LogVar("Total Energy Peak width",  m_etot_sigma)
           << LogVar("Caltime Peak",  m_caltime_mean)
           << LogVar("Caltime Peak width",  m_caltime_sigma)
           << LogVar("Cluster Energy Peak",  clusterE_vm)
           << LogVar("Cluster Energy Peak width",  clusterE_vs));
 
 
    B2DEBUG(100, "TRGECLQAMModule ... endRun called ");
 
  }

endRun() [7/9]

void endRun ( void  )

overridevirtual

Called when run ended.

Reimplemented from Module.

Definition at line 184 of file TRGECLRawdataAnalysisModule.cc.

  {
    if (TRGDebug::level()) {
      std::cout << "TRGECLRawdataAnalysisModule ... endRun called " << std::endl;
    }
  }

endRun() [8/9]

void endRun ( void  )

overridevirtual

Called when run ended.

Reimplemented from HistoModule.

Definition at line 212 of file TRGGDLModule.cc.

  {
    B2DEBUG(200, "TRGGDLModule ... endRun called ");
  }

endRun() [9/9]

void endRun ( void  )

overridevirtual

Called when run ended.

Reimplemented from Module.

Definition at line 120 of file TRGGRLModule.cc.

  {
    B2DEBUG(300, "TRGGRLModule ... endRun called ");
  }

event() [1/11]

void event ( void  )

overridevirtual

Called once for each event.

Prepare input and target for each track and store it.

Reimplemented from Module.

Definition at line 134 of file CDCTriggerNeuroDataModule.cc.

  {
    StoreObjPtr<EventMetaData> evtmetadata;
    for (int itrack = 0; itrack < m_tracks.getEntries(); ++itrack) {
      // get related MCParticle/RecoTrack for target
      // and retrieve track parameters
      //float phi0Target = 0;
      //float invptTarget = 0;
      float thetaTarget = 0;
      float zTarget = 0;
      if (m_trainOnRecoTracks) {
        RecoTrack* recoTrack =
          m_tracks[itrack]->getRelatedTo<RecoTrack>(m_targetCollectionName);
        if (!recoTrack) {
          B2DEBUG(150, "Skipping CDCTriggerTrack without relation to RecoTrack.");
          continue;
        }
        // a RecoTrack has multiple representations for different particle hypothesis
        // -> just take the first one that does not give errors.
        const std::vector<genfit::AbsTrackRep*>& reps = recoTrack->getRepresentations();
        bool foundValidRep = false;
        for (unsigned irep = 0; irep < reps.size() && !foundValidRep; ++irep) {
          if (!recoTrack->wasFitSuccessful(reps[irep])) {
            continue;
          }
 
          // get state (position, momentum etc.) from hit closest to IP and
          // extrapolate to z-axis (may throw an exception -> continue to next representation)
          try {
            genfit::MeasuredStateOnPlane state =
              recoTrack->getMeasuredStateOnPlaneClosestTo(ROOT::Math::XYZVector(0, 0, 0), reps[irep]);
            reps[irep]->extrapolateToLine(state, TVector3(0, 0, -1000), TVector3(0, 0, 2000));
            // flip tracks if necessary, such that trigger tracks and reco tracks
            // point in the same direction
            if (state.getMom().Dot(B2Vector3D(m_tracks[itrack]->getDirection())) < 0) {
              state.setPosMom(state.getPos(), -state.getMom());
              state.setChargeSign(-state.getCharge());
            }
            // get track parameters
            //phi0Target = state.getMom().Phi();
            //invptTarget = state.getCharge() / state.getMom().Pt();
            thetaTarget = state.getMom().Theta();
            zTarget = state.getPos().Z();
          } catch (...) {
            continue;
          }
          // break loop
          foundValidRep = true;
        }
        if (!foundValidRep) {
          B2DEBUG(150, "No valid representation found for RecoTrack, skipping.");
          continue;
        }
      } else {
        MCParticle* mcTrack =
          m_tracks[itrack]->getRelatedTo<MCParticle>(m_targetCollectionName);
        if (not mcTrack) {
          B2DEBUG(150, "Skipping CDCTriggerTrack without relation to MCParticle.");
          continue;
        }
        //phi0Target = mcTrack->getMomentum().Phi();
        //invptTarget = mcTrack->getCharge() / mcTrack->getMomentum().Pt();
        thetaTarget = mcTrack->getMomentum().Theta();
        zTarget = mcTrack->getProductionVertex().Z();
      }
      m_NeuroTrigger.updateTrack(*m_tracks[itrack]);
 
      // find all matching sectors
      float phi0 = m_tracks[itrack]->getPhi0();
      float invpt = m_tracks[itrack]->getKappa(1.5);
      float theta = atan2(1., m_tracks[itrack]->getCotTheta());
      std::vector<int> sectors = m_NeuroTrigger.selectMLPsTrain(phi0, invpt, theta);
      if (sectors.size() == 0) continue;
 
      // get target values
      std::vector<float> targetRaw = {};
      if (m_neuroParameters.targetZ)
        targetRaw.push_back(zTarget);
      if (m_neuroParameters.targetTheta)
        targetRaw.push_back(thetaTarget);
      for (unsigned i = 0; i < sectors.size(); ++i) {
        int isector = sectors[i];
        std::vector<float> target = m_NeuroTrigger[isector].scaleTarget(targetRaw);
        // skip out of range targets or rescale them
        bool outOfRange = false;
        for (unsigned itarget = 0; itarget < target.size(); ++itarget) {
          if (fabs(target[itarget]) > 1.) {
            outOfRange = true;
            target[itarget] /= fabs(target[itarget]);
          }
        }
        if (!m_neuroParameters.rescaleTarget && outOfRange) {
          continue;
        }
        //
        // read out or determine event time
        m_NeuroTrigger.getEventTime(isector, *m_tracks[itrack], m_neuroParameters.et_option(), m_neuroTrackInputMode);
        // check hit pattern
        unsigned long hitPattern = m_NeuroTrigger.getCompleteHitPattern(isector, *m_tracks[itrack], m_neuroTrackInputMode); // xxxxx0xxx
        // sectorpattern holds the absolut necessary SLs for the expert
        unsigned long sectorPattern = m_NeuroTrigger[isector].getSLpattern(); // 010100010
        // sectorpatternmask holds the SLs, which are generally used to determine the right expert
        unsigned long sectorPatternMask = m_NeuroTrigger[isector].getSLpatternMask(); // 010101010
        B2DEBUG(250, "hitPattern " << hitPattern << " sectorPattern " << sectorPattern);
        // if multiple experts should be trained with one track, the necessary SLs have to be in the hitpattern
        if (!m_singleUse && sectorPattern > 0 && (sectorPattern & hitPattern) != sectorPattern) {
          B2DEBUG(250, "hitPattern not matching " << (sectorPattern & hitPattern));
          continue;
          // if one track should only be used for one expert, the absolute necessary SLs for this expert should match exactly the set within the hitpattern of the generally used SLs for determining the expert.
        } else if (m_singleUse && sectorPattern > 0 && (sectorPattern) != (hitPattern & sectorPatternMask)) {
          // 010100010 != 0x0x0x0x0
          B2DEBUG(250, "hitPattern not matching " << (sectorPatternMask & hitPattern));
          continue;
        }
        // check, if enough axials are there. first, we select the axial bits from the
        // hitpattern (341 = int('101010101',base=2)) and then check if the number of
        // ones is equal or greater than 4.
        if ((hitPattern & 341) != 341 && // this is an ugly workaround, because popcount is only
            (hitPattern & 341) != 340 && // available with c++20 and newer
            (hitPattern & 341) != 337 &&
            (hitPattern & 341) != 325 &&
            (hitPattern & 341) != 277 &&
            (hitPattern & 341) != 85) {
          B2DEBUG(250, "Not enough axial hits (<4), skipping!");
          continue;
        }
        // get training data
        std::vector<unsigned> hitIds;
        if (m_neuroTrackInputMode) {
          hitIds = m_NeuroTrigger.selectHitsHWSim(isector, *m_tracks[itrack]);
        } else {
          hitIds = m_NeuroTrigger.selectHits(isector, *m_tracks[itrack]);
        }
        // add a "sample" it is a full training set including some zeroes for the neurotrigger values (raw and scaled z and theta). Those are sometimes used in the training for reference, but cannot be added here without running *some* neuotrigger instance.
        CDCTriggerMLPData::NeuroSet<27, 2> sample(m_NeuroTrigger.getInputVector(isector, hitIds).data(), target.data(),
                                                  evtmetadata->getExperiment(), evtmetadata->getRun(), evtmetadata->getSubrun(), evtmetadata->getEvent(), itrack, i,
                                                  m_tracks.getEntries(), 0, 0, 0, 0, phi0, theta, invpt);
        //check whether we already have enough samples
        m_trainSet[isector].addSample(sample);
        if ((m_trainSet)[isector].nSamples() % 1000 == 0) {
          B2DEBUG(50, m_trainSet[isector].nSamples() << " samples for training collected for sector " << isector);
        }
        std::ofstream gzipfile4(m_filename, std::ios_base::app | std::ios_base::binary);
        boost::iostreams::filtering_ostream outStream;
        outStream.push(boost::iostreams::gzip_compressor());
        outStream.push(gzipfile4);
        outStream << sample << std::endl;
        m_trackcounter++;
      }
 
    }
  }

event() [2/11]

void event ( void  )

overridevirtual

Called once for each event.

Prepare input and target for each track and store it.

Reimplemented from Module.

Definition at line 126 of file CDCTriggerNeuroIDHistModule.cc.

  {
    StoreObjPtr<EventMetaData> evtmetadata;
    for (int itrack = 0; itrack < m_tracks.getEntries(); ++itrack) {
      // get related MCParticle/RecoTrack for target
      // and retrieve track parameters
 
      std::vector<float> targetvector = NeuroTrainer::getTrainTargets(m_trainOnRecoTracks, m_tracks[itrack], m_targetCollectionName);
      if (targetvector[4] == 0) {
        continue;
      } // no valid representation found
      // float phi0Target = targetvector[0];  // currently unused
      // float invptTarget = targetvector[1]; // currently unused
      float thetaTarget = targetvector[2];
      float zTarget = targetvector[3];
 
 
      // update 2D track variables
      m_NeuroTrigger.updateTrack(*m_tracks[itrack]);
 
      // find all matching sectors
      float phi0 = m_tracks[itrack]->getPhi0();
      float invpt = m_tracks[itrack]->getKappa(1.5);
      float theta = atan2(1., m_tracks[itrack]->getCotTheta());
      std::vector<int> sectors = m_NeuroTrigger.selectMLPsTrain(phi0, invpt, theta);
      if (sectors.size() == 0) continue;
      // get target values
      std::vector<float> targetRaw = {};
      if (m_neuroParameters.targetZ)
        targetRaw.push_back(zTarget);
      if (m_neuroParameters.targetTheta)
        targetRaw.push_back(thetaTarget);
      for (unsigned i = 0; i < sectors.size(); ++i) {
        int isector = sectors[i];
        std::vector<float> target = m_NeuroTrigger[isector].scaleTarget(targetRaw);
        // skip out of range targets or rescale them
        bool outOfRange = false;
        for (unsigned itarget = 0; itarget < target.size(); ++itarget) {
          if (fabs(target[itarget]) > 1.) {
            outOfRange = true;
            target[itarget] /= fabs(target[itarget]);
          }
        }
        if (!m_neuroParameters.rescaleTarget && outOfRange) continue;
        if (m_nPrepare == 0 || m_trainSets_prepare[isector].getTrackCounter() < m_nPrepare) {
          // get relative ids for all hits related to the MCParticle / RecoTrack
          // and count them to find relevant id range
          // using only related hits suppresses background EXCEPT for curling tracks
          if (m_trainOnRecoTracks) {
            RecoTrack* recoTrack =
              m_tracks[itrack]->getRelatedTo<RecoTrack>(m_targetCollectionName);
            for (const CDCTriggerSegmentHit& hit :
                 recoTrack->getRelationsTo<CDCTriggerSegmentHit>(m_hitCollectionName)) {
              // get relative id
              double relId = m_NeuroTrigger.getRelId(hit);
              m_trainSets_prepare[isector].addHit(hit.getISuperLayer(), round(relId));
            }
          } else {
            MCParticle* mcTrack =
              m_tracks[itrack]->getRelatedTo<MCParticle>(m_targetCollectionName);
            for (const CDCTriggerSegmentHit& hit :
                 mcTrack->getRelationsTo<CDCTriggerSegmentHit>(m_hitCollectionName)) {
              // get relative id
              double relId = m_NeuroTrigger.getRelId(hit);
              m_trainSets_prepare[isector].addHit(hit.getISuperLayer(), round(relId));
            }
          }
          m_trainSets_prepare[isector].countTrack();
        }
      }
    }
    bool stop = true;
    for (unsigned isector = 0; isector < m_trainSets_prepare.size(); ++isector) {
      if (m_nPrepare == 0 || m_trainSets_prepare[isector].getTrackCounter() < m_nPrepare) {
        stop = false;
        break;
      }
    }
    if (stop) {
      B2INFO("Training sample preparation for NeuroTrigger finished, stopping event loop.");
      // if required hit number is reached, get relevant ids
      StoreObjPtr<EventMetaData> eventMetaData;
      eventMetaData->setEndOfData();
    }
  }

event() [3/11]

void event ( void  )

overridevirtual

Called event by event.

Reimplemented from Module.

Definition at line 360 of file TRGCDCModule.cc.

  {
    // initialize returnValue
    _cdc->setReturnValue(0);
 
    //...CDC trigger simulation...
    _cdc->update();
    _cdc->simulate();
 
    // save the results to DataStore
    _cdc->storeSimulationResults(m_2DfinderCollectionName,
                                 m_2DfitterCollectionName,
                                 m_3DfitterCollectionName);
 
    // Set return value
    int returnValue = _cdc->getReturnValue();
    int mask = 0;
    for (vector<string>::iterator it = _returnValueModuleNames.begin(); it != _returnValueModuleNames.end(); ++it) {
      string const& moduleName = *it;
      if (moduleName == "TSF") mask |= TRGCDC::EReturnValueType::TSF;
      else if (moduleName == "ETF") mask |= TRGCDC::EReturnValueType::ETF;
      else if (moduleName == "find2D") mask |= TRGCDC::EReturnValueType::find2D;
      else if (moduleName == "fit2D") mask |= TRGCDC::EReturnValueType::fit2D;
      else if (moduleName == "find3D") mask |= TRGCDC::EReturnValueType::find3D;
      else if (moduleName == "fit3D") mask |= TRGCDC::EReturnValueType::fit3D;
    }
    if (mask != 0) returnValue &= mask;
    setReturnValue(returnValue);
  }

event() [4/11]

void event ( void  )

overridevirtual

Called event by event.

Reimplemented from Module.

Definition at line 135 of file TRGCDCTSStreamModule.cc.

  {
 
    //...To dump wire hits...
    if (TRGDebug::level())
      _cdc->dump("trgWireCentralHits");
 
    if (_out) {
      unsigned val = TRGBSRecord_BeginEvent;
      _out->write((char*) & val, 4);
      val = 0;
      _out->write((char*) & val, 4);
    }
 
    //...Clock loop (from 0 to 99 cycles, about 800 ns)...
    for (unsigned c = 0; c < 100; c++) {
 
      if (_out) {
        unsigned val = TRGBSRecord_Clock;
        _out->write((char*) & val, 4);
        val = 32;
        _out->write((char*) & val, 4);
        val = c;
        _out->write((char*) & val, 4);
      }
 
      //...Super layer loop...
      for (unsigned l = 0; l < 9; l++) {
 
        //...Bit stream for this super layer...
        TRGBitStream stream;
 
        //...Wire loop...
        for (unsigned i = 0; i < _wires[l].size(); i++) {
          const TRGSignal& s = _wires[l][i]->signal();
          bool hit = s.state(c);
          stream.append(hit);
        }
 
        if (_out) {
          unsigned val = TRGBSRecord_SegmentSL0;
          val += l;
          _out->write((char*) & val, 4);
          val = stream.size();
          _out->write((char*) & val, 4);
          for (unsigned i = 0; i < stream.sizeInChar(); i++) {
            char cs = stream.c(i);
            _out->write(& cs, 1);
          }
        }
 
        if (TRGDebug::level()) {
          cout << "Super layer " << l << ", clock " << c << endl;
          stream.dump();
        }
      }
    }
 
    if (_out) {
      unsigned val = TRGBSRecord_EndEvent;
      _out->write((char*) & val, 4);
      val = 0;
      _out->write((char*) & val, 4);
    }
 
  }

event() [5/11]

void event ( void  )

overridevirtual

Called event by event.

Reimplemented from Module.

Definition at line 149 of file TRGECLFAMModule.cc.

  {
 
    B2DEBUG(200, "TRGECLFAMMoudle ... event called");
    //
    //
    if (m_nEvent < 1e2) {if (m_nEvent %    10 == 0) {B2DEBUG(200, "TRGECLFAMModule::event> evtno= " << m_nEvent);}}
    else if (m_nEvent < 1e3) {if (m_nEvent %   100 == 0) {B2DEBUG(200, "TRGECLFAMModule::event> evtno= " << m_nEvent);}}
    else if (m_nEvent < 1e4) {if (m_nEvent %  1000 == 0) {B2DEBUG(200, "TRGECLFAMModule::event> evtno= " << m_nEvent);}}
    else if (m_nEvent < 1e5) {if (m_nEvent % 10000 == 0) {B2DEBUG(200, "TRGECLFAMModule::event> evtno= " << m_nEvent);}}
    else if (m_nEvent < 1e6) {if (m_nEvent % 100000 == 0) {B2DEBUG(200, "TRGECLFAMModule::event> evtno= " << m_nEvent);}}
 
 
    //
    //
    //
    // FAM Digitizer
    TrgEclDigitizer* obj_trgeclDigi = new TrgEclDigitizer();
    obj_trgeclDigi->setWaveform(_waveform);
    obj_trgeclDigi->setFADC(_FADC);
    obj_trgeclDigi->setup(m_SourceOfTC);
    if (_famMethod == 1 || _famMethod == 2) {
      // no-fit method = backup method 1
      obj_trgeclDigi->digitization01(TCDigiE, TCDigiT);
    } else if (_famMethod == 3) {
      // orignal method = backup method 2
      obj_trgeclDigi->digitization02(TCDigiE, TCDigiT);
    }
    obj_trgeclDigi-> save(m_nEvent);
 
 
    // FAM Fitter
    TrgEclFAMFit* obj_trgeclfit = new TrgEclFAMFit();
    obj_trgeclfit-> SetBeamBkgTagFlag(_beambkgtag);
    obj_trgeclfit-> SetAnaTagFlag(_famana);
    obj_trgeclfit-> setup(m_nEvent);
    obj_trgeclfit-> SetThreshold(Threshold);
 
    if (_famMethod == 1) {obj_trgeclfit->  FAMFit01(TCDigiE, TCDigiT); } // fitting method
    else if (_famMethod == 2) {obj_trgeclfit->  FAMFit02(TCDigiE, TCDigiT); } // no-fit method = backup method 1
    else if (_famMethod == 3) { obj_trgeclfit-> FAMFit03(TCDigiE, TCDigiT); } // orignal method = backup method 2
    obj_trgeclfit-> save(m_nEvent);
 
 
    // Count number of trigger cells in each ECL region for EventLevelClusteringInfo
    uint16_t nTCsPerRegion[3] = {};
    const double absTimeRequirement = 9999.; // Selection on time to reproduce data
    const int firstBarrelId = 81; // First TCId in the barrel
    const int lastBarrelId = 512; // Last TCId in the barrel
    for (auto& trgeclhit : m_TRGECLHit) {
      const int tcId = trgeclhit.getTCId();
      const double tcTime = trgeclhit.getTimeAve();
      if (std::abs(tcTime) < absTimeRequirement) {
        if (tcId < firstBarrelId) {
          nTCsPerRegion[0]++;
        } else if (tcId > lastBarrelId) {
          nTCsPerRegion[2]++;
        } else {
          nTCsPerRegion[1]++;
        }
      }
    }
 
    // Store
    if (!m_eventLevelClusteringInfo) { m_eventLevelClusteringInfo.create();}
    m_eventLevelClusteringInfo->setNECLTriggerCellsFWD(nTCsPerRegion[0]);
    m_eventLevelClusteringInfo->setNECLTriggerCellsBarrel(nTCsPerRegion[1]);
    m_eventLevelClusteringInfo->setNECLTriggerCellsBWD(nTCsPerRegion[2]);
 
 
    //
    //
    //
    m_nEvent++;
    delete obj_trgeclDigi;
    delete obj_trgeclfit;
    //
    //
    //
  }

event() [6/11]

void event ( void  )

overridevirtual

Called event by event.

Reimplemented from Module.

Definition at line 363 of file TRGECLModule.cc.

  {
 
    if (TRGDebug::level()) {
      std::cout << "TRGECLMoudle ... event called" << std::endl;
    }
    //
    // simulation
    //
    TrgEclMaster* _ecl = new TrgEclMaster();
 
    _ecl->initialize(m_nEvent);
    _ecl->setClusterMethod(_Clustering);
    _ecl->setClusterLimit(_ClusterLimit);
    _ecl->setBhabhaMethod(_Bhabha);
    _ecl->setEventTimingMethod(_EventTiming);
    _ecl->setTimeWindow(_TimeWindow);
    _ecl->setOverlapWindow(_OverlapWindow);
    _ecl->setNofTopTC(_NofTopTC);
    _ecl->set2DBhabhaThreshold(_2DBhabhaThresholdFWD, _2DBhabhaThresholdBWD);
    _ecl->set3DBhabhaSelectionThreshold(_3DBhabhaSelectionThreshold);
    _ecl->set3DBhabhaVetoThreshold(_3DBhabhaVetoThreshold);
    _ecl->set3DBhabhaSelectionAngle(_3DBhabhaSelectionAngle);
    _ecl->set3DBhabhaVetoAngle(_3DBhabhaVetoAngle);
    _ecl->setmumuThreshold(_mumuThreshold);
    _ecl->setmumuAngle(_mumuAngle);
    _ecl->set3DBhabhaAddAngleCut(m_3DBhabhaAddAngleCut);
    _ecl->setTaub2bAngleCut(m_taub2bAngleCut);
    _ecl->setTaub2bEtotCut(m_taub2bEtotCut);
    _ecl->setTaub2bClusterECut(m_taub2bClusterECut1,
                               m_taub2bClusterECut2);
    _ecl->setTaub2b2Cut(m_taub2b2AngleCut,
                        m_taub2b2EtotCut,
                        m_taub2b2CLEEndcapCut,
                        m_taub2b2CLECut);
    _ecl->setTaub2b3Cut(m_taub2b3AngleCut,
                        m_taub2b3EtotCut,
                        m_taub2b3CLEb2bCut,
                        m_taub2b3CLELowCut,
                        m_taub2b3CLEHighCut);
    _ecl->setTotalEnergyThreshold(_TotalEnergy);
    _ecl->setLowMultiplicityThreshold(_LowMultiThreshold);
    _ecl->setn300MeVClusterThreshold(_n300MeVCluster);
    _ecl->setECLBurstThreshold(_ECLBurstThreshold);
    _ecl->set3DBhabhaVetoInTrackThetaRegion(m_3DBhabhaVetoInTrackThetaRegion);
    _ecl->setEventTimingQualityThresholds(m_EventTimingQualityThresholds);
 
    if (_SelectEvent == 0) {
      _ecl->simulate01(m_nEvent);
    } else if (_SelectEvent == 1) {
      _ecl->simulate02(m_nEvent);
    }
    // printf("TRGECLModule> eventId = %d \n", m_nEvent);
    //
    //
    m_nEvent++;
    //
    //
    //
    delete _ecl;
  }

event() [7/11]

void event ( void  )

overridevirtual

Called event by event.

Reimplemented from Module.

Definition at line 111 of file TRGECLBGTCHitModule.cc.

  {
    B2DEBUG(200, "TRGECLBGTCHitMoudle ... event called");
 
    if (m_eclHits.getEntries() > 0) {
      genSimulationObj();
    } else if (m_trgeclUnpackerStores.getEntries() > 0) {
      genRandomTrgObj();
    }
  }

event() [8/11]

void event ( void  )

overridevirtual

Event.

Reimplemented from Module.

Definition at line 142 of file TRGECLQAMModule.cc.

  {
 
    B2DEBUG(200, "TRGECLFAMMoudle ... event called");
 
 
 
    //QAM ECL Hit Map
 
    //    StoreArray<TRGECLUnpackerStore> m_TRGECLUnpackerStore;
    for (int ii = 0; ii < m_TRGECLUnpackerStore.getEntries(); ii++) {
      TRGECLUnpackerStore* m_TRGECLUnpacker = m_TRGECLUnpackerStore[ii];
 
 
      int TCId = (m_TRGECLUnpacker->getTCId());
      int HitEnergy = (m_TRGECLUnpacker->getTCEnergy());
 
      m_etot += HitEnergy;
      if (TCId < 1 || TCId > 576 || HitEnergy == 0) {continue;}
      TCID[TCId - 1]++;
      int Caltime = m_TRGECLUnpacker->getTCCALTime();
      if (Caltime != 0) {
        caltime.push_back(Caltime);
      }
    }
 
    etot.push_back(m_etot);
 
    m_etot = 0;
 
    //Unpacker EvtStore
    //for ( int ii = 0; ii < m_TRGECLUnpackerSumStore.getEntries(); ii++){
    TRGECLUnpackerEvtStore* m_TRGECLUnpackerEvt = m_TRGECLUnpackerEvtStore[0];
 
    if (m_TRGECLUnpackerEvt->get3DBhabhaS() == 1) {
      //Cluster Energy 2D Check
      const int* clE = m_TRGECLUnpackerEvt->getCLEnergy();
 
      int clER[6] = {0};
 
      for (int i = 0; i < 6; i++) {
        clER[i] = clE[i];
      }
 
      sort(clER, clER + 6);
 
      clusterE = clER[5] + clER[4]; // Ecl1+Ecl2
      cluster.push_back(clusterE);
    }
 
    StoreObjPtr<EventMetaData> bevt;
    m_nRun = bevt->getRun();
    m_nevent++;
 
  }

event() [9/11]

void event ( void  )

overridevirtual

Called event by event.

Reimplemented from Module.

Definition at line 110 of file TRGECLRawdataAnalysisModule.cc.

  {
 
    if (TRGDebug::level()) {
      std::cout << "TRGECLMoudle ... event called" << std::endl;
    }
 
    std::vector<int> TCId;
    std::vector<double> TCTiming;
    std::vector<double> TCEnergy;
 
    TCId.clear();
    TCTiming.clear();
    TCEnergy.clear();
 
    int HitRevoFAM = 0;
    int HitFineTime = 0;
 
 
    //
    StoreArray<TRGECLUnpackerStore> trgeclUnpackerStoreArray;
    for (int ii = 0; ii < trgeclUnpackerStoreArray.getEntries(); ii++) {
 
      TRGECLUnpackerStore* TCHit = trgeclUnpackerStoreArray[ii];
      int  iTCID = (TCHit->getTCId() - 1);
      int HitTiming    = TCHit ->getTCTime();
      int HitEnergy =  TCHit -> getTCEnergy();
      HitRevoFAM = TCHit -> getRevoFAM();
      HitFineTime = TCHit -> getTCTime();
 
      if (iTCID == -1) {continue;}
 
      TCId.push_back(iTCID + 1);
      TCTiming.push_back(HitTiming);
      TCEnergy.push_back(HitEnergy);
 
 
 
    }
    TrgEclCluster obj_cluster;
    if (TCId.size() > 0) {
 
      obj_cluster.setClusteringMethod(_Clustering);
      obj_cluster.setEventId(m_nEvent);
      obj_cluster.setICN(TCId, TCEnergy, TCTiming); // Make Cluster
      obj_cluster.save(m_nEvent); // Save Clusters to TRGECLCluster
 
    }
 
    int Timing = ((HitFineTime >> 3) & 0xF) + ((HitRevoFAM & 0x7F) << 4);
 
    m_hitNum = 0;
    StoreArray<TRGECLTiming> TimingArray;
    TimingArray.appendNew();
    m_hitNum = TimingArray.getEntries() - 1;
    TimingArray[m_hitNum]->setEventId(m_nEvent);
    TimingArray[m_hitNum]->setTiming(Timing);
 
 
    //
 
    //
    //
    m_nEvent++;
    //
    //
    //
 
 
  }

event() [10/11]

void event ( void  )

overridevirtual

Called event by event.

Reimplemented from HistoModule.

Definition at line 185 of file TRGGDLModule.cc.

  {
 
    if (_debugLevel > 9) printf("TRGGDLModule::event() starts.\n");
 
    newDir->cd();
 
    //...GDL simulation...
    _gdl->update(true);
    _gdl->simulate();
    _gdl->accumulateInp(h_inp);
    _gdl->accumulateFtd(h_ftd);
    _gdl->accumulatePsn(h_psn);
 
    //StoreObjPtr<TRGSummary> m_TRGSummary; /**< output for TRGSummary */
    bool result_summary = false;
    if (m_TRGSummary) {
      result_summary = m_TRGSummary->test();
    } else {
      B2WARNING("TRGGDLModule.cc: TRGSummary not found. Check it!!!!");
    }
    setReturnValue(result_summary);
 
    oldDir->cd();
  }

event() [11/11]

void event ( void  )

overridevirtual

Called event by event.

Reimplemented from Module.

Definition at line 110 of file TRGGRLModule.cc.

  {
 
    //...GRL simulation...
    _grl->update(true);
    _grl->simulate();
 
  }

genRandomTrgObj()

void genRandomTrgObj ( )

private

Definition at line 181 of file TRGECLBGTCHitModule.cc.

  {
 
    for (const TRGECLUnpackerStore& ttt : m_trgeclUnpackerStores) {
      int tcID = ttt.getTCId();
      // TC energy(ADC)
      double tcEnergyADC = (double) ttt.getTCEnergy();
      // TC timing(ns)
      double tcTime   = (double) ttt.getTCTime();
      if (tcID <= 0) {continue;}
      // ADC to Energy(MeV) conversion
      double adc2energy = 5.231;
      // TC energy (GeV)
      double tcEnergy = tcEnergyADC * (adc2energy / 1000);
      // TC energy cut to reduce object size
      if (tcEnergy < m_TCEnergyCut) { continue; }
      // TC timing cut to reduce object size
      if (tcTime < m_TCTimingCutLow || tcTime > m_TCTimingCutHigh) { continue; }
      // store data
      m_trgeclBGTCHits.appendNew(tcID,
                                 tcEnergy,
                                 tcTime);
    }
  }

genSimulationObj()

void genSimulationObj ( )

private

Parameters.

Definition at line 141 of file TRGECLBGTCHitModule.cc.

  {
 
    struct hit_t { double e, t; };
    map<int, hit_t> a;
 
    for (const ECLHit& t : m_eclHits) {
      double txtal = t.getTimeAve();
      // timing cut, +-8ns, (almost) same as ECL in SensitiveDetector.cc
      if (txtal >= 8000.0 || txtal < -8000.0) continue;
      double edep = t.getEnergyDep();
      int TimeIndex = (txtal + 8000.0) * (1. / 100);
      int cellId = t.getCellId() - 1;
      int iTCIdm = m_TCMap->getTCIdFromXtalId(cellId + 1) - 1;
 
      int key = iTCIdm * 160 + TimeIndex;
      hit_t& h = a[key];
      double old_edep = h.e, old_txtal = h.t;
      double new_edep = old_edep + edep;
      h.e = new_edep;
      h.t = (old_edep * old_txtal + edep * txtal) / new_edep;
    }
    // save TC data to dataobject
    for (pair<int, hit_t> t : a) {
      int key  = t.first;
      int tcid = key / 160;
      const hit_t& h = t.second;
      // TC energy cut to reduce object size
      if (h.e < m_TCEnergyCut) { continue; }
      // TC timing cut to reduce object size
      if (h.t < m_TCTimingCutLow || h.t > m_TCTimingCutHigh) { continue; }
      // store data
      m_trgeclBGTCHits.appendNew(tcid + 1, h.e, h.t);
    }
 
  }

get_BKLM_back_to_back_flag()

int get_BKLM_back_to_back_flag

(

int

subdetector

)

Definition at line 157 of file KLM_Trig.cc.

  {
    for (const auto& e : m_summary1) {
      if (e.subdetector == subdetector) {
        return e.back2back;
      }
    }
    return 0;
  }

get_index()

constexpr int64_t get_index ( const T &  e, int  track_id = 0  )

constexpr

Definition at line 43 of file klm_trig_linear_fit.cc.

  {
    int64_t ret = 0;
    int64_t shift = 0;
    bitshift(ret, shift, track_id, 5);
    bitshift(ret, shift, e.layer,  6);
    bitshift(ret, shift, e.plane, 2);
    bitshift(ret, shift, e.sector, 6);
    bitshift(ret, shift, e.section, 3);
    bitshift(ret, shift, e.subdetector, 3);
 
 
    return ret;
  };

get_n_sections_trig()

int get_n_sections_trig

(

int

subdetector

)

Definition at line 137 of file KLM_Trig.cc.

  {
    for (const auto& e : m_summary1) {
      if (e.subdetector == subdetector) {
        return e.n_sections_trig;
      }
    }
    return 0;
  }

get_result()

const KLM_trig_linear_fits & get_result

(

)

const

Definition at line 275 of file klm_trig_linear_fit.cc.

  {
    return m_linear_fits;
  }

get_triggermask() [1/2]

int get_triggermask	(	int	subdetector,
		int	section
	)

Definition at line 147 of file KLM_Trig.cc.

  {
    for (const auto& e : m_sections_trig) {
      if (e.subdetector == subdetector && e.section == section) {
        return e.sector_mask;
      }
    }
    return 0;
  }

get_triggermask() [2/2]

int get_triggermask	(	int	subdetector,
		int	section
	)

Definition at line 282 of file klm_trig_linear_fit.cc.

  {
    for (const auto& e : m_sections_trig) {
      if (e.subdetector == subdetector && e.section == section) {
        return e.sector_mask;
      }
    }
    return 0;
  }

get_triggermask_or()

int get_triggermask_or	(	int	subdetector,
		int	section
	)

Definition at line 292 of file klm_trig_linear_fit.cc.

  {
    for (const auto& e : m_sections_trig) {
      if (e.subdetector == subdetector && e.section == section) {
        return e.sector_mask_or;
      }
    }
    return 0;
  }

initialize() [1/11]

void initialize ( void  )

overridevirtual

Initialize the module.

Initialize the networks and register datastore objects.

Reimplemented from Module.

Definition at line 83 of file CDCTriggerNeuroDataModule.cc.

  {
    // register store objects
    m_tracks.isRequired(m_inputCollectionName);
    // decided wether to load mcparticles or recotracks as training targets.
    // This has to be done because the train target values are obtained
    // for both classes in a different way.
    if (m_trainOnRecoTracks) {
      StoreArray<RecoTrack> targets(m_targetCollectionName);
      targets.isRequired(m_targetCollectionName);
    } else {
      StoreArray<MCParticle> targets(m_targetCollectionName);
      targets.isRequired(m_targetCollectionName);
    }
    // initialize the neurotrigger object, but use the parameters given in the module
    if (m_configFileName != "") {
      m_neuroParameters.loadconfigtxt(m_configFileName);
      m_NeuroTrigger.initialize(m_neuroParameters);
    } else {
      B2ERROR("The Neurotrigger needs to be initialized by a configuration file! Make sure to give the configuration file as a parameter.");
    }
    // in this version, we first need an idhistfile in prior to collect the training data.
    // this idhistfile is created by running another steering file before
    // and stored in the config file loaded in the previous step.
    m_trainSet.clear();
    CDC::CDCGeometryPar& cdc = CDC::CDCGeometryPar::Instance();
    // create an empty dataset of training data for each expert network
    for (unsigned iMLP = 0; iMLP < m_NeuroTrigger.nSectors(); ++iMLP) {
      m_trainSet.push_back(CDCTriggerMLPData());
      int layerId = 3;
      for (int iSL = 0; iSL < 9; ++iSL) {
        m_trainSet[iMLP].addCounters(cdc.nWiresInLayer(layerId));
        layerId += (iSL > 0 ? 6 : 7);
      }
    }
    // this one sets up the other root store arrays needed to collect training data
    m_NeuroTrigger.initializeCollections(m_hitCollectionName, m_EventTimeName, m_neuroParameters.et_option());
    // consistency check of training parameters
    if (m_NeuroTrigger.nSectors() != m_trainSet.size()) {
      B2ERROR("Number of ID sets (" << m_trainSet.size() << ") should match " <<
              "number of sectors (" << m_NeuroTrigger.nSectors() << ")");
    }
    // overwrite previous file with empty file, in case it already exists
    std::ofstream gzipfile4(m_filename, std::ios_base::trunc | std::ios_base::binary);
    boost::iostreams::filtering_ostream outStream;
    outStream.push(boost::iostreams::gzip_compressor());
    outStream.push(gzipfile4);
    CDCTriggerMLPData::NeuroSet<27, 2> sample;
    outStream << sample.headline << std::endl;
  }

initialize() [2/11]

void initialize ( void  )

overridevirtual

Initialize the module.

Initialize the networks and register datastore objects.

First Initialize the neurotrigger object, which is needed because of the original structure of the training procedure in the old module. Second, the dataset for the idhist data will be initialized and set up.

Reimplemented from Module.

Definition at line 70 of file CDCTriggerNeuroIDHistModule.cc.

  {
    if (m_configFileName != "") {
      m_neuroParameters.loadconfigtxt(m_configFileName);
      m_NeuroTrigger.initialize(m_neuroParameters);
    } else {
      m_neuroParameters.saveconfigtxt("neuroconfig_example.conf");
      B2ERROR("Configuration file is missing! Make sure to give the configuration file as a parameter. \
              An example file neuroconfig_example.conf has been saved.");
    }
    m_trainSets_prepare.clear();
    CDC::CDCGeometryPar& cdc = CDC::CDCGeometryPar::Instance();
    for (unsigned iMLP = 0; iMLP < m_NeuroTrigger.nSectors(); ++iMLP) {
      // sectors means experts here; this is the old naming within the neurotrigger
      // class which originates from times were we had different expert networks for
      // different geometrical sectors. In future, this might come back because of
      // the 3DFinder and its rough information about the theta angle of the tracks.
      m_trainSets_prepare.push_back(CDCTriggerMLPData());
      // layerid is the layer number where the priority hits are.
      // for every 1st priority wire there is a corresponding track segment.
      int layerId = 3;
      for (int iSL = 0; iSL < 9; ++iSL) {
 
        m_trainSets_prepare[iMLP].addCounters(cdc.nWiresInLayer(layerId));
        // the first superlayer has 2 layers extra compared to the rest
        layerId += (iSL > 0 ? 6 : 7);
      }
      m_tracks.isRequired(m_inputCollectionName);
      if (m_trainOnRecoTracks) {
        StoreArray<RecoTrack> targets(m_targetCollectionName);
        targets.isRequired(m_targetCollectionName);
      } else {
        StoreArray<MCParticle> targets(m_targetCollectionName);
        targets.isRequired(m_targetCollectionName);
      }
 
      m_NeuroTrigger.initializeCollections(m_hitCollectionName);
 
 
 
 
    }
    if (m_NeuroTrigger.nSectors() == 0) {
      B2ERROR("No networks defined, please make sure to have a proper configuration file! Example file will be created here: ./neurotrigger_default.conf");
      m_neuroParameters.saveconfigtxt("neurotrigger_default.conf");
    }
  }

initialize() [3/11]

void initialize ( void  )

overridevirtual

Initilizes TRGCDCModule.

Reimplemented from Module.

Definition at line 257 of file TRGCDCModule.cc.

  {
 
    //...Stop curl buck...
    if (_curlBackStop) {
      G4RunManager* g4rm = G4RunManager::GetRunManager();
      if (g4rm != 0) {
        if (g4rm->GetUserPhysicsList() != 0) {
          _sa = new TCSAction();
          g4rm->SetUserAction(_sa);
        }
      }
    }
    TRGDebug::level(_debugLevel);
 
    if (TRGDebug::level()) {
      cout << "TRGCDCModule::initialize ... options" << endl;
      cout << TRGDebug::tab(4) << "debug level = " << TRGDebug::level()
           << endl;
      cout << TRGDebug::tab(4) << "back stop = " << _curlBackStop << endl;
      cout << TRGDebug::tab(4) << "2D finder perfect = " << _perfect2DFinder
           << endl;
      cout << TRGDebug::tab(4) << "Hough finder Min. Peak = "
           << _hFinderPeakMin
           << endl;
    }
 
    // register DataStore elements
    StoreArray<CDCTriggerSegmentHit> segmentHits;
    segmentHits.registerInDataStore();
    StoreArray<CDCHit> cdcHits(m_inputCollectionName);
    cdcHits.isRequired();
    StoreArray<MCParticle> mcparticles;
    segmentHits.registerRelationTo(cdcHits);
    mcparticles.registerRelationTo(segmentHits);
    // register tracks only for full simulation (not for TSF simulation)
    if (!(_fastSimulationMode & 1)) {
      StoreArray<CDCTriggerTrack> tracks2Dfinder(m_2DfinderCollectionName);
      StoreArray<CDCTriggerTrack> tracks2Dfitter(m_2DfitterCollectionName);
      StoreArray<CDCTriggerTrack> tracks3Dfitter(m_3DfitterCollectionName);
      tracks2Dfinder.registerInDataStore(m_2DfinderCollectionName);
      tracks2Dfitter.registerInDataStore(m_2DfitterCollectionName);
      tracks3Dfitter.registerInDataStore(m_3DfitterCollectionName);
      tracks2Dfinder.registerRelationTo(segmentHits); // hits related over Hough cell
      tracks2Dfitter.registerRelationTo(segmentHits); // hits used for the 2D fit
      tracks3Dfitter.registerRelationTo(segmentHits); // hits used for the 2D and 3D fit
      tracks2Dfinder.registerRelationTo(tracks2Dfitter);
      tracks2Dfinder.registerRelationTo(tracks3Dfitter);
    }
    StoreObjPtr<TRGTiming> eventTime("CDCTriggerEventTime");
    eventTime.registerInDataStore("CDCTriggerEventTime");
 
    //...CDC trigger config. name...
    static string cfn = _configFilename;
 
    //...CDC trigger...
    if ((cfn != _configFilename) || (_cdc == 0)) {
      Belle2::CDC::CDCGeometryPar& cdc = Belle2::CDC::CDCGeometryPar::Instance();
      _cdc = TRGCDC::getTRGCDC(_configFilename,
                               _simulationMode,
                               _fastSimulationMode,
                               _firmwareSimulationMode,
                               _firmwareSimulationStart,
                               _firmwareSimulationStop,
                               _fileTRGCDC,
                               _perfect2DFinder,
                               _perfect3DFinder,
                               _innerTSLUTFilename,
                               _outerTSLUTFilename,
                               _rootTRGCDCFilename,
                               _rootFitter3DFilename,
                               _hFinderPeakMin,
                               _hFinderMappingFilePlus,
                               _hFinderMappingFileMinus,
                               _hFinderDoit,
                               _fLogicLUTTSF,
                               _fLRLUT,
                               _fFitter3Dsmclr,
                               _fFitter3Ds2DFit,
                               _fFitter3Ds2DFitDrift,
                               _wireHitInefficiency,
                               _fileTSF,
                               _fileETF,
                               _fverETF,
                               _fprintFirmETF,
                               _fileHough3D,
                               _finder3DMode,
                               _fileFitter3D,
                               _fXtSimpleFitter3D,
                               cdc.getTdcBinWidth(),
                               _trgCDCDataInputMode,
                               m_inputCollectionName);
    }
  }

initialize() [4/11]

void initialize ( void  )

overridevirtual

Initilizes TRGCDCTSStreamModule.

Reimplemented from Module.

Definition at line 81 of file TRGCDCTSStreamModule.cc.

  {
    if (TRGDebug::level()) {
      cout << "TRGCDCTSStreamModule::initialize ... options" << endl;
      cout << TRGDebug::tab(4) << "debug level = " << _debugLevel << endl;
      cout << TRGDebug::tab(4) << "       mode = " << _mode << endl;
      cout << TRGDebug::tab(4) << "output file = " << _streamFilename << endl;
    }
 
    if (_streamFilename != "unknown") {
      _out = new ofstream(_streamFilename.c_str(), ios::out | ios::binary);
      unsigned val = TRGBSRecord_Comment;
      _out->write((char*) & val, 4);
      const string cmt = "test data ";
      val = cmt.size() * 8;
      _out->write((char*) & val, 4);
      _out->write(cmt.c_str(), cmt.size());
    }
  }

initialize() [5/11]

void initialize ( void  )

overridevirtual

Initilizes TRGECLFAMModule.

Reimplemented from Module.

Definition at line 101 of file TRGECLFAMModule.cc.

  {
 
    B2DEBUG(100, "TRGECLFAMModule::initialize ... options");
    B2DEBUG(100, "TRGECLFAMModule::initialize> FAM Fit Method = " << _famMethod
            << "  ; Bin of Time Interval = " << _binTimeInterval
            << " ;output TC waveforml = " << _waveform);
    if (m_SourceOfTC <= 0 ||
        m_SourceOfTC >= 4) {
      B2FATAL("TRGECLFAMModule::initialize> SourceOfTC must be 1 or 2 or 3");
    }
 
    m_nRun   = 0;
    m_nEvent = 1;
 
    m_TRGECLDigi0.registerInDataStore();
    m_TRGECLWaveform.registerInDataStore();
    m_TRGECLHit.registerInDataStore();
    m_TRGECLFAMAna.registerInDataStore();
    // This object is registered by few packages. Let's be agnostic about the
    // execution order of the modules: the first package run registers the module
    m_eventLevelClusteringInfo.isOptional() ? m_eventLevelClusteringInfo.isRequired() :
    m_eventLevelClusteringInfo.registerInDataStore();
 
    //    m_FAMPara = new DBObjPtr<TRGECLFAMPara>;
  }

initialize() [6/11]

void initialize ( void  )

overridevirtual

Initilizes TRGECLModule.

Reimplemented from Module.

Definition at line 228 of file TRGECLModule.cc.

  {
 
    TRGDebug::level(_debugLevel);
 
    if (TRGDebug::level()) {
      std::cout << "TRGECLModule::initialize ... options" << std::endl;
      std::cout << TRGDebug::tab(4) << "debug level = " << TRGDebug::level()
                << std::endl;
    }
    //
    //
    //
    m_nRun   = 0 ;
    m_nEvent = 0 ;
    m_hitNum = 0;
    m_hitTCNum = 0;
 
    m_TRGECLHit.registerInDataStore();
    m_TRGECLTrg.registerInDataStore();
    m_TRGECLCluster.registerInDataStore();
 
  }

initialize() [7/11]

void initialize ( void  )

overridevirtual

Initilizes TRGECLBGTCHitModule.

Reimplemented from Module.

Definition at line 84 of file TRGECLBGTCHitModule.cc.

  {
    B2DEBUG(100, "TRGECLBGTCHitModule::initialize>");
 
    B2INFO("[TRGECLBGTCHitModule] TCEnergyCut     : TC(E) GeV = "
           << m_TCEnergyCut);
    B2INFO("[TRGECLBGTCHitModule] TCTimingCutLow  : TC(T) ns  = "
           << m_TCTimingCutLow);
    B2INFO("[TRGECLBGTCHitModule] TCTimingCutHigh : TC(T) ns  = "
           << m_TCTimingCutHigh);
 
    m_eclHits.registerInDataStore();
    m_trgeclUnpackerStores.registerInDataStore();
    m_trgeclBGTCHits.registerInDataStore();
  }

initialize() [8/11]

void initialize ( void  )

overridevirtual

initialize

Reimplemented from Module.

Definition at line 114 of file TRGECLQAMModule.cc.

  {
 
    B2DEBUG(100, "TRGECLQAMModule::initialize ... options");
 
 
    m_TRGECLUnpackerStore.registerInDataStore();
    m_TRGECLUnpackerEvtStore.registerInDataStore();
 
    m_TRGSummary.registerInDataStore();
    //  EvtMeta.registterInDataStore();
 
 
  }

initialize() [9/11]

void initialize ( void  )

overridevirtual

Initilizes TRGECLModule.

Reimplemented from Module.

Definition at line 70 of file TRGECLRawdataAnalysisModule.cc.

  {
 
    TRGDebug::level(_debugLevel);
 
    if (TRGDebug::level()) {
      std::cout << "TRGECLRawdataAnalysisModule::initialize ... options" << std::endl;
      std::cout << TRGDebug::tab(4) << "debug level = " << TRGDebug::level()
                << std::endl;
    }
    //
    //
    //
    m_nRun   = 0 ;
    m_nEvent = 0 ;
    m_hitNum = 0;
    m_hitTCNum = 0;
 
    m_TRGECLUnpackerStore.registerInDataStore();
    m_TRGECLCluster.registerInDataStore();
    m_TRGECLTiming.registerInDataStore();
 
  }

initialize() [10/11]

void initialize ( void  )

overridevirtual

Initilizes TRGGDLModule.

Reimplemented from HistoModule.

Definition at line 126 of file TRGGDLModule.cc.

  {
 
    REG_HISTOGRAM
    defineHisto();
 
//  m_TRGSummary.isRequired();
    TRGDebug::level(_debugLevel);
 
    B2INFO("TRGGDLModule::initialize. simulationMode=" << _simulationMode);
    if (_simulationMode != 3) {
      m_TRGGRLInfo.isRequired("TRGGRLObjects");
    }
 
    if (m_simulateT0Jitter)
      m_simClockState.isRequired();
 
    m_TRGSummary.registerInDataStore();
  }

initialize() [11/11]

void initialize ( void  )

overridevirtual

Initilizes TRGGRLModule.

Reimplemented from Module.

Definition at line 77 of file TRGGRLModule.cc.

  {
 
//  TRGDebug::level(_debugLevel);
 
    B2DEBUG(100, "TRGGRLModule::initialize ... options");
  }

operator==()

bool operator== ( const CDCTriggerHoughCand &  a, const CDCTriggerHoughCand &  b  )

inline

Two cells are identical if they have the same coordinates.

Definition at line 77 of file CDCTriggerHoughETFModule.h.

  {
    return (a.getCoord().first.X() == b.getCoord().first.X() &&
            a.getCoord().first.Y() == b.getCoord().first.Y());
  }

run() [1/2]

void run

(

const KLM_TRG_definitions::KLM_Digit_compact_ts &

hits

)

Definition at line 54 of file KLM_Trig.cc.

  {
 
 
    m_sections_trig.clear();
    m_summary1.clear();
 
    auto hits1 = nt::algorithms::join_vectors(
                   hits,
                   m_layersUsed,
                   nt::comparators::on_common_args,
    [&](auto e1, const auto&) {
      e1.layer = e1.layer * 2 + e1.plane;
      return e1;
    }
                 );
 
    __CSV__WRITE__(hits1);
 
 
    auto grouped = nt_group(
                     hits1[0].subdetector,
                     hits1[0].section,
                     hits1[0].sector
                   ).apply_append(
                     hits,
    [](const auto & e1) {
      auto bit_mask = to_bit_mask<layer>(e1);
      return nt::ntuple(
               ax_maker(layer_count) = countBits(bit_mask),
               ax_maker(layer_mask) = bit_mask
             );
    });
 
 
    __CSV__WRITE__(grouped);
 
    nt::algorithms::filter(grouped, [&](const auto & e) { return e.layer_count >= m_NlayerTrigger;  });
 
 
 
 
 
    m_sections_trig = nt_group(
                        grouped[0].subdetector,
                        grouped[0].section
                      ).apply_append(
                        grouped,
    [&](const auto & e1) {
      auto bit_mask = to_bit_mask<sector>(e1);
      n_sections_trig  nsector = countBits(bit_mask);
      return  nt::ntuple(
                sector_mask(bit_mask),
                nsector
              );
    }
                      );
 
    __CSV__WRITE__(m_sections_trig);
 
 
 
    m_summary1 = nt_group(
                   m_sections_trig[0].subdetector
                 ).apply_append(
                   m_sections_trig,
    [](const auto & e1) {
      auto bit_mask = countBits(to_bit_mask<section>(e1));
      n_sections_trig  nsector = nt::algorithms::sum(e1, [](auto&& e) {return e.n_sections_trig; });
 
      return  nt::ntuple(
                nsector,
                back2back(bit_mask == 2)
 
              );
    }
                 );
 
    __CSV__WRITE__(m_summary1);
 
 
  }

run() [2/2]

void run

(

const KLM_TRG_definitions::KLM_Digit_compact_ts &

hits

)

Definition at line 167 of file klm_trig_linear_fit.cc.

  {
 
    __CSV__WRITE__(hits);
    m_linear_fits.clear();
 
 
    auto hits_w_geo_fit = nt::algorithms::add_column(
                            hits,
    [&](const auto & e1) {
      auto ret = nt::ntuple(
                   Belle2::KLM_TRG_definitions::x_pos(0),
                   Belle2::KLM_TRG_definitions::y_pos(0)
                 );
 
      auto&& e2 = m_KLMgeomap[get_index(e1)];
      ret.x_pos.v = e1.layer * e2.slopeX + e2.offsetX;
      ret.y_pos.v = e1.strip * e2.slopeY + e2.offsetY;
      return ret;
 
    }
                          );
 
 
    __CSV__WRITE__(hits_w_geo_fit);
    nt::algorithms::sort(hits_w_geo_fit);
 
 
    auto track_maker_ = track_maker_t{};
    track_maker_.y_cutoff = y_cutoff;
 
 
    auto hits_w_geo_fit_w_tracks = nt::algorithms::add_column(
                                     hits_w_geo_fit,
    [&](const auto & e) {
      return nt::ntuple{
        track_id(track_maker_(e))
      };
 
 
    }
 
                                   );
 
 
    __CSV__WRITE__(hits_w_geo_fit_w_tracks);
 
 
 
    m_linear_fits = nt_group(
                      hits_w_geo_fit_w_tracks[0].subdetector,
                      hits_w_geo_fit_w_tracks[0].section,
                      hits_w_geo_fit_w_tracks[0].sector,
                      hits_w_geo_fit_w_tracks[0].plane,
                      hits_w_geo_fit_w_tracks[0].track_id
                    ).apply_append(
                      hits_w_geo_fit_w_tracks,
                      Linear_fit_of_Hits
                    );
    __CSV__WRITE__(m_linear_fits);
 
 
    nt::algorithms::filter(m_linear_fits, [&](const auto & e) { return std::abs(e.interceptXY) <= m_intercept_cutoff; });
 
 
    auto m_linear_fits1 = nt_group(
                            m_linear_fits[0].subdetector,
                            m_linear_fits[0].section,
                            m_linear_fits[0].sector
                          ).apply_append(
                            m_linear_fits,
    [](const auto & es) {
      bool plane0 = false;
      bool plane1 = false;
      for (const auto& e : es) {
        if (e.plane == 0) { plane0 = true; }
        if (e.plane == 1) { plane1 = true; }
      }
      return nt::ntuple{
        ax_maker(trigger_or) = (int)(plane0 || plane1),
        ax_maker(trigger_and) = (int)(plane0 && plane1)
      };
    }
                          );
    __CSV__WRITE__(m_linear_fits1);
 
    m_sections_trig = nt_group(
                        m_linear_fits1[0].subdetector,
                        m_linear_fits1[0].section
                      ).apply_append(
                        m_linear_fits1,
    [](const auto & e) {
      auto ret = nt::ntuple{
        KLM_TRG_definitions::sector_mask {},
        KLM_TRG_definitions::sector_mask_or{}
      };
      auto trig_and = nt::algorithms::filter_copy(e, [](const auto & x) { return x.trigger_and > 0; });
      ret.sector_mask.v = Belle2::to_bit_mask< sector >(trig_and);
      ret.sector_mask_or.v = Belle2::to_bit_mask< sector >(e);
      return ret;
 
    }
                      );
    __CSV__WRITE__(m_sections_trig);
 
 
  }

set_intercept_cutoff()

void set_intercept_cutoff

(

int

cutoff

)

Definition at line 307 of file klm_trig_linear_fit.cc.

  {
    m_intercept_cutoff = cutoff;
  }

set_NLayerTrigger()

void set_NLayerTrigger

(

int

NlayerTrigger

)

Definition at line 49 of file KLM_Trig.cc.

  {
    m_NlayerTrigger = NlayerTrigger;
  }

set_y_cutoff()

void set_y_cutoff

(

int

cutoff

)

Definition at line 302 of file klm_trig_linear_fit.cc.

  {
    y_cutoff = cutoff;
  }

terminate() [1/11]

void terminate ( void  )

overridevirtual

Do the training for all sectors.

Reimplemented from Module.

Definition at line 287 of file CDCTriggerNeuroDataModule.cc.

  {
    std::stringstream ss;
    for (unsigned int i = 0; i < m_trainSet.size(); ++i) {
      ss << "expert " << i << " : " << m_trainSet[i].nSamples() << ", ";
    }
    // lock the parameters, which were used to obtain the training data, so they cannot be altered in the next steps.
    m_neuroParameters.ETOption.lock();
    m_neuroParameters.rescaleTarget.lock();
    m_neuroParameters.targetZ.lock();
    m_neuroParameters.targetTheta.lock();
    for (auto x : m_neuroParameters.maxHitsPerSL) {
      x.lock();
    }
    if (m_writeconfigFileName == "") {
      B2INFO("write writeconfig now: " + m_configFileName);
      m_neuroParameters.saveconfigtxt(m_configFileName);
    } else {
      B2INFO("write writeconfig now: " + m_writeconfigFileName);
      m_neuroParameters.saveconfigtxt(m_writeconfigFileName);
    }
 
    B2DEBUG(10, "Collected events: " << ss.str());
  }

terminate() [2/11]

void terminate ( void  )

overridevirtual

Do the training for all sectors.

Reimplemented from Module.

Definition at line 213 of file CDCTriggerNeuroIDHistModule.cc.

  {
 
    if (m_neuroParameters.IDRanges.size() > 0) {
      if (m_neuroParameters.IDRanges[0].size() > 0) {
        if (m_neuroParameters.IDRanges[0][0].isSet()) {
          // the idranges are already set, print warning:
          if (!m_neuroParameters.IDRanges[0][0].isLocked()) {
            B2WARNING("ID ranges are already set in the config file, they will be updated now!");
          } else {
            B2ERROR("The ID ranges in the config file are already locked and cannot be updated!");
            return;
          }
        }
      }
    }
    m_neuroParameters.IDRanges.clear();
    for (unsigned isector = 0; isector < m_trainSets_prepare.size(); ++isector) {
      CDCTriggerMLPData::HeaderSet hset(isector, NeuroTrainer::getRelevantID(
                                          m_trainSets_prepare[isector],
                                          m_neuroParameters.cutSum,
                                          m_neuroParameters.relevantCut));
      std::vector<NNTParam<float>> expertline;
      expertline.push_back(float(isector));
      expertline.back().lock();
      for (auto x : hset.relID) {
        expertline.push_back(x);
        expertline.back().lock();
      }
      m_neuroParameters.IDRanges.push_back(expertline);
      B2DEBUG(15, hset);
    }
    // lock the variables used in this module, that are not supposed be changed
    // further down the training chain because of the danger of implications or
    // wrong assumptions.
    m_neuroParameters.relevantCut.lock();
    m_neuroParameters.cutSum.lock();
    // the IDRanges are set here; however, they can be altered manually in the
    // configuration file to achieve potentially better results, eg. widen
    // the range of the axial phi acceptance. This is why they are not to be "locked".
    m_neuroParameters.nInput.lock();
    m_neuroParameters.nOutput.lock();
    m_neuroParameters.nMLP.lock();
    for (auto x : m_neuroParameters.SLpattern) {
      x.lock();
    }
    for (auto x : m_neuroParameters.SLpatternMask) {
      x.lock();
    }
    for (auto x : m_neuroParameters.maxHitsPerSL) {
      x.lock();
    }
    m_neuroParameters.multiplyHidden.lock();
 
    if (m_writeconfigFileName == "") {
      m_neuroParameters.saveconfigtxt(m_configFileName);
    } else {
      m_neuroParameters.saveconfigtxt(m_writeconfigFileName);
    }
 
 
    // the *rangeTrain variables are used here, but just for obtaining the idranges.
    // because they only have a very minor effect on those, they are not locked here.
    //TODO: also write the config file to be directly able to start the training
  }

terminate() [3/11]

void terminate ( void  )

overridevirtual

Called when processing ended.

Reimplemented from Module.

Definition at line 398 of file TRGCDCModule.cc.

  {
 
    _cdc->terminate();
 
    if (TRGDebug::level())
      cout << "TRGCDCModule ... terminate called " << endl;
  }

terminate() [4/11]

void terminate ( void  )

overridevirtual

Called when processing ended.

Reimplemented from Module.

Definition at line 217 of file TRGCDCTSStreamModule.cc.

  {
    if (_out) {
      _out->flush();
      _out->close();
    }
 
    if (TRGDebug::level())
      cout << "TRGCDCTSStreamModule ... terminate called " << endl;
  }

terminate() [5/11]

void terminate ( void  )

overridevirtual

Called when processing ended.

Reimplemented from Module.

Definition at line 241 of file TRGECLFAMModule.cc.

  {
    B2DEBUG(100, "TRGECLFAMModule ... terminate called ");
  }

terminate() [6/11]

void terminate ( void  )

overridevirtual

Called when processing ended.

Reimplemented from Module.

Definition at line 438 of file TRGECLModule.cc.

  {
    if (TRGDebug::level()) {
      std::cout << "TRGECLModule ... terminate called " << std::endl;
    }
  }

terminate() [7/11]

void terminate ( void  )

overridevirtual

Called when processing ended.

Reimplemented from Module.

Definition at line 133 of file TRGECLBGTCHitModule.cc.

  {
    B2DEBUG(100, "TRGECLBGTCHitModule ... terminate called ");
  }

terminate() [8/11]

void terminate ( void  )

overridevirtual

terminate

Reimplemented from Module.

Definition at line 410 of file TRGECLQAMModule.cc.

  {
    B2DEBUG(100, "TRGECLQAMModule ... terminate called ");
 
  }

terminate() [9/11]

void terminate ( void  )

overridevirtual

Called when processing ended.

Reimplemented from Module.

Definition at line 194 of file TRGECLRawdataAnalysisModule.cc.

  {
    if (TRGDebug::level()) {
      std::cout << "TRGECLRawdataAnalysisModule ... terminate called " << std::endl;
    }
  }

terminate() [10/11]

void terminate ( void  )

overridevirtual

Called when processing ended.

Reimplemented from HistoModule.

Definition at line 218 of file TRGGDLModule.cc.

  {
 
    _gdl->terminate();
 
    B2DEBUG(100, "TRGGDLModule ... terminate called ");
  }

terminate() [11/11]

void terminate ( void  )

overridevirtual

Called when processing ended.

Reimplemented from Module.

Definition at line 127 of file TRGGRLModule.cc.

  {
    _grl->terminate();
    B2DEBUG(200, "TRGGRLModule ... terminate called ");
  }

to_bit_mask()

uint64_t to_bit_mask

(

const CONTAINER_T &

container

)

Definition at line 30 of file bit_operations.h.

  {
 
    return std::accumulate(container.begin(), container.end(), uint64_t(0),
    [](const auto & lhs, const auto & rhs) {
      const auto bitshift = uint64_t(AXIS_NAME_T::get(rhs));
      if (bitshift > 32) {
        throw std::runtime_error("from:\nuint64_t to_bit_mask(const CONTAINER_T& container)\ninput number to large.\n\n");
      }
      return lhs | (uint64_t(1) << bitshift);
    });
 
 
 
  }

TRGCDCModule()

TRGCDCModule

(

)

Constructor.

Definition at line 47 of file TRGCDCModule.cc.

    : Module::Module(),
      _debugLevel(0),
      _configFilename("TRGCDCConfig.dat"),
      _innerTSLUTFilename("undefined"),
      _outerTSLUTFilename("undefined"),
      _rootTRGCDCFilename("undefined"),
      _rootFitter3DFilename("undefined"),
      _curlBackStop(0),
      _simulationMode(1),
      _fastSimulationMode(0),
      _firmwareSimulationMode(0),
      _firmwareSimulationStart(0),
      _firmwareSimulationStop(32 * 32 - 1),
      _fileTRGCDC(0),
      _perfect2DFinder(false),
      _perfect3DFinder(false),
      _hFinderPeakMin(5),
      _hFinderMappingFilePlus("HoughMappingPlus20160223.dat"),
      _hFinderMappingFileMinus("HoughMappingMinus20160223.dat"),
      _hFinderDoit(2),
      _fLogicLUTTSF(1),
      _fLRLUT(1),
      _fFitter3Dsmclr(0),
      _fFitter3Ds2DFit(1),
      _fFitter3Ds2DFitDrift(0),
      _wireHitInefficiency(0.),
      _fileTSF(0),
      _fileETF(0),
      _fverETF(0),
      _fprintFirmETF(0),
      _fileHough3D(0),
      _finder3DMode(2),
      _fileFitter3D(0),
      _fXtSimpleFitter3D(0),
      _trgCDCDataInputMode(0),
      _cdc(0),
      _sa(0)
  {
 
    string desc = "TRGCDCModule(" + version() + ")";
    setDescription(desc);
    setPropertyFlags(c_ParallelProcessingCertified);
 
    vector<string> emptyvector;
 
    addParam("DebugLevel", _debugLevel, "TRGCDC debug level", _debugLevel);
    addParam("ConfigFile",
             _configFilename,
             "The filename of CDC trigger config file",
             _configFilename);
    addParam("ReturnValueModuleNames", _returnValueModuleNames,
             "Chooses modules for return value. Can be a list of [TSF, ETF, 2DFind, 2DFit, 3DFind, 3DFit]. If empty list then returnValue will be for all modules.",
             emptyvector);
    addParam("InnerTSLUTFile",
             _innerTSLUTFilename,
             "The filename of LUT for inner-most track segments",
             _innerTSLUTFilename);
    addParam("OuterTSLUTFile",
             _outerTSLUTFilename,
             "The filename of LUT for outer track segments",
             _outerTSLUTFilename);
    addParam("RootTRGCDCFile",
             _rootTRGCDCFilename,
             "The filename of root file for TRGCDC",
             string("TRGCDC.root"));
    addParam("RootFitter3DFile",
             _rootFitter3DFilename,
             "The filename of root file for Fitter3D",
             string("Fitter3D.root"));
    addParam("CurlBackStop",
             _curlBackStop,
             "Curl back stop parameter",
             _curlBackStop);
    addParam("SimulationMode",
             _simulationMode,
             "TRGCDC simulation switch",
             _simulationMode);
    addParam("FastSimulationMode",
             _fastSimulationMode,
             "TRGCDC fast simulation mode",
             _fastSimulationMode);
    addParam("FirmwareSimulationMode",
             _firmwareSimulationMode,
             "TRGCDC firmware simulation mode",
             _firmwareSimulationMode);
    addParam("FirmwareSimulationStart",
             _firmwareSimulationStart,
             "TRGCDC firmware simulation start clock of CDC FE",
             _firmwareSimulationStart);
    addParam("FirmwareSimulationStop",
             _firmwareSimulationStop,
             "TRGCDC firmware simulation stop clock of CDC FE",
             _firmwareSimulationStop);
    addParam("TRGCDCRootFile",
             _fileTRGCDC,
             "Flag for making TRGCDC.root",
             _fileTRGCDC);
    addParam("TSFLogicLUT",
             _fLogicLUTTSF,
             "Use logic or LUT for TSF",
             _fLogicLUTTSF);
    addParam("2DFinderPerfect",
             _perfect2DFinder,
             "2D finder perfect option",
             _perfect2DFinder);
    addParam("3DFinderPerfect",
             _perfect3DFinder,
             "3D finder perfect option",
             _perfect3DFinder);
    addParam("HoughFinderPeakMin",
             _hFinderPeakMin,
             "Hough finder min. peak height",
             _hFinderPeakMin);
    addParam("HoughFinderMappingFilePlus",
             _hFinderMappingFilePlus,
             "Hough finder cell mapping file for charge plus",
             _hFinderMappingFilePlus);
    addParam("HoughFinderMappingFileMinus",
             _hFinderMappingFileMinus,
             "Hough finder cell mapping file for charge minus",
             _hFinderMappingFileMinus);
    addParam("HoughFinderDoit",
             _hFinderDoit,
             "Hough finder doit() version",
             _hFinderDoit);
    addParam("Fitter3DsLRLUT",
             _fLRLUT,
             "Using L/R LUT in fitter3D",
             _fLRLUT);
    addParam("Fitter3DsMcLR",
             _fFitter3Dsmclr,
             "Using MC L/R information for Fitter3D",
             _fFitter3Dsmclr);
    addParam("Fitter3Ds2DFit",
             _fFitter3Ds2DFit,
             "Using 1: 2D fit or 0: Hough Finder results for Fitter3D",
             _fFitter3Ds2DFit);
    addParam("Fitter3Ds2DFitDrift",
             _fFitter3Ds2DFitDrift,
             "Using 1: drift 2D fit or 0: wire 2D fit for Fitter3D",
             _fFitter3Ds2DFitDrift);
    addParam("wireHitInefficiency",
             _wireHitInefficiency,
             "wire hit inefficiency",
             _wireHitInefficiency);
    addParam("TSFRootFile",
             _fileTSF,
             "Flag for making TSF.root",
             _fileTSF);
    addParam("ETFRootFile",
             _fileETF,
             "Flag for making ETF.root",
             _fileETF);
    addParam("verETF",
             _fverETF,
             "Select ETF ver, 0:histogram, 1:old Version, 2:True T0",
             _fverETF);
    addParam("printFirmETF",
             _fprintFirmETF,
             "Flag for printing Firmware inputs for ETF",
             _fprintFirmETF);
    addParam("Hough3DRootFile",
             _fileHough3D,
             "Flag for making Hough3D.root",
             _fileHough3D);
    addParam("Finder3DMode",
             _finder3DMode,
             "Mode for finder3D",
             _finder3DMode);
    addParam("Fitter3DRootFile",
             _fileFitter3D,
             "Flag for making Fitter3D.root",
             _fileFitter3D);
    addParam("Fitter3DsXtSimple",
             _fXtSimpleFitter3D,
             "0: Use simple xt curve table. 1: use non-linear xt curve table.",
             _fXtSimpleFitter3D);
    addParam("TRGCDCDataInputMode",
             _trgCDCDataInputMode,
             "Flag for TRG CDC input mode",
             _trgCDCDataInputMode);
    addParam("inputCollection", m_inputCollectionName,
             "Name of the CDCHit StoreArray to be used as input.",
             string(""));
    addParam("2DfinderCollection", m_2DfinderCollectionName,
             "Name of the StoreArray holding the tracks made by the 2D finder.",
             string("Trg2DFinderTracks"));
    addParam("2DfitterCollection", m_2DfitterCollectionName,
             "Name of the StoreArray holding the tracks made by the 2D fitter.",
             string("Trg2DFitterTracks"));
    addParam("3DfitterCollection", m_3DfitterCollectionName,
             "Name of the StoreArray holding the tracks made by the 3D fitter.",
             string("Trg3DFitterTracks"));
 
 
    if (TRGDebug::level())
      cout << "TRGCDCModule ... created" << endl;
  }

TRGCDCTSStreamModule()

TRGCDCTSStreamModule

(

)

Constructor.

Definition at line 42 of file TRGCDCTSStreamModule.cc.

    : Module::Module(),
      _debugLevel(0),
      _mode(0),
      _streamFilename("unknown"),
      _cdc(0),
      _out(0)
  {
 
    string desc = "TRGCDCTSStreamModule(" + version() + ")";
    setDescription(desc);
 
    addParam("DebugLevel",
             _debugLevel,
             "TRGCDCTSStream debug level",
             _debugLevel);
 
    addParam("Mode", _mode, "TRGCDCTSStream mode", _mode);
 
    addParam("OutputStreamFile",
             _streamFilename,
             "The filename of bit stream",
             _streamFilename);
 
    if (TRGDebug::level())
      cout << "TRGCDCTSStreamModule ... created" << endl;
  }

TRGECLBGTCHitModule()

TRGECLBGTCHitModule

(

)

Constructor.

Definition at line 41 of file TRGECLBGTCHitModule.cc.

    : Module::Module(),
      m_TCEnergyCut(0.005),
      m_TCTimingCutLow(-1500),
      m_TCTimingCutHigh(1500),
      m_debugLevel(0)
  {
    string desc = "TRGECLBGTCHitModule(" + version() + ")";
    setDescription(desc);
    setPropertyFlags(c_ParallelProcessingCertified);
 
    addParam("TCEnergyCut",
             m_TCEnergyCut,
             "TC energy cut: TC(E)<X (GeV), X=",
             m_TCEnergyCut);
    addParam("TCTimingCutLow",
             m_TCTimingCutLow,
             "TC Timing cut low: TC(T)>X (ns), X=",
             m_TCTimingCutLow);
    addParam("TCTimingCutHigh",
             m_TCTimingCutHigh,
             "TC Timing cut high: TC(T)<X (ns), X=",
             m_TCTimingCutHigh);
    addParam("DebugLevel",
             m_debugLevel,
             "TRGECL debug level",
             m_debugLevel);
 
    m_TCMap = new TrgEclMapping();
 
    B2DEBUG(100, "TRGECLBGTCHitModule ... created");
  }

TRGECLFAMModule()

TRGECLFAMModule

(

)

Constructor.

Definition at line 45 of file TRGECLFAMModule.cc.

    : Module::Module(),
      _debugLevel(0),
      _famMethod(1),
      _binTimeInterval(125),
      _waveform(0),
      _beambkgtag(0),
      _famana(0),
      _threshold(100.0),
      _FADC(1),
      _ConditionDB(0),
      m_SourceOfTC(3)
  {
 
    string desc = "TRGECLFAMModule(" + version() + ")";
    setDescription(desc);
    setPropertyFlags(c_ParallelProcessingCertified);
 
    addParam("DebugLevel", _debugLevel, "TRGECL debug level", _debugLevel);
    addParam("FAMFitMethod", _famMethod, "TRGECLFAM fit method", _famMethod);
    addParam("FAMBinTimeInterval", _binTimeInterval, "TRGECLFAM binTimeInterval",
             _binTimeInterval);
    addParam("TCWaveform", _waveform, "Output the TC waveform ",
             _waveform);
    addParam("FAMAnaTable", _famana, "Save FAM ana table ",
             _famana);
    addParam("BeamBkgTag", _beambkgtag, "Save beambackground tag in TRGECLHit table ",
             _beambkgtag);
    addParam("TCThreshold", _threshold, "Set FAM TC threshold ",
             _threshold);
    addParam("ShapingFunction", _FADC, "Set function of shaper ",  _FADC);
    addParam("ConditionDB", _ConditionDB, "Use conditionDB ",  _ConditionDB);
    addParam("SourceOfTC", m_SourceOfTC,
             "Select source of TC data(1:=ECLHit or 2:=ECLSimHit or 3:=ECLHit+TRGECLBGTCHit)",
             m_SourceOfTC);
 
    if (_ConditionDB == 1) { //Use global tag
      m_FAMPara.addCallback(this, &TRGECLFAMModule::beginRun);
    }
    B2DEBUG(100, "TRGECLFAMModule ... created");
    Threshold.clear();
 
  }

TRGECLModule()

TRGECLModule

(

)

Constructor.

Definition at line 44 of file TRGECLModule.cc.

                            :
    Module::Module(), _debugLevel(0), _Bhabha(0),
    _Clustering(1), _ClusterLimit(6), _EventTiming(1),
    _TimeWindow(250.0), _OverlapWindow(125.0), _NofTopTC(3),
    _SelectEvent(1), _ConditionDB(0), _mumuThreshold(20),
    _n300MeVCluster(1), _ECLBurstThreshold(200)
  {
 
    string desc = "TRGECLModule(" + version() + ")";
    setDescription(desc);
    setPropertyFlags(c_ParallelProcessingCertified);
 
    addParam("DebugLevel",
             _debugLevel,
             "TRGECL debug level",
             _debugLevel);
    addParam("Bhabha",
             _Bhabha,
             "TRGECL Bhabha method  0 : Belle I, 1 :belle II(defult)",
             _Bhabha);
    addParam("Clustering",
             _Clustering,
             "TRGECL Clustering method  0 : use only ICN, 1 : ICN + Energy(Defult)",
             _Clustering);
    addParam("ClusterLimit",
             _ClusterLimit,
             "The Limit number of cluster (Defult:6)",
             _ClusterLimit);
    addParam("EventTiming",
             _EventTiming,
             "TRGECL EventTiming method  0 : Belle I, 1 : Energetic TC, 2 : Energy Weighted timing (defult)",
             _EventTiming);
    addParam("NofTopTC",
             _NofTopTC,
             "TRGECL # of considered TC in energy weighted Timing method(Only work for EvenTiming Method 2)",
             _NofTopTC);
    addParam("TimeWindow",
             _TimeWindow,
             "TRGECL Trigger decision Time Window",
             _TimeWindow);
    addParam("OverlapWindow",
             _OverlapWindow,
             "TRGECL Trigger decision Time Window",
             _OverlapWindow);
    addParam("EventSelect",
             _SelectEvent,
             "TRGECL Select one trigger window for logic study",
             _SelectEvent);
    addParam("ConditionDB",
             _ConditionDB,
             "Flag to use Condition Database(Defult 0)",
             _ConditionDB);
    addParam("3DBhabhaVetoInTrackThetaRegion",
             m_3DBhabhaVetoInTrackThetaRegion,
             "set 3DBhabhaVetoInTrackThetaRegion (low TCID and high TCID)",
             m_3DBhabhaVetoInTrackThetaRegion);
    addParam("EventTimingQualityThresholds",
             m_EventTimingQualityThresholds,
             "set EventTimingQualityThresholds (low energy(GeV) and high energy(GeV))",
             m_EventTimingQualityThresholds);
    addParam("3DBhabhaAddAngleCut",
             m_3DBhabhaAddAngleCut,
             "set Angle Cut of Bhbaha Additional veto for hie (dphi low and high, theta sum low and high(degree)",
             m_3DBhabhaAddAngleCut);
    addParam("Taub2bAngleCut",
             m_taub2bAngleCut,
             "set Taub2bAngleCut(dphi low and high, theta sum low and high(degree)",
             m_taub2bAngleCut);
    addParam("Taub2bEtotCut",
             m_taub2bEtotCut,
             "set Taub2b total energy cut in TC ThetaID 1-17(GeV)",
             m_taub2bEtotCut);
    addParam("Taub2bClusterECut1",
             m_taub2bClusterECut1,
             "set Taub2b 1 cluster energy Cut (GeV) for CL1 : ",
             m_taub2bClusterECut1);
    addParam("Taub2bClusterECut2",
             m_taub2bClusterECut2,
             "set Taub2b 1 cluster energy Cut (GeV) for CL2 : ",
             m_taub2bClusterECut2);
 
    addParam("Taub2b2AngleCut",
             m_taub2b2AngleCut,
             "set Taub2b2AngleCut(dphi low and high, theta sum low and high(degree))",
             m_taub2b2AngleCut);
    addParam("Taub2b2EtotCut",
             m_taub2b2EtotCut,
             "set Taub2b2 total energy cut in TC ThetaID 1-17(GeV)",
             m_taub2b2EtotCut);
    addParam("Taub2b2CLEEndcapCut",
             m_taub2b2CLEEndcapCut,
             "set Taub2b2 energy Cut (GeV) for cluseter in endcap (GeV) : ",
             m_taub2b2CLEEndcapCut);
    addParam("Taub2b2CLECut",
             m_taub2b2CLECut,
             "set Taub2b2 cluster energy Cut (GeV) : ",
             m_taub2b2CLECut);
 
    addParam("Taub2b3AngleCut",
             m_taub2b3AngleCut,
             "set Taub2b3AngleCut(dphi low and high, theta sum low and high in cms (degree))",
             m_taub2b3AngleCut);
    addParam("Taub2b3EtotCut",
             m_taub2b3EtotCut,
             "set Taub2b3 total energy cut in lab in ThetaID 1-17(GeV)",
             m_taub2b3EtotCut);
    addParam("Taub2b3CLEb2bCut",
             m_taub2b3CLEb2bCut,
             "set Taub2b3 cluster energy Cut in lab for one of b2b clusters (GeV) : ",
             m_taub2b3CLEb2bCut);
    addParam("Taub2b3CLELowCut",
             m_taub2b3CLELowCut,
             "set Taub2b3 cluster energy low Cut in lab for all clusters (GeV) : ",
             m_taub2b3CLELowCut);
    addParam("Taub2b3CLEHighCut",
             m_taub2b3CLEHighCut,
             "set Taub2b3 cluster energy high Cut in lab for all clusters (GeV) : ",
             m_taub2b3CLEHighCut);
 
    _2DBhabhaThresholdFWD.clear();
    _2DBhabhaThresholdBWD.clear();
    _3DBhabhaSelectionThreshold.clear();
    _3DBhabhaVetoThreshold.clear();
    _3DBhabhaSelectionAngle.clear();
    _3DBhabhaVetoAngle.clear();
    _mumuAngle.clear();
    m_3DBhabhaAddAngleCut.clear();
    m_taub2bAngleCut.clear();
    m_taub2b2AngleCut.clear();
    m_taub2b3AngleCut.clear();
 
    _TotalEnergy.clear();
    _LowMultiThreshold.clear();
 
    _TotalEnergy = {5, 10, 30}; // /100 MeV
    _2DBhabhaThresholdFWD = {40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 30, 35}; // /100 MeV
    _2DBhabhaThresholdBWD  = {25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 30, 30}; // /100 MeV
    _3DBhabhaVetoThreshold = {30, 45}; //  /100 MeV
    _3DBhabhaSelectionThreshold = {20, 40}; //  /100 MeV
    _3DBhabhaVetoAngle = {160, 200, 165, 190}; //  /100 MeV
    _3DBhabhaSelectionAngle = {140, 220, 160, 200}; //  /100 MeV
    _mumuThreshold = 20; //100 MeV
    _mumuAngle = {160, 200, 165, 190}; //  degree
    m_3DBhabhaAddAngleCut   = {150, 210, 160, 200}; // degree
    _LowMultiThreshold = {10, 20, 25, 30}; // degree
    _n300MeVCluster = 1;
    _ECLBurstThreshold = 20;
    m_3DBhabhaVetoInTrackThetaRegion = {3, 15};
    m_EventTimingQualityThresholds = {0.5, 10}; // GeV
    // taub2b
    m_taub2bAngleCut = {110, 250, 130, 230}; // degree
    m_taub2bEtotCut = 7; // GeV
    m_taub2bClusterECut1 = 1.9; // GeV
    m_taub2bClusterECut2 = 999; // GeV
    // taub2b2
    m_taub2b2AngleCut = {120, 240, 140, 220}; // degree
    m_taub2b2EtotCut = 7.0; // GeV
    m_taub2b2CLEEndcapCut = 3.0; // GeV
    m_taub2b2CLECut = 0.162; // GeV
    // taub2b3
    m_taub2b3AngleCut = {120, 240, 140, 220}; // degree
    m_taub2b3EtotCut = 7.0; // GeV
    m_taub2b3CLEb2bCut  = 0.14; // GeV
    m_taub2b3CLELowCut  = 0.12; // GeV
    m_taub2b3CLEHighCut = 4.50; // GeV
 
    if (TRGDebug::level()) {
      std::cout << "TRGECLModule ... created" << std::endl;
    }
  }

TRGECLQAMModule()

TRGECLQAMModule

(

)

Costructor.

Definition at line 56 of file TRGECLQAMModule.cc.

    : Module::Module(),
      m_nevent(1), m_outputfile("")
  {
 
    string desc = "TRGECLQAMModule(" + version() + ")";
    setDescription(desc);
    setPropertyFlags(c_ParallelProcessingCertified);
 
    addParam("outputFileName", m_outputfile, "TRGECL QAM file", m_outputfile);
 
    m_nRun   = 0;
    m_nevent    = 1;
 
    //Problem TC Check
    m_FWD = 0;
    m_BAR = 0;
    m_BWD = 0;
    m_ALL = 0;
    memset(TCID, 0, sizeof(TCID));
 
 
    //Total Energy
 
    m_etot = 0;
    m_etot_mean = 0;
    m_etot_error = 0;
    m_etot_sigma = 0;
 
    //  std::vector<int> etot;
    etot.clear();
 
    //Caltime
 
    m_caltime_mean = 0;
    m_caltime_error = 0;
    m_caltime_sigma = 0;
    caltime.clear();
 
    //Cluster E
 
    clusterE = 0;
    cluster.clear();
 
    B2DEBUG(100, "TRGECLQAMModule ... created");
  }

TRGECLRawdataAnalysisModule()

TRGECLRawdataAnalysisModule

(

)

Constructor.

Definition at line 37 of file TRGECLRawdataAnalysisModule.cc.

    : Module::Module(),
      _debugLevel(0),  _Clustering(1)
  {
 
    string desc = "TRGECLRawdataAnalysisModule(" + version() + ")";
    setDescription(desc);
    //   setPropertyFlags(c_ParallelProcessingCertified | c_InitializeInProcess);
 
    addParam("DebugLevel", _debugLevel, "TRGECL debug level", _debugLevel);
    addParam("Clustering", _Clustering, "TRGECL Clustering method  0 : use only ICN, 1 : ICN + Energy(Defult)", _Clustering);
 
 
 
    if (TRGDebug::level()) {
      std::cout << "TRGECLRawdataAnalysisModule ... created" << std::endl;
    }
  }

TRGGDLModule()

TRGGDLModule

(

)

Constructor.

Definition at line 37 of file TRGGDLModule.cc.

                             : HistoModule(),
//TRGGDLModule::TRGGDLModule()
//  : Module::Module(),
    _debugLevel(0),
    _configFilename("TRGGDLConfig.dat"),
    _simulationMode(1),
    _fastSimulationMode(0),
    _firmwareSimulationMode(0),
    _Phase("Phase2"),
    _algFilePath("ftd.alg")
  {
 
    string desc = "TRGGDLModule(" + version() + ")";
    setDescription(desc);
    setPropertyFlags(c_ParallelProcessingCertified);
 
    addParam("debugLevel", _debugLevel, "Debug Level", _debugLevel);
    addParam("Belle2Phase", _Phase, "Phase2 or Phase3", _Phase);
    addParam("ConfigFile",
             _configFilename,
             "The filename of CDC trigger config file",
             _configFilename);
    addParam("SimulationMode",
             _simulationMode,
             "TRGGDL simulation switch",
             _simulationMode);
    addParam("FastSimulationMode",
             _fastSimulationMode,
             "TRGGDL fast simulation mode",
             _fastSimulationMode);
    addParam("FirmwareSimulationMode",
             _firmwareSimulationMode,
             "TRGGDL firmware simulation mode",
             _firmwareSimulationMode);
 
    addParam("algFromDB",
             _algFromDB,
             "Set false when alg is taken from local file.",
             true);
 
    addParam("algFilePath",
             _algFilePath,
             ".alg file path",
             _algFilePath);
 
    addParam("timquality_threshold_sfin",
             _timquality_threshold_sfin,
             "Threshold to determine timing quality flag with MC truth: superfine",
             -1.0);
 
    addParam("timquality_threshold_fine",
             _timquality_threshold_fine,
             "Threshold to determine timing quality flag with MC truth: fine",
             10.0);
 
    addParam("simulateT0jitter",
             m_simulateT0Jitter,
             "if True, L1 jitter is simulated by EventT0Generator.",
             false);
 
    B2DEBUG(100, "TRGGDLModule ... created");
  }

TRGGRLModule()

TRGGRLModule

(

)

Constructor.

Definition at line 36 of file TRGGRLModule.cc.

    : Module::Module(),
      _debugLevel(0),
      _configFilename("TRGGDLConfig.dat"),
      _simulationMode(2),
      _fastSimulationMode(0),
      _firmwareSimulationMode(0)
  {
 
    string desc = "TRGGRLModule(" + version() + ")";
    setDescription(desc);
    setPropertyFlags(c_ParallelProcessingCertified);
 
    addParam("DebugLevel", _debugLevel, "TRGGRL debug level", _debugLevel);
    addParam("ConfigFile",
             _configFilename,
             "The filename of CDC trigger config file",
             _configFilename);
    addParam("SimulationMode",
             _simulationMode,
             "TRGGRL simulation switch",
             _simulationMode);
    addParam("FastSimulationMode",
             _fastSimulationMode,
             "TRGGRL fast simulation mode",
             _fastSimulationMode);
    addParam("FirmwareSimulationMode",
             _firmwareSimulationMode,
             "TRGGRL firmware simulation mode",
             _firmwareSimulationMode);
  }

version() [1/9]

string version

(

void

)

const

version

Definition at line 49 of file TRGECLQAMModule.cc.

  {
    return string("TRGECLQAMModule 1.00");
  }

version() [2/9]

string version

(

void

)

const

returns version of TRGCDCModule.

Definition at line 42 of file TRGCDCModule.cc.

  {
    return string("TRGCDCModule 5.07");
  }

version() [3/9]

string version

(

void

)

const

returns version of TRGCDCTSStreamModule.

Definition at line 37 of file TRGCDCTSStreamModule.cc.

  {
    return string("TRGCDCTSStreamModule 0.00");
  }

version() [4/9]

string version

(

void

)

const

returns version of TRGECLFAMModule.

Definition at line 38 of file TRGECLFAMModule.cc.

  {
    return string("TRGECLFAMModule 1.00");
  }

version() [5/9]

string version

(

void

)

const

returns version of TRGECLModule.

Definition at line 37 of file TRGECLModule.cc.

  {
    return string("TRGECLModule 1.00");
  }

version() [6/9]

string version

(

void

)

const

returns version of TRGECLBGTCHitModule.

Definition at line 34 of file TRGECLBGTCHitModule.cc.

  {
    return string("TRGECLBGTCHitModule 1.00");
  }

version() [7/9]

string version

(

void

)

const

returns version of TRGECLModule.

Definition at line 30 of file TRGECLRawdataAnalysisModule.cc.

  {
    return string("TRGECLAnalysisModule 1.00");
  }

version() [8/9]

string version

(

void

)

const

returns version of TRGGDLModule.

Definition at line 32 of file TRGGDLModule.cc.

  {
    return string("TRGGDLModule 0.00");
  }

version() [9/9]

string version

(

void

)

const

returns version of TRGGRLModule.

Definition at line 31 of file TRGGRLModule.cc.

  {
    return string("TRGGRLModule 0.00");
  }

~TRGCDCModule()

~TRGCDCModule ( )

virtual

Destructor.

Definition at line 247 of file TRGCDCModule.cc.

  {
    if (_cdc)
      TRGCDC::getTRGCDC("good-bye");
 
    if (TRGDebug::level())
      cout << "TRGCDCModule ... destructed " << endl;
  }

~TRGCDCTSStreamModule()

~TRGCDCTSStreamModule ( )

virtual

Destructor.

Definition at line 70 of file TRGCDCTSStreamModule.cc.

  {
 
    if (_cdc)
      TRGCDC::getTRGCDC("good-bye");
 
    if (TRGDebug::level())
      cout << "TRGCDCTSStreamModule ... destructed " << endl;
  }

~TRGECLBGTCHitModule()

~TRGECLBGTCHitModule ( )

virtual

Destructor.

Definition at line 76 of file TRGECLBGTCHitModule.cc.

  {
    B2DEBUG(100, "TRGECLBGTCHitModule ... destructed ");
  }

~TRGECLFAMModule()

~TRGECLFAMModule ( )

virtual

Destructor.

Definition at line 91 of file TRGECLFAMModule.cc.

  {
 
    B2DEBUG(100, "TRGECLFAMModule ... destructed ");
 
  }

~TRGECLModule()

~TRGECLModule ( )

virtual

Destructor

Definition at line 217 of file TRGECLModule.cc.

  {
    if (TRGDebug::level()) {
      std::cout << "TRGECLModule ... destructed " << std::endl;
 
    }
  }

~TRGECLQAMModule()

~TRGECLQAMModule ( )

virtual

Destrunctor.

Definition at line 105 of file TRGECLQAMModule.cc.

  {
 
    B2DEBUG(100, "TRGECLQAMModule ... destructed ");
  }

~TRGECLRawdataAnalysisModule()

~TRGECLRawdataAnalysisModule ( )

virtual

Destructor

Definition at line 58 of file TRGECLRawdataAnalysisModule.cc.

  {
 
    if (TRGDebug::level()) {
      std::cout << "TRGECLRawdataAnalysisModule ... destructed " << std::endl;
 
    }
  }

~TRGGDLModule()

~TRGGDLModule ( )

virtual

Destructor.

Definition at line 100 of file TRGGDLModule.cc.

  {
 
    if (_gdl)
      B2DEBUG(100, "good-bye");
 
    B2DEBUG(100,  "TRGGDLModule ... destructed ");
  }

~TRGGRLModule()

~TRGGRLModule ( )

virtual

Destructor.

Definition at line 68 of file TRGGRLModule.cc.

  {
 
    if (_grl)
      TRGGRL::getTRGGRL("good-bye");
 
  }

Variable Documentation

_gdl

TRGGDL * _gdl = 0

staticprivate

A pointer to a TRGGDL;.

Definition at line 113 of file TRGGDLModule.h.

_grl

TRGGRL * _grl = 0

staticprivate

A pointer to a TRGGRL;.

Definition at line 79 of file TRGGRLModule.h.

intNaN

const int intNaN = std::numeric_limits<int>::quiet_NaN()

static

Definition at line 53 of file trgrawdataModule.h.

Linear_fit_of_Hits

constexpr Linear_fit_of_Hits_t Linear_fit_of_Hits

constexpr

Definition at line 140 of file klm_trig_linear_fit.cc.

mergerWidth

constexpr int mergerWidth = 256

staticconstexpr

Merger data width.

Definition at line 45 of file CDCTriggerUnpackerModule.h.

nAllMergers

constexpr int nAllMergers = 146

staticconstexpr

Number of Mergers.

Definition at line 47 of file CDCTriggerUnpackerModule.h.

nFinesse

constexpr int nFinesse = 48

staticconstexpr

Number of FINESSE in a PCIe40.

Definition at line 51 of file CDCTriggerUnpackerModule.h.

wordWidth

constexpr int wordWidth = 32

staticconstexpr

width of a single word in the raw int buffer

Definition at line 49 of file CDCTriggerUnpackerModule.h.