diff --git a/analyzers/dataframe/myUtils.cc b/analyzers/dataframe/myUtils.cc index a477265ce6..387a7e72e9 100644 --- a/analyzers/dataframe/myUtils.cc +++ b/analyzers/dataframe/myUtils.cc @@ -279,8 +279,9 @@ myUtils::get_VertexObject(ROOT::VecOps::RVec seltracks; seltracks.push_back(recop.at(pions.at(i))); seltracks.push_back(recop.at(kaons.at(j))); - VertexingUtils::FCCAnalysesVertex TheVertex = VertexFitterSimple::VertexFitter(0,seltracks, tracks ); + VertexingUtils::FCCAnalysesVertex TheVertex = VertexFitterSimple::VertexFitter(0,seltracks, tracks); float chi2 = TheVertex.vertex.chi2; if (chi2<0.01)continue; @@ -1630,7 +1631,7 @@ ROOT::VecOps::RVec myUtils::awkwardtest(ROOT::VecOps::RVec tk_ind; ROOT::VecOps::RVec seltracks = VertexingUtils::selTracks(0.,3.,0.,3.)( recop, tracks); - VertexingUtils::FCCAnalysesVertex ThePVertex = VertexFitterSimple::VertexFitter(0,seltracks, tracks ); + VertexingUtils::FCCAnalysesVertex ThePVertex = VertexFitterSimple::VertexFitter(0,seltracks, tracks); int PV_ntrk = ThePVertex.ntracks; float PV_chi2 = ThePVertex.vertex.chi2; @@ -1750,7 +1751,7 @@ ROOT::VecOps::RVec myUtils::awkwardtest(ROOT::VecOps::RVec bdt("Bc2TauNu_BDT", "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/xgb_bdt_vtx.root"); - computeModel = TMVA::Experimental::Compute<18, float>(bdt); ''') @@ -42,7 +41,7 @@ def __init__(self, inputlist, outname, ncpu): print (" init done, about to run") #__________________________________________________________ def run(self): - #df2 = (self.df.Range(10) + #df2 = (self.df.Range(100) df2 = (self.df ############################################# ## Aliases for # in python ## @@ -376,7 +375,7 @@ def run(self): df2.Snapshot("events", self.outname, branchList) # example call for standalone file -# python examples/FCCee/flavour/Bc2TauNu/analysis_stage1.py p8_ee_Zbb_Bc2TauNu_stage1.root /eos/experiment/fcc/ee/generation/DelphesEvents/spring2021/IDEA/p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU/events_003834121.root +# python examples/FCCee/flavour/Bc2TauNu/analysis_stage1.py p8_ee_Zbb_Bc2TauNu_stage1.root /eos/experiment/fcc/ee/generation/DelphesEvents/spring2021/IDEA/p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU/events_003907469.root # python examples/FCCee/flavour/Bc2TauNu/analysis_stage1.py p8_ee_Zbb_Bu2TauNu_stage1.root /eos/experiment/fcc/ee/generation/DelphesEvents/fcc_tmp_v03/p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU/events_026079857.root diff --git a/examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py b/examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py index 54781734ca..cc56a10ccc 100644 --- a/examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py +++ b/examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py @@ -22,7 +22,12 @@ ROOT.gInterpreter.ProcessLine(''' TMVA::Experimental::RBDT<> bdt("Bc2TauNu_BDT2", "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/xgb_bdt_stage2.root"); -computeModel = TMVA::Experimental::Compute<20, float>(bdt); +computeModel = TMVA::Experimental::Compute<21, float>(bdt); +''') + +ROOT.gInterpreter.ProcessLine(''' +TMVA::Experimental::RBDT<> bdt("Bc2TauNu_BDT", "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/xgb_bdt_vtx.root"); +computeModelBis = TMVA::Experimental::Compute<18, float>(bdt); ''') @@ -58,7 +63,7 @@ def run(self): .Define("LOCAL_CandVtxInd", "Tau23PiCandidates_vertex.at(CUT_CandInd)") ##LOCAL INDEX screwed up in prod02!!!! need -1 because PV is removed - .Define("CUT_CandVtxThrustEmin", "Vertex_thrusthemis_emin.at(LOCAL_CandVtxInd-1)") + .Define("CUT_CandVtxThrustEmin", "Vertex_thrusthemis_emin.at(LOCAL_CandVtxInd)") .Filter("CUT_CandVtxThrustEmin>0") .Define("EVT_CandN", "float(Tau23PiCandidates_vertex.size())") @@ -103,12 +108,30 @@ def run(self): .Define("EVT_DVmass_ave", "myFinalSel::get_ave(Vertex_mass, Vertex_isPV, LOCAL_CandVtxInd)") .Define("EVT_PVmass", "Vertex_mass.at(0)") - .Define("MVAVec", ROOT.computeModel, ("EVT_CandMass","EVT_CandRho1Mass","EVT_CandRho2Mass","EVT_CandN","EVT_CandVtxFD","EVT_CandVtxChi2","EVT_CandPx","EVT_CandPy","EVT_CandPz","EVT_CandP","EVT_CandD0","EVT_CandZ0","EVT_CandAngleThrust","EVT_DVd0_min","EVT_DVd0_max","EVT_DVd0_ave","EVT_DVz0_min","EVT_DVz0_max","EVT_DVz0_ave","EVT_Nominal_B_E")) - + .Define("MVAVec", ROOT.computeModel, ("EVT_CandMass", "EVT_CandRho1Mass", "EVT_CandRho2Mass", + "EVT_CandN", "EVT_CandVtxFD", "EVT_CandVtxChi2", + "EVT_CandPx", "EVT_CandPy", "EVT_CandPz", + "EVT_CandP", "EVT_CandD0", "EVT_CandZ0", + "EVT_CandAngleThrust", "EVT_DVd0_min", "EVT_DVd0_max", + "EVT_DVd0_ave", "EVT_DVz0_min", "EVT_DVz0_max", + "EVT_DVz0_ave", "EVT_PVmass", "EVT_Nominal_B_E")) + + .Define("EVT_MVA2", "MVAVec.at(0)") .Filter(MVAFilter) - + .Define("MVAVecBis", ROOT.computeModelBis, ("EVT_ThrustEmin_E", "EVT_ThrustEmax_E", + "EVT_ThrustEmin_Echarged", "EVT_ThrustEmax_Echarged", + "EVT_ThrustEmin_Eneutral", "EVT_ThrustEmax_Eneutral", + "EVT_ThrustEmin_Ncharged", "EVT_ThrustEmax_Ncharged", + "EVT_ThrustEmin_Nneutral", "EVT_ThrustEmax_Nneutral", + "EVT_NtracksPV", "EVT_NVertex", + "EVT_NTau23Pi", "EVT_ThrustEmin_NDV", + "EVT_ThrustEmax_NDV", "EVT_dPV2DVmin", + "EVT_dPV2DVmax", "EVT_dPV2DVave")) + .Define("EVT_MVA1Bis", "MVAVecBis.at(0)") + + .Define("EVT_minRhoMass", "if (EVT_CandRho1MassEVT_CandRho2Mass) return EVT_CandRho1Mass; else return EVT_CandRho2Mass;") .Define("EVT_ThrustDiff_E", "EVT_ThrustEmax_E-EVT_ThrustEmin_E") @@ -156,7 +179,7 @@ def run(self): "EVT_CandPx","EVT_CandP","EVT_CandPz","EVT_CandPy", "EVT_CandD0","EVT_CandZ0","EVT_CandAngleThrust", "EVT_minRhoMass", "EVT_maxRhoMass", - "EVT_MVA1", "EVT_MVA2", + "EVT_MVA1", "EVT_MVA2","EVT_MVA1Bis", "EVT_CandPion1P","EVT_CandPion1D0","EVT_CandPion1Z0", "EVT_CandPion2P","EVT_CandPion2D0","EVT_CandPion2Z0", @@ -181,28 +204,11 @@ def run(self): #df2.Snapshot("events", self.outname, branchList, opts) df2.Snapshot("events", self.outname, branchList) -# example call for standalone file -# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/21042021/Batch_Analysis_stage1/p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU/*.root" - -# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/21042021/Batch_Analysis_stage1/p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU/*.root" - -# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/21042021/Batch_Analysis_stage1/p8_ee_Zbb_ecm91/*.root" - -# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zcc_ecm91.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/21042021/Batch_Analysis_stage1/p8_ee_Zcc_ecm91/*.root" - -# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zuds_ecm91.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/21042021/Batch_Analysis_stage1/p8_ee_Zuds_ecm91/*.root" ######################################### -# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Batch_Training_4stage1/p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU/*.root" - -# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Batch_Training_4stage1/p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU/*.root" - -# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Batch_Training_4stage1/p8_ee_Zbb_ecm91/*.root" - -# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zcc_ecm91.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Batch_Training_4stage1/p8_ee_Zcc_ecm91/*.root" +# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU.root /eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_03/Batch_Training_4stage1/p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU/flat_chunk_0.root -# python examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py p8_ee_Zuds_ecm91.root "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Batch_Training_4stage1/p8_ee_Zuds_ecm91/*.root" if __name__ == "__main__": diff --git a/examples/FCCee/flavour/Bc2TauNu/finalSel.py b/examples/FCCee/flavour/Bc2TauNu/finalSel.py index 0c08729197..a040879d76 100644 --- a/examples/FCCee/flavour/Bc2TauNu/finalSel.py +++ b/examples/FCCee/flavour/Bc2TauNu/finalSel.py @@ -5,18 +5,15 @@ import ROOT ###Input directory where the files produced at the pre-selection level are -baseDir = "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Analysis_stage2/" +baseDir = "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_04/Analysis_stage2/" ###Link to the dictonary that contains all the cross section informations etc... procDict = os.path.join(os.getenv('FCCDICTSDIR', deffccdicts), '') + "FCCee_procDict_spring2021_IDEA.json" -#baseDir = "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/21042021/Training_4stage2/" -#procDict = os.path.join(os.getenv('FCCDICTSDIR', deffccdicts), '') + "FCCee_procDict_fcc_tmp_training.json" - - process_list=['p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU', 'p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU', + 'p8_ee_Zbb_ecm91_EvtGen', 'p8_ee_Zbb_ecm91', 'p8_ee_Zcc_ecm91', 'p8_ee_Zuds_ecm91', @@ -54,7 +51,6 @@ 'p8_ee_Zbb_ecm91_EvtGen_Lb2LcstTauNu', ] -#process_list=['p8_ee_Zbb_ecm91_EvtGen_Bs2DsstDsst'] define_list={ # "EVT_minRhoMass":"if (EVT_CandRho1Mass0.8 && EVT_MVA2>0.8", +cut_list = {#"sel0":"EVT_MVA1>0.8 && EVT_MVA2>0.8", #"sel1":"EVT_MVA1>0.9 && EVT_MVA2>0.9", #"sel2":"EVT_MVA1>0.95 && EVT_MVA2>0.95", #"sel3":"EVT_MVA1>0.98 && EVT_MVA2>0.98", #"sel4":"EVT_MVA1>0.99 && EVT_MVA2>0.99", + "sel5":"EVT_MVA1>0.9998979591836735 && EVT_MVA2>0.9956122448979592 && EVT_ThrustDiff_E > 10", } + ###Dictionary for the ouput variable/hitograms. The key is the name of the variable in the output files. "name" is the name of the variable in the input file, "title" is the x-axis label of the histogram, "bin" the number of bins of the histogram, "xmin" the minimum x-axis value and "xmax" the maximum x-axis value. variables = { @@ -97,11 +95,23 @@ "EVT_CandAngleThrust_2" :{"name":"EVT_CandAngleThrust","title":"Angle between candidate and thrust","bin":50,"xmin":0.,"xmax":0.5}, "EVT_MVA1" :{"name":"EVT_MVA1","title":"MVA1 score","bin":100,"xmin":0.6,"xmax":1.}, "EVT_MVA2" :{"name":"EVT_MVA2","title":"MVA2 score","bin":100,"xmin":0.6,"xmax":1.}, - + "EVT_Nominal_B_E" :{"name":"EVT_Nominal_B_E","title":"Nominal B Energy [GeV]","bin":100,"xmin":0.,"xmax":50.}, + "EVT_PVmass" :{"name":"EVT_PVmass","title":"Primary Vertex Mass [GeV]","bin":100,"xmin":0.,"xmax":10.}, + "EVT_DVmass_min" :{"name":"EVT_DVmass_min","title":"Displaced Vertex Mass Min [GeV]","bin":100,"xmin":0.,"xmax":20.}, + "EVT_DVmass_max" :{"name":"EVT_DVmass_max","title":"Displaced Vertex Mass Max [GeV]","bin":100,"xmin":0.,"xmax":20.}, + "EVT_DVmass_ave" :{"name":"EVT_DVmass_ave","title":"Displaced Vertex Mass Ave [GeV]","bin":100,"xmin":0.,"xmax":20.}, + "EVT_DVd0_min" :{"name":"EVT_DVd0_min","title":"Displaced Vertex d0 Min [mm]","bin":100,"xmin":-10.,"xmax":10.}, + "EVT_DVd0_max" :{"name":"EVT_DVd0_max","title":"Displaced Vertex d0 Max [mm]","bin":100,"xmin":-10.,"xmax":10.}, + "EVT_DVd0_ave" :{"name":"EVT_DVd0_ave","title":"Displaced Vertex d0 Ave [mm]","bin":100,"xmin":-10.,"xmax":10.}, + "EVT_DVz0_min" :{"name":"EVT_DVz0_min","title":"Displaced Vertex z0 Min [mm]","bin":100,"xmin":-10.,"xmax":10.}, + "EVT_DVz0_max" :{"name":"EVT_DVz0_max","title":"Displaced Vertex z0 Max [mm]","bin":100,"xmin":-10.,"xmax":10.}, + "EVT_DVz0_ave" :{"name":"EVT_DVz0_ave","title":"Displaced Vertex z0 Ave [mm]","bin":100,"xmin":-10.,"xmax":10.}, + + } ###Number of CPUs to use -NUM_CPUS = 4 +NUM_CPUS = 8 ###This part is standard to all analyses import config.runDataFrameFinal as rdf diff --git a/examples/FCCee/flavour/Bc2TauNu/plots.py b/examples/FCCee/flavour/Bc2TauNu/plots.py index e4ad83959c..82195371d1 100644 --- a/examples/FCCee/flavour/Bc2TauNu/plots.py +++ b/examples/FCCee/flavour/Bc2TauNu/plots.py @@ -6,11 +6,11 @@ delphesVersion = "3.4.2" energy = 91.0 collider = "FCC-ee" -inputDir = "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_01/Analysis_stage2/" +inputDir = "/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_04/Analysis_stage2/" formats = ['png','pdf'] yaxis = ['lin','log'] stacksig = ['nostack'] -outdir = '/eos/home-h/helsens/www/FCC/ee/flavour/Bc2TauNu/spring2021/prod_01/' +outdir = '/eos/home-h/helsens/www/FCC/ee/flavour/Bc2TauNu/spring2021/prod_04/plots/' variables = [ "EVT_CandMass", "EVT_ThrustEmin_E","EVT_ThrustEmax_E", "EVT_ThrustEmin_Echarged", "EVT_ThrustEmax_Echarged", @@ -25,8 +25,21 @@ "EVT_CandVtxFD", "EVT_CandAngleThrust", "EVT_CandAngleThrust_2", - "EVT_MVA1","EVT_MVA2" + "EVT_MVA1","EVT_MVA2", + + "EVT_Nominal_B_E", + "EVT_PVmass", + "EVT_DVmass_min", + "EVT_DVmass_max", + "EVT_DVmass_ave", + "EVT_DVd0_min", + "EVT_DVd0_max", + "EVT_DVd0_ave", + "EVT_DVz0_min", + "EVT_DVz0_max", + "EVT_DVz0_ave", + ] scaleSig=1. @@ -34,8 +47,10 @@ selections = {} selections['Inclusive'] = ["sel1","sel2","sel3","sel4"] selections['ExclusiveMerged'] = ["sel1","sel2","sel3","sel4"] -selections['Inclusive'] = ["sel0"] -selections['ExclusiveMerged'] = ["sel0"] + +selections['Inclusive'] = ["sel5"] +selections['InclusiveEvtGen'] = ["sel5"] +selections['ExclusiveMerged'] = ["sel5"] extralabel = {} extralabel['sel0'] = "Selection: MVA1>0.8, MVA2>0.8" @@ -43,10 +58,12 @@ extralabel['sel2'] = "Selection: MVA1>0.95, MVA2>0.95" extralabel['sel3'] = "Selection: MVA1>0.98, MVA2>0.98" extralabel['sel4'] = "Selection: MVA1>0.99, MVA2>0.99" +extralabel['sel5'] = "Selection: Donal's cuts" colors = {} colors['Z_Bc'] = ROOT.kRed colors['Z_bb'] = ROOT.kBlue +colors['Z_bb_EvtGen'] = ROOT.kBlue colors['Z_bb_B'] = ROOT.kBlue colors['Z_bb_Bu'] = ROOT.kBlue+1 colors['Z_bb_Bd'] = ROOT.kBlue+2 @@ -65,6 +82,13 @@ } } +plots['InclusiveEvtGen'] = {'signal':{'Z_Bc':['p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU']}, + 'backgrounds':{'Z_bb_EvtGen':['p8_ee_Zbb_ecm91_EvtGen'], + 'Z_cc':['p8_ee_Zcc_ecm91'], + 'Z_uds':['p8_ee_Zuds_ecm91'], + 'Z_Bu':['p8_ee_Zbb_ecm91_EvtGen_Bu2TauNuTAUHADNU'] + } + } plots['ExclusiveMerged'] = {'signal':{'Z_Bc':['p8_ee_Zbb_ecm91_EvtGen_Bc2TauNuTAUHADNU']}, @@ -109,6 +133,7 @@ legend = {} legend['Z_Bc'] = 'B_{c}#rightarrow #tau#nu (#tau#rightarrow 3#pi)' legend['Z_bb'] = 'Z#rightarrow b#bar{b}' +legend['Z_bb_EvtGen'] = 'Z#rightarrow b#bar{b} EvtGen' legend['Z_bb_B'] = 'Z#rightarrow b#bar{b} B_{u}+B_{d}+B_{s}+L_{b}' legend['Z_cc'] = 'Z#rightarrow c#bar{c}' legend['Z_uds'] = 'Z#rightarrow q#bar{q}' diff --git a/examples/FCCee/flavour/Bc2TauNu/preSel_Batch.py b/examples/FCCee/flavour/Bc2TauNu/preSel_Batch.py index 29c3098877..cbd8b9382f 100644 --- a/examples/FCCee/flavour/Bc2TauNu/preSel_Batch.py +++ b/examples/FCCee/flavour/Bc2TauNu/preSel_Batch.py @@ -5,7 +5,7 @@ import os basedir=os.path.join(os.getenv('FCCDICTSDIR', deffccdicts), '') + "yaml/FCCee/spring2021/IDEA/" -outdir="/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_02/Batch_Analysis_stage1/" +outdir="/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_04/Batch_Analysis_stage1/" NUM_CPUS=8 output_list=[] @@ -13,7 +13,8 @@ inputana="analysis_stage1.py" -process_list=['p8_ee_Zbb_ecm91', +process_list=['p8_ee_Zbb_ecm91_EvtGen', + 'p8_ee_Zbb_ecm91', 'p8_ee_Zcc_ecm91', 'p8_ee_Zuds_ecm91', ] diff --git a/examples/FCCee/flavour/Bc2TauNu/preSel_BatchTraining.py b/examples/FCCee/flavour/Bc2TauNu/preSel_BatchTraining.py index 49c1f92966..662bfab2d1 100644 --- a/examples/FCCee/flavour/Bc2TauNu/preSel_BatchTraining.py +++ b/examples/FCCee/flavour/Bc2TauNu/preSel_BatchTraining.py @@ -11,7 +11,7 @@ ##To produce variables for MVA stage 1 training inputana="analysis_stage1.py" -outdir="/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_02/Batch_Training_4stage1/" +outdir="/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_03/Batch_Training_4stage1/" fraction=1. diff --git a/examples/FCCee/flavour/Bc2TauNu/send_bg_stage2.sh b/examples/FCCee/flavour/Bc2TauNu/send_bg_stage2.sh index bb3232e8b8..eadd9bfe61 100755 --- a/examples/FCCee/flavour/Bc2TauNu/send_bg_stage2.sh +++ b/examples/FCCee/flavour/Bc2TauNu/send_bg_stage2.sh @@ -1,13 +1,13 @@ # example call for standalone file -input=/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_02/Batch_Analysis_stage1 -output=/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_02/Analysis_stage2 +input=/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_04/Batch_Analysis_stage1 +output=/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_04/Analysis_stage2 ana=examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py -#for process in p8_ee_Zbb_ecm91 p8_ee_Zcc_ecm91 p8_ee_Zuds_ecm91; -for process in p8_ee_Zcc_ecm91; +for process in p8_ee_Zbb_ecm91 p8_ee_Zbb_ecm91_EvtGen p8_ee_Zcc_ecm91 p8_ee_Zuds_ecm91 ; +#for process in p8_ee_Zcc_ecm91; do python $ana $output/$process.root "$input/$process/*.root" done diff --git a/examples/FCCee/flavour/Bc2TauNu/send_bgexdec_stage2.sh b/examples/FCCee/flavour/Bc2TauNu/send_bgexdec_stage2.sh index 24ebc54362..f935eb2515 100755 --- a/examples/FCCee/flavour/Bc2TauNu/send_bgexdec_stage2.sh +++ b/examples/FCCee/flavour/Bc2TauNu/send_bgexdec_stage2.sh @@ -1,14 +1,14 @@ # example call for standalone file -input=/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_02/Batch_Analysis_stage1 -output=/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_02/Analysis_stage2 +input=/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_03/Batch_Analysis_stage1 +output=/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_03/Analysis_stage2 ana=examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py -#for process in p8_ee_Zbb_ecm91_EvtGen_Bd2D3Pi p8_ee_Zbb_ecm91_EvtGen_Bd2DDs p8_ee_Zbb_ecm91_EvtGen_Bd2DTauNu p8_ee_Zbb_ecm91_EvtGen_Bd2Dst3Pi p8_ee_Zbb_ecm91_EvtGen_Bd2DstDs p8_ee_Zbb_ecm91_EvtGen_Bd2DstDsst p8_ee_Zbb_ecm91_EvtGen_Bd2DstTauNu; -#do -# python $ana $output/$process.root "$input/$process/*.root" -#done +for process in p8_ee_Zbb_ecm91_EvtGen_Bd2D3Pi p8_ee_Zbb_ecm91_EvtGen_Bd2DDs p8_ee_Zbb_ecm91_EvtGen_Bd2DTauNu p8_ee_Zbb_ecm91_EvtGen_Bd2Dst3Pi p8_ee_Zbb_ecm91_EvtGen_Bd2DstDs p8_ee_Zbb_ecm91_EvtGen_Bd2DstDsst p8_ee_Zbb_ecm91_EvtGen_Bd2DstTauNu; +do + python $ana $output/$process.root "$input/$process/*.root" +done #for process in p8_ee_Zbb_ecm91_EvtGen_Bs2Ds3Pi p8_ee_Zbb_ecm91_EvtGen_Bs2DsDs p8_ee_Zbb_ecm91_EvtGen_Bs2DsTauNu p8_ee_Zbb_ecm91_EvtGen_Bs2Dsst3Pi p8_ee_Zbb_ecm91_EvtGen_Bs2DsstDs p8_ee_Zbb_ecm91_EvtGen_Bs2DsstDsst p8_ee_Zbb_ecm91_EvtGen_Bs2DsstTauNu; @@ -23,7 +23,7 @@ ana=examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py #done -for process in p8_ee_Zbb_ecm91_EvtGen_Lb2Lc3Pi p8_ee_Zbb_ecm91_EvtGen_Lb2LcDs p8_ee_Zbb_ecm91_EvtGen_Lb2LcTauNu p8_ee_Zbb_ecm91_EvtGen_Lb2Lcst3Pi p8_ee_Zbb_ecm91_EvtGen_Lb2LcstDs p8_ee_Zbb_ecm91_EvtGen_Lb2LcstDsst p8_ee_Zbb_ecm91_EvtGen_Lb2LcstTauNu; -do - python $ana $output/$process.root "$input/$process/*.root" -done +#for process in p8_ee_Zbb_ecm91_EvtGen_Lb2Lc3Pi p8_ee_Zbb_ecm91_EvtGen_Lb2LcDs p8_ee_Zbb_ecm91_EvtGen_Lb2LcTauNu p8_ee_Zbb_ecm91_EvtGen_Lb2Lcst3Pi p8_ee_Zbb_ecm91_EvtGen_Lb2LcstDs p8_ee_Zbb_ecm91_EvtGen_Lb2LcstDsst p8_ee_Zbb_ecm91_EvtGen_Lb2LcstTauNu; +#do +# python $ana $output/$process.root "$input/$process/*.root" +#done diff --git a/examples/FCCee/flavour/Bc2TauNu/send_signal_stage2.sh b/examples/FCCee/flavour/Bc2TauNu/send_signal_stage2.sh index 6f60d8257b..3ce707c78d 100755 --- a/examples/FCCee/flavour/Bc2TauNu/send_signal_stage2.sh +++ b/examples/FCCee/flavour/Bc2TauNu/send_signal_stage2.sh @@ -1,6 +1,6 @@ # example call for standalone file -input=/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_02/Batch_Analysis_stage1 -output=/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_02/Analysis_stage2 +input=/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_04/Batch_Analysis_stage1 +output=/eos/experiment/fcc/ee/analyses/case-studies/flavour/Bc2TauNu/flatNtuples/spring2021/prod_04/Analysis_stage2 ana=examples/FCCee/flavour/Bc2TauNu/analysis_stage2.py