Processing run_reco.C... PSaid instance created... access via gSaid->f() - RTDB container factory CbmBaseContFact - RTDB container factory CbmFieldContFact - RTDB container factory CbmPassiveContFact - RTDB container factory CbmMvdContFact - RTDB container factory CbmStsContFact - RTDB container factory CbmRichContFact - RTDB container factory CbmTrdContFact - RTDB container factory CbmTofContFact - RTDB container factory CbmEcalContFact -I- CbmRunAna: Opening Input file: data/test.mc.root -I- Map Correction input file: /u/alla/cbmroot/DEC08/parameters/rich/muon_radius_correction_map.root -I- CbmRunAna::Init : data/test.mc.root is connected with: create CbmFieldPar container CbmFieldPar -I container name CbmGeoStsPar -I container name CbmStsDigiPar -I container name CbmMvdGeoPar -I container name CbmGeoTofPar * HitProducer * :: SetParContainers() ************************************************************* initialisation for run id 1813359249 ************************************************************* Container CbmBaseParSet initialized from ROOT file. Container CbmFieldPar initialized from ROOT file. Container CbmGeoStsPar initialized from ROOT file. -I- Ascii Io init() CbmStsDigiPar CbmStsDigiPar initialized from Ascii file Container CbmGeoPassivePar initialized from ROOT file. Container CbmMvdGeoPar initialized from ROOT file. Container CbmGeoTrdPar initialized from ROOT file. Container CbmGeoRichPar initialized from ROOT file. Container CbmGeoTofPar initialized from ROOT file. New field at 0xb452330, type 3 -I- CbmFieldMap: Reading field map from ROOT file /u/alla/cbmroot/DEC08/input/FieldMuonMagnet.root ====================================================== ---- CBM Magnetic Field : FieldMuonMagnet ---- ---- Field type : Map sym3 ---- ---- Field map grid : ---- x = 0.000 to 150.0 cm, 151 grid points, dx = 1.000 cm ---- y = 0.000 to 150.0 cm, 151 grid points, dy = 1.000 cm ---- z = 0.000 to 250.0 cm, 251 grid points, dz = 1.000 cm ---- Field centre position: ( 0.000, 0.000, 50.00) cm ---- Field scaling factor: 1.000 ---- ---- Field at origin is ( 0.000, -7.272, -2.724e-09) kG ====================================================== -I- STSDigitize::Init: STS digitisation scheme succesfully initialised Stations: 8, Sectors: 930, Channels: 123573544 -I- STSFindHits::Init: STS digitisation scheme succesfully initialised Stations: 8, Sectors: 930, Channels: 123573544 -I- STSMatchHits::Init: STS digitisation scheme succesfully initialised Stations: 8, Sectors: 930, Channels: 123573544 KALMAN FILTER : === INIT MAGNETIC FIELD === Magnetic field done KALMAN FILTER : === READ MVD MATERIAL === Mvd detector tube ( id, x,y,z, dz, r, R, RadL )= ( 1101, 0.000000, 0.000000, 10.000000, 0.015000, 0.550000, 5.000000, 9.349529 ) Mvd detector tube ( id, x,y,z, dz, r, R, RadL )= ( 1102, 0.000000, 0.000000, 20.000000, 0.015000, 1.050000, 10.000000, 9.349529 ) KALMAN FILTER : === READ STS MATERIAL === Sts material ( id, z, dz, r, R, RadL )= ( 1001, 30.00, 0.07500, 0.000, 100.0, 9.350 ) Sts material ( id, z, dz, r, R, RadL )= ( 1002, 35.00, 0.07500, 0.000, 100.0, 9.350 ) Sts material ( id, z, dz, r, R, RadL )= ( 1003, 40.00, 0.04000, 0.000, 100.0, 9.350 ) Sts material ( id, z, dz, r, R, RadL )= ( 1004, 50.00, 0.04000, 0.000, 100.0, 9.350 ) Sts material ( id, z, dz, r, R, RadL )= ( 1005, 60.00, 0.04000, 0.000, 100.0, 9.350 ) Sts material ( id, z, dz, r, R, RadL )= ( 1006, 75.00, 0.04000, 0.000, 100.0, 9.350 ) Sts material ( id, z, dz, r, R, RadL )= ( 1007, 95.00, 0.04000, 0.000, 100.0, 9.350 ) Sts material ( id, z, dz, r, R, RadL )= ( 1008, 100.0, 0.04000, 0.000, 100.0, 9.350 ) KALMAN FILTER : === READ RICH MATERIAL === Rich material tube ( id, x,y,z, dz, r, R, RadL )= ( 0, 0.000000, 0.000000, 160.012500, 0.025000, 0.000000, 1000.000000, 28.409874 ) Rich material tube ( id, x,y,z, dz, r, R, RadL )= ( 0, 0.000000, 0.000000, 245.012500, 169.975000, 0.000000, 1000.000000, 18234.587393 ) Rich material tube ( id, x,y,z, dz, r, R, RadL )= ( 1005, 0.000000, 92.974000, 330.300000, 0.600000, 0.000000, 1000.000000, 12.263924 ) Rich material tube ( id, x,y,z, dz, r, R, RadL )= ( 0, 0.000000, 0.000000, 340.287500, 19.375000, 0.000000, 1000.000000, 18234.587393 ) Rich material tube ( id, x,y,z, dz, r, R, RadL )= ( 0, 0.000000, 0.000000, 349.987500, 0.025000, 0.000000, 1000.000000, 28.409874 ) KALMAN FILTER : === READ TRD PASSIVE MATERIAL === printout skipped for Verbose mode 1 KALMAN FILTER : === READ TRD DETECTORS === Trd detector trd1mod1gas: ( id, z, dz, RadL )= ( 1014, 503.300000, 0.600000, 1792.871610 ), station 0 Trd detector trd1mod2gas: ( id, z, dz, RadL )= ( 1022, 503.300000, 0.600000, 1792.871610 ), station 0 Trd detector trd1mod3gas: ( id, z, dz, RadL )= ( 1030, 503.300000, 0.600000, 1792.871610 ), station 0 Trd detector trd2mod1gas: ( id, z, dz, RadL )= ( 1040, 728.300000, 0.600000, 1792.871610 ), station 1 Trd detector trd2mod2gas: ( id, z, dz, RadL )= ( 1048, 728.300000, 0.600000, 1792.871610 ), station 1 Trd detector trd2mod3gas: ( id, z, dz, RadL )= ( 1056, 728.300000, 0.600000, 1792.871610 ), station 1 Trd detector trd3mod3gas: ( id, z, dz, RadL )= ( 1066, 953.300000, 0.600000, 1792.871610 ), station 2 KALMAN FILTER : === READ TOF PASSIVE MATERIAL === printout skipped for Verbose mode 1 KALMAN FILTER : === READ TOF DETECTORS === KALMAN FILTER : === READ BEAM PIPE MATERIAL === Pipe material ( id, {z1, r1, R1}, {z2, r2, R2}, RL )= ( 0, {-5.000, 2.500, 2.550}, {2.500, 2.500, 2.550}, 18.76 ) Pipe material ( id, {z1, r1, R1}, {z2, r2, R2}, RL )= ( 0, {2.500, 2.500, 2.550}, {3.500, 13.00, 13.05}, 18.76 ) Pipe material ( id, {z1, r1, R1}, {z2, r2, R2}, RL )= ( 0, {3.500, 13.00, 13.05}, {24.00, 13.00, 13.05}, 18.76 ) Pipe material ( id, {z1, r1, R1}, {z2, r2, R2}, RL )= ( 0, {24.00, 13.00, 13.05}, {27.00, 1.000, 1.050}, 18.76 ) Pipe material ( id, {z1, r1, R1}, {z2, r2, R2}, RL )= ( 0, {27.00, 1.000, 1.050}, {160.0, 3.200, 3.250}, 18.76 ) KALMAN FILTER : === READ TARGET MATERIAL === Target material tube ( id, x,y,z, dz, r, R, RadL )= ( 0, 0.000000, 0.000000, 0.000000, 0.025000, 0.000000, 0.250000, 0.335086 ) L1Algo Input Magnetic field:0.000 0.000 -7.272 -2.724e-09 L1Algo Input Magnetic field:2.500 6.327e-30 -7.595 -2.722e-09 L1Algo Input Magnetic field:5.000 0.000 -7.917 -2.709e-09 L1Algo Input 10 Stations: 10.00 0.01500 9.350, 10 field coeff. 0.000 0.0005000 1.571 0.0005000 20.00 0.01500 9.350, 10 field coeff. 0.000 0.0005000 1.571 0.0005000 30.00 0.07500 9.350, 10 field coeff. 0.000 0.001443 1.571 0.001443 35.00 0.07500 9.350, 10 field coeff. 0.000 0.001443 1.571 0.001443 40.00 0.04000 9.350, 10 field coeff. 0.000 0.001732 0.2618 0.001673 50.00 0.04000 9.350, 10 field coeff. 0.000 0.001732 0.2618 0.001673 60.00 0.04000 9.350, 10 field coeff. 0.000 0.001732 0.2618 0.001673 75.00 0.04000 9.350, 10 field coeff. 0.000 0.001732 0.2618 0.001673 95.00 0.04000 9.350, 10 field coeff. 0.000 0.001732 0.2618 0.001673 100.0 0.04000 9.350, 10 field coeff. 0.000 0.001732 0.2618 0.001673 L1Algo initialized -I- STSFindTracks::Init: STS digitisation scheme succesfully initialised Stations: 8, Sectors: 930, Channels: 123573544 -I- STSFindTracks intialised -I- Track finder engine selected -I- No MVD hits available! -E- CbmRootManager Branch: MVDHit not found in Tree -W- CbmStsKFTrackFitter::Init: no MVD hits array * HitProducer * :: Init() ----- CbmGeoTrdPar ----- Context/Purpose: TestDefaultContext first input version: 1 second input version: -1 has changed is not static -I version is with 12 layers (3 stations x 4 layers each) and new segmented geometry *********** Initilization of Trd Radiator ************** Nr. of foils : 130 Foil thickness : 0.001300 cm Gap thickness : 0.02000 cm Simple TR production: 1 ************* End of Trd Radiator Init ************** ********** Initilization of TRD Gas ************ -I- Detector type : standard GSI geometry (2) -I- Gas thickness : 0.60 cm -I- Percent (Xe) : 85. -I- Percent (CO2) : 15. ********** End of TRD Gas Init ************ -I- CbmL1TrdTrackFinderSts::TrdLayout : TRD geometry : 3x4 standard. -E- CbmTrdSetTracksPidANN::Init: wrong geometry type. -I- CbmTrdSetTracksPidANN::Init : get NeuralNet weight parameters from: /u/alla/cbmroot/DEC08/parameters/trd/Neural_Net_Weights_El_ID_12.txt , /u/alla/cbmroot/DEC08/parameters/trd/Neural_Net_Weights_El_ID_11.txt , /u/alla/cbmroot/DEC08/parameters/trd/Neural_Net_Weights_El_ID_10.txt , /u/alla/cbmroot/DEC08/parameters/trd/Neural_Net_Weights_El_ID_9.txt , /u/alla/cbmroot/DEC08/parameters/trd/Neural_Net_Weights_El_ID_8.txt , /u/alla/cbmroot/DEC08/parameters/trd/Neural_Net_Weights_El_ID_7.txt , /u/alla/cbmroot/DEC08/parameters/trd/Neural_Net_Weights_El_ID_6.txt , -I- CbmTofHitProducer: Intialization successfull -W- CbmFindGlobalTracks::Init: No RICH merger selected! -E- CbmRootManager Branch: RICHRing not found in Tree -W- CbmStsFindTracks::Init: No RichRing array! -I- CbmTrdSetTracksPidLike::ReadData -I- InputFile: /u/alla/cbmroot/DEC08/parameters/trd/Likelihood_Xenon_85_GSI_6.root Array detector 50. Array detector 30. Array detector 0.025 Global TrackArray used: extrapolation to z-plane at z = 3.0e+02 cm -------------------------------------------------------------------------------- CbmRichRingFitter::Init() -I- Map Correction input file: /u/alla/cbmroot/DEC08/parameters/rich/muon_radius_correction_map.root CbmRichRingFitter::Init() CbmRichRingSelect::Init() -I- CbmRichRingSelectNeuralNet: get NeuralNet weight parameters from: /u/alla/cbmroot/DEC08/parameters/rich/NeuralNet_RingSelection_Weights.txt CbmRichRingFitter::Init() CbmRichRingSelect::Init() -I- CbmRichRingSelectNeuralNet: get NeuralNet weight parameters from: /u/alla/cbmroot/DEC08/parameters/rich/NeuralNet_RingSelection_Weights.txt Starting run No primary vertex !!! L1 ACCUMULATED STAT : 1 EVENTS RefPrim efficiency : 0.000 / 0.000 | 0 RefSec efficiency : 0.000 / 0.000 | 0 Refset efficiency : 0.000 / 0.000 | 0 Allset efficiency : 0.000 / 0.000 | 0 ExtraPrim efficiency : 0.000 / 0.000 | 0 ExtraSec efficiency : 0.000 / 0.000 | 0 Extra efficiency : 0.000 / 0.000 | 0 Clone probability : 0.000 | 0 Ghost probability : 0.000 | 0 MC tracks/event found : 0 CA Track Finder: 0.000 s/ev + STSFitTracks : 0.0000 s, tracks fitted 0 * No. of Trds per Station = 4 ** 13438 Trd hits to be created in this event** -I- CbmL1TrdTrackFinderSts::SortTrdHits : TRD plane no. 0 has 887 hits TRD plane no. 1 has 900 hits TRD plane no. 2 has 932 hits TRD plane no. 3 has 983 hits TRD plane no. 4 has 1100 hits TRD plane no. 5 has 1129 hits TRD plane no. 6 has 1184 hits TRD plane no. 7 has 1231 hits TRD plane no. 8 has 1232 hits TRD plane no. 9 has 1260 hits TRD plane no. 10 has 1277 hits TRD plane no. 11 has 1323 hits -I- CbmL1TrdTrackFinderSts::DoFind : 0 TRD tracks found. -I- CbmL1TrdTrackFinder : 1 events processed -I- CbmTrdMatchTracks: rec. 0, quota 0.0000 % -I- CbmTofHitProducer : 14187 points in Tof for this event -I- CbmTofHitProducer : 2363 tracks in Tof -I- CbmTofHitProducer : 434 tracks in Tof from vertex -I- CbmTofHitProducer : 960 tracks in Tof produced locally (Z>996 cm) -I- CbmTofHitProducer : 1403 tracks in Tof able to produce a hit -I- CbmTofHitProducer : 1656 hits in Tof created -I- CbmFindGlobalTracks: 0 global tracks, 0 merged, 0 with RICH ring, 0 with TOF hit ------------------------------------------------------- -I- Primary Vertex Finder -I- Vertex coord. (0.0000,0.0000,0.0000) cm, chi2/ndf = -0.0000, 0 tracks used ------------------------------------------------------- 0 global tracks -------------------------------------------------------------------------------------- -I- Track Extrapolation -I- Global tracks: 0 Points written to Point Array for Projection to RICH photodetector plane: 0 -------------------------------------------------------------------------------------- =====> Event no. 1 Number of found ring = 0 TIME CPU time 0.0300 s -I- CbmRichMatchRings: rec. 0, quota 0.0000 % -I- CbmRichAssignTrack: call concrete assign class =====> Event no. 1 Number of input MC points: 4364 ============================================================ ===== STSMatchHits: Run summary ===== ===== Good events : 1 ===== Failed events : 0 ===== Average time : 0.0000 s ===== ===== Hits per event : 0. ===== Matched hits : nan % ===== Distant hits : nan % ===== Background hits : nan % ============================================================ ============================================================ ===== STSFindTracks: Run summary ===== ===== Good events : 1 ===== Failed events : 0 ===== Average time : 0.0000 s ===== ===== Found tracks per event : 0. ============================================================ ============================================================ ===== STSMatchTracks: Run summary ===== ===== Good events : 1 ===== Failed events : 0 ===== Average time : 0.0000 s ===== ===== Tracks per event : 0. ===== True hits : nan % ============================================================ ============================================================ ===== STSFitTracks: Run summary ===== ===== Good events : 1 ===== Failed events : 0 ===== Average time : 0.0000 s ===== ===== Fitted tracks per event : 0. ============================================================ Macro finished succesfully. Output file is data/test.eds.root Parameter file is data/params.root Real time 21. s, CPU time 16. s Test passed All ok