NueVars.cxx

Go to the documentation of this file.

#include "3FlavorAna/Vars/NueVars.h"

#include "StandardRecord/Proxy/SRProxy.h"

namespace ana
{
  const Var kRawE    = SIMPLEVAR(energy.nue.rawgevslice);
  const Var kRawEshw = SIMPLEVAR(energy.nue.rawgevshw);
  const Var kRawEhad = SIMPLEVAR(energy.nue.rawgevhad);

  const Var kDistAllTop (
                         [](const caf::SRProxy* sr)
                         {
                           if(std::isnan(1.*sr->sel.nuecosrej.distallpngtop))
                             return -1000.0f;
                           return float(sr->sel.nuecosrej.distallpngtop);
                         });


  const Var kDistAllBottom (
                            [](const caf::SRProxy* sr)
                            {
                              if(std::isnan(1.*sr->sel.nuecosrej.distallpngbottom))
                                return -1000.0f;
                              return float(sr->sel.nuecosrej.distallpngbottom);
                            });


  const Var kDistAllWest (
                          [](const caf::SRProxy* sr)
                          {
                            if(std::isnan(1.*sr->sel.nuecosrej.distallpngwest))
                              return -1000.0f;
                            return float(sr->sel.nuecosrej.distallpngwest);
                          });


  const Var kDistAllEast (
                          [](const caf::SRProxy* sr)
                          {
                            if(std::isnan(1.*sr->sel.nuecosrej.distallpngeast))
                              return -1000.0f;
                            return float(sr->sel.nuecosrej.distallpngeast);
                          });


  const Var kDistAllBack (
                          [](const caf::SRProxy* sr)
                          {
                            if(std::isnan(1.*sr->sel.nuecosrej.distallpngback))
                              return -1000.0f;
                            return float(sr->sel.nuecosrej.distallpngback);
                          });


  const Var kDistAllFront (
                           [](const caf::SRProxy* sr)
                           {
                             if(std::isnan(1.*sr->sel.nuecosrej.distallpngfront))
                               return -1000.0f;
                             return float(sr->sel.nuecosrej.distallpngfront);
                           });


  const Var kDistNotTop (
                         [](const caf::SRProxy* sr)
                         {
                           float distnottop = std::min(kDistAllBottom(sr),kDistAllFront(sr));
                           distnottop = std::min(float(kDistAllBack(sr)),distnottop);
                           distnottop = std::min(float(kDistAllEast(sr)),distnottop);
                           distnottop = std::min(float(kDistAllWest(sr)),distnottop);

                           return distnottop;
                         });


  const Var kMaxY(
                  [](const caf::SRProxy* sr)
                  {
                    if(!sr->vtx.elastic.IsValid) return -1000.0f;
                    if(sr->vtx.elastic.fuzzyk.nshwlid == 0) return -1000.0f;
                    float maxyall = -999.0;
                    for(unsigned int i=0;i<sr->vtx.elastic.fuzzyk.nshwlid;i++){
                      maxyall = std::max(std::max(sr->vtx.elastic.fuzzyk.png[i].shwlid.start.Y(),
                                                  sr->vtx.elastic.fuzzyk.png[i].shwlid.stop.Y()),maxyall);
                    }
                    return maxyall;
                  });

  const Var kPtP = SIMPLEVAR(sel.nuecosrej.partptp);

  const Var kSparsenessAsymm = SIMPLEVAR(sel.nuecosrej.sparsenessasymm);


  const Var kMaxZ(
                  [](const caf::SRProxy* sr)
                  {
                    if(!sr->vtx.elastic.IsValid) return -1000.0f;
                    if(sr->vtx.elastic.fuzzyk.nshwlid == 0) return -1000.0f;
                    float maxzall = -999.0;
                    for(unsigned int i=0;i<sr->vtx.elastic.fuzzyk.nshwlid;i++){
                      maxzall = std::max(std::max(sr->vtx.elastic.fuzzyk.png[i].shwlid.start.Z(),
                                                  sr->vtx.elastic.fuzzyk.png[i].shwlid.stop.Z()),maxzall);
                    }
                    return maxzall;
                  });


  const Var kCosPID        = SIMPLEVAR(sel.nuecosrej.cospidlight);
  const Var kCosPIDContain = SIMPLEVAR(sel.nuecosrej.cospidcontain);

  const Var kCosPIDCoreBDT = SIMPLEVAR(sel.nuecosrej.cospidcorebdt);
  const Var kCosPIDPeriBDT = SIMPLEVAR(sel.nuecosrej.cospidperibdt);
}