#include "SkimTuple.hh" using namespace std; typedef int Int_t; TStopwatch timer; int main() { timer.Start(); cout << "Started timer" << endl; inpfile >> run_range_min >> run_range_max >> nevmax >> ChainFileName; cout << "Run range used: " << run_range_min << " " << run_range_max << endl; cout << "Chain file = " << ChainFileName << endl; TChain *ntp100 = new TChain("ntp100"); fin = new ifstream(ChainFileName); cout << "Setup chain" << endl; // Setup file to hold runNumber and number of events // Do this early so that if it's going to crash, it does so soon! char cfn[250]; cout << "Setup Chain" << endl; while(fin->getline(cfn, 250)){ ntp100->Add(cfn); //cout << "Setup Chain" << endl; } cout << "Setup Chain" << endl; cout << "Files added to chain" << endl; ntp100->GetEntry(0); cout << "Found first entry" << endl; //Set Branch addresses cout << "Setting Branch Addresses" << endl; ntp100->SetBranchAddress("rndmNum",&rndmNum); ntp100->SetBranchAddress("DstD0pi_num",&DstD0pi_num); ntp100->SetBranchAddress("DKKpi_num",&DKKpi_num); ntp100->SetBranchAddress("DKpipi_num",&DKpipi_num); ntp100->SetBranchAddress("DKpipi_unbias_num",&DKpipi_unbias_num); ntp100->SetBranchAddress("runNumber",&runNumber); ntp100->SetBranchAddress("platform",&platform); //ntp100->SetBranchAddress("partition",&partition); ntp100->SetBranchAddress("upperID",&upperID); ntp100->SetBranchAddress("lowerID",&lowerID); ntp100->SetBranchAddress("energyCM",&energyCM); ntp100->SetBranchAddress("ePlusMomentumX",&ePlusMomentumX); ntp100->SetBranchAddress("ePlusMomentumY",&ePlusMomentumY); ntp100->SetBranchAddress("ePlusMomentumZ",&ePlusMomentumZ); ntp100->SetBranchAddress("eMinusMomentumX",&eMinusMomentumX); ntp100->SetBranchAddress("eMinusMomentumY",&eMinusMomentumY); ntp100->SetBranchAddress("eMinusMomentumZ",&eMinusMomentumZ); ntp100->SetBranchAddress("boostCMtoLabX",&boostCMtoLabX); ntp100->SetBranchAddress("boostCMtoLabY",&boostCMtoLabY); ntp100->SetBranchAddress("boostCMtoLabZ",&boostCMtoLabZ); ntp100->SetBranchAddress("IPX",&IPX); ntp100->SetBranchAddress("IPY",&IPY); ntp100->SetBranchAddress("IPZ",&IPZ); ntp100->SetBranchAddress("xxpv",&xxpv); ntp100->SetBranchAddress("xypv",&xypv); ntp100->SetBranchAddress("xzpv",&xzpv); ntp100->SetBranchAddress("yypv",&yypv); ntp100->SetBranchAddress("yzpv",&yzpv); ntp100->SetBranchAddress("zzpv",&zzpv); ntp100->SetBranchAddress("ThrustMagCTrk",&ThrustMagCTrk); ntp100->SetBranchAddress("ThrustMagAll",&ThrustMagAll); ntp100->SetBranchAddress("ThrustCosThCTrk",&ThrustCosThCTrk); ntp100->SetBranchAddress("ThrustCosThAll",&ThrustCosThAll); ntp100->SetBranchAddress("ThrustPhiCTrk",&ThrustPhiCTrk); ntp100->SetBranchAddress("ThrustPhiAll",&ThrustPhiAll); ntp100->SetBranchAddress("R2CTrk",&R2CTrk); ntp100->SetBranchAddress("R2All",&R2All); ntp100->SetBranchAddress("SphericityAll",&SphericityAll); ntp100->SetBranchAddress("xPrimaryVtx",&xPrimaryVtx); ntp100->SetBranchAddress("yPrimaryVtx",&yPrimaryVtx); ntp100->SetBranchAddress("zPrimaryVtx",&zPrimaryVtx); ntp100->SetBranchAddress("probPrimaryVtx",&probPrimaryVtx); ntp100->SetBranchAddress("mc_num",&mc_num); ntp100->SetBranchAddress("trk_num",&trk_num); ntp100->SetBranchAddress("trk_GenWord",trk_GenWord); ntp100->SetBranchAddress("trk_PIDLongWord",trk_PIDLongWord); ntp100->SetBranchAddress("trk_eLHbit",trk_eLHbit); ntp100->SetBranchAddress("trk_charge",trk_charge); ntp100->SetBranchAddress("trk_px",trk_px); ntp100->SetBranchAddress("trk_py",trk_py); ntp100->SetBranchAddress("trk_pz",trk_pz); ntp100->SetBranchAddress("trk_p4ErrCovXX",trk_p4ErrCovXX); ntp100->SetBranchAddress("trk_p4ErrCovXY",trk_p4ErrCovXY); ntp100->SetBranchAddress("trk_p4ErrCovXZ",trk_p4ErrCovXZ); ntp100->SetBranchAddress("trk_p4ErrCovYY",trk_p4ErrCovYY); ntp100->SetBranchAddress("trk_p4ErrCovYZ",trk_p4ErrCovYZ); ntp100->SetBranchAddress("trk_p4ErrCovZZ",trk_p4ErrCovZZ); ntp100->SetBranchAddress("trk_ndof",trk_ndof); ntp100->SetBranchAddress("trk_chisq",trk_chisq); ntp100->SetBranchAddress("trk_svtpat",trk_svtpat); ntp100->SetBranchAddress("trk_phiatemc",trk_phiatemc); ntp100->SetBranchAddress("trk_dedxsvt",trk_dedxsvt); ntp100->SetBranchAddress("trk_nsvt_dedx",trk_nsvt_dedx); ntp100->SetBranchAddress("trk_nsvt_mom",trk_nsvt_mom); ntp100->SetBranchAddress("trk_dedxdch",trk_dedxdch); ntp100->SetBranchAddress("trk_ndch_dedx",trk_ndch_dedx); ntp100->SetBranchAddress("trk_ndch_mom",trk_ndch_mom); ntp100->SetBranchAddress("trk_tht",trk_tht); ntp100->SetBranchAddress("trk_thterr",trk_thterr); ntp100->SetBranchAddress("trk_nph",trk_nph); ntp100->SetBranchAddress("trk_nphbkg",trk_nphbkg); ntp100->SetBranchAddress("trk_phicluster",trk_phicluster); ntp100->SetBranchAddress("trk_ecal",trk_ecal); ntp100->SetBranchAddress("trk_lat",trk_lat); ntp100->SetBranchAddress("trk_z42",trk_z42); ntp100->SetBranchAddress("trk_z20",trk_z20); ntp100->SetBranchAddress("trk_ncrys",trk_ncrys); ntp100->SetBranchAddress("trk_s1s9",trk_s1s9); ntp100->SetBranchAddress("trk_s9s25",trk_s9s25); ntp100->SetBranchAddress("trk_ifrexp",trk_ifrexp); ntp100->SetBranchAddress("trk_ifrmeas",trk_ifrmeas); ntp100->SetBranchAddress("trk_ifrbefore",trk_ifrbefore); ntp100->SetBranchAddress("trk_ifrchisq",trk_ifrchisq); ntp100->SetBranchAddress("DKKpi_pipx",DKKpi_pipx); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_pipy",DKKpi_pipy); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_pipz",DKKpi_pipz); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_K1px",DKKpi_K1px); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_K1py",DKKpi_K1py); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_K1pz",DKKpi_K1pz); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_K2px",DKKpi_K2px); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_K2py",DKKpi_K2py); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_K2pz", DKKpi_K2pz); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_DVtxPosx",DKKpi_DVtxPosx); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_DVtxPosy",DKKpi_DVtxPosy); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_DVtxPosz",DKKpi_DVtxPosz); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_DChiSquared",DKKpi_DChiSquared); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_DnDof", DKKpi_DnDof); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_piIdx", DKKpi_piIdx); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_K1Idx",DKKpi_K1Idx); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_K2Idx",DKKpi_K2Idx); //[DKKpi_num] ntp100->SetBranchAddress("DKKpi_DVtxSepChiSq", DKKpi_DVtxSepChiSq); //[DKKpi_num] ntp100->SetBranchAddress("DKpipi_Kpx",DKpipi_Kpx); ntp100->SetBranchAddress("DKpipi_Kpy",DKpipi_Kpy); ntp100->SetBranchAddress("DKpipi_Kpz",DKpipi_Kpz); ntp100->SetBranchAddress("DKpipi_pi1px",DKpipi_pi1px); ntp100->SetBranchAddress("DKpipi_pi1py",DKpipi_pi1py); ntp100->SetBranchAddress("DKpipi_pi1pz",DKpipi_pi1pz); ntp100->SetBranchAddress("DKpipi_pi2px",DKpipi_pi2px); ntp100->SetBranchAddress("DKpipi_pi2py",DKpipi_pi2py); ntp100->SetBranchAddress("DKpipi_pi2pz",DKpipi_pi2pz); ntp100->SetBranchAddress("DKpipi_DVtxPosx",DKpipi_DVtxPosx); ntp100->SetBranchAddress("DKpipi_DVtxPosy",DKpipi_DVtxPosy); ntp100->SetBranchAddress("DKpipi_DVtxPosz",DKpipi_DVtxPosz); ntp100->SetBranchAddress("DKpipi_DChiSquared",DKpipi_DChiSquared); ntp100->SetBranchAddress("DKpipi_DnDof",DKpipi_DnDof); ntp100->SetBranchAddress("DKpipi_KIdx",DKpipi_KIdx); ntp100->SetBranchAddress("DKpipi_pi1Idx",DKpipi_pi1Idx); ntp100->SetBranchAddress("DKpipi_pi2Idx",DKpipi_pi2Idx); ntp100->SetBranchAddress("DKpipi_DVtxSepChiSq", DKpipi_DVtxSepChiSq); ntp100->SetBranchAddress("DKpipi_KTrkpx",DKpipi_KTrkpx); ntp100->SetBranchAddress("DKpipi_KTrkpy",DKpipi_KTrkpy); ntp100->SetBranchAddress("DKpipi_KTrkpz",DKpipi_KTrkpz); ntp100->SetBranchAddress("DKpipi_pi1Trkpx",DKpipi_pi1Trkpx); ntp100->SetBranchAddress("DKpipi_pi1Trkpy",DKpipi_pi1Trkpy); ntp100->SetBranchAddress("DKpipi_pi1Trkpz",DKpipi_pi1Trkpz); ntp100->SetBranchAddress("DKpipi_pi2Trkpx",DKpipi_pi2Trkpx); ntp100->SetBranchAddress("DKpipi_pi2Trkpy",DKpipi_pi2Trkpy); ntp100->SetBranchAddress("DKpipi_pi2Trkpz",DKpipi_pi2Trkpz); ntp100->SetBranchAddress("DKpipi_DTrkVtxPosx",DKpipi_DTrkVtxPosx); ntp100->SetBranchAddress("DKpipi_DTrkVtxPosy",DKpipi_DTrkVtxPosy); ntp100->SetBranchAddress("DKpipi_DTrkVtxPosz",DKpipi_DTrkVtxPosz); ntp100->SetBranchAddress("DKpipi_DTrkChiSquared",DKpipi_DTrkChiSquared); ntp100->SetBranchAddress("DKpipi_DTrknDof",DKpipi_DTrknDof); ntp100->SetBranchAddress("DKpipi_KTrkIdx",DKpipi_KTrkIdx); ntp100->SetBranchAddress("DKpipi_pi1TrkIdx",DKpipi_pi1TrkIdx); ntp100->SetBranchAddress("DKpipi_pi2TrkIdx",DKpipi_pi2TrkIdx); ntp100->SetBranchAddress("DKpipi_DTrkVtxSepChiSq", DKpipi_DTrkVtxSepChiSq); ntp100->SetBranchAddress("D0Kpi_Kpx",D0Kpi_Kpx); ntp100->SetBranchAddress("D0Kpi_Kpy",D0Kpi_Kpy); ntp100->SetBranchAddress("D0Kpi_Kpz",D0Kpi_Kpz); ntp100->SetBranchAddress("D0Kpi_pipx",D0Kpi_pipx); ntp100->SetBranchAddress("D0Kpi_pipy",D0Kpi_pipy); ntp100->SetBranchAddress("D0Kpi_pipz",D0Kpi_pipz); ntp100->SetBranchAddress("D0Kpi_D0VtxPosx",D0Kpi_D0VtxPosx); ntp100->SetBranchAddress("D0Kpi_D0VtxPosy",D0Kpi_D0VtxPosy); ntp100->SetBranchAddress("D0Kpi_D0VtxPosz",D0Kpi_D0VtxPosz); ntp100->SetBranchAddress("D0Kpi_D0ChiSquared",D0Kpi_D0ChiSquared); ntp100->SetBranchAddress("D0Kpi_D0nDof",D0Kpi_D0nDof); ntp100->SetBranchAddress("D0Kpi_KIdx",D0Kpi_KIdx); ntp100->SetBranchAddress("D0Kpi_piIdx",D0Kpi_piIdx); ntp100->SetBranchAddress("DstD0pi_DstVtxPosx",DstD0pi_DstVtxPosx); ntp100->SetBranchAddress("DstD0pi_DstVtxPosy",DstD0pi_DstVtxPosy); ntp100->SetBranchAddress("DstD0pi_DstVtxPosz",DstD0pi_DstVtxPosz); ntp100->SetBranchAddress("DstD0pi_DstChiSquared",DstD0pi_DstChiSquared); ntp100->SetBranchAddress("DstD0pi_DstnDof",DstD0pi_DstnDof); ntp100->SetBranchAddress("DstD0pi_D0px",DstD0pi_D0px); ntp100->SetBranchAddress("DstD0pi_D0py",DstD0pi_D0py); ntp100->SetBranchAddress("DstD0pi_D0pz",DstD0pi_D0pz); ntp100->SetBranchAddress("DstD0pi_spipx",DstD0pi_spipx); ntp100->SetBranchAddress("DstD0pi_spipy",DstD0pi_spipy); ntp100->SetBranchAddress("DstD0pi_spipz",DstD0pi_spipz); ntp100->SetBranchAddress("DstD0pi_Kpx",DstD0pi_Kpx); ntp100->SetBranchAddress("DstD0pi_Kpy",DstD0pi_Kpy); ntp100->SetBranchAddress("DstD0pi_Kpz",DstD0pi_Kpz); ntp100->SetBranchAddress("DstD0pi_pipx",DstD0pi_pipx); ntp100->SetBranchAddress("DstD0pi_pipy",DstD0pi_pipy); ntp100->SetBranchAddress("DstD0pi_pipz",DstD0pi_pipz); ntp100->SetBranchAddress("DstD0pi_spiIdx",DstD0pi_spiIdx); TFile *outfile = new TFile("Test.root","RECREATE"); TTree *ntp1 = new TTree("ntp100","A new Tree"); ntp1->Branch("rndmNum",&rndmNum,"rndmNum/D"); ntp1->Branch("DstD0pi_num",&DstD0pi_num,"DstD0pi_num/I"); ntp1->Branch("DKKpi_num",&DKKpi_num,"DKKpi_num/I"); ntp1->Branch("DKpipi_num",&DKpipi_num,"DKpipi_num/I"); ntp1->Branch("DKpipi_unbias_num",&DKpipi_unbias_num,"DKpipi_unbias_num/I"); ntp1->Branch("runNumber",&runNumber,"runNumber/I"); ntp1->Branch("platform",&platform,"platform/I"); //ntp1->Branch("partition",&partition,"partition/I"); ntp1->Branch("upperID",&upperID,"upperID/I"); ntp1->Branch("lowerID",&lowerID,"lowerID/I"); ntp1->Branch("energyCM",&energyCM,"energyCM/D"); ntp1->Branch("ePlusMomentumX",&ePlusMomentumX,"ePlusMomentumX/D"); ntp1->Branch("ePlusMomentumY",&ePlusMomentumY,"ePlusMomentumY/D"); ntp1->Branch("ePlusMomentumZ",&ePlusMomentumZ,"ePlusMomentumZ/D"); ntp1->Branch("eMinusMomentumX",&eMinusMomentumX,"eMinusMomentumX/D"); ntp1->Branch("eMinusMomentumY",&eMinusMomentumY,"eMinusMomentumY/D"); ntp1->Branch("eMinusMomentumZ",&eMinusMomentumZ,"eMinusMomentumZ/D"); ntp1->Branch("boostCMtoLabX",&boostCMtoLabX,"boostCMtoLabX/D"); ntp1->Branch("boostCMtoLabY",&boostCMtoLabY,"boostCMtoLabY/D"); ntp1->Branch("boostCMtoLabZ",&boostCMtoLabZ,"boostCMtoLabZ/D"); ntp1->Branch("IPX",&IPX,"IPX/D"); ntp1->Branch("IPY",&IPY,"IPY/D"); ntp1->Branch("IPZ",&IPZ,"IPZ/D"); ntp1->Branch("xxpv",&xxpv,"xxpv/D"); ntp1->Branch("xypv",&xypv,"xypv/D"); ntp1->Branch("xzpv",&xzpv,"xzpv/D"); ntp1->Branch("yypv",&yypv,"yypv/D"); ntp1->Branch("yzpv",&yzpv,"yzpv/D"); ntp1->Branch("zzpv",&zzpv,"zzpv/D"); ntp1->Branch("ThrustMagCTrk",&ThrustMagCTrk,"ThrustMagCTrk/D"); ntp1->Branch("ThrustMagAll",&ThrustMagAll,"ThrustMagAll/D"); ntp1->Branch("ThrustCosThCTrk",&ThrustCosThCTrk,"ThrustCosThCTrk/D"); ntp1->Branch("ThrustCosThAll",&ThrustCosThAll,"ThrustCosThAll/D"); ntp1->Branch("ThrustPhiCTrk",&ThrustPhiCTrk,"ThrustPhiCTrk/D"); ntp1->Branch("ThrustPhiAll",&ThrustPhiAll,"ThrustPhiAll/D"); ntp1->Branch("R2CTrk",&R2CTrk,"R2CTrk/D"); ntp1->Branch("R2All",&R2All,"R2All/D"); ntp1->Branch("SphericityAll",&SphericityAll,"SphericityAll/D"); ntp1->Branch("xPrimaryVtx",&xPrimaryVtx,"xPrimaryVtx/D"); ntp1->Branch("yPrimaryVtx",&yPrimaryVtx,"yPrimaryVtx/D"); ntp1->Branch("zPrimaryVtx",&zPrimaryVtx,"zPrimaryVtx/D"); ntp1->Branch("probPrimaryVtx",&probPrimaryVtx,"probPrimaryVtx/D"); ntp1->Branch("mc_num",&mc_num,"mc_num/I"); ntp1->Branch("trk_num",&trk_num,"trk_num/I"); ntp1->Branch("trk_GenWord",trk_GenWord,"trk_GenWord/I"); ntp1->Branch("trk_PIDLongWord",trk_PIDLongWord,"trk_PIDLongWord/I"); ntp1->Branch("trk_eLHbit",trk_eLHbit,"trk_eLHbit/I"); ntp1->Branch("trk_charge",trk_charge,"trk_charge/D"); ntp1->Branch("trk_px",trk_px,"trk_px/D"); ntp1->Branch("trk_py",trk_py,"trk_py/D"); ntp1->Branch("trk_pz",trk_pz,"trk_pz/D"); ntp1->Branch("trk_p4ErrCovXX",trk_p4ErrCovXX,"trk_p4ErrCovXX/D"); ntp1->Branch("trk_p4ErrCovXY",trk_p4ErrCovXY,"trk_p4ErrCovXY/D"); ntp1->Branch("trk_p4ErrCovXZ",trk_p4ErrCovXZ,"trk_p4ErrCovXZ/D"); ntp1->Branch("trk_p4ErrCovYY",trk_p4ErrCovYY,"trk_p4ErrCovYY/D"); ntp1->Branch("trk_p4ErrCovYZ",trk_p4ErrCovYZ,"trk_p4ErrCovYZ/D"); ntp1->Branch("trk_p4ErrCovZZ",trk_p4ErrCovZZ,"trk_p4ErrCovZZ/D"); ntp1->Branch("trk_ndof",trk_ndof,"trk_ndof/I"); ntp1->Branch("trk_chisq",trk_chisq,"trk_chisq/D"); ntp1->Branch("trk_svtpat",trk_svtpat,"trk_svtpat/I"); ntp1->Branch("trk_phiatemc",trk_phiatemc,"trk_phiatemc/D"); ntp1->Branch("trk_dedxsvt",trk_dedxsvt,"trk_dedxsvt/D"); ntp1->Branch("trk_nsvt_dedx",trk_nsvt_dedx,"trk_nsvt_dedx/I"); ntp1->Branch("trk_nsvt_mom",trk_nsvt_mom,"trk_nsvt_mom/I"); ntp1->Branch("trk_dedxdch",trk_dedxdch,"trk_dedxdch/D"); ntp1->Branch("trk_ndch_dedx",trk_ndch_dedx,"trk_ndch_dedx/I"); ntp1->Branch("trk_ndch_mom",trk_ndch_mom,"trk_ndch_mom/I"); ntp1->Branch("trk_tht",trk_tht,"trk_tht/D"); ntp1->Branch("trk_thterr",trk_thterr,"trk_thterr/D"); ntp1->Branch("trk_nph",trk_nph,"trk_nph/I"); ntp1->Branch("trk_nphbkg",trk_nphbkg,"trk_nphbkg/I"); ntp1->Branch("trk_phicluster",trk_phicluster,"trk_phicluster/D"); ntp1->Branch("trk_ecal",trk_ecal,"trk_ecal/D"); ntp1->Branch("trk_lat",trk_lat,"trk_lat/D"); ntp1->Branch("trk_z42",trk_z42,"trk_z42/D"); ntp1->Branch("trk_z20",trk_z20,"trk_z20/D"); ntp1->Branch("trk_ncrys",trk_ncrys,"trk_ncrys/D"); ntp1->Branch("trk_s1s9",trk_s1s9,"trk_s1s9/D"); ntp1->Branch("trk_s9s25",trk_s9s25,"trk_s9s25/D"); ntp1->Branch("trk_ifrexp",trk_ifrexp,"trk_ifrexp/D"); ntp1->Branch("trk_ifrmeas",trk_ifrmeas,"trk_ifrmeas/D"); ntp1->Branch("trk_ifrbefore",trk_ifrbefore,"trk_ifrbefore/D"); ntp1->Branch("trk_ifrchisq",trk_ifrchisq,"trk_ifrchisq/D"); ntp1->Branch("DKKpi_pipx",DKKpi_pipx,"DKKpi_pipx/D"); //[DKKpi_num] ntp1->Branch("DKKpi_pipy",DKKpi_pipy,"DKKpi_pipy/D"); //[DKKpi_num] ntp1->Branch("DKKpi_pipz",DKKpi_pipz,"DKKpi_pipz/D"); //[DKKpi_num] ntp1->Branch("DKKpi_K1px",DKKpi_K1px,"DKKpi_K1px/D"); //[DKKpi_num] ntp1->Branch("DKKpi_K1py",DKKpi_K1py,"DKKpi_K1py/D"); //[DKKpi_num] ntp1->Branch("DKKpi_K1pz",DKKpi_K1pz,"DKKpi_K1pz/D"); //[DKKpi_num] ntp1->Branch("DKKpi_K2px",DKKpi_K2px,"DKKpi_K2px/D"); //[DKKpi_num] ntp1->Branch("DKKpi_K2py",DKKpi_K2py,"DKKpi_K2py/D"); //[DKKpi_num] ntp1->Branch("DKKpi_K2pz", DKKpi_K2pz,"DKKpi_K2pz/D"); //[DKKpi_num] ntp1->Branch("DKKpi_DVtxPosx",DKKpi_DVtxPosx,"DKKpi_DVtxPosx/D"); //[DKKpi_num] ntp1->Branch("DKKpi_DVtxPosy",DKKpi_DVtxPosy,"DKKpi_DVtxPosy/D"); //[DKKpi_num] ntp1->Branch("DKKpi_DVtxPosz",DKKpi_DVtxPosz,"DKKpi_DVtxPosz/D"); //[DKKpi_num] ntp1->Branch("DKKpi_DChiSquared",DKKpi_DChiSquared,"DKKpi_DChiSquared/D"); //[DKKpi_num] ntp1->Branch("DKKpi_DnDof", DKKpi_DnDof,"DKKpi_DnDof/I"); //[DKKpi_num] ntp1->Branch("DKKpi_piIdx", DKKpi_piIdx,"DKKpi_piIdx/I"); //[DKKpi_num] ntp1->Branch("DKKpi_K1Idx",DKKpi_K1Idx,"DKKpi_K1Idx/I"); //[DKKpi_num] ntp1->Branch("DKKpi_K2Idx",DKKpi_K2Idx,"DKKpi_K2Idx/I"); //[DKKpi_num] ntp1->Branch("DKKpi_DVtxSepChiSq", DKKpi_DVtxSepChiSq,"DKKpi_DVtxSepChiSq/D"); //[DKKpi_num] ntp1->Branch("DKpipi_Kpx",DKpipi_Kpx,"DKpipi_Kpx/D"); ntp1->Branch("DKpipi_Kpy",DKpipi_Kpy,"DKpipi_Kpy/D"); ntp1->Branch("DKpipi_Kpz",DKpipi_Kpz,"DKpipi_Kpz/D"); ntp1->Branch("DKpipi_pi1px",DKpipi_pi1px,"DKpipi_pi1px/D"); ntp1->Branch("DKpipi_pi1py",DKpipi_pi1py,"DKpipi_pi1py/D"); ntp1->Branch("DKpipi_pi1pz",DKpipi_pi1pz,"DKpipi_pi1pz/D"); ntp1->Branch("DKpipi_pi2px",DKpipi_pi2px,"DKpipi_pi2px/D"); ntp1->Branch("DKpipi_pi2py",DKpipi_pi2py,"DKpipi_pi2py/D"); ntp1->Branch("DKpipi_pi2pz",DKpipi_pi2pz,"DKpipi_pi2pz/D"); ntp1->Branch("DKpipi_DVtxPosx",DKpipi_DVtxPosx,"DKpipi_DVtxPosx/D"); ntp1->Branch("DKpipi_DVtxPosy",DKpipi_DVtxPosy,"DKpipi_DVtxPosy/D"); ntp1->Branch("DKpipi_DVtxPosz",DKpipi_DVtxPosz,"DKpipi_DVtxPosz/D"); ntp1->Branch("DKpipi_DChiSquared",DKpipi_DChiSquared,"DKpipi_DChiSquared/D"); ntp1->Branch("DKpipi_DnDof",DKpipi_DnDof,"DKpipi_DnDof/I"); ntp1->Branch("DKpipi_KIdx",DKpipi_KIdx,"DKpipi_KIdx/I"); ntp1->Branch("DKpipi_pi1Idx",DKpipi_pi1Idx,"DKpipi_pi1Idx/I"); ntp1->Branch("DKpipi_pi2Idx",DKpipi_pi2Idx,"DKpipi_pi2Idx/I"); ntp1->Branch("DKpipi_DVtxSepChiSq", DKpipi_DVtxSepChiSq,"DKpipi_DVtxSepChiSq/D"); ntp1->Branch("DKpipi_KTrkpx",DKpipi_KTrkpx,"DKpipi_KTrkpx/D"); ntp1->Branch("DKpipi_KTrkpy",DKpipi_KTrkpy,"DKpipi_KTrkpy/D"); ntp1->Branch("DKpipi_KTrkpz",DKpipi_KTrkpz,"DKpipi_KTrkpz/D"); ntp1->Branch("DKpipi_pi1Trkpx",DKpipi_pi1Trkpx,"DKpipi_pi1Trkpx/D"); ntp1->Branch("DKpipi_pi1Trkpy",DKpipi_pi1Trkpy,"DKpipi_pi1Trkpy/D"); ntp1->Branch("DKpipi_pi1Trkpz",DKpipi_pi1Trkpz,"DKpipi_pi1Trkpz/D"); ntp1->Branch("DKpipi_pi2Trkpx",DKpipi_pi2Trkpx,"DKpipi_pi2Trkpx/D"); ntp1->Branch("DKpipi_pi2Trkpy",DKpipi_pi2Trkpy,"DKpipi_pi2Trkpy/D"); ntp1->Branch("DKpipi_pi2Trkpz",DKpipi_pi2Trkpz,"DKpipi_pi2Trkpz/D"); ntp1->Branch("DKpipi_DTrkVtxPosx",DKpipi_DTrkVtxPosx,"DKpipi_DTrkVtxPosx/D"); ntp1->Branch("DKpipi_DTrkVtxPosy",DKpipi_DTrkVtxPosy,"DKpipi_DTrkVtxPosy/D"); ntp1->Branch("DKpipi_DTrkVtxPosz",DKpipi_DTrkVtxPosz,"DKpipi_DTrkVtxPosz/D"); ntp1->Branch("DKpipi_DTrkChiSquared",DKpipi_DTrkChiSquared,"DKpipi_DTrkChiSquared/D"); ntp1->Branch("DKpipi_DTrknDof",DKpipi_DTrknDof,"DKpipi_DTrknDof/I"); ntp1->Branch("DKpipi_KTrkIdx",DKpipi_KTrkIdx,"DKpipi_KTrkIdx/I"); ntp1->Branch("DKpipi_pi1TrkIdx",DKpipi_pi1TrkIdx,"DKpipi_pi1TrkIdx/I"); ntp1->Branch("DKpipi_pi2TrkIdx",DKpipi_pi2TrkIdx,"DKpipi_pi2TrkIdx/I"); ntp1->Branch("DKpipi_DTrkVtxSepChiSq", DKpipi_DTrkVtxSepChiSq,"DKpipi_DTrkVtxSepChiSq/D"); ntp1->Branch("D0Kpi_Kpx",D0Kpi_Kpx,"D0Kpi_Kpx/D"); ntp1->Branch("D0Kpi_Kpy",D0Kpi_Kpy,"D0Kpi_Kpy/D"); ntp1->Branch("D0Kpi_Kpz",D0Kpi_Kpz,"D0Kpi_Kpz/D"); ntp1->Branch("D0Kpi_pipx",D0Kpi_pipx,"D0Kpi_pipx/D"); ntp1->Branch("D0Kpi_pipy",D0Kpi_pipy,"D0Kpi_pipy/D"); ntp1->Branch("D0Kpi_pipz",D0Kpi_pipz,"D0Kpi_pipz/D"); ntp1->Branch("D0Kpi_D0VtxPosx",D0Kpi_D0VtxPosx,"D0Kpi_D0VtxPosx/D"); ntp1->Branch("D0Kpi_D0VtxPosy",D0Kpi_D0VtxPosy,"D0Kpi_D0VtxPosy/D"); ntp1->Branch("D0Kpi_D0VtxPosz",D0Kpi_D0VtxPosz,"D0Kpi_D0VtxPosz/D"); ntp1->Branch("D0Kpi_D0ChiSquared",D0Kpi_D0ChiSquared,"D0Kpi_D0ChiSquared/D"); ntp1->Branch("D0Kpi_D0nDof",D0Kpi_D0nDof,"D0Kpi_D0nDof/I"); ntp1->Branch("D0Kpi_KIdx",D0Kpi_KIdx,"D0Kpi_KIdx/I"); ntp1->Branch("D0Kpi_piIdx",D0Kpi_piIdx,"D0Kpi_piIdx/I"); ntp1->Branch("DstD0pi_DstVtxPosx",DstD0pi_DstVtxPosx,"DstD0pi_DstVtxPosx/D"); ntp1->Branch("DstD0pi_DstVtxPosy",DstD0pi_DstVtxPosy,"DstD0pi_DstVtxPosy/D"); ntp1->Branch("DstD0pi_DstVtxPosz",DstD0pi_DstVtxPosz,"DstD0pi_DstVtxPosz/D"); ntp1->Branch("DstD0pi_DstChiSquared",DstD0pi_DstChiSquared,"DstD0pi_DstChiSquared/D"); ntp1->Branch("DstD0pi_DstnDof",DstD0pi_DstnDof,"DstD0pi_DstnDof/I"); ntp1->Branch("DstD0pi_D0px",DstD0pi_D0px,"DstD0pi_D0px/D"); ntp1->Branch("DstD0pi_D0py",DstD0pi_D0py,"DstD0pi_D0py/D"); ntp1->Branch("DstD0pi_D0pz",DstD0pi_D0pz,"DstD0pi_D0pz/D"); ntp1->Branch("DstD0pi_spipx",DstD0pi_spipx,"DstD0pi_spipx/D"); ntp1->Branch("DstD0pi_spipy",DstD0pi_spipy,"DstD0pi_spipy/D"); ntp1->Branch("DstD0pi_spipz",DstD0pi_spipz,"DstD0pi_spipz/D"); ntp1->Branch("DstD0pi_Kpx",DstD0pi_Kpx,"DstD0pi_Kpx/D"); ntp1->Branch("DstD0pi_Kpy",DstD0pi_Kpy,"DstD0pi_Kpy/D"); ntp1->Branch("DstD0pi_Kpz",DstD0pi_Kpz, "DstD0pi_Kpz"); ntp1->Branch("DstD0pi_pipx",DstD0pi_pipx,"DstD0pi_pipx/D"); ntp1->Branch("DstD0pi_pipy",DstD0pi_pipy,"DstD0pi_pipy/D"); ntp1->Branch("DstD0pi_pipz",DstD0pi_pipz,"DstD0pi_pipz/D"); ntp1->Branch("DstD0pi_spiIdx",DstD0pi_spiIdx,"DstD0pi_spiIdx/I"); Int_t nentries = ntp100->GetEntries(); int NewEntries = 0; cout << "Total Entries: " << nentries << endl; for (Int_t i=0; iGetEntry(i); if((DstD0pi_num > 0) || (DKpipi_num > 0) || (DKpipi_unbias_num > 0) || (DKKpi_num > 0)){ NewEntries++; ntp1->Fill(); } }//for i ntp1->Write(); outfile->Write(); outfile->Close(); cout << "Total new entries: " << NewEntries << endl; timer.Stop(); cout << "Stopped timer" << endl; timer.Print(); }//main