#include "TROOT.h" #include "TApplication.h" #include "TClass.h" #include "TClassTable.h" #include "TCollection.h" #include "TGraph.h" #include "TPostScript.h" #include "TObject.h" #include "TFile.h" #include "TH1.h" #include "TH2.h" #include "TCanvas.h" #include "TF1.h" #include "TRandom.h" #include "TStopwatch.h" #include "TStyle.h" #include "TChain.h" #include "TPaveStats.h" #include "TLegend.h" #include "TPaveLabel.h" #include "TArrayD.h" // #include #include #include #include #include #include #include char ChainFileName[250]; ifstream inpfile("SkimTuple.inp"); Int_t run_range_min, run_range_max, nevmax; ifstream * fin; //Declaration of leaf types Double_t rndmNum; Int_t DstD0pi_num; Int_t DKpipi_num; Int_t DKpipi_unbias_num; Int_t DKKpi_num; Int_t runNumber; Int_t platform; //Int_t partition; Int_t upperID; Int_t lowerID; Double_t energyCM; Double_t ePlusMomentumX; Double_t ePlusMomentumY; Double_t ePlusMomentumZ; Double_t eMinusMomentumX; Double_t eMinusMomentumY; Double_t eMinusMomentumZ; Double_t boostCMtoLabX; Double_t boostCMtoLabY; Double_t boostCMtoLabZ; Double_t IPX; Double_t IPY; Double_t IPZ; Double_t xxpv; Double_t xypv; Double_t xzpv; Double_t yypv; Double_t yzpv; Double_t zzpv; Double_t ThrustMagCTrk; Double_t ThrustMagAll; Double_t ThrustCosThCTrk; Double_t ThrustCosThAll; Double_t ThrustPhiCTrk; Double_t ThrustPhiAll; Double_t R2CTrk; Double_t R2All; Double_t SphericityAll; Double_t xPrimaryVtx; Double_t yPrimaryVtx; Double_t zPrimaryVtx; Double_t probPrimaryVtx; Int_t mc_num; Int_t trk_num; Int_t trk_GenWord[29]; //[trk_num] Int_t trk_PIDLongWord[29]; //[trk_num] Int_t trk_eLHbit[29]; //[trk_num] Double_t trk_charge[29]; //[trk_num] Double_t trk_px[29]; //[trk_num] Double_t trk_py[29]; //[trk_num] Double_t trk_pz[29]; //[trk_num] Double_t trk_p4ErrCovXX[29]; //[trk_num] Double_t trk_p4ErrCovXY[29]; //[trk_num] Double_t trk_p4ErrCovXZ[29]; //[trk_num] Double_t trk_p4ErrCovYY[29]; //[trk_num] Double_t trk_p4ErrCovYZ[29]; //[trk_num] Double_t trk_p4ErrCovZZ[29]; //[trk_num] Int_t trk_ndof[29]; //[trk_num] Double_t trk_chisq[29]; //[trk_num] Int_t trk_svtpat[29]; //[trk_num] Double_t trk_phiatemc[29]; //[trk_num] Double_t trk_dedxsvt[29]; //[trk_num] Int_t trk_nsvt_dedx[29]; //[trk_num] Int_t trk_nsvt_mom[29]; //[trk_num] Double_t trk_dedxdch[29]; //[trk_num] Int_t trk_ndch_dedx[29]; //[trk_num] Int_t trk_ndch_mom[29]; //[trk_num] Double_t trk_tht[29]; //[trk_num] Double_t trk_thterr[29]; //[trk_num] Int_t trk_nph[29]; //[trk_num] Int_t trk_nphbkg[29]; //[trk_num] Double_t trk_phicluster[29]; //[trk_num] Double_t trk_ecal[29]; //[trk_num] Double_t trk_lat[29]; //[trk_num] Double_t trk_z42[29]; //[trk_num] Double_t trk_z20[29]; //[trk_num] Double_t trk_ncrys[29]; //[trk_num] Double_t trk_s1s9[29]; //[trk_num] Double_t trk_s9s25[29]; //[trk_num] Double_t trk_ifrexp[29]; //[trk_num] Double_t trk_ifrmeas[29]; //[trk_num] Double_t trk_ifrbefore[29]; //[trk_num] Double_t trk_ifrchisq[29]; //[trk_num] Double_t DKKpi_pipx[6]; //[DKKpi_num] Double_t DKKpi_pipy[6]; //[DKKpi_num] Double_t DKKpi_pipz[6]; //[DKKpi_num] Double_t DKKpi_K1px[6]; //[DKKpi_num] Double_t DKKpi_K1py[6]; //[DKKpi_num] Double_t DKKpi_K1pz[6]; //[DKKpi_num] Double_t DKKpi_K2px[6]; //[DKKpi_num] Double_t DKKpi_K2py[6]; //[DKKpi_num] Double_t DKKpi_K2pz[6]; //[DKKpi_num] Double_t DKKpi_DVtxPosx[6]; //[DKKpi_num] Double_t DKKpi_DVtxPosy[6]; //[DKKpi_num] Double_t DKKpi_DVtxPosz[6]; //[DKKpi_num] Double_t DKKpi_DChiSquared[6]; //[DKKpi_num] Int_t DKKpi_DnDof[6]; //[DKKpi_num] Int_t DKKpi_piIdx[6]; //[DKKpi_num] Int_t DKKpi_K1Idx[6]; //[DKKpi_num] Int_t DKKpi_K2Idx[6]; //[DKKpi_num] Double_t DKKpi_DVtxSepChiSq[6]; //[DKKpi_num] Double_t DKpipi_Kpx[6]; //[DKpipi_num] Double_t DKpipi_Kpy[6]; //[DKpipi_num] Double_t DKpipi_Kpz[6]; //[DKpipi_num] Double_t DKpipi_pi1px[6]; //[DKpipi_num] Double_t DKpipi_pi1py[6]; //[DKpipi_num] Double_t DKpipi_pi1pz[6]; //[DKpipi_num] Double_t DKpipi_pi2px[6]; //[DKpipi_num] Double_t DKpipi_pi2py[6]; //[DKpipi_num] Double_t DKpipi_pi2pz[6]; //[DKpipi_num] Double_t DKpipi_DVtxPosx[6]; //[DKpipi_num] Double_t DKpipi_DVtxPosy[6]; //[DKpipi_num] Double_t DKpipi_DVtxPosz[6]; //[DKpipi_num] Double_t DKpipi_DChiSquared[6]; //[DKpipi_num] Int_t DKpipi_DnDof[6]; //[DKpipi_num] Int_t DKpipi_KIdx[6]; //[DKpipi_num] Int_t DKpipi_pi1Idx[6]; //[DKpipi_num] Int_t DKpipi_pi2Idx[6]; //[DKpipi_num] Double_t DKpipi_DVtxSepChiSq[4]; //[DKpipi_num] Double_t DKpipi_KTrkpx[5]; //[DKpipi_unbias_num] Double_t DKpipi_KTrkpy[5]; //[DKpipi_unbias_num] Double_t DKpipi_KTrkpz[5]; //[DKpipi_unbias_num] Double_t DKpipi_pi1Trkpx[5]; //[DKpipi_unbias_num] Double_t DKpipi_pi1Trkpy[5]; //[DKpipi_unbias_num] Double_t DKpipi_pi1Trkpz[5]; //[DKpipi_unbias_num] Double_t DKpipi_pi2Trkpx[5]; //[DKpipi_unbias_num] Double_t DKpipi_pi2Trkpy[5]; //[DKpipi_unbias_num] Double_t DKpipi_pi2Trkpz[5]; //[DKpipi_unbias_num] Double_t DKpipi_DTrkVtxPosx[5]; //[DKpipi_unbias_num] Double_t DKpipi_DTrkVtxPosy[5]; //[DKpipi_unbias_num] Double_t DKpipi_DTrkVtxPosz[5]; //[DKpipi_unbias_num] Double_t DKpipi_DTrkChiSquared[5]; //[DKpipi_unbias_num] Int_t DKpipi_DTrknDof[5]; //[DKpipi_unbias_num] Int_t DKpipi_KTrkIdx[5]; //[DKpipi_unbias_num] Int_t DKpipi_pi1TrkIdx[5]; //[DKpipi_unbias_num] Int_t DKpipi_pi2TrkIdx[5]; //[DKpipi_unbias_num] Double_t DKpipi_DTrkVtxSepChiSq[7]; //[DKpipi_unbias_num] Double_t D0Kpi_Kpx[2]; //[DstD0pi_num] Double_t D0Kpi_Kpy[2]; //[DstD0pi_num] Double_t D0Kpi_Kpz[2]; //[DstD0pi_num] Double_t D0Kpi_pipx[2]; //[DstD0pi_num] Double_t D0Kpi_pipy[2]; //[DstD0pi_num] Double_t D0Kpi_pipz[2]; //[DstD0pi_num] Double_t D0Kpi_D0VtxPosx[2]; //[DstD0pi_num] Double_t D0Kpi_D0VtxPosy[2]; //[DstD0pi_num] Double_t D0Kpi_D0VtxPosz[2]; //[DstD0pi_num] Double_t D0Kpi_D0ChiSquared[2]; //[DstD0pi_num] Int_t D0Kpi_D0nDof[2]; //[DstD0pi_num] Int_t D0Kpi_KIdx[2]; //[DstD0pi_num] Int_t D0Kpi_piIdx[2]; //[DstD0pi_num] Double_t DstD0pi_DstVtxPosx[2]; //[DstD0pi_num] Double_t DstD0pi_DstVtxPosy[2]; //[DstD0pi_num] Double_t DstD0pi_DstVtxPosz[2]; //[DstD0pi_num] Double_t DstD0pi_DstChiSquared[2]; //[DstD0pi_num] Int_t DstD0pi_DstnDof[2]; //[DstD0pi_num] Double_t DstD0pi_D0px[2]; //[DstD0pi_num] Double_t DstD0pi_D0py[2]; //[DstD0pi_num] Double_t DstD0pi_D0pz[2]; //[DstD0pi_num] Double_t DstD0pi_spipx[2]; //[DstD0pi_num] Double_t DstD0pi_spipy[2]; //[DstD0pi_num] Double_t DstD0pi_spipz[2]; //[DstD0pi_num] Double_t DstD0pi_Kpx[2]; //[DstD0pi_num] Double_t DstD0pi_Kpy[2]; //[DstD0pi_num] Double_t DstD0pi_Kpz[2]; //[DstD0pi_num] Double_t DstD0pi_pipx[2]; //[DstD0pi_num] Double_t DstD0pi_pipy[2]; //[DstD0pi_num] Double_t DstD0pi_pipz[2]; //[DstD0pi_num] Int_t DstD0pi_spiIdx[2]; //[DstD0pi_num]