Belle II Software development
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Static Private Member Functions | Private Attributes | Static Private Attributes | List of all members
B2BIIConvertMdstModule Class Reference
Packages » b2bii | Modules » b2bii modules

Module converts Belle MDST objects (Panther records) to Belle II MDST objects. More...

#include <B2BIIConvertMdstModule.h>

Inheritance diagram for B2BIIConvertMdstModule:
Module PathElement

Public Types

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

Public Member Functions

 B2BIIConvertMdstModule ()
 Constructor.
 
virtual ~B2BIIConvertMdstModule () override
 Destructor.
 
virtual void initialize () override
 Initialize the module.
 
virtual void beginRun () override
 Module functions to be called from event process.
 
virtual void event () override
 Called for each event.
 
virtual void endRun () override
 Called when the current run is finished.
 
virtual void terminate () override
 Terminates the module.
 
virtual std::vector< std::string > getFileNames (bool outputFiles)
 Return a list of output filenames for this modules.
 
const std::string & getName () const
 Returns the name of the module.
 
const std::string & getType () const
 Returns the type of the module (i.e.
 
const std::string & getPackage () const
 Returns the package this module is in.
 
const std::string & getDescription () const
 Returns the description of the module.
 
void setName (const std::string &name)
 Set the name of the module.
 
void setPropertyFlags (unsigned int propertyFlags)
 Sets the flags for the module properties.
 
LogConfig & getLogConfig ()
 Returns the log system configuration.
 
void setLogConfig (const LogConfig &logConfig)
 Set the log system configuration.
 
void setLogLevel (int logLevel)
 Configure the log level.
 
void setDebugLevel (int debugLevel)
 Configure the debug messaging level.
 
void setAbortLevel (int abortLevel)
 Configure the abort log level.
 
void setLogInfo (int logLevel, unsigned int logInfo)
 Configure the printed log information for the given level.
 
void if_value (const std::string &expression, const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
 Add a condition to the module.
 
void if_false (const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
 A simplified version to add a condition to the module.
 
void if_true (const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
 A simplified version to set the condition of the module.
 
bool hasCondition () const
 Returns true if at least one condition was set for the module.
 
const ModuleCondition * getCondition () const
 Return a pointer to the first condition (or nullptr, if none was set)
 
const std::vector< ModuleCondition > & getAllConditions () const
 Return all set conditions for this module.
 
bool evalCondition () const
 If at least one condition was set, it is evaluated and true returned if at least one condition returns true.
 
std::shared_ptr< Path > getConditionPath () const
 Returns the path of the last true condition (if there is at least one, else reaturn a null pointer).
 
Module::EAfterConditionPath getAfterConditionPath () const
 What to do after the conditional path is finished.
 
std::vector< std::shared_ptr< Path > > getAllConditionPaths () const
 Return all condition paths currently set (no matter if the condition is true or not).
 
bool hasProperties (unsigned int propertyFlags) const
 Returns true if all specified property flags are available in this module.
 
bool hasUnsetForcedParams () const
 Returns true and prints error message if the module has unset parameters which the user has to set in the steering file.
 
const ModuleParamList & getParamList () const
 Return module param list.
 
template<typename T >
ModuleParam< T > & getParam (const std::string &name) const
 Returns a reference to a parameter.
 
bool hasReturnValue () const
 Return true if this module has a valid return value set.
 
int getReturnValue () const
 Return the return value set by this module.
 
std::shared_ptr< PathElement > clone () const override
 Create an independent copy of this module.
 
std::shared_ptr< boost::python::list > getParamInfoListPython () const
 Returns a python list of all parameters.
 

Static Public Member Functions

static void exposePythonAPI ()
 Exposes methods of the Module class to Python.
 

Protected Member Functions

virtual void def_initialize ()
 Wrappers to make the methods without "def_" prefix callable from Python.
 
virtual void def_beginRun ()
 Wrapper method for the virtual function beginRun() that has the implementation to be used in a call from Python.
 
virtual void def_event ()
 Wrapper method for the virtual function event() that has the implementation to be used in a call from Python.
 
virtual void def_endRun ()
 This method can receive that the current run ends as a call from the Python side.
 
virtual void def_terminate ()
 Wrapper method for the virtual function terminate() that has the implementation to be used in a call from Python.
 
void setDescription (const std::string &description)
 Sets the description of the module.
 
void setType (const std::string &type)
 Set the module type.
 
template<typename T >
void addParam (const std::string &name, T &paramVariable, const std::string &description, const T &defaultValue)
 Adds a new parameter to the module.
 
template<typename T >
void addParam (const std::string &name, T &paramVariable, const std::string &description)
 Adds a new enforced parameter to the module.
 
void setReturnValue (int value)
 Sets the return value for this module as integer.
 
void setReturnValue (bool value)
 Sets the return value for this module as bool.
 
void setParamList (const ModuleParamList &params)
 Replace existing parameter list.
 

Private Types

enum  MCMatchingMode {
  c_Direct ,
  c_GeneratorLevel
}
 MC matching mode. More...
 

Private Member Functions

void initializeDataStore ()
 Initializes Belle II DataStore.
 
void convertRecTrgTable ()
 Reads and converts m_final from rectrg_summary3.
 
void convertEvtclsTable ()
 Reads and converts all entries of evtcls Panther table.
 
void convertGenHepEvtTable ()
 Reads and converts all entries of Gen_hepevt Panther table to MCParticle dataobjects and adds them to StoreArray<MCParticle>.
 
void convertMdstECLTable ()
 Reads and converts all entries of Mdst_ecl(_aux) Panther table to ECLCluster dataobjects and adds them to StoreArray<ECLCluster>.
 
void convertMdstKLMTable ()
 Reads and converts all entries of Mdst_klm_cluster Panther table to KLMCluster dataobjects and adds them to StoreArray<KLMCluster>.
 
void convertMdstChargedTable ()
 Reads and converts all entries of Mdst_charged (Mdst_trk and Mdst_trk_fit) Panther table to Track (TrackFitResult) dataobjects and adds them to StoreArray<Track> (StoreArray<TrackFitResult>).
 
void convertMdstGammaTable ()
 Reads all entries of Mdst_Gamma Panther table, creates a particle list 'gamma:mdst' and adds them to StoreArray<Particles>.
 
void copyNbarFromGamma ()
 Copies Particles in 'gamma:mdst' with energy > 0.5 GeV to be anti-n0:mdst.
 
void convertMdstKLongTable ()
 Reads all entries of Mdst_Klong Panther table, creates a particle list 'K_L0:mdst' and adds them to StoreArray<Particles>.
 
void convertMdstPi0Table ()
 Reads all entries of Mdst_Pi0 Panther table, creates a particle list 'pi0:mdst' and adds them to StoreArray<Particles>.
 
void convertMdstVee2Table ()
 Reads and converts all entries of Mdst_vee2 Panther table to V0 dataobjects and adds them to StoreArray<V0>.
 
void convertBeamEnergy ()
 Stores beam parameters (energy, angles) in CollisionInvariantMass and CollisionBoostVector (currently in the DataStore).
 
void convertIPProfile (bool beginRun=false)
 Stores the IPProfiles in BeamSpot (currently in DataStore)
 
void convertGenHepevtObject (const Belle::Gen_hepevt &genHepevt, MCParticleGraph::GraphParticle *mcParticle)
 Converts Gen_hepevt record to MCParticleGraph object.
 
void convertMdstECLObject (const Belle::Mdst_ecl &ecl, const Belle::Mdst_ecl_aux &eclAux, ECLCluster *eclCluster)
 Converts Mdst_ecl(_aux) record to ECLCluster object.
 
double computeTrkMinDistanceBelle (ECLCluster *eclCluster)
 calculate the minimal distance between a cluster and a set of tracks on the ECL surface.
 
void convertMdstKLMObject (const Belle::Mdst_klm_cluster &klm, KLMCluster *klmCluster)
 Converts Mdst_klm_cluster record to KLMCluster object.
 
void convertMdstChargedObject (const Belle::Mdst_charged &belleTrack, Track *track)
 Converts Mdst_charged (Mdst_trk(_fit)) record to Track (TrackFitResult) object.
 
TrackFitResult createTrackFitResult (const CLHEP::HepLorentzVector &momentum, const HepPoint3D &position, const CLHEP::HepSymMatrix &error, const short int charge, const Const::ParticleType &pType, const float pValue, const uint64_t hitPatternCDCInitializer, const uint32_t hitPatternVXDInitializer, const uint16_t ndf)
 Creates TrackFitResult and fills it.
 
int getHelixParameters (const Belle::Mdst_trk_fit &trk_fit, const double mass, const HepPoint3D &newPivot, std::vector< float > &helixParams, CLHEP::HepSymMatrix &error5x5, CLHEP::HepLorentzVector &momentum, HepPoint3D &position, CLHEP::HepSymMatrix &error7x7, const double dPhi=0.0)
 Fills Helix parameters (converted to Belle II version), 5x5 error matrix, 4-momentum, position and 7x7 error matrix from Belle Helix stored in Mdst_trk_fit.
 
void convertHelix (Belle::Helix &helix, std::vector< float > &helixParams, CLHEP::HepSymMatrix &error5x5)
 Converts Belle's Helix parameters and it's covariance matrix to Belle II's version.
 
void convertHelix (const Belle::Mdst_trk_fit &trk_fit, const HepPoint3D &newPivot, std::vector< float > &helixParams, std::vector< float > &helixError)
 Converts Belle's Helix parameters and it's covariance matrix to Belle II's version.
 
void belleVeeDaughterToCartesian (const Belle::Mdst_vee2 &vee, const int charge, const Const::ParticleType &pType, CLHEP::HepLorentzVector &momentum, HepPoint3D &position, CLHEP::HepSymMatrix &error)
 Fills 4-momentum, position and 7x7 error matrix from Belle Vee daughter.
 
void belleVeeDaughterHelix (const Belle::Mdst_vee2 &vee, const int charge, std::vector< float > &helixParam, std::vector< float > &helixError)
 obtains the helix parameters of the vee daughters
 
void convertPIDData (const Belle::Mdst_charged &belleTrack, const Track *track)
 Get PID information for belleTrack and add it to PIDLikelihood (with relation from track).
 
void setLikelihoods (PIDLikelihood *pid, Const::EDetector det, double likelihoods[c_nHyp], bool discard_allzero=false)
 Add given Belle likelihoods (not log-likelihoods, in Belle hypothesis order) for given detector to pid.
 
double cdc_pid (const Belle::Mdst_charged &chg, int idp)
 Returns CDC likelihood for given hypothesis idp.
 
double atcPID (const PIDLikelihood *pid, int sigHyp, int bkgHyp)
 calculates atc_pid(3,1,5,sigHyp,bkgHyp).prob() from converted PIDLikelihood
 
void setTracksToECLClustersRelations ()
 Sets Track -> ECLCluster relations.
 
void setKLMClustersRelations ()
 Sets KLMCluster -> Track and ECLCluster relations.
 
int recoverMoreThan24bitIDHEP (int id)
 Helper function to recover falsely set idhep info in GenHepEvt list.
 
void testMCRelation (const Belle::Gen_hepevt &belleMC, const MCParticle *mcP, const std::string &objectName)
 Checks if the reconstructed object (Track, ECLCluster, ...) was matched to the same MCParticle.
 
std::list< ModulePtr > getModules () const override
 no submodules, return empty list
 
std::string getPathString () const override
 return the module name.
 
void setParamPython (const std::string &name, const boost::python::object &pyObj)
 Implements a method for setting boost::python objects.
 
void setParamPythonDict (const boost::python::dict &dictionary)
 Implements a method for reading the parameter values from a boost::python dictionary.
 

Static Private Member Functions

static double acc_pid (const Belle::Mdst_charged &chg, int idp)
 Returns ACC likelihood for given hypothesis idp.
 

Private Attributes

bool m_realData
 flag that tells whether given data sample is for real data or MC
 
bool m_convertBeamParameters
 Convert beam parameters or use information stored in Belle II database.
 
bool m_use6x6CovarianceMatrix4Tracks
 flag that tells which form of covariance matrix should be used in the conversion of charged tracks
 
std::string m_mcMatchingModeString
 MC matching mode.
 
MCMatchingMode m_mcMatchingMode
 C matching mode.
 
bool m_convertEvtcls
 Flag to switch on conversion of Evtcls table.
 
bool m_nisEnable
 Flag to switch on conversion of nisKsFinder info.
 
bool m_convertRecTrg
 Flag to switch on conversion of rectrg_summary3.
 
bool m_convertTrkExtra
 Flag to switch on conversion of first(last)_{x,y,z} of mdst_trk_fit.
 
bool m_convertNbar
 Flag to create anti-n0:mdst list from gamma:mdst.
 
double m_matchType2E9oE25Threshold
 E9/E25 threshold value clusters with a value above this threshold are classified as neutral even if tracks are matched to their connected region (matchType == 2)
 
int m_lastIPProfileBin { -1}
 variable to tell us which IPProfile bin was active last time we looked
 
Belle2::MCParticleGraph m_particleGraph
 MCParticle Graph to build Belle2 MC Particles.
 
std::map< int, int > genHepevtToMCParticle
 map of Gen_hepevt Panther IDs and corresponding MCParticle StoreArray indices
 
std::map< int, int > mdstEclToECLCluster
 map of Mdst_ecl Panther IDs and corresponding ECLCluster StoreArray indices
 
std::map< int, int > mdstGammaToParticle
 map of gamma Panther IDs and corresponding Particle StoreArray indices
 
std::map< int, int > mdstKlmToKLMCluster
 map of Mdst_klm Panther IDs and corresponding KLMCluster StoreArray indices
 
std::map< int, int > mdstKlongToParticle
 map of Klong Panther IDs and corresponding Particle StoreArray indices
 
StoreArray< MCParticle > m_mcParticles
 MC particles.
 
StoreArray< ECLCluster > m_eclClusters
 ECL clusters.
 
StoreArray< KLMCluster > m_klmClusters
 KLM clusters.
 
StoreArray< Track > m_tracks
 Tracks.
 
StoreArray< TrackFitResult > m_trackFitResults
 Track fir results.
 
StoreArray< V0 > m_v0s
 V0-particles.
 
StoreArray< Particle > m_particles
 Particles.
 
StoreArray< BelleTrkExtra > m_belleTrkExtra
 Belle CDC extra information.
 
StoreArray< PIDLikelihood > m_pidLikelihoods
 output PIDLikelihood array.
 
StoreObjPtr< EventExtraInfo > m_evtInfo
 Event Extra Info.
 
OptionalDBObjPtr< BeamSpot > m_beamSpotDB
 BeamSpot for IP.
 
BeamSpot m_beamSpot
 Interaction Point of the beam.
 
OptionalDBObjPtr< CollisionBoostVector > m_collisionBoostVectorDB
 CollisionBoostVector for boost vector.
 
CollisionBoostVector m_collisionBoostVector
 CollisionBoostVector for bosst vector of the beam.
 
OptionalDBObjPtr< CollisionInvariantMass > m_collisionInvMDB
 CollisionInvariantMass for Invariant Mass of Beam.
 
CollisionInvariantMass m_collisionInvM
 CollisionInvariantMass for the invariant mass of the beam.
 
OptionalDBObjPtr< CollisionAxisCMS > m_collisionAxisCMSDB
 CollisionAxisCMS.
 
CollisionAxisCMS m_collisionAxisCMS
 CollisionAxisCMS of the beam.
 
const int ERRMCONV [7] = {3, 4, 5, -1, 0, 1, 2}
 CONVERSION OF TRACK ERROR MATRIX ELEMENTS.
 
const double BFIELD = 1.5
 B filed in TESLA.
 
const double KAPPA2OMEGA = 1.5 * TMath::C() * 1E-11
 Conversion factor for Kappa -> Omega helix parameters.
 
std::string m_name
 The name of the module, saved as a string (user-modifiable)
 
std::string m_type
 The type of the module, saved as a string.
 
std::string m_package
 Package this module is found in (may be empty).
 
std::string m_description
 The description of the module.
 
unsigned int m_propertyFlags
 The properties of the module as bitwise or (with |) of EModulePropFlags.
 
LogConfig m_logConfig
 The log system configuration of the module.
 
ModuleParamList m_moduleParamList
 List storing and managing all parameter of the module.
 
bool m_hasReturnValue
 True, if the return value is set.
 
int m_returnValue
 The return value.
 
std::vector< ModuleCondition > m_conditions
 Module condition, only non-null if set.
 

Static Private Attributes

static const int c_nHyp = 5
 Number of Belle track hypotheses (see c_belleHyp_to_chargedStable).
 
static const Const::ChargedStable c_belleHyp_to_chargedStable [c_nHyp] = { Const::electron, Const::muon, Const::pion, Const::kaon, Const::proton }
 maps Belle hypotheses to Const::ChargedStable (from http://belle.kek.jp/secured/wiki/doku.php?id=software:atc_pid).
 

Detailed Description

Module converts Belle MDST objects (Panther records) to Belle II MDST objects.

The module performs conversion of the following Belle Panther records:

o) Gen_hepevt -> MCParticle

o) Mdst_ecl and Mdst_ecl_aux -> ECLCluster

The following relations are created:

o) ECLCluster -> MCParticle

Definition at line 92 of file B2BIIConvertMdstModule.h.

Member Typedef Documentation

◆ EAfterConditionPath

typedef ModuleCondition::EAfterConditionPath EAfterConditionPath
inherited

Forward the EAfterConditionPath definition from the ModuleCondition.

Definition at line 88 of file Module.h.

Member Enumeration Documentation

◆ EModulePropFlags

enum EModulePropFlags
inherited

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

Enumerator
c_Input 

This module is an input module (reads data).

c_Output 

This module is an output module (writes data).

c_ParallelProcessingCertified 

This module can be run in parallel processing mode safely (All I/O must be done through the data store, in particular, the module must not write any files.)

c_HistogramManager 

This module is used to manage histograms accumulated by other modules.

c_InternalSerializer 

This module is an internal serializer/deserializer for parallel processing.

c_TerminateInAllProcesses 

When using parallel processing, call this module's terminate() function in all processes().

This will also ensure that there is exactly one process (single-core if no parallel modules found) or at least one input, one main and one output process.

c_DontCollectStatistics 

No statistics is collected for this module.

Definition at line 77 of file Module.h.

77 {
78 c_Input = 1,
79 c_Output = 2,
80 c_ParallelProcessingCertified = 4,
81 c_HistogramManager = 8,
82 c_InternalSerializer = 16,
83 c_TerminateInAllProcesses = 32,
84 c_DontCollectStatistics = 64,
85 };
Belle2::Module::c_HistogramManager
@ c_HistogramManager
This module is used to manage histograms accumulated by other modules.
Definition: Module.h:81
Belle2::Module::c_Input
@ c_Input
This module is an input module (reads data).
Definition: Module.h:78
Belle2::Module::c_DontCollectStatistics
@ c_DontCollectStatistics
No statistics is collected for this module.
Definition: Module.h:84
Belle2::Module::c_ParallelProcessingCertified
@ c_ParallelProcessingCertified
This module can be run in parallel processing mode safely (All I/O must be done through the data stor...
Definition: Module.h:80
Belle2::Module::c_InternalSerializer
@ c_InternalSerializer
This module is an internal serializer/deserializer for parallel processing.
Definition: Module.h:82
Belle2::Module::c_Output
@ c_Output
This module is an output module (writes data).
Definition: Module.h:79
Belle2::Module::c_TerminateInAllProcesses
@ c_TerminateInAllProcesses
When using parallel processing, call this module's terminate() function in all processes().
Definition: Module.h:83

◆ MCMatchingMode

enum MCMatchingMode
private

MC matching mode.

Enumerator
c_Direct 

Direct matching.

c_GeneratorLevel 

Match to generator-level particles.

Definition at line 97 of file B2BIIConvertMdstModule.h.

97 {
98
100 c_Direct,
101
103 c_GeneratorLevel,
104
105 };
Belle2::B2BIIConvertMdstModule::c_GeneratorLevel
@ c_GeneratorLevel
Match to generator-level particles.
Definition: B2BIIConvertMdstModule.h:103

Constructor & Destructor Documentation

◆ B2BIIConvertMdstModule()

B2BIIConvertMdstModule ( )

Constructor.

Definition at line 164 of file B2BIIConvertMdstModule.cc.

164 : Module(),
165 m_mcMatchingMode(c_Direct)
166{
167 //Set module properties
168 setDescription("Converts Belle mDST objects (Panther tables and records) to Belle II mDST objects.");
169
170 addParam("convertBeamParameters", m_convertBeamParameters,
171 "Convert beam parameters or use information stored in "
172 "Belle II database.", true);
173 addParam("use6x6CovarianceMatrix4Tracks", m_use6x6CovarianceMatrix4Tracks,
174 "Use 6x6 (position, momentum) covariance matrix for charged tracks instead of 5x5 (helix parameters) covariance matrix", false);
175 addParam("mcMatchingMode", m_mcMatchingModeString,
176 "MC matching mode: 'Direct', or 'GeneratorLevel'",
177 std::string("Direct"));
178 addParam("matchType2E9oE25Threshold", m_matchType2E9oE25Threshold,
179 "clusters with a E9/E25 value above this threshold are classified as neutral even if tracks are matched to their connected region (matchType == 2)",
180 -1.1);
181
182 addParam("convertEvtcls", m_convertEvtcls, "Flag to switch on conversion of Mdst_evtcls", true);
183 addParam("nisKsInfo", m_nisEnable, "Flag to switch on conversion of nisKsFinder info", true);
184 addParam("RecTrg", m_convertRecTrg, "Flag to switch on conversion of rectrg_summary3", false);
185 addParam("TrkExtra", m_convertTrkExtra, " Flag to switch on conversion of first_x,y,z and last_x,y,z from Mdst_trk_fit", true);
186 addParam("convertNbar", m_convertNbar, " Flag to switch on conversion of nbar:mdst (copy from gamma:mdst)", false);
187
188 m_realData = false;
189
190 B2DEBUG(1, "B2BIIConvertMdst: Constructor done.");
191}
Belle2::B2BIIConvertMdstModule::m_convertEvtcls
bool m_convertEvtcls
Flag to switch on conversion of Evtcls table.
Definition: B2BIIConvertMdstModule.h:151
Belle2::B2BIIConvertMdstModule::m_nisEnable
bool m_nisEnable
Flag to switch on conversion of nisKsFinder info.
Definition: B2BIIConvertMdstModule.h:153
Belle2::B2BIIConvertMdstModule::m_convertNbar
bool m_convertNbar
Flag to create anti-n0:mdst list from gamma:mdst.
Definition: B2BIIConvertMdstModule.h:159
Belle2::B2BIIConvertMdstModule::m_convertBeamParameters
bool m_convertBeamParameters
Convert beam parameters or use information stored in Belle II database.
Definition: B2BIIConvertMdstModule.h:138
Belle2::B2BIIConvertMdstModule::m_realData
bool m_realData
flag that tells whether given data sample is for real data or MC
Definition: B2BIIConvertMdstModule.h:135
Belle2::B2BIIConvertMdstModule::m_mcMatchingModeString
std::string m_mcMatchingModeString
MC matching mode.
Definition: B2BIIConvertMdstModule.h:146
Belle2::B2BIIConvertMdstModule::m_convertRecTrg
bool m_convertRecTrg
Flag to switch on conversion of rectrg_summary3.
Definition: B2BIIConvertMdstModule.h:155
Belle2::B2BIIConvertMdstModule::m_matchType2E9oE25Threshold
double m_matchType2E9oE25Threshold
E9/E25 threshold value clusters with a value above this threshold are classified as neutral even if t...
Definition: B2BIIConvertMdstModule.h:165
Belle2::B2BIIConvertMdstModule::m_mcMatchingMode
MCMatchingMode m_mcMatchingMode
C matching mode.
Definition: B2BIIConvertMdstModule.h:149
Belle2::B2BIIConvertMdstModule::m_use6x6CovarianceMatrix4Tracks
bool m_use6x6CovarianceMatrix4Tracks
flag that tells which form of covariance matrix should be used in the conversion of charged tracks
Definition: B2BIIConvertMdstModule.h:143
Belle2::B2BIIConvertMdstModule::m_convertTrkExtra
bool m_convertTrkExtra
Flag to switch on conversion of first(last)_{x,y,z} of mdst_trk_fit.
Definition: B2BIIConvertMdstModule.h:157
Belle2::Module::setDescription
void setDescription(const std::string &description)
Sets the description of the module.
Definition: Module.cc:214
Belle2::Module::Module
Module()
Constructor.
Definition: Module.cc:30
Belle2::Module::addParam
void addParam(const std::string &name, T &paramVariable, const std::string &description, const T &defaultValue)
Adds a new parameter to the module.
Definition: Module.h:559

◆ ~B2BIIConvertMdstModule()

~B2BIIConvertMdstModule ( )
overridevirtual

Destructor.

Definition at line 194 of file B2BIIConvertMdstModule.cc.

195{
196}

Member Function Documentation

◆ acc_pid()

static double acc_pid ( const Belle::Mdst_charged &  chg,
int  idp 
)
staticprivate

Returns ACC likelihood for given hypothesis idp.

Copied from atc_pid::acc_pid().

◆ atcPID()

double atcPID ( const PIDLikelihood *  pid,
int  sigHyp,
int  bkgHyp 
)
private

calculates atc_pid(3,1,5,sigHyp,bkgHyp).prob() from converted PIDLikelihood

Definition at line 2303 of file B2BIIConvertMdstModule.cc.

2304{
2305 if (!pid) return 0.5;
2306
2307 // ACC = ARICH
2308 Const::PIDDetectorSet set = Const::ARICH;
2309 double acc_sig = exp(pid->getLogL(c_belleHyp_to_chargedStable[sigHyp], set));
2310 double acc_bkg = exp(pid->getLogL(c_belleHyp_to_chargedStable[bkgHyp], set));
2311 double acc = 0.5;
2312 if (acc_sig + acc_bkg > 0.0)
2313 acc = acc_sig / (acc_sig + acc_bkg);
2314
2315 // TOF = TOP
2316 set = Const::TOP;
2317 double tof_sig = exp(pid->getLogL(c_belleHyp_to_chargedStable[sigHyp], set));
2318 double tof_bkg = exp(pid->getLogL(c_belleHyp_to_chargedStable[bkgHyp], set));
2319 double tof = 0.5;
2320 double tof_all = tof_sig + tof_bkg;
2321 if (tof_all != 0) {
2322 tof = tof_sig / tof_all;
2323 if (tof < 0.001) tof = 0.001;
2324 if (tof > 0.999) tof = 0.999;
2325 }
2326
2327 // dE/dx = CDC
2328 set = Const::CDC;
2329 double cdc_sig = exp(pid->getLogL(c_belleHyp_to_chargedStable[sigHyp], set));
2330 double cdc_bkg = exp(pid->getLogL(c_belleHyp_to_chargedStable[bkgHyp], set));
2331 double cdc = 0.5;
2332 double cdc_all = cdc_sig + cdc_bkg;
2333 if (cdc_all != 0) {
2334 cdc = cdc_sig / cdc_all;
2335 if (cdc < 0.001) cdc = 0.001;
2336 if (cdc > 0.999) cdc = 0.999;
2337 }
2338
2339 // Combined
2340 double pid_sig = acc * tof * cdc;
2341 double pid_bkg = (1. - acc) * (1. - tof) * (1. - cdc);
2342
2343 return pid_sig / (pid_sig + pid_bkg);
2344}
Belle2::B2BIIConvertMdstModule::c_belleHyp_to_chargedStable
static const Const::ChargedStable c_belleHyp_to_chargedStable[c_nHyp]
maps Belle hypotheses to Const::ChargedStable (from http://belle.kek.jp/secured/wiki/doku....
Definition: B2BIIConvertMdstModule.h:323
Belle2::Const::RestrictedDetectorSet
A class for sets of detector IDs whose content is limited to restricted set of valid detector IDs.
Definition: Const.h:296

◆ beginRun()

void beginRun ( void  )
overridevirtual

Module functions to be called from event process.

Called when the current run begins.

Reimplemented from Module.

Definition at line 269 of file B2BIIConvertMdstModule.cc.

270{
271 B2DEBUG(99, "B2BIIConvertMdst: beginRun called.");
272
273 if (m_convertBeamParameters) {
274 //BeamEnergy class updated by fixmdst module in beginRun()
275 Belle::BeamEnergy::begin_run();
276 convertBeamEnergy();
277 Belle::BeamEnergy::dump();
278
279 // load IP data from DB server
280 Belle::IpProfile::begin_run();
281 convertIPProfile(true);
282 Belle::IpProfile::dump();
283 bool usableIP = Belle::IpProfile::usable();
284 B2DEBUG(99, "B2BIIConvertMdst: IpProfile is usable = " << usableIP);
285 }
286
287 //init eID
288#ifdef HAVE_EID
289 Belle::eid::init_data();
290 Belle::eid::show_use("ALL");
291#endif
292}
Belle2::B2BIIConvertMdstModule::convertIPProfile
void convertIPProfile(bool beginRun=false)
Stores the IPProfiles in BeamSpot (currently in DataStore)
Definition: B2BIIConvertMdstModule.cc:424
Belle2::B2BIIConvertMdstModule::convertBeamEnergy
void convertBeamEnergy()
Stores beam parameters (energy, angles) in CollisionInvariantMass and CollisionBoostVector (currently...
Definition: B2BIIConvertMdstModule.cc:386

◆ belleVeeDaughterHelix()

void belleVeeDaughterHelix ( const Belle::Mdst_vee2 &  vee,
const int  charge,
std::vector< float > &  helixParam,
std::vector< float > &  helixError 
)
private

obtains the helix parameters of the vee daughters

Definition at line 2211 of file B2BIIConvertMdstModule.cc.

2213{
2214 const HepPoint3D pivot(vee.vx(), vee.vy(), vee.vz());
2215 CLHEP::HepVector a(5);
2216 CLHEP::HepSymMatrix Ea(5, 0);
2217 if (charge > 0) {
2218 a[0] = vee.daut().helix_p(0); a[1] = vee.daut().helix_p(1);
2219 a[2] = vee.daut().helix_p(2); a[3] = vee.daut().helix_p(3);
2220 a[4] = vee.daut().helix_p(4);
2221 Ea[0][0] = vee.daut().error_p(0);
2222 Ea[1][0] = vee.daut().error_p(1);
2223 Ea[1][1] = vee.daut().error_p(2);
2224 Ea[2][0] = vee.daut().error_p(3);
2225 Ea[2][1] = vee.daut().error_p(4);
2226 Ea[2][2] = vee.daut().error_p(5);
2227 Ea[3][0] = vee.daut().error_p(6);
2228 Ea[3][1] = vee.daut().error_p(7);
2229 Ea[3][2] = vee.daut().error_p(8);
2230 Ea[3][3] = vee.daut().error_p(9);
2231 Ea[4][0] = vee.daut().error_p(10);
2232 Ea[4][1] = vee.daut().error_p(11);
2233 Ea[4][2] = vee.daut().error_p(12);
2234 Ea[4][3] = vee.daut().error_p(13);
2235 Ea[4][4] = vee.daut().error_p(14);
2236 } else {
2237 a[0] = vee.daut().helix_m(0); a[1] = vee.daut().helix_m(1);
2238 a[2] = vee.daut().helix_m(2); a[3] = vee.daut().helix_m(3);
2239 a[4] = vee.daut().helix_m(4);
2240 Ea[0][0] = vee.daut().error_m(0);
2241 Ea[1][0] = vee.daut().error_m(1);
2242 Ea[1][1] = vee.daut().error_m(2);
2243 Ea[2][0] = vee.daut().error_m(3);
2244 Ea[2][1] = vee.daut().error_m(4);
2245 Ea[2][2] = vee.daut().error_m(5);
2246 Ea[3][0] = vee.daut().error_m(6);
2247 Ea[3][1] = vee.daut().error_m(7);
2248 Ea[3][2] = vee.daut().error_m(8);
2249 Ea[3][3] = vee.daut().error_m(9);
2250 Ea[4][0] = vee.daut().error_m(10);
2251 Ea[4][1] = vee.daut().error_m(11);
2252 Ea[4][2] = vee.daut().error_m(12);
2253 Ea[4][3] = vee.daut().error_m(13);
2254 Ea[4][4] = vee.daut().error_m(14);
2255 }
2256
2257 Belle::Helix helix(pivot, a, Ea);
2258
2259 // go to the new pivot
2260 helix.pivot(HepPoint3D(0., 0., 0.));
2261
2262 CLHEP::HepSymMatrix error5x5(5, 0);
2263 convertHelix(helix, helixParam, error5x5);
2264
2265 unsigned int size = 5;
2266 unsigned int counter = 0;
2267 for (unsigned int i = 0; i < size; i++)
2268 for (unsigned int j = i; j < size; j++)
2269 helixError[counter++] = error5x5[i][j];
2270}
Belle2::B2BIIConvertMdstModule::convertHelix
void convertHelix(Belle::Helix &helix, std::vector< float > &helixParams, CLHEP::HepSymMatrix &error5x5)
Converts Belle's Helix parameters and it's covariance matrix to Belle II's version.
Definition: B2BIIConvertMdstModule.cc:1654

◆ belleVeeDaughterToCartesian()

void belleVeeDaughterToCartesian ( const Belle::Mdst_vee2 &  vee,
const int  charge,
const Const::ParticleType &  pType,
CLHEP::HepLorentzVector &  momentum,
HepPoint3D &  position,
CLHEP::HepSymMatrix &  error 
)
private

Fills 4-momentum, position and 7x7 error matrix from Belle Vee daughter.

Definition at line 2172 of file B2BIIConvertMdstModule.cc.

2175{
2176 const HepPoint3D pivot(vee.vx(), vee.vy(), vee.vz());
2177 CLHEP::HepVector a(5);
2178 CLHEP::HepSymMatrix Ea(5, 0);
2179 if (charge > 0) {
2180 a[0] = vee.daut().helix_p(0); a[1] = vee.daut().helix_p(1);
2181 a[2] = vee.daut().helix_p(2); a[3] = vee.daut().helix_p(3);
2182 a[4] = vee.daut().helix_p(4);
2183 Ea[0][0] = vee.daut().error_p(0); Ea[1][0] = vee.daut().error_p(1);
2184 Ea[1][1] = vee.daut().error_p(2); Ea[2][0] = vee.daut().error_p(3);
2185 Ea[2][1] = vee.daut().error_p(4); Ea[2][2] = vee.daut().error_p(5);
2186 Ea[3][0] = vee.daut().error_p(6); Ea[3][1] = vee.daut().error_p(7);
2187 Ea[3][2] = vee.daut().error_p(8); Ea[3][3] = vee.daut().error_p(9);
2188 Ea[4][0] = vee.daut().error_p(10); Ea[4][1] = vee.daut().error_p(11);
2189 Ea[4][2] = vee.daut().error_p(12); Ea[4][3] = vee.daut().error_p(13);
2190 Ea[4][4] = vee.daut().error_p(14);
2191 } else {
2192 a[0] = vee.daut().helix_m(0); a[1] = vee.daut().helix_m(1);
2193 a[2] = vee.daut().helix_m(2); a[3] = vee.daut().helix_m(3);
2194 a[4] = vee.daut().helix_m(4);
2195 Ea[0][0] = vee.daut().error_m(0); Ea[1][0] = vee.daut().error_m(1);
2196 Ea[1][1] = vee.daut().error_m(2); Ea[2][0] = vee.daut().error_m(3);
2197 Ea[2][1] = vee.daut().error_m(4); Ea[2][2] = vee.daut().error_m(5);
2198 Ea[3][0] = vee.daut().error_m(6); Ea[3][1] = vee.daut().error_m(7);
2199 Ea[3][2] = vee.daut().error_m(8); Ea[3][3] = vee.daut().error_m(9);
2200 Ea[4][0] = vee.daut().error_m(10); Ea[4][1] = vee.daut().error_m(11);
2201 Ea[4][2] = vee.daut().error_m(12); Ea[4][3] = vee.daut().error_m(13);
2202 Ea[4][4] = vee.daut().error_m(14);
2203 }
2204
2205 Belle::Helix helix(pivot, a, Ea);
2206
2207 // this is Vee daughter momentum/position/error at pivot = V0 Decay Vertex
2208 momentum = helix.momentum(0., pType.getMass(), position, error);
2209}
Belle2::Const::ParticleType::getMass
double getMass() const
Particle mass.
Definition: UnitConst.cc:353

◆ cdc_pid()

double cdc_pid ( const Belle::Mdst_charged &  chg,
int  idp 
)
private

Returns CDC likelihood for given hypothesis idp.

Copied from atc_pid::cdd_pid().

◆ clone()

std::shared_ptr< PathElement > clone ( ) const
overridevirtualinherited

Create an independent copy of this module.

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

Implements PathElement.

Definition at line 179 of file Module.cc.

180{
181 ModulePtr newModule = ModuleManager::Instance().registerModule(getType());
182 newModule->m_moduleParamList.setParameters(getParamList());
183 newModule->setName(getName());
184 newModule->m_package = m_package;
185 newModule->m_propertyFlags = m_propertyFlags;
186 newModule->m_logConfig = m_logConfig;
187 newModule->m_conditions = m_conditions;
188
189 return newModule;
190}
Belle2::ModuleManager::registerModule
std::shared_ptr< Module > registerModule(const std::string &moduleName, std::string sharedLibPath="") noexcept(false)
Creates an instance of a module and registers it to the ModuleManager.
Definition: ModuleManager.cc:80
Belle2::ModuleManager::Instance
static ModuleManager & Instance()
Exception is thrown if the requested module could not be created by the ModuleManager.
Definition: ModuleManager.cc:24
Belle2::Module::getParamList
const ModuleParamList & getParamList() const
Return module param list.
Definition: Module.h:362
Belle2::Module::getName
const std::string & getName() const
Returns the name of the module.
Definition: Module.h:186
Belle2::Module::getType
const std::string & getType() const
Returns the type of the module (i.e.
Definition: Module.cc:41
Belle2::Module::m_propertyFlags
unsigned int m_propertyFlags
The properties of the module as bitwise or (with |) of EModulePropFlags.
Definition: Module.h:511
Belle2::Module::m_logConfig
LogConfig m_logConfig
The log system configuration of the module.
Definition: Module.h:513
Belle2::Module::m_conditions
std::vector< ModuleCondition > m_conditions
Module condition, only non-null if set.
Definition: Module.h:520
Belle2::Module::m_package
std::string m_package
Package this module is found in (may be empty).
Definition: Module.h:509
Belle2::ModulePtr
std::shared_ptr< Module > ModulePtr
Defines a pointer to a module object as a boost shared pointer.
Definition: Module.h:43

◆ computeTrkMinDistanceBelle()

double computeTrkMinDistanceBelle ( ECLCluster *  eclCluster)
private

calculate the minimal distance between a cluster and a set of tracks on the ECL surface.

Definition at line 1941 of file B2BIIConvertMdstModule.cc.

1942{
1943 const double m_extRadius(125.0);
1944 const double m_extZFWD(196.0);
1945 const double m_extZBWD(-102.2);
1946 double minDist(10000);
1947
1948 // get cluster info
1949 const int reg = eclCluster->getDetectorRegion();
1950 double eclClusterR_surface = m_extRadius / sin(eclCluster->getTheta());
1951 if (reg == 1) {eclClusterR_surface = m_extZFWD / cos(eclCluster->getTheta());}
1952 else if (reg == 3) {eclClusterR_surface = m_extZBWD / cos(eclCluster->getTheta());}
1953
1954 ROOT::Math::XYZVector eclCluster_surface_position(0, 0, 0);
1955 VectorUtil::setMagThetaPhi(eclCluster_surface_position, eclClusterR_surface, eclCluster->getTheta(), eclCluster->getPhi());
1956
1957 for (const auto& track : m_tracks) {
1958 const TrackFitResult* trackFit = track.getTrackFitResultWithClosestMass(Const::ChargedStable(Const::pion));
1959
1960 if (trackFit == NULL) {continue;}
1961 // get the track parameters
1962 const double z0 = trackFit->getZ0();
1963 const double tanlambda = trackFit->getTanLambda();
1964
1965 // use the helix class
1966 Helix h = trackFit->getHelix();
1967
1968 // extrapolate to radius
1969 const double lHelixRadius = h.getArcLength2DAtCylindricalR(m_extRadius) > 0 ? h.getArcLength2DAtCylindricalR(m_extRadius) : 999999.;
1970
1971 // extrapolate to FWD z
1972 const double lFWD = (m_extZFWD - z0) / tanlambda > 0 ? (m_extZFWD - z0) / tanlambda : 999999.;
1973
1974 // extrapolate to backward z
1975 const double lBWD = (m_extZBWD - z0) / tanlambda > 0 ? (m_extZBWD - z0) / tanlambda : 999999.;
1976
1977 // pick smallest arclength
1978 const double l = std::min(std::min(lHelixRadius, lFWD), lBWD);
1979
1980 B2DEBUG(50, lHelixRadius << " " << lFWD << " " << lBWD << " -> " << l);
1981
1982 ROOT::Math::XYZVector ext_helix = h.getPositionAtArcLength2D(l);
1983 double helixExtR_surface = m_extRadius / sin(ext_helix.Theta());
1984 if (l == lFWD) { helixExtR_surface = m_extZFWD / cos(ext_helix.Theta());}
1985 else if (l == lBWD) { helixExtR_surface = m_extZBWD / cos(ext_helix.Theta());}
1986
1987 ROOT::Math::XYZVector helixExt_surface_position(0, 0, 0);
1988 VectorUtil::setMagThetaPhi(helixExt_surface_position, helixExtR_surface, ext_helix.Theta(), ext_helix.Phi());
1989
1990 double distance = (eclCluster_surface_position - helixExt_surface_position).R();
1991 if (distance < minDist) {minDist = distance;}
1992 }
1993 if (minDist > 9999) minDist = -1;
1994 return minDist;
1995}
R
double R
typedef autogenerated by FFTW
Definition: DiscreteCosineTransform_31points.cc:35
Belle2::B2BIIConvertMdstModule::m_tracks
StoreArray< Track > m_tracks
Tracks.
Definition: B2BIIConvertMdstModule.h:390
Belle2::CDC::Helix
Helix parameter class.
Definition: Helix.h:48
Belle2::Const::ChargedStable
Provides a type-safe way to pass members of the chargedStableSet set.
Definition: Const.h:589
Belle2::Const::pion
static const ChargedStable pion
charged pion particle
Definition: Const.h:661
Belle2::ECLCluster::getPhi
double getPhi() const
Return Corrected Phi of Shower (radian).
Definition: ECLCluster.h:304
Belle2::ECLCluster::getDetectorRegion
int getDetectorRegion() const
Return detector region: 0: below acceptance, 1: FWD, 2: BRL, 3: BWD, 11: FWDGAP, 13: BWDGAP.
Definition: ECLCluster.cc:70
Belle2::ECLCluster::getTheta
double getTheta() const
Return Corrected Theta of Shower (radian).
Definition: ECLCluster.h:307
Belle2::TrackFitResult
Values of the result of a track fit with a given particle hypothesis.
Definition: TrackFitResult.h:49
Belle2::TrackFitResult::getHelix
Helix getHelix() const
Conversion to framework Helix (without covariance).
Definition: TrackFitResult.h:231
Belle2::TrackFitResult::getTanLambda
double getTanLambda() const
Getter for tanLambda.
Definition: TrackFitResult.h:207
Belle2::TrackFitResult::getZ0
double getZ0() const
Getter for z0.
Definition: TrackFitResult.h:201

◆ convertBeamEnergy()

void convertBeamEnergy ( )
private

Stores beam parameters (energy, angles) in CollisionInvariantMass and CollisionBoostVector (currently in the DataStore).

Definition at line 386 of file B2BIIConvertMdstModule.cc.

387{
388 const double Eher = Belle::BeamEnergy::E_HER();
389 const double Eler = Belle::BeamEnergy::E_LER();
390 const double crossingAngle = Belle::BeamEnergy::Cross_angle();
391 const double angleLer = M_PI; //parallel to negative z axis (different from Belle II!)
392 const double angleHer = crossingAngle; //in positive z and x direction, verified to be consistent with Upsilon(4S) momentum
393 const double mass_e = Const::electronMass; //mass of electron: 0.0 in basf, 0.000510998902 in basf2;
394 TMatrixDSym covariance(3); //0 entries = no error
395 HepLorentzVector p_beam = Belle::BeamEnergy::p_beam(); // Testing only
396
397 // Get four momentum of LER and HER
398 ROOT::Math::PxPyPzEVector P_her(0.0, 0.0, TMath::Sqrt(Eher * Eher - mass_e * mass_e), Eher);
399 ROOT::Math::RotationY rotateAroundYAxis(angleHer);
400 P_her = rotateAroundYAxis(P_her);
401 ROOT::Math::PxPyPzEVector P_ler(0.0, 0.0, TMath::Sqrt(Eler * Eler - mass_e * mass_e), Eler);
402 rotateAroundYAxis.SetAngle(angleLer);
403 P_ler = rotateAroundYAxis(P_ler);
404
405 // Get four momentum of beam
406 ROOT::Math::PxPyPzEVector P_beam = P_her + P_ler;
407
408 ROOT::Math::PxPyPzEVector momentumHER = P_her;
409 ROOT::Math::VectorUtil::boost(momentumHER, -P_beam.BoostToCM());
410 double angleXZ = std::atan(momentumHER.X() / momentumHER.Z());
411 double angleYZ = std::atan(momentumHER.Y() / momentumHER.Z());
412
413 m_collisionBoostVector.setBoost(B2Vector3D(-P_beam.BoostToCM()), covariance);
414 m_collisionInvM.setMass(P_beam.M(), 0.0, 0.0);
415 m_collisionAxisCMS.setAngles(angleXZ, angleYZ, TMatrixTSym<double>(2));
416 m_collisionAxisCMS.setSpread(TMatrixTSym<double>(2), 0, 0, 0);
417
418 B2DEBUG(99, "Beam Energy: E_HER = " << Eher << "; E_LER = " << Eler << "; angle = " << crossingAngle);
419 B2DEBUG(99, "Beam Momentum (pre-convert) : P_X = " << p_beam.px() << "; P_Y = " << p_beam.py() << "; P_Z = " << p_beam.pz());
420 B2DEBUG(99, "Beam Momentum (post-convert) : P_X = " << P_beam.Px() << "; P_Y = " << P_beam.Py() << "; P_Z = " << P_beam.Pz());
421 B2DEBUG(99, "CollisionAxisCMS: angleXZ = " << m_collisionAxisCMS.getAngleXZ() << "; angleYZ = " << m_collisionAxisCMS.getAngleYZ());
422}
Belle2::B2BIIConvertMdstModule::m_collisionAxisCMS
CollisionAxisCMS m_collisionAxisCMS
CollisionAxisCMS of the beam.
Definition: B2BIIConvertMdstModule.h:424
Belle2::B2BIIConvertMdstModule::m_collisionInvM
CollisionInvariantMass m_collisionInvM
CollisionInvariantMass for the invariant mass of the beam.
Definition: B2BIIConvertMdstModule.h:420
Belle2::B2BIIConvertMdstModule::m_collisionBoostVector
CollisionBoostVector m_collisionBoostVector
CollisionBoostVector for bosst vector of the beam.
Definition: B2BIIConvertMdstModule.h:416
Belle2::CollisionAxisCMS::getAngleYZ
double getAngleYZ() const
Get the mean YZ angle of the HER beam in the CM frame obtained by pure boost.
Definition: CollisionAxisCMS.h:60
Belle2::CollisionAxisCMS::setSpread
void setSpread(const TMatrixDSym &spreadCovariance, double spreadXZunc, double spreadYZunc, double spreadPhiUnc)
Set spread covariance and uncertainties of the eigenvalues of this matrix.
Definition: CollisionAxisCMS.h:47
Belle2::CollisionAxisCMS::getAngleXZ
double getAngleXZ() const
Get the mean XZ angle of the HER beam in the CM frame obtained by pure boost.
Definition: CollisionAxisCMS.h:57
Belle2::CollisionAxisCMS::setAngles
void setAngles(double angleXZ, double angleYZ, const TMatrixDSym &centerCovariance)
Set the central values and uncertainty of them.
Definition: CollisionAxisCMS.h:39
Belle2::CollisionBoostVector::setBoost
void setBoost(const TVector3 &boost, const TMatrixDSym &covariance)
Set the boost vector and its error matrix.
Definition: CollisionBoostVector.h:36
Belle2::CollisionInvariantMass::setMass
void setMass(double mass, double error, double spread)
Set the CMS energy and its uncertainty.
Definition: CollisionInvariantMass.h:35
Belle2::Const::electronMass
static const double electronMass
electron mass
Definition: Const.h:685
Belle2::B2Vector3D
B2Vector3< double > B2Vector3D
typedef for common usage with double
Definition: B2Vector3.h:516

◆ convertEvtclsTable()

void convertEvtclsTable ( )
private

Reads and converts all entries of evtcls Panther table.

Definition at line 1266 of file B2BIIConvertMdstModule.cc.

1267{
1268 // Create StoreObj if it is not valid
1269 if (not m_evtInfo.isValid()) {
1270 m_evtInfo.create();
1271 }
1272 // Pull Evtcls_flag(2) from manager
1273 Belle::Evtcls_flag_Manager& EvtFlagMgr = Belle::Evtcls_flag_Manager::get_manager();
1274 Belle::Evtcls_flag2_Manager& EvtFlag2Mgr = Belle::Evtcls_flag2_Manager::get_manager();
1275
1276 // Pull Evtcls_hadronic_flag from manager
1277 Belle::Evtcls_hadronic_flag_Manager& EvtHadFlagMgr = Belle::Evtcls_hadronic_flag_Manager::get_manager();
1278
1279 std::string name = "evtcls_flag";
1280 std::string name_had = "evtcls_hadronic_flag";
1281 // Only one entry in each event
1282 std::vector<Belle::Evtcls_flag>::iterator eflagIterator = EvtFlagMgr.begin();
1283 std::vector<Belle::Evtcls_flag2>::iterator eflag2Iterator = EvtFlag2Mgr.begin();
1284 std::vector<Belle::Evtcls_hadronic_flag>::iterator ehadflagIterator = EvtHadFlagMgr.begin();
1285
1286 // Converting evtcls_flag(2)
1287 std::vector<int> flag(20);
1288 for (int index = 0; index < 20; ++index) {
1289 // flag(14, 16): not filled
1290 if (index == 14 || index == 16) continue;
1291 std::string iVar = name + std::to_string(index);
1292 // 0-9 corresponding to evtcls_flag.flag(0-9)
1293 if (index < 10) {
1294 m_evtInfo->addExtraInfo(iVar, (*eflagIterator).flag(index));
1295 } else {
1296 // 10-19 corresponding to evtcls_flag2.flag(0-9)
1297 m_evtInfo->addExtraInfo(iVar, (*eflag2Iterator).flag(index - 10));
1298 }
1299 B2DEBUG(99, "evtcls_flag(" << index << ") = " << m_evtInfo->getExtraInfo(iVar));
1300 }
1301
1302 // Converting evtcls_hadronic_flag
1303 for (int index = 0; index < 6; ++index) {
1304 std::string iVar = name_had + std::to_string(index);
1305 m_evtInfo->addExtraInfo(iVar, (*ehadflagIterator).hadronic_flag(index));
1306 B2DEBUG(99, "evtcls_hadronic_flag(" << index << ") = " << m_evtInfo->getExtraInfo(iVar));
1307 }
1308
1309}
Belle2::B2BIIConvertMdstModule::m_evtInfo
StoreObjPtr< EventExtraInfo > m_evtInfo
Event Extra Info.
Definition: B2BIIConvertMdstModule.h:408

◆ convertGenHepevtObject()

void convertGenHepevtObject ( const Belle::Gen_hepevt &  genHepevt,
MCParticleGraph::GraphParticle *  mcParticle 
)
private

Converts Gen_hepevt record to MCParticleGraph object.

Definition at line 1856 of file B2BIIConvertMdstModule.cc.

1857{
1858 //B2DEBUG(80, "Gen_ehepevt: idhep " << genHepevt.idhep() << " (" << genHepevt.isthep() << ") with ID = " << genHepevt.get_ID());
1859
1860 // updating the GraphParticle information from the Gen_hepevt information
1861 const int idHep = recoverMoreThan24bitIDHEP(genHepevt.idhep());
1862 if (idHep == 0) {
1863 B2WARNING("Trying to convert Gen_hepevt with idhep = " << idHep <<
1864 ". This should never happen.");
1865 mcParticle->setPDG(Const::photon.getPDGCode());
1866 } else {
1867 mcParticle->setPDG(idHep);
1868 }
1869
1870 if (genHepevt.isthep() > 0) {
1871 mcParticle->setStatus(MCParticle::c_PrimaryParticle);
1872 }
1873
1874 mcParticle->setMass(genHepevt.M());
1875
1876 ROOT::Math::PxPyPzEVector p4(genHepevt.PX(), genHepevt.PY(), genHepevt.PZ(), genHepevt.E());
1877 mcParticle->set4Vector(p4);
1878
1879 mcParticle->setProductionVertex(genHepevt.VX()*Unit::mm, genHepevt.VY()*Unit::mm, genHepevt.VZ()*Unit::mm);
1880 mcParticle->setProductionTime(genHepevt.T()*Unit::mm / Const::speedOfLight);
1881
1882 // decay time of this particle is production time of the daughter particle
1883 if (genHepevt.daFirst() > 0) {
1884 Belle::Gen_hepevt_Manager& genMgr = Belle::Gen_hepevt_Manager::get_manager();
1885 Belle::Gen_hepevt daughterParticle = genMgr(Belle::Panther_ID(genHepevt.daFirst()));
1886 mcParticle->setDecayTime(daughterParticle.T()*Unit::mm / Const::speedOfLight);
1887 mcParticle->setDecayVertex(daughterParticle.VX()*Unit::mm, daughterParticle.VY()*Unit::mm, daughterParticle.VZ()*Unit::mm);
1888 } else {
1889 //otherwise, assume it's stable
1890 mcParticle->setDecayTime(std::numeric_limits<float>::infinity());
1891 }
1892
1893 mcParticle->setValidVertex(true);
1894}
Belle2::B2BIIConvertMdstModule::recoverMoreThan24bitIDHEP
int recoverMoreThan24bitIDHEP(int id)
Helper function to recover falsely set idhep info in GenHepEvt list.
Definition: B2BIIConvertMdstModule.cc:2077
Belle2::Const::speedOfLight
static const double speedOfLight
[cm/ns]
Definition: Const.h:695
Belle2::Const::photon
static const ParticleType photon
photon particle
Definition: Const.h:673
Belle2::MCParticle::c_PrimaryParticle
@ c_PrimaryParticle
bit 0: Particle is primary particle.
Definition: MCParticle.h:47
Belle2::MCParticle::setDecayTime
void setDecayTime(float time)
Set decay time.
Definition: MCParticle.h:379
Belle2::MCParticle::setMass
void setMass(float mass)
Set particle mass.
Definition: MCParticle.h:355
Belle2::MCParticle::setDecayVertex
void setDecayVertex(const ROOT::Math::XYZVector &vertex)
Set decay vertex.
Definition: MCParticle.h:436
Belle2::MCParticle::setValidVertex
void setValidVertex(bool valid)
Set indication whether vertex and time information is valid or just default.
Definition: MCParticle.h:367
Belle2::MCParticle::setProductionVertex
void setProductionVertex(const ROOT::Math::XYZVector &vertex)
Set production vertex position.
Definition: MCParticle.h:385
Belle2::MCParticle::setPDG
void setPDG(int pdg)
Set PDG code of the particle.
Definition: MCParticle.h:324
Belle2::MCParticle::set4Vector
void set4Vector(const ROOT::Math::PxPyPzEVector &p4)
Sets the 4Vector of particle.
Definition: MCParticle.h:427
Belle2::MCParticle::setStatus
void setStatus(unsigned short int status)
Set Status code for the particle.
Definition: MCParticle.h:335
Belle2::MCParticle::setProductionTime
void setProductionTime(float time)
Set production time.
Definition: MCParticle.h:373
Belle2::Unit::mm
static const double mm
[millimeters]
Definition: Unit.h:70

◆ convertGenHepEvtTable()

void convertGenHepEvtTable ( )
private

Reads and converts all entries of Gen_hepevt Panther table to MCParticle dataobjects and adds them to StoreArray<MCParticle>.

Definition at line 884 of file B2BIIConvertMdstModule.cc.

885{
886 if (m_realData)
887 return;
888
889 // clear the Gen_hepevt_ID <-> MCParticleGraphPosition map
890 genHepevtToMCParticle.clear();
891
892 // check if the Gen_hepevt table has any entries
893 Belle::Gen_hepevt_Manager& genMgr = Belle::Gen_hepevt_Manager::get_manager();
894 if (genMgr.count() == 0)
895 return;
896
897 typedef std::pair<MCParticleGraph::GraphParticle*, Belle::Gen_hepevt> halfFamily;
898 halfFamily currFamily;
899 halfFamily family;
900 std::queue < halfFamily > heritancesQueue;
901
902 // Add motherless particles. The root particle is often the first one,
903 // but this is not correct in general case; thus, all particles are added,
904 // including the beam-background ones.
905 m_particleGraph.clear();
906 for (Belle::Gen_hepevt_Manager::iterator genIterator = genMgr.begin();
907 genIterator != genMgr.end(); ++genIterator) {
908 Belle::Gen_hepevt hep = *genIterator;
909 // Select particles without mother.
910 if (!(hep.moFirst() == 0 && hep.moLast() == 0))
911 continue;
912
913 int position = m_particleGraph.size();
914 m_particleGraph.addParticle();
915 genHepevtToMCParticle[hep.get_ID()] = position;
916 MCParticleGraph::GraphParticle* graphParticle = &m_particleGraph[position];
917 convertGenHepevtObject(hep, graphParticle);
918 for (int iDaughter = hep.daFirst(); iDaughter <= hep.daLast();
919 ++iDaughter) {
920 if (iDaughter == 0) {
921 B2DEBUG(95, "Trying to access generated daughter with Panther ID == 0");
922 continue;
923 }
924 currFamily.first = graphParticle;
925 currFamily.second = genMgr(Belle::Panther_ID(iDaughter));
926 heritancesQueue.push(currFamily);
927 }
928 }
929
930 //now we can go through the queue:
931 while (!heritancesQueue.empty()) {
932 currFamily = heritancesQueue.front(); //get the first entry from the queue
933 heritancesQueue.pop(); //remove the entry.
934
935 MCParticleGraph::GraphParticle* currMother = currFamily.first;
936 Belle::Gen_hepevt& currDaughter = currFamily.second;
937
938 // Skip particles with idhep = 0.
939 if (currDaughter.idhep() == 0)
940 continue;
941
942 //putting the daughter in the graph:
943 int position = m_particleGraph.size();
944 m_particleGraph.addParticle();
945 genHepevtToMCParticle[currDaughter.get_ID()] = position;
946
947 MCParticleGraph::GraphParticle* graphDaughter = &m_particleGraph[position];
948 convertGenHepevtObject(currDaughter, graphDaughter);
949
950 //add relation between mother and daughter to graph:
951 currMother->decaysInto((*graphDaughter));
952
953 int nGrandChildren = currDaughter.daLast() - currDaughter.daFirst() + 1;
954
955 if (nGrandChildren > 0 && currDaughter.daFirst() != 0) {
956 for (int igrandchild = currDaughter.daFirst(); igrandchild <= currDaughter.daLast(); ++igrandchild) {
957 if (igrandchild == 0) {
958 B2DEBUG(95, "Trying to access generated daughter with Panther ID == 0");
959 continue;
960 }
961
962 family.first = graphDaughter;
963 family.second = genMgr(Belle::Panther_ID(igrandchild));
964 heritancesQueue.push(family);
965 }
966 }
967 }
968
969 m_particleGraph.generateList();
970}
Belle2::B2BIIConvertMdstModule::genHepevtToMCParticle
std::map< int, int > genHepevtToMCParticle
map of Gen_hepevt Panther IDs and corresponding MCParticle StoreArray indices
Definition: B2BIIConvertMdstModule.h:370
Belle2::B2BIIConvertMdstModule::m_particleGraph
Belle2::MCParticleGraph m_particleGraph
MCParticle Graph to build Belle2 MC Particles.
Definition: B2BIIConvertMdstModule.h:367
Belle2::B2BIIConvertMdstModule::convertGenHepevtObject
void convertGenHepevtObject(const Belle::Gen_hepevt &genHepevt, MCParticleGraph::GraphParticle *mcParticle)
Converts Gen_hepevt record to MCParticleGraph object.
Definition: B2BIIConvertMdstModule.cc:1856
Belle2::MCParticleGraph::GraphParticle
Class to represent Particle data in graph.
Definition: MCParticleGraph.h:75
Belle2::MCParticleGraph::GraphParticle::decaysInto
void decaysInto(GraphParticle &daughter)
Tells the graph that this particle decays into daughter.
Definition: MCParticleGraph.h:99
Belle2::MCParticleGraph::size
size_t size() const
Return the number of particles in the graph.
Definition: MCParticleGraph.h:242
Belle2::MCParticleGraph::generateList
void generateList(const std::string &name="", int options=c_setNothing)
Generates the MCParticle list and stores it in the StoreArray with the given name.
Definition: MCParticleGraph.cc:208
Belle2::MCParticleGraph::clear
void clear()
Reset particles and decay information to make the class reusable.
Definition: MCParticleGraph.h:287
Belle2::MCParticleGraph::addParticle
GraphParticle & addParticle()
Add new particle to the graph.
Definition: MCParticleGraph.h:294

◆ convertHelix() [1/2]

void convertHelix ( Belle::Helix &  helix,
std::vector< float > &  helixParams,
CLHEP::HepSymMatrix &  error5x5 
)
private

Converts Belle's Helix parameters and it's covariance matrix to Belle II's version.

Definition at line 1654 of file B2BIIConvertMdstModule.cc.

1655{
1656 CLHEP::HepVector a(5);
1657 CLHEP::HepSymMatrix Ea(5, 0);
1658
1659 a = helix.a();
1660 Ea = helix.Ea();
1661
1662 // param 0: d_0 = d_rho
1663 helixParams[0] = a[0];
1664
1665 // param 1: phi = phi_0 + pi/2
1666 helixParams[1] = adjustAngleRange(a[1] + TMath::Pi() / 2.0);
1667
1668 // param 2: omega = Kappa * alpha = Kappa * B[Tesla] * speed_of_light[m/s] * 1e-11
1669 helixParams[2] = a[2] * KAPPA2OMEGA;
1670
1671 // param 3: d_z = z0
1672 helixParams[3] = a[3];
1673
1674 // param 4: tan(Lambda) = tanLambda
1675 helixParams[4] = a[4];
1676
1677 unsigned int size = 5;
1678 for (unsigned int i = 0; i < size; i++) {
1679 for (unsigned int j = 0; j < size; j++) {
1680 error5x5[i][j] = Ea[i][j];
1681 if (i == 2)
1682 error5x5[i][j] *= KAPPA2OMEGA;
1683 if (j == 2)
1684 error5x5[i][j] *= KAPPA2OMEGA;
1685
1686 if (std::isinf(error5x5[i][j])) {
1687 B2DEBUG(99, "Helix covariance matrix element found to be infinite. Setting value to DBL_MAX/2.0.");
1688 error5x5[i][j] = DBL_MAX / 2.0;
1689 }
1690 }
1691 }
1692}
Belle2::B2BIIConvertMdstModule::KAPPA2OMEGA
const double KAPPA2OMEGA
Conversion factor for Kappa -> Omega helix parameters.
Definition: B2BIIConvertMdstModule.h:441

◆ convertHelix() [2/2]

void convertHelix ( const Belle::Mdst_trk_fit &  trk_fit,
const HepPoint3D &  newPivot,
std::vector< float > &  helixParams,
std::vector< float > &  helixError 
)
private

Converts Belle's Helix parameters and it's covariance matrix to Belle II's version.

Definition at line 1603 of file B2BIIConvertMdstModule.cc.

1606{
1607 const HepPoint3D pivot(trk_fit.pivot_x(),
1608 trk_fit.pivot_y(),
1609 trk_fit.pivot_z());
1610
1611 CLHEP::HepVector a(5);
1612 a[0] = trk_fit.helix(0);
1613 a[1] = trk_fit.helix(1);
1614 a[2] = trk_fit.helix(2);
1615 a[3] = trk_fit.helix(3);
1616 a[4] = trk_fit.helix(4);
1617 CLHEP::HepSymMatrix Ea(5, 0);
1618 Ea[0][0] = trk_fit.error(0);
1619 Ea[1][0] = trk_fit.error(1);
1620 Ea[1][1] = trk_fit.error(2);
1621 Ea[2][0] = trk_fit.error(3);
1622 Ea[2][1] = trk_fit.error(4);
1623 Ea[2][2] = trk_fit.error(5);
1624 Ea[3][0] = trk_fit.error(6);
1625 Ea[3][1] = trk_fit.error(7);
1626 Ea[3][2] = trk_fit.error(8);
1627 Ea[3][3] = trk_fit.error(9);
1628 Ea[4][0] = trk_fit.error(10);
1629 Ea[4][1] = trk_fit.error(11);
1630 Ea[4][2] = trk_fit.error(12);
1631 Ea[4][3] = trk_fit.error(13);
1632 Ea[4][4] = trk_fit.error(14);
1633
1634 Belle::Helix helix(pivot, a, Ea);
1635
1636 if (newPivot.x() != 0. || newPivot.y() != 0. || newPivot.z() != 0.) {
1637 helix.pivot(newPivot);
1638 } else {
1639 if (pivot.x() != 0. || pivot.y() != 0. || pivot.z() != 0.) {
1640 helix.pivot(HepPoint3D(0., 0., 0.));
1641 }
1642 }
1643
1644 CLHEP::HepSymMatrix error5x5(5, 0);
1645 convertHelix(helix, helixParams, error5x5);
1646
1647 unsigned int size = 5;
1648 unsigned int counter = 0;
1649 for (unsigned int i = 0; i < size; i++)
1650 for (unsigned int j = i; j < size; j++)
1651 helixError[counter++] = error5x5[i][j];
1652}

◆ convertIPProfile()

void convertIPProfile ( bool  beginRun = false)
private

Stores the IPProfiles in BeamSpot (currently in DataStore)

Definition at line 424 of file B2BIIConvertMdstModule.cc.

425{
426 if (!Belle::IpProfile::usable()) {
427 // No IPProfile for this run ...
428 if (beginRun) {
429 // no IPProfile, set vertex to NaN without errors for the full run
430 m_beamSpot.setIP(
431 TVector3(std::numeric_limits<double>::quiet_NaN(),
432 std::numeric_limits<double>::quiet_NaN(),
433 std::numeric_limits<double>::quiet_NaN()
434 ), TMatrixTSym<double>()
435 );
436 }
437 return;
438 }
439 HepPoint3D ip;
440 CLHEP::HepSymMatrix ipErr;
441 if (beginRun) {
442 // use event independent average in begin run
443 ip = Belle::IpProfile::position();
444 ipErr = Belle::IpProfile::position_err();
445 } else {
446 // update evtbin
447 Belle::IpProfile::set_evtbin_number();
448 // check if it changed, if not there's nothing to do
449 if (Belle::IpProfile::EvtBinNo() == m_lastIPProfileBin) return;
450 // get event dependent position and error
451 ip = Belle::IpProfile::e_position();
452 ipErr = Belle::IpProfile::e_position_err();
453 }
454 // reset last ipprofile bin
455 m_lastIPProfileBin = Belle::IpProfile::EvtBinNo();
456
457 TMatrixDSym cov(ipErr.num_col());
458 for (int i = 0; i < ipErr.num_row(); ++i) {
459 for (int j = 0; j < ipErr.num_col(); ++j) {
460 cov(i, j) = ipErr(i + 1, j + 1);
461 }
462 }
463 m_beamSpot.setIP(TVector3(ip.x(), ip.y(), ip.z()), cov);
464}
Belle2::B2BIIConvertMdstModule::m_beamSpot
BeamSpot m_beamSpot
Interaction Point of the beam.
Definition: B2BIIConvertMdstModule.h:412
Belle2::B2BIIConvertMdstModule::beginRun
virtual void beginRun() override
Module functions to be called from event process.
Definition: B2BIIConvertMdstModule.cc:269
Belle2::B2BIIConvertMdstModule::m_lastIPProfileBin
int m_lastIPProfileBin
variable to tell us which IPProfile bin was active last time we looked
Definition: B2BIIConvertMdstModule.h:168
Belle2::BeamSpot::setIP
void setIP(const TVector3 &ipPosition, const TMatrixDSym &covariance)
Set the IP position and its error matrix.
Definition: BeamSpot.h:59

◆ convertMdstChargedObject()

void convertMdstChargedObject ( const Belle::Mdst_charged &  belleTrack,
Track *  track 
)
private

Converts Mdst_charged (Mdst_trk(_fit)) record to Track (TrackFitResult) object.

If running on MC, the Track -> MCParticle relation is set as well.

Definition at line 1694 of file B2BIIConvertMdstModule.cc.

1695{
1696 Belle::Mdst_trk& trk = belleTrack.trk();
1697
1698 for (int mhyp = 0 ; mhyp < c_nHyp; ++mhyp) {
1699 Belle::Mdst_trk_fit& trk_fit = trk.mhyp(mhyp);
1700
1701 double pValue = TMath::Prob(trk_fit.chisq(), trk_fit.ndf());
1702 if (std::isnan(pValue)) continue;
1703
1704 const Const::ChargedStable& pType = c_belleHyp_to_chargedStable[mhyp];
1705 double thisMass = pType.getMass();
1706
1707 // Converted helix parameters
1708 std::vector<float> helixParam(5);
1709 // Converted 5x5 error matrix
1710 CLHEP::HepSymMatrix error5x5(5, 0);
1711 // 4-momentum
1712 CLHEP::HepLorentzVector momentum;
1713 // 7x7 (momentum, position) error matrix
1714 CLHEP::HepSymMatrix error7x7(7, 0);
1715 // position
1716 HepPoint3D position;
1717
1718 getHelixParameters(trk_fit, thisMass, HepPoint3D(0., 0., 0.),
1719 helixParam, error5x5,
1720 momentum, position, error7x7, 0.0);
1721
1722 std::vector<float> helixError(15);
1723 unsigned int size = 5;
1724 unsigned int counter = 0;
1725 for (unsigned int i = 0; i < size; i++)
1726 for (unsigned int j = i; j < size; j++)
1727 helixError[counter++] = error5x5[i][j];
1728
1729 // Create an empty cdc hitpattern and set the number of total hits
1730 // use hits from 0: axial-wire, 1:stereo-wire, 2:cathode
1731 // the actual cdc hitpattern is not converted
1732
1733 int cdcNHits = 0;
1734 for (unsigned int i = 0; i < 3; i++)
1735 cdcNHits += trk_fit.nhits(i);
1736
1737 HitPatternCDC patternCdc;
1738 patternCdc.setNHits(cdcNHits);
1739
1740 // conversion of track position in CDC layers
1741 if (m_convertTrkExtra) {
1742 double tof = 0;
1743 double path_length = 0;
1744 double tof_sigma = 0;
1745 short tof_qual = 0;
1746 int acc_ph = 0;
1747 short acc_qual = 0;
1748 double dedx = trk.dEdx();
1749 short dedx_qual = trk.quality_dedx();
1750
1751 const Belle::Mdst_tof& tof_obj = belleTrack.tof();
1752 if (tof_obj) {
1753 tof = tof_obj.tof();
1754 path_length = tof_obj.path_length();
1755 tof_qual = tof_obj.quality();
1756 tof_sigma = tof_obj.sigma_tof();
1757 }
1758
1759 const Belle::Mdst_acc& acc_obj = belleTrack.acc();
1760 if (acc_obj) {
1761 acc_ph = acc_obj.photo_electron();
1762 acc_qual = acc_obj.quality();
1763 }
1764
1765
1766 auto cdcExtraInfo = m_belleTrkExtra.appendNew(trk_fit.first_x(), trk_fit.first_y(), trk_fit.first_z(),
1767 trk_fit.last_x(), trk_fit.last_y(), trk_fit.last_z(),
1768 tof, path_length, tof_qual, tof_sigma,
1769 acc_ph, acc_qual, dedx, dedx_qual);
1770 track->addRelationTo(cdcExtraInfo);
1771 }
1772 // conversion of the SVD hit pattern
1773 int svdHitPattern = trk_fit.hit_svd();
1774 // use hits from 3: SVD-rphi, 4: SVD-z
1775 // int svdNHits = trk_fit.nhits(3) + trk_fit.nhits(4);
1776
1777 std::bitset<32> svdBitSet(svdHitPattern);
1778
1779 HitPatternVXD patternVxd;
1780
1781 unsigned short svdLayers;
1782 // taken from: http://belle.kek.jp/group/indirectcp/cpfit/cpfit-festa/2004/talks/Apr.14/CPfesta-2005-Higuchi(3).pdf
1783 StoreObjPtr<EventMetaData> event;
1784 // mask for the rphi hits, first 6 (8) bits/ 2 bits per layer
1785 std::bitset<32> svdUMask(static_cast<std::string>("00000000000000000000000000000011"));
1786 // mask for the z hits, second 6 (8) bits/ 2 bits per layer
1787 std::bitset<32> svdVMask;
1788
1789 // find out if the SVD has 3 (4) layers; if exp <= (>) exp 27
1790 if (event->getExperiment() <= 27) {
1791 svdVMask = svdUMask << 6;
1792 svdLayers = 3;
1793 } else {
1794 svdVMask = svdUMask << 8;
1795 svdLayers = 4;
1796 }
1797
1798 // loop over all svd layers (layer index is shifted + 3 for basf2)
1799 for (unsigned short layerId = 0; layerId < svdLayers; layerId++) {
1800 unsigned short uHits = (svdBitSet & svdUMask).count();
1801 unsigned short vHits = (svdBitSet & svdVMask).count();
1802 patternVxd.setSVDLayer(layerId + 3, uHits, vHits);
1803 // shift masks to the left
1804 svdUMask <<= 2;
1805 svdVMask <<= 2;
1806 }
1807
1808 TrackFitResult helixFromHelix(helixParam, helixError, pType, pValue, -1, patternVxd.getInteger(), 0);
1809
1810 if (m_use6x6CovarianceMatrix4Tracks) {
1811 TMatrixDSym cartesianCovariance(6);
1812 for (unsigned i = 0; i < 7; i++) {
1813 if (i == 3)
1814 continue;
1815 for (unsigned j = 0; j < 7; j++) {
1816 if (j == 3)
1817 continue;
1818
1819 cartesianCovariance(ERRMCONV[i], ERRMCONV[j]) = error7x7[i][j];
1820 }
1821 }
1822 UncertainHelix helixFromCartesian(helixFromHelix.getPosition(), helixFromHelix.getMomentum(), helixFromHelix.getChargeSign(),
1823 BFIELD, cartesianCovariance, pValue);
1824
1825 TMatrixDSym helixCovariance = helixFromCartesian.getCovariance();
1826
1827 counter = 0;
1828 for (unsigned int i = 0; i < 5; ++i)
1829 for (unsigned int j = i; j < 5; ++j)
1830 helixError[counter++] = helixCovariance(i, j);
1831 }
1832
1833 auto trackFit = m_trackFitResults.appendNew(helixParam, helixError, pType, pValue, patternCdc.getInteger(),
1834 patternVxd.getInteger(), trk_fit.ndf());
1835 track->setTrackFitResultIndex(pType, trackFit->getArrayIndex());
1836 /*
1837 B2INFO("--- B1 Track: ");
1838 std::cout << "Momentum = " << momentum << std::endl;
1839 std::cout << "Position = " << position << std::endl;
1840 std::cout << "7x7 error matrix = " << error7x7 << std::endl;
1841 B2INFO("--- B2 Track: ");
1842 std::cout << "Momentum = " << std::endl;
1843 trackFit->get4Momentum().Print();
1844 std::cout << "Position = " << std::endl;
1845 trackFit->getPosition().Print();
1846 std::cout << "6x6 error matrix = " << std::endl;
1847 trackFit->getCovariance6().Print();
1848 TMatrixDSym b2Error7x7(7);
1849 fill7x7ErrorMatrix(trackFit, b2Error7x7, thisMass, 1.5);
1850 std::cout << "7x7 error matrix = " << std::endl;
1851 b2Error7x7.Print();
1852 */
1853 }
1854}
Belle2::B2BIIConvertMdstModule::ERRMCONV
const int ERRMCONV[7]
CONVERSION OF TRACK ERROR MATRIX ELEMENTS.
Definition: B2BIIConvertMdstModule.h:435
Belle2::B2BIIConvertMdstModule::event
virtual void event() override
Called for each event.
Definition: B2BIIConvertMdstModule.cc:295
Belle2::B2BIIConvertMdstModule::m_trackFitResults
StoreArray< TrackFitResult > m_trackFitResults
Track fir results.
Definition: B2BIIConvertMdstModule.h:393
Belle2::B2BIIConvertMdstModule::BFIELD
const double BFIELD
B filed in TESLA.
Definition: B2BIIConvertMdstModule.h:438
Belle2::B2BIIConvertMdstModule::m_belleTrkExtra
StoreArray< BelleTrkExtra > m_belleTrkExtra
Belle CDC extra information.
Definition: B2BIIConvertMdstModule.h:402
Belle2::B2BIIConvertMdstModule::getHelixParameters
int getHelixParameters(const Belle::Mdst_trk_fit &trk_fit, const double mass, const HepPoint3D &newPivot, std::vector< float > &helixParams, CLHEP::HepSymMatrix &error5x5, CLHEP::HepLorentzVector &momentum, HepPoint3D &position, CLHEP::HepSymMatrix &error7x7, const double dPhi=0.0)
Fills Helix parameters (converted to Belle II version), 5x5 error matrix, 4-momentum,...
Definition: B2BIIConvertMdstModule.cc:1542
Belle2::B2BIIConvertMdstModule::c_nHyp
static const int c_nHyp
Number of Belle track hypotheses (see c_belleHyp_to_chargedStable).
Definition: B2BIIConvertMdstModule.h:321
Belle2::HitPatternCDC
Hit pattern of CDC hits within a track.
Definition: HitPatternCDC.h:35
Belle2::HitPatternCDC::setNHits
void setNHits(unsigned short nHits)
Sets the 8 MSBs to the total number of hits in the CDC.
Definition: HitPatternCDC.cc:53
Belle2::HitPatternCDC::getInteger
ULong64_t getInteger() const
Getter for underlying integer type.
Definition: HitPatternCDC.cc:69
Belle2::HitPatternVXD
Hit pattern of the VXD within a track.
Definition: HitPatternVXD.h:37
Belle2::HitPatternVXD::getInteger
unsigned int getInteger() const
Getter for the underlying integer.
Definition: HitPatternVXD.h:58
Belle2::HitPatternVXD::setSVDLayer
void setSVDLayer(const unsigned short layerId, unsigned short uHits, unsigned short vHits)
Set the number of hits in a specific layer of the SVD.
Definition: HitPatternVXD.cc:51
Belle2::StoreObjPtr
Type-safe access to single objects in the data store.
Definition: StoreObjPtr.h:95
Belle2::UncertainHelix
This class represents an ideal helix in perigee parameterization including the covariance matrix of t...
Definition: UncertainHelix.h:36

◆ convertMdstChargedTable()

void convertMdstChargedTable ( )
private

Reads and converts all entries of Mdst_charged (Mdst_trk and Mdst_trk_fit) Panther table to Track (TrackFitResult) dataobjects and adds them to StoreArray<Track> (StoreArray<TrackFitResult>).

Definition at line 466 of file B2BIIConvertMdstModule.cc.

467{
468 // Relations
469 RelationArray tracksToMCParticles(m_tracks, m_mcParticles);
470
471 // Loop over all Belle charged tracks
472 Belle::Mdst_charged_Manager& m = Belle::Mdst_charged_Manager::get_manager();
473 for (Belle::Mdst_charged_Manager::iterator chargedIterator = m.begin(); chargedIterator != m.end(); ++chargedIterator) {
474 Belle::Mdst_charged belleTrack = *chargedIterator;
475
476 Belle::Mdst_trk& trk = belleTrack.trk();
477 bool foundValidTrackFitResult = false;
478 for (int mhyp = 0 ; mhyp < c_nHyp; ++mhyp) {
479 Belle::Mdst_trk_fit& trk_fit = trk.mhyp(mhyp);
480
481 double pValue = TMath::Prob(trk_fit.chisq(), trk_fit.ndf());
482 if (std::isnan(pValue)) continue;
483 foundValidTrackFitResult = true;
484 }
485 if (not foundValidTrackFitResult) continue;
486
487 auto track = m_tracks.appendNew();
488
489 // convert MDST_Charged -> Track
490 convertMdstChargedObject(belleTrack, track);
491
492 convertPIDData(belleTrack, track);
493
494 if (m_realData)
495 continue;
496
497 // create Track -> MCParticle relation
498 // step 1: MDSTCharged -> Gen_hepevt
499 const Belle::Gen_hepevt& hep0 = get_hepevt(belleTrack);
500 if (hep0 == 0)
501 continue;
502 const Belle::Gen_hepevt* hep = nullptr;
503 switch (m_mcMatchingMode) {
504 case c_Direct:
505 hep = &hep0;
506 break;
507 case c_GeneratorLevel:
508 hep = &gen_level(hep0);
509 break;
510 }
511 // step 2: Gen_hepevt -> MCParticle
512 if (genHepevtToMCParticle.count(hep->get_ID()) > 0) {
513 int matchedMCParticle = genHepevtToMCParticle[hep->get_ID()];
514
515 // step 3: set the relation
516 tracksToMCParticles.add(track->getArrayIndex(), matchedMCParticle);
517
518 testMCRelation(*hep, m_mcParticles[matchedMCParticle], "Track");
519 } else {
520 B2DEBUG(99, "Can not find MCParticle corresponding to this gen_hepevt (Panther ID = " << hep->get_ID() << ")");
521 B2DEBUG(99, "Gen_hepevt: Panther ID = " << hep->get_ID() << "; idhep = " << hep->idhep() << "; isthep = " << hep->isthep());
522 }
523 }
524}
Belle2::B2BIIConvertMdstModule::testMCRelation
void testMCRelation(const Belle::Gen_hepevt &belleMC, const MCParticle *mcP, const std::string &objectName)
Checks if the reconstructed object (Track, ECLCluster, ...) was matched to the same MCParticle.
Definition: B2BIIConvertMdstModule.cc:2145
Belle2::B2BIIConvertMdstModule::convertPIDData
void convertPIDData(const Belle::Mdst_charged &belleTrack, const Track *track)
Get PID information for belleTrack and add it to PIDLikelihood (with relation from track).
Definition: B2BIIConvertMdstModule.cc:1414
Belle2::B2BIIConvertMdstModule::convertMdstChargedObject
void convertMdstChargedObject(const Belle::Mdst_charged &belleTrack, Track *track)
Converts Mdst_charged (Mdst_trk(_fit)) record to Track (TrackFitResult) object.
Definition: B2BIIConvertMdstModule.cc:1694
Belle2::B2BIIConvertMdstModule::m_mcParticles
StoreArray< MCParticle > m_mcParticles
MC particles.
Definition: B2BIIConvertMdstModule.h:381
Belle2::RelationArray
Low-level class to create/modify relations between StoreArrays.
Definition: RelationArray.h:62

◆ convertMdstECLObject()

void convertMdstECLObject ( const Belle::Mdst_ecl &  ecl,
const Belle::Mdst_ecl_aux &  eclAux,
ECLCluster *  eclCluster 
)
private

Converts Mdst_ecl(_aux) record to ECLCluster object.

If running on MC, the ECLCluster -> MCParticle relation is set as well.

Definition at line 1896 of file B2BIIConvertMdstModule.cc.

1898{
1899 if (eclAux.e9oe25() < m_matchType2E9oE25Threshold)
1900 eclCluster->setIsTrack(ecl.match() > 0);
1901 else
1902 eclCluster->setIsTrack(ecl.match() == 1);
1903
1904 eclCluster->setEnergy(ecl.energy()); //must happen before setCovarianceMatrix()!
1905 if (eclCluster->getEnergy(ECLCluster::EHypothesisBit::c_nPhotons) > 0.5)
1906 eclCluster->addHypothesis(ECLCluster::EHypothesisBit::c_neutralHadron);
1907 eclCluster->setPhi(ecl.phi());
1908 eclCluster->setTheta(ecl.theta());
1909 eclCluster->setR(ecl.r());
1910 eclCluster->setdeltaL(ecl.quality());
1911
1912 double covarianceMatrix[6];
1913 covarianceMatrix[0] = ecl.error(0); // error on energy
1914 covarianceMatrix[1] = ecl.error(1);
1915 covarianceMatrix[2] = ecl.error(2); // error on phi
1916 covarianceMatrix[3] = ecl.error(3);
1917 covarianceMatrix[4] = ecl.error(4);
1918 covarianceMatrix[5] = ecl.error(5); // error on theta
1919 eclCluster->setCovarianceMatrix(covarianceMatrix);
1920
1921 eclCluster->setLAT(eclAux.width());
1922 eclCluster->setEnergyRaw(eclAux.mass());
1923 eclCluster->setE9oE21(eclAux.e9oe25());
1924 eclCluster->setEnergyHighestCrystal(eclAux.seed());
1925 // The property 2 of eclAux contains the timing information
1926 // in a bit encoded format.
1927 // The 16 bits: 0-15 contain tdc count
1928 float prop2 = eclAux.property(2);
1929 // a float to int conversion
1930 int property2;
1931 std::memcpy(&property2, &prop2, sizeof(int));
1932 //decode the bit encoded variables
1933 int tdccount;
1934 tdccount = property2 & 0xffff;
1935 eclCluster->setTime(tdccount);
1936 eclCluster->setNumberOfCrystals(eclAux.nhits());
1937 double dist = computeTrkMinDistanceBelle(eclCluster);
1938 eclCluster->setMinTrkDistance(dist);
1939}
Belle2::B2BIIConvertMdstModule::computeTrkMinDistanceBelle
double computeTrkMinDistanceBelle(ECLCluster *eclCluster)
calculate the minimal distance between a cluster and a set of tracks on the ECL surface.
Definition: B2BIIConvertMdstModule.cc:1941
Belle2::ECLCluster::setE9oE21
void setE9oE21(double E9oE21)
Set E9/E21 energy ratio.
Definition: ECLCluster.h:193
Belle2::ECLCluster::addHypothesis
void addHypothesis(EHypothesisBit bitmask)
Add bitmask to current hypothesis.
Definition: ECLCluster.h:129
Belle2::ECLCluster::setTheta
void setTheta(double theta)
Set Theta of Shower (radian).
Definition: ECLCluster.h:217
Belle2::ECLCluster::setPhi
void setPhi(double phi)
Set Phi of Shower (radian).
Definition: ECLCluster.h:220
Belle2::ECLCluster::setdeltaL
void setdeltaL(double deltaL)
Set deltaL for shower shape.
Definition: ECLCluster.h:172
Belle2::ECLCluster::setEnergyHighestCrystal
void setEnergyHighestCrystal(double energyhighestcrystal)
Set energy of highest energetic crystal (GeV).
Definition: ECLCluster.h:232
Belle2::ECLCluster::setEnergyRaw
void setEnergyRaw(double energyraw)
Set Incorrect Energy deposited (GeV).
Definition: ECLCluster.h:229
Belle2::ECLCluster::setTime
void setTime(double time)
Set time information.
Definition: ECLCluster.h:211
Belle2::ECLCluster::setCovarianceMatrix
void setCovarianceMatrix(double covArray[6])
Set covariance matrix (3x3), i.e.
Definition: ECLCluster.h:152
Belle2::ECLCluster::setMinTrkDistance
void setMinTrkDistance(double distance)
Set distance between cluster COG and track extrapolation to ECL.
Definition: ECLCluster.h:175
Belle2::ECLCluster::setLAT
void setLAT(double LAT)
Set Lateral distribution parameter.
Definition: ECLCluster.h:205
Belle2::ECLCluster::getEnergy
double getEnergy(EHypothesisBit hypothesis) const
Return Energy (GeV).
Definition: ECLCluster.cc:23
Belle2::ECLCluster::setNumberOfCrystals
void setNumberOfCrystals(double noc)
Set number of crystals (sum of weights).
Definition: ECLCluster.h:208
Belle2::ECLCluster::setEnergy
void setEnergy(double energy)
Set Corrected Energy (GeV).
Definition: ECLCluster.h:226
Belle2::ECLCluster::setIsTrack
void setIsTrack(bool istrack)
Set m_isTrack true if the cluster matches with a track.
Definition: ECLCluster.h:104
Belle2::ECLCluster::EHypothesisBit::c_nPhotons
@ c_nPhotons
CR is split into n photons (N1)
Belle2::ECLCluster::EHypothesisBit::c_neutralHadron
@ c_neutralHadron
CR is reconstructed as a neutral hadron (N2)
Belle2::ECLCluster::setR
void setR(double r)
Set R (in cm).
Definition: ECLCluster.h:223

◆ convertMdstECLTable()

void convertMdstECLTable ( )
private

Reads and converts all entries of Mdst_ecl(_aux) Panther table to ECLCluster dataobjects and adds them to StoreArray<ECLCluster>.

Definition at line 973 of file B2BIIConvertMdstModule.cc.

974{
975 // Relations
976 RelationArray eclClustersToMCParticles(m_eclClusters, m_mcParticles);
977
978 // Clear the mdstEcl <-> ECLCluster map
979 mdstEclToECLCluster.clear();
980
981 // Loop over all Belle Mdst_ecl
982 Belle::Mdst_ecl_Manager& ecl_manager = Belle::Mdst_ecl_Manager::get_manager();
983 Belle::Mdst_ecl_aux_Manager& ecl_aux_manager = Belle::Mdst_ecl_aux_Manager::get_manager();
984
985 for (Belle::Mdst_ecl_Manager::iterator eclIterator = ecl_manager.begin(); eclIterator != ecl_manager.end(); ++eclIterator) {
986
987 // Pull Mdst_ecl from manager
988 Belle::Mdst_ecl mdstEcl = *eclIterator;
989 Belle::Mdst_ecl_aux mdstEclAux(ecl_aux_manager(mdstEcl.get_ID()));
990
991 // Create Belle II ECLCluster
992 auto B2EclCluster = m_eclClusters.appendNew();
993
994 // Convert Mdst_ecl -> ECLCluster and create map of indices
995 convertMdstECLObject(mdstEcl, mdstEclAux, B2EclCluster);
996 mdstEclToECLCluster[mdstEcl.get_ID()] = B2EclCluster->getArrayIndex();
997
998 // set ConnectedRegionID and ClusterID to
999 // cluster's array index + 1 and 1, respectively
1000 B2EclCluster->setConnectedRegionId(B2EclCluster->getArrayIndex() + 1);
1001 B2EclCluster->setClusterId(1);
1002
1003 if (m_realData)
1004 continue;
1005
1006 // Create ECLCluster -> MCParticle relation
1007 // Step 1: MDST_ECL -> Gen_hepevt
1008 const Belle::Gen_hepevt& hep0 = get_hepevt(mdstEcl);
1009 if (hep0 == 0)
1010 continue;
1011 const Belle::Gen_hepevt* hep = nullptr;
1012 switch (m_mcMatchingMode) {
1013 case c_Direct:
1014 hep = &hep0;
1015 break;
1016 case c_GeneratorLevel:
1017 hep = &gen_level(hep0);
1018 break;
1019 }
1020 // Step 2: Gen_hepevt -> MCParticle
1021 if (genHepevtToMCParticle.count(hep->get_ID()) > 0) {
1022 int matchedMCParticleID = genHepevtToMCParticle[hep->get_ID()];
1023 // Step 3: set the relation
1024 eclClustersToMCParticles.add(B2EclCluster->getArrayIndex(), matchedMCParticleID);
1025 testMCRelation(*hep, m_mcParticles[matchedMCParticleID], "ECLCluster");
1026 } else {
1027 B2DEBUG(79, "Cannot find MCParticle corresponding to this gen_hepevt (Panther ID = " << hep->get_ID() << ")");
1028 B2DEBUG(79, "Gen_hepevt: Panther ID = " << hep->get_ID() << "; idhep = " << hep->idhep() << "; isthep = " << hep->isthep());
1029 }
1030 }
1031}
Belle2::B2BIIConvertMdstModule::convertMdstECLObject
void convertMdstECLObject(const Belle::Mdst_ecl &ecl, const Belle::Mdst_ecl_aux &eclAux, ECLCluster *eclCluster)
Converts Mdst_ecl(_aux) record to ECLCluster object.
Definition: B2BIIConvertMdstModule.cc:1896
Belle2::B2BIIConvertMdstModule::m_eclClusters
StoreArray< ECLCluster > m_eclClusters
ECL clusters.
Definition: B2BIIConvertMdstModule.h:384
Belle2::B2BIIConvertMdstModule::mdstEclToECLCluster
std::map< int, int > mdstEclToECLCluster
map of Mdst_ecl Panther IDs and corresponding ECLCluster StoreArray indices
Definition: B2BIIConvertMdstModule.h:372

◆ convertMdstGammaTable()

void convertMdstGammaTable ( )
private

Reads all entries of Mdst_Gamma Panther table, creates a particle list 'gamma:mdst' and adds them to StoreArray<Particles>.

Definition at line 1060 of file B2BIIConvertMdstModule.cc.

1061{
1062 // Relations
1063 RelationArray particlesToMCParticles(m_particles, m_mcParticles);
1064
1065 // Clear the mdstGamma <-> Particle map
1066 mdstGammaToParticle.clear();
1067
1068 // Create and initialize particle list
1069 StoreObjPtr<ParticleList> plist("gamma:mdst");
1070 plist.create();
1071 plist->initialize(22, "gamma:mdst");
1072
1073 // Loop over all Belle Mdst_gamma
1074 Belle::Mdst_gamma_Manager& gamma_manager = Belle::Mdst_gamma_Manager::get_manager();
1075
1076 for (Belle::Mdst_gamma_Manager::iterator gammaIterator = gamma_manager.begin(); gammaIterator != gamma_manager.end();
1077 ++gammaIterator) {
1078
1079 // Pull Mdst_gamma from manager and Mdst_ecl from pointer to Mdst_ecl
1080 Belle::Mdst_gamma mdstGamma = *gammaIterator;
1081 Belle::Mdst_ecl mdstEcl = mdstGamma.ecl();
1082 if (!mdstEcl)
1083 continue;
1084
1085 // Get ECLCluster from map
1086 ECLCluster* B2EclCluster = m_eclClusters[mdstEclToECLCluster[mdstEcl.get_ID()]];
1087 if (!B2EclCluster)
1088 continue;
1089
1090 // Create Particle from ECLCluster, add to StoreArray, create gamma map entry
1091 Particle* B2Gamma = m_particles.appendNew(B2EclCluster);
1092 mdstGammaToParticle[mdstGamma.get_ID()] = B2Gamma->getArrayIndex();
1093
1094 // Add particle to particle list
1095 plist->addParticle(B2Gamma);
1096
1097 if (m_realData)
1098 continue;
1099
1100 // Relation to MCParticle
1101 MCParticle* matchedMCParticle = B2EclCluster->getRelated<MCParticle>();
1102 if (matchedMCParticle)
1103 B2Gamma->addRelationTo(matchedMCParticle);
1104 }
1105}
Belle2::B2BIIConvertMdstModule::m_particles
StoreArray< Particle > m_particles
Particles.
Definition: B2BIIConvertMdstModule.h:399
Belle2::B2BIIConvertMdstModule::mdstGammaToParticle
std::map< int, int > mdstGammaToParticle
map of gamma Panther IDs and corresponding Particle StoreArray indices
Definition: B2BIIConvertMdstModule.h:374
Belle2::ECLCluster
ECL cluster data.
Definition: ECLCluster.h:27
Belle2::MCParticle
A Class to store the Monte Carlo particle information.
Definition: MCParticle.h:32
Belle2::Particle
Class to store reconstructed particles.
Definition: Particle.h:76
Belle2::RelationsInterface::addRelationTo
void addRelationTo(const RelationsInterface< BASE > *object, float weight=1.0, const std::string &namedRelation="") const
Add a relation from this object to another object (with caching).
Definition: RelationsObject.h:142
Belle2::RelationsInterface::getArrayIndex
int getArrayIndex() const
Returns this object's array index (in StoreArray), or -1 if not found.
Definition: RelationsObject.h:385
Belle2::RelationsInterface::getRelated
T * getRelated(const std::string &name="", const std::string &namedRelation="") const
Get the object to or from which this object has a relation.
Definition: RelationsObject.h:278
Belle2::StoreArray::appendNew
T * appendNew()
Construct a new T object at the end of the array.
Definition: StoreArray.h:246

◆ convertMdstKLMObject()

void convertMdstKLMObject ( const Belle::Mdst_klm_cluster &  klm,
KLMCluster *  klmCluster 
)
private

Converts Mdst_klm_cluster record to KLMCluster object.

No MCRelation is set, since there was no matching in Belle.

Definition at line 1997 of file B2BIIConvertMdstModule.cc.

1998{
1999 // note: Belle quality flag is not saved (no free int variable in Belle2 KLMCluster)
2000 klmCluster->setLayers(klm_cluster.layers());
2001 klmCluster->setInnermostLayer(klm_cluster.first_layer());
2002}
Belle2::KLMCluster::setLayers
void setLayers(int layers)
Set number of layers with hits.
Definition: KLMCluster.h:145
Belle2::KLMCluster::setInnermostLayer
void setInnermostLayer(int innermostLayer)
Set number of the innermost layer with hits.
Definition: KLMCluster.h:152

◆ convertMdstKLMTable()

void convertMdstKLMTable ( )
private

Reads and converts all entries of Mdst_klm_cluster Panther table to KLMCluster dataobjects and adds them to StoreArray<KLMCluster>.

Definition at line 1033 of file B2BIIConvertMdstModule.cc.

1034{
1035 // There was no MC matching in Belle for KLM Clusters
1036
1037 // Clear the mdstKlm <-> KLMCluster map
1038 mdstKlmToKLMCluster.clear();
1039
1040 // Loop over all Belle Mdst_klm_cluster
1041 Belle::Mdst_klm_cluster_Manager& klm_cluster_manager = Belle::Mdst_klm_cluster_Manager::get_manager();
1042
1043 for (Belle::Mdst_klm_cluster_Manager::iterator klmC_Ite = klm_cluster_manager.begin(); klmC_Ite != klm_cluster_manager.end();
1044 ++klmC_Ite) {
1045
1046 // Pull Mdst_ecl from manager
1047 Belle::Mdst_klm_cluster mdstKlm_cluster = *klmC_Ite;
1048
1049 // Create Belle II ECLCluster
1050 auto B2KlmCluster = m_klmClusters.appendNew();
1051
1052 // Convert Mdst_klm_cluster -> KLMCluster and create map of indices
1053 convertMdstKLMObject(mdstKlm_cluster, B2KlmCluster);
1054 mdstKlmToKLMCluster[mdstKlm_cluster.get_ID()] = B2KlmCluster->getArrayIndex();
1055
1056 }
1057}
Belle2::B2BIIConvertMdstModule::mdstKlmToKLMCluster
std::map< int, int > mdstKlmToKLMCluster
map of Mdst_klm Panther IDs and corresponding KLMCluster StoreArray indices
Definition: B2BIIConvertMdstModule.h:376
Belle2::B2BIIConvertMdstModule::m_klmClusters
StoreArray< KLMCluster > m_klmClusters
KLM clusters.
Definition: B2BIIConvertMdstModule.h:387
Belle2::B2BIIConvertMdstModule::convertMdstKLMObject
void convertMdstKLMObject(const Belle::Mdst_klm_cluster &klm, KLMCluster *klmCluster)
Converts Mdst_klm_cluster record to KLMCluster object.
Definition: B2BIIConvertMdstModule.cc:1997

◆ convertMdstKLongTable()

void convertMdstKLongTable ( )
private

Reads all entries of Mdst_Klong Panther table, creates a particle list 'K_L0:mdst' and adds them to StoreArray<Particles>.

Definition at line 1181 of file B2BIIConvertMdstModule.cc.

1182{
1183 // Relations
1184 RelationArray particlesToMCParticles(m_particles, m_mcParticles);
1185
1186
1187 // Create and initialize particle list
1188 StoreObjPtr<ParticleList> plist("K_L0:mdst");
1189 plist.create();
1190 plist->initialize(Const::Klong.getPDGCode(), "K_L0:mdst");
1191
1192 Belle::Mdst_klong_Manager& klong_manager = Belle::Mdst_klong_Manager::get_manager();
1193 for (Belle::Mdst_klong_Manager::iterator klong_Ite = klong_manager.begin(); klong_Ite != klong_manager.end(); ++klong_Ite) {
1194
1195 // Pull Mdst_klong from manager and Mdst_klm from pointer to Mdst_klm
1196 Belle::Mdst_klong mdstKlong = *klong_Ite;
1197 Belle::Mdst_klm_cluster mdstKlm = mdstKlong.klmc();
1198
1199 if (!mdstKlm)
1200 continue;
1201
1202
1203 // Get KLMCluster from map
1204 KLMCluster* B2KlmCluster = m_klmClusters[mdstKlmToKLMCluster[mdstKlm.get_ID()]];
1205 if (!B2KlmCluster)
1206 continue;
1207
1208 // Extract cluster position from Klong and save it in KLMCluster
1209 B2KlmCluster->setClusterPosition(mdstKlong.cos_x(), mdstKlong.cos_y(), mdstKlong.cos_z());
1210
1211 // Create Particle from KLMCluster, add to StoreArray, create Klong map entry
1212 Particle* B2Klong = m_particles.appendNew(B2KlmCluster);
1213 mdstKlongToParticle[mdstKlong.get_ID()] = B2Klong->getArrayIndex();
1214
1215 // Add particle to particle list
1216 plist->addParticle(B2Klong);
1217 }
1218
1219 // (Vague) MC Matching
1220 // There was no MC matching for KLongs in Belle , but a hack:
1221 // Check if MC KLong and reconstructed KLong (only without ecl) are within 15 degree for phi and theta, we set a relation
1222 // for the best reconstructed KLong to the MC KLong.
1223 // Taken and adapted from http://belle.kek.jp/secured/wiki/doku.php?id=physics:ckm:kleff
1224
1225 if (!m_realData) {
1226
1227 Belle::Gen_hepevt_Manager& GenMgr = Belle::Gen_hepevt_Manager::get_manager();
1228 const double dang(15. / 180.*M_PI); // check reconstructed candidates within 15 degrees
1229
1230 for (Belle::Gen_hepevt_Manager::iterator klong_hep_it = GenMgr.begin(); klong_hep_it != GenMgr.end(); ++klong_hep_it) {
1231
1232 if (abs((*klong_hep_it).idhep()) == Const::Klong.getPDGCode() && klong_hep_it->isthep() > 0) {
1233
1234 CLHEP::HepLorentzVector gp4(klong_hep_it->PX(), klong_hep_it->PY(), klong_hep_it->PZ(), klong_hep_it->E());
1235 double sum(0.0);
1236 int bestRecKlongID(0);
1237
1238 for (Belle::Mdst_klong_Manager::iterator klong_rec_it = klong_manager.begin(); klong_rec_it != klong_manager.end();
1239 ++klong_rec_it) {
1240
1241 // if((*klong_rec_it).klmc().ecl())continue; // check only klm cand.
1242 if ((*klong_rec_it).ecl())
1243 continue; // check only klm cand.
1244 CLHEP::Hep3Vector klp3(klong_rec_it->cos_x(), klong_rec_it->cos_y(), klong_rec_it->cos_z());
1245
1246 if (cos(gp4.theta() - klp3.theta()) > cos(dang) && cos(gp4.phi() - klp3.phi()) > cos(dang)) {
1247
1248 double tmp_sum = cos(gp4.theta() - klp3.theta()) + cos(gp4.phi() - klp3.phi());
1249 if (tmp_sum > sum) {
1250 bestRecKlongID = mdstKlongToParticle[(*klong_rec_it).get_ID()];
1251 sum = tmp_sum;
1252 }
1253 }
1254
1255 }
1256 if (sum > 0.0) {
1257 int matchedMCParticleID = genHepevtToMCParticle[(*klong_hep_it).get_ID()];
1258 particlesToMCParticles.add(bestRecKlongID, matchedMCParticleID);
1259 testMCRelation((*klong_hep_it), m_mcParticles[matchedMCParticleID], "m_particles");
1260 }
1261 }
1262 }
1263 }
1264}
Belle2::B2BIIConvertMdstModule::mdstKlongToParticle
std::map< int, int > mdstKlongToParticle
map of Klong Panther IDs and corresponding Particle StoreArray indices
Definition: B2BIIConvertMdstModule.h:378
Belle2::Const::ParticleType::getPDGCode
int getPDGCode() const
PDG code.
Definition: Const.h:473
Belle2::Const::Klong
static const ParticleType Klong
K^0_L particle.
Definition: Const.h:678
Belle2::KLMCluster
KLM cluster data.
Definition: KLMCluster.h:28
Belle2::KLMCluster::setClusterPosition
void setClusterPosition(float globalX, float globalY, float globalZ)
Set global position.
Definition: KLMCluster.h:161

◆ convertMdstPi0Table()

void convertMdstPi0Table ( )
private

Reads all entries of Mdst_Pi0 Panther table, creates a particle list 'pi0:mdst' and adds them to StoreArray<Particles>.

Definition at line 1133 of file B2BIIConvertMdstModule.cc.

1134{
1135 // Create and initialize particle list
1136 StoreObjPtr<ParticleList> plist("pi0:mdst");
1137 plist.create();
1138 plist->initialize(111, "pi0:mdst");
1139
1140 // Loop over all Mdst_pi0
1141 Belle::Mdst_pi0_Manager& pi0_manager = Belle::Mdst_pi0_Manager::get_manager();
1142 for (Belle::Mdst_pi0_Manager::iterator pi0Iterator = pi0_manager.begin(); pi0Iterator != pi0_manager.end(); ++pi0Iterator) {
1143
1144 // Pull Mdst_pi0 from manager and Mdst_gammas from pointers to Mdst_gammas
1145 Belle::Mdst_pi0 mdstPi0 = *pi0Iterator;
1146 Belle::Mdst_gamma mdstGamma1 = mdstPi0.gamma(0);
1147 Belle::Mdst_gamma mdstGamma2 = mdstPi0.gamma(1);
1148 if (!mdstGamma1 || !mdstGamma2)
1149 continue;
1150
1151 ROOT::Math::PxPyPzEVector p4(mdstPi0.px(), mdstPi0.py(), mdstPi0.pz(), mdstPi0.energy());
1152
1153 // Create Particle from LorentzVector and PDG code, add to StoreArray
1154 Particle* B2Pi0 = m_particles.appendNew(p4, 111);
1155
1156 // Get Belle II photons from map
1157 Particle* B2Gamma1 = m_particles[mdstGammaToParticle[mdstGamma1.get_ID()]];
1158 Particle* B2Gamma2 = m_particles[mdstGammaToParticle[mdstGamma2.get_ID()]];
1159 if (!B2Gamma1 || !B2Gamma2)
1160 continue;
1161
1162 // Append photons as pi0 daughters
1163 B2Pi0->appendDaughter(B2Gamma1);
1164 B2Pi0->appendDaughter(B2Gamma2);
1165
1166 // Add chisq of mass-constrained Kfit
1167 B2Pi0->addExtraInfo("chiSquared", mdstPi0.chisq());
1168
1169 // Ndf of a pi0 mass-constrained kinematic fit is 1
1170 B2Pi0->addExtraInfo("ndf", 1);
1171
1172 // Add p-value to extra Info
1173 double prob = TMath::Prob(mdstPi0.chisq(), 1);
1174 B2Pi0->setPValue(prob);
1175
1176 // Add particle to particle list
1177 plist->addParticle(B2Pi0);
1178 }
1179}
Belle2::Particle::appendDaughter
void appendDaughter(const Particle *daughter, const bool updateType=true, const int daughterProperty=c_Ordinary)
Appends index of daughter to daughters index array.
Definition: Particle.cc:707
Belle2::Particle::addExtraInfo
void addExtraInfo(const std::string &name, double value)
Sets the user-defined data of given name to the given value.
Definition: Particle.cc:1421
Belle2::Particle::setPValue
void setPValue(double pValue)
Sets chi^2 probability of fit.
Definition: Particle.h:377

◆ convertMdstVee2Table()

void convertMdstVee2Table ( )
private

Reads and converts all entries of Mdst_vee2 Panther table to V0 dataobjects and adds them to StoreArray<V0>.

V0 daughters are converted from Mdst_Vee_Daughters to TrackFitResults, which are appended to the TrackFitResult StoreArray.

Definition at line 526 of file B2BIIConvertMdstModule.cc.

527{
528 // Create and initialize K_S0 particle list
529 StoreObjPtr<ParticleList> ksPList("K_S0:mdst");
530 ksPList.create();
531 ksPList->initialize(310, ksPList.getName());
532
533 // Create and initialize Lambda0 and anti-Lamda0 particle list
534 StoreObjPtr<ParticleList> lambda0PList("Lambda0:mdst");
535 lambda0PList.create();
536 lambda0PList->initialize(3122, lambda0PList.getName());
537
538 StoreObjPtr<ParticleList> antiLambda0PList("anti-Lambda0:mdst");
539 antiLambda0PList.create();
540 antiLambda0PList->initialize(-3122, antiLambda0PList.getName());
541
542 antiLambda0PList->bindAntiParticleList(*lambda0PList);
543
544 // Create and initialize converted gamma particle list
545 StoreObjPtr<ParticleList> convGammaPList("gamma:v0mdst");
546 convGammaPList.create();
547 convGammaPList->initialize(22, convGammaPList.getName());
548
549 // Loop over all Belle Vee2 candidates
550 Belle::Mdst_vee2_Manager& m = Belle::Mdst_vee2_Manager::get_manager();
551 for (Belle::Mdst_vee2_Manager::iterator vee2Iterator = m.begin(); vee2Iterator != m.end(); ++vee2Iterator) {
552 Belle::Mdst_vee2 belleV0 = *vee2Iterator;
553
554 // +ve track
555 Belle::Mdst_charged belleTrackP = belleV0.chgd(0);
556 // -ve track
557 Belle::Mdst_charged belleTrackM = belleV0.chgd(1);
558
559 // type of V0
560 Const::ChargedStable pTypeP(Const::pion);
561 Const::ChargedStable pTypeM(Const::pion);
562 int belleHypP = -1;
563 int belleHypM = -1;
564
565 switch (belleV0.kind()) {
566 case 1 : // K0s -> pi+ pi-
567 pTypeP = Const::pion;
568 pTypeM = Const::pion;
569 belleHypP = 2;
570 belleHypM = 2;
571 break;
572 case 2 : // Lambda -> p+ pi-
573 pTypeP = Const::proton;
574 pTypeM = Const::pion;
575 belleHypP = 4;
576 belleHypM = 2;
577 break;
578 case 3 : // anti-Lambda -> pi+ anti-p-
579 pTypeP = Const::pion;
580 pTypeM = Const::proton;
581 belleHypP = 2;
582 belleHypM = 4;
583 break;
584 case 4 : // gamma -> e+ e-
585 pTypeP = Const::electron;
586 pTypeM = Const::electron;
587 belleHypP = 0;
588 belleHypM = 0;
589 break;
590 default :
591 B2WARNING("Conversion of vee2 candidate of unknown kind! kind = " << belleV0.kind());
592 }
593
594 // This part is copied from Relation.cc in BASF
595 int trackID[2] = {0, 0};
596 unsigned nTrack = 0;
597 Belle::Mdst_charged_Manager& charged_mag = Belle::Mdst_charged_Manager::get_manager();
598 for (std::vector<Belle::Mdst_charged>::iterator chgIterator = charged_mag.begin(); chgIterator != charged_mag.end();
599 ++chgIterator) {
600 if (belleV0.chgd(0).get_ID() >= 1 && trackID[0] == 0 && belleV0.chgd(0).get_ID() == chgIterator->get_ID()) {
601 trackID[0] = (int)(chgIterator->get_ID()); //+ve trac
602 ++nTrack;
603 }
604 if (belleV0.chgd(1).get_ID() >= 1 && trackID[1] == 0 && belleV0.chgd(1).get_ID() == chgIterator->get_ID()) {
605 trackID[1] = (int)(chgIterator->get_ID()); //-ve trac
606 ++nTrack;
607 }
608 if (nTrack == 2)
609 break;
610 }
611 if (std::max(trackID[0], trackID[1]) > m_tracks.getEntries()) continue;
612
613 HepPoint3D dauPivot(belleV0.vx(), belleV0.vy(), belleV0.vz());
614 int trackFitPIndex = -1;
615 int trackFitMIndex = -1;
616 Particle daughterP, daughterM;
617 CLHEP::HepLorentzVector momentumP;
618 CLHEP::HepSymMatrix error7x7P(7, 0);
619 HepPoint3D positionP;
620 TMatrixFSym errMatrixP(7);
621 CLHEP::HepLorentzVector momentumM;
622 CLHEP::HepSymMatrix error7x7M(7, 0);
623 HepPoint3D positionM;
624 TMatrixFSym errMatrixM(7);
625 CLHEP::HepSymMatrix error5x5(5, 0);
626 if (trackID[0] >= 1) {
627 if (belleV0.daut()) {
628 std::vector<float> helixParam(5);
629 std::vector<float> helixError(15);
630 belleVeeDaughterHelix(belleV0, 1, helixParam, helixError);
631
632 auto trackFitP = m_trackFitResults.appendNew(helixParam, helixError, pTypeP, 0.5, -1, -1, 0);
633 trackFitPIndex = trackFitP->getArrayIndex();
634
635 belleVeeDaughterToCartesian(belleV0, 1, pTypeP, momentumP, positionP, error7x7P);
636 TrackFitResult* tmpTFR = new TrackFitResult(createTrackFitResult(momentumP, positionP, error7x7P, 1, pTypeP, 0.5, -1, -1, 0));
637 // TrackFitResult internally stores helix parameters at pivot = (0,0,0) so the momentum of the Particle will be wrong again.
638 // Overwrite it.
639
640 for (unsigned i = 0; i < 7; i++)
641 for (unsigned j = 0; j < 7; j++)
642 errMatrixP(i, j) = error7x7P[i][j];
643
644 daughterP = Particle(trackID[0] - 1, tmpTFR, pTypeP);
645 daughterP.updateMomentum(ROOT::Math::PxPyPzEVector(momentumP.px(), momentumP.py(), momentumP.pz(), momentumP.e()),
646 ROOT::Math::XYZVector(positionP.x(), positionP.y(), positionP.z()),
647 errMatrixP, 0.5);
648 delete tmpTFR;
649 } else {
650 Belle::Mdst_trk_fit& trk_fit = charged_mag[trackID[0] - 1].trk().mhyp(belleHypP);
651 double pValue = TMath::Prob(trk_fit.chisq(), trk_fit.ndf());
652 if (std::isnan(pValue)) continue;
653
654 std::vector<float> helixParam(5);
655 std::vector<float> helixError(15);
656 convertHelix(trk_fit, HepPoint3D(0., 0., 0.), helixParam, helixError);
657
658 // Checking for invalid helix curvature with parameter 2 equal to 0:
659 if (helixParam[2] == 0) {
660 B2WARNING("Helix parameter for curvature == 0. Skipping Track! The parameter is: " << helixParam[2] << "...");
661 continue;
662 }
663
664 auto trackFitP = m_trackFitResults.appendNew(helixParam, helixError, pTypeP, pValue, -1, -1, 0);
665
666 trackFitPIndex = trackFitP->getArrayIndex();
667
668 daughterP = Particle(trackID[0] - 1, trackFitP, pTypeP);
669 // set momentum/positions at pivot = V0 decay vertex
670 getHelixParameters(trk_fit, pTypeP.getMass(), dauPivot,
671 helixParam, error5x5,
672 momentumP, positionP, error7x7P);
673
674 for (unsigned i = 0; i < 7; i++)
675 for (unsigned j = 0; j < 7; j++)
676 errMatrixP(i, j) = error7x7P[i][j];
677
678 daughterP.updateMomentum(ROOT::Math::PxPyPzEVector(momentumP.px(), momentumP.py(), momentumP.pz(), momentumP.e()),
679 ROOT::Math::XYZVector(positionP.x(), positionP.y(), positionP.z()),
680 errMatrixP, pValue);
681 }
682 }
683 if (trackID[1] >= 1) {
684 if (belleV0.daut()) {
685 std::vector<float> helixParam(5);
686 std::vector<float> helixError(15);
687 belleVeeDaughterHelix(belleV0, -1, helixParam, helixError);
688
689 auto trackFitM = m_trackFitResults.appendNew(helixParam, helixError, pTypeM, 0.5, -1, -1, 0);
690 trackFitMIndex = trackFitM->getArrayIndex();
691
692 belleVeeDaughterToCartesian(belleV0, -1, pTypeM, momentumM, positionM, error7x7M);
693 TrackFitResult* tmpTFR = new TrackFitResult(createTrackFitResult(momentumM, positionM, error7x7M, -1, pTypeM, 0.5, -1, -1, 0));
694 // TrackFitResult internally stores helix parameters at pivot = (0,0,0) so the momentum of the Particle will be wrong again.
695 // Overwrite it.
696 for (unsigned i = 0; i < 7; i++)
697 for (unsigned j = 0; j < 7; j++)
698 errMatrixM(i, j) = error7x7M[i][j];
699
700 daughterM = Particle(trackID[1] - 1, tmpTFR, pTypeM);
701 daughterM.updateMomentum(ROOT::Math::PxPyPzEVector(momentumM.px(), momentumM.py(), momentumM.pz(), momentumM.e()),
702 ROOT::Math::XYZVector(positionM.x(), positionM.y(), positionM.z()),
703 errMatrixM, 0.5);
704 delete tmpTFR;
705 } else {
706 Belle::Mdst_trk_fit& trk_fit = charged_mag[trackID[1] - 1].trk().mhyp(belleHypM);
707 double pValue = TMath::Prob(trk_fit.chisq(), trk_fit.ndf());
708 if (std::isnan(pValue)) continue;
709
710 std::vector<float> helixParam(5);
711 std::vector<float> helixError(15);
712 convertHelix(trk_fit, HepPoint3D(0., 0., 0.), helixParam, helixError);
713
714 // Checking for invalid helix curvature with parameter 2 equal to 0:
715 if (helixParam[2] == 0) {
716 B2WARNING("Helix parameter for curvature == 0. Skipping Track! The parameter is: " << helixParam[2] << "...");
717 continue;
718 }
719
720 auto trackFitM = m_trackFitResults.appendNew(helixParam, helixError, pTypeM, pValue, -1, -1, 0);
721
722 trackFitMIndex = trackFitM->getArrayIndex();
723
724 daughterM = Particle(trackID[1] - 1, trackFitM, pTypeM);
725 // set momentum/positions at pivot = V0 decay vertex
726 getHelixParameters(trk_fit, pTypeM.getMass(), dauPivot,
727 helixParam, error5x5,
728 momentumM, positionM, error7x7M);
729
730 for (unsigned i = 0; i < 7; i++)
731 for (unsigned j = 0; j < 7; j++)
732 errMatrixM(i, j) = error7x7M[i][j];
733
734 daughterM.updateMomentum(ROOT::Math::PxPyPzEVector(momentumM.px(), momentumM.py(), momentumM.pz(), momentumM.e()),
735 ROOT::Math::XYZVector(positionM.x(), positionM.y(), positionM.z()),
736 errMatrixM, pValue);
737 }
738 }
739
740 Track* trackP = m_tracks[trackID[0] - 1];
741 Track* trackM = m_tracks[trackID[1] - 1];
742
743 TrackFitResult* trackFitP = m_trackFitResults[trackFitPIndex];
744 TrackFitResult* trackFitM = m_trackFitResults[trackFitMIndex];
745
746 m_v0s.appendNew(std::make_pair(trackP, trackFitP), std::make_pair(trackM, trackFitM), belleV0.vx(), belleV0.vy(), belleV0.vz());
747
748 // create Ks Particle and add it to the 'K_S0:mdst' ParticleList
749 const PIDLikelihood* pidP = trackP->getRelated<PIDLikelihood>();
750 const PIDLikelihood* pidM = trackM->getRelated<PIDLikelihood>();
751 const MCParticle* mcParticleP = trackP->getRelated<MCParticle>();
752 const MCParticle* mcParticleM = trackM->getRelated<MCParticle>();
753
754 Particle* newDaugP = m_particles.appendNew(daughterP);
755 if (pidP)
756 newDaugP->addRelationTo(pidP);
757 if (mcParticleP)
758 newDaugP->addRelationTo(mcParticleP);
759 Particle* newDaugM = m_particles.appendNew(daughterM);
760 if (pidM)
761 newDaugM->addRelationTo(pidM);
762 if (mcParticleM)
763 newDaugM->addRelationTo(mcParticleM);
764
765 ROOT::Math::PxPyPzEVector v0Momentum(belleV0.px(), belleV0.py(), belleV0.pz(), belleV0.energy());
766 ROOT::Math::XYZVector v0Vertex(belleV0.vx(), belleV0.vy(), belleV0.vz());
767
768 /*
769 * Documentation of Mdst_vee2 vertex fit:
770 * /sw/belle/belle/b20090127_0910/share/tables/mdst.tdf (L96-L125)
771 */
772 auto appendVertexFitInfo = [](Belle::Mdst_vee2 & _belle_V0, Particle & _belle2_V0) {
773 // Add chisq of vertex fit. chiSquared=10^10 means the fit fails.
774 _belle2_V0.addExtraInfo("chiSquared", _belle_V0.chisq());
775 // Ndf of the vertex kinematic fit is 1
776 _belle2_V0.addExtraInfo("ndf", 1);
777 // Add p-value to extra Info
778 double prob = TMath::Prob(_belle_V0.chisq(), 1);
779 _belle2_V0.setPValue(prob);
780 };
781
782 Particle* newV0 = nullptr;
783 if (belleV0.kind() == 1) { // K0s -> pi+ pi-
784 Particle KS(v0Momentum, 310);
785 KS.appendDaughter(newDaugP);
786 KS.appendDaughter(newDaugM);
787 KS.setVertex(v0Vertex);
788 appendVertexFitInfo(belleV0, KS);
789 newV0 = m_particles.appendNew(KS);
790 ksPList->addParticle(newV0);
791
792 // append extra info: goodKs flag
793 Belle::FindKs belleKSFinder;
794 belleKSFinder.candidates(belleV0, Belle::IpProfile::position(1));
795 newV0->addExtraInfo("goodKs", belleKSFinder.goodKs());
796
797 /*
798 std::cout << " ---- B1 Ks ---- " << std::endl;
799 std::cout << " momentum = " << std::endl;
800 v0Momentum.Print();
801 std::cout << " position = " << std::endl;
802 v0Vertex.Print();
803 std::cout << " ---- B2 Ks ---- " << std::endl;
804 std::cout << " momentum = " << std::endl;
805 newKS->get4Vector().Print();
806 std::cout << " position = " << std::endl;
807 newKS->getVertex().Print();
808 std::cout << " ---- B1 Ks.child(0) ---- " << std::endl;
809 std::cout << " momentum = " << momentumP << std::endl;
810 std::cout << " position = " << positionP << std::endl;
811 std::cout << " error7x7 = " << error7x7P << std::endl;
812 std::cout << " ---- B2 Ks.child(0) ---- " << std::endl;
813 std::cout << " momentum = " << std::endl;
814 newKS->getDaughter(0)->get4Vector().Print();
815 std::cout << " position = " << std::endl;
816 newKS->getDaughter(0)->getVertex().Print();
817 std::cout << " error7x7 = " << std::endl;
818 newKS->getDaughter(0)->getMomentumVertexErrorMatrix().Print();
819 std::cout << " ---- B1 Ks.child(1) ---- " << std::endl;
820 std::cout << " momentum = " << momentumM << std::endl;
821 std::cout << " position = " << positionM << std::endl;
822 std::cout << " error7x7 = " << error7x7M << std::endl;
823 std::cout << " ---- B2 Ks.child(1) ---- " << std::endl;
824 std::cout << " momentum = " << std::endl;
825 newKS->getDaughter(1)->get4Vector().Print();
826 std::cout << " position = " << std::endl;
827 newKS->getDaughter(1)->getVertex().Print();
828 std::cout << " error7x7 = " << std::endl;
829 newKS->getDaughter(1)->getMomentumVertexErrorMatrix().Print();
830 */
831 } else if (belleV0.kind() == 2) { // Lambda -> p+ pi-
832 Particle Lambda0(v0Momentum, 3122);
833 Lambda0.appendDaughter(newDaugP);
834 Lambda0.appendDaughter(newDaugM);
835 Lambda0.setVertex(v0Vertex);
836 appendVertexFitInfo(belleV0, Lambda0);
837 newV0 = m_particles.appendNew(Lambda0);
838 lambda0PList->addParticle(newV0);
839
840 // GoodLambda flag as extra info
841 Belle::FindLambda lambdaFinder;
842 lambdaFinder.candidates(belleV0, Belle::IpProfile::position(1));
843 newV0->addExtraInfo("goodLambda", lambdaFinder.goodLambda());
844 } else if (belleV0.kind() == 3) { // anti-Lambda -> pi+ anti-p
845 Particle antiLambda0(v0Momentum, -3122);
846 antiLambda0.appendDaughter(newDaugM);
847 antiLambda0.appendDaughter(newDaugP);
848 antiLambda0.setVertex(v0Vertex);
849 appendVertexFitInfo(belleV0, antiLambda0);
850 newV0 = m_particles.appendNew(antiLambda0);
851 antiLambda0PList->addParticle(newV0);
852
853 // GoodLambda flag as extra info
854 Belle::FindLambda lambdaFinder;
855 lambdaFinder.candidates(belleV0, Belle::IpProfile::position(1));
856 newV0->addExtraInfo("goodLambda", lambdaFinder.goodLambda());
857 } else if (belleV0.kind() == 4) { // gamma -> e+ e-
858 Particle gamma(v0Momentum, 22);
859 gamma.appendDaughter(newDaugP);
860 gamma.appendDaughter(newDaugM);
861 gamma.setVertex(v0Vertex);
862 appendVertexFitInfo(belleV0, gamma);
863 newV0 = m_particles.appendNew(gamma);
864 convGammaPList->addParticle(newV0);
865 }
866 // append extra info: nisKsFinder quality indicators
867 if (m_nisEnable) {
868 if (belleV0.kind() > 0 and belleV0.kind() <= 3) { // K_S0, Lambda, anti-Lambda
869 Belle::nisKsFinder ksnb;
870 double protIDP = atcPID(pidP, 2, 4);
871 double protIDM = atcPID(pidM, 2, 4);
872 ksnb.candidates(belleV0, Belle::IpProfile::position(1), momentumP, protIDP, protIDM);
873 // K_S0 and Lambda (inverse cut on ksnbNoLam for Lambda selection).
874 newV0->addExtraInfo("ksnbVLike", ksnb.nb_vlike());
875 newV0->addExtraInfo("ksnbNoLam", ksnb.nb_nolam());
876 // K_S0 only
877 if (belleV0.kind() == 1)
878 newV0->addExtraInfo("ksnbStandard", ksnb.standard());
879 }
880 }
881 }
882}
Belle2::B2BIIConvertMdstModule::belleVeeDaughterHelix
void belleVeeDaughterHelix(const Belle::Mdst_vee2 &vee, const int charge, std::vector< float > &helixParam, std::vector< float > &helixError)
obtains the helix parameters of the vee daughters
Definition: B2BIIConvertMdstModule.cc:2211
Belle2::B2BIIConvertMdstModule::atcPID
double atcPID(const PIDLikelihood *pid, int sigHyp, int bkgHyp)
calculates atc_pid(3,1,5,sigHyp,bkgHyp).prob() from converted PIDLikelihood
Definition: B2BIIConvertMdstModule.cc:2303
Belle2::B2BIIConvertMdstModule::belleVeeDaughterToCartesian
void belleVeeDaughterToCartesian(const Belle::Mdst_vee2 &vee, const int charge, const Const::ParticleType &pType, CLHEP::HepLorentzVector &momentum, HepPoint3D &position, CLHEP::HepSymMatrix &error)
Fills 4-momentum, position and 7x7 error matrix from Belle Vee daughter.
Definition: B2BIIConvertMdstModule.cc:2172
Belle2::B2BIIConvertMdstModule::m_v0s
StoreArray< V0 > m_v0s
V0-particles.
Definition: B2BIIConvertMdstModule.h:396
Belle2::B2BIIConvertMdstModule::createTrackFitResult
TrackFitResult createTrackFitResult(const CLHEP::HepLorentzVector &momentum, const HepPoint3D &position, const CLHEP::HepSymMatrix &error, const short int charge, const Const::ParticleType &pType, const float pValue, const uint64_t hitPatternCDCInitializer, const uint32_t hitPatternVXDInitializer, const uint16_t ndf)
Creates TrackFitResult and fills it.
Definition: B2BIIConvertMdstModule.cc:2272
Belle2::Const::proton
static const ChargedStable proton
proton particle
Definition: Const.h:663
Belle2::Const::electron
static const ChargedStable electron
electron particle
Definition: Const.h:659
Belle2::PIDLikelihood
Class to collect log likelihoods from TOP, ARICH, dEdx, ECL and KLM aimed for output to mdst includes...
Definition: PIDLikelihood.h:29
Belle2::Particle::updateMomentum
void updateMomentum(const ROOT::Math::PxPyPzEVector &p4, const ROOT::Math::XYZVector &vertex, const TMatrixFSym &errMatrix, double pValue)
Sets Lorentz vector, position, 7x7 error matrix and p-value.
Definition: Particle.h:397
Belle2::Track
Class that bundles various TrackFitResults.
Definition: Track.h:25

◆ convertPIDData()

void convertPIDData ( const Belle::Mdst_charged &  belleTrack,
const Track *  track 
)
private

Get PID information for belleTrack and add it to PIDLikelihood (with relation from track).

Definition at line 1414 of file B2BIIConvertMdstModule.cc.

1415{
1416 PIDLikelihood* pid = m_pidLikelihoods.appendNew();
1417 track->addRelationTo(pid);
1418
1419 //convert data handled by atc_pid: dE/dx (-> CDC), TOF (-> TOP), ACC ( -> ARICH)
1420 //this should result in the same likelihoods used when creating atc_pid(3, 1, 5, ..., ...)
1421 //and calling prob(const Mdst_charged & chg).
1422
1423 double likelihoods[c_nHyp];
1424 double accL[c_nHyp];
1425 double tofL[c_nHyp];
1426 double cdcL[c_nHyp];
1427 for (int i = 0; i < c_nHyp; i++) {
1428 accL[i] = tofL[i] = cdcL[i] = 1.0; // cppcheck-suppress unreadVariable
1429 }
1430#ifdef HAVE_KID_ACC
1431 //accq0 = 3, as implemented in acc_prob3()
1432 const auto& acc = belleTrack.acc();
1433 if (acc and acc.quality() == 0) {
1434 for (int i = 0; i < c_nHyp; i++)
1435 accL[i] = likelihoods[i] = acc_pid(belleTrack, i); // cppcheck-suppress unreadVariable
1436 setLikelihoods(pid, Const::ARICH, likelihoods, true);
1437 }
1438#endif
1439
1440 //tofq0 = 1, as implemented in tof_prob1()
1441 //uses p1 / (p1 + p2) to create probability, so this should map directly to likelihoods
1442 const Belle::Mdst_tof& tof = belleTrack.tof();
1443 if (tof and tof.quality() == 0) {
1444 for (int i = 0; i < c_nHyp; i++)
1445 tofL[i] = likelihoods[i] = tof.pid(i); // cppcheck-suppress unreadVariable
1446 setLikelihoods(pid, Const::TOP, likelihoods, true);
1447 }
1448
1449 // cdcq0 = 5, as implemented in cdc_prob0() (which is used for all values of cdcq0!)
1450 //uses p1 / (p1 + p2) to create probability, so this should map directly to likelihoods
1451 // eID actually uses cdc_pid (cdc_prob5)
1452 const Belle::Mdst_trk& trk = belleTrack.trk();
1453 if (trk.dEdx() > 0) {
1454 for (int i = 0; i < c_nHyp; i++) {
1455 likelihoods[i] = trk.pid(i);
1456 cdcL[i] = cdc_pid(belleTrack, i); // cppcheck-suppress unreadVariable
1457 }
1458 setLikelihoods(pid, Const::CDC, likelihoods, true);
1459 }
1460
1461
1462 // eid
1463 // eid is combination of atc_pid and ecl related information
1464 // since atc_pid part is already converted above only the ECL part
1465 // is converted
1466 // ECL pdfs are available only for electrons and hadrons (assumed to be pions)
1467 // likelihoods for others are set to 0
1468
1469#ifdef HAVE_EID
1470 Belle::eid electronID(belleTrack);
1471 float eclID_e_pdf = electronID.pdf_e_ecl();
1472 float eclID_h_pdf = electronID.pdf_h_ecl();
1473 float atcID_e_pdf = electronID.atc_pid_pdf(true, accL, tofL, cdcL);
1474 float atcID_h_pdf = electronID.atc_pid_pdf(false, accL, tofL, cdcL);
1475
1476 // eID
1477 float eclProb = eclID_e_pdf / (eclID_e_pdf + eclID_h_pdf);
1478 float atcProb = atcID_e_pdf / (atcID_e_pdf + atcID_h_pdf);
1479
1480 if (atcProb > 0.999999) atcProb = 0.999999;
1481 // combine the two probabilities.
1482 double eidCombinedSig = eclProb * atcProb;
1483 double eidCombinedBkg = (1. - eclProb) * (1. - atcProb);
1484
1485 likelihoods[0] = eidCombinedSig;
1486 likelihoods[1] = 0; // no muons
1487 likelihoods[2] = eidCombinedBkg;
1488 likelihoods[3] = 0; // no kaons
1489 likelihoods[4] = 0; // no protons
1490
1491 setLikelihoods(pid, Const::ECL, likelihoods, true);
1492
1493 //Hep3Vector mom(belleTrack.px(), belleTrack.py(), belleTrack.pz());
1494 //B2INFO(" p = " << mom.mag() << " le_ecl = " << electronID.le_ecl());
1495#endif
1496
1497 //muid
1498 //Note that though it says "_likelihood()" on the label, those are
1499 //actually likelihood ratios of the type L(hyp) / (L(mu) + L(pi) + L(K)),
1500 //which are set in the FixMdst module.
1501 int muid_trackid = belleTrack.muid_ID();
1502 if (muid_trackid) {
1503 //Using approach 2. from http://belle.kek.jp/secured/muid/usage_muid.html since
1504 //it's much simpler than what Muid_mdst does.
1505 Belle::Mdst_klm_mu_ex_Manager& ex_mgr = Belle::Mdst_klm_mu_ex_Manager::get_manager();
1506 Belle::Mdst_klm_mu_ex& ex = ex_mgr(Belle::Panther_ID(muid_trackid));
1507
1508 //filter out tracks with insufficient #hits (equal to cut on Muid_mdst::Chi_2())
1509 if (ex.Chi_2() > 0) {
1510 likelihoods[0] = 0; //no electrons
1511 likelihoods[1] = ex.Muon_likelihood();
1512 likelihoods[2] = ex.Pion_likelihood();
1513 likelihoods[3] = ex.Kaon_likelihood();
1514 likelihoods[4] = 0; //no protons
1515 //Miss_likelihood should only be != 0 for tracks that do not pass the Chi_2 cut.
1516
1517 // in some cases the ex.XYZ_likelihood() < 0; Set it to 0 in these cases.
1518 for (int i = 0; i < 5; i++)
1519 if (likelihoods[i] < 0)
1520 likelihoods[i] = 0;
1521
1522 //note: discard_allzero = false since all likelihoods = 0 usually means that Junk_likelihood is 1
1523 // PIDLikelihood::getProbability(hyp) will correctly return 0 then.
1524 setLikelihoods(pid, Const::KLM, likelihoods);
1525
1526 /*
1527 const double tolerance = 1e-7;
1528 if (fabs(pid->getProbability(Const::muon, nullptr, Const::KLM) - ex.Muon_likelihood()) > tolerance ||
1529 fabs(pid->getProbability(Const::pion, nullptr, Const::KLM) - ex.Pion_likelihood()) > tolerance ||
1530 fabs(pid->getProbability(Const::kaon, nullptr, Const::KLM) - ex.Kaon_likelihood()) > tolerance) {
1531
1532 B2INFO("muons: " << pid->getProbability(Const::muon, nullptr, Const::KLM) << " " << ex.Muon_likelihood());
1533 B2INFO("pion: " << pid->getProbability(Const::pion, nullptr, Const::KLM) << " " << ex.Pion_likelihood());
1534 B2INFO("kaon: " << pid->getProbability(Const::kaon, nullptr, Const::KLM) << " " << ex.Kaon_likelihood());
1535 B2INFO("miss/junk: " << ex.Miss_likelihood() << " " << ex.Junk_likelihood());
1536 }
1537 */
1538 }
1539 }
1540}
Belle2::B2BIIConvertMdstModule::setLikelihoods
void setLikelihoods(PIDLikelihood *pid, Const::EDetector det, double likelihoods[c_nHyp], bool discard_allzero=false)
Add given Belle likelihoods (not log-likelihoods, in Belle hypothesis order) for given detector to pi...
Definition: B2BIIConvertMdstModule.cc:1396
Belle2::B2BIIConvertMdstModule::cdc_pid
double cdc_pid(const Belle::Mdst_charged &chg, int idp)
Returns CDC likelihood for given hypothesis idp.
Belle2::B2BIIConvertMdstModule::acc_pid
static double acc_pid(const Belle::Mdst_charged &chg, int idp)
Returns ACC likelihood for given hypothesis idp.
Belle2::B2BIIConvertMdstModule::m_pidLikelihoods
StoreArray< PIDLikelihood > m_pidLikelihoods
output PIDLikelihood array.
Definition: B2BIIConvertMdstModule.h:405

◆ convertRecTrgTable()

void convertRecTrgTable ( )
private

Reads and converts m_final from rectrg_summary3.

Definition at line 1311 of file B2BIIConvertMdstModule.cc.

1312{
1313
1314 // Create StoreObj if it is not valid
1315 if (not m_evtInfo.isValid()) {
1316 m_evtInfo.create();
1317 }
1318 // Load event info
1319 StoreObjPtr<EventMetaData> event;
1320
1321 if (event->getExperiment() <= 27) { // Check if it's SVD 1
1322 // Pull rectrg_summary from namager
1323 Belle::Rectrg_summary_Manager& RecTrgSummaryMgr = Belle::Rectrg_summary_Manager::get_manager();
1324 std::vector<Belle::Rectrg_summary>::iterator eflagIterator = RecTrgSummaryMgr.begin();
1325 std::string name_summary = "rectrg_summary_m_final";
1326
1327 // Converting m_final(2) from rectrg_summary
1328 for (int index = 0; index < 2; ++index) {
1329 std::string iVar = name_summary + std::to_string(index);
1330 m_evtInfo->addExtraInfo(iVar, (*eflagIterator).final(index));
1331 B2DEBUG(99, "m_final(" << index << ") = " << m_evtInfo->getExtraInfo(iVar));
1332 }
1333 } else { // For SVD2
1334 // Pull rectrg_summary3 from manager
1335 Belle::Rectrg_summary3_Manager& RecTrgSummary3Mgr = Belle::Rectrg_summary3_Manager::get_manager();
1336
1337 std::string name_summary3 = "rectrg_summary3_m_final";
1338 // Only one entry in each event
1339 std::vector<Belle::Rectrg_summary3>::iterator eflagIterator3 = RecTrgSummary3Mgr.begin();
1340
1341 // Converting m_final(3)
1342 for (int index = 0; index < 3; ++index) {
1343 std::string iVar = name_summary3 + std::to_string(index);
1344 m_evtInfo->addExtraInfo(iVar, (*eflagIterator3).final(index));
1345 B2DEBUG(99, "m_final(" << index << ") = " << m_evtInfo->getExtraInfo(iVar));
1346 }
1347 }
1348
1349}

◆ copyNbarFromGamma()

void copyNbarFromGamma ( )
private

Copies Particles in 'gamma:mdst' with energy > 0.5 GeV to be anti-n0:mdst.

Definition at line 1107 of file B2BIIConvertMdstModule.cc.

1108{
1109 StoreObjPtr<ParticleList> plist("anti-n0:mdst");
1110 plist.create();
1111 plist->initialize(Const::antiNeutron.getPDGCode(), "anti-n0:mdst");
1112
1113 B2DEBUG(99, "Getting gamma:mdst in copyNbarFromGamma");
1114 StoreObjPtr<ParticleList> plist_gamma("gamma:mdst");
1115 for (const Particle& gamma : *plist_gamma) {
1116 auto* eclCluster = gamma.getECLCluster();
1117 // Pre-select energetic gamma
1118 if (eclCluster->getEnergy(ECLCluster::EHypothesisBit::c_nPhotons) <= 0.5) continue;
1119 B2DEBUG(99, "Copying anti-n0:mdst from gamma:mdst");
1120 Particle* nbar = m_particles.appendNew(eclCluster, Const::antiNeutron);
1121 plist->addParticle(nbar);
1122
1123 if (m_realData)
1124 continue;
1125
1126 // Relation to MCParticle
1127 MCParticle* matchedMCParticle = eclCluster->getRelated<MCParticle>();
1128 if (matchedMCParticle)
1129 nbar->addRelationTo(matchedMCParticle);
1130 }
1131}
Belle2::Const::antiNeutron
static const ParticleType antiNeutron
Anti-neutron particle.
Definition: Const.h:676

◆ createTrackFitResult()

TrackFitResult createTrackFitResult ( const CLHEP::HepLorentzVector &  momentum,
const HepPoint3D &  position,
const CLHEP::HepSymMatrix &  error,
const short int  charge,
const Const::ParticleType &  pType,
const float  pValue,
const uint64_t  hitPatternCDCInitializer,
const uint32_t  hitPatternVXDInitializer,
const uint16_t  ndf 
)
private

Creates TrackFitResult and fills it.

Definition at line 2272 of file B2BIIConvertMdstModule.cc.

2281{
2282 ROOT::Math::XYZVector pos(position.x(), position.y(), position.z());
2283 ROOT::Math::XYZVector mom(momentum.px(), momentum.py(), momentum.pz());
2284
2285 TMatrixDSym errMatrix(6);
2286 for (unsigned i = 0; i < 7; i++) {
2287 if (i == 3)
2288 continue;
2289 for (unsigned j = 0; j < 7; j++) {
2290 if (j == 3)
2291 continue;
2292
2293 if (i == j)
2294 errMatrix(ERRMCONV[i], ERRMCONV[i]) = error[i][i];
2295 else
2296 errMatrix(ERRMCONV[i], ERRMCONV[j]) = errMatrix(ERRMCONV[j], ERRMCONV[i]) = error[i][j];
2297 }
2298 }
2299
2300 return TrackFitResult(pos, mom, errMatrix, charge, pType, pValue, BFIELD, hitPatternCDCInitializer, hitPatternVXDInitializer, ndf);
2301}

◆ def_beginRun()

virtual void def_beginRun ( )
inlineprotectedvirtualinherited

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

Reimplemented in PyModule.

Definition at line 425 of file Module.h.

425{ beginRun(); }
Belle2::Module::beginRun
virtual void beginRun()
Called when entering a new run.
Definition: Module.h:146

◆ def_endRun()

virtual void def_endRun ( )
inlineprotectedvirtualinherited

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

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

Reimplemented in PyModule.

Definition at line 438 of file Module.h.

438{ endRun(); }
Belle2::Module::endRun
virtual void endRun()
This method is called if the current run ends.
Definition: Module.h:165

◆ def_event()

virtual void def_event ( )
inlineprotectedvirtualinherited

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

Reimplemented in PyModule.

Definition at line 431 of file Module.h.

431{ event(); }
Belle2::Module::event
virtual void event()
This method is the core of the module.
Definition: Module.h:156

◆ def_initialize()

virtual void def_initialize ( )
inlineprotectedvirtualinherited

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

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

Reimplemented in PyModule.

Definition at line 419 of file Module.h.

419{ initialize(); }
Belle2::Module::initialize
virtual void initialize()
Initialize the Module.
Definition: Module.h:109

◆ def_terminate()

virtual void def_terminate ( )
inlineprotectedvirtualinherited

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

Reimplemented in PyModule.

Definition at line 444 of file Module.h.

444{ terminate(); }
Belle2::Module::terminate
virtual void terminate()
This method is called at the end of the event processing.
Definition: Module.h:175

◆ endRun()

void endRun ( void  )
overridevirtual

Called when the current run is finished.

Reimplemented from Module.

Definition at line 2347 of file B2BIIConvertMdstModule.cc.

2348{
2349 B2DEBUG(99, "B2BIIConvertMdst: endRun done.");
2350}

◆ evalCondition()

bool evalCondition ( ) const
inherited

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

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

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

Definition at line 96 of file Module.cc.

97{
98 if (m_conditions.empty()) return false;
99
100 //okay, a condition was set for this Module...
101 if (!m_hasReturnValue) {
102 B2FATAL("A condition was set for '" << getName() << "', but the module did not set a return value!");
103 }
104
105 for (const auto& condition : m_conditions) {
106 if (condition.evaluate(m_returnValue)) {
107 return true;
108 }
109 }
110 return false;
111}
Belle2::Module::m_returnValue
int m_returnValue
The return value.
Definition: Module.h:518
Belle2::Module::m_hasReturnValue
bool m_hasReturnValue
True, if the return value is set.
Definition: Module.h:517

◆ event()

void event ( void  )
overridevirtual

Called for each event.

Reimplemented from Module.

Definition at line 295 of file B2BIIConvertMdstModule.cc.

296{
297 if (m_convertBeamParameters) {
298 // Are we running on MC or DATA?
299 Belle::Belle_event_Manager& evman = Belle::Belle_event_Manager::get_manager();
300 Belle::Belle_event& evt = evman[0];
301
302 if (evt.ExpMC() == 2)
303 m_realData = false; // <- this is MC sample
304 else
305 m_realData = true; // <- this is real data sample
306
307 // 0. Convert IPProfile to BeamSpot
308 convertIPProfile();
309
310 // Make sure beam parameters are correct: if they are not found in the
311 // database or different from the ones in the database we need to override them
312 if (!m_beamSpotDB || !(m_beamSpot == *m_beamSpotDB) ||
313 !m_collisionBoostVectorDB || !m_collisionInvMDB || !m_collisionAxisCMSDB) {
314 if ((!m_beamSpotDB || !m_collisionBoostVectorDB || !m_collisionInvMDB || !m_collisionAxisCMSDB) && !m_realData) {
315 B2INFO("No database entry for this run yet, create one");
316 StoreObjPtr<EventMetaData> event;
317 IntervalOfValidity iov(event->getExperiment(), event->getRun(), event->getExperiment(), event->getRun());
318 Database::Instance().storeData("CollisionBoostVector", &m_collisionBoostVector, iov);
319 Database::Instance().storeData("CollisionInvariantMass", &m_collisionInvM, iov);
320 Database::Instance().storeData("CollisionAxisCMS", &m_collisionAxisCMS, iov);
321 Database::Instance().storeData("BeamSpot", &m_beamSpot, iov);
322 }
323 if (m_realData) {
324 B2ERROR("BeamParameters from condition database are different from converted "
325 "ones, overriding database. Did you make sure the globaltag B2BII is used?");
326 } else {
327 B2INFO("BeamSpot, BoostVector, and InvariantMass from condition database are different from converted "
328 "ones, overriding database");
329 }
330 if (ProcHandler::parallelProcessingUsed()) {
331 B2FATAL("Cannot reliably override the Database content in parallel processing "
332 "mode, please run the conversion in single processing mode");
333 }
334 DBStore::Instance().addConstantOverride("CollisionBoostVector", new CollisionBoostVector(m_collisionBoostVector), true);
335 DBStore::Instance().addConstantOverride("CollisionInvariantMass", new CollisionInvariantMass(m_collisionInvM), true);
336 DBStore::Instance().addConstantOverride("CollisionAxisCMS", new CollisionAxisCMS(m_collisionAxisCMS), true);
337 DBStore::Instance().addConstantOverride("BeamSpot", new BeamSpot(m_beamSpot), true);
338 }
339 }
340
341 // 1. Convert MC information
342 convertGenHepEvtTable();
343
344 // 2. Convert Tracking information
345 convertMdstChargedTable();
346
347 // 3. Convert ECL information
348 convertMdstECLTable();
349
350 // 4. Convert KLM information
351 convertMdstKLMTable();
352
353 // 5. Set Track -> ECLCluster relations
354 setTracksToECLClustersRelations();
355
356 // 6. Set KLMCluster -> Track, ECLCluster relations
357 setKLMClustersRelations();
358
359 // 7. Convert Gamma information
360 convertMdstGammaTable();
361
362 // 8. Convert Pi0 information
363 convertMdstPi0Table();
364
365 // 9. Convert V0s
366 convertMdstVee2Table();
367
368 // 10. Convert KLong information
369 convertMdstKLongTable();
370
371 // 11. Convert Evtcls panther table information
372 if (m_convertEvtcls) convertEvtclsTable();
373
374 // 12. Convert trigger information from rectrg_summary3
375 if (m_convertRecTrg) convertRecTrgTable();
376
377 // 13. Copy nbar from Gamma with the cut E > 0.5 GeV
378 if (m_convertNbar) copyNbarFromGamma();
379
380}
Belle2::B2BIIConvertMdstModule::convertMdstGammaTable
void convertMdstGammaTable()
Reads all entries of Mdst_Gamma Panther table, creates a particle list 'gamma:mdst' and adds them to ...
Definition: B2BIIConvertMdstModule.cc:1060
Belle2::B2BIIConvertMdstModule::copyNbarFromGamma
void copyNbarFromGamma()
Copies Particles in 'gamma:mdst' with energy > 0.5 GeV to be anti-n0:mdst.
Definition: B2BIIConvertMdstModule.cc:1107
Belle2::B2BIIConvertMdstModule::convertMdstChargedTable
void convertMdstChargedTable()
Reads and converts all entries of Mdst_charged (Mdst_trk and Mdst_trk_fit) Panther table to Track (Tr...
Definition: B2BIIConvertMdstModule.cc:466
Belle2::B2BIIConvertMdstModule::setTracksToECLClustersRelations
void setTracksToECLClustersRelations()
Sets Track -> ECLCluster relations.
Definition: B2BIIConvertMdstModule.cc:2008
Belle2::B2BIIConvertMdstModule::m_collisionBoostVectorDB
OptionalDBObjPtr< CollisionBoostVector > m_collisionBoostVectorDB
CollisionBoostVector for boost vector.
Definition: B2BIIConvertMdstModule.h:415
Belle2::B2BIIConvertMdstModule::m_collisionAxisCMSDB
OptionalDBObjPtr< CollisionAxisCMS > m_collisionAxisCMSDB
CollisionAxisCMS.
Definition: B2BIIConvertMdstModule.h:423
Belle2::B2BIIConvertMdstModule::convertGenHepEvtTable
void convertGenHepEvtTable()
Reads and converts all entries of Gen_hepevt Panther table to MCParticle dataobjects and adds them to...
Definition: B2BIIConvertMdstModule.cc:884
Belle2::B2BIIConvertMdstModule::convertRecTrgTable
void convertRecTrgTable()
Reads and converts m_final from rectrg_summary3.
Definition: B2BIIConvertMdstModule.cc:1311
Belle2::B2BIIConvertMdstModule::setKLMClustersRelations
void setKLMClustersRelations()
Sets KLMCluster -> Track and ECLCluster relations.
Definition: B2BIIConvertMdstModule.cc:2051
Belle2::B2BIIConvertMdstModule::m_collisionInvMDB
OptionalDBObjPtr< CollisionInvariantMass > m_collisionInvMDB
CollisionInvariantMass for Invariant Mass of Beam.
Definition: B2BIIConvertMdstModule.h:419
Belle2::B2BIIConvertMdstModule::convertMdstPi0Table
void convertMdstPi0Table()
Reads all entries of Mdst_Pi0 Panther table, creates a particle list 'pi0:mdst' and adds them to Stor...
Definition: B2BIIConvertMdstModule.cc:1133
Belle2::B2BIIConvertMdstModule::convertMdstKLongTable
void convertMdstKLongTable()
Reads all entries of Mdst_Klong Panther table, creates a particle list 'K_L0:mdst' and adds them to S...
Definition: B2BIIConvertMdstModule.cc:1181
Belle2::B2BIIConvertMdstModule::m_beamSpotDB
OptionalDBObjPtr< BeamSpot > m_beamSpotDB
BeamSpot for IP.
Definition: B2BIIConvertMdstModule.h:411
Belle2::B2BIIConvertMdstModule::convertMdstECLTable
void convertMdstECLTable()
Reads and converts all entries of Mdst_ecl(_aux) Panther table to ECLCluster dataobjects and adds the...
Definition: B2BIIConvertMdstModule.cc:973
Belle2::B2BIIConvertMdstModule::convertMdstKLMTable
void convertMdstKLMTable()
Reads and converts all entries of Mdst_klm_cluster Panther table to KLMCluster dataobjects and adds t...
Definition: B2BIIConvertMdstModule.cc:1033
Belle2::B2BIIConvertMdstModule::convertMdstVee2Table
void convertMdstVee2Table()
Reads and converts all entries of Mdst_vee2 Panther table to V0 dataobjects and adds them to StoreArr...
Definition: B2BIIConvertMdstModule.cc:526
Belle2::B2BIIConvertMdstModule::convertEvtclsTable
void convertEvtclsTable()
Reads and converts all entries of evtcls Panther table.
Definition: B2BIIConvertMdstModule.cc:1266
Belle2::BeamSpot
This class contains the beam spot position and size modeled as a gaussian distribution in space.
Definition: BeamSpot.h:22
Belle2::CollisionAxisCMS
This class contains the measured values of the orientation of the collision axis in the CM system obt...
Definition: CollisionAxisCMS.h:27
Belle2::CollisionBoostVector
This class contains the measured average boost vector vec(beta) = (beta_x, beta_y,...
Definition: CollisionBoostVector.h:23
Belle2::CollisionInvariantMass
This class contains the measured average center-of-mass energy, which is equal to the invariant mass ...
Definition: CollisionInvariantMass.h:21
Belle2::IntervalOfValidity
A class that describes the interval of experiments/runs for which an object in the database is valid.
Definition: IntervalOfValidity.h:25
Belle2::ProcHandler::parallelProcessingUsed
static bool parallelProcessingUsed()
Returns true if multiple processes have been spawned, false in single-core mode.
Definition: ProcHandler.cc:226
Belle2::Database::Instance
static Database & Instance()
Instance of a singleton Database.
Definition: Database.cc:41
Belle2::DBStore::Instance
static DBStore & Instance()
Instance of a singleton DBStore.
Definition: DBStore.cc:26
Belle2::Database::storeData
bool storeData(const std::string &name, TObject *object, const IntervalOfValidity &iov)
Store an object in the database.
Definition: Database.cc:140
Belle2::DBStore::addConstantOverride
void addConstantOverride(const std::string &name, TObject *obj, bool oneRun=false)
Add constant override payload.
Definition: DBStore.cc:202

◆ exposePythonAPI()

void exposePythonAPI ( )
staticinherited

Exposes methods of the Module class to Python.

Definition at line 325 of file Module.cc.

326{
327 // to avoid confusion between std::arg and boost::python::arg we want a shorthand namespace as well
328 namespace bp = boost::python;
329
330 docstring_options options(true, true, false); //userdef, py sigs, c++ sigs
331
332 void (Module::*setReturnValueInt)(int) = &Module::setReturnValue;
333
334 enum_<Module::EAfterConditionPath>("AfterConditionPath",
335 R"(Determines execution behaviour after a conditional path has been executed:
336
337.. attribute:: END
338
339 End processing of this path after the conditional path. (this is the default for if_value() etc.)
340
341.. attribute:: CONTINUE
342
343 After the conditional path, resume execution after this module.)")
344 .value("END", Module::EAfterConditionPath::c_End)
345 .value("CONTINUE", Module::EAfterConditionPath::c_Continue)
346 ;
347
348 /* Do not change the names of >, <, ... we use them to serialize conditional paths */
349 enum_<Belle2::ModuleCondition::EConditionOperators>("ConditionOperator")
350 .value(">", Belle2::ModuleCondition::EConditionOperators::c_GT)
351 .value("<", Belle2::ModuleCondition::EConditionOperators::c_ST)
352 .value(">=", Belle2::ModuleCondition::EConditionOperators::c_GE)
353 .value("<=", Belle2::ModuleCondition::EConditionOperators::c_SE)
354 .value("==", Belle2::ModuleCondition::EConditionOperators::c_EQ)
355 .value("!=", Belle2::ModuleCondition::EConditionOperators::c_NE)
356 ;
357
358 enum_<Module::EModulePropFlags>("ModulePropFlags",
359 R"(Flags to indicate certain low-level features of modules, see :func:`Module.set_property_flags()`, :func:`Module.has_properties()`. Most useful flags are:
360
361.. attribute:: PARALLELPROCESSINGCERTIFIED
362
363 This module can be run in parallel processing mode safely (All I/O must be done through the data store, in particular, the module must not write any files.)
364
365.. attribute:: HISTOGRAMMANAGER
366
367 This module is used to manage histograms accumulated by other modules
368
369.. attribute:: TERMINATEINALLPROCESSES
370
371 When using parallel processing, call this module's terminate() function in all processes. This will also ensure that there is exactly one process (single-core if no parallel modules found) or at least one input, one main and one output process.
372)")
373 .value("INPUT", Module::EModulePropFlags::c_Input)
374 .value("OUTPUT", Module::EModulePropFlags::c_Output)
375 .value("PARALLELPROCESSINGCERTIFIED", Module::EModulePropFlags::c_ParallelProcessingCertified)
376 .value("HISTOGRAMMANAGER", Module::EModulePropFlags::c_HistogramManager)
377 .value("INTERNALSERIALIZER", Module::EModulePropFlags::c_InternalSerializer)
378 .value("TERMINATEINALLPROCESSES", Module::EModulePropFlags::c_TerminateInAllProcesses)
379 ;
380
381 //Python class definition
382 class_<Module, PyModule> module("Module", R"(
383Base class for Modules.
384
385A module is the smallest building block of the framework.
386A typical event processing chain consists of a Path containing
387modules. By inheriting from this base class, various types of
388modules can be created. To use a module, please refer to
389:func:`Path.add_module()`. A list of modules is available by running
390``basf2 -m`` or ``basf2 -m package``, detailed information on parameters is
391given by e.g. ``basf2 -m RootInput``.
392
393The 'Module Development' section in the manual provides detailed information
394on how to create modules, setting parameters, or using return values/conditions:
395https://xwiki.desy.de/xwiki/rest/p/f4fa4/#HModuleDevelopment
396
397)");
398 module
399 .def("__str__", &Module::getPathString)
400 .def("name", &Module::getName, return_value_policy<copy_const_reference>(),
401 "Returns the name of the module. Can be changed via :func:`set_name() <Module.set_name()>`, use :func:`type() <Module.type()>` for identifying a particular module class.")
402 .def("type", &Module::getType, return_value_policy<copy_const_reference>(),
403 "Returns the type of the module (i.e. class name minus 'Module')")
404 .def("set_name", &Module::setName, args("name"), R"(
405Set custom name, e.g. to distinguish multiple modules of the same type.
406
407>>> path.add_module('EventInfoSetter')
408>>> ro = path.add_module('RootOutput', branchNames=['EventMetaData'])
409>>> ro.set_name('RootOutput_metadata_only')
410>>> print(path)
411[EventInfoSetter -> RootOutput_metadata_only]
412
413)")
414 .def("description", &Module::getDescription, return_value_policy<copy_const_reference>(),
415 "Returns the description of this module.")
416 .def("package", &Module::getPackage, return_value_policy<copy_const_reference>(),
417 "Returns the package this module belongs to.")
418 .def("available_params", &_getParamInfoListPython,
419 "Return list of all module parameters as `ModuleParamInfo` instances")
420 .def("has_properties", &Module::hasProperties, (bp::arg("properties")),
421 R"DOCSTRING(Allows to check if the module has the given properties out of `ModulePropFlags` set.
422
423>>> if module.has_properties(ModulePropFlags.PARALLELPROCESSINGCERTIFIED):
424>>> ...
425
426Parameters:
427 properties (int): bitmask of `ModulePropFlags` to check for.
428)DOCSTRING")
429 .def("set_property_flags", &Module::setPropertyFlags, args("property_mask"),
430 "Set module properties in the form of an OR combination of `ModulePropFlags`.");
431 {
432 // python signature is too crowded, make ourselves
433 docstring_options subOptions(true, false, false); //userdef, py sigs, c++ sigs
434 module
435 .def("if_value", &Module::if_value,
436 (bp::arg("expression"), bp::arg("condition_path"), bp::arg("after_condition_path")= Module::EAfterConditionPath::c_End),
437 R"DOCSTRING(if_value(expression, condition_path, after_condition_path=AfterConditionPath.END)
438
439Sets a conditional sub path which will be executed after this
440module if the return value set in the module passes the given ``expression``.
441
442Modules can define a return value (int or bool) using ``setReturnValue()``,
443which can be used in the steering file to split the Path based on this value, for example
444
445>>> module_with_condition.if_value("<1", another_path)
446
447In case the return value of the ``module_with_condition`` for a given event is
448less than 1, the execution will be diverted into ``another_path`` for this event.
449
450You could for example set a special return value if an error occurs, and divert
451the execution into a path containing :b2:mod:`RootOutput` if it is found;
452saving only the data producing/produced by the error.
453
454After a conditional path has executed, basf2 will by default stop processing
455the path for this event. This behaviour can be changed by setting the
456``after_condition_path`` argument.
457
458Parameters:
459 expression (str): Expression to determine if the conditional path should be executed.
460 This should be one of the comparison operators ``<``, ``>``, ``<=``,
461 ``>=``, ``==``, or ``!=`` followed by a numerical value for the return value
462 condition_path (Path): path to execute in case the expression is fulfilled
463 after_condition_path (AfterConditionPath): What to do once the ``condition_path`` has been executed.
464)DOCSTRING")
465 .def("if_false", &Module::if_false,
466 (bp::arg("condition_path"), bp::arg("after_condition_path")= Module::EAfterConditionPath::c_End),
467 R"DOC(if_false(condition_path, after_condition_path=AfterConditionPath.END)
468
469Sets a conditional sub path which will be executed after this module if
470the return value of the module evaluates to False. This is equivalent to
471calling `if_value` with ``expression=\"<1\"``)DOC")
472 .def("if_true", &Module::if_true,
473 (bp::arg("condition_path"), bp::arg("after_condition_path")= Module::EAfterConditionPath::c_End),
474 R"DOC(if_true(condition_path, after_condition_path=AfterConditionPath.END)
475
476Sets a conditional sub path which will be executed after this module if
477the return value of the module evaluates to True. It is equivalent to
478calling `if_value` with ``expression=\">=1\"``)DOC");
479 }
480 module
481 .def("has_condition", &Module::hasCondition,
482 "Return true if a conditional path has been set for this module "
483 "using `if_value`, `if_true` or `if_false`")
484 .def("get_all_condition_paths", &_getAllConditionPathsPython,
485 "Return a list of all conditional paths set for this module using "
486 "`if_value`, `if_true` or `if_false`")
487 .def("get_all_conditions", &_getAllConditionsPython,
488 "Return a list of all conditional path expressions set for this module using "
489 "`if_value`, `if_true` or `if_false`")
490 .add_property("logging", make_function(&Module::getLogConfig, return_value_policy<reference_existing_object>()),
491 &Module::setLogConfig)
Belle2::ModuleCondition::c_GE
@ c_GE
Greater or equal than: ">=".
Definition: ModuleCondition.h:35
Belle2::ModuleCondition::c_SE
@ c_SE
Smaller or equal than: "<=".
Definition: ModuleCondition.h:36
Belle2::ModuleCondition::c_GT
@ c_GT
Greater than: ">"
Definition: ModuleCondition.h:33
Belle2::ModuleCondition::c_NE
@ c_NE
Not equal: "!=".
Definition: ModuleCondition.h:38
Belle2::ModuleCondition::c_EQ
@ c_EQ
Equal: "=" or "=="
Definition: ModuleCondition.h:37
Belle2::ModuleCondition::c_ST
@ c_ST
Smaller than: "<"
Definition: ModuleCondition.h:34
Belle2::Module
Base class for Modules.
Definition: Module.h:72
Belle2::Module::getLogConfig
LogConfig & getLogConfig()
Returns the log system configuration.
Definition: Module.h:224
Belle2::Module::if_value
void if_value(const std::string &expression, const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
Add a condition to the module.
Definition: Module.cc:79
Belle2::Module::setPropertyFlags
void setPropertyFlags(unsigned int propertyFlags)
Sets the flags for the module properties.
Definition: Module.cc:208
Belle2::Module::if_true
void if_true(const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
A simplified version to set the condition of the module.
Definition: Module.cc:90
Belle2::Module::setReturnValue
void setReturnValue(int value)
Sets the return value for this module as integer.
Definition: Module.cc:220
Belle2::Module::setLogConfig
void setLogConfig(const LogConfig &logConfig)
Set the log system configuration.
Definition: Module.h:229
Belle2::Module::getDescription
const std::string & getDescription() const
Returns the description of the module.
Definition: Module.h:201
Belle2::Module::if_false
void if_false(const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
A simplified version to add a condition to the module.
Definition: Module.cc:85
Belle2::Module::hasCondition
bool hasCondition() const
Returns true if at least one condition was set for the module.
Definition: Module.h:310
Belle2::Module::getPackage
const std::string & getPackage() const
Returns the package this module is in.
Definition: Module.h:196
Belle2::Module::setName
void setName(const std::string &name)
Set the name of the module.
Definition: Module.h:213
Belle2::Module::hasProperties
bool hasProperties(unsigned int propertyFlags) const
Returns true if all specified property flags are available in this module.
Definition: Module.cc:160
Belle2::Module::getPathString
std::string getPathString() const override
return the module name.
Definition: Module.cc:192

◆ getAfterConditionPath()

Module::EAfterConditionPath getAfterConditionPath ( ) const
inherited

What to do after the conditional path is finished.

(defaults to c_End if no condition is set)

Definition at line 133 of file Module.cc.

134{
135 if (m_conditions.empty()) return EAfterConditionPath::c_End;
136
137 //okay, a condition was set for this Module...
138 if (!m_hasReturnValue) {
139 B2FATAL("A condition was set for '" << getName() << "', but the module did not set a return value!");
140 }
141
142 for (const auto& condition : m_conditions) {
143 if (condition.evaluate(m_returnValue)) {
144 return condition.getAfterConditionPath();
145 }
146 }
147
148 return EAfterConditionPath::c_End;
149}

◆ getAllConditionPaths()

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

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

Definition at line 150 of file Module.cc.

151{
152 std::vector<std::shared_ptr<Path>> allConditionPaths;
153 for (const auto& condition : m_conditions) {
154 allConditionPaths.push_back(condition.getPath());
155 }
156
157 return allConditionPaths;
158}

◆ getAllConditions()

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

Return all set conditions for this module.

Definition at line 323 of file Module.h.

324 {
325 return m_conditions;
326 }

◆ getCondition()

const ModuleCondition * getCondition ( ) const
inlineinherited

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

Definition at line 313 of file Module.h.

314 {
315 if (m_conditions.empty()) {
316 return nullptr;
317 } else {
318 return &m_conditions.front();
319 }
320 }

◆ getConditionPath()

std::shared_ptr< Path > getConditionPath ( ) const
inherited

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


Definition at line 113 of file Module.cc.

114{
115 PathPtr p;
116 if (m_conditions.empty()) return p;
117
118 //okay, a condition was set for this Module...
119 if (!m_hasReturnValue) {
120 B2FATAL("A condition was set for '" << getName() << "', but the module did not set a return value!");
121 }
122
123 for (const auto& condition : m_conditions) {
124 if (condition.evaluate(m_returnValue)) {
125 return condition.getPath();
126 }
127 }
128
129 // if none of the conditions were true, return a null pointer.
130 return p;
131}
Belle2::PathPtr
std::shared_ptr< Path > PathPtr
Defines a pointer to a path object as a boost shared pointer.
Definition: Path.h:35

◆ getDescription()

const std::string & getDescription ( ) const
inlineinherited

Returns the description of the module.

Definition at line 201 of file Module.h.

201{return m_description;}
Belle2::Module::m_description
std::string m_description
The description of the module.
Definition: Module.h:510

◆ getFileNames()

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

Return a list of output filenames for this modules.

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

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

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

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

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

Reimplemented in RootInputModule, StorageRootOutputModule, and RootOutputModule.

Definition at line 133 of file Module.h.

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

◆ getHelixParameters()

int getHelixParameters ( const Belle::Mdst_trk_fit &  trk_fit,
const double  mass,
const HepPoint3D &  newPivot,
std::vector< float > &  helixParams,
CLHEP::HepSymMatrix &  error5x5,
CLHEP::HepLorentzVector &  momentum,
HepPoint3D &  position,
CLHEP::HepSymMatrix &  error7x7,
const double  dPhi = 0.0 
)
private

Fills Helix parameters (converted to Belle II version), 5x5 error matrix, 4-momentum, position and 7x7 error matrix from Belle Helix stored in Mdst_trk_fit.

Definition at line 1542 of file B2BIIConvertMdstModule.cc.

1550{
1551 const HepPoint3D pivot(trk_fit.pivot_x(),
1552 trk_fit.pivot_y(),
1553 trk_fit.pivot_z());
1554
1555 CLHEP::HepVector a(5);
1556 a[0] = trk_fit.helix(0);
1557 a[1] = trk_fit.helix(1);
1558 a[2] = trk_fit.helix(2);
1559 a[3] = trk_fit.helix(3);
1560 a[4] = trk_fit.helix(4);
1561 CLHEP::HepSymMatrix Ea(5, 0);
1562 Ea[0][0] = trk_fit.error(0);
1563 Ea[1][0] = trk_fit.error(1);
1564 Ea[1][1] = trk_fit.error(2);
1565 Ea[2][0] = trk_fit.error(3);
1566 Ea[2][1] = trk_fit.error(4);
1567 Ea[2][2] = trk_fit.error(5);
1568 Ea[3][0] = trk_fit.error(6);
1569 Ea[3][1] = trk_fit.error(7);
1570 Ea[3][2] = trk_fit.error(8);
1571 Ea[3][3] = trk_fit.error(9);
1572 Ea[4][0] = trk_fit.error(10);
1573 Ea[4][1] = trk_fit.error(11);
1574 Ea[4][2] = trk_fit.error(12);
1575 Ea[4][3] = trk_fit.error(13);
1576 Ea[4][4] = trk_fit.error(14);
1577
1578 Belle::Helix helix(pivot, a, Ea);
1579
1580 int charge = 0;
1581 if (helix.kappa() > 0)
1582 charge = 1;
1583 else
1584 charge = -1;
1585
1586 if (newPivot.x() != 0. || newPivot.y() != 0. || newPivot.z() != 0.) {
1587 helix.pivot(newPivot);
1588 momentum = helix.momentum(dPhi, mass, position, error7x7);
1589 } else {
1590 if (pivot.x() != 0. || pivot.y() != 0. || pivot.z() != 0.) {
1591 helix.pivot(HepPoint3D(0., 0., 0.));
1592 momentum = helix.momentum(dPhi, mass, position, error7x7);
1593 } else {
1594 momentum = helix.momentum(dPhi, mass, position, error7x7);
1595 }
1596 }
1597
1598 convertHelix(helix, helixParams, error5x5);
1599
1600 return charge;
1601}
Belle2::EvtPDLUtil::charge
double charge(int pdgCode)
Returns electric charge of a particle with given pdg code.
Definition: EvtPDLUtil.cc:44

◆ getLogConfig()

LogConfig & getLogConfig ( )
inlineinherited

Returns the log system configuration.

Definition at line 224 of file Module.h.

224{return m_logConfig;}

◆ getModules()

std::list< ModulePtr > getModules ( ) const
inlineoverrideprivatevirtualinherited

no submodules, return empty list

Implements PathElement.

Definition at line 505 of file Module.h.

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

◆ getName()

const std::string & getName ( ) const
inlineinherited

Returns the name of the module.

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

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

Definition at line 186 of file Module.h.

186{return m_name;}
Belle2::Module::m_name
std::string m_name
The name of the module, saved as a string (user-modifiable)
Definition: Module.h:507

◆ getPackage()

const std::string & getPackage ( ) const
inlineinherited

Returns the package this module is in.

Definition at line 196 of file Module.h.

196{return m_package;}

◆ getParamInfoListPython()

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

Returns a python list of all parameters.

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

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

Definition at line 279 of file Module.cc.

280{
281 return m_moduleParamList.getParamInfoListPython();
282}
Belle2::ModuleParamList::getParamInfoListPython
std::shared_ptr< boost::python::list > getParamInfoListPython() const
Returns a python list of all parameters.
Definition: ModuleParamList.cc:64
Belle2::Module::m_moduleParamList
ModuleParamList m_moduleParamList
List storing and managing all parameter of the module.
Definition: Module.h:515

◆ getParamList()

const ModuleParamList & getParamList ( ) const
inlineinherited

Return module param list.

Definition at line 362 of file Module.h.

362{ return m_moduleParamList; }

◆ getPathString()

std::string getPathString ( ) const
overrideprivatevirtualinherited

return the module name.

Implements PathElement.

Definition at line 192 of file Module.cc.

193{
194
195 std::string output = getName();
196
197 for (const auto& condition : m_conditions) {
198 output += condition.getString();
199 }
200
201 return output;
202}

◆ getReturnValue()

int getReturnValue ( ) const
inlineinherited

Return the return value set by this module.

This value is only meaningful if hasReturnValue() is true

Definition at line 380 of file Module.h.

380{ return m_returnValue; }

◆ getType()

const std::string & getType ( ) const
inherited

Returns the type of the module (i.e.

class name minus 'Module')

Definition at line 41 of file Module.cc.

42{
43 if (m_type.empty())
44 B2FATAL("Module type not set for " << getName());
45 return m_type;
46}
Belle2::Module::m_type
std::string m_type
The type of the module, saved as a string.
Definition: Module.h:508

◆ hasCondition()

bool hasCondition ( ) const
inlineinherited

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

Definition at line 310 of file Module.h.

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

◆ hasProperties()

bool hasProperties ( unsigned int  propertyFlags) const
inherited

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

Parameters
propertyFlagsOred EModulePropFlags which should be compared with the module flags.

Definition at line 160 of file Module.cc.

161{
162 return (propertyFlags & m_propertyFlags) == propertyFlags;
163}

◆ hasReturnValue()

bool hasReturnValue ( ) const
inlineinherited

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

Definition at line 377 of file Module.h.

377{ return m_hasReturnValue; }

◆ hasUnsetForcedParams()

bool hasUnsetForcedParams ( ) const
inherited

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

Definition at line 166 of file Module.cc.

167{
168 auto missing = m_moduleParamList.getUnsetForcedParams();
169 std::string allMissing = "";
170 for (const auto& s : missing)
171 allMissing += s + " ";
172 if (!missing.empty())
173 B2ERROR("The following required parameters of Module '" << getName() << "' were not specified: " << allMissing <<
174 "\nPlease add them to your steering file.");
175 return !missing.empty();
176}
Belle2::ModuleParamList::getUnsetForcedParams
std::vector< std::string > getUnsetForcedParams() const
Returns list of unset parameters (if they are required to have a value.
Definition: ModuleParamList.cc:54

◆ if_false()

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

A simplified version to add a condition to the module.

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

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

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

Parameters
pathShared pointer to the Path which will be executed if the return value is false.
afterConditionPathWhat to do after executing 'path'.

Definition at line 85 of file Module.cc.

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

◆ if_true()

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

A simplified version to set the condition of the module.

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

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

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

Parameters
pathShared pointer to the Path which will be executed if the return value is true.
afterConditionPathWhat to do after executing 'path'.

Definition at line 90 of file Module.cc.

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

◆ if_value()

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

Add a condition to the module.

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

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

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

Parameters
expressionThe expression of the condition.
pathShared pointer to the Path which will be executed if the condition is evaluated to true.
afterConditionPathWhat to do after executing 'path'.

Definition at line 79 of file Module.cc.

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

◆ initialize()

void initialize ( void  )
overridevirtual

Initialize the module.

Reimplemented from Module.

Definition at line 198 of file B2BIIConvertMdstModule.cc.

199{
200 // Initialize Belle II DataStore
201 initializeDataStore();
202 if (m_mcMatchingModeString == "Direct")
203 m_mcMatchingMode = c_Direct;
204 else if (m_mcMatchingModeString == "GeneratorLevel")
205 m_mcMatchingMode = c_GeneratorLevel;
206 else
207 B2FATAL("Unknown MC matching mode: " << m_mcMatchingModeString);
208 B2INFO("B2BIIConvertMdst: initialized.");
209 if (!m_nisEnable)
210 B2WARNING("nisKsFinder output has been disabled. ksnbVLike, ksnbNoLam, ksnbStandard will not be converted.");
211}
Belle2::B2BIIConvertMdstModule::initializeDataStore
void initializeDataStore()
Initializes Belle II DataStore.
Definition: B2BIIConvertMdstModule.cc:213

◆ initializeDataStore()

void initializeDataStore ( )
private

Initializes Belle II DataStore.

Definition at line 213 of file B2BIIConvertMdstModule.cc.

214{
215 B2DEBUG(99, "[B2BIIConvertMdstModule::initializeDataStore] initialization of DataStore started");
216
217 // list here all converted Belle2 objects
218 m_eclClusters.registerInDataStore();
219 m_klmClusters.registerInDataStore();
220 m_tracks.registerInDataStore();
221 m_trackFitResults.registerInDataStore();
222 m_v0s.registerInDataStore();
223 m_particles.registerInDataStore();
224
225 StoreObjPtr<ParticleExtraInfoMap> extraInfoMap;
226 extraInfoMap.registerInDataStore();
227
228 if (m_convertEvtcls || m_convertRecTrg) m_evtInfo.isRequired();
229
230 if (m_convertTrkExtra) m_belleTrkExtra.registerInDataStore();
231
232 StoreObjPtr<ParticleList> gammaParticleList("gamma:mdst");
233 gammaParticleList.registerInDataStore();
234 StoreObjPtr<ParticleList> nbarParticleList("anti-n0:mdst");
235 nbarParticleList.registerInDataStore();
236 StoreObjPtr<ParticleList> pi0ParticleList("pi0:mdst");
237 pi0ParticleList.registerInDataStore();
238 StoreObjPtr<ParticleList> kShortParticleList("K_S0:mdst");
239 kShortParticleList.registerInDataStore();
240 StoreObjPtr<ParticleList> kLongParticleList("K_L0:mdst");
241 kLongParticleList.registerInDataStore();
242 StoreObjPtr<ParticleList> lambdaParticleList("Lambda0:mdst");
243 lambdaParticleList.registerInDataStore();
244 StoreObjPtr<ParticleList> antiLambdaParticleList("anti-Lambda0:mdst");
245 antiLambdaParticleList.registerInDataStore();
246 StoreObjPtr<ParticleList> gammaConversionsParticleList("gamma:v0mdst");
247 gammaConversionsParticleList.registerInDataStore();
248
249 m_pidLikelihoods.registerInDataStore();
250
251 // needs to be registered, even if running over data, since this information is available only at the begin_run function
252 m_mcParticles.registerInDataStore();
253
254 //list here all Relations between Belle2 objects
255 m_tracks.registerRelationTo(m_mcParticles);
256 m_tracks.registerRelationTo(m_pidLikelihoods);
257 if (m_convertTrkExtra) m_tracks.registerRelationTo(m_belleTrkExtra);
258 m_eclClusters.registerRelationTo(m_mcParticles);
259 m_tracks.registerRelationTo(m_eclClusters);
260 m_klmClusters.registerRelationTo(m_tracks);
261 m_klmClusters.registerRelationTo(m_eclClusters);
262 m_particles.registerRelationTo(m_mcParticles);
263 m_particles.registerRelationTo(m_pidLikelihoods);
264
265 B2DEBUG(99, "[B2BIIConvertMdstModule::initializeDataStore] initialization of DataStore ended");
266}
Belle2::StoreAccessorBase::registerInDataStore
bool registerInDataStore(DataStore::EStoreFlags storeFlags=DataStore::c_WriteOut)
Register the object/array in the DataStore.
Definition: StoreAccessorBase.h:52
Belle2::StoreArray::registerRelationTo
bool registerRelationTo(const StoreArray< TO > &toArray, DataStore::EDurability durability=DataStore::c_Event, DataStore::EStoreFlags storeFlags=DataStore::c_WriteOut, const std::string &namedRelation="") const
Register a relation to the given StoreArray.
Definition: StoreArray.h:140

◆ recoverMoreThan24bitIDHEP()

int recoverMoreThan24bitIDHEP ( int  id)
private

Helper function to recover falsely set idhep info in GenHepEvt list.

Definition at line 2077 of file B2BIIConvertMdstModule.cc.

2078{
2079 /*
2080 QUICK CHECK: most of the normal particles are smaller than
2081 0x100000, while all the corrupt id has some of the high bits on.
2082
2083 This bit check has to be revised when the table below is updated.
2084 */
2085 const int mask = 0x00f00000;
2086 int high_bits = id & mask;
2087 if (high_bits == 0 || high_bits == mask) return id;
2088
2089 switch (id) {
2090 case 7114363:
2091 return 91000443; // X(3940)
2092 case 6114363:
2093 return 90000443; // Y(3940)
2094 case 6114241:
2095 return 90000321; // K_0*(800)+
2096 case 6114231:
2097 return 90000311; // K_0*(800)0
2098 case -6865004:
2099 return 9912212; // p_diff+
2100 case -6865104:
2101 return 9912112; // n_diffr
2102 case -6866773:
2103 return 9910443; // psi_diff
2104 case -6866883:
2105 return 9910333; // phi_diff
2106 case -6866993:
2107 return 9910223; // omega_diff
2108 case -6867005:
2109 return 9910211; // pi_diff+
2110 case -6867103:
2111 return 9910113; // rho_diff0
2112 case -7746995:
2113 return 9030221; // f_0(1500)
2114 case -7756773:
2115 return 9020443; // psi(4415)
2116 case -7756995:
2117 return 9020221; // eta(1405)
2118 case -7766773:
2119 return 9010443; // psi(4160)
2120 case -7776663:
2121 return 9000553; // Upsilon(5S)
2122 case -7776773:
2123 return 9000443; // psi(4040)
2124 case -7776783:
2125 return 9000433; // D_sj(2700)+
2126 case -7776995:
2127 return 9000221; // f_0(600)
2128 case -6114241:
2129 return -90000321; // K_0*(800)-
2130 case -6114231:
2131 return -90000311; // anti-K_0*(800)0
2132 case 6865004:
2133 return -9912212; // anti-p_diff-
2134 case 6865104:
2135 return -9912112; // anti-n_diffr
2136 case 6867005:
2137 return -9910211; // pi_diff-
2138 case 7776783:
2139 return -9000433; // D_sj(2700)-
2140 default:
2141 return id;
2142 }
2143}

◆ setAbortLevel()

void setAbortLevel ( int  abortLevel)
inherited

Configure the abort log level.

Definition at line 67 of file Module.cc.

68{
69 m_logConfig.setAbortLevel(static_cast<LogConfig::ELogLevel>(abortLevel));
70}
Belle2::LogConfig::ELogLevel
ELogLevel
Definition of the supported log levels.
Definition: LogConfig.h:26
Belle2::LogConfig::setAbortLevel
void setAbortLevel(ELogLevel abortLevel)
Configure the abort level.
Definition: LogConfig.h:112

◆ setDebugLevel()

void setDebugLevel ( int  debugLevel)
inherited

Configure the debug messaging level.

Definition at line 61 of file Module.cc.

62{
63 m_logConfig.setDebugLevel(debugLevel);
64}
Belle2::LogConfig::setDebugLevel
void setDebugLevel(int debugLevel)
Configure the debug messaging level.
Definition: LogConfig.h:98

◆ setDescription()

void setDescription ( const std::string &  description)
protectedinherited

Sets the description of the module.

Parameters
descriptionA description of the module.

Definition at line 214 of file Module.cc.

215{
216 m_description = description;
217}

◆ setKLMClustersRelations()

void setKLMClustersRelations ( )
private

Sets KLMCluster -> Track and ECLCluster relations.

Definition at line 2051 of file B2BIIConvertMdstModule.cc.

2052{
2053 // Relations
2054 RelationArray klmClustersToTracks(m_klmClusters, m_tracks);
2055 RelationArray klmClustersToEclClusters(m_klmClusters, m_eclClusters);
2056
2057 Belle::Mdst_klm_cluster_Manager& klm_cluster_manager = Belle::Mdst_klm_cluster_Manager::get_manager();
2058
2059
2060 for (Belle::Mdst_klm_cluster_Manager::iterator klmC_Ite = klm_cluster_manager.begin(); klmC_Ite != klm_cluster_manager.end();
2061 ++klmC_Ite) {
2062
2063 Belle::Mdst_klm_cluster mdstKlm_cluster = *klmC_Ite;
2064 Belle::Mdst_trk mTRK = mdstKlm_cluster.trk();
2065 Belle::Mdst_ecl mECL = mdstKlm_cluster.ecl();
2066
2067 if (mTRK) klmClustersToTracks.add(mdstKlm_cluster.get_ID() - 1, mTRK.get_ID() - 1);
2068 if (mECL) klmClustersToEclClusters.add(mdstKlm_cluster.get_ID() - 1, mECL.get_ID() - 1);
2069 }
2070}

◆ setLikelihoods()

void setLikelihoods ( PIDLikelihood *  pid,
Const::EDetector  det,
double  likelihoods[c_nHyp],
bool  discard_allzero = false 
)
private

Add given Belle likelihoods (not log-likelihoods, in Belle hypothesis order) for given detector to pid.

Definition at line 1396 of file B2BIIConvertMdstModule.cc.

1398{
1399 if (discard_allzero) {
1400 const double max_l = *std::max_element(likelihoods, likelihoods + c_nHyp);
1401 if (max_l <= 0.0) {
1402 return; //likelihoods broken, ignore
1403 }
1404 }
1405
1406 for (int i = 0; i < c_nHyp; i++) {
1407 float logl = log(likelihoods[i]);
1408 pid->setLogLikelihood(det, c_belleHyp_to_chargedStable[i], logl);
1409 }
1410 //copy proton likelihood to deuterons
1411 pid->setLogLikelihood(det, Const::deuteron, pid->getLogL(Const::proton, det));
1412}
Belle2::Const::deuteron
static const ChargedStable deuteron
deuteron particle
Definition: Const.h:664

◆ setLogConfig()

void setLogConfig ( const LogConfig &  logConfig)
inlineinherited

Set the log system configuration.

Definition at line 229 of file Module.h.

229{m_logConfig = logConfig;}

◆ setLogInfo()

void setLogInfo ( int  logLevel,
unsigned int  logInfo 
)
inherited

Configure the printed log information for the given level.

Parameters
logLevelThe log level (one of LogConfig::ELogLevel)
logInfoWhat kind of info should be printed? ORed combination of LogConfig::ELogInfo flags.

Definition at line 73 of file Module.cc.

74{
75 m_logConfig.setLogInfo(static_cast<LogConfig::ELogLevel>(logLevel), logInfo);
76}
Belle2::LogConfig::setLogInfo
void setLogInfo(ELogLevel logLevel, unsigned int logInfo)
Configure the printed log information for the given level.
Definition: LogConfig.h:127

◆ setLogLevel()

void setLogLevel ( int  logLevel)
inherited

Configure the log level.

Definition at line 55 of file Module.cc.

56{
57 m_logConfig.setLogLevel(static_cast<LogConfig::ELogLevel>(logLevel));
58}
Belle2::LogConfig::setLogLevel
void setLogLevel(ELogLevel logLevel)
Configure the log level.
Definition: LogConfig.cc:25

◆ setName()

void setName ( const std::string &  name)
inlineinherited

Set the name of the module.

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

Definition at line 213 of file Module.h.

213{ m_name = name; };

◆ setParamList()

void setParamList ( const ModuleParamList &  params)
inlineprotectedinherited

Replace existing parameter list.

Definition at line 500 of file Module.h.

500{ m_moduleParamList = params; }

◆ setParamPython()

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

Implements a method for setting boost::python objects.

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

Parameters
nameThe unique name of the parameter.
pyObjThe object which should be converted and stored as the parameter value.

Definition at line 234 of file Module.cc.

235{
236 LogSystem& logSystem = LogSystem::Instance();
237 logSystem.updateModule(&(getLogConfig()), getName());
238 try {
239 m_moduleParamList.setParamPython(name, pyObj);
240 } catch (std::runtime_error& e) {
241 throw std::runtime_error("Cannot set parameter '" + name + "' for module '"
242 + m_name + "': " + e.what());
243 }
244
245 logSystem.updateModule(nullptr);
246}
Belle2::LogSystem
Class for logging debug, info and error messages.
Definition: LogSystem.h:46
Belle2::LogSystem::updateModule
void updateModule(const LogConfig *moduleLogConfig=nullptr, const std::string &moduleName="")
Sets the log configuration to the given module log configuration and sets the module name This method...
Definition: LogSystem.h:191
Belle2::LogSystem::Instance
static LogSystem & Instance()
Static method to get a reference to the LogSystem instance.
Definition: LogSystem.cc:28
Belle2::ModuleParamList::setParamPython
void setParamPython(const std::string &name, const PythonObject &pyObj)
Implements a method for setting boost::python objects.
Definition: ModuleParamList.templateDetails.h:104

◆ setParamPythonDict()

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

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

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

Parameters
dictionaryThe python dictionary from which the parameter values are read.

Definition at line 249 of file Module.cc.

250{
251
252 LogSystem& logSystem = LogSystem::Instance();
253 logSystem.updateModule(&(getLogConfig()), getName());
254
255 boost::python::list dictKeys = dictionary.keys();
256 int nKey = boost::python::len(dictKeys);
257
258 //Loop over all keys in the dictionary
259 for (int iKey = 0; iKey < nKey; ++iKey) {
260 boost::python::object currKey = dictKeys[iKey];
261 boost::python::extract<std::string> keyProxy(currKey);
262
263 if (keyProxy.check()) {
264 const boost::python::object& currValue = dictionary[currKey];
265 setParamPython(keyProxy, currValue);
266 } else {
267 B2ERROR("Setting the module parameters from a python dictionary: invalid key in dictionary!");
268 }
269 }
270
271 logSystem.updateModule(nullptr);
272}
Belle2::Module::setParamPython
void setParamPython(const std::string &name, const boost::python::object &pyObj)
Implements a method for setting boost::python objects.
Definition: Module.cc:234

◆ setPropertyFlags()

void setPropertyFlags ( unsigned int  propertyFlags)
inherited

Sets the flags for the module properties.

Parameters
propertyFlagsbitwise OR of EModulePropFlags

Definition at line 208 of file Module.cc.

209{
210 m_propertyFlags = propertyFlags;
211}

◆ setReturnValue() [1/2]

void setReturnValue ( bool  value)
protectedinherited

Sets the return value for this module as bool.

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

Parameters
valueThe value of the return value.

Definition at line 227 of file Module.cc.

228{
229 m_hasReturnValue = true;
230 m_returnValue = value;
231}

◆ setReturnValue() [2/2]

void setReturnValue ( int  value)
protectedinherited

Sets the return value for this module as integer.

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

Parameters
valueThe value of the return value.

Definition at line 220 of file Module.cc.

221{
222 m_hasReturnValue = true;
223 m_returnValue = value;
224}

◆ setTracksToECLClustersRelations()

void setTracksToECLClustersRelations ( )
private

Sets Track -> ECLCluster relations.

Definition at line 2008 of file B2BIIConvertMdstModule.cc.

2009{
2010 // Relations
2011 RelationArray tracksToECLClusters(m_tracks, m_eclClusters);
2012
2013 Belle::Mdst_ecl_trk_Manager& m = Belle::Mdst_ecl_trk_Manager::get_manager();
2014 Belle::Mdst_charged_Manager& chgMg = Belle::Mdst_charged_Manager::get_manager();
2015
2016 // We first insert relations to tracks which are directly matched (type == 1)
2017 // secondly we had CR matched tracks (connected region) (type == 2)
2018 // finally tracks which are geometrically matched (type == 0)
2019 std::vector<int> insert_order_types = {1, 2, 0};
2020 for (auto& insert_type : insert_order_types) {
2021 for (Belle::Mdst_ecl_trk_Manager::iterator ecltrkIterator = m.begin(); ecltrkIterator != m.end(); ++ecltrkIterator) {
2022 Belle::Mdst_ecl_trk mECLTRK = *ecltrkIterator;
2023
2024 if (mECLTRK.type() != insert_type)
2025 continue;
2026
2027 Belle::Mdst_ecl mdstEcl = mECLTRK.ecl();
2028 Belle::Mdst_trk mTRK = mECLTRK.trk();
2029
2030 if (!mdstEcl)
2031 continue;
2032
2033 // the numbering in mdst_charged
2034 // not necessarily the same as in mdst_trk
2035 // therefore have to find corresponding mdst_charged
2036 for (Belle::Mdst_charged_Manager::iterator chgIterator = chgMg.begin(); chgIterator != chgMg.end(); ++chgIterator) {
2037 Belle::Mdst_charged mChar = *chgIterator;
2038 Belle::Mdst_trk mTRK_in_charged = mChar.trk();
2039
2040 if (mTRK_in_charged.get_ID() == mTRK.get_ID()) {
2041 // found the correct mdst_charged
2042 tracksToECLClusters.add(mChar.get_ID() - 1, mdstEcl.get_ID() - 1, 1.0);
2043 break;
2044 }
2045 }
2046 }
2047 }
2048}

◆ setType()

void setType ( const std::string &  type)
protectedinherited

Set the module type.

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

Definition at line 48 of file Module.cc.

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

◆ terminate()

void terminate ( void  )
overridevirtual

Terminates the module.

Reimplemented from Module.

Definition at line 2353 of file B2BIIConvertMdstModule.cc.

2354{
2355 B2DEBUG(99, "B2BIIConvertMdst: terminate called");
2356}

◆ testMCRelation()

void testMCRelation ( const Belle::Gen_hepevt &  belleMC,
const MCParticle *  mcP,
const std::string &  objectName 
)
private

Checks if the reconstructed object (Track, ECLCluster, ...) was matched to the same MCParticle.

Definition at line 2145 of file B2BIIConvertMdstModule.cc.

2146{
2147 int bellePDGCode = belleMC.idhep();
2148 int belleIIPDGCode = mcP->getPDG();
2149
2150 if (bellePDGCode == 0)
2151 B2WARNING("[B2BIIConvertMdstModule] " << objectName << " matched to Gen_hepevt with idhep = 0.");
2152
2153 if (bellePDGCode != belleIIPDGCode)
2154 B2WARNING("[B2BIIConvertMdstModule] " << objectName << " matched to different MCParticle! " << bellePDGCode << " vs. " <<
2155 belleIIPDGCode);
2156
2157 const double belleMomentum[] = { belleMC.PX(), belleMC.PY(), belleMC.PZ() };
2158 const double belle2Momentum[] = { mcP->get4Vector().Px(), mcP->get4Vector().Py(), mcP->get4Vector().Pz() };
2159
2160 for (unsigned i = 0; i < 3; i++) {
2161 double relDev = (belle2Momentum[i] - belleMomentum[i]) / belleMomentum[i];
2162
2163 if (relDev > 1e-3) {
2164 B2WARNING("[B2BIIConvertMdstModule] " << objectName << " matched to different MCParticle!");
2165 B2INFO(" - Gen_hepevt [" << bellePDGCode << "] px/py/pz = " << belleMC.PX() << "/" << belleMC.PY() << "/" << belleMC.PZ());
2166 B2INFO(" - TrackFitResult [" << belleIIPDGCode << "] px/py/pz = " << mcP->get4Vector().Px() << "/" << mcP->get4Vector().Py() << "/"
2167 << mcP->get4Vector().Pz());
2168 }
2169 }
2170}
Belle2::MCParticle::get4Vector
ROOT::Math::PxPyPzEVector get4Vector() const
Return 4Vector of particle.
Definition: MCParticle.h:196
Belle2::MCParticle::getPDG
int getPDG() const
Return PDG code of particle.
Definition: MCParticle.h:101

Member Data Documentation

◆ BFIELD

const double BFIELD = 1.5
private

B filed in TESLA.

Definition at line 438 of file B2BIIConvertMdstModule.h.

◆ c_belleHyp_to_chargedStable

const Const::ChargedStable c_belleHyp_to_chargedStable = { Const::electron, Const::muon, Const::pion, Const::kaon, Const::proton }
staticprivate

maps Belle hypotheses to Const::ChargedStable (from http://belle.kek.jp/secured/wiki/doku.php?id=software:atc_pid).

Definition at line 323 of file B2BIIConvertMdstModule.h.

◆ c_nHyp

const int c_nHyp = 5
staticprivate

Number of Belle track hypotheses (see c_belleHyp_to_chargedStable).

Definition at line 321 of file B2BIIConvertMdstModule.h.

◆ ERRMCONV

const int ERRMCONV[7] = {3, 4, 5, -1, 0, 1, 2}
private

CONVERSION OF TRACK ERROR MATRIX ELEMENTS.

Belle error matrix elements are in the following order px/py/pz/E/x/y/z or in indices are (0,1,2,3,4,5,6) Belle II error matrix elements are in the following order x/y/z/px/py/pz or in indices (0,1,2,3,4,5). The conversion of indices is therefore: 4 -> 0, 5 -> 1, 6 -> 2, 0 -> 3, 1 -> 4, 2 -> 5, and 3 -> nothing.

Definition at line 435 of file B2BIIConvertMdstModule.h.

◆ genHepevtToMCParticle

std::map<int, int> genHepevtToMCParticle
private

map of Gen_hepevt Panther IDs and corresponding MCParticle StoreArray indices

Definition at line 370 of file B2BIIConvertMdstModule.h.

◆ KAPPA2OMEGA

const double KAPPA2OMEGA = 1.5 * TMath::C() * 1E-11
private

Conversion factor for Kappa -> Omega helix parameters.

Definition at line 441 of file B2BIIConvertMdstModule.h.

◆ m_beamSpot

BeamSpot m_beamSpot
private

Interaction Point of the beam.

Definition at line 412 of file B2BIIConvertMdstModule.h.

◆ m_beamSpotDB

OptionalDBObjPtr<BeamSpot> m_beamSpotDB
private

BeamSpot for IP.

Definition at line 411 of file B2BIIConvertMdstModule.h.

◆ m_belleTrkExtra

StoreArray<BelleTrkExtra> m_belleTrkExtra
private

Belle CDC extra information.

Definition at line 402 of file B2BIIConvertMdstModule.h.

◆ m_collisionAxisCMS

CollisionAxisCMS m_collisionAxisCMS
private

CollisionAxisCMS of the beam.

Definition at line 424 of file B2BIIConvertMdstModule.h.

◆ m_collisionAxisCMSDB

OptionalDBObjPtr<CollisionAxisCMS> m_collisionAxisCMSDB
private

CollisionAxisCMS.

Definition at line 423 of file B2BIIConvertMdstModule.h.

◆ m_collisionBoostVector

CollisionBoostVector m_collisionBoostVector
private

CollisionBoostVector for bosst vector of the beam.

Definition at line 416 of file B2BIIConvertMdstModule.h.

◆ m_collisionBoostVectorDB

OptionalDBObjPtr<CollisionBoostVector> m_collisionBoostVectorDB
private

CollisionBoostVector for boost vector.

Definition at line 415 of file B2BIIConvertMdstModule.h.

◆ m_collisionInvM

CollisionInvariantMass m_collisionInvM
private

CollisionInvariantMass for the invariant mass of the beam.

Definition at line 420 of file B2BIIConvertMdstModule.h.

◆ m_collisionInvMDB

OptionalDBObjPtr<CollisionInvariantMass> m_collisionInvMDB
private

CollisionInvariantMass for Invariant Mass of Beam.

Definition at line 419 of file B2BIIConvertMdstModule.h.

◆ m_conditions

std::vector<ModuleCondition> m_conditions
privateinherited

Module condition, only non-null if set.

Definition at line 520 of file Module.h.

◆ m_convertBeamParameters

bool m_convertBeamParameters
private

Convert beam parameters or use information stored in Belle II database.

Definition at line 138 of file B2BIIConvertMdstModule.h.

◆ m_convertEvtcls

bool m_convertEvtcls
private

Flag to switch on conversion of Evtcls table.

Definition at line 151 of file B2BIIConvertMdstModule.h.

◆ m_convertNbar

bool m_convertNbar
private

Flag to create anti-n0:mdst list from gamma:mdst.

Definition at line 159 of file B2BIIConvertMdstModule.h.

◆ m_convertRecTrg

bool m_convertRecTrg
private

Flag to switch on conversion of rectrg_summary3.

Definition at line 155 of file B2BIIConvertMdstModule.h.

◆ m_convertTrkExtra

bool m_convertTrkExtra
private

Flag to switch on conversion of first(last)_{x,y,z} of mdst_trk_fit.

Definition at line 157 of file B2BIIConvertMdstModule.h.

◆ m_description

std::string m_description
privateinherited

The description of the module.

Definition at line 510 of file Module.h.

◆ m_eclClusters

StoreArray<ECLCluster> m_eclClusters
private

ECL clusters.

Definition at line 384 of file B2BIIConvertMdstModule.h.

◆ m_evtInfo

StoreObjPtr<EventExtraInfo> m_evtInfo
private

Event Extra Info.

Definition at line 408 of file B2BIIConvertMdstModule.h.

◆ m_hasReturnValue

bool m_hasReturnValue
privateinherited

True, if the return value is set.

Definition at line 517 of file Module.h.

◆ m_klmClusters

StoreArray<KLMCluster> m_klmClusters
private

KLM clusters.

Definition at line 387 of file B2BIIConvertMdstModule.h.

◆ m_lastIPProfileBin

int m_lastIPProfileBin { -1}
private

variable to tell us which IPProfile bin was active last time we looked

Definition at line 168 of file B2BIIConvertMdstModule.h.

◆ m_logConfig

LogConfig m_logConfig
privateinherited

The log system configuration of the module.

Definition at line 513 of file Module.h.

◆ m_matchType2E9oE25Threshold

double m_matchType2E9oE25Threshold
private

E9/E25 threshold value clusters with a value above this threshold are classified as neutral even if tracks are matched to their connected region (matchType == 2)

Definition at line 165 of file B2BIIConvertMdstModule.h.

◆ m_mcMatchingMode

MCMatchingMode m_mcMatchingMode
private

C matching mode.

Definition at line 149 of file B2BIIConvertMdstModule.h.

◆ m_mcMatchingModeString

std::string m_mcMatchingModeString
private

MC matching mode.

Definition at line 146 of file B2BIIConvertMdstModule.h.

◆ m_mcParticles

StoreArray<MCParticle> m_mcParticles
private

MC particles.

Definition at line 381 of file B2BIIConvertMdstModule.h.

◆ m_moduleParamList

ModuleParamList m_moduleParamList
privateinherited

List storing and managing all parameter of the module.

Definition at line 515 of file Module.h.

◆ m_name

std::string m_name
privateinherited

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

Definition at line 507 of file Module.h.

◆ m_nisEnable

bool m_nisEnable
private

Flag to switch on conversion of nisKsFinder info.

Definition at line 153 of file B2BIIConvertMdstModule.h.

◆ m_package

std::string m_package
privateinherited

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

Definition at line 509 of file Module.h.

◆ m_particleGraph

Belle2::MCParticleGraph m_particleGraph
private

MCParticle Graph to build Belle2 MC Particles.

Definition at line 367 of file B2BIIConvertMdstModule.h.

◆ m_particles

StoreArray<Particle> m_particles
private

Particles.

Definition at line 399 of file B2BIIConvertMdstModule.h.

◆ m_pidLikelihoods

StoreArray<PIDLikelihood> m_pidLikelihoods
private

output PIDLikelihood array.

Definition at line 405 of file B2BIIConvertMdstModule.h.

◆ m_propertyFlags

unsigned int m_propertyFlags
privateinherited

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

Definition at line 511 of file Module.h.

◆ m_realData

bool m_realData
private

flag that tells whether given data sample is for real data or MC

Definition at line 135 of file B2BIIConvertMdstModule.h.

◆ m_returnValue

int m_returnValue
privateinherited

The return value.

Definition at line 518 of file Module.h.

◆ m_trackFitResults

StoreArray<TrackFitResult> m_trackFitResults
private

Track fir results.

Definition at line 393 of file B2BIIConvertMdstModule.h.

◆ m_tracks

StoreArray<Track> m_tracks
private

Tracks.

Definition at line 390 of file B2BIIConvertMdstModule.h.

◆ m_type

std::string m_type
privateinherited

The type of the module, saved as a string.

Definition at line 508 of file Module.h.

◆ m_use6x6CovarianceMatrix4Tracks

bool m_use6x6CovarianceMatrix4Tracks
private

flag that tells which form of covariance matrix should be used in the conversion of charged tracks

Definition at line 143 of file B2BIIConvertMdstModule.h.

◆ m_v0s

StoreArray<V0> m_v0s
private

V0-particles.

Definition at line 396 of file B2BIIConvertMdstModule.h.

◆ mdstEclToECLCluster

std::map<int, int> mdstEclToECLCluster
private

map of Mdst_ecl Panther IDs and corresponding ECLCluster StoreArray indices

Definition at line 372 of file B2BIIConvertMdstModule.h.

◆ mdstGammaToParticle

std::map<int, int> mdstGammaToParticle
private

map of gamma Panther IDs and corresponding Particle StoreArray indices

Definition at line 374 of file B2BIIConvertMdstModule.h.

◆ mdstKlmToKLMCluster

std::map<int, int> mdstKlmToKLMCluster
private

map of Mdst_klm Panther IDs and corresponding KLMCluster StoreArray indices

Definition at line 376 of file B2BIIConvertMdstModule.h.

◆ mdstKlongToParticle

std::map<int, int> mdstKlongToParticle
private

map of Klong Panther IDs and corresponding Particle StoreArray indices

Definition at line 378 of file B2BIIConvertMdstModule.h.

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