Skip to content
Snippets Groups Projects
Commit 13ca9fd4 authored by Alberica Toia's avatar Alberica Toia
Browse files

fixed clang format

parent 542b4ed9
No related branches found
No related tags found
1 merge request!153adding sts reco with abs and rel path
...@@ -24,10 +24,11 @@ void build_event_win_kronos(UInt_t uRunIdx = 0, ...@@ -24,10 +24,11 @@ void build_event_win_kronos(UInt_t uRunIdx = 0,
*/ */
}; };
if (vuListRunId.size() <= uRunIdx) return kFALSE; if (vuListRunId.size() <= uRunIdx) return kFALSE;
uRunId = vuListRunId[uRunIdx]; uRunId = vuListRunId[uRunIdx];
//fileName = Form("data/unp_mcbm_%03u.root", uRunId); //fileName = Form("data/unp_mcbm_%03u.root", uRunId);
fileName = Form("/lustre/cbm/users/ploizeau/mcbm2020/unp_evt_data_7f229b3f_20201103/unp_mcbm_%i.root", fileName = Form("/lustre/cbm/users/ploizeau/mcbm2020/"
uRunId); "unp_evt_data_7f229b3f_20201103/unp_mcbm_%i.root",
uRunId);
} // if( 99999 != uRunIdx ) } // if( 99999 != uRunIdx )
if (uRunId < 692 && 0 != uRunId) return kFALSE; if (uRunId < 692 && 0 != uRunId) return kFALSE;
...@@ -105,18 +106,20 @@ void build_event_win_kronos(UInt_t uRunIdx = 0, ...@@ -105,18 +106,20 @@ void build_event_win_kronos(UInt_t uRunIdx = 0,
fRun->AddTask(eventBuilder); fRun->AddTask(eventBuilder);
// ----- Parameter database -------------------------------------------- // ----- Parameter database --------------------------------------------
TString parFile = TString parFile =
Form("/lustre/cbm/users/ploizeau/mcbm2020/unp_evt_data_7f229b3f_20201103/unp_mcbm_params_%i.root", Form("/lustre/cbm/users/ploizeau/mcbm2020/unp_evt_data_7f229b3f_20201103/"
uRunId); "unp_mcbm_params_%i.root",
FairRuntimeDb* rtdb = fRun->GetRuntimeDb(); uRunId);
FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
FairParRootFileIo* parIo1 = new FairParRootFileIo(); FairParRootFileIo* parIo1 = new FairParRootFileIo();
parIo1->open(parFile.Data(),"UPDATE"); parIo1->open(parFile.Data(), "UPDATE");
rtdb->setFirstInput(parIo1); rtdb->setFirstInput(parIo1);
// ------------------------------------------------------------------------ // ------------------------------------------------------------------------
TString geoFileSts = TString geoFileSts =
"/lustre/cbm/users/alberica/cbmroot/macro/beamtime/mcbm2020/data/test.geo.root"; // to be created by a simulation run "/lustre/cbm/users/alberica/cbmroot/macro/beamtime/mcbm2020/data/"
"test.geo.root"; // to be created by a simulation run
fRun->SetGeomFile(geoFileSts); fRun->SetGeomFile(geoFileSts);
// ----- Local reconstruction in STS ---------------------------------- // ----- Local reconstruction in STS ----------------------------------
...@@ -125,16 +128,17 @@ void build_event_win_kronos(UInt_t uRunIdx = 0, ...@@ -125,16 +128,17 @@ void build_event_win_kronos(UInt_t uRunIdx = 0,
//recoSts->SetTimeCutDigisAbs( 20 );// cluster finder: time cut in ns //recoSts->SetTimeCutDigisAbs( 20 );// cluster finder: time cut in ns
//recoSts->SetTimeCutClustersAbs(20.); // hit finder: time cut in ns //recoSts->SetTimeCutClustersAbs(20.); // hit finder: time cut in ns
// ASIC params: #ADC channels, dyn. range, threshold, time resol., dead time, // ASIC params: #ADC channels, dyn. range, threshold, time resol., dead time,
// noise RMS, zero-threshold crossing rate // noise RMS, zero-threshold crossing rate
auto parAsic = new CbmStsParAsic(32, 75000., 3000., 5., 800., 1000., 3.9789e-3); auto parAsic =
new CbmStsParAsic(32, 75000., 3000., 5., 800., 1000., 3.9789e-3);
// Module params: number of channels, number of channels per ASIC // Module params: number of channels, number of channels per ASIC
auto parMod = new CbmStsParModule(2048, 128); auto parMod = new CbmStsParModule(2048, 128);
parMod->SetAllAsics(*parAsic); parMod->SetAllAsics(*parAsic);
recoSts->UseModulePar(parMod); recoSts->UseModulePar(parMod);
// Sensor params // Sensor params
auto sensorPar = new CbmStsParSensor(CbmStsSensorClass::kDssdStereo); auto sensorPar = new CbmStsParSensor(CbmStsSensorClass::kDssdStereo);
sensorPar->SetPar(0, 6.2092); // Extension in x sensorPar->SetPar(0, 6.2092); // Extension in x
...@@ -143,18 +147,18 @@ void build_event_win_kronos(UInt_t uRunIdx = 0, ...@@ -143,18 +147,18 @@ void build_event_win_kronos(UInt_t uRunIdx = 0,
sensorPar->SetPar(3, 5.9692); // Active size in y sensorPar->SetPar(3, 5.9692); // Active size in y
sensorPar->SetPar(4, 1024.); // Number of strips front side sensorPar->SetPar(4, 1024.); // Number of strips front side
sensorPar->SetPar(5, 1024.); // Number of strips back side sensorPar->SetPar(5, 1024.); // Number of strips back side
sensorPar->SetPar(6, 0.0058); // Strip pitch front side sensorPar->SetPar(6, 0.0058); // Strip pitch front side
sensorPar->SetPar(7, 0.0058); // Strip pitch back side sensorPar->SetPar(7, 0.0058); // Strip pitch back side
sensorPar->SetPar(8, 7.5); // Stereo angle front side sensorPar->SetPar(8, 7.5); // Stereo angle front side
sensorPar->SetPar(9, 0.0); // Stereo angle back side sensorPar->SetPar(9, 0.0); // Stereo angle back side
recoSts->UseSensorPar(sensorPar); recoSts->UseSensorPar(sensorPar);
// Sensor conditions: full depletion voltage, bias voltage, temperature, // Sensor conditions: full depletion voltage, bias voltage, temperature,
// coupling capacitance, inter-strip capacitance // coupling capacitance, inter-strip capacitance
auto sensorCond = new CbmStsParSensorCond(70., 140., 268., 17.5, 1.); auto sensorCond = new CbmStsParSensorCond(70., 140., 268., 17.5, 1.);
recoSts->UseSensorCond(sensorCond); recoSts->UseSensorCond(sensorCond);
fRun->AddTask(recoSts); fRun->AddTask(recoSts);
std::cout << "-I- : Added task " << recoSts->GetName() << std::endl; std::cout << "-I- : Added task " << recoSts->GetName() << std::endl;
// ------------------------------------------------------------------------ // ------------------------------------------------------------------------
......
...@@ -21,15 +21,18 @@ void mcbm_reco(Int_t runId = 812, Int_t nTimeslices = 0) { ...@@ -21,15 +21,18 @@ void mcbm_reco(Int_t runId = 812, Int_t nTimeslices = 0) {
// ----- In- and output file names ------------------------------------ // ----- In- and output file names ------------------------------------
TString inFile = Form("/lustre/cbm/users/ploizeau/mcbm2020/" TString inFile = Form("/lustre/cbm/users/ploizeau/mcbm2020/"
"unp_evt_data_7f229b3f_20201103/unp_mcbm_%i.root", "unp_evt_data_7f229b3f_20201103/unp_mcbm_%i.root",
runId); runId);
TString parFile = TString parFile =
Form("/lustre/cbm/users/ploizeau/mcbm2020/unp_evt_data_7f229b3f_20201103/unp_mcbm_params_%i.root", Form("/lustre/cbm/users/ploizeau/mcbm2020/unp_evt_data_7f229b3f_20201103/"
runId); "unp_mcbm_params_%i.root",
runId);
TString geoFile = TString geoFile =
"/lustre/cbm/users/alberica/cbmroot/macro/beamtime/mcbm2020/data/test.geo.root"; // to be created by a simulation run "/lustre/cbm/users/alberica/cbmroot/macro/beamtime/mcbm2020/data/"
TString outFile = "test.geo.root"; // to be created by a simulation run
Form("/lustre/cbm/users/ploizeau/mcbm2020/unp_evt_data_7f229b3f_20201103/reco_mcbm_%i.root", runId); TString outFile = Form("/lustre/cbm/users/ploizeau/mcbm2020/"
"unp_evt_data_7f229b3f_20201103/reco_mcbm_%i.root",
runId);
// ------------------------------------------------------------------------ // ------------------------------------------------------------------------
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment