TPCdEdxMCCorrection Class Reference

#include <TPCdEdxMCCorrection.hxx>

Inheritance diagram for TPCdEdxMCCorrection:

Public Member Functions
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)

Additional Inherited Members
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

This class defines a correction that affects the measured dE/dx of an AnaTpcSegment. This in turn affects the computed pulls (how electron-like a track is, for example).

This correction is applied to MC, and currently has NO EFFECT.

Definition at line 12 of file TPCdEdxMCCorrection.hxx.

Member Function Documentation

Apply()

void TPCdEdxMCCorrection::Apply ( AnaSpillC &  spill )

virtual

Apply the dE/dx correction to all the AnaTpcSegment members of all tracks in the spill. See TPCdEdxMCCorrection class documentation for details.

Reimplemented from CorrectionBase.

Definition at line 30 of file TPCdEdxMCCorrection.cxx.

                                                 {
//********************************************************************

  AnaSpill& spill = *static_cast<AnaSpill*>(&spillBB);

  // No correction for data
  if (!spill.GetIsMC())
    return;

  for (unsigned int i = 0; i < spill.Bunches.size(); i++) {
    AnaTrackB* allTpcTracks[100];
    AnaBunch* bunch = static_cast<AnaBunch*>(spill.Bunches[i]);
    int nTPC = anaUtils::GetAllTracksUsingDet(*bunch, SubDetId::kTPC,  allTpcTracks);
    for (Int_t j = 0; j < nTPC; j++) {
      for (int k = 0; k < allTpcTracks[j]->nTPCSegments; k++) {

        // The raw TPC track
        AnaTPCParticleB* original = static_cast<const AnaTrackB*>(allTpcTracks[j]->Original)->TPCSegments[k];
        // The corrected TPC track
        AnaTPCParticle* tpcTrack = static_cast<AnaTPCParticle*>(allTpcTracks[j]->TPCSegments[k]);

        // Get the raw CT
        double CT0 = original->dEdxMeas;

        // Apply the correction
        if (CT0 != 999999) {
          // Apply correction only if the CT is a valid number.
          int pdg = 0;
          if (tpcTrack->TrueObject) {
            pdg = tpcTrack->GetTrueParticle()->PDG;
          }

          tpcTrack->dEdxMeas = CT0 / GetCorrection(pdg);
 
          //recompute TPC pulls
          Float_t pulls[4];
          anaUtils::ComputeTPCPull(*tpcTrack, pulls);
          
          tpcTrack->Pullele = pulls[1];
          tpcTrack->Pullmu  = pulls[0];
          tpcTrack->Pullp   = pulls[2];
          tpcTrack->Pullpi  = pulls[3];
          
          //for the kaon need to recompute explicitely   
          tpcTrack->Pullk = ((tpcTrack->dEdxMeas - tpcTrack->dEdxexpKaon) / tpcTrack->dEdxSigmaKaon); 
        
        }
      }
    }
  }
}

---

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

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