getCosmicList.C File Reference

#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Cuts/TimingCuts.h"
#include "3FlavorAna/Cuts/NumuCuts2017.h"
#include "3FlavorAna/Vars/NumuVars.h"
#include "CAFAna/Core/EventList.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "CAFAna/Vars/GenieWeights.h"

Go to the source code of this file.

Functions
void	getCosmicList ()

Function Documentation

void getCosmicList

(

)

Definition at line 14 of file getCosmicList.C.

References caf::Proxy< caf::SRCosRej >::anglekal, caf::Proxy< caf::SRIDBranch >::cosrej, energy, ana::kCCE, ana::kHadEFrac, ana::kInTimingSideband, ana::kNumuContainFD, ana::kNumuCutFD2017, ana::kNumuNCRej, ana::kNumuQuality, ana::kStandardSpillCuts, ana::MakeTextListFile(), caf::Proxy< caf::SRSlice >::nhit, caf::Proxy< caf::StandardRecord >::sel, SIMPLEVAR, caf::Proxy< caf::StandardRecord >::slc, sr, and string.

{
  // files
  const std::string fd_data = "prod_sumrestricteddecaf_R17-03-01-prod3reco.k_fd_numi_fhc_full_v1_goodruns_numu_sa_2017";

  // cuts
  const Cut kNumuCosmicRejSA(
                           [](const caf::SRProxy* sr)
                           { return ( sr->sel.cosrej.anglekal  > 0.5 &&
                                      sr->sel.cosrej.numuSAcontpid > 0.535 &&
                                      //sr->sel.cosrej.numucontpid > 0.535 &&
                                      sr->slc.nhit < 400);
                           }
                           );

  const Cut kNumuFDSA = kNumuQuality && kNumuContainFD && kNumuNCRej && kNumuCosmicRejSA;

  // Vars

  const Var kIsSA(
                      [kNumuFDSA](const caf::SRProxy *sr)
                      {
                        return kNumuFDSA(sr);
                      });

  const Var kIsA17(
                      [kNumuCutFD2017](const caf::SRProxy *sr)
                      {
                        return kNumuCutFD2017(sr);
                      });

  const Var kCCESA = SIMPLEVAR(energy.numu.trkccE);

  
  MakeTextListFile(fd_data, { (kNumuFDSA || kNumuCutFD2017) && kInTimingSideband}, {"CosmicList.txt"} , {&kIsSA, &kIsA17, &kHadEFrac, &kCCE, &kCCESA}, &kStandardSpillCuts);

} // End of function