Calibration_module.cc File Reference

#include <iostream>
#include "art/Framework/Core/EDAnalyzer.h"
#include "canvas/Persistency/Common/FindOneP.h"
#include "art/Framework/Core/ModuleMacros.h"
#include "art/Framework/Principal/Event.h"
#include "art/Framework/Principal/Handle.h"
#include "art/Framework/Principal/Run.h"
#include "art/Framework/Principal/SubRun.h"
#include "canvas/Utilities/InputTag.h"
#include "art/Framework/Services/Optional/TFileService.h"
#include <TH2I.h>
#include <TProfile.h>
#include <TMath.h>
#include "fhiclcpp/ParameterSet.h"
#include "messagefacility/MessageLogger/MessageLogger.h"
#include "DDTBaseDataProducts/BaseProducts.h"
#include "DDTBaseDataProducts/HitList.h"
#include "DDTBaseDataProducts/DAQHit.h"
#include "DDTBaseDataProducts/Track3D.h"
#include "DDTBaseDataProducts/HitSet.h"
#include "OnlineCalibration/ChannelMapService.h"

Go to the source code of this file.

Classes
class	novaddt::Calibration

Namespaces
	novaddt

Functions
double	FindW (const novaddt::Track3D &t, const novaddt::DAQHit &h)
double	FindE (const novaddt::Track3D &t, const novaddt::DAQHit &h)

Function Documentation

double FindE	(	const novaddt::Track3D &	t,
		const novaddt::DAQHit &	h
	)

Definition at line 90 of file Calibration_module.cc.

References novaddt::DAQHit::ADC(), dir, novaddt::Track3D::End(), ana::Sqrt(), novaddt::Track3D::Start(), Unit(), novaddt::ADC::val, novaddt::View::val, and novaddt::DAQHit::View().

Referenced by novaddt::Calibration::analyze().

                                                            {
    TVector3 dir=(t.End()-t.Start()).Unit();
    //geometric factor = sin(angle between dir and cell orientation)
    double Factor=1;
    if(h.View().val==1)Factor=TMath::Sqrt(1-dir.X()*dir.X());
    if(h.View().val==2)Factor=TMath::Sqrt(1-dir.Y()*dir.Y());
    return h.ADC().val*Factor;
}

double FindW	(	const novaddt::Track3D &	t,
		const novaddt::DAQHit &	h
	)

Definition at line 78 of file Calibration_module.cc.

References dir, novaddt::Track3D::End(), novaddt::DAQHit::Plane(), elec2geo::pos, novaddt::Track3D::Start(), novaddt::Plane::val, novaddt::View::val, and novaddt::DAQHit::View().

Referenced by novaddt::Calibration::analyze().

                                                            {
    TVector3 dir=(t.End()-t.Start());
    dir*=1./dir.Z();
    TVector3 pos=t.Start()+dir*(h.Plane().val-t.Start().Z());
/*    std::cout<<"Z: "<<t.Start().Z()<<" - "<<t.End().Z()
             <<"pos["<<h.Plane().val<<"] = "
             <<pos.X()<<", "<<pos.Y()<<", "<<pos.Z()
    <<std::endl;*/
    if(h.View().val==1)return pos.Y();
    if(h.View().val==2)return pos.X();
    return -999;
}