//Function for printing the main screen (everything on the main screne) void lifeTracker() { setup(); for(i = 0; i < 20; i++){ xt_par2(XT_SET_ROW_COL_POS, row = i+7, col = 1); printf("|"); xt_par2(XT_SET_ROW_COL_POS, row = i+7, col = 16); printf("|"); xt_par2(XT_SET_ROW_COL_POS, row = i+7, col = 19); printf("|"); xt_par2(XT_SET_ROW_COL_POS, row = i+7, col = 91); printf("|"); xt_par2(XT_SET_ROW_COL_POS, row = i+7, col = 94); printf("|"); xt_par2(XT_SET_ROW_COL_POS, row = i+7, col = 120); printf("|"); } xt_par2(XT_SET_ROW_COL_POS, row = 7, col = 2); fillCategory(); xt_par0(XT_CH_GREEN); printf("--Categories--"); xt_par0(XT_CH_WHITE); if(sort == 0) bubbleSortTime(); else if (sort == 1) bubbleSort(); xt_par2(XT_SET_ROW_COL_POS, row = 26, col = 2); xt_par0(XT_CH_BLUE); printf("---F6-Sort----"); updateRecords(0); xt_par2(XT_SET_ROW_COL_POS, row = 26, col = 20); xt_par0(XT_CH_BLUE); printf("-----------------------F2-Add F3-Delete F4-Edit------------------------"); xt_par0(XT_CH_GREEN); xt_par2(XT_SET_ROW_COL_POS, row = 7, col = 95); printf("---------Search----------"); clearSearch(); //also delete the code in between these lines xt_par0(XT_CH_BLUE); xt_par2(XT_SET_ROW_COL_POS, row = 26, col = 95); printf("--------F7-Search--------"); xt_par2(XT_SET_ROW_COL_POS, row = 28, col = 1); xt_par0(XT_CH_BLUE); printf("---------------------------------------------------------F9-EXIT--------------------------------------------------------"); xt_par0(XT_CH_WHITE); xt_par2(XT_SET_ROW_COL_POS,row=CATEGORY,col=2); }
int main(int argc, char **argv){ if( argc != 4 ){ printf( "Must specify a Database File, book order file, and a category names. \n Feel Free to refer to the readme for usage.\n" ); _exit( 1 ); // crash and burn } FILE *fp; if( (fp = fopen(argv[1], "r") ) == NULL ) { printf("fopen didn't work in %s line %d.\n", __FILE__ , __LINE__ ); printf("Please Enter Valid Database File \n"); _exit(1); // crash and burn } else if( (fp = fopen(argv[2], "r") ) == NULL ) { printf("fopen didn't work in %s line %d.\n", __FILE__ , __LINE__ ); printf("Please Enter Valid Book Order File \n"); _exit(1); // crash and burn } else if( (fp = fopen(argv[3], "r") ) == NULL ) { printf("fopen didn't work in %s line %d.\n", __FILE__ , __LINE__ ); printf("Please Enter Valid Category File \n"); _exit(1); // crash and burn } /* Done with error checking. */ fp = fopen(argv[1], "r"); int fSize; fseek(fp, 0L, SEEK_END); fSize = ftell(fp); fseek(fp, 0 , SEEK_SET); if(fSize == 0){ printf("The database file is empty, cannot procceed \n"); _exit(1); } customerList = readDatabase(fp); if( (error = fclose( fp ) ) > 0 ){ printf("fclose died in %s line %d: %s \n", __FILE__ , __LINE__ , strerror(error)); _exit(1); // crash and burn } totalProfit = 0; /* this gives us the full set of orders that we can then separate into the array of smaller queues based on the category. */ fp = fopen(argv[3], "r"); fseek(fp, 0L, SEEK_END); fSize = ftell(fp); fseek(fp , 0 , SEEK_SET); if(fSize == 0){ printf("The category file is empty, cannot procceed \n"); _exit(1); } numcats = readCategories(fp); categoryList = (category*)malloc(sizeof(category)*numcats); fillCategory(fp); if(( error = fclose( fp ) )> 0 ){ printf("fclose stuck it's nose where it didn't belong in %s line %d: %s\n", __FILE__ , __LINE__ , strerror( error ) ); _exit(1); // crash and burn } fp = fopen(argv[2], "r"); pthread_t tid; if ((error = pthread_create( &tid, NULL , readBookOrders , fp ))){ printf( "pthread_create() sparked the flames of revolution in %s line %d: %s\n", __FILE__ , __LINE__ , strerror( error ) ); _exit( 1 ); // crash and burn } createThreads(); /*create the consumer threads to wait for the producer, force a join on the producer after we create the consumers that way all threads are guaranteed to finish. */ if ((error = pthread_join(tid, NULL))) { printf( "pthread_join()ed the dark side in %s line %d: %s \n", __FILE__, __LINE__, strerror( error ) ); _exit(1); // crash and burn } if( (error = fclose(fp) ) > 0 ) { printf("fclose was beaten down by the mob in %s line %d: %s \n", __FILE__ , __LINE__ , strerror( error ) ); _exit(1); // crash and burn } finalReport(); free(categoryList); printf("Total Revenue Gained: \x1b[32m$\x1b[0m%.2lf\n", totalProfit); return 0; //success }
int main(int ac, char** av){ if(ac < 2){ std::cout << "usage: ./signalAcceptance inputFile[s]" << std::endl; return -1; } TH1F* allCategory = new TH1F("allCategory","all Category",11,0.5,11.5); TH1F* preselCategory = new TH1F("preselCategory","presel Category",11,0.5,11.5); TH1F* photonCategory = new TH1F("photonCategory","reco photon Category",11,0.5,11.5); TH1F* VisAllCategory = new TH1F("VisAllCategory","all Category, Vis",11,0.5,11.5); TH1F* VisPreselCategory = new TH1F("VisPreselCategory","presel Category, Vis",11,0.5,11.5); TH1F* VisPhotonCategory = new TH1F("VisPhotonCategory","reco photon Category, Vis",11,0.5,11.5); TH1F* dROtherGen = new TH1F("dROtherGen", "dROtherGen", 800, 0.0, 4.0); TH1F* parentage = new TH1F("parentage","parentage",30, 0, 30); TH1F* dptOverpt = new TH1F("dptOverpt","dptOverpt", 400, -2.0, 2.0); TH1F* dRrecoGen = new TH1F("dRrecoGen","dRrecoGen", 200, 0.0, 0.2); TH1F* dPhiRecoGen = new TH1F("dPhiRecoGen","dPhiRecoGen", 400, 0.0, 0.2); TH1F* dEtaRecoGen = new TH1F("dEtaRecoGen","dEtaRecoGen", 800, -0.2, 0.2); TH1F* dRGenNearJet = new TH1F("dRGenNearJet","dRGenNearJet", 200, 0.0, 1.0); TH1F* dPhiGenNearJet = new TH1F("dPhiGenNearJet","dPhiGenNearJet", 100, 0.0, 0.5); TH1F* dEtaGenNearJet = new TH1F("dEtaGenNearJet","dEtaGenNearJet", 200, -0.5, 0.5); //TH1F* dRGenNextNearJet = new TH1F("dRGenNextNearJet","dRGenNextNearJet", 600, 0.0, 6.0); //TH1F* dPhiGenNextNearJet = new TH1F("dPhiGenNextNearJet","dPhiGenNextNearJet", 300, 0.0, 3.0); //TH1F* dEtaGenNextNearJet = new TH1F("dEtaGenNextNearJet","dEtaGenNextNearJet", 600, -3.0, 3.0); // object selector Selector* selectorLoose = new Selector(); // create event selectors here EventPick* evtPickLoose = new EventPick("LoosePhotonID"); // do not do jet to photon dR cleaning evtPickLoose->veto_pho_jet_dR = 0.0; EventTree* tree = new EventTree(ac-1, av+1); double PUweight = 1.0; Long64_t nEntr = tree->GetEntries(); for(Long64_t entry=0; entry<nEntr; entry++){ if(entry%10000 == 0) std::cout << "processing entry " << entry << " out of " << nEntr << std::endl; tree->GetEntry(entry); doJER(tree); selectorLoose->process_objects(tree); evtPickLoose->process_event(tree, selectorLoose, PUweight); // fill the histograms fillCategory(tree, allCategory, PUweight); if(evtPickLoose->passPreSel) fillCategory(tree, preselCategory, PUweight); if(evtPickLoose->passAll) fillCategory(tree, photonCategory, PUweight); // fill histograms for gen photon passing the acceptance cuts defined in analysis bool inAcc = false; for(int mcInd=0; mcInd<tree->nMC_; ++mcInd){ if(tree->mcPID->at(mcInd) == 22 && (tree->mcParentage->at(mcInd)==2 || tree->mcParentage->at(mcInd)==10 || tree->mcParentage->at(mcInd)==26) && tree->mcPt->at(mcInd) > 25 && fabs(tree->mcEta->at(mcInd)) < 1.4442){ inAcc = true; } } if(inAcc){ fillCategory(tree, VisAllCategory, PUweight); if(evtPickLoose->passPreSel) fillCategory(tree, VisPreselCategory, PUweight); if(evtPickLoose->passAll) fillCategory(tree, VisPhotonCategory, PUweight); } // have at least one good photon if(!evtPickLoose->passAll) continue; // test //if(overlapMadGraph(tree)) continue; int phoInd = evtPickLoose->Photons.at(0); // experiment with delta R cuts for photons for(int mcInd=0; mcInd<tree->nMC_; ++mcInd){ bool etetamatch = dR(tree->mcEta->at(mcInd),tree->mcPhi->at(mcInd),tree->phoEta_->at(phoInd),tree->phoPhi_->at(phoInd)) < 0.2 && (fabs(tree->phoEt_->at(phoInd) - tree->mcPt->at(mcInd)) / tree->mcPt->at(mcInd)) < 1.0; if( etetamatch && tree->mcPID->at(mcInd) == 22){ // test if(!(tree->mcParentage->at(mcInd)==2 || tree->mcParentage->at(mcInd)==10 || tree->mcParentage->at(mcInd)==26)) continue; // fill histograms for mathced photon candidate parentage->Fill( tree->mcParentage->at(mcInd) ); dptOverpt->Fill( (tree->phoEt_->at(phoInd) - tree->mcPt->at(mcInd)) / tree->mcPt->at(mcInd)); dRrecoGen->Fill( dR(tree->mcEta->at(mcInd),tree->mcPhi->at(mcInd),tree->phoEta_->at(phoInd),tree->phoPhi_->at(phoInd)) ); dPhiRecoGen->Fill( dPhi( tree->phoPhi_->at(phoInd) - tree->mcPhi->at(mcInd) ) ); dEtaRecoGen->Fill( tree->phoEta_->at(phoInd) - tree->mcEta->at(mcInd) ); int closestGenInd = secondMinDrIndex( mcInd, tree ); if(dR(tree->mcEta->at(mcInd), tree->mcPhi->at(mcInd), tree->mcEta->at(closestGenInd), tree->mcPhi->at(closestGenInd)) < 0.01){ std::cout << "closest PID " << tree->mcPID->at(closestGenInd) << " MomPID " << tree->mcMomPID->at(closestGenInd) << std::endl; std::cout << "photon mother PID " << tree->mcMomPID->at(mcInd) << std::endl; } dROtherGen->Fill( dR(tree->mcEta->at(mcInd), tree->mcPhi->at(mcInd), tree->mcEta->at(closestGenInd), tree->mcPhi->at(closestGenInd)) ); int closestJetInd = minDrIndex( tree->mcEta->at(mcInd), tree->mcPhi->at(mcInd), tree->jetEta_, tree->jetPhi_ ); dRGenNearJet->Fill( dR(tree->mcEta->at(mcInd), tree->mcPhi->at(mcInd), tree->jetEta_->at(closestJetInd), tree->jetPhi_->at(closestJetInd) ) ); dPhiGenNearJet->Fill( dPhi( tree->jetPhi_->at(closestJetInd) - tree->mcPhi->at(mcInd) ) ); dEtaGenNearJet->Fill( tree->jetEta_->at(closestJetInd) - tree->mcEta->at(mcInd) ); //closestJetInd = secondMinDrIndex( tree->mcEta->at(mcInd), tree->mcPhi->at(mcInd), tree->jetEta_, tree->jetPhi_ ); //dRGenNextNearJet->Fill( dR(tree->mcEta->at(mcInd), tree->mcPhi->at(mcInd), tree->jetEta_->at(closestJetInd), tree->jetPhi_->at(closestJetInd) ) ); //dPhiGenNextNearJet->Fill( dPhi( tree->jetPhi_->at(closestJetInd) - tree->mcPhi->at(mcInd) ) ); //dEtaGenNextNearJet->Fill( tree->jetEta_->at(closestJetInd) - tree->mcEta->at(mcInd) ); } } } evtPickLoose->print_cutflow(); // write histograms TFile outFile("signalAcc.root","RECREATE"); saveHist(allCategory, &outFile); saveHist(preselCategory, &outFile); saveHist(photonCategory, &outFile); saveHist(VisAllCategory, &outFile); saveHist(VisPreselCategory, &outFile); saveHist(VisPhotonCategory, &outFile); saveHist(dROtherGen, &outFile); saveHist(parentage, &outFile); saveHist(dptOverpt, &outFile); saveHist(dRrecoGen, &outFile); saveHist(dPhiRecoGen, &outFile); saveHist(dEtaRecoGen, &outFile); saveHist(dRGenNearJet, &outFile); saveHist(dPhiGenNearJet, &outFile); saveHist(dEtaGenNearJet, &outFile); //saveHist(dRGenNextNearJet, &outFile); //saveHist(dPhiGenNextNearJet, &outFile); //saveHist(dEtaGenNextNearJet, &outFile); outFile.Close(); delete tree; return 0; }