void calc_SBRatio_DecMod_DchKpipi() {

#include <iostream.h>
#include <fstream.h>

  gROOT->Reset();

  // Root macro to calculate uncertainty of signal
  //#######################################################
  // Change these values to meet your condition
  //#######################################################
 
  // Info of input data/MC root files
  TString baseDir("/home1x/OtherMounts/hep03/khamano/R16b/Run3");
  //TString baseDir("/home1x/OtherMounts/hep03/khamano/R16b/kCut1");
  TString nameBase_BpBm("SP-1235-BToDlnu-Run3-R16b");
  int nStart_BpBm = 1;
  int nStop_BpBm  = 10; //1193 //87
  TString nameBase_B0B0bar("SP-1237-BToDlnu-Run3-R16b");
  int nStart_B0B0bar = 100;
  int nStop_B0B0bar  = 10; //855 //88
  TString nameBase_ccbar("SP-1005-BToDlnu-Run3-R16b");
  int nStart_ccbar = 100;
  int nStop_ccbar  = 10; //208 //83
  TString nameBase_uds("SP-998-BToDlnu-Run3-R16b");
  int nStart_uds = 1000;
  int nStop_uds  = 5; //143 //58
  TString nameBase_tautau("SP-3429-BToDlnu-Run3-R16b");
  int nStart_tautau = 100;
  int nStop_tautau  = 2; //11 //5
  // Name of the root tree in the input root files
  TString ntpName("ntp2");


  // Info of output files
  TString outFileName("test_DchKpipi_Run3_BpBm_DecMod.txt");

  // D mass peak region
  float DMassMin = 1.845;  //for Dch
  float DMassMax = 1.893; 

  // D mass off peak (side band)region
  float DSideMin1 = 1.821;  //for Dch
  float DSideMax1 = 1.845; 
  float DSideMin2 = 1.893;  //for Dch
  float DSideMax2 = 1.917; 

  // Select D decay mode
  int setDDecMod = 5; // DchKpipi
  //int setDDecMod = 6; // DchKspipi0
  //int setDDecMod = 7; // DchKspi

  // Vertex cut
  float DprobMin = 0.01; //0.001
  float BprobMin = 0.01; //0.001

  // Thrust cut
  float cosThrMax = 0.80;

  // Select a bin
  float DMomMin = 0.2;
  float DMomMax = 2.8;
  float LMomMin = 0.8;
  float LMomMax = 2.8;
  float cosBYMin = -10;
  float cosBYMax = 5;

  // uncertainty on background
  float rho = 0;

  //#######################################################################

  // Chain input root files and define trees.
  cout << "start chain\n";

  TChain * chain_MC = new TChain(ntpName);
  for (Int_t i=nStart_BpBm; i<=nStop_BpBm; i++) {
  chain_MC->Add(baseDir+"/"+nameBase_BpBm+"-"+i+".root");
  }
  for (Int_t i=nStart_B0B0bar; i<=nStop_B0B0bar; i++) {
  chain_MC->Add(baseDir+"/"+nameBase_B0B0bar+"-"+i+".root");
  }
  for (Int_t i=nStart_ccbar; i<=nStop_ccbar; i++) {
  chain_MC->Add(baseDir+"/"+nameBase_ccbar+"-"+i+".root");
  }
  for (Int_t i=nStart_uds; i<=nStop_uds; i++) {
  chain_MC->Add(baseDir+"/"+nameBase_uds+"-"+i+".root");
  }
  for (Int_t i=nStart_tautau; i<=nStop_tautau; i++) {
  chain_MC->Add(baseDir+"/"+nameBase_tautau+"-"+i+".root");
  }
  TTree * tree_MC = chain_MC;
  
  int numDl(0);
  int DDecMod[200], DType[200];
  double fitDprob[200], fitBprob[200];
  float cosThr[200];
  int isABEvent[200], isDirectFromB[200], isFromSameB[200];
  int BDecMod[200], BDecMod2[200], BDecMod3[200]; 
  int mcDDecMod[200], mcD[200];
  int noExtKPi[200], KMissID[200], PiMissID[200], KPiSameD[200];
  float Dmass[200];
  float pD[200], pLep[200], cosBY[200];
  tree_MC->SetBranchAddress("nDl", &numDl);
  tree_MC->SetBranchAddress("DDecayMode", &DDecMod);
  tree_MC->SetBranchAddress("DhadType", &DType);
  tree_MC->SetBranchAddress("FitDprob", &fitDprob);
  tree_MC->SetBranchAddress("FitBprob", &fitBprob);
  tree_MC->SetBranchAddress("cosDlnonDl", &cosThr);
  tree_MC->SetBranchAddress("trueLepFromB", &isABEvent);
  tree_MC->SetBranchAddress("mcDirectLepFromB", &isDirectFromB);
  tree_MC->SetBranchAddress("trueLepFromSameB", &isFromSameB);
  tree_MC->SetBranchAddress("trueBDecMod", &BDecMod);
  tree_MC->SetBranchAddress("trueBDecMod2", &BDecMod2);
  tree_MC->SetBranchAddress("trueBDecMod3", &BDecMod3);
  tree_MC->SetBranchAddress("trueDDecMod", &mcDDecMod);
  tree_MC->SetBranchAddress("trueD", &mcD);
  tree_MC->SetBranchAddress("noExtraKorPi", &noExtKPi);
  tree_MC->SetBranchAddress("mcKaonMissID", &KMissID);
  tree_MC->SetBranchAddress("mcPiMissID", &PiMissID);
  tree_MC->SetBranchAddress("mcKPifromSameD", &KPiSameD);
  tree_MC->SetBranchAddress("DMass", &Dmass);
  tree_MC->SetBranchAddress("DMom", &pD);
  tree_MC->SetBranchAddress("DlLeptMom", &pLep);
  tree_MC->SetBranchAddress("DLcosBY", &cosBY);
  
  int realD_Peak(0), CascL_SameB_Peak(0), CascL_DiffB_Peak(0), DiffB_Peak;
  int realD_nDToKenu_Peak(0), realD_nDToKstarenu_Peak(0);
  int CascL_SameB_nDToKenu_Peak(0), CascL_SameB_nDToKstarenu_Peak(0);
  int CascL_SameB_nDDX_Peak(0);
  int CascL_SameB_nDD_Peak(0), CascL_SameB_nDDs_Peak(0);
  int CascL_SameB_nDDK_Peak(0), CascL_SameB_nDDPi_Peak(0);
  int CascL_SameB_nDDsK_Peak(0), CascL_SameB_nDDsPi_Peak(0);
  int CascL_DiffB_nDToKenu_Peak(0), CascL_DiffB_nDToKstarenu_Peak(0);
  int CascL_DiffB_nDDX_Peak(0);
  int CascL_DiffB_nDD_Peak(0), CascL_DiffB_nDDs_Peak(0);
  int CascL_DiffB_nDDK_Peak(0), CascL_DiffB_nDDPi_Peak(0);
  int CascL_DiffB_nDDsK_Peak(0), CascL_DiffB_nDDsPi_Peak(0);
  int DiffB_nDDX_Peak(0);
  int DiffB_nDD_Peak(0), DiffB_nDDs_Peak(0);
  int DiffB_nDDK_Peak(0), DiffB_nDDPi_Peak(0);
  int DiffB_nDDsK_Peak(0), DiffB_nDDsPi_Peak(0);

  Int_t nEvent =(Int_t)tree_MC->GetEntries(); // get number of events
  cout << "nEvent = " << nEvent << endl;
  int nSelEvent(0);
  for (Int_t i=0;i<nEvent;i++) { // loop over events
    if (i==1) cout << "event 1" << endl;
    if (i==100) cout << "event 100" << endl;
    if (i==1000) cout << "event 1000" << endl;
    if (i==10000) cout << "event 10000" << endl;
    if (i==100000) cout << "event 100000" << endl;
    if (i==1000000) cout << "event 1000000" << endl;
    if (i==2000000) cout << "event 2000000" << endl;
    if (i==3000000) cout << "event 3000000" << endl;
    if (i==4000000) cout << "event 4000000" << endl;

    tree_MC->GetEntry(i); //read complete event in memory
    int nDl_total(0), nDl_sig(0), nDl_bkg(0);
    int nDl_side(0);

    for (Int_t j=0;j<numDl;j++) { // loop over Dl candidates
      int DDecMod1 = DDecMod[j];
      int DType1 = DType[j];
      int BDecMod1 = BDecMod[j];
      int BDecMod21 = BDecMod2[j];
      int BDecMod31 = BDecMod3[j];
      double fitDprob1 = fitDprob[j]; 
      double fitBprob1 = fitBprob[j]; 
      float cosThr1 = cosThr[j];
      //cout << "cosThr1 = " << cosThr1 << endl;
      int isABEvent1 = isABEvent[j];
      int isDirectFromB1 = isDirectFromB[j];
      int isFromSameB1 = isFromSameB[j];
      int mcDDecMod1 = mcDDecMod[j];
      int mcD1 = mcD[j];
      int noExtKPi1 = noExtKPi[j];
      int KMissID1 = KMissID[j];
      int PiMissID1 = PiMissID[j];
      int KPiSameD1 = KPiSameD[j];
      float Dmass1 = Dmass[j];
      float pD1 = pD[j];
      float pLep1 = pLep[j];
      float cosBY1 = cosBY[j];

      // D decay mode selection
      if (DDecMod1 != setDDecMod) continue;
      if (DType1 == 1) continue;

      // Vertex cut
      //if (fitDprob1 < DprobMin) continue;
      //if (fitBprob1 < BprobMin) continue;
      
      // Thrust cut
      //if (cosThr1<-cosThrMax || cosThr1>cosThrMax) continue;
      
      // DMass cut
      if (Dmass1<DSideMin1 || DSideMax2<Dmass1) continue;
      if (Dmass1<DMassMin || DMassMax<Dmass1) continue; // DMass side band

      // true D cut
      //if (mcD1 == 0) continue;
      
      // ccbar trueD
      if (isABEvent1 == 0) {
	if (KMissID1 == 0 && PiMissID1 == 0 && KPiSameD1 == 1
	    && noExtKPi1 == 1) {
	  realD_Peak++;
	  if (BDecMod1 == 6111) {realD_nDToKenu_Peak++;}
	  if (BDecMod1 == 6131) {realD_nDToKstarenu_Peak++;}
	}
      }

      // CascL_SameB or DiffB  
      if (isABEvent1 == 1) {
	if (KMissID1 == 0 && PiMissID1 == 0 && KPiSameD1 == 1 
	    && noExtKPi1 == 1 && BDecMod1 != 0) {
	  if (isDirectFromB1 == 0 && isFromSameB1 == 1) { // CascL_SameB
	    CascL_SameB_Peak++;
	    if (BDecMod1 == 6111) {CascL_SameB_nDToKenu_Peak++;}
	    if (BDecMod1 == 6131) {CascL_SameB_nDToKstarenu_Peak++;}
	    if (BDecMod21 == 61110 || BDecMod21 == 61220 ||
		BDecMod21 == 61120 || BDecMod21 == 61210 ||
		BDecMod21 == 61330 || BDecMod21 == 61440 ||
		BDecMod21 == 61130 || BDecMod21 == 61310 ||
		BDecMod21 == 61140 || BDecMod21 == 61410 ||
		BDecMod21 == 61230 || BDecMod21 == 61320 ||
		BDecMod21 == 61240 || BDecMod21 == 61420 ||
		BDecMod21 == 61340 || BDecMod21 == 61430 ) {
	      CascL_SameB_nDD_Peak++;
	    }
	    if (BDecMod21 == 62110 || BDecMod21 == 62220 ||
		BDecMod21 == 62120 || BDecMod21 == 62210 ||
		BDecMod21 == 62330 || BDecMod21 == 62440 ||
		BDecMod21 == 62130 || BDecMod21 == 62310 ||
		BDecMod21 == 62140 || BDecMod21 == 62410 ||
		BDecMod21 == 62230 || BDecMod21 == 62320 ||
		BDecMod21 == 62240 || BDecMod21 == 62420 ||
		BDecMod21 == 62340 || BDecMod21 == 62430 ) {
	      CascL_SameB_nDD_Peak++;
	    }
	    if (BDecMod21 == 63110 || BDecMod21 == 63220 ||
		BDecMod21 == 63120 || BDecMod21 == 63210 ||
		BDecMod21 == 63330 || BDecMod21 == 63440 ||
		BDecMod21 == 63130 || BDecMod21 == 63310 ||
		BDecMod21 == 63140 || BDecMod21 == 63410 ||
		BDecMod21 == 63230 || BDecMod21 == 63320 ||
		BDecMod21 == 63240 || BDecMod21 == 63420 ||
		BDecMod21 == 63340 || BDecMod21 == 63430 ) {
	      CascL_SameB_nDD_Peak++;
	    }
	    if (BDecMod21 == 64110 || BDecMod21 == 64220 ||
		BDecMod21 == 64120 || BDecMod21 == 64210 ||
		BDecMod21 == 64330 || BDecMod21 == 64440 ||
		BDecMod21 == 64130 || BDecMod21 == 64310 ||
		BDecMod21 == 64140 || BDecMod21 == 64410 ||
		BDecMod21 == 64230 || BDecMod21 == 64320 ||
		BDecMod21 == 64240 || BDecMod21 == 64420 ||
		BDecMod21 == 64340 || BDecMod21 == 64430 ) {
	      CascL_SameB_nDD_Peak++;
	    }

	    if (BDecMod21 == 61160 || BDecMod21 == 61610 ||
		BDecMod21 == 61170 || BDecMod21 == 61710 ||
		BDecMod21 == 61260 || BDecMod21 == 61620 ||
		BDecMod21 == 61270 || BDecMod21 == 61720 ||
		BDecMod21 == 61360 || BDecMod21 == 61630 ||
		BDecMod21 == 61370 || BDecMod21 == 61730 ||
		BDecMod21 == 61460 || BDecMod21 == 61640 ||
		BDecMod21 == 61470 || BDecMod21 == 61740 ) {
	      CascL_SameB_nDDs_Peak++;
	    }
	    if (BDecMod21 == 62160 || BDecMod21 == 62610 ||
		BDecMod21 == 62170 || BDecMod21 == 62710 ||
		BDecMod21 == 62260 || BDecMod21 == 62620 ||
		BDecMod21 == 62270 || BDecMod21 == 62720 ||
		BDecMod21 == 62360 || BDecMod21 == 62630 ||
		BDecMod21 == 62370 || BDecMod21 == 62730 ||
		BDecMod21 == 62460 || BDecMod21 == 62640 ||
		BDecMod21 == 62470 || BDecMod21 == 62740 ) {
	      CascL_SameB_nDDs_Peak++;
	    }
	    if (BDecMod21 == 63160 || BDecMod21 == 63610 ||
		BDecMod21 == 63170 || BDecMod21 == 63710 ||
		BDecMod21 == 63260 || BDecMod21 == 63620 ||
		BDecMod21 == 63270 || BDecMod21 == 63720 ||
		BDecMod21 == 63360 || BDecMod21 == 63630 ||
		BDecMod21 == 63370 || BDecMod21 == 63730 ||
		BDecMod21 == 63460 || BDecMod21 == 63640 ||
		BDecMod21 == 63470 || BDecMod21 == 63740 ) {
	      CascL_SameB_nDDs_Peak++;
	    }
	    if (BDecMod21 == 64160 || BDecMod21 == 64610 ||
		BDecMod21 == 64170 || BDecMod21 == 64710 ||
		BDecMod21 == 64260 || BDecMod21 == 64620 ||
		BDecMod21 == 64270 || BDecMod21 == 64720 ||
		BDecMod21 == 64360 || BDecMod21 == 64630 ||
		BDecMod21 == 64370 || BDecMod21 == 64730 ||
		BDecMod21 == 64460 || BDecMod21 == 64640 ||
		BDecMod21 == 64470 || BDecMod21 == 64740 ) {
	      CascL_SameB_nDDs_Peak++;
	    }

	    if ((BDecMod21 >= 61111 && BDecMod21 <= 61116) ||
		(BDecMod21 >= 61221 && BDecMod21 <= 61226) ||
		(BDecMod21 >= 61121 && BDecMod21 <= 61126) ||
		(BDecMod21 >= 61211 && BDecMod21 <= 61216) ||
		(BDecMod21 >= 61331 && BDecMod21 <= 61336) ||
		(BDecMod21 >= 61441 && BDecMod21 <= 61446) ||
		(BDecMod21 >= 61131 && BDecMod21 <= 61136) ||
		(BDecMod21 >= 61311 && BDecMod21 <= 61316) ||
		(BDecMod21 >= 61141 && BDecMod21 <= 61146) ||
		(BDecMod21 >= 61411 && BDecMod21 <= 61416) ||
		(BDecMod21 >= 61231 && BDecMod21 <= 61236) ||
		(BDecMod21 >= 61321 && BDecMod21 <= 61326) ||
		(BDecMod21 >= 61241 && BDecMod21 <= 61246) ||
		(BDecMod21 >= 61421 && BDecMod21 <= 61426) ||
		(BDecMod21 >= 61341 && BDecMod21 <= 61346) ||
		 (BDecMod21 >= 61431 && BDecMod21 <= 61436)) {
	      CascL_SameB_nDDK_Peak++;
	    }	    
	    if ((BDecMod21 >= 61117 && BDecMod21 <= 61119) ||
		(BDecMod21 >= 61227 && BDecMod21 <= 61229) ||
		(BDecMod21 >= 61127 && BDecMod21 <= 61129) ||
		(BDecMod21 >= 61217 && BDecMod21 <= 61219) ||
		(BDecMod21 >= 61337 && BDecMod21 <= 61339) ||
		(BDecMod21 >= 61447 && BDecMod21 <= 61449) ||
		(BDecMod21 >= 61137 && BDecMod21 <= 61139) ||
		(BDecMod21 >= 61317 && BDecMod21 <= 61319) ||
		(BDecMod21 >= 61147 && BDecMod21 <= 61149) ||
		(BDecMod21 >= 61417 && BDecMod21 <= 61419) ||
		(BDecMod21 >= 61237 && BDecMod21 <= 61239) ||
		(BDecMod21 >= 61327 && BDecMod21 <= 61329) ||
		(BDecMod21 >= 61247 && BDecMod21 <= 61249) ||
		(BDecMod21 >= 61427 && BDecMod21 <= 61429) ||
		(BDecMod21 >= 61347 && BDecMod21 <= 61349) ||
		(BDecMod21 >= 61437 && BDecMod21 <= 61439)) {
	      CascL_SameB_nDDPi_Peak++;
	    }	    
	    if ((BDecMod21 >= 61161 && BDecMod21 <= 61166) || 
		(BDecMod21 >= 61611 && BDecMod21 <= 61616) ||
		(BDecMod21 >= 61171 && BDecMod21 <= 61176) || 
		(BDecMod21 >= 61711 && BDecMod21 <= 61716) ||
		(BDecMod21 >= 61261 && BDecMod21 <= 61266) || 
		(BDecMod21 >= 61621 && BDecMod21 <= 61626) ||
		(BDecMod21 >= 61271 && BDecMod21 <= 61276) || 
		(BDecMod21 >= 61721 && BDecMod21 <= 61726) ||
		(BDecMod21 >= 61361 && BDecMod21 <= 61366) || 
		(BDecMod21 >= 61631 && BDecMod21 <= 61636) ||
		(BDecMod21 >= 61371 && BDecMod21 <= 61376) || 
		(BDecMod21 >= 61731 && BDecMod21 <= 61736) ||
		(BDecMod21 >= 61461 && BDecMod21 <= 61466) || 
		(BDecMod21 >= 61641 && BDecMod21 <= 61646) ||
		(BDecMod21 >= 61471 && BDecMod21 <= 61476) || 
		(BDecMod21 >= 61741 && BDecMod21 <= 61746)) {
	      CascL_SameB_nDDsK_Peak++;
	    }
	    if ((BDecMod21 >= 61167 && BDecMod21 <= 61169) || 
		(BDecMod21 >= 61617 && BDecMod21 <= 61619) ||
		(BDecMod21 >= 61177 && BDecMod21 <= 61179) || 
		(BDecMod21 >= 61717 && BDecMod21 <= 61719) ||
		(BDecMod21 >= 61267 && BDecMod21 <= 61269) || 
		(BDecMod21 >= 61627 && BDecMod21 <= 61629) ||
		(BDecMod21 >= 61277 && BDecMod21 <= 61279) || 
		(BDecMod21 >= 61727 && BDecMod21 <= 61729) ||
		(BDecMod21 >= 61367 && BDecMod21 <= 61369) || 
		(BDecMod21 >= 61637 && BDecMod21 <= 61639) ||
		(BDecMod21 >= 61377 && BDecMod21 <= 61379) || 
		(BDecMod21 >= 61737 && BDecMod21 <= 61739) ||
		(BDecMod21 >= 61467 && BDecMod21 <= 61469) || 
		(BDecMod21 >= 61647 && BDecMod21 <= 61649) ||
		(BDecMod21 >= 61477 && BDecMod21 <= 61479) || 
		(BDecMod21 >= 61747 && BDecMod21 <= 61749)) {
	      CascL_SameB_nDDsPi_Peak++;
	    }

	    if ((BDecMod21 >= 62111 && BDecMod21 <= 62116) ||
		(BDecMod21 >= 62221 && BDecMod21 <= 62226) ||
		(BDecMod21 >= 62121 && BDecMod21 <= 62126) ||
		(BDecMod21 >= 62211 && BDecMod21 <= 62216) ||
		(BDecMod21 >= 62331 && BDecMod21 <= 62336) ||
		(BDecMod21 >= 62441 && BDecMod21 <= 62446) ||
		(BDecMod21 >= 62131 && BDecMod21 <= 62136) ||
		(BDecMod21 >= 62311 && BDecMod21 <= 62316) ||
		(BDecMod21 >= 62141 && BDecMod21 <= 62146) ||
		(BDecMod21 >= 62411 && BDecMod21 <= 62416) ||
		(BDecMod21 >= 62231 && BDecMod21 <= 62236) ||
		(BDecMod21 >= 62321 && BDecMod21 <= 62326) ||
		(BDecMod21 >= 62241 && BDecMod21 <= 62246) ||
		(BDecMod21 >= 62421 && BDecMod21 <= 62426) ||
		(BDecMod21 >= 62341 && BDecMod21 <= 62346) ||
		(BDecMod21 >= 62431 && BDecMod21 <= 62436)) {
		  CascL_SameB_nDDK_Peak++;
		}	    
	    if ((BDecMod21 >= 62117 && BDecMod21 <= 62119) ||
		(BDecMod21 >= 62227 && BDecMod21 <= 62229) ||
		(BDecMod21 >= 62127 && BDecMod21 <= 62129) ||
		(BDecMod21 >= 62217 && BDecMod21 <= 62219) ||
		(BDecMod21 >= 62337 && BDecMod21 <= 62339) ||
		(BDecMod21 >= 62447 && BDecMod21 <= 62449) ||
		(BDecMod21 >= 62137 && BDecMod21 <= 62139) ||
		(BDecMod21 >= 62317 && BDecMod21 <= 62319) ||
		(BDecMod21 >= 62147 && BDecMod21 <= 62149) ||
		(BDecMod21 >= 62417 && BDecMod21 <= 62419) ||
		(BDecMod21 >= 62237 && BDecMod21 <= 62239) ||
		(BDecMod21 >= 62327 && BDecMod21 <= 62329) ||
		(BDecMod21 >= 62247 && BDecMod21 <= 62249) ||
		(BDecMod21 >= 62427 && BDecMod21 <= 62429) ||
		(BDecMod21 >= 62347 && BDecMod21 <= 62349) ||
		(BDecMod21 >= 62437 && BDecMod21 <= 62439)) {
		  CascL_SameB_nDDPi_Peak++;
		}	    
	    if ((BDecMod21 >= 62161 && BDecMod21 <= 62166) || 
		(BDecMod21 >= 62611 && BDecMod21 <= 62616) ||
		(BDecMod21 >= 62171 && BDecMod21 <= 62176) || 
		(BDecMod21 >= 62711 && BDecMod21 <= 62716) ||
		(BDecMod21 >= 62261 && BDecMod21 <= 62266) || 
		(BDecMod21 >= 62621 && BDecMod21 <= 62626) ||
		(BDecMod21 >= 62271 && BDecMod21 <= 62276) || 
		(BDecMod21 >= 62721 && BDecMod21 <= 62726) ||
		(BDecMod21 >= 62361 && BDecMod21 <= 62366) || 
		(BDecMod21 >= 62631 && BDecMod21 <= 62636) ||
		(BDecMod21 >= 62371 && BDecMod21 <= 62376) || 
		(BDecMod21 >= 62731 && BDecMod21 <= 62736) ||
		(BDecMod21 >= 62461 && BDecMod21 <= 62466) || 
		(BDecMod21 >= 62641 && BDecMod21 <= 62646) ||
		(BDecMod21 >= 62471 && BDecMod21 <= 62476) || 
		(BDecMod21 >= 62741 && BDecMod21 <= 62746)) {
	      CascL_SameB_nDDsK_Peak++;
	    }
	    if ((BDecMod21 >= 62167 && BDecMod21 <= 62169) || 
		(BDecMod21 >= 62617 && BDecMod21 <= 62619) ||
		(BDecMod21 >= 62177 && BDecMod21 <= 62179) || 
		(BDecMod21 >= 62717 && BDecMod21 <= 62719) ||
		(BDecMod21 >= 62267 && BDecMod21 <= 62269) || 
		(BDecMod21 >= 62627 && BDecMod21 <= 62629) ||
		(BDecMod21 >= 62277 && BDecMod21 <= 62279) || 
		(BDecMod21 >= 62727 && BDecMod21 <= 62729) ||
		(BDecMod21 >= 62367 && BDecMod21 <= 62369) || 
		(BDecMod21 >= 62637 && BDecMod21 <= 62639) ||
		(BDecMod21 >= 62377 && BDecMod21 <= 62379) || 
		(BDecMod21 >= 62737 && BDecMod21 <= 62739) ||
		(BDecMod21 >= 62467 && BDecMod21 <= 62469) || 
		(BDecMod21 >= 62647 && BDecMod21 <= 62649) ||
		(BDecMod21 >= 62477 && BDecMod21 <= 62479) || 
		(BDecMod21 >= 62747 && BDecMod21 <= 62749)) {
	      CascL_SameB_nDDsPi_Peak++;
	    }

	    if ((BDecMod21 >= 63111 && BDecMod21 <= 63116) ||
		(BDecMod21 >= 63221 && BDecMod21 <= 63226) ||
		(BDecMod21 >= 63121 && BDecMod21 <= 63126) ||
		(BDecMod21 >= 63211 && BDecMod21 <= 63216) ||
		(BDecMod21 >= 63331 && BDecMod21 <= 63336) ||
		(BDecMod21 >= 63441 && BDecMod21 <= 63446) ||
		(BDecMod21 >= 63131 && BDecMod21 <= 63136) ||
		(BDecMod21 >= 63311 && BDecMod21 <= 63316) ||
		(BDecMod21 >= 63141 && BDecMod21 <= 63146) ||
		(BDecMod21 >= 63411 && BDecMod21 <= 63416) ||
		(BDecMod21 >= 63231 && BDecMod21 <= 63236) ||
		(BDecMod21 >= 63321 && BDecMod21 <= 63326) ||
		(BDecMod21 >= 63241 && BDecMod21 <= 63246) ||
		(BDecMod21 >= 63421 && BDecMod21 <= 63426) ||
		(BDecMod21 >= 63341 && BDecMod21 <= 63346) ||
		(BDecMod21 >= 63431 && BDecMod21 <= 63436)) {
	      CascL_SameB_nDDK_Peak++;
	    }	    
	    if ((BDecMod21 >= 63117 && BDecMod21 <= 63119) ||
		(BDecMod21 >= 63227 && BDecMod21 <= 63229) ||
		(BDecMod21 >= 63127 && BDecMod21 <= 63129) ||
		(BDecMod21 >= 63217 && BDecMod21 <= 63219) ||
		(BDecMod21 >= 63337 && BDecMod21 <= 63339) ||
		(BDecMod21 >= 63447 && BDecMod21 <= 63449) ||
		(BDecMod21 >= 63137 && BDecMod21 <= 63139) ||
		(BDecMod21 >= 63317 && BDecMod21 <= 63319) ||
		(BDecMod21 >= 63147 && BDecMod21 <= 63149) ||
		(BDecMod21 >= 63417 && BDecMod21 <= 63419) ||
		(BDecMod21 >= 63237 && BDecMod21 <= 63239) ||
		(BDecMod21 >= 63327 && BDecMod21 <= 63329) ||
		(BDecMod21 >= 63247 && BDecMod21 <= 63249) ||
		(BDecMod21 >= 63427 && BDecMod21 <= 63429) ||
		(BDecMod21 >= 63347 && BDecMod21 <= 63349) ||
		(BDecMod21 >= 63437 && BDecMod21 <= 63439)) {
	      CascL_SameB_nDDPi_Peak++;
	    }	    
	    if ((BDecMod21 >= 63161 && BDecMod21 <= 63166) || 
		(BDecMod21 >= 63611 && BDecMod21 <= 63616) ||
		(BDecMod21 >= 63171 && BDecMod21 <= 63176) || 
		(BDecMod21 >= 63711 && BDecMod21 <= 63716) ||
		(BDecMod21 >= 63261 && BDecMod21 <= 63266) || 
		(BDecMod21 >= 63621 && BDecMod21 <= 63626) ||
		(BDecMod21 >= 63271 && BDecMod21 <= 63276) || 
		(BDecMod21 >= 63721 && BDecMod21 <= 63726) ||
		(BDecMod21 >= 63361 && BDecMod21 <= 63366) || 
		(BDecMod21 >= 63631 && BDecMod21 <= 63636) ||
		(BDecMod21 >= 63371 && BDecMod21 <= 63376) || 
		(BDecMod21 >= 63731 && BDecMod21 <= 63736) ||
		(BDecMod21 >= 63461 && BDecMod21 <= 63466) || 
		(BDecMod21 >= 63641 && BDecMod21 <= 63646) ||
		(BDecMod21 >= 63471 && BDecMod21 <= 63476) || 
		(BDecMod21 >= 63741 && BDecMod21 <= 63746)) {
	      CascL_SameB_nDDsK_Peak++;
	    }
	    if ((BDecMod21 >= 63167 && BDecMod21 <= 63169) || 
		(BDecMod21 >= 63617 && BDecMod21 <= 63619) ||
		(BDecMod21 >= 63177 && BDecMod21 <= 63179) || 
		(BDecMod21 >= 63717 && BDecMod21 <= 63719) ||
		(BDecMod21 >= 63267 && BDecMod21 <= 63269) || 
		(BDecMod21 >= 63627 && BDecMod21 <= 63629) ||
		(BDecMod21 >= 63277 && BDecMod21 <= 63279) || 
		(BDecMod21 >= 63727 && BDecMod21 <= 63729) ||
		(BDecMod21 >= 63367 && BDecMod21 <= 63369) || 
		(BDecMod21 >= 63637 && BDecMod21 <= 63639) ||
		(BDecMod21 >= 63377 && BDecMod21 <= 63379) || 
		(BDecMod21 >= 63737 && BDecMod21 <= 63739) ||
		(BDecMod21 >= 63467 && BDecMod21 <= 63469) || 
		(BDecMod21 >= 63647 && BDecMod21 <= 63649) ||
		(BDecMod21 >= 63477 && BDecMod21 <= 63479) || 
		(BDecMod21 >= 63747 && BDecMod21 <= 63749)) {
	      CascL_SameB_nDDsPi_Peak++;
	    }

	    if ((BDecMod21 >= 64111 && BDecMod21 <= 64116) ||
		(BDecMod21 >= 64221 && BDecMod21 <= 64226) ||
		(BDecMod21 >= 64121 && BDecMod21 <= 64126) ||
		(BDecMod21 >= 64211 && BDecMod21 <= 64216) ||
		(BDecMod21 >= 64331 && BDecMod21 <= 64336) ||
		(BDecMod21 >= 64441 && BDecMod21 <= 64446) ||
		(BDecMod21 >= 64131 && BDecMod21 <= 64136) ||
		(BDecMod21 >= 64311 && BDecMod21 <= 64316) ||
		(BDecMod21 >= 64141 && BDecMod21 <= 64146) ||
		(BDecMod21 >= 64411 && BDecMod21 <= 64416) ||
		(BDecMod21 >= 64231 && BDecMod21 <= 64236) ||
		(BDecMod21 >= 64321 && BDecMod21 <= 64326) ||
		(BDecMod21 >= 64241 && BDecMod21 <= 64246) ||
		(BDecMod21 >= 64421 && BDecMod21 <= 64426) ||
		(BDecMod21 >= 64341 && BDecMod21 <= 64346) ||
		(BDecMod21 >= 64431 && BDecMod21 <= 64436)) {
	      CascL_SameB_nDDK_Peak++;
	    }	    
	    if ((BDecMod21 >= 64117 && BDecMod21 <= 64119) ||
		(BDecMod21 >= 64227 && BDecMod21 <= 64229) ||
		(BDecMod21 >= 64127 && BDecMod21 <= 64129) ||
		(BDecMod21 >= 64217 && BDecMod21 <= 64219) ||
		(BDecMod21 >= 64337 && BDecMod21 <= 64339) ||
		(BDecMod21 >= 64447 && BDecMod21 <= 64449) ||
		(BDecMod21 >= 64137 && BDecMod21 <= 64139) ||
		(BDecMod21 >= 64317 && BDecMod21 <= 64319) ||
		(BDecMod21 >= 64147 && BDecMod21 <= 64149) ||
		(BDecMod21 >= 64417 && BDecMod21 <= 64419) ||
		(BDecMod21 >= 64237 && BDecMod21 <= 64239) ||
		(BDecMod21 >= 64327 && BDecMod21 <= 64329) ||
		(BDecMod21 >= 64247 && BDecMod21 <= 64249) ||
		(BDecMod21 >= 64427 && BDecMod21 <= 64429) ||
		(BDecMod21 >= 64347 && BDecMod21 <= 64349) ||
		(BDecMod21 >= 64437 && BDecMod21 <= 64439)) {
	      CascL_SameB_nDDPi_Peak++;
	    }	    
	    if ((BDecMod21 >= 64161 && BDecMod21 <= 64166) || 
		(BDecMod21 >= 64611 && BDecMod21 <= 64616) ||
		(BDecMod21 >= 64171 && BDecMod21 <= 64176) || 
		(BDecMod21 >= 64711 && BDecMod21 <= 64716) ||
		(BDecMod21 >= 64261 && BDecMod21 <= 64266) || 
		(BDecMod21 >= 64621 && BDecMod21 <= 64626) ||
		(BDecMod21 >= 64271 && BDecMod21 <= 64276) || 
		(BDecMod21 >= 64721 && BDecMod21 <= 64726) ||
		(BDecMod21 >= 64361 && BDecMod21 <= 64366) || 
		(BDecMod21 >= 64631 && BDecMod21 <= 64636) ||
		(BDecMod21 >= 64371 && BDecMod21 <= 64376) || 
		(BDecMod21 >= 64731 && BDecMod21 <= 64736) ||
		(BDecMod21 >= 64461 && BDecMod21 <= 64466) || 
		(BDecMod21 >= 64641 && BDecMod21 <= 64646) ||
		(BDecMod21 >= 64471 && BDecMod21 <= 64476) || 
		(BDecMod21 >= 64741 && BDecMod21 <= 64746)) {
	      CascL_SameB_nDDsK_Peak++;
	    }
	    if ((BDecMod21 >= 64167 && BDecMod21 <= 64169) || 
		(BDecMod21 >= 64617 && BDecMod21 <= 64619) ||
		(BDecMod21 >= 64177 && BDecMod21 <= 64179) || 
		(BDecMod21 >= 64717 && BDecMod21 <= 64719) ||
		(BDecMod21 >= 64267 && BDecMod21 <= 64269) || 
		(BDecMod21 >= 64627 && BDecMod21 <= 64629) ||
		(BDecMod21 >= 64277 && BDecMod21 <= 64279) || 
		(BDecMod21 >= 64727 && BDecMod21 <= 64729) ||
		(BDecMod21 >= 64367 && BDecMod21 <= 64369) || 
		(BDecMod21 >= 64637 && BDecMod21 <= 64639) ||
		(BDecMod21 >= 64377 && BDecMod21 <= 64379) || 
		(BDecMod21 >= 64737 && BDecMod21 <= 64739) ||
		(BDecMod21 >= 64467 && BDecMod21 <= 64469) || 
		(BDecMod21 >= 64647 && BDecMod21 <= 64649) ||
		(BDecMod21 >= 64477 && BDecMod21 <= 64479) || 
		(BDecMod21 >= 64747 && BDecMod21 <= 64749)) {
	      CascL_SameB_nDDsPi_Peak++;
	    }
	  }
	  
	  if (isDirectFromB1 == 0 && isFromSameB1 == 0) { // CascL_DiffB
	    CascL_DiffB_Peak++;
	    if (BDecMod1 == 6111) CascL_DiffB_nDToKenu_Peak++;
	    if (BDecMod1 == 6131) CascL_DiffB_nDToKstarenu_Peak++;
	    if (BDecMod31 == 1110 || BDecMod31 == 1220 ||
		BDecMod31 == 1120 || BDecMod31 == 1210 ||
		BDecMod31 == 1330 || BDecMod31 == 1440 ||
		BDecMod31 == 1130 || BDecMod31 == 1310 ||
		BDecMod31 == 1140 || BDecMod31 == 1410 ||
		BDecMod31 == 1230 || BDecMod31 == 1320 ||
		BDecMod31 == 1240 || BDecMod31 == 1420 ||
		BDecMod31 == 1340 || BDecMod31 == 1430 ) {
	      CascL_DiffB_nDD_Peak++;
	    }
	    if (BDecMod31 == 2110 || BDecMod31 == 2220 ||
		BDecMod31 == 2120 || BDecMod31 == 2210 ||
		BDecMod31 == 2330 || BDecMod31 == 2440 ||
		BDecMod31 == 2130 || BDecMod31 == 2310 ||
		BDecMod31 == 2140 || BDecMod31 == 2410 ||
		BDecMod31 == 2230 || BDecMod31 == 2320 ||
		BDecMod31 == 2240 || BDecMod31 == 2420 ||
		BDecMod31 == 2340 || BDecMod31 == 2430 ) {
	      CascL_DiffB_nDD_Peak++;
	    }
	    if (BDecMod31 == 3110 || BDecMod31 == 3220 ||
		BDecMod31 == 3120 || BDecMod31 == 3210 ||
		BDecMod31 == 3330 || BDecMod31 == 3440 ||
		BDecMod31 == 3130 || BDecMod31 == 3310 ||
		BDecMod31 == 3140 || BDecMod31 == 3410 ||
		BDecMod31 == 3230 || BDecMod31 == 3320 ||
		BDecMod31 == 3240 || BDecMod31 == 3420 ||
		BDecMod31 == 3340 || BDecMod31 == 3430 ) {
	      CascL_DiffB_nDD_Peak++;
	    }
	    if (BDecMod31 == 4110 || BDecMod31 == 4220 ||
		BDecMod31 == 4120 || BDecMod31 == 4210 ||
		BDecMod31 == 4330 || BDecMod31 == 4440 ||
		BDecMod31 == 4130 || BDecMod31 == 4310 ||
		BDecMod31 == 4140 || BDecMod31 == 4410 ||
		BDecMod31 == 4230 || BDecMod31 == 4320 ||
		BDecMod31 == 4240 || BDecMod31 == 4420 ||
		BDecMod31 == 4340 || BDecMod31 == 4430 ) {
	      CascL_DiffB_nDD_Peak++;
	    }

	    if (BDecMod31 == 1160 || BDecMod31 == 1610 ||
		BDecMod31 == 1170 || BDecMod31 == 1710 ||
		BDecMod31 == 1260 || BDecMod31 == 1620 ||
		BDecMod31 == 1270 || BDecMod31 == 1720 ||
		BDecMod31 == 1360 || BDecMod31 == 1630 ||
		BDecMod31 == 1370 || BDecMod31 == 1730 ||
		BDecMod31 == 1460 || BDecMod31 == 1640 ||
		BDecMod31 == 1470 || BDecMod31 == 1740 ) {
	      CascL_DiffB_nDDs_Peak++;
	    }
	    if (BDecMod31 == 2160 || BDecMod31 == 2610 ||
		BDecMod31 == 2170 || BDecMod31 == 2710 ||
		BDecMod31 == 2260 || BDecMod31 == 2620 ||
		BDecMod31 == 2270 || BDecMod31 == 2720 ||
		BDecMod31 == 2360 || BDecMod31 == 2630 ||
		BDecMod31 == 2370 || BDecMod31 == 2730 ||
		BDecMod31 == 2460 || BDecMod31 == 2640 ||
		BDecMod31 == 2470 || BDecMod31 == 2740 ) {
	      CascL_DiffB_nDDs_Peak++;
	    }
	    if (BDecMod31 == 3160 || BDecMod31 == 3610 ||
		BDecMod31 == 3170 || BDecMod31 == 3710 ||
		BDecMod31 == 3260 || BDecMod31 == 3620 ||
		BDecMod31 == 3270 || BDecMod31 == 3720 ||
		BDecMod31 == 3360 || BDecMod31 == 3630 ||
		BDecMod31 == 3370 || BDecMod31 == 3730 ||
		BDecMod31 == 3460 || BDecMod31 == 3640 ||
		BDecMod31 == 3470 || BDecMod31 == 3740 ) {
	      CascL_DiffB_nDDs_Peak++;
	    }
	    if (BDecMod31 == 4160 || BDecMod31 == 4610 ||
		BDecMod31 == 4170 || BDecMod31 == 4710 ||
		BDecMod31 == 4260 || BDecMod31 == 4620 ||
		BDecMod31 == 4270 || BDecMod31 == 4720 ||
		BDecMod31 == 4360 || BDecMod31 == 4630 ||
		BDecMod31 == 4370 || BDecMod31 == 4730 ||
		BDecMod31 == 4460 || BDecMod31 == 4640 ||
		BDecMod31 == 4470 || BDecMod31 == 4740 ) {
	      CascL_DiffB_nDDs_Peak++;
	    }

	    if ((BDecMod31 >= 1111 && BDecMod31 <= 1116) ||
		(BDecMod31 >= 1221 && BDecMod31 <= 1226) ||
		(BDecMod31 >= 1121 && BDecMod31 <= 1126) ||
		(BDecMod31 >= 1211 && BDecMod31 <= 1216) ||
		(BDecMod31 >= 1331 && BDecMod31 <= 1336) ||
		(BDecMod31 >= 1441 && BDecMod31 <= 1446) ||
		(BDecMod31 >= 1131 && BDecMod31 <= 1136) ||
		(BDecMod31 >= 1311 && BDecMod31 <= 1316) ||
		(BDecMod31 >= 1141 && BDecMod31 <= 1146) ||
		(BDecMod31 >= 1411 && BDecMod31 <= 1416) ||
		(BDecMod31 >= 1231 && BDecMod31 <= 1236) ||
		(BDecMod31 >= 1321 && BDecMod31 <= 1326) ||
		(BDecMod31 >= 1241 && BDecMod31 <= 1246) ||
		(BDecMod31 >= 1421 && BDecMod31 <= 1426) ||
		(BDecMod31 >= 1341 && BDecMod31 <= 1346) ||
		(BDecMod31 >= 1431 && BDecMod31 <= 1436)) {
	      CascL_DiffB_nDDK_Peak++;
	    }	    
	    if ((BDecMod31 >= 1117 && BDecMod31 <= 1119) ||
		(BDecMod31 >= 1227 && BDecMod31 <= 1229) ||
		(BDecMod31 >= 1127 && BDecMod31 <= 1129) ||
		(BDecMod31 >= 1217 && BDecMod31 <= 1219) ||
		(BDecMod31 >= 1337 && BDecMod31 <= 1339) ||
		(BDecMod31 >= 1447 && BDecMod31 <= 1449) ||
		(BDecMod31 >= 1137 && BDecMod31 <= 1139) ||
		(BDecMod31 >= 1317 && BDecMod31 <= 1319) ||
		(BDecMod31 >= 1147 && BDecMod31 <= 1149) ||
		(BDecMod31 >= 1417 && BDecMod31 <= 1419) ||
		(BDecMod31 >= 1237 && BDecMod31 <= 1239) ||
		(BDecMod31 >= 1327 && BDecMod31 <= 1329) ||
		(BDecMod31 >= 1247 && BDecMod31 <= 1249) ||
		(BDecMod31 >= 1427 && BDecMod31 <= 1429) ||
		(BDecMod31 >= 1347 && BDecMod31 <= 1349) ||
		(BDecMod31 >= 1437 && BDecMod31 <= 1439)) {
	      CascL_DiffB_nDDPi_Peak++;
	    }	    
	    if ((BDecMod31 >= 1161 && BDecMod31 <= 1166) || 
		(BDecMod31 >= 1611 && BDecMod31 <= 1616) ||
		(BDecMod31 >= 1171 && BDecMod31 <= 1176) || 
		(BDecMod31 >= 1711 && BDecMod31 <= 1716) ||
		(BDecMod31 >= 1261 && BDecMod31 <= 1266) || 
		(BDecMod31 >= 1621 && BDecMod31 <= 1626) ||
		(BDecMod31 >= 1271 && BDecMod31 <= 1276) || 
		(BDecMod31 >= 1721 && BDecMod31 <= 1726) ||
		(BDecMod31 >= 1361 && BDecMod31 <= 1366) || 
		(BDecMod31 >= 1631 && BDecMod31 <= 1636) ||
		(BDecMod31 >= 1371 && BDecMod31 <= 1376) || 
		(BDecMod31 >= 1731 && BDecMod31 <= 1736) ||
		(BDecMod31 >= 1461 && BDecMod31 <= 1466) || 
		(BDecMod31 >= 1641 && BDecMod31 <= 1646) ||
		(BDecMod31 >= 1471 && BDecMod31 <= 1476) || 
		(BDecMod31 >= 1741 && BDecMod31 <= 1746)) {
	      CascL_DiffB_nDDsK_Peak++;
	    }
	    if ((BDecMod31 >= 1167 && BDecMod31 <= 1169) || 
		(BDecMod31 >= 1617 && BDecMod31 <= 1619) ||
		(BDecMod31 >= 1177 && BDecMod31 <= 1179) || 
		(BDecMod31 >= 1717 && BDecMod31 <= 1719) ||
		(BDecMod31 >= 1267 && BDecMod31 <= 1269) || 
		(BDecMod31 >= 1627 && BDecMod31 <= 1629) ||
		(BDecMod31 >= 1277 && BDecMod31 <= 1279) || 
		(BDecMod31 >= 1727 && BDecMod31 <= 1729) ||
		(BDecMod31 >= 1367 && BDecMod31 <= 1369) || 
		(BDecMod31 >= 1637 && BDecMod31 <= 1639) ||
		(BDecMod31 >= 1377 && BDecMod31 <= 1379) || 
		(BDecMod31 >= 1737 && BDecMod31 <= 1739) ||
		(BDecMod31 >= 1467 && BDecMod31 <= 1469) || 
		(BDecMod31 >= 1647 && BDecMod31 <= 1649) ||
		(BDecMod31 >= 1477 && BDecMod31 <= 1479) || 
		(BDecMod31 >= 1747 && BDecMod31 <= 1749)) {
	      CascL_DiffB_nDDsPi_Peak++;
	    }

	    if ((BDecMod31 >= 2111 && BDecMod31 <= 2116) ||
		(BDecMod31 >= 2221 && BDecMod31 <= 2226) ||
		(BDecMod31 >= 2121 && BDecMod31 <= 2126) ||
		(BDecMod31 >= 2211 && BDecMod31 <= 2216) ||
		(BDecMod31 >= 2331 && BDecMod31 <= 2336) ||
		(BDecMod31 >= 2441 && BDecMod31 <= 2446) ||
		(BDecMod31 >= 2131 && BDecMod31 <= 2136) ||
		(BDecMod31 >= 2311 && BDecMod31 <= 2316) ||
		(BDecMod31 >= 2141 && BDecMod31 <= 2146) ||
		(BDecMod31 >= 2411 && BDecMod31 <= 2416) ||
		(BDecMod31 >= 2231 && BDecMod31 <= 2236) ||
		(BDecMod31 >= 2321 && BDecMod31 <= 2326) ||
		(BDecMod31 >= 2241 && BDecMod31 <= 2246) ||
		(BDecMod31 >= 2421 && BDecMod31 <= 2426) ||
		(BDecMod31 >= 2341 && BDecMod31 <= 2346) ||
		(BDecMod31 >= 2431 && BDecMod31 <= 2436)) {
	      CascL_DiffB_nDDK_Peak++;
	    }	    
	    if ((BDecMod31 >= 2117 && BDecMod31 <= 2119) ||
		(BDecMod31 >= 2227 && BDecMod31 <= 2229) ||
		(BDecMod31 >= 2127 && BDecMod31 <= 2129) ||
		(BDecMod31 >= 2217 && BDecMod31 <= 2219) ||
		(BDecMod31 >= 2337 && BDecMod31 <= 2339) ||
		(BDecMod31 >= 2447 && BDecMod31 <= 2449) ||
		(BDecMod31 >= 2137 && BDecMod31 <= 2139) ||
		(BDecMod31 >= 2317 && BDecMod31 <= 2319) ||
		(BDecMod31 >= 2147 && BDecMod31 <= 2149) ||
		(BDecMod31 >= 2417 && BDecMod31 <= 2419) ||
		(BDecMod31 >= 2237 && BDecMod31 <= 2239) ||
		(BDecMod31 >= 2327 && BDecMod31 <= 2329) ||
		(BDecMod31 >= 2247 && BDecMod31 <= 2249) ||
		(BDecMod31 >= 2427 && BDecMod31 <= 2429) ||
		(BDecMod31 >= 2347 && BDecMod31 <= 2349) ||
		(BDecMod31 >= 2437 && BDecMod31 <= 2439)) {
	      CascL_DiffB_nDDPi_Peak++;
	    }	    
	    if ((BDecMod31 >= 2161 && BDecMod31 <= 2166) || 
		(BDecMod31 >= 2611 && BDecMod31 <= 2616) ||
		(BDecMod31 >= 2171 && BDecMod31 <= 2176) || 
		(BDecMod31 >= 2711 && BDecMod31 <= 2716) ||
		(BDecMod31 >= 2261 && BDecMod31 <= 2266) || 
		(BDecMod31 >= 2621 && BDecMod31 <= 2626) ||
		(BDecMod31 >= 2271 && BDecMod31 <= 2276) || 
		(BDecMod31 >= 2721 && BDecMod31 <= 2726) ||
		(BDecMod31 >= 2361 && BDecMod31 <= 2366) || 
		(BDecMod31 >= 2631 && BDecMod31 <= 2636) ||
		(BDecMod31 >= 2371 && BDecMod31 <= 2376) || 
		(BDecMod31 >= 2731 && BDecMod31 <= 2736) ||
		(BDecMod31 >= 2461 && BDecMod31 <= 2466) || 
		(BDecMod31 >= 2641 && BDecMod31 <= 2646) ||
		(BDecMod31 >= 2471 && BDecMod31 <= 2476) || 
		(BDecMod31 >= 2741 && BDecMod31 <= 2746)) {
	      CascL_DiffB_nDDsK_Peak++;
	    }
	    if ((BDecMod31 >= 2167 && BDecMod31 <= 2169) || 
		(BDecMod31 >= 2617 && BDecMod31 <= 2619) ||
		(BDecMod31 >= 2177 && BDecMod31 <= 2179) || 
		(BDecMod31 >= 2717 && BDecMod31 <= 2719) ||
		(BDecMod31 >= 2267 && BDecMod31 <= 2269) || 
		(BDecMod31 >= 2627 && BDecMod31 <= 2629) ||
		(BDecMod31 >= 2277 && BDecMod31 <= 2279) || 
		(BDecMod31 >= 2727 && BDecMod31 <= 2729) ||
		(BDecMod31 >= 2367 && BDecMod31 <= 2369) || 
		(BDecMod31 >= 2637 && BDecMod31 <= 2639) ||
		(BDecMod31 >= 2377 && BDecMod31 <= 2379) || 
		(BDecMod31 >= 2737 && BDecMod31 <= 2739) ||
		(BDecMod31 >= 2467 && BDecMod31 <= 2469) || 
		(BDecMod31 >= 2647 && BDecMod31 <= 2649) ||
		(BDecMod31 >= 2477 && BDecMod31 <= 2479) || 
		(BDecMod31 >= 2747 && BDecMod31 <= 2749)) {
	      CascL_DiffB_nDDsPi_Peak++;
	    }

	    if ((BDecMod31 >= 3111 && BDecMod31 <= 3116) ||
		(BDecMod31 >= 3221 && BDecMod31 <= 3226) ||
		(BDecMod31 >= 3121 && BDecMod31 <= 3126) ||
		(BDecMod31 >= 3211 && BDecMod31 <= 3216) ||
		(BDecMod31 >= 3331 && BDecMod31 <= 3336) ||
		(BDecMod31 >= 3441 && BDecMod31 <= 3446) ||
		(BDecMod31 >= 3131 && BDecMod31 <= 3136) ||
		(BDecMod31 >= 3311 && BDecMod31 <= 3316) ||
		(BDecMod31 >= 3141 && BDecMod31 <= 3146) ||
		(BDecMod31 >= 3411 && BDecMod31 <= 3416) ||
		(BDecMod31 >= 3231 && BDecMod31 <= 3236) ||
		(BDecMod31 >= 3321 && BDecMod31 <= 3326) ||
		(BDecMod31 >= 3241 && BDecMod31 <= 3246) ||
		(BDecMod31 >= 3421 && BDecMod31 <= 3426) ||
		(BDecMod31 >= 3341 && BDecMod31 <= 3346) ||
		(BDecMod31 >= 3431 && BDecMod31 <= 3436)) {
	      CascL_DiffB_nDDK_Peak++;
	    }	    
	    if ((BDecMod31 >= 3117 && BDecMod31 <= 3119) ||
		(BDecMod31 >= 3227 && BDecMod31 <= 3229) ||
		(BDecMod31 >= 3127 && BDecMod31 <= 3129) ||
		(BDecMod31 >= 3217 && BDecMod31 <= 3219) ||
		(BDecMod31 >= 3337 && BDecMod31 <= 3339) ||
		(BDecMod31 >= 3447 && BDecMod31 <= 3449) ||
		(BDecMod31 >= 3137 && BDecMod31 <= 3139) ||
		(BDecMod31 >= 3317 && BDecMod31 <= 3319) ||
		(BDecMod31 >= 3147 && BDecMod31 <= 3149) ||
		(BDecMod31 >= 3417 && BDecMod31 <= 3419) ||
		(BDecMod31 >= 3237 && BDecMod31 <= 3239) ||
		(BDecMod31 >= 3327 && BDecMod31 <= 3329) ||
		(BDecMod31 >= 3247 && BDecMod31 <= 3249) ||
		(BDecMod31 >= 3427 && BDecMod31 <= 3429) ||
		(BDecMod31 >= 3347 && BDecMod31 <= 3349) ||
		(BDecMod31 >= 3437 && BDecMod31 <= 3439)) {
	      CascL_DiffB_nDDPi_Peak++;
	    }	    
	    if ((BDecMod31 >= 3161 && BDecMod31 <= 3166) || 
		(BDecMod31 >= 3611 && BDecMod31 <= 3616) ||
		(BDecMod31 >= 3171 && BDecMod31 <= 3176) || 
		(BDecMod31 >= 3711 && BDecMod31 <= 3716) ||
		(BDecMod31 >= 3261 && BDecMod31 <= 3266) || 
		(BDecMod31 >= 3621 && BDecMod31 <= 3626) ||
		(BDecMod31 >= 3271 && BDecMod31 <= 3276) || 
		(BDecMod31 >= 3721 && BDecMod31 <= 3726) ||
		(BDecMod31 >= 3361 && BDecMod31 <= 3366) || 
		(BDecMod31 >= 3631 && BDecMod31 <= 3636) ||
		(BDecMod31 >= 3371 && BDecMod31 <= 3376) || 
		(BDecMod31 >= 3731 && BDecMod31 <= 3736) ||
		(BDecMod31 >= 3461 && BDecMod31 <= 3466) || 
		(BDecMod31 >= 3641 && BDecMod31 <= 3646) ||
		(BDecMod31 >= 3471 && BDecMod31 <= 3476) || 
		(BDecMod31 >= 3741 && BDecMod31 <= 3746)) {
	      CascL_DiffB_nDDsK_Peak++;
	    }
	    if ((BDecMod31 >= 3167 && BDecMod31 <= 3169) || 
		(BDecMod31 >= 3617 && BDecMod31 <= 3619) ||
		(BDecMod31 >= 3177 && BDecMod31 <= 3179) || 
		(BDecMod31 >= 3717 && BDecMod31 <= 3719) ||
		(BDecMod31 >= 3267 && BDecMod31 <= 3269) || 
		(BDecMod31 >= 3627 && BDecMod31 <= 3629) ||
		(BDecMod31 >= 3277 && BDecMod31 <= 3279) || 
		(BDecMod31 >= 3727 && BDecMod31 <= 3729) ||
		(BDecMod31 >= 3367 && BDecMod31 <= 3369) || 
		(BDecMod31 >= 3637 && BDecMod31 <= 3639) ||
		(BDecMod31 >= 3377 && BDecMod31 <= 3379) || 
		(BDecMod31 >= 3737 && BDecMod31 <= 3739) ||
		(BDecMod31 >= 3467 && BDecMod31 <= 3469) || 
		(BDecMod31 >= 3647 && BDecMod31 <= 3649) ||
		(BDecMod31 >= 3477 && BDecMod31 <= 3479) || 
		(BDecMod31 >= 3747 && BDecMod31 <= 3749)) {
	      CascL_DiffB_nDDsPi_Peak++;
	    }

	    if ((BDecMod31 >= 4111 && BDecMod31 <= 4116) ||
		(BDecMod31 >= 4221 && BDecMod31 <= 4226) ||
		(BDecMod31 >= 4121 && BDecMod31 <= 4126) ||
		(BDecMod31 >= 4211 && BDecMod31 <= 4216) ||
		(BDecMod31 >= 4331 && BDecMod31 <= 4336) ||
		(BDecMod31 >= 4441 && BDecMod31 <= 4446) ||
		(BDecMod31 >= 4131 && BDecMod31 <= 4136) ||
		(BDecMod31 >= 4311 && BDecMod31 <= 4316) ||
		(BDecMod31 >= 4141 && BDecMod31 <= 4146) ||
		(BDecMod31 >= 4411 && BDecMod31 <= 4416) ||
		(BDecMod31 >= 4231 && BDecMod31 <= 4236) ||
		(BDecMod31 >= 4321 && BDecMod31 <= 4326) ||
		(BDecMod31 >= 4241 && BDecMod31 <= 4246) ||
		(BDecMod31 >= 4421 && BDecMod31 <= 4426) ||
		(BDecMod31 >= 4341 && BDecMod31 <= 4346) ||
		(BDecMod31 >= 4431 && BDecMod31 <= 4436)) {
	      CascL_DiffB_nDDK_Peak++;
	    }	    
	    if ((BDecMod31 >= 4117 && BDecMod31 <= 4119) ||
		(BDecMod31 >= 4227 && BDecMod31 <= 4229) ||
		(BDecMod31 >= 4127 && BDecMod31 <= 4129) ||
		(BDecMod31 >= 4217 && BDecMod31 <= 4219) ||
		(BDecMod31 >= 4337 && BDecMod31 <= 4339) ||
		(BDecMod31 >= 4447 && BDecMod31 <= 4449) ||
		(BDecMod31 >= 4137 && BDecMod31 <= 4139) ||
		(BDecMod31 >= 4317 && BDecMod31 <= 4319) ||
		(BDecMod31 >= 4147 && BDecMod31 <= 4149) ||
		(BDecMod31 >= 4417 && BDecMod31 <= 4419) ||
		(BDecMod31 >= 4237 && BDecMod31 <= 4239) ||
		(BDecMod31 >= 4327 && BDecMod31 <= 4329) ||
		(BDecMod31 >= 4247 && BDecMod31 <= 4249) ||
		(BDecMod31 >= 4427 && BDecMod31 <= 4429) ||
		(BDecMod31 >= 4347 && BDecMod31 <= 4349) ||
		(BDecMod31 >= 4437 && BDecMod31 <= 4439)) {
	      CascL_DiffB_nDDPi_Peak++;
	    }	    
	    if ((BDecMod31 >= 4161 && BDecMod31 <= 4166) || 
		(BDecMod31 >= 4611 && BDecMod31 <= 4616) ||
		(BDecMod31 >= 4171 && BDecMod31 <= 4176) || 
		(BDecMod31 >= 4711 && BDecMod31 <= 4716) ||
		(BDecMod31 >= 4261 && BDecMod31 <= 4266) || 
		(BDecMod31 >= 4621 && BDecMod31 <= 4626) ||
		(BDecMod31 >= 4271 && BDecMod31 <= 4276) || 
		(BDecMod31 >= 4721 && BDecMod31 <= 4726) ||
		(BDecMod31 >= 4361 && BDecMod31 <= 4366) || 
		(BDecMod31 >= 4631 && BDecMod31 <= 4636) ||
		(BDecMod31 >= 4371 && BDecMod31 <= 4376) || 
		(BDecMod31 >= 4731 && BDecMod31 <= 4736) ||
		(BDecMod31 >= 4461 && BDecMod31 <= 4466) || 
		(BDecMod31 >= 4641 && BDecMod31 <= 4646) ||
		(BDecMod31 >= 4471 && BDecMod31 <= 4476) || 
		(BDecMod31 >= 4741 && BDecMod31 <= 4746)) {
	      CascL_DiffB_nDDsK_Peak++;
	    }
	    if ((BDecMod31 >= 4167 && BDecMod31 <= 4169) || 
		(BDecMod31 >= 4617 && BDecMod31 <= 4619) ||
		(BDecMod31 >= 4177 && BDecMod31 <= 4179) || 
		(BDecMod31 >= 4717 && BDecMod31 <= 4719) ||
		(BDecMod31 >= 4267 && BDecMod31 <= 4269) || 
		(BDecMod31 >= 4627 && BDecMod31 <= 4629) ||
		(BDecMod31 >= 4277 && BDecMod31 <= 4279) || 
		(BDecMod31 >= 4727 && BDecMod31 <= 4729) ||
		(BDecMod31 >= 4367 && BDecMod31 <= 4369) || 
		(BDecMod31 >= 4637 && BDecMod31 <= 4639) ||
		(BDecMod31 >= 4377 && BDecMod31 <= 4379) || 
		(BDecMod31 >= 4737 && BDecMod31 <= 4739) ||
		(BDecMod31 >= 4467 && BDecMod31 <= 4469) || 
		(BDecMod31 >= 4647 && BDecMod31 <= 4649) ||
		(BDecMod31 >= 4477 && BDecMod31 <= 4479) || 
		(BDecMod31 >= 4747 && BDecMod31 <= 4749)) {
	      CascL_DiffB_nDDsPi_Peak++;
	    }
	  } 
	       	  
	  if (isDirectFromB1 == 1 && isFromSameB1 == 0) { // DiffB
	    DiffB_Peak++;
	    if (BDecMod31 == 1110 || BDecMod31 == 1220 ||
		BDecMod31 == 1120 || BDecMod31 == 1210 ||
		BDecMod31 == 1330 || BDecMod31 == 1440 ||
		BDecMod31 == 1130 || BDecMod31 == 1310 ||
		BDecMod31 == 1140 || BDecMod31 == 1410 ||
		BDecMod31 == 1230 || BDecMod31 == 1320 ||
		BDecMod31 == 1240 || BDecMod31 == 1420 ||
		BDecMod31 == 1340 || BDecMod31 == 1430 ) DiffB_nDD_Peak++;
	    if (BDecMod31 == 2110 || BDecMod31 == 2220 ||
		BDecMod31 == 2120 || BDecMod31 == 2210 ||
		BDecMod31 == 2330 || BDecMod31 == 2440 ||
		BDecMod31 == 2130 || BDecMod31 == 2310 ||
		BDecMod31 == 2140 || BDecMod31 == 2410 ||
		BDecMod31 == 2230 || BDecMod31 == 2320 ||
		BDecMod31 == 2240 || BDecMod31 == 2420 ||
		BDecMod31 == 2340 || BDecMod31 == 2430 ) DiffB_nDD_Peak++;
	    if (BDecMod31 == 3110 || BDecMod31 == 3220 ||
		BDecMod31 == 3120 || BDecMod31 == 3210 ||
		BDecMod31 == 3330 || BDecMod31 == 3440 ||
		BDecMod31 == 3130 || BDecMod31 == 3310 ||
		BDecMod31 == 3140 || BDecMod31 == 3410 ||
		BDecMod31 == 3230 || BDecMod31 == 3320 ||
		BDecMod31 == 3240 || BDecMod31 == 3420 ||
		BDecMod31 == 3340 || BDecMod31 == 3430 ) DiffB_nDD_Peak++;
	    if (BDecMod31 == 4110 || BDecMod31 == 4220 ||
		BDecMod31 == 4120 || BDecMod31 == 4210 ||
		BDecMod31 == 4330 || BDecMod31 == 4440 ||
		BDecMod31 == 4130 || BDecMod31 == 4310 ||
		BDecMod31 == 4140 || BDecMod31 == 4410 ||
		BDecMod31 == 4230 || BDecMod31 == 4320 ||
		BDecMod31 == 4240 || BDecMod31 == 4420 ||
		BDecMod31 == 4340 || BDecMod31 == 4430 ) DiffB_nDD_Peak++;

	    if (BDecMod31 == 1160 || BDecMod31 == 1610 ||
		BDecMod31 == 1170 || BDecMod31 == 1710 ||
		BDecMod31 == 1260 || BDecMod31 == 1620 ||
		BDecMod31 == 1270 || BDecMod31 == 1720 ||
		BDecMod31 == 1360 || BDecMod31 == 1630 ||
		BDecMod31 == 1370 || BDecMod31 == 1730 ||
		BDecMod31 == 1460 || BDecMod31 == 1640 ||
		BDecMod31 == 1470 || BDecMod31 == 1740 ) DiffB_nDDs_Peak++;
	    if (BDecMod31 == 2160 || BDecMod31 == 2610 ||
		BDecMod31 == 2170 || BDecMod31 == 2710 ||
		BDecMod31 == 2260 || BDecMod31 == 2620 ||
		BDecMod31 == 2270 || BDecMod31 == 2720 ||
		BDecMod31 == 2360 || BDecMod31 == 2630 ||
		BDecMod31 == 2370 || BDecMod31 == 2730 ||
		BDecMod31 == 2460 || BDecMod31 == 2640 ||
		BDecMod31 == 2470 || BDecMod31 == 2740 ) DiffB_nDDs_Peak++;
	    if (BDecMod31 == 3160 || BDecMod31 == 3610 ||
		BDecMod31 == 3170 || BDecMod31 == 3710 ||
		BDecMod31 == 3260 || BDecMod31 == 3620 ||
		BDecMod31 == 3270 || BDecMod31 == 3720 ||
		BDecMod31 == 3360 || BDecMod31 == 3630 ||
		BDecMod31 == 3370 || BDecMod31 == 3730 ||
		BDecMod31 == 3460 || BDecMod31 == 3640 ||
		BDecMod31 == 3470 || BDecMod31 == 3740 ) DiffB_nDDs_Peak++;
	    if (BDecMod31 == 4160 || BDecMod31 == 4610 ||
		BDecMod31 == 4170 || BDecMod31 == 4710 ||
		BDecMod31 == 4260 || BDecMod31 == 4620 ||
		BDecMod31 == 4270 || BDecMod31 == 4720 ||
		BDecMod31 == 4360 || BDecMod31 == 4630 ||
		BDecMod31 == 4370 || BDecMod31 == 4730 ||
		BDecMod31 == 4460 || BDecMod31 == 4640 ||
		BDecMod31 == 4470 || BDecMod31 == 4740 ) DiffB_nDDs_Peak++;

	    if ((BDecMod31 >= 1111 && BDecMod31 <= 1116) ||
		(BDecMod31 >= 1221 && BDecMod31 <= 1226) ||
		(BDecMod31 >= 1121 && BDecMod31 <= 1126) ||
		(BDecMod31 >= 1211 && BDecMod31 <= 1216) ||
		(BDecMod31 >= 1331 && BDecMod31 <= 1336) ||
		(BDecMod31 >= 1441 && BDecMod31 <= 1446) ||
		(BDecMod31 >= 1131 && BDecMod31 <= 1136) ||
		(BDecMod31 >= 1311 && BDecMod31 <= 1316) ||
		(BDecMod31 >= 1141 && BDecMod31 <= 1146) ||
		(BDecMod31 >= 1411 && BDecMod31 <= 1416) ||
		(BDecMod31 >= 1231 && BDecMod31 <= 1236) ||
		(BDecMod31 >= 1321 && BDecMod31 <= 1326) ||
		(BDecMod31 >= 1241 && BDecMod31 <= 1246) ||
		(BDecMod31 >= 1421 && BDecMod31 <= 1426) ||
		(BDecMod31 >= 1341 && BDecMod31 <= 1346) ||
		(BDecMod31 >= 1431 && BDecMod31 <= 1436)) {
	      DiffB_nDDK_Peak++;
	    }	    
	    if ((BDecMod31 >= 1117 && BDecMod31 <= 1119) ||
		(BDecMod31 >= 1227 && BDecMod31 <= 1229) ||
		(BDecMod31 >= 1127 && BDecMod31 <= 1129) ||
		(BDecMod31 >= 1217 && BDecMod31 <= 1219) ||
		(BDecMod31 >= 1337 && BDecMod31 <= 1339) ||
		(BDecMod31 >= 1447 && BDecMod31 <= 1449) ||
		(BDecMod31 >= 1137 && BDecMod31 <= 1139) ||
		(BDecMod31 >= 1317 && BDecMod31 <= 1319) ||
		(BDecMod31 >= 1147 && BDecMod31 <= 1149) ||
		(BDecMod31 >= 1417 && BDecMod31 <= 1419) ||
		(BDecMod31 >= 1237 && BDecMod31 <= 1239) ||
		(BDecMod31 >= 1327 && BDecMod31 <= 1329) ||
		(BDecMod31 >= 1247 && BDecMod31 <= 1249) ||
		(BDecMod31 >= 1427 && BDecMod31 <= 1429) ||
		(BDecMod31 >= 1347 && BDecMod31 <= 1349) ||
		(BDecMod31 >= 1437 && BDecMod31 <= 1439)) {
	      DiffB_nDDPi_Peak++;
	    }	    
	    if ((BDecMod31 >= 1161 && BDecMod31 <= 1166) || 
		(BDecMod31 >= 1611 && BDecMod31 <= 1616) ||
		(BDecMod31 >= 1171 && BDecMod31 <= 1176) || 
		(BDecMod31 >= 1711 && BDecMod31 <= 1716) ||
		(BDecMod31 >= 1261 && BDecMod31 <= 1266) || 
		(BDecMod31 >= 1621 && BDecMod31 <= 1626) ||
		(BDecMod31 >= 1271 && BDecMod31 <= 1276) || 
		(BDecMod31 >= 1721 && BDecMod31 <= 1726) ||
		(BDecMod31 >= 1361 && BDecMod31 <= 1366) || 
		(BDecMod31 >= 1631 && BDecMod31 <= 1636) ||
		(BDecMod31 >= 1371 && BDecMod31 <= 1376) || 
		(BDecMod31 >= 1731 && BDecMod31 <= 1736) ||
		(BDecMod31 >= 1461 && BDecMod31 <= 1466) || 
		(BDecMod31 >= 1641 && BDecMod31 <= 1646) ||
		(BDecMod31 >= 1471 && BDecMod31 <= 1476) || 
		(BDecMod31 >= 1741 && BDecMod31 <= 1746)) {
	      DiffB_nDDsK_Peak++;
	    }
	    if ((BDecMod31 >= 1167 && BDecMod31 <= 1169) || 
		(BDecMod31 >= 1617 && BDecMod31 <= 1619) ||
		(BDecMod31 >= 1177 && BDecMod31 <= 1179) || 
		(BDecMod31 >= 1717 && BDecMod31 <= 1719) ||
		(BDecMod31 >= 1267 && BDecMod31 <= 1269) || 
		(BDecMod31 >= 1627 && BDecMod31 <= 1629) ||
		(BDecMod31 >= 1277 && BDecMod31 <= 1279) || 
		(BDecMod31 >= 1727 && BDecMod31 <= 1729) ||
		(BDecMod31 >= 1367 && BDecMod31 <= 1369) || 
		(BDecMod31 >= 1637 && BDecMod31 <= 1639) ||
		(BDecMod31 >= 1377 && BDecMod31 <= 1379) || 
		(BDecMod31 >= 1737 && BDecMod31 <= 1739) ||
		(BDecMod31 >= 1467 && BDecMod31 <= 1469) || 
		(BDecMod31 >= 1647 && BDecMod31 <= 1649) ||
		(BDecMod31 >= 1477 && BDecMod31 <= 1479) || 
		(BDecMod31 >= 1747 && BDecMod31 <= 1749)) {
	      DiffB_nDDsPi_Peak++;
	    }

	    if ((BDecMod31 >= 2111 && BDecMod31 <= 2116) ||
		(BDecMod31 >= 2221 && BDecMod31 <= 2226) ||
		(BDecMod31 >= 2121 && BDecMod31 <= 2126) ||
		(BDecMod31 >= 2211 && BDecMod31 <= 2216) ||
		(BDecMod31 >= 2331 && BDecMod31 <= 2336) ||
		(BDecMod31 >= 2441 && BDecMod31 <= 2446) ||
		(BDecMod31 >= 2131 && BDecMod31 <= 2136) ||
		(BDecMod31 >= 2311 && BDecMod31 <= 2316) ||
		(BDecMod31 >= 2141 && BDecMod31 <= 2146) ||
		(BDecMod31 >= 2411 && BDecMod31 <= 2416) ||
		(BDecMod31 >= 2231 && BDecMod31 <= 2236) ||
		(BDecMod31 >= 2321 && BDecMod31 <= 2326) ||
		(BDecMod31 >= 2241 && BDecMod31 <= 2246) ||
		(BDecMod31 >= 2421 && BDecMod31 <= 2426) ||
		(BDecMod31 >= 2341 && BDecMod31 <= 2346) ||
		(BDecMod31 >= 2431 && BDecMod31 <= 2436)) {
	      DiffB_nDDK_Peak++;
	    }	    
	    if ((BDecMod31 >= 2117 && BDecMod31 <= 2119) ||
		(BDecMod31 >= 2227 && BDecMod31 <= 2229) ||
		(BDecMod31 >= 2127 && BDecMod31 <= 2129) ||
		(BDecMod31 >= 2217 && BDecMod31 <= 2219) ||
		(BDecMod31 >= 2337 && BDecMod31 <= 2339) ||
		(BDecMod31 >= 2447 && BDecMod31 <= 2449) ||
		(BDecMod31 >= 2137 && BDecMod31 <= 2139) ||
		(BDecMod31 >= 2317 && BDecMod31 <= 2319) ||
		(BDecMod31 >= 2147 && BDecMod31 <= 2149) ||
		(BDecMod31 >= 2417 && BDecMod31 <= 2419) ||
		(BDecMod31 >= 2237 && BDecMod31 <= 2239) ||
		(BDecMod31 >= 2327 && BDecMod31 <= 2329) ||
		(BDecMod31 >= 2247 && BDecMod31 <= 2249) ||
		(BDecMod31 >= 2427 && BDecMod31 <= 2429) ||
		(BDecMod31 >= 2347 && BDecMod31 <= 2349) ||
		(BDecMod31 >= 2437 && BDecMod31 <= 2439)) {
	      DiffB_nDDPi_Peak++;
	    }	    
	    if ((BDecMod31 >= 2161 && BDecMod31 <= 2166) || 
		(BDecMod31 >= 2611 && BDecMod31 <= 2616) ||
		(BDecMod31 >= 2171 && BDecMod31 <= 2176) || 
		(BDecMod31 >= 2711 && BDecMod31 <= 2716) ||
		(BDecMod31 >= 2261 && BDecMod31 <= 2266) || 
		(BDecMod31 >= 2621 && BDecMod31 <= 2626) ||
		(BDecMod31 >= 2271 && BDecMod31 <= 2276) || 
		(BDecMod31 >= 2721 && BDecMod31 <= 2726) ||
		(BDecMod31 >= 2361 && BDecMod31 <= 2366) || 
		(BDecMod31 >= 2631 && BDecMod31 <= 2636) ||
		(BDecMod31 >= 2371 && BDecMod31 <= 2376) || 
		(BDecMod31 >= 2731 && BDecMod31 <= 2736) ||
		(BDecMod31 >= 2461 && BDecMod31 <= 2466) || 
		(BDecMod31 >= 2641 && BDecMod31 <= 2646) ||
		(BDecMod31 >= 2471 && BDecMod31 <= 2476) || 
		(BDecMod31 >= 2741 && BDecMod31 <= 2746)) {
	      DiffB_nDDsK_Peak++;
	    }
	    if ((BDecMod31 >= 2167 && BDecMod31 <= 2169) || 
		(BDecMod31 >= 2617 && BDecMod31 <= 2619) ||
		(BDecMod31 >= 2177 && BDecMod31 <= 2179) || 
		(BDecMod31 >= 2717 && BDecMod31 <= 2719) ||
		(BDecMod31 >= 2267 && BDecMod31 <= 2269) || 
		(BDecMod31 >= 2627 && BDecMod31 <= 2629) ||
		(BDecMod31 >= 2277 && BDecMod31 <= 2279) || 
		(BDecMod31 >= 2727 && BDecMod31 <= 2729) ||
		(BDecMod31 >= 2367 && BDecMod31 <= 2369) || 
		(BDecMod31 >= 2637 && BDecMod31 <= 2639) ||
		(BDecMod31 >= 2377 && BDecMod31 <= 2379) || 
		(BDecMod31 >= 2737 && BDecMod31 <= 2739) ||
		(BDecMod31 >= 2467 && BDecMod31 <= 2469) || 
		(BDecMod31 >= 2647 && BDecMod31 <= 2649) ||
		(BDecMod31 >= 2477 && BDecMod31 <= 2479) || 
		(BDecMod31 >= 2747 && BDecMod31 <= 2749)) {
	      DiffB_nDDsPi_Peak++;
	    }

	    if ((BDecMod31 >= 3111 && BDecMod31 <= 3116) ||
		(BDecMod31 >= 3221 && BDecMod31 <= 3226) ||
		(BDecMod31 >= 3121 && BDecMod31 <= 3126) ||
		(BDecMod31 >= 3211 && BDecMod31 <= 3216) ||
		(BDecMod31 >= 3331 && BDecMod31 <= 3336) ||
		(BDecMod31 >= 3441 && BDecMod31 <= 3446) ||
		(BDecMod31 >= 3131 && BDecMod31 <= 3136) ||
		(BDecMod31 >= 3311 && BDecMod31 <= 3316) ||
		(BDecMod31 >= 3141 && BDecMod31 <= 3146) ||
		(BDecMod31 >= 3411 && BDecMod31 <= 3416) ||
		(BDecMod31 >= 3231 && BDecMod31 <= 3236) ||
		(BDecMod31 >= 3321 && BDecMod31 <= 3326) ||
		(BDecMod31 >= 3241 && BDecMod31 <= 3246) ||
		(BDecMod31 >= 3421 && BDecMod31 <= 3426) ||
		(BDecMod31 >= 3341 && BDecMod31 <= 3346) ||
		(BDecMod31 >= 3431 && BDecMod31 <= 3436)) {
	      DiffB_nDDK_Peak++;
	    }	    
	    if ((BDecMod31 >= 3117 && BDecMod31 <= 3119) ||
		(BDecMod31 >= 3227 && BDecMod31 <= 3229) ||
		(BDecMod31 >= 3127 && BDecMod31 <= 3129) ||
		(BDecMod31 >= 3217 && BDecMod31 <= 3219) ||
		(BDecMod31 >= 3337 && BDecMod31 <= 3339) ||
		(BDecMod31 >= 3447 && BDecMod31 <= 3449) ||
		(BDecMod31 >= 3137 && BDecMod31 <= 3139) ||
		(BDecMod31 >= 3317 && BDecMod31 <= 3319) ||
		(BDecMod31 >= 3147 && BDecMod31 <= 3149) ||
		(BDecMod31 >= 3417 && BDecMod31 <= 3419) ||
		(BDecMod31 >= 3237 && BDecMod31 <= 3239) ||
		(BDecMod31 >= 3327 && BDecMod31 <= 3329) ||
		(BDecMod31 >= 3247 && BDecMod31 <= 3249) ||
		(BDecMod31 >= 3427 && BDecMod31 <= 3429) ||
		(BDecMod31 >= 3347 && BDecMod31 <= 3349) ||
		(BDecMod31 >= 3437 && BDecMod31 <= 3439)) {
	      DiffB_nDDPi_Peak++;
	    }	    
	    if ((BDecMod31 >= 3161 && BDecMod31 <= 3166) || 
		(BDecMod31 >= 3611 && BDecMod31 <= 3616) ||
		(BDecMod31 >= 3171 && BDecMod31 <= 3176) || 
		(BDecMod31 >= 3711 && BDecMod31 <= 3716) ||
		(BDecMod31 >= 3261 && BDecMod31 <= 3266) || 
		(BDecMod31 >= 3621 && BDecMod31 <= 3626) ||
		(BDecMod31 >= 3271 && BDecMod31 <= 3276) || 
		(BDecMod31 >= 3721 && BDecMod31 <= 3726) ||
		(BDecMod31 >= 3361 && BDecMod31 <= 3366) || 
		(BDecMod31 >= 3631 && BDecMod31 <= 3636) ||
		(BDecMod31 >= 3371 && BDecMod31 <= 3376) || 
		(BDecMod31 >= 3731 && BDecMod31 <= 3736) ||
		(BDecMod31 >= 3461 && BDecMod31 <= 3466) || 
		(BDecMod31 >= 3641 && BDecMod31 <= 3646) ||
		(BDecMod31 >= 3471 && BDecMod31 <= 3476) || 
		(BDecMod31 >= 3741 && BDecMod31 <= 3746)) {
	      DiffB_nDDsK_Peak++;
	    }
	    if ((BDecMod31 >= 3167 && BDecMod31 <= 3169) || 
		(BDecMod31 >= 3617 && BDecMod31 <= 3619) ||
		(BDecMod31 >= 3177 && BDecMod31 <= 3179) || 
		(BDecMod31 >= 3717 && BDecMod31 <= 3719) ||
		(BDecMod31 >= 3267 && BDecMod31 <= 3269) || 
		(BDecMod31 >= 3627 && BDecMod31 <= 3629) ||
		(BDecMod31 >= 3277 && BDecMod31 <= 3279) || 
		(BDecMod31 >= 3727 && BDecMod31 <= 3729) ||
		(BDecMod31 >= 3367 && BDecMod31 <= 3369) || 
		(BDecMod31 >= 3637 && BDecMod31 <= 3639) ||
		(BDecMod31 >= 3377 && BDecMod31 <= 3379) || 
		(BDecMod31 >= 3737 && BDecMod31 <= 3739) ||
		(BDecMod31 >= 3467 && BDecMod31 <= 3469) || 
		(BDecMod31 >= 3647 && BDecMod31 <= 3649) ||
		(BDecMod31 >= 3477 && BDecMod31 <= 3479) || 
		(BDecMod31 >= 3747 && BDecMod31 <= 3749)) {
	      DiffB_nDDsPi_Peak++;
	    }

	    if ((BDecMod31 >= 4111 && BDecMod31 <= 4116) ||
		(BDecMod31 >= 4221 && BDecMod31 <= 4226) ||
		(BDecMod31 >= 4121 && BDecMod31 <= 4126) ||
		(BDecMod31 >= 4211 && BDecMod31 <= 4216) ||
		(BDecMod31 >= 4331 && BDecMod31 <= 4336) ||
		(BDecMod31 >= 4441 && BDecMod31 <= 4446) ||
		(BDecMod31 >= 4131 && BDecMod31 <= 4136) ||
		(BDecMod31 >= 4311 && BDecMod31 <= 4316) ||
		(BDecMod31 >= 4141 && BDecMod31 <= 4146) ||
		(BDecMod31 >= 4411 && BDecMod31 <= 4416) ||
		(BDecMod31 >= 4231 && BDecMod31 <= 4236) ||
		(BDecMod31 >= 4321 && BDecMod31 <= 4326) ||
		(BDecMod31 >= 4241 && BDecMod31 <= 4246) ||
		(BDecMod31 >= 4421 && BDecMod31 <= 4426) ||
		(BDecMod31 >= 4341 && BDecMod31 <= 4346) ||
		(BDecMod31 >= 4431 && BDecMod31 <= 4436)) {
	      DiffB_nDDK_Peak++;
	    }	    
	    if ((BDecMod31 >= 4117 && BDecMod31 <= 4119) ||
		(BDecMod31 >= 4227 && BDecMod31 <= 4229) ||
		(BDecMod31 >= 4127 && BDecMod31 <= 4129) ||
		(BDecMod31 >= 4217 && BDecMod31 <= 4219) ||
		(BDecMod31 >= 4337 && BDecMod31 <= 4339) ||
		(BDecMod31 >= 4447 && BDecMod31 <= 4449) ||
		(BDecMod31 >= 4137 && BDecMod31 <= 4139) ||
		(BDecMod31 >= 4317 && BDecMod31 <= 4319) ||
		(BDecMod31 >= 4147 && BDecMod31 <= 4149) ||
		(BDecMod31 >= 4417 && BDecMod31 <= 4419) ||
		(BDecMod31 >= 4237 && BDecMod31 <= 4239) ||
		(BDecMod31 >= 4327 && BDecMod31 <= 4329) ||
		(BDecMod31 >= 4247 && BDecMod31 <= 4249) ||
		(BDecMod31 >= 4427 && BDecMod31 <= 4429) ||
		(BDecMod31 >= 4347 && BDecMod31 <= 4349) ||
		(BDecMod31 >= 4437 && BDecMod31 <= 4439)) {
	      DiffB_nDDPi_Peak++;
	    }	    
	    if ((BDecMod31 >= 4161 && BDecMod31 <= 4166) || 
		(BDecMod31 >= 4611 && BDecMod31 <= 4616) ||
		(BDecMod31 >= 4171 && BDecMod31 <= 4176) || 
		(BDecMod31 >= 4711 && BDecMod31 <= 4716) ||
		(BDecMod31 >= 4261 && BDecMod31 <= 4266) || 
		(BDecMod31 >= 4621 && BDecMod31 <= 4626) ||
		(BDecMod31 >= 4271 && BDecMod31 <= 4276) || 
		(BDecMod31 >= 4721 && BDecMod31 <= 4726) ||
		(BDecMod31 >= 4361 && BDecMod31 <= 4366) || 
		(BDecMod31 >= 4631 && BDecMod31 <= 4636) ||
		(BDecMod31 >= 4371 && BDecMod31 <= 4376) || 
		(BDecMod31 >= 4731 && BDecMod31 <= 4736) ||
		(BDecMod31 >= 4461 && BDecMod31 <= 4466) || 
		(BDecMod31 >= 4641 && BDecMod31 <= 4646) ||
		(BDecMod31 >= 4471 && BDecMod31 <= 4476) || 
		(BDecMod31 >= 4741 && BDecMod31 <= 4746)) {
	      DiffB_nDDsK_Peak++;
	    }
	    if ((BDecMod31 >= 4167 && BDecMod31 <= 4169) || 
		(BDecMod31 >= 4617 && BDecMod31 <= 4619) ||
		(BDecMod31 >= 4177 && BDecMod31 <= 4179) || 
		(BDecMod31 >= 4717 && BDecMod31 <= 4719) ||
		(BDecMod31 >= 4267 && BDecMod31 <= 4269) || 
		(BDecMod31 >= 4627 && BDecMod31 <= 4629) ||
		(BDecMod31 >= 4277 && BDecMod31 <= 4279) || 
		(BDecMod31 >= 4727 && BDecMod31 <= 4729) ||
		(BDecMod31 >= 4367 && BDecMod31 <= 4369) || 
		(BDecMod31 >= 4637 && BDecMod31 <= 4639) ||
		(BDecMod31 >= 4377 && BDecMod31 <= 4379) || 
		(BDecMod31 >= 4737 && BDecMod31 <= 4739) ||
		(BDecMod31 >= 4467 && BDecMod31 <= 4469) || 
		(BDecMod31 >= 4647 && BDecMod31 <= 4649) ||
		(BDecMod31 >= 4477 && BDecMod31 <= 4479) || 
		(BDecMod31 >= 4747 && BDecMod31 <= 4749)) {
	      DiffB_nDDsPi_Peak++;
	    }
	  }
	}
      }
    } // end of Dl loop
  } //end of event loop

  CascL_SameB_nDDX_Peak 
    = CascL_SameB_nDD_Peak + CascL_SameB_nDDs_Peak 
    + CascL_SameB_nDDK_Peak + CascL_SameB_nDDPi_Peak
    + CascL_SameB_nDDsK_Peak + CascL_SameB_nDDsPi_Peak;
  CascL_DiffB_nDDX_Peak 
    = CascL_DiffB_nDD_Peak + CascL_DiffB_nDDs_Peak 
    + CascL_DiffB_nDDK_Peak + CascL_DiffB_nDDPi_Peak
    + CascL_DiffB_nDDsK_Peak + CascL_DiffB_nDDsPi_Peak;
  DiffB_nDDX_Peak 
    = DiffB_nDD_Peak + DiffB_nDDs_Peak 
    + DiffB_nDDK_Peak + DiffB_nDDPi_Peak
    + DiffB_nDDsK_Peak + DiffB_nDDsPi_Peak;

  cout << "realD_Peak = " << realD_Peak << endl;
  cout << "realD_nDToKenu_Peak = " << realD_nDToKenu_Peak << endl;
  cout << "realD_nDToKstarenu_Peak = " << realD_nDToKstarenu_Peak << endl;
  cout << "CascL_SameB_Peak = " << CascL_SameB_Peak << endl;
  cout << "CascL_SameB_nDToKenu_Peak = " << CascL_SameB_nDToKenu_Peak << endl;
  cout << "CascL_SameB_nDToKstarenu_Peak = " << CascL_SameB_nDToKstarenu_Peak << endl;
  cout << "CascL_SameB_nDD_Peak = " << CascL_SameB_nDD_Peak << endl;
  cout << "CascL_SameB_nDDs_Peak = " << CascL_SameB_nDDs_Peak << endl;
  cout << "CascL_SameB_nDDK_Peak = " << CascL_SameB_nDDK_Peak << endl;
  cout << "CascL_SameB_nDDPi_Peak = " << CascL_SameB_nDDPi_Peak << endl;
  cout << "CascL_SameB_nDDsK_Peak = " << CascL_SameB_nDDsK_Peak << endl;
  cout << "CascL_SameB_nDDsPi_Peak = " << CascL_SameB_nDDsPi_Peak << endl;
  cout << "CascL_SameB_nDDX_Peak = " << CascL_SameB_nDDX_Peak << endl;
  cout << "CascL_DiffB_Peak = " << CascL_DiffB_Peak << endl;
  cout << "CascL_DiffB_nDToKenu_Peak = " << CascL_DiffB_nDToKenu_Peak << endl;
  cout << "CascL_DiffB_nDToKstarenu_Peak = " << CascL_DiffB_nDToKstarenu_Peak << endl;
  cout << "CascL_DiffB_nDD_Peak = " << CascL_DiffB_nDD_Peak << endl;
  cout << "CascL_DiffB_nDDs_Peak = " << CascL_DiffB_nDDs_Peak << endl;
  cout << "CascL_DiffB_nDDK_Peak = " << CascL_DiffB_nDDK_Peak << endl;
  cout << "CascL_DiffB_nDDPi_Peak = " << CascL_DiffB_nDDPi_Peak << endl;
  cout << "CascL_DiffB_nDDsK_Peak = " << CascL_DiffB_nDDsK_Peak << endl;
  cout << "CascL_DiffB_nDDsPi_Peak = " << CascL_DiffB_nDDsPi_Peak << endl;
  cout << "CascL_DiffB_nDDX_Peak = " << CascL_DiffB_nDDX_Peak << endl;
  cout << "DiffB_Peak = " << DiffB_Peak << endl;
  cout << "DiffB_nDD_Peak = " << DiffB_nDD_Peak << endl;
  cout << "DiffB_nDDs_Peak = " << DiffB_nDDs_Peak << endl;
  cout << "DiffB_nDDK_Peak = " << DiffB_nDDK_Peak << endl;
  cout << "DiffB_nDDPi_Peak = " << DiffB_nDDPi_Peak << endl;
  cout << "DiffB_nDDsK_Peak = " << DiffB_nDDsK_Peak << endl;
  cout << "DiffB_nDDsPi_Peak = " << DiffB_nDDsPi_Peak << endl;
  cout << "DiffB_nDDX_Peak = " << DiffB_nDDX_Peak << endl;

  // Open output file and save results. #####
  ofstream outFile(outFileName); //open output file
  if (!outFile) {
    cout << "Cannot open file " << outFile << endl;
  }
  outFile << "realD_Peak = " << realD_Peak << endl;
  outFile << "realD_nDToKenu_Peak = " << realD_nDToKenu_Peak << endl;
  outFile << "realD_nDToKstarenu_Peak = " << realD_nDToKstarenu_Peak << endl;
  outFile << "CascL_SameB_Peak = " << CascL_SameB_Peak << endl;
  outFile << "CascL_SameB_nDToKenu_Peak = " << CascL_SameB_nDToKenu_Peak << endl;
  outFile << "CascL_SameB_nDToKstarenu_Peak = " << CascL_SameB_nDToKstarenu_Peak << endl;
  outFile << "CascL_SameB_nDD_Peak = " << CascL_SameB_nDD_Peak << endl;
  outFile << "CascL_SameB_nDDs_Peak = " << CascL_SameB_nDDs_Peak << endl;
  outFile << "CascL_SameB_nDDK_Peak = " << CascL_SameB_nDDK_Peak << endl;
  outFile << "CascL_SameB_nDDPi_Peak = " << CascL_SameB_nDDPi_Peak << endl;
  outFile << "CascL_SameB_nDDsK_Peak = " << CascL_SameB_nDDsK_Peak << endl;
  outFile << "CascL_SameB_nDDsPi_Peak = " << CascL_SameB_nDDsPi_Peak << endl;
  outFile << "CascL_SameB_nDDX_Peak = " << CascL_SameB_nDDX_Peak << endl;
  outFile << "DiffB_Peak = " << DiffB_Peak << endl;
  outFile << "CascL_DiffB_Peak = " << CascL_DiffB_Peak << endl;
  outFile << "CascL_DiffB_nDToKenu_Peak = " << CascL_DiffB_nDToKenu_Peak << endl;
  outFile << "CascL_DiffB_nDToKstarenu_Peak = " << CascL_DiffB_nDToKstarenu_Peak << endl;
  outFile << "CascL_DiffB_nDD_Peak = " << CascL_DiffB_nDD_Peak << endl;
  outFile << "CascL_DiffB_nDDs_Peak = " << CascL_DiffB_nDDs_Peak << endl;
  outFile << "CascL_DiffB_nDDsK_Peak = " << CascL_DiffB_nDDsK_Peak << endl;
  outFile << "CascL_DiffB_nDDsPi_Peak = " << CascL_DiffB_nDDsPi_Peak << endl;
  outFile << "CascL_DiffB_nDDK_Peak = " << CascL_DiffB_nDDK_Peak << endl;
  outFile << "CascL_DiffB_nDDPi_Peak = " << CascL_DiffB_nDDPi_Peak << endl;
  outFile << "CascL_DiffB_nDDX_Peak = " << CascL_DiffB_nDDX_Peak << endl;
  outFile << "DiffB_nDD_Peak = " << DiffB_nDD_Peak << endl;
  outFile << "DiffB_nDDs_Peak = " << DiffB_nDDs_Peak << endl;
  outFile << "DiffB_nDDK_Peak = " << DiffB_nDDK_Peak << endl;
  outFile << "DiffB_nDDPi_Peak = " << DiffB_nDDPi_Peak << endl;
  outFile << "DiffB_nDDsK_Peak = " << DiffB_nDDsK_Peak << endl;
  outFile << "DiffB_nDDsPi_Peak = " << DiffB_nDDsPi_Peak << endl;
  outFile << "DiffB_nDDX_Peak = " << DiffB_nDDX_Peak << endl;
  outFile.close();


}