The (almost) Complete Code

So to find the calibrated data, and store it in a BrTileRdo object for each input event, one could have code like

    BrModuleContainer* mainCont =
      new BrModuleContainer("mainCont", "Main Cont");
    BrRdoModuleZDC* zdcModule =
      new BrRdoModuleZDC("ZDC", "ZDC reducer");
    mainCont->AddModule(zdcModule);
    BrRdoModuleBB* bbModule = 
      new BrRdoModuleBB("BB", "Beam-Beam Reducer");
    mainCont->AddModule(bbModule);
    // ... more vertex modules (?)
    BrTileRdoModule* tileModule =
      new BrTileRdoModule("MultTile", "Tile Reducer");
    mainCont->AddModule(tileModule);
    mainCont->Init();
    // From here on normal event processing

To find the functions $S_{cut}(V_z)$, one would take the output from above and pass it through code like the below:

    BrModuleContainer* mainCont = 
      new BrModuleContainer("mainCont", "Main Cont");
    BrTriggerFilter* trigFilt = 
      new BrTriggerFilter("trigFilt", "Trigger filter");
    mainCont->AddModule(trigFilt);
    trigFilt->AddTrigger(4);
    BrRdoModuleZDC* zdcModule = 
      new BrRdoModuleZDC("ZDC", "ZDC reducer");
    mainCont->AddModule(zdcModule);
    BrVertexFilter* vtxFilter = 
      new BrVertexFilter("VTX", "VerTeX Filter");
    mainCont->AddModule(vtxFilter);
    vtxFilter->SetMaxZdcVertexZ(50);
    BrTileCentCalModule* centcalModule = 
      new BrTileCentCalModule("centTile", "Centrality Calib");
    mainCont->AddModule(centcalModule);
    mainCont->Init();
    // From here on normal event processing

To find the centrality limits of a given event, one should have:

    BrModuleContainer* mainCont = 
      new BrModuleContainer("mainCont", "Main Cont");
    BrRdoModuleZDC* zdcModule = 
      new BrRdoModuleZDC("ZDC", "ZDC reducer");
    mainCont->AddModule(zdcModule);
    BrRdoModuleBB* bbModule = 
      new BrRdoModuleBB("BB", "Beam-Beam Reducer");
    mainCont->AddModule(bbModule);
    // ... more vertex modules (?)
    BrTileCentModule* centModule = 
      new BrTileCalModule("centTile", "Centrality finder");
    mainCont->AddModule(centModule);
    BrCentFilter* centfilterModule =
      new BrCentFilter("centFilter", "Centrality filter");
    mainCont->AddModule(centfilterModule);
    centfilterModule->SetLowCut(0); // From 0 to 10%
    centfilterModule->SetHighCut(.1); // From 0 to 10%
    mainCont->Init();
    // From here on normal event processing