TPCdEdxDataCorrection Class Reference

#include <TPCdEdxDataCorrection.hxx>

Inheritance diagram for TPCdEdxDataCorrection:

Public Member Functions
	TPCdEdxDataCorrection ()
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 data, and is based on the run/subrun of this spill.

The corrections are defined in the baseAnalysis/data/TPCdEdxCorrection.dat. The format of this file is

run subrun corr_tpc1 err_tpc1 corr_tpc2 err_tpc2 corr_tpc3 err_tpc3

where corr_tpcN is the CT correction to apply to TPC N, and err_tpcN is the error on that correction.

Definition at line 21 of file TPCdEdxDataCorrection.hxx.

Constructor & Destructor Documentation

TPCdEdxDataCorrection()

TPCdEdxDataCorrection::TPCdEdxDataCorrection

(

)

Constructor, which reads the corrections from the text file, and fills the member variables with the details.

Definition at line 7 of file TPCdEdxDataCorrection.cxx.

                                            :CorrectionBase(){
//********************************************************************

  char filename[300];
  if (versionUtils::prod6_corrections)
    sprintf(filename, "%s/data/TPCdEdxCorrection_p6B.dat", getenv("HIGHLANDCORRECTIONSROOT"));
  else
    sprintf(filename, "%s/data/TPCdEdxCorrection_p5F.dat", getenv("HIGHLANDCORRECTIONSROOT"));
  

  std::cout << " CT correction data " << filename << std::endl;
  FILE *pFile = fopen(filename, "r");

  if (pFile == NULL) {
    printf("Cannot open file.\n");
    exit(1);
  }

  int run;
  int subrun;

  float c[3]; // CT_expected -CT_measured
  float ec[3]; // error on CT_expected - CT_measured

  while (fscanf(pFile, "%d%d%f%f%f%f%f%f", &run, &subrun, &c[0], &ec[0], &c[1], &ec[1], &c[2], &ec[2]) == 8) {
    comb = std::make_pair(run, subrun);
    corrTPC1[comb] = c[0];
    corrTPC2[comb] = c[1];
    corrTPC3[comb] = c[2];
  }

  itr = corrTPC1.begin();
  fclose(pFile);
}

Member Function Documentation

Apply()

void TPCdEdxDataCorrection::Apply ( AnaSpillC &  spill )

virtual

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

Reimplemented from CorrectionBase.

Definition at line 80 of file TPCdEdxDataCorrection.cxx.

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

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

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

  Int_t subrun = spill.EventInfo->SubRun;

  for (unsigned int i = 0; i < spill.Bunches.size(); i++) {
    AnaBunch* bunch = static_cast<AnaBunch*>(spill.Bunches[i]);
    AnaTrackB* allTpcTracks[100];
    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.
          tpcTrack->dEdxMeas = CT0 / GetCorrection(spill.EventInfo->Run, subrun, tpcTrack->Detector);
          
          //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/TPCdEdxDataCorrection.hxx
• /hep/T2K/nd280rep/ANALYSIS/HIGHLAND2/ANAPARTICLE/highland2/highlandCorrections/v2r8/src/TPCdEdxDataCorrection.cxx