FGDHybridTrackEffSystematics Class Reference

Inheritance diagram for FGDHybridTrackEffSystematics:

Public Member Functions
	FGDHybridTrackEffSystematics (bool computecounters=false)
Weight_h	ComputeWeight (const ToyExperiment &toy, const AnaEventC &event, const ToyBoxB &box)
Public Member Functions inherited from EventWeightBase
	EventWeightBase ()
	Create the systematic, with one parameter.
	EventWeightBase (UInt_t npar)
	Create the systematic, specifying the number of systematic parameters.
virtual	~EventWeightBase ()
	Everyone should have a destructor.
virtual void	InitializeEvent (const AnaEventC &event, const SelectionBase &sel, Int_t ibranch)
	Initialize the SystBox for this event.
virtual Weight_h	ComputeWeight (const ToyExperiment &toy, const AnaEventC &event, const ToyBoxB &box, const SelectionBase &sel)
	This is now the actual method called by SystematicManager, which allows further selection tunning of relevant objects once the ToyBox is filled.
Public Member Functions inherited from SystematicBase
	SystematicBase ()
	Create the systematic, with one parameter.
	SystematicBase (UInt_t npar)
	Create the systematic, specifying the number of systematic parameters.
virtual	~SystematicBase ()
	Everyone should have a destructor.
UInt_t	GetNParameters () const
	Returns the number of systematic parameters associated to this systematic.
void	SetNParameters (int N)
	Set the number of systematic parameters associated to this systematic.
virtual const char *	GetName () const
	Return the name of this systematic. This overrides the TObject::GetName() interface.
const std::string &	Name () const
	Return the name of this systematic.
void	SetName (const std::string &name)
	Set the name of this systematic.
Int_t	GetIndex () const
	Return the index of this systematic.
void	SetIndex (Int_t index)
	Set the index of this systematic.
TypeEnum	Type () const
	Return the type of this systematic.
std::string	ConvertType ()
	Return the type of this systematic.
void	SetType (TypeEnum type)
	Set the type of this systematic.
PDFEnum	PDF () const
	Return the PDF of this systematic.
std::string	ConvertPDF ()
	Return the PDF of this systematic.
void	SetPDF (PDFEnum PDF)
	Set the PDF of this systematic.
void	SetEnabled (bool ena)
	Set the enable bit.
bool	IsEnabled () const
	Returns the enable bit.
virtual void	FinalizeEvent (const AnaEventC &event)
	Delete the SystBox for this event.
SystBoxB *	GetSystBox (const AnaEventC &event, Int_t isel=0, Int_t ibranch=0) const
	Get the SystBox corresponding to a selection, branch and event.
virtual void	Initialize (Int_t nsel, Int_t isel, Int_t nbranch, Int_t nevents)
	Create the array of SystBox.
	ClassDef (SystematicBase, 2)

Public Attributes
BinnedParams *	_fgd1
BinnedParams *	_fgd2

Protected Member Functions
bool	IsRelevantTrueObject (const AnaEventC &event, const AnaTrueObjectC &trueTrack) const
	Is this true track relevant for this systematic ?
Int_t	GetRelevantRecObjectGroups (const SelectionBase &sel, Int_t ibranch, Int_t *IDs) const
	Get the TrackGroup IDs array for this systematic.
Int_t	GetRelevantTrueObjectGroups (const SelectionBase &sel, Int_t ibranch, Int_t *IDs) const
	Get the TrueTrackGroup IDs array for this systematic.
Protected Member Functions inherited from EventWeightBase
virtual bool	IsRelevantRecObject (const AnaEventC &, const AnaRecObjectC &) const
	Check whether a AnaRecObject is relevant for this systematic or not.
virtual void	FillSystBox (const AnaEventC &event, const SelectionBase &sel, Int_t ibranch)
	Fills the SystBox.
virtual bool	CheckTrueRecoAssociation (const AnaTrueObjectC &trueRecObject, const AnaRecObjectC &track) const
	Criteria for true-reco association for this systematic.
void	FillTrueRecoAssociation (SystBoxB &box)
	Fll the true-reco association array.

Protected Attributes
bool	_computecounters
Protected Attributes inherited from SystematicBase
std::string	_name
	The name of this systematic.
Int_t	_index
	The index of this systematic (needed by SystematicsManager);.
TypeEnum	_type
	The type of this systematic (variation, weight or flux)
bool	_enabled
	Is this systematic enabled ?
PDFEnum	_PDF
	The PDF use for the systematic parameter scan.
UInt_t	_nParameters
	the number of systematic parameters
SystBoxB ****	_systBoxes
	----—— Relevant objects for this systematic ------------—— More...
Int_t	_systBoxesNSel
Int_t	_systBoxesNBranches [NMAXSELECTIONS]
Int_t	_systBoxesNEvents [NMAXSELECTIONS][NMAXBRANCHES]

Additional Inherited Members
Public Types inherited from SystematicBase
enum	TypeEnum { kVariation = 0, kWeight, kFlux, kUnknown }
	Enumerator describing the values that _type member can take.
enum	PDFEnum {   kGaussian = 0, kUniform, kBinomial, kMultinomial,   kUnknownPDF }
	Enumerator describing the values that _PDF member can take.

Detailed Description

Definition at line 8 of file FGDHybridTrackEffSystematics.hxx.

Member Function Documentation

ComputeWeight()

Weight_h FGDHybridTrackEffSystematics::ComputeWeight ( const ToyExperiment &  toy, const AnaEventC &  event, const ToyBoxB &  box  )

virtual

This is now the actual method called by SystematicManager, which allows further selection tunning of relevant objects once the ToyBox is filled This MUST be overridden in the derived class.

Implements EventWeightBase.

Definition at line 26 of file FGDHybridTrackEffSystematics.cxx.

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

  Weight_h eventWeight=1; 
  BinnedParamsParams params;

  // Cast the ToyBox to the appropriate type
  const ToyBoxTracker& box = *static_cast<const ToyBoxTracker*>(&boxB); 

  // Get the SystBox for this event, and the appropriate selection and branch
  SystBoxB* SystBox = GetSystBox(event,box.SelectionEnabledIndex,box.SuccessfulBranch);

  if(_computecounters){
    _fgd1->InitializeEfficiencyCounter();
    _fgd2->InitializeEfficiencyCounter();
  }
  // Get the direction of the main track
  AnaTrackB *maintrack = box.MainTrack; 
  if(!maintrack)  return eventWeight;

  for (Int_t itrue=0;itrue< SystBox->nRelevantTrueObjects; itrue++){      
    AnaTrueParticleB* trueTrack = static_cast<AnaTrueParticleB*>(SystBox->RelevantTrueObjects[itrue]);            

    // TODO by now the reco track associated to the true is the first one found, we should use instead a criteria to decide.
    // by now since it's very rare it's like that.
    // retrieve the reconstructed track associated
    AnaTrackB* recoTrack = static_cast<AnaTrackB*>(SystBox->RelevantTrueObjectsReco[itrue]);

    // Get the cos (theta), where theta is the angle between the main track and the FGD track
    Float_t costheta = 0;
    if (recoTrack)
      costheta = anaUtils::ScalarProduct(maintrack->DirectionStart, recoTrack->DirectionStart);
    else
      costheta = anaUtils::ScalarProduct(maintrack->DirectionStart, trueTrack->Direction);

    Float_t pdg=(Float_t)trueTrack->PDG;
    bool found = (recoTrack);
    
    // Get the FGD tracking efficiency
    Int_t index;
    if(box.DetectorFV==SubDetId::kFGD1){
      if(!_fgd1->GetBinValues(fabs(pdg),costheta, params, index)) continue;
      // save the index of the varied parameters into the systematics box
      if(_computecounters)
        _fgd1->UpdateEfficiencyCounter(index,found);    
    }
    else if(box.DetectorFV==SubDetId::kFGD2){
      if(!_fgd2->GetBinValues(fabs(pdg),costheta, params, index)) continue;
      if(_computecounters)
        _fgd2->UpdateEfficiencyCounter(index,found);    
      
    }
    else return eventWeight;
    
    eventWeight *= systUtils::ComputeEffLikeWeight(true, toy, GetIndex(), index, params);
  }

  return eventWeight;

}

---

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

• /hep/T2K/nd280rep/ANALYSIS/HIGHLAND2/ANAPARTICLE/psyche/psycheSystematics/v3r12/src/FGDHybridTrackEffSystematics.hxx
• /hep/T2K/nd280rep/ANALYSIS/HIGHLAND2/ANAPARTICLE/psyche/psycheSystematics/v3r12/src/FGDHybridTrackEffSystematics.cxx