attach_mrcc_path_far.h

Go to the documentation of this file.

JobCPath& attach_mrcc_path_far(JobC &jc,
                           JobCPath& attach_to,
                           int spill_m,
                           TString config,
                           bool blinding,
                           SimFlag::SimFlag_t simflag,
                           TString output_filename = "ntupleMRNT.root")
{
  assert(simflag == SimFlag::kData || simflag == SimFlag::kMC);

  // For use in making unique job path names
  const TString blindStr = blinding ? "Blind" : "";

  JobCPath& removal = jc.Path.Create("muonremoval"+blindStr,
                                     "FiltTriggerPrescale::Ana ");

  JobCPath& muonreco = jc.Path.Create("muonreco"+blindStr);

  muonreco.PushBack("RemoveMuon", "Reco");
  muonreco.PushBack("StripCands", "Reco");
  muonreco.PushBack("ChopModule", "Reco");
  muonreco.PushBack("BiggestChopModule", "Reco");
  muonreco.PushBack("DeMuxDigitListModule", "Reco");
  muonreco.PushBack("StripSRListModule", "Reco");
  muonreco.PushBack("SliceSRListModule", "Reco");
  muonreco.PushBack("TrackCamListModule", "Reco");
  muonreco.PushBack("FitTrackCamListModule", "Reco");
  muonreco.PushBack("ClusterSRListModule", "Reco");
  muonreco.PushBack("SubShowerSRListModule", "Reco");
  muonreco.PushBack("ShowerSRListModule", "Reco");
  muonreco.PushBack("EventSRListModule", "Reco");
  muonreco.PushBack("RecordSetupModule", "Reco");

  if(simflag == SimFlag::kData)
    muonreco.PushBack("NtpBDLite", "Reco");

  muonreco.PushBack("NtpStModule", "Get");
  muonreco.PushBack("NtpSRModule", "Reco");

  if(simflag == SimFlag::kMC){
    muonreco.PushBack("NtpMCModule", "Reco"); //Necessary to write the MC info tree
    JobCModule& ntpmc = muonreco.Mod("NtpMCModule");
    ntpmc.Set("UseStandard=1");
    ntpmc.Set("RecordName=MuonRemoved");

    muonreco.PushBack("NtpTHModule", "Reco"); //Necessary to write the MC th tree
    JobCModule& ntpth = muonreco.Mod("NtpTHModule");
    ntpth.Set("UseStandard=1");
    ntpth.Set("RecordName=MuonRemoved");
  }

  muonreco.PushBack("FitTrackSAListModule", "Reco");
  muonreco.PushBack("NtpFitSAModule", "Reco");
  muonreco.PushBack("NtpMRModule", "Reco");
  muonreco.PushBack("Output", "Put");

  removal.Attach(&muonreco);

  removal.SetAllFilters(1);

  removal.Mod("FiltTriggerPrescale").Set(TString::Format("Spill_N=1    Spill_M=%d", spill_m));

  muonreco.SetAllFilters(0);

  //configure my code
  muonreco.Mod("RemoveMuon").Set("StripListOut=stripdigitlist");
  muonreco.Mod("StripCands").Set("keepdigitlist=stripdigitlist");

  //Post Muon removal reconstruction
  JobCModule& stripsrlist = muonreco.Mod("StripSRListModule");
  stripsrlist.Set("ListIn=canddigitlist");

  JobCModule& slicesrlist = muonreco.Mod("SliceSRListModule");
  slicesrlist.Set(("SliceListAlgConfig="+config).Data());
  slicesrlist.Set("ListIn=CandStripList");

  JobCModule& trackcamlist = muonreco.Mod("TrackCamListModule");
  trackcamlist.Set("TrackListAlgConfig=default");
  trackcamlist.Set("NameListIn=CandSliceList");
  trackcamlist.Set("NameListOut=CandTrackSRList");

  JobCModule& fittrackcam = muonreco.Mod("FitTrackCamListModule");
  fittrackcam.Set("NameListIn=CandTrackSRList");
  fittrackcam.Set("UseGeoSwimmer=1");

  muonreco.Mod("ClusterSRListModule").Set(("ClusterListAlgConfig="+config).Data());
  muonreco.Mod("SubShowerSRListModule").Set("SubShowerSRListAlgConfig=default");
  JobCModule& showersrlist = muonreco.Mod("ShowerSRListModule");
  showersrlist.Set("ShowerListAlgorithm=AlgShowerSSList");
  showersrlist.Set("ShowerListAlgConfig=default");
  showersrlist.Set("ListOut=CandShowerSRList");

  JobCModule& eventsrlist = muonreco.Mod("EventSRListModule");
  eventsrlist.Set("EventListAlgorithm=AlgEventSSList");
  eventsrlist.Set(("EventListAlgConfig="+config).Data());







  muonreco.Mod("ChopModule").Set("ChopAlgorithm=AlgChopListFar");
  muonreco.Mod("BiggestChopModule").Set("MinEnergy=0.0"); //No cut.
  muonreco.Mod("BiggestChopModule").Set("OmitLiEvents=1"); //kill LI
  muonreco.Mod("BiggestChopModule").Set("LiVetoWindowLow=-1.0e-6");
  muonreco.Mod("BiggestChopModule").Set("LiVetoWindowHigh=31.0e-6");


  muonreco.Mod("DeMuxDigitListModule").Set("SwitchPersToTemp=1");
  muonreco.Mod("DeMuxDigitListModule").Set("NameListOut=altdemux");
  muonreco.Mod("DeMuxDigitListModule").Set("DeMuxDigitListAlgConfig=devel");
  muonreco.Mod("StripSRListModule").Set("ListIn=altdemux");





  //configure output
  JobCModule& output = muonreco.Mod("Output");
  output.Set("FileName="+output_filename);

  muonreco.Mod("NtpStModule").Set("RecordName=MuonRemoved");
  muonreco.Mod("NtpSRModule").Set("UseStandard=1");
  muonreco.Mod("NtpSRModule").Set("RecordName=MuonRemoved");

  output.Cmd("DefineStream NtpSt NtpStRecord MuonRemoved");
  output.Cmd("DefineStream NtpMR NtpMRRecord");
  output.Set("Streams=NtpSt,NtpMR");


  if(blinding){
    JobCPath& mrblindpath = jc.Path.Create("mrblindpath",
                                           "Blinder::Ana "
                                           );

    //Explicitly turn on the Blinder Filtering
    mrblindpath.Node("Blinder::Ana").FilterOn();

    JobCPath& mrblindSR = jc.Path.Create("mrblindSR",
                                         "Output::Put ");

    mrblindSR.SetAllFilters(false);

    mrblindpath.Attach(&mrblindSR);
    muonreco.Attach(&mrblindpath);

    JobCModule& op = mrblindSR.Mod("Output");
    op.Set("FileName=ntupleMRNT.root");
    op.Cmd("DefineStream NtpSt NtpStRecord MuonRemoved");
    op.Cmd("DefineStream NtpMR NtpMRRecord");
    op.Set("Streams=NtpSt,NtpMR");
  } // end if blinding


  attach_to.Attach(&removal);

  return muonreco;
}

---