Selection criteria

[maxwellcui]

Event should only be selected when b-jet present, but the code has an ambiguous criteria.

for(Int_t i=0;i<nentries;i++)
    {
      oldTree->GetEntry(i);
      bjet=0;
      if(SSee_2016 || SSem_2016 || SSmm_2016 || eee_2016 || eem_2016 || emm_2016 || mmm_2016)
        {
          for(unsigned int ibjet=0;ibjet<jet_mv2c10->size();ibjet++)
            {
              ///     if  (jet_mv2c10->at(ibjet) > 0.1758475) {  // 85% WP 
          if  (jet_mv2c10->at(ibjet) > 0.645925) 
        {  // 77% WP
          //if (jet_mv2c10->at(ibjet) > 0.8244273) {  // 70% WP
        bjet++;
        }      
            }
      bjetBranch->Fill();
        }
    }

[maxwellcui]

The code is now changed to the following, where b-jet is a hard requirement.

  for(Long64_t i=0;i<nentries;i++)
    {
      //Get entry
      oldTree->GetEntry(i);
      bjet=0;                    //Initilization

      //B-tagging
      if(SSee_2016 || SSem_2016 || SSmm_2016 || eee_2016 || eem_2016 || emm_2016 || mmm_2016)
        {
          for(unsigned int ibjet=0;ibjet<jet_mv2c10->size();ibjet++)
            {
              //0.1758475 85% WP
          //0.8244273 70% WP
          if  (jet_mv2c10->at(ibjet) > 0.645925)  // 77% WP          
        { 
          bjet++;
        }      
            }
        }

      //Event weight calculation
      //Calculate event weight
      if(bjet!=0)
    {
      evtWeight=weight_mc*weight_jvt*(weight_leptonSF_tightLeps/weight_indiv_SF_MU_TTVA)*weight_pileup*weight_bTagSF_77*lumi/mcnorm;
      newTree->Fill();
    }

    }