Belle II Software  release-08-01-10
CDCCRTestModule Class Reference

CDC Cosmic test calibration module. More...

#include <CDCCRTestModule.h>

Inheritance diagram for CDCCRTestModule:
Collaboration diagram for CDCCRTestModule:

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

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

Static Public Member Functions

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

Protected Member Functions

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

Private Member Functions

TH1 * getHist (const char *name, const char *title, int nBins, double x0, double x1)
 Create 1D histogram.
 
TProfile * getHistProfile (const char *name, const char *title, int nBins, double x0, double x1)
 Create profile plot.
 
TH2 * getHist (const char *name, const char *title, int nBinsX, double x0, double x1, int nBinsY, double y0, double y1)
 Create 2d-histogram.
 
TH1 * getHist (const std::string &name, const std::string &title, int nBins, double x0, double x1)
 Create 1d-histogram.
 
TProfile * getHistProfile (const std::string &name, const std::string &title, int nBins, double x0, double x1)
 Create profile plot.
 
TH2 * getHist (const std::string &name, const std::string &title, int nBinsX, double x0, double x1, int nBinsY, double y0, double y1)
 Create 2d-histogram.
 
const genfit::SharedPlanePtr constructPlane (const genfit::MeasuredStateOnPlane &state, WireID m_wireID)
 Construct a plane for the hit.
 
void getResidualOfUnFittedLayer (Belle2::RecoTrack *track)
 Calculate residual for Layers which didn't use int fitting.
 
void plotResults (Belle2::RecoTrack *track)
 Plot track parameters and related variables.
 
void getHitDistInTrackCand (const RecoTrack *track)
 Make hit distribution from track candidate.
 
B2Vector3D getTriggerHitPosition (Belle2::RecoTrack *track)
 extrapolation track to trigger counter plane (y position).
 
void HitEfficiency (const Belle2::RecoTrack *track)
 Cal Hit eff.
 
int getICLayer (int slayer, int ilayer)
 Convert slayer and ilayer to iclayer.
 
std::list< ModulePtrgetModules () const override
 no submodules, return empty list
 
std::string getPathString () const override
 return the module name.
 
void setParamPython (const std::string &name, const boost::python::object &pyObj)
 Implements a method for setting boost::python objects. More...
 
void setParamPythonDict (const boost::python::dict &dictionary)
 Implements a method for reading the parameter values from a boost::python dictionary. More...
 

Private Attributes

const Belle2::TrackFitResultfitresult
 Track fit result.
 
StoreObjPtr< EventT0m_eventTimeStoreObject
 Event timing. More...
 
StoreArray< Trackm_Tracks
 Tracks.
 
StoreArray< RecoTrackm_RecoTracks
 Tracks.
 
StoreArray< TrackFitResultm_TrackFitResults
 Track fit results.
 
StoreArray< CDCHitm_CDCHits
 CDC hits.
 
std::string m_trackArrayName
 Belle2::Track StoreArray name.
 
std::string m_cdcHitArrayName
 Belle2::CDCHit StoreArray name.
 
std::string m_recoTrackArrayName
 Belle2::RecoTrack StoreArray nam.e.
 
std::string m_trackFitResultArrayName
 Belle2::TrackFitResult StoreArray name.
 
std::string m_relRecoTrackTrackName
 Relation between RecoTrack and Belle2:Track.
 
std::string m_histogramDirectoryName
 subdir where to place the histograms.
 
std::vector< TH1 * > m_allHistos
 A list of 1d histograms.
 
std::string m_treeName
 Name of tree for the output file.
 
TTree * m_tree
 Output tree recording the information of each hit.
 
TH1 * m_hNTracks
 Number of track fitted, Convergence, not conv, not fit.
 
TH1 * m_hNTracksPerEvent
 Number of TrackCand per Event.
 
TH1 * m_hNTracksPerEventFitted
 Number of TrackCand per Event.
 
TH1 * m_hNDF
 Number of Degree Freedom.
 
TH1 * m_hNHits
 Number of Hits per track.
 
TH1 * m_hNHits_trackcand
 Number of Hits per trackCand.
 
TH1 * m_hChi2
 Chi2 histo.
 
TH1 * m_hPval
 Fit Probability histo.
 
TH2 * m_hNDFChi2
 Chi2 vs degree-of-freedom histo.
 
TH2 * m_hNDFPval
 Degree-of-freedom vs Probability histo.
 
TH1 * m_hAlpha
 Alpha of each Hit.
 
TH1 * m_hPhi0
 Phi0 of ttrack, see Helix.
 
TH1 * m_hTheta
 Theta of each Hit.
 
TH1 * m_hHitDistInCDCHit [56]
 Hit Dist. More...
 
TH1 * m_hHitDistInTrCand [56]
 Hit Dist. More...
 
TH1 * m_hHitDistInTrack [56]
 Hit Dist. More...
 
TH1 * m_hResidualU [56]
 Residual distribution (in cm)

 
TH1 * m_hEvtT0
 Event T0.
 
TH2 * m_hNDFResidualU [56]
 Residual vs. More...
 
TH1 * m_hNormalizedResidualU [56]
 Residual distribution normalized with tracking error. More...
 
TH2 * m_hNDFNormalizedResidualU [56]
 Normalized residual vs. More...
 
TH2 * m_hDxDt [56]
 Unbiased x_fit vs. More...
 
TProfile * m_hHitEff_soft [56]
 Hit efficience of each layer, software.
 
TH2 * m_h2DHitDistInCDCHit
 2D Hit Dist. More...
 
TH2 * m_h2DHitDistInTrCand
 2D Hit Dist. More...
 
TH2 * m_h2DHitDistInTrack
 2D Hit Dist..(ICLay vs IWire) have weight>0.5 after fit with DAF
 
TH2 * m_hTriggerHitZX
 Trigger hit image.
 
double res_b
 Biased residual.
 
double res_u
 Unbiased residual.
 
double res_b_err
 Biased residual error.
 
double res_u_err
 Unbiased residual error.
 
double weight
 Weight of hit.
 
double absRes_b
 absolute value of biased residual.
 
double absRes_u
 absolute value of unbiased residual.
 
double alpha
 Entrance Azimuthal angle of hit (degree).
 
double theta
 Entrance Polar angle of hit (degree).
 
unsigned short adc
 adc value.
 
short tdc
 tdc value.
 
double t
 Measurement Drift time.
 
double t_fit
 Drift time calculated from x_fit.
 
double dt_flight
 Time of flight.
 
double dt_flight_sim
 Time of flight (Simulation).
 
double dt_prop
 Time of propagation.
 
double evtT0
 Event time.
 
double Pt
 Transverse momentum.
 
double x_mea
 measure drift length (signed by left right).
 
double x_u
 X_fit for unbiased track fit.
 
double x_b
 X_fit for biased track fit.
 
double x_sim
 Simulation DriftLength .
 
double z
 Z of hit on wire.
 
double z_sim
 Z of hit on wire (simulation).
 
double z_prop
 Propagation Length along the sense wire.
 
int lay
 Layer ID.
 
int IWire
 Wire ID.
 
int lr
 Left or right.
 
int numhits
 Number of hits.
 
int boardID
 Electrical Board ID.
 
double Pval
 P-value of fitted track. More...
 
double TrPval
 P-value of fitted track. More...
 
double ndf
 degree of freedom.
 
double d0
 Track Parameter, d0.
 
double z0
 Track Parameter, z0.
 
double phi0
 Track Parameter, phi0.
 
double tanL
 Track Parameter, tanL.
 
double omega
 Track Parameter, omega.
 
double m_MinimumPt
 Minimum Transverse momentum of tracks.
 
B2Vector3D m_trigHitPos
 Trigger position.
 
std::vector< double > m_TriggerPos
 Nominal center position of trigger counter.
 
std::vector< double > m_TriggerPlaneDirection
 Nominal center position of trigger counter.
 
std::vector< double > m_TriggerSize
 Size of trigger counter (Width x length).
 
std::vector< int > m_up
 upper channel list for each board.
 
std::vector< int > m_low
 lower channel list for each board.
 
double trigHitPos_x
 X-position of track at trigger counter.
 
double trigHitPos_z
 Z-position of track at trigger counter.
 
int trighit
 Trigger hit information. More...
 
bool m_fillExpertHistos
 Fill some histogram for monitoring fit quality.
 
bool m_plotResidual
 Process track to get the hit information of fitted track.
 
bool m_hitEfficiency
 calculate hit eff or not, Haven't finished.
 
bool m_calExpectedDriftTime
 Calculate expected drift time from x_fit or not.
 
bool m_noBFit
 fit incase no magnetic Field of not, if true, NDF=4 in cal P-value
 
bool m_ToP
 Enable to correct ToP if true.
 
bool m_ToF
 Enable to correct ToF if true.
 
bool m_StoreCDCSimHitInfo
 Store CDCSimHit Information.
 
bool m_EstimateResultForUnFittedLayer
 Calculate residual for layer that we do not use in track fitting.
 
bool m_SmallerOutput
 make output smaller by ignore some variable.
 
bool m_StoreTrackParams
 Store Track parameter or not.
 
bool m_MakeHitDist
 Switch to make histograms of hit distribution.
 
bool m_EventT0Extraction
 use Event T0 extract t0 or not.
 
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< ModuleConditionm_conditions
 Module condition, only non-null if set.
 

Detailed Description

CDC Cosmic test calibration module.

Definition at line 42 of file CDCCRTestModule.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.

Member Function Documentation

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

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

◆ def_initialize()

virtual void def_initialize ( )
inlineprotectedvirtualinherited

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

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

Reimplemented in PyModule.

Definition at line 420 of file Module.h.

◆ defineHisto()

void defineHisto ( )
overridevirtual

Histogram definitions such as TH1(), TH2(), TNtuple(), TTree()....

are supposed to be placed in this function.

Reimplemented from HistoModule.

Definition at line 80 of file CDCCRTestModule.cc.

81 {
82  m_tree = new TTree(m_treeName.c_str(), "tree");
83  m_tree->Branch("x_mea", &x_mea, "x_mea/D");
84  m_tree->Branch("x_u", &x_u, "x_u/D");
85  m_tree->Branch("x_b", &x_b, "x_b/D");
86  m_tree->Branch("z", &z, "z/D");
87  m_tree->Branch("alpha", &alpha, "alpha/D");
88  m_tree->Branch("theta", &theta, "theta/D");
89  m_tree->Branch("t", &t, "t/D");
90  m_tree->Branch("evtT0", &evtT0, "evtT0/D");
91  m_tree->Branch("adc", &adc, "adc/s");
92  m_tree->Branch("boardID", &boardID, "boardID/I");
93  m_tree->Branch("lay", &lay, "lay/I");
94  m_tree->Branch("weight", &weight, "weight/D");
95  m_tree->Branch("IWire", &IWire, "IWire/I");
96  m_tree->Branch("Pval", &Pval, "Pval/D");
97  m_tree->Branch("ndf", &ndf, "ndf/D");
98  // m_tree->Branch("trighit", &trighit, "trighit/I");
99  if (m_StoreTrackParams) {
100  m_tree->Branch("d0", &d0, "d0/D");
101  m_tree->Branch("z0", &z0, "z0/D");
102  m_tree->Branch("phi0", &phi0, "phi0/D");
103  m_tree->Branch("tanL", &tanL, "tanL/D");
104  m_tree->Branch("omega", &omega, "omega/D");
105  m_tree->Branch("Pt", &Pt, "Pt/D");
106  }
107  if (m_StoreCDCSimHitInfo) {
108  m_tree->Branch("z_sim", &z_sim, "z_sim/D");
109  m_tree->Branch("x_sim", &x_sim, "x_sim/D");
110  m_tree->Branch("dt_flight_sim", &dt_flight_sim, "dt_flight_sim/D");
111  }
112  if (m_calExpectedDriftTime) { // expected drift time, calculated form xfit
113  m_tree->Branch("t_fit", &t_fit, "t_fit/D");
114  }
115  if (!m_SmallerOutput) {
116  m_tree->Branch("tdc", &tdc, "tdc/I");
117  m_tree->Branch("z_prop", &z_prop, "z_prop/D");
118  m_tree->Branch("res_b", &res_b, "res_b/D");
119  m_tree->Branch("res_u", &res_u, "res_u/D");
120  m_tree->Branch("lr", &lr, "lr/I");
121  m_tree->Branch("trigHitPos_x", &trigHitPos_x, "trigHitPos_x/D");
122  m_tree->Branch("trigHitPos_z", &trigHitPos_z, "trigHitPos_z/D");
123  m_tree->Branch("numhits", &numhits, "numhits/I");
124  m_tree->Branch("res_b_err", &res_b_err, "res_b_err/D");
125  m_tree->Branch("res_u_err", &res_u_err, "res_u_err/D");
126  m_tree->Branch("absRes_u", &absRes_u, "absRes_u/D");
127  m_tree->Branch("absRes_b", &absRes_b, "absRes_b/D");
128  m_tree->Branch("dt_prop", &dt_prop, "dt_prop/D");
129  m_tree->Branch("dt_flight", &dt_flight, "dt_flight/D");
130  }
131 
132  // int N =m_Nchannel;//Number of Wire per Layer used;
133  TDirectory* oldDir = gDirectory;
134  TDirectory* histDir = oldDir->mkdir(m_histogramDirectoryName.c_str());
135  histDir->cd();
136  m_hNTracks = getHist("hNTracks", "number of tracks", 3, 0, 3);
137  m_hNTracks->GetXaxis()->SetBinLabel(1, "fitted, converged");
138  m_hNTracks->GetXaxis()->SetBinLabel(2, "fitted, not converged");
139  m_hNTracks->GetXaxis()->SetBinLabel(3, "TrackCand, but no Track");
140 
141  m_hNDF = getHist("hNDF", "NDF of fitted track;NDF;Tracks", 71, -1, 150);
142  m_hNHits = getHist("hNHits", "#hit of fitted track;#hit;Tracks", 61, -1, 150);
143  m_hNHits_trackcand = getHist("hNHits_trackcand", "#hit of track candidate;#hit;Tracks", 71, -1, 150);
144  m_hNTracksPerEvent = getHist("hNTracksPerEvent", "#tracks/Event;#Tracks;Event", 20, 0, 20);
145  m_hNTracksPerEventFitted = getHist("hNTracksPerEventFitted", "#tracks/Event After Fit;#Tracks;Event", 20, 0, 20);
146  m_hChi2 = getHist("hChi2", "#chi^{2} of tracks;#chi^{2};Tracks", 400, 0, 400);
147  m_hPhi0 = getHist("hPhi0", "#Phi_{0} of tracks;#phi_{0} (Degree);Tracks", 400, -190, 190);
148  m_hAlpha = getHist("hAlpha", "#alpha Dist.;#alpha (Degree);Hits", 360, -90, 90);
149  m_hTheta = getHist("hTheta", "#theta Dist.;#theta (Degree);Hits", 360, 0, 180);
150  m_hPval = getHist("hPval", "p-values of tracks;pVal;Tracks", 1000, 0, 1);
151  m_hEvtT0 = getHist("hEvtT0", "Event T0; EvtT0 (ns); #event", 200, -100, 100);
152 
153  m_hTriggerHitZX = getHist("TriggerHitZX", "Hit Position on trigger counter;z(cm);x(cm)", 300, -100, 100, 120, -15, 15);
154  if (m_MakeHitDist) {
155  m_h2DHitDistInCDCHit = getHist("2DHitDistInCDCHit", " CDCHit;WireID;LayerID",
156  m_up[8] - m_low[0], m_low[0], m_up[8], 56, 0, 56);
157  m_h2DHitDistInTrCand = getHist("2DHitDistInTrCand", "Track Cand ;WireID;LayerID",
158  m_up[8] - m_low[0], m_low[0], m_up[8], 56, 0, 56);
159  m_h2DHitDistInTrack = getHist("2DHitDistInTrack", "Fitted Track ;WireID;LayerID",
160  m_up[8] - m_low[0], m_low[0], m_up[8], 56, 0, 56);
161  }
162  if (m_fillExpertHistos) {
163  m_hNDFChi2 = getHist("hNDFChi2", "#chi^{2} of tracks;NDF;#chi^{2};Tracks", 8, 0, 8, 800, 0, 200);
164  m_hNDFPval = getHist("hNDFPval", "p-values of tracks;NDF;pVal;Tracks", 8, 0, 8, 100, 0, 1);
165  }
166  int sl;
167  for (int i = 0; i < 56; ++i) {
168  if (m_hitEfficiency) {
169  m_hHitEff_soft[i] = getHistProfile(Form("hHitEff_soft_L%d", i),
170  Form("hit efficiency(soft) of Layer %d ;Drift distance;Software Efficiency", i), 200, -1, 1);
171  }
172  if (m_MakeHitDist) {
173  if (i < 8) {sl = 0;} else { sl = floor((i - 8) / 6) + 1;}
174  m_hHitDistInCDCHit[i] = getHist(Form("hHitDistInCDCHit_layer%d", i), Form("Hit Dist. ICLayer_%d;WireID;#Hits", i),
175  m_up.at(sl) - m_low.at(sl), m_low.at(sl), m_up.at(sl));
176  m_hHitDistInCDCHit[i]->SetLineColor(kGreen);
177  m_hHitDistInTrCand[i] = getHist(Form("hHitDistInTrCand_layer%d", i), Form("Hit Dist. ICLayer_%d;WireID;#Hits", i),
178  m_up.at(sl) - m_low.at(sl), m_low.at(sl), m_up.at(sl));
179  m_hHitDistInTrCand[i]->SetLineColor(kRed);
180  m_hHitDistInTrack[i] = getHist(Form("hHitDistInTrack_layer%d", i), Form("Hit Dist. ICLayer_%d;WireID;#Hits", i),
181  m_up.at(sl) - m_low.at(sl), m_low.at(sl), m_up.at(sl));
182  }
183  const double normResRange = 20;
184  const double residualRange = 0.3;
185  std::string title, name;
186  if (m_plotResidual) {
187  name = (boost::format("hist_ResidualsU%1%") % i).str();
188  title = (boost::format("unnormalized, unbiased residuals in layer %1%;cm;Tracks") % i).str();
189  m_hResidualU[i] = getHist(name, title, 500, -residualRange, residualRange);
190 
191  name = (boost::format("hNormalizedResidualsU%1%") % i).str();
192  title = (boost::format("normalized, unbiased residuals in layer %1%;NDF;#sigma (cm);Tracks") % i).str();
193  m_hNormalizedResidualU[i] = getHist(name, title, 500, -normResRange, normResRange);
194 
195  name = (boost::format("DxDt%1%") % i).str();
196  title = (boost::format("Drift Length vs Drift time at Layer_%1%;Drift Length (cm);Drift time (ns)") % i).str();
197  m_hDxDt[i] = getHist(name, title, 200, -1, 1, 450, -50, 400);
198  }
199  if (m_fillExpertHistos) {
200  name = (boost::format("hNDFResidualsU%1%") % i).str();
201  title = (boost::format("unnormalized, unbiased residuals along U in layer %1%;NDF;cm;Tracks") % i).str();
202  m_hNDFResidualU[i] = getHist(name, title, 8, 0, 8, 1000, -residualRange, residualRange);
203 
204  name = (boost::format("hNDFNormalizedResidualsU%1%") % i).str();
205  title = (boost::format("normalized, unbiased residuals in layer %1%;NDF;#sigma (cm);Tracks") % i).str();
206  m_hNDFNormalizedResidualU[i] = getHist(name, title, 8, 0, 8, 1000, -normResRange, normResRange);
207  }
208  }
209  oldDir->cd();
210 }
bool m_calExpectedDriftTime
Calculate expected drift time from x_fit or not.
double x_sim
Simulation DriftLength .
TTree * m_tree
Output tree recording the information of each hit.
double res_b
Biased residual.
double alpha
Entrance Azimuthal angle of hit (degree).
TH1 * m_hNHits
Number of Hits per track.
bool m_plotResidual
Process track to get the hit information of fitted track.
double res_u_err
Unbiased residual error.
TH1 * m_hNHits_trackcand
Number of Hits per trackCand.
double absRes_b
absolute value of biased residual.
bool m_StoreTrackParams
Store Track parameter or not.
double dt_flight
Time of flight.
bool m_SmallerOutput
make output smaller by ignore some variable.
TH1 * m_hHitDistInTrack[56]
Hit Dist.
double res_u
Unbiased residual.
TH1 * m_hPhi0
Phi0 of ttrack, see Helix.
TH1 * getHist(const char *name, const char *title, int nBins, double x0, double x1)
Create 1D histogram.
TH1 * m_hNormalizedResidualU[56]
Residual distribution normalized with tracking error.
TH1 * m_hNTracks
Number of track fitted, Convergence, not conv, not fit.
TH1 * m_hHitDistInCDCHit[56]
Hit Dist.
double z0
Track Parameter, z0.
bool m_hitEfficiency
calculate hit eff or not, Haven't finished.
std::vector< int > m_low
lower channel list for each board.
TH2 * m_hNDFNormalizedResidualU[56]
Normalized residual vs.
double Pval
P-value of fitted track.
double z_prop
Propagation Length along the sense wire.
double res_b_err
Biased residual error.
double t
Measurement Drift time.
TH2 * m_hNDFPval
Degree-of-freedom vs Probability histo.
TH1 * m_hTheta
Theta of each Hit.
TH1 * m_hNTracksPerEvent
Number of TrackCand per Event.
std::string m_histogramDirectoryName
subdir where to place the histograms.
double omega
Track Parameter, omega.
double weight
Weight of hit.
double Pt
Transverse momentum.
double x_b
X_fit for biased track fit.
double t_fit
Drift time calculated from x_fit.
bool m_fillExpertHistos
Fill some histogram for monitoring fit quality.
TH1 * m_hAlpha
Alpha of each Hit.
double z_sim
Z of hit on wire (simulation).
bool m_MakeHitDist
Switch to make histograms of hit distribution.
TH2 * m_hNDFChi2
Chi2 vs degree-of-freedom histo.
int boardID
Electrical Board ID.
double z
Z of hit on wire.
TProfile * m_hHitEff_soft[56]
Hit efficience of each layer, software.
TH2 * m_h2DHitDistInTrack
2D Hit Dist..(ICLay vs IWire) have weight>0.5 after fit with DAF
double trigHitPos_x
X-position of track at trigger counter.
TProfile * getHistProfile(const char *name, const char *title, int nBins, double x0, double x1)
Create profile plot.
double dt_prop
Time of propagation.
double x_u
X_fit for unbiased track fit.
double theta
Entrance Polar angle of hit (degree).
std::vector< int > m_up
upper channel list for each board.
TH1 * m_hNTracksPerEventFitted
Number of TrackCand per Event.
TH1 * m_hPval
Fit Probability histo.
std::string m_treeName
Name of tree for the output file.
double d0
Track Parameter, d0.
TH1 * m_hHitDistInTrCand[56]
Hit Dist.
TH2 * m_h2DHitDistInTrCand
2D Hit Dist.
TH2 * m_hDxDt[56]
Unbiased x_fit vs.
TH2 * m_hNDFResidualU[56]
Residual vs.
double tanL
Track Parameter, tanL.
unsigned short adc
adc value.
bool m_StoreCDCSimHitInfo
Store CDCSimHit Information.
double phi0
Track Parameter, phi0.
double trigHitPos_z
Z-position of track at trigger counter.
double ndf
degree of freedom.
double absRes_u
absolute value of unbiased residual.
TH1 * m_hNDF
Number of Degree Freedom.
TH1 * m_hResidualU[56]
Residual distribution (in cm)
TH2 * m_hTriggerHitZX
Trigger hit image.
TH2 * m_h2DHitDistInCDCHit
2D Hit Dist.
double dt_flight_sim
Time of flight (Simulation).
double x_mea
measure drift length (signed by left right).

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

◆ event()

void event ( void  )
overridevirtual

Event action (main routine).

find results in track fit results

Reimplemented from HistoModule.

Definition at line 238 of file CDCCRTestModule.cc.

◆ getAfterConditionPath()

Module::EAfterConditionPath getAfterConditionPath ( ) const
inherited

What to do after the conditional path is finished.

(defaults to c_End if no condition is set)

Definition at line 133 of file Module.cc.

◆ getConditionPath()

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

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


Definition at line 113 of file Module.cc.

◆ 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 RootOutputModule, StorageRootOutputModule, and RootInputModule.

Definition at line 134 of file Module.h.

◆ getName()

const std::string& getName ( ) const
inlineinherited

Returns the name of the module.

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

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

Definition at line 187 of file Module.h.

◆ getParamInfoListPython()

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

Returns a python list of all parameters.

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

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

Definition at line 279 of file Module.cc.

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

◆ getType()

const std::string & getType ( ) const
inherited

Returns the type of the module (i.e.

class name minus 'Module')

Definition at line 41 of file Module.cc.

◆ hasProperties()

bool hasProperties ( unsigned int  propertyFlags) const
inherited

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

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

Definition at line 160 of file Module.cc.

◆ if_false()

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

A simplified version to add a condition to the module.

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

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

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

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

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

◆ if_value()

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

Add a condition to the module.

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

See https://confluence.desy.de/display/BI/Software+ModCondTut or ModuleCondition for a description of the syntax.

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

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

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

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

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

◆ setParamPython()

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

Implements a method for setting boost::python objects.

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

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

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

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

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

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

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

Member Data Documentation

◆ m_eventTimeStoreObject

StoreObjPtr<EventT0> m_eventTimeStoreObject
private

Event timing.

The event time is fetched from the data store using this pointer.

Definition at line 197 of file CDCCRTestModule.h.

◆ m_h2DHitDistInCDCHit

TH2* m_h2DHitDistInCDCHit
private

2D Hit Dist.

(ICLay vs IWire) from CDCHit.

Definition at line 244 of file CDCCRTestModule.h.

◆ m_h2DHitDistInTrCand

TH2* m_h2DHitDistInTrCand
private

2D Hit Dist.

(ICLay vs IWire) of Track candidates.

Definition at line 245 of file CDCCRTestModule.h.

◆ m_hDxDt

TH2* m_hDxDt[56]
private

Unbiased x_fit vs.

drift time.

Definition at line 242 of file CDCCRTestModule.h.

◆ m_hHitDistInCDCHit

TH1* m_hHitDistInCDCHit[56]
private

Hit Dist.

from CDCHit.

Definition at line 234 of file CDCCRTestModule.h.

◆ m_hHitDistInTrack

TH1* m_hHitDistInTrack[56]
private

Hit Dist.

after fit (Weight of Hit >0.5).

Definition at line 236 of file CDCCRTestModule.h.

◆ m_hHitDistInTrCand

TH1* m_hHitDistInTrCand[56]
private

Hit Dist.

Before Fit.

Definition at line 235 of file CDCCRTestModule.h.

◆ m_hNDFNormalizedResidualU

TH2* m_hNDFNormalizedResidualU[56]
private

Normalized residual vs.

ndf.

Definition at line 241 of file CDCCRTestModule.h.

◆ m_hNDFResidualU

TH2* m_hNDFResidualU[56]
private

Residual vs.

ndf.

Definition at line 239 of file CDCCRTestModule.h.

◆ m_hNormalizedResidualU

TH1* m_hNormalizedResidualU[56]
private

Residual distribution normalized with tracking error.


Definition at line 240 of file CDCCRTestModule.h.

◆ Pval

double Pval
private

P-value of fitted track.


Definition at line 281 of file CDCCRTestModule.h.

◆ trighit

int trighit
private

Trigger hit information.

1 if track hits trigger counter, otherwise 0.

Definition at line 299 of file CDCCRTestModule.h.

◆ TrPval

double TrPval
private

P-value of fitted track.


Definition at line 282 of file CDCCRTestModule.h.


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