Dear colleagues,
I have observed a strangeness with number of entries in TTree for simulation, digitization, reconstruction and pid, i.e. these number are different and I expect them to be equal.
I do simulation with eta_c macros from /macro/run/tdrct/eta_c.
With 2000 events simulation with trunk rev.12565
Running simple macro on produced data
void test_data()
{
TString inSimFileName = "evt_points_stt.root";
TString inDigiFileName = "evt_digi_stt.root";
TString inRecoFileName = "evt_reco_stt.root";
TString inPidFileName = "evt_pid_stt.root";
TFile *inSimFile = TFile::Open(inSimFileName,"READ");
TTree *treeSim=(TTree *) inSimFile->Get("cbmsim");
TFile *inDigiFile = TFile::Open(inDigiFileName,"READ");
TTree *treeDigi=(TTree *) inDigiFile->Get("cbmsim");
TFile *inRecoFile = TFile::Open(inRecoFileName,"READ");
TTree *treeReco=(TTree *) inRecoFile->Get("cbmsim");
TFile *inPidFile = TFile::Open(inPidFileName,"READ");
TTree *treePid=(TTree *) inPidFile->Get("cbmsim");
std::cout<<"treeSim->GetEntriesFast()="<<treeSim->GetEntriesFast()<<std::endl;
std::cout<<"treeDigi->GetEntriesFast()="<<treeDigi->GetEntriesFast()<<std::endl;
std::cout<<"treeReco->GetEntriesFast()="<<treeReco->GetEntriesFast()<<std::endl;
std::cout<<"treePid->GetEntriesFast()="<<treePid->GetEntriesFast()<<std::endl;
}
I have the following output
treeSim->GetEntriesFast()=2000
treeDigi->GetEntriesFast()=2001
treeReco->GetEntriesFast()=2002
treePid->GetEntriesFast()=2001
And those number are the same for stt and tpc simulation.
Should those numbers be equal?
Dima
GSI Forum
Search
Help
Members
Login
Home
