TPCExpecteddEdxCorrection Class Reference

#include <TPCExpecteddEdxCorrection.hxx>

Inheritance diagram for TPCExpecteddEdxCorrection:

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 correction uses a new parameterisation of the TPC dEdx lines to re-calculate the TPC pulls. The parameterisation is different to that use in tpcRecon, and comes from work Claudio Giganti did to improve the Data/MC comparisons for Production 5.

Definition at line 11 of file TPCExpecteddEdxCorrection.hxx.

Member Function Documentation

Apply()

void TPCExpecteddEdxCorrection::Apply ( AnaSpillC &  spill )

virtual

Apply the TPCExpecteddEdxCorrection. Updates the expected dEdx for all particle hypotheses, and recomputes the pulls based on these values.

Reimplemented from CorrectionBase.

Definition at line 11 of file TPCExpecteddEdxCorrection.cxx.

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

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

  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) {
          tpcTrack->dEdxexpEle = anaUtils::ExpectedTPCdEdx(*tpcTrack, "electron");
          tpcTrack->dEdxexpMuon = anaUtils::ExpectedTPCdEdx(*tpcTrack, "muon");
          tpcTrack->dEdxexpProton = anaUtils::ExpectedTPCdEdx(*tpcTrack, "proton");
          tpcTrack->dEdxexpPion = anaUtils::ExpectedTPCdEdx(*tpcTrack, "pion");
                tpcTrack->dEdxexpKaon = anaUtils::ExpectedTPCdEdx(*tpcTrack, "kaon");

        
          //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/TPCExpecteddEdxCorrection.hxx
• /hep/T2K/nd280rep/ANALYSIS/HIGHLAND2/ANAPARTICLE/highland2/highlandCorrections/v2r8/src/TPCExpecteddEdxCorrection.cxx