 |
Belle II Software
release-05-02-19
|
11 #ifndef CDCCRTestModule_H_
12 #define CDCCRTestModule_H_
14 #include <framework/core/HistoModule.h>
16 #include <cdc/dataobjects/WireID.h>
20 #include <mdst/dataobjects/TrackFitResult.h>
21 #include <tracking/dataobjects/RecoTrack.h>
22 #include <framework/dataobjects/EventT0.h>
42 class CDCCRTestModule :
public HistoModule {
70 void event()
override;
93 TH1*
getHist(
const char* name,
const char* title,
94 int nBins,
double x0,
double x1)
96 TH1* h =
new TH1D(name, title, nBins, x0, x1);
105 int nBins,
double x0,
double x1)
107 TProfile* h =
new TProfile(name, title, nBins, x0, x1);
115 TH2*
getHist(
const char* name,
const char* title,
116 int nBinsX,
double x0,
double x1,
117 int nBinsY,
double y0,
double y1)
119 TH2* h =
new TH2D(name, title, nBinsX, x0, x1, nBinsY, y0, y1);
127 TH1*
getHist(
const std::string& name,
const std::string& title,
128 int nBins,
double x0,
double x1)
130 return getHist(name.c_str(), title.c_str(), nBins, x0, x1);
136 TProfile*
getHistProfile(
const std::string& name,
const std::string& title,
137 int nBins,
double x0,
double x1)
139 return getHistProfile(name.c_str(), title.c_str(), nBins, x0, x1);
145 TH2*
getHist(
const std::string& name,
const std::string& title,
146 int nBinsX,
double x0,
double x1,
147 int nBinsY,
double y0,
double y1)
149 return getHist(name.c_str(), title.c_str(), nBinsX, x0, x1, nBinsY, y0, y1);
188 if (slayer == 0) {
return ilayer;}
189 else {
return 8 + (slayer - 1) * 6 + ilayer;}
TTree * m_tree
Output tree recording the information of each hit.
double res_u
Unbiased residual.
double dt_prop
Time of propagation.
const genfit::SharedPlanePtr constructPlane(const genfit::MeasuredStateOnPlane &state, WireID m_wireID)
Construct a plane for the hit.
double res_b
Biased residual.
void terminate() override
Termination action.
TH2 * m_hNDFNormalizedResidualU[56]
Normalized residual vs.
std::shared_ptr< genfit::DetPlane > SharedPlanePtr
Shared Pointer to a DetPlane.
double phi0
Track Parameter, phi0.
TH1 * m_hHitDistInTrack[56]
Hit Dist.
std::string m_recoTrackArrayName
Belle2::RecoTrack StoreArray nam.e.
TH1 * m_hHitDistInTrCand[56]
Hit Dist.
void initialize() override
Initializes the Module.
double alpha
Entrance Azimuthal angle of hit (degree).
TVector3 getTriggerHitPosition(Belle2::RecoTrack *track)
extrapolation track to trigger counter plane (y position).
double z_prop
Propagation Length along the sense wire.
TH1 * m_hNTracks
Number of track fitted, Convergence, not conv, not fit.
bool m_noBFit
fit incase no magnetic Field of not, if true, NDF=4 in cal P-value
TH1 * getHist(const char *name, const char *title, int nBins, double x0, double x1)
Create 1D histogram.
std::string m_treeName
Name of tree for the output file.
bool m_MakeHitDist
Switch to make histograms of hit distribution.
#StateOnPlane with additional covariance matrix.
double Pval
P-value of fitted track.
double dt_flight_sim
Time of flight (Simulation).
TH1 * m_hNTracksPerEvent
Number of TrackCand per Event.
int numhits
Number of hits.
double d0
Track Parameter, d0.
int getICLayer(int slayer, int ilayer)
Convert slayer and ilayer to iclayer.
double m_MinimumPt
Minimum Transverse momentum of tracks.
void getResidualOfUnFittedLayer(Belle2::RecoTrack *track)
Calculate residual for Layers which didn't use int fitting.
TH2 * m_h2DHitDistInCDCHit
2D Hit Dist.
double weight
Weight of hit.
double Pt
Transverse momentum.
double trigHitPos_x
X-position of track at trigger counter.
TVector3 m_trigHitPos
Trigger position.
bool m_plotResidual
Process track to get the hit information of fitted track.
bool m_SmallerOutput
make output smaller by ignore some variable.
bool m_EventT0Extraction
use Event T0 extract t0 or not.
bool m_fillExpertHistos
Fill some histogram for monitoring fit quality.
virtual ~CDCCRTestModule()
Destructor.
int trighit
Trigger hit information.
Values of the result of a track fit with a given particle hypothesis.
bool m_StoreTrackParams
Store Track parameter or not.
TH1 * m_hPval
Fit Probability histo.
StoreObjPtr< EventT0 > m_eventTimeStoreObject
Event timing.
std::vector< double > m_TriggerPos
Nominal center position of trigger counter.
TH1 * m_hNDF
Number of Degree Freedom.
This is the Reconstruction Event-Data Model Track.
double t_fit
Drift time calculated from x_fit.
double x_sim
Simulation DriftLength .
std::string m_trackArrayName
Belle2::Track StoreArray name.
unsigned short adc
adc value.
TH2 * m_h2DHitDistInTrack
2D Hit Dist..(ICLay vs IWire) have weight>0.5 after fit with DAF
std::string m_relRecoTrackTrackName
Relation between RecoTrack and Belle2:Track.
Abstract base class for different kinds of events.
double z0
Track Parameter, z0.
double x_b
X_fit for biased track fit.
TH2 * m_hNDFChi2
Chi2 vs degree-of-freedom histo.
TProfile * m_hHitEff_soft[56]
Hit efficience of each layer, software.
const Belle2::TrackFitResult * fitresult
Track fit result.
bool m_EstimateResultForUnFittedLayer
Calculate residual for layer that we do not use in track fitting.
TH1 * m_hAlpha
Alpha of each Hit.
int boardID
Electrical Board ID.
void endRun() override
End run action.
TH2 * m_hNDFPval
Degree-of-freedom vs Probability histo.
double absRes_u
absolute value of unbiased residual.
std::vector< TH1 * > m_allHistos
A list of 1d histograms.
std::string m_cdcHitArrayName
Belle2::CDCHit StoreArray name.
std::vector< int > m_up
upper channel list for each board.
double theta
Entrance Polar angle of hit (degree).
TH1 * m_hNHits
Number of Hits per track.
double z_sim
Z of hit on wire (simulation).
TH1 * m_hNHits_trackcand
Number of Hits per trackCand.
TH2 * m_h2DHitDistInTrCand
2D Hit Dist.
void defineHisto() override
Histogram definitions such as TH1(), TH2(), TNtuple(), TTree()....
double t
Measurement Drift time.
TH1 * m_hNormalizedResidualU[56]
Residual distribution normalized with tracking error.
void event() override
Event action (main routine).
double ndf
degree of freedom.
double x_mea
measure drift length (signed by left right).
void HitEfficiency(const Belle2::RecoTrack *track)
Cal Hit eff.
TH2 * m_hTriggerHitZX
Trigger hit image.
TH1 * m_hNTracksPerEventFitted
Number of TrackCand per Event.
std::string m_trackFitResultArrayName
Belle2::TrackFitResult StoreArray name.
double TrPval
P-value of fitted track.
CDCCRTestModule()
Constructor.
double z
Z of hit on wire.
std::string m_histogramDirectoryName
subdir where to place the histograms.
void beginRun() override
Begin run action.
TH1 * m_hHitDistInCDCHit[56]
Hit Dist.
TH1 * m_hResidualU[56]
Residual distribution (in cm)
TH1 * m_hPhi0
Phi0 of ttrack, see Helix.
void plotResults(Belle2::RecoTrack *track)
Plot track parameters and related variables.
double res_b_err
Biased residual error.
bool m_calExpectedDriftTime
Calculate expected drift time from x_fit or not.
double omega
Track Parameter, omega.
std::vector< double > m_TriggerSize
Size of trigger counter (Width x length).
bool m_StoreCDCSimHitInfo
Store CDCSimHit Information.
bool m_hitEfficiency
calculate hit eff or not, Haven't finished.
bool m_ToP
Enable to correct ToP if true.
double x_u
X_fit for unbiased track fit.
void getHitDistInTrackCand(const RecoTrack *track)
Make hit distribution from track candidate.
double res_u_err
Unbiased residual error.
bool m_ToF
Enable to correct ToF if true.
double absRes_b
absolute value of biased residual.
TProfile * getHistProfile(const char *name, const char *title, int nBins, double x0, double x1)
Create profile plot.
std::vector< int > m_low
lower channel list for each board.
double trigHitPos_z
Z-position of track at trigger counter.
TH1 * m_hTheta
Theta of each Hit.
TH2 * m_hNDFResidualU[56]
Residual vs.
double tanL
Track Parameter, tanL.
TH2 * m_hDxDt[56]
Unbiased x_fit vs.
std::vector< double > m_TriggerPlaneDirection
Nominal center position of trigger counter.
double dt_flight
Time of flight.