//--------------------------------------------------------- bool CCMD_Module::_Load_Input(CSG_Parameter *pParameter) { wxString FileName; if( !pParameter->is_Input() || !pParameter->is_Enabled() ) { return( true ); } if( !m_CMD.Found(_Get_ID(pParameter), &FileName) ) { return( pParameter->is_Optional() ); } if( pParameter->is_DataObject() ) { if( !SG_Get_Data_Manager().Find(&FileName) && !SG_Get_Data_Manager().Add(&FileName) && !pParameter->is_Optional() ) { CMD_Print_Error(_TL("input file"), &FileName); return( false ); } return( pParameter->Set_Value(SG_Get_Data_Manager().Find(&FileName)) ); } else if( pParameter->is_DataObject_List() ) { pParameter->asList()->Del_Items(); wxString FileNames(FileName); do { FileName = FileNames.BeforeFirst (';'); FileNames = FileNames.AfterFirst (';'); if( !SG_Get_Data_Manager().Find(&FileName) ) { SG_Get_Data_Manager().Add(&FileName); } pParameter->asList()->Add_Item(SG_Get_Data_Manager().Find(&FileName)); } while( FileNames.Length() > 0 ); } return( true ); }
// Main function int POTCheck(std::string GeneratorName, unsigned int ThreadNumber, unsigned int NumberOfThreads) { std::string Version = "v05_08_00"; std::string FileNumberStr; if(NumberOfThreads > 1) { FileNumberStr = "_" + std::to_string(ThreadNumber); } else { FileNumberStr = ""; } std::cout << "Data Sample : " << GeneratorName << std::endl; TChain *treenc = new TChain("analysistree/pottree"); if(GeneratorName == "data_onbeam_bnb") { // treenc -> Add( ("/pnfs/uboone/persistent/users/aschu/onbeam_data_bnbSWtrigger/"+GeneratorName+"_"+Version+"_anatree.root").c_str() ); std::ifstream FileNames("/pnfs/uboone/persistent/users/aschu/devel/v05_11_01/hadd/GOODBNBBEAM/filesana.list"); std::string FileName; while(std::getline(FileNames,FileName)) { std::cout << FileName << std::endl; treenc -> Add((FileName).c_str()); } } else if(GeneratorName == "data_offbeam_bnbext") { // treenc -> Add( ("/pnfs/uboone/persistent/users/aschu/offbeam_data_bnbSWtrigger/"+GeneratorName+"_"+Version+"_anatree.root").c_str() ); std::ifstream FileNames("/pnfs/uboone/persistent/users/aschu/devel/v05_11_01/hadd/GOODEXTBNB/filesana.list"); std::string FileName; while(std::getline(FileNames,FileName)) { std::cout << FileName << std::endl; treenc -> Add((FileName).c_str()); } } else if(GeneratorName == "prodgenie_bnb_nu_cosmic_uboone_field") { // treenc -> Add( ("/pnfs/uboone/persistent/users/aschu/offbeam_data_bnbSWtrigger/"+GeneratorName+"_"+Version+"_anatree.root").c_str() ); std::ifstream FileNames("/pnfs/uboone/persistent/users/sowjanya/v05_08_00/bnbpluscosmics_nominal_scOn_100k/ana/filesana.list"); std::string FileName; while(std::getline(FileNames,FileName)) { std::cout << FileName << std::endl; treenc -> Add((FileName).c_str()); } } else if(GeneratorName == "prodgenie_bnb_nu_cosmic_uboone") { treenc -> Add( ("/pnfs/uboone/persistent/users/aschu/MC_BNB_Cosmic/"+GeneratorName+"_"+Version+"_anatree.root").c_str() ); } else if(GeneratorName == "TEM" || GeneratorName == "MEC" || GeneratorName == "MA" || GeneratorName == "BITE") { treenc -> Add( ("/pnfs/uboone/persistent/users/aschu/"+GeneratorName+"/"+GeneratorName+"merge.root").c_str() ); } else if(GeneratorName == "test") { for(unsigned int index = 1; index < 8; index++) { treenc -> Add( ("/uboone/data/users/aschu/FilterTestFiles/ana_hist_"+std::to_string(index)+".root").c_str() ); } } else { treenc -> Add( ("/lheppc46/data/uBData/anatrees/"+GeneratorName+"_"+Version+"_anatree.root").c_str() ); } // treenc -> Add( ("/media/christoph/200EFBDA63AA160B/anatrees/"+GeneratorName+"_"+Version+"_anatree.root").c_str() ); Double_t pot; Double_t potbnbETOR860; Double_t potbnbETOR875; Double_t potnumiETORTGT; //define cut variables double TotPOT = 0; //cm. Distance flash-track double TotPOTBNB860 = 0; //cm. Distance track start/end to vertex double TotPOTBNB875 = 0; //cm. Length of longest track double TotPOTBNBTGT = 0; unsigned long int Size = treenc -> GetEntries(); // Event Loop for(int i = 0; i < Size; i++) { // Asign Branches treenc -> SetBranchAddress("pot", &pot); treenc -> SetBranchAddress("potbnbETOR860", &potbnbETOR860); treenc -> SetBranchAddress("potbnbETOR875", &potbnbETOR875); treenc -> SetBranchAddress("potnumiETORTGT", &potnumiETORTGT); if(pot > .1) std::cout << i << " " << pot << std::endl; // Calculate total POTs TotPOT += pot; TotPOTBNB860 += potbnbETOR860; TotPOTBNB875 += potbnbETOR875; TotPOTBNBTGT += potnumiETORTGT; }//loop over all events std::cout << "TotPOT " << TotPOT << std::endl; std::cout << "TotPOTBNB860 " << TotPOTBNB860 << std::endl; std::cout << "TotPOTBNB875 " << TotPOTBNB875 << std::endl; std::cout << "TotPOTBNBTGT " << TotPOTBNBTGT << std::endl; return 0; } // end main function