DriftAna Class Reference

Inheritance diagram for DriftAna:

Public Member Functions
	DriftAna (fhicl::ParameterSet const &p)
	DriftAna (DriftAna const &)=delete
	DriftAna (DriftAna &&)=delete
DriftAna &	operator= (DriftAna const &)=delete
DriftAna &	operator= (DriftAna &&)=delete
void	analyze (art::Event const &e) override
void	beginJob () override
Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDAnalyzer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept
	Analyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	Analyzer (Table< Config > const &config)
void	doBeginJob (SharedResources const &resources)
void	doEndJob ()
void	doRespondToOpenInputFile (FileBlock const &fb)
void	doRespondToCloseInputFile (FileBlock const &fb)
void	doRespondToOpenOutputFiles (FileBlock const &fb)
void	doRespondToCloseOutputFiles (FileBlock const &fb)
bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
Public Member Functions inherited from art::Observer
	~Observer () noexcept
	Observer (Observer const &)=delete
	Observer (Observer &&)=delete
Observer &	operator= (Observer const &)=delete
Observer &	operator= (Observer &&)=delete
void	registerProducts (ProductDescriptions &, ModuleDescription const &)
void	fillDescriptions (ModuleDescription const &)
fhicl::ParameterSetID	selectorConfig () const
Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept
	ModuleBase ()
ModuleDescription const &	moduleDescription () const
void	setModuleDescription (ModuleDescription const &)
std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const
void	sortConsumables (std::string const &current_process_name)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Attributes
TH1D *	drifttimenosce
TH1D *	drifttimeave
TH1D *	drifttimeint
TH1D *	fEfield
TH1D *	fdeltatnosce
TH1D *	fdeltatave
TH1D *	fdeltatint
TH1D *	fdriftvnosce
TH1D *	fdriftvave
TH1D *	fdriftvint
TH1D *	fefieldnosce
TH1D *	fefieldave
TH1D *	fefieldint

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer
Protected Member Functions inherited from art::Observer
std::string const &	processName () const
bool	wantAllEvents () const noexcept
bool	wantEvent (ScheduleID id, Event const &e) const
Handle< TriggerResults >	getTriggerResults (Event const &e) const
	Observer (fhicl::ParameterSet const &config)
	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)
Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 30 of file DriftAna_module.cc.

Constructor & Destructor Documentation

DriftAna::DriftAna ( fhicl::ParameterSet const &  p )

explicit

Definition at line 69 of file DriftAna_module.cc.

  : EDAnalyzer{p}  // ,
  // More initializers here.
{
  // Call appropriate consumes<>() for any products to be retrieved by this module.
}

DriftAna::DriftAna ( DriftAna const &  )

delete

DriftAna::DriftAna ( DriftAna &&  )

delete

Member Function Documentation

void DriftAna::analyze ( art::Event const &  e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 76 of file DriftAna_module.cc.

{
  auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService>()->DataFor(e);
  double driftv = detProp.DriftVelocity(detProp.Efield(),
                                        detProp.Temperature());
  double efield = detProp.Efield();
  std::cout<<"E field = "<<efield<<" Nominal drift velocity = "<<driftv<<std::endl;

  const int n = 100;
  const double e0 = 0.4;
  const double e1 = 0.7;
  std::vector<double> vecv;
  std::vector<double> vece;
  for (int i = 0; i<n; ++i){
    double e_field = e0 + i*(e1-e0)/n;
    double drift_v = detProp.DriftVelocity(e_field,
                                           detProp.Temperature());
    vecv.push_back(drift_v);
    vece.push_back(e_field);
  }
  TGraph *gr = new TGraph(n, &vecv[0], &vece[0]);
  for (int i = 1; i<=drifttimeave->GetNbinsX(); ++i){
    double xyz[3] = {0};
    xyz[0] = drifttimeave->GetBinCenter(i);
    xyz[1] = 300;
    xyz[2] = 350;
    std::cout<<xyz[0]<<" "<<xyz[1]<<" "<<xyz[2]<<std::endl;
    double driftdis = 358.6 - std::abs(xyz[0]);
    double drifttime = driftdis/driftv;
    drifttimenosce->SetBinContent(i, drifttime);
    geo::Vector_t posOffsets;
    auto const* SCE = lar::providerFrom<spacecharge::SpaceChargeService>();
    if (SCE->EnableSimSpatialSCE() == true){
      posOffsets = SCE->GetPosOffsets({ xyz[0], xyz[1], xyz[2] });
    }
    posOffsets.SetX(-posOffsets.X());
    driftdis += posOffsets.X();
    drifttime = driftdis/driftv;
    drifttimeave->SetBinContent(i, drifttime);
    
    double dir = 1;
    if (xyz[0]<0) dir = -1;
    double xyz1[3];
    xyz1[0] = xyz[0];
    xyz1[1] = 300;
    xyz1[2] = 350;
    double deltax = 1;
    double dt = 0;
    while (std::abs(xyz1[0])<358.6){
      geo::Vector_t fEfieldOffsets = SCE->GetEfieldOffsets({ xyz1[0], xyz1[1], xyz1[2] });
      double EField = std::sqrt( (efield + efield*fEfieldOffsets.X())*(efield + efield*fEfieldOffsets.X()) +
                                 (efield*fEfieldOffsets.Y()*efield*fEfieldOffsets.Y()) +
                                 (efield*fEfieldOffsets.Z()*efield*fEfieldOffsets.Z()) );
      double v = detProp.DriftVelocity(EField,
                                       detProp.Temperature());
      dt += deltax/v;
      xyz1[0] += deltax*dir;
    }
    drifttimeint->SetBinContent(i, dt);

    geo::Vector_t fEfieldOffsets = SCE->GetEfieldOffsets({ xyz[0], xyz[1], xyz[2] });
    double EField = std::sqrt( (efield + efield*fEfieldOffsets.X())*(efield + efield*fEfieldOffsets.X()) +
                               (efield*fEfieldOffsets.Y()*efield*fEfieldOffsets.Y()) +
                               (efield*fEfieldOffsets.Z()*efield*fEfieldOffsets.Z()) );
    fEfield->SetBinContent(i, EField);
  }

  for (int i = 1; i<=drifttimeave->GetNbinsX() - 1; ++i){
    double x = drifttimeave->GetBinCenter(i);
    if (x<0){
      fdeltatnosce->SetBinContent(i, drifttimenosce->GetBinContent(i) - drifttimenosce->GetBinContent(i-1));
      fdeltatave->SetBinContent(i, drifttimeave->GetBinContent(i) - drifttimeave->GetBinContent(i-1));
      fdeltatint->SetBinContent(i, drifttimeint->GetBinContent(i) - drifttimeint->GetBinContent(i-1));
    }
    else{
      fdeltatnosce->SetBinContent(i, drifttimenosce->GetBinContent(i) - drifttimenosce->GetBinContent(i+1));
      fdeltatave->SetBinContent(i, drifttimeave->GetBinContent(i) - drifttimeave->GetBinContent(i+1));
      fdeltatint->SetBinContent(i, drifttimeint->GetBinContent(i) - drifttimeint->GetBinContent(i+1));
    }
    fdriftvnosce->SetBinContent(i, fdeltatnosce->GetBinWidth(i)/fdeltatnosce->GetBinContent(i));
    fdriftvave->SetBinContent(i, fdeltatave->GetBinWidth(i)/fdeltatave->GetBinContent(i));
    fdriftvint->SetBinContent(i, fdeltatint->GetBinWidth(i)/fdeltatint->GetBinContent(i));
    fefieldnosce->SetBinContent(i, gr->Eval(fdeltatnosce->GetBinWidth(i)/fdeltatnosce->GetBinContent(i)));
    fefieldave->SetBinContent(i, gr->Eval(fdeltatave->GetBinWidth(i)/fdeltatave->GetBinContent(i)));
    fefieldint->SetBinContent(i, gr->Eval(fdeltatint->GetBinWidth(i)/fdeltatint->GetBinContent(i)));
  }
}

void DriftAna::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 164 of file DriftAna_module.cc.

{
    art::ServiceHandle<art::TFileService> tfs;
    drifttimenosce = tfs->make<TH1D>("drifttimenosce",";x (cm);Drift time (#mus)",120,-360,360);
    drifttimeave = tfs->make<TH1D>("drifttimeave",";x (cm);Drift time (#mus)",120,-360,360);
    drifttimeint = tfs->make<TH1D>("drifttimeint",";x (cm);Drift time (#mus)",120,-360,360);
    fEfield = tfs->make<TH1D>("fEfield","E field; x (cm);E field (kV/cm)",120,-360,360);

    fdeltatnosce = tfs->make<TH1D>("fdeltatnosce",";x (cm);#Delta t (#mus)", 120, -360, 360);
    fdeltatave = tfs->make<TH1D>("fdeltatave",";x (cm);#Delta t (#mus)", 120, -360, 360);
    fdeltatint = tfs->make<TH1D>("fdeltatint",";x (cm);#Delta t (#mus)", 120, -360, 360);

    fdriftvnosce = tfs->make<TH1D>("fdriftvnosce",";x (cm);Drift velocity (cm/#mus)", 120, -360, 360);
    fdriftvave = tfs->make<TH1D>("fdriftvave",";x (cm);Drift velocity (cm/#mus)", 120, -360, 360);
    fdriftvint = tfs->make<TH1D>("fdriftvint",";x (cm);Drift velocity (cm/#mus)", 120, -360, 360);

    fefieldnosce = tfs->make<TH1D>("fefieldnosce",";x (cm);E field (kV/cm)", 120, -360, 360);
    fefieldave = tfs->make<TH1D>("fefieldave",";x (cm);E field (kV/cm)", 120, -360, 360);
    fefieldint = tfs->make<TH1D>("fefieldint",";x (cm);E field (kV/cm)", 120, -360, 360);

}

DriftAna& DriftAna::operator= ( DriftAna const &  )

delete

DriftAna& DriftAna::operator= ( DriftAna &&  )

delete

Member Data Documentation

TH1D* DriftAna::drifttimeave

private

Definition at line 51 of file DriftAna_module.cc.

TH1D* DriftAna::drifttimeint

private

Definition at line 52 of file DriftAna_module.cc.

TH1D* DriftAna::drifttimenosce

private

Definition at line 50 of file DriftAna_module.cc.

TH1D* DriftAna::fdeltatave

private

Definition at line 56 of file DriftAna_module.cc.

TH1D* DriftAna::fdeltatint

private

Definition at line 57 of file DriftAna_module.cc.

TH1D* DriftAna::fdeltatnosce

private

Definition at line 55 of file DriftAna_module.cc.

TH1D* DriftAna::fdriftvave

private

Definition at line 60 of file DriftAna_module.cc.

TH1D* DriftAna::fdriftvint

private

Definition at line 61 of file DriftAna_module.cc.

TH1D* DriftAna::fdriftvnosce

private

Definition at line 59 of file DriftAna_module.cc.

TH1D* DriftAna::fEfield

private

Definition at line 53 of file DriftAna_module.cc.

TH1D* DriftAna::fefieldave

private

Definition at line 64 of file DriftAna_module.cc.

TH1D* DriftAna::fefieldint

private

Definition at line 65 of file DriftAna_module.cc.

TH1D* DriftAna::fefieldnosce

private

Definition at line 63 of file DriftAna_module.cc.

---

The documentation for this class was generated from the following file:

• protoduneana/protoduneana/singlephase/driftvel/DriftAna_module.cc