From cac041c03dbf3d3044379b2691578d8810421ad0 Mon Sep 17 00:00:00 2001 From: Kamil Sedlak Date: Fri, 4 Nov 2011 13:26:02 +0000 Subject: [PATCH] 4.11.2011 - Kamil Sedlak Increased maximum size of the simulation output file to 100 GB. --- musrSimAna/musrCounter.cxx | 27 +++++++++++++++------------ src/musrRootOutput.cc | 2 ++ 2 files changed, 17 insertions(+), 12 deletions(-) diff --git a/musrSimAna/musrCounter.cxx b/musrSimAna/musrCounter.cxx index ee64acf..25a60bb 100644 --- a/musrSimAna/musrCounter.cxx +++ b/musrSimAna/musrCounter.cxx @@ -174,13 +174,14 @@ Bool_t musrCounter::GetNextGoodMuon(Int_t evtID, Long64_t timeBinMin, Long64_t& Long64_t timeBinOfCount_tmp = it->first; timeBinOfNextHit = timeBinOfCount_tmp; + Int_t modulo = timeBinOfNextHit % 524288; if (timeBinOfCount_tmp <= timeBinMin) continue; // This hit was already processed previously ==> skip it Int_t eventNumber = (it->second)->eventIDnumber; if (eventNumber!=evtID) continue; // This trigger hit does not correspond to the currently processed event // ==> skip it, because it was already proceesed or will be processed in future numberOfMuonCandidates++; - if (bool_debugingRequired) {if (debugEventMap[evtID]>3) std::cout<<"GetNextGoodMuon: muon candidate found ("<3) std::cout<<"GetNextGoodMuon: muon candidate found ("<3) std::cout<<"GetNextGoodMuon: muon candidate not in koincidence ("<3) std::cout<<"GetNextGoodMuon: muon candidate not in koincidence ("<3) std::cout<<"GetNextGoodMuon: muon candidate vetoed ("<3) std::cout<<"GetNextGoodMuon: muon candidate vetoed ("<3) std::cout<<"GetNextGoodMuon: muon candidate after VK ("<3) std::cout<<"GetNextGoodMuon: muon candidate after VK ("<3) std::cout<<"GetNextGoodMuon: muon candidate killed by pileup muon ("<3) std::cout<<"GetNextGoodMuon: muon candidate killed by pileup muon ("<second)->eventEntry; @@ -218,7 +219,7 @@ Bool_t musrCounter::GetNextGoodMuon(Int_t evtID, Long64_t timeBinMin, Long64_t& idetID = (it->second)->det_id; idetEdep = (it->second)->det_edep; if (bool_debugingRequired) {if (debugEventMap[evtID]>1) std::cout<<"GetNextGoodMuon: GOOD muon candidate found (" - <::iterator itt = it_muon_hits_to_be_deleted.begin(); itt != it_muon_hits_to_be_deleted.end(); ++itt) { hitMap.erase(*itt); } @@ -278,6 +279,7 @@ Int_t musrCounter::GetNextGoodPositron(Int_t evtID, Long64_t timeBinMin, Long64_ for (hitMap_TYPE::iterator it = hitMap.begin(); it != hitMap.end(); ++it) { // Int_t eventNumber = (it->second)->eventIDnumber; Long64_t timeBinOfCount_tmp = it->first; + Int_t modulo = timeBinOfCount_tmp % 524288; if ((timeBinOfCount_tmp <= timeBinMin) || (timeBinOfCount_tmp > timeBinMax)) { if (bool_debugingRequired) { if (debugEventMap[evtID]>3) {std::cout<<"GetNextGoodPositron: Hit out of data interval"<second)->IsInCoincidence(timeBinOfCount_tmp,'P') )) { if (bool_debugingRequired) { - if (debugEventMap[evtID]>3) {std::cout<<"GetNextGoodPositron: Coincidence required but not found (timeBin="<3) {std::cout<<"GetNextGoodPositron: Coincidence required but not found (timeBin="< remove the candidate. if (bool_ignoreUnperfectPositrons) it_positron_hits_to_be_deleted.push_back(it); @@ -305,7 +307,7 @@ Int_t musrCounter::GetNextGoodPositron(Int_t evtID, Long64_t timeBinMin, Long64_ for (counterMapType::const_iterator itCounter = vetoCounterMap.begin(); itCounter!=vetoCounterMap.end(); ++itCounter) { if ( (itCounter->second)->IsInCoincidence(timeBinOfCount_tmp,'P') ) { if (bool_debugingRequired) { - if (debugEventMap[evtID]>3) {std::cout<<"GetNextGoodPositron: Coincidence vith veto detector found (timeBin="<3) {std::cout<<"GetNextGoodPositron: Coincidence vith veto detector found (timeBin="< remove the candidate. if (bool_ignoreUnperfectPositrons) it_positron_hits_to_be_deleted.push_back(it); @@ -332,7 +334,7 @@ Int_t musrCounter::GetNextGoodPositron(Int_t evtID, Long64_t timeBinMin, Long64_ timeBinOfNextGoodHit = timeBinOfCount_tmp; timeBinOfNextGoodHit_phaseShifted = (it->second) -> timeBin2; if (bool_debugingRequired) { - if (debugEventMap[evtID]>3) {std::cout<<"GetNextGoodPositron: Good positron candidate found in this counter. (timeBin="<second)->eventIDnumber<<")"<3) {std::cout<<"GetNextGoodPositron: Good positron candidate found in this counter. (timeBin="<second)->eventIDnumber<<")"<::iterator itt = it_positron_hits_to_be_deleted.begin(); itt != it_positron_hits_to_be_deleted.end(); ++itt) { @@ -376,10 +378,11 @@ void musrCounter::SetCoincidenceTimeWindowOfAllVetoDetectors(Long64_t maxCoinc, void musrCounter::myPrintThisCounter(Int_t evtID, Int_t detail) { Bool_t eventMixing=false; if ((hitMap.begin()==hitMap.end()) && (detail<=1) ) return; - if (detail>1) std::cout<<"musrCounter::myPrintThisCounter: counterNr = "<1) std::cout<<"musrCounter::myPrintThisCounter: counterNr = "<first<<" "<<(it->second)->eventIDnumber<<","; + std::cout<<"\t"<first<<" ("<<(it->first)%524288 <<") "<<(it->second)->eventIDnumber<<","; if (evtID != (it->second)->eventIDnumber) {eventMixing=true;} } if (eventMixing) {std::cout<<" Potential event mixing";} diff --git a/src/musrRootOutput.cc b/src/musrRootOutput.cc index 76e2431..446ad10 100644 --- a/src/musrRootOutput.cc +++ b/src/musrRootOutput.cc @@ -30,6 +30,8 @@ //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... musrRootOutput::musrRootOutput() { + TTree::SetMaxTreeSize(100000000000LL); // Set maximum size of the tree file + // to 100 GB (instead of 1.9 GB). pointerToRoot=this; boolIsAnySpecialSaveVolumeDefined=false; nFieldNomVal=0;