#include <FlipKinematicsCorrection.hxx>

Inheritance diagram for FlipKinematicsCorrection:

Public Member Functions
	FlipKinematicsCorrection (SubDetId::SubDetEnum det=SubDetId::kFGD)

void	Apply (AnaSpillC &spill)

Public Member Functions inherited from CorrectionBase
	CorrectionBase (const CorrectionBase &corr)
	Copy constructor.

virtual	~CorrectionBase ()
	Everyone should have a destructor.

virtual const char *	GetName () const
	Return the name of this correction. This overrides the TObject::GetName() interface.

const std::string &	Name () const
	Return the name of this correction.

void	SetName (const std::string &name)
	Set the name of this correction.

void	Enable ()
	Enable the correction.

void	Disable ()
	Disable the correction.

bool	IsEnabled () const
	Is the correction enabled.

void	SetAppliedInInput (bool ap)
	Set the correction as applied in the input file.

bool	IsAppliedInInput () const
	Is the correction already applied in the input file ?

Int_t	GetIndex () const
	Return the index of this correction.

void	SetIndex (Int_t index)
	Set the index of this correction.

	ClassDef (CorrectionBase, 2)

Protected Member Functions
bool	IsRelevantTrack (const AnaTrackB &track) const

Protected Attributes
SubDetId::SubDetEnum	_det

Protected Attributes inherited from CorrectionBase
std::string	_name
	The name of this correction.

bool	_enabled
	Is the correction enabled?

bool	_appliedInInput
	Is the correction already applied in the input file ?

Int_t	_index
	The index of the correction.

Detailed Description

Each track that has a TPC component present together with a segment in other detector should have an "alternate" re-fit with a reversed charge and direction Sometimes it happens that this reverse fit fails on oaRecon level, hence the info is not available, this correction "fixes" the problem by taking the closest (to the end point) TPC momentum and applying a "reversed" dEdX correction up to the end point

Definition at line 17 of file FlipKinematicsCorrection.hxx.

Member Function Documentation

§ Apply()

void FlipKinematicsCorrection::Apply ( AnaSpillC & spill )

virtual

This is the function that applies the correction to the input data. This MUST be overridden in the derived class.

Reimplemented from CorrectionBase.

Definition at line 19 of file FlipKinematicsCorrection.cxx.

                                                       {
 //********************************************************************
 
 
   AnaSpill& spill = *static_cast<AnaSpill*>(&spillC);
 
   for (UInt_t i = 0; i < spill.Bunches.size(); i++) {
     AnaBunch* bunch = static_cast<AnaBunch*>(spill.Bunches[i]);
     for (UInt_t j = 0; j < bunch->Particles.size(); j++) {
 
       AnaTrack* track = static_cast<AnaTrack*>(bunch->Particles[j]);
 
       if (!track) continue;
 
       // Check whether a track is relevant to apply the correction
       if (!IsRelevantTrack(*track)) continue;
 
       AnaTPCParticle* tpcTrack = static_cast<AnaTPCParticle*>(anaUtils::GetSegmentWithMostNodesInClosestTpc(*track, track->PositionEnd));
 
       if (!tpcTrack) continue;
 
       if (!cutUtils::TPCTrackQualityCut(*tpcTrack)) continue;
 
 #ifdef  UseRecPack
       track->MomentumFlip = -999.;
       double momentum_tmp = tpcTrack->MomentumEnd; 
       if (ND::tman().CorrectEnergyLossLinear(tpcTrack->PositionEnd, track->PositionEnd, momentum_tmp, 
             ParticleId::GetParticle(track->ReconPDG), AnaRecPackManager::kBwd)) // from closest TPC endpoint and to track end backwards
         track->MomentumFlip = momentum_tmp;
     
 
 #else
       std::cout << " WARNING! FlipKinematicsCorrection is run w/o highlandRecPack usage. \n"
         " It is highly recommended to use RecPack to restore possible fit failures of the \n"
         " global recon re-fits with opposite track sense hypotheis assumes (end to start) " << std::endl;
 
 #endif
 
     } 
   }       
 }

§ IsRelevantTrack()

bool FlipKinematicsCorrection::IsRelevantTrack ( const AnaTrackB & track ) const

protected

A correction is supposed to be relevant for the tracks that have a potential to be reversed (should have a TOF) and for which B field + eloss momentum estimation is crucial: TPC should be present

Definition at line 62 of file FlipKinematicsCorrection.cxx.

                                                                           {
   //********************************************************************
 
   if (!SubDetId::GetDetectorUsed(track.Detector, _det)) return false;
 
   if (!SubDetId::GetDetectorUsed(track.Detector, SubDetId::kTPC)) return false;
 
   if (ToFSenseCorrector::GetToFTopology(track) == ToFSenseCorrector::kToFInvalid) return false;
         
   return true;
 
 }

Member Data Documentation

§ _det

SubDetId::SubDetEnum FlipKinematicsCorrection::_det

protected

A detector bit that defines which ToFs should be present to consider the object for correction

Definition at line 35 of file FlipKinematicsCorrection.hxx.

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

/hep/T2K/nd280rep/ANALYSIS/HIGHLAND2/ANAPARTICLE/highland2/highlandCorrections/v2r8/src/FlipKinematicsCorrection.hxx
/hep/T2K/nd280rep/ANALYSIS/HIGHLAND2/ANAPARTICLE/highland2/highlandCorrections/v2r8/src/FlipKinematicsCorrection.cxx

Public Member Functions