Skip to content
Snippets Groups Projects
Commit e1719e4c authored by Valentina Akishina's avatar Valentina Akishina Committed by Florian Uhlig
Browse files

L1: add option for mc-point-based reconstruction

parent bd33f3dc
Branches
Tags
1 merge request!742L1: add mc-point based reconstruction
Expand all Show whitespace changes
Inline Side-by-side
reco/L1/CbmL1.cxx
+ 22
− 10
View file @ e1719e4c
......@@ -195,7 +195,7 @@ InitStatus CbmL1::Init()
fUseTOF = 0;
if (fTrackingMode == L1Algo::TrackingMode::kMcbm) {
fUseMUCH = 0;
fUseMUCH = 1;
fUseTRD = 1;
fUseTOF = 1;
}
......@@ -433,6 +433,8 @@ InitStatus CbmL1::Init()
if (NULL == fChannelInfo) break;
float z = fChannelInfo->GetZ();
float x = fChannelInfo->GetX();
float y = fChannelInfo->GetY();
if (fMissingHits) {
if ((x > 20) && (z > 270) && (station == 1)) station = 2;
if (z > 400) continue;
......@@ -534,11 +536,12 @@ InitStatus CbmL1::Init()
geo.push_back(t.RadLength);
fscal f_phi = 0, f_sigma = mvdStationPar->GetXRes(ist) / 10000, b_phi = 3.14159265358 / 2.,
b_sigma = mvdStationPar->GetYRes(ist) / 10000;
b_sigma = mvdStationPar->GetYRes(ist) / 10000, dt = 1000;
geo.push_back(f_phi);
geo.push_back(f_sigma);
geo.push_back(b_phi);
geo.push_back(b_sigma);
geo.push_back(dt);
z = t.z;
Xmax = Ymax = t.R;
......@@ -564,6 +567,7 @@ InitStatus CbmL1::Init()
b_phi += station->GetSensorStereoAngle(1) * Pi / 180.;
f_sigma = station->GetSensorPitch(0) / TMath::Sqrt(12);
b_sigma = f_sigma;
fscal dt = 5;
//f_sigma *= cos(f_phi); // TODO: think about this
//b_sigma *= cos(b_phi);
......@@ -571,6 +575,7 @@ InitStatus CbmL1::Init()
geo.push_back(f_sigma);
geo.push_back(b_phi);
geo.push_back(b_sigma);
geo.push_back(dt);
z = station->GetZ();
Xmax = station->GetXmax();
......@@ -597,15 +602,17 @@ InitStatus CbmL1::Init()
geo.push_back(2);
geo.push_back(z);
geo.push_back(layer->GetDz());
geo.push_back(0);
geo.push_back(10);
geo.push_back(100); //station->GetRmax()
geo.push_back(0);
fscal f_phi = 0, f_sigma = 0.1, b_phi = 3.14159265358 / 2., b_sigma = 0.1;
// fscal f_phi = 0, f_sigma = 0.1, b_phi = 3.14159265358 / 2., b_sigma = 0.1;
fscal f_phi = 0, f_sigma = 0.35, b_phi = 3.14159265358 / 2., b_sigma = 0.35, dt = 3.9;
geo.push_back(f_phi);
geo.push_back(f_sigma);
geo.push_back(b_phi);
geo.push_back(b_sigma);
geo.push_back(dt);
Xmax = 100; //station->GetRmax();
......@@ -644,15 +651,17 @@ InitStatus CbmL1::Init()
geo.push_back(stationZ);
geo.push_back(2 * module->GetSizeZ());
geo.push_back(0);
geo.push_back(10);
geo.push_back(2 * module->GetSizeX());
geo.push_back(10);
fscal f_phi = 0, f_sigma = 1., b_phi = 3.14159265358 / 2., b_sigma = 1.;
// fscal f_phi = 0, f_sigma = 1., b_phi = 3.14159265358 / 2., b_sigma = 1.;
fscal f_phi = 0, f_sigma = 0.15, b_phi = 3.14159265358 / 2., b_sigma = 0.15, dt = 10;
geo.push_back(f_phi);
geo.push_back(f_sigma);
geo.push_back(b_phi);
geo.push_back(b_sigma);
geo.push_back(dt);
Xmax = module->GetSizeX();
Ymax = module->GetSizeY();
LOG(info) << "L1: Trd station " << num << " at z " << stationZ << endl;
......@@ -672,11 +681,13 @@ InitStatus CbmL1::Init()
geo.push_back(150); /// TODO: add Tof max radius
geo.push_back(10);
fscal f_phi = 0, f_sigma = 1 / 10000, b_phi = 3.14159265358 / 2., b_sigma = 1 / 10000;
// fscal f_phi = 0, f_sigma = 1 / 10000, b_phi = 3.14159265358 / 2., b_sigma = 1 / 10000;
fscal f_phi = 0, f_sigma = 0.42, b_phi = 3.14159265358 / 2., b_sigma = 0.23, dt = 0.075;
geo.push_back(f_phi);
geo.push_back(f_sigma);
geo.push_back(b_phi);
geo.push_back(b_sigma);
geo.push_back(dt);
Xmax = Ymax = 20;
LOG(info) << "L1: Tof station " << num << " at z " << z << endl;
}
......@@ -875,7 +886,7 @@ InitStatus CbmL1::Init()
// MVD // TODO: to be exchanged with specific flags (timeInfo, fieldInfo etc.) (S.Zh.)
stationInfo.SetTimeInfo(1);
stationInfo.SetZ(layer->GetZ());
stationInfo.SetMaterial(layer->GetDz(), 0); // TODO: Why rad len is 0????? (S.Zh.)
stationInfo.SetMaterial(layer->GetDz(), 10); // TODO: Why rad len is 0????? (S.Zh.)
stationInfo.SetXmax(100.);
stationInfo.SetYmax(100.);
stationInfo.SetRmin(0.);
......@@ -1074,9 +1085,10 @@ InitStatus CbmL1::Init()
algo->fRadThick[iSta].SetBins(NBins, RMax);
algo->fRadThick[iSta].table.resize(NBins);
float hole = 0.15;
for (int iB = 0; iB < NBins; iB++) {
algo->fRadThick[iSta].table[iB].resize(NBins);
float hole = 0.15;
for (int iB2 = 0; iB2 < NBins; iB2++) {
algo->fRadThick[iSta].table[iB][iB2] = 0.01 * hStaRadLen->GetBinContent(iB, iB2);
// Correction for holes in material map
......@@ -1137,10 +1149,10 @@ InitStatus CbmL1::Init()
const float RMax = hStaRadLen->GetXaxis()->GetXmax(); // should be same as min
algo->fRadThick[iSta].SetBins(NBins, RMax);
algo->fRadThick[iSta].table.resize(NBins);
float hole = 0.15;
for (int iB = 0; iB < NBins; iB++) {
algo->fRadThick[iSta].table[iB].resize(NBins);
float hole = 0.15;
for (int iB2 = 0; iB2 < NBins; iB2++) {
algo->fRadThick[iSta].table[iB][iB2] = 0.01 * hStaRadLen->GetBinContent(iB, iB2);
// Correction for holes in material map
......
reco/L1/CbmL1.h
+ 14
− 0
View file @ e1719e4c
......@@ -149,6 +149,14 @@ public:
/// Sets a vector of detectors used in tracking
void SetActiveTrackingDetectorIDs(const std::set<L1DetectorID>& init) { fActiveTrackingDetectorIDs = init; }
void SetUseMcHit(int StsUseMcHit = 0, int MuchUseMcHit = 0, int TrdUseMcHit = 0, int TofUseMcHit = 0)
{
fStsUseMcHit = StsUseMcHit;
fMuchUseMcHit = MuchUseMcHit;
fTrdUseMcHit = TrdUseMcHit;
fTofUseMcHit = TofUseMcHit;
}
void SetStsMaterialBudgetFileName(TString fileName) { fStsMatBudgetFileName = fileName; }
void SetMvdMaterialBudgetFileName(TString fileName) { fMvdMatBudgetFileName = fileName; }
......@@ -278,6 +286,11 @@ private:
Double_t fMomentumCutOff {0.1}; // currently not used
Bool_t fGhostSuppression {true}; // currently not used
Int_t fStsUseMcHit {-1}; // if STS data should be processed
Int_t fMuchUseMcHit {-1}; // if Much data should be processed
Int_t fTrdUseMcHit {-1}; // if Trd data should be processed
Int_t fTofUseMcHit {-1}; // if Tof data should be processed
Bool_t fUseMVD {false}; // if Mvd data should be processed
Bool_t fUseMUCH {false}; // if Much data should be processed
Bool_t fUseTRD {false}; // if Trd data should be processed
......@@ -326,6 +339,7 @@ private:
TClonesArray* fTofHits {nullptr}; // CbmMatches array
CbmTofDigiPar* fDigiPar {nullptr};
CbmTofDigiBdfPar* fTofDigiBdfPar {nullptr};
vector<vector<int>> TofPointToTrack;
TFile* fPerfFile {nullptr};
TDirectory* fHistoDir {nullptr};
......
reco/L1/CbmL1MCPoint.h
+ 1
− 0
View file @ e1719e4c
......@@ -66,6 +66,7 @@ struct CbmL1MCPoint {
int pointId = -1;
int file = -1;
int event = -1;
int trackId = -1;
L1Vector<int> hitIds {"CbmL1MCPoint::hitIds"}; // indices of CbmL1Hits in L1->vStsHits array
};
......
reco/L1/CbmL1ReadEvent.cxx
+ 221
− 166
View file @ e1719e4c
This diff is collapsed.
reco/L1/L1Algo/L1Algo.cxx
+ 2
− 0
View file @ e1719e4c
......@@ -132,6 +132,7 @@ void L1Algo::Init(const L1Vector<fscal>& geo, const bool UseHitErrors, const Tra
double f_sigma = geo[ind++];
double b_phi = geo[ind++];
double b_sigma = geo[ind++];
double dt = geo[ind++];
double c_f = cos(f_phi);
double s_f = sin(f_phi);
double c_b = cos(b_phi);
......@@ -144,6 +145,7 @@ void L1Algo::Init(const L1Vector<fscal>& geo, const bool UseHitErrors, const Tra
st.backInfo.cos_phi = c_b;
st.backInfo.sin_phi = s_b;
st.backInfo.sigma2 = b_sigma * b_sigma;
st.dt = dt;
if (fabs(b_phi - f_phi) < .1) {
double th = b_phi - f_phi;
......
reco/L1/L1Algo/L1Station.h
+ 1
− 1
View file @ e1719e4c
......@@ -12,7 +12,6 @@
class L1Station {
public:
// TODO: test this constructors
//L1Station(const L1Station &) = default;
//L1Station & operator=(const L1Station &) = default;
......@@ -25,6 +24,7 @@ public:
fvec Rmin {0};
fvec Rmax {0};
fvec Sy {0};
fvec dt {0};
L1MaterialInfo materialInfo {};
L1FieldSlice fieldSlice {};
L1UMeasurementInfo frontInfo {};
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment