---

HardPhotonBarConfigurer.cxx

Go to the documentation of this file.

#include "FastShowerUtils/HardPhotonBarConfigurer.h"


#include "FastShowerUtils/Samplers/S0EmEcalBar1.h"
#include "FastShowerUtils/Samplers/PeakMeanS0EmEcalBar1.h"
#include "FastShowerUtils/Samplers/PeakSigmaS0EmEcalBar1.h"
#include "FastShowerUtils/Samplers/TailFracS0EmEcalBar1.h"
#include "FastShowerUtils/Samplers/TailAuxS0EmEcalBar1.h"
//Cell0
#include "FastShowerUtils/Samplers/C0EmEcalBar1.h"
#include "FastShowerUtils/Samplers/PeakMeanC0EmEcalBar1.h"
#include "FastShowerUtils/Samplers/PeakSigmaC0EmEcalBar1.h"
#include "FastShowerUtils/Samplers/TailFracC0EmEcalBar1.h"
#include "FastShowerUtils/Samplers/TailAuxC0EmEcalBar1.h"
//SN1
//#include "FastShowerUtils/SamplerRouterC0Sample.h"
#include "FastShowerUtils/SNSamplerRouterC0Sample.h"
//
#include "FastShowerUtils/Samplers/SNEmEcalBar1.h"
#include "FastShowerUtils/Samplers/PeakMeanSNEmEcalBar1.h"
#include "FastShowerUtils/Samplers/PeakSigmaSNEmEcalBar1.h"
#include "FastShowerUtils/Samplers/TailFracSNEmEcalBar1.h"
#include "FastShowerUtils/Samplers/TailAuxSNEmEcalBar1.h"
//SN2
#include "FastShowerUtils/Samplers/SNEmEcalBar2.h"
#include "FastShowerUtils/Samplers/TailFracSNEmEcalBar2.h"
#include "FastShowerUtils/Samplers/TailAuxSNEmEcalBar2.h"
//SN3 - as SN1
#include "FastShowerUtils/Samplers/SNEmEcalBar3.h"
//AN1
//#include "FastShowerUtils/SamplerRouterSNSample.h"
#include "FastShowerUtils/ANSamplerRouterSNSample.h"
#
#include "FastShowerUtils/Samplers/ANEmEcalBar1.h"
#include "FastShowerUtils/Samplers/PeakMeanANEmEcalBar1.h"
#include "FastShowerUtils/Samplers/PeakSigmaANEmEcalBar1.h"
//AN2
#include "FastShowerUtils/Samplers/ANEmEcalBar2.h"
#include "FastShowerUtils/Samplers/PeakMeanANEmEcalBar2.h"
#include "FastShowerUtils/Samplers/PeakSigmaANEmEcalBar2.h"
//
//Normalisers
#include "FastShowerUtils/Normalisers/EmEtaResponse.h"
#include "FastShowerUtils/Normalisers/EmHcalLeakageBar.h"
//
#include "FastShowerUtils/PtrConverter.h"
#include "FastShowerUtils/ISampler.h"
#include "FastShowerUtils/SamplerPtr.h"
#include "FastShowerUtils/INormaliser.h"
#include "FastShowerUtils/NormaliserPtr.h"
#include "FastShowerUtils/IDepositor.h"
#include "FastShowerUtils/DepositorPtr.h"
#include "FastShowerUtils/CoreDepositor.h"
#include "FastShowerUtils/PencilDepositor.h"
//
#include "FastShowerUtils/PtrAdapter.h"
//
#include "FastShowerUtils/Normalisers/EmEcalNormaliser.h"
#include "FastShowerUtils/Normalisers/EmHcalNormaliser.h"
//
#include <string>
#include <algorithm>

namespace FastShower{

  //
  HardPhotonBarConfigurer::HardPhotonBarConfigurer(Moni m):ConfigurerBase(m){
    this->initialise();
  }
  HardPhotonBarConfigurer::HardPhotonBarConfigurer(NoMoni nm):ConfigurerBase(nm){
    this->initialise();
  }
  void HardPhotonBarConfigurer::initialise(){
    //Slice0
    this->addFn("S0PeakMean",        new  PeakMeanS0EmEcalBar1);
    this->addFn("S0PeakSigma",       new  PeakSigmaS0EmEcalBar1);
    this->addFn("S0TailFrac",        new  TailFracS0EmEcalBar1);
    this->addFn("S0TailLowerBound",  new  TailAuxS0EmEcalBar1);
    //Cell0
    this->addFn("C0PeakMean",        new  PeakMeanC0EmEcalBar1);
    this->addFn("C0PeakSigma",       new  PeakSigmaC0EmEcalBar1);
    this->addFn("C0TailFrac",        new  TailFracC0EmEcalBar1);
    this->addFn("C0TailLowerBound",  new  TailAuxC0EmEcalBar1);
    //SN1
    this->addFn("SN1PeakMean",       new  PeakMeanSNEmEcalBar1);
    this->addFn("SN1PeakSigma",      new  PeakSigmaSNEmEcalBar1);
    this->addFn("SN1TailFrac",       new  TailFracSNEmEcalBar1);
    this->addFn("SN1TailLowerBound", new  TailAuxSNEmEcalBar1);
    //SN2 - uses SN1 functions for mean and sigma
    this->addFn("SN2PeakMean",       new  PeakMeanSNEmEcalBar1);
    this->addFn("SN2PeakSigma",      new  PeakSigmaSNEmEcalBar1);
    this->addFn("SN2TailFrac",       new  TailFracSNEmEcalBar2);
    this->addFn("SN2TailLowerBound", new  TailAuxSNEmEcalBar2);
    //SN3 - uses SN1 functions for mean and sigma
    this->addFn("SN3PeakMean",       new  PeakMeanSNEmEcalBar1);
    this->addFn("SN3PeakSigma",      new  PeakSigmaSNEmEcalBar1);
    //AN1 
    this->addFn("AN1PeakMean",       new  PeakMeanANEmEcalBar1);
    this->addFn("AN1PeakSigma",      new  PeakSigmaANEmEcalBar1);
    //AN2 
    this->addFn("AN2PeakMean",       new  PeakMeanANEmEcalBar2);
    this->addFn("AN2PeakSigma",      new  PeakSigmaANEmEcalBar2);
    //
    //*******************************************************
    //*              Functions for Normalisers              *
    //*******************************************************
    this->addFnPP("EtaResponse",    new  EmEtaResponse);
    //*******************************************************
    //*              Functions2 for Normalisers             *
    //*******************************************************
    this->addFnPP2("E/HcalFractions", new  EmHcalLeakageBar);
    //
    //*******************************************************
    //*       IUpdating Gaussians (Peaks for Samplers)      *
    //*******************************************************
    this->addIUG("AN1Peak", this->makeUpdatingGaussian0LP("AN1") );
    this->addIUG("AN2Peak", this->makeUpdatingGaussian0LP("AN2") );
    this->addIUG("C0Peak",  this->makeUpdatingGaussian0LP("C0")  );
    this->addIUG("SN1Peak", this->makeUpdatingGaussian0LP("SN1") );
    this->addIUG("SN2Peak", this->makeUpdatingGaussian0LP("SN2") );
    this->addIUG("SN3Peak", this->makeUpdatingGaussian0LP("SN3") );
    this->addIUG("S0Peak",  this->makeUpdatingGaussian0LP("S0")  );
    //
    //*******************************************************
    //*       Samplers                                      *
    //*******************************************************
    this->addSampler("S0",  new S0EmEcalBar1(this, "")  );
    this->addSampler("C0",  new C0EmEcalBar1(this, "")  );
    this->addSampler("SN1", new SNEmEcalBar1(this, "") );
    this->addSampler("SN2", new SNEmEcalBar2(this, "") );
    this->addSampler("SN3", new SNEmEcalBar3(this, "") );
    this->addSampler("AN1", new ANEmEcalBar1(this, "") );
    this->addSampler("AN2", new ANEmEcalBar2(this, "") );
    this->addSampler("SN",  this->makeSN() );
    this->addSampler("AN",  this->makeAN() );
    //
    //*******************************************************
    //*       Depositors                                    *
    //*******************************************************
    this->addDepositorsE();
    this->addDepositorsH();
    //*******************************************************
    //*       Normalisers                                   *
    //*******************************************************
    this->normalisers();
   }
  //
  ISampler* HardPhotonBarConfigurer::makeSN() const {
    
    std::vector< std::pair<double, ISampler*> > v;
    v.push_back(std::pair<double, ISampler*>(0.50, this->findSampler("SN1"))); 
    v.push_back(std::pair<double, ISampler*>(0.95, this->findSampler("SN2"))); 
    v.push_back(std::pair<double, ISampler*>(1.00, this->findSampler("SN3"))); 
    //return new SamplerRouterC0Sample(v, 0.0, "SNRouter");
    return new SNSamplerRouterC0Sample(v, 0.0, "SNRouter");
    
  }
  ISampler* HardPhotonBarConfigurer::makeAN() const {

    std::vector< std::pair<double, ISampler*> > v;
    v.push_back(std::pair<double, ISampler*>(0.05, this->findSampler("AN1"))); 
    v.push_back(std::pair<double, ISampler*>(1.00, this->findSampler("AN2"))); 
    //return new SamplerRouterSNSample(v, 0.0, "ANRouter");
    return new ANSamplerRouterSNSample(v, 0.0, "ANRouter");

  }
  void HardPhotonBarConfigurer::addDepositorsE(){
    
    this->addEdepositor( "EcalCore",
                         new CoreDepositor(
                                           this->findSampler("S0"),
                                           this->findSampler("C0"),
                                           this->findSampler("SN"),
                                           this->findSampler("AN")
                                           )
                         );
    
  }
  void HardPhotonBarConfigurer::addDepositorsH(){
    
    this->addHdepositor( "HCalPencil", new PencilDepositor);

  }
  void HardPhotonBarConfigurer::normalisers(){
    this->addNormaliser("EcalNormaliser",   
                        new EmEcalNormaliser(this,"HardPhotonBarEcalNormaliser"));
    this->addNormaliser("HcalNormaliser",   
                        new EmHcalNormaliser(this,"HardPhotonBarHcalNormaliser"));
  }
}//namespace

---