| #include "hpidalgstandcuts.h" |
HPidAlgStandCuts
class description - source file - inheritance tree (.pdf)
public:
HPidAlgStandCuts(TString outFileName, Option_t* opt_algorithms = "RICHRINGCUTS,TOFPVSBETA,TOFINOPVSBETA,SHOWERMAXFVSP,SHOWERSUM0VSP")
HPidAlgStandCuts(Float_t fWeight, TString outFileName, Option_t* opt_algorithms = "RICHRINGCUTS,TOFPVSBETA,TOFINOPVSBETA,SHOWERMAXFVSP,SHOWERSUM0VSP")
HPidAlgStandCuts(const HPidAlgStandCuts&)
~HPidAlgStandCuts()
Bool_t buildOutNtuple()
virtual Bool_t calculateCL(HPidTrackCand* pTrack, Float_t* afReturn, Short_t& nCatIndex)
virtual Bool_t calculatePDF(HPidTrackCand* pTrack, Float_t* afReturn, Short_t& nCatIndex)
static TClass* Class()
void clear()
Bool_t fillOutNtuple(Int_t iPid, HPidTrackCandSim* pTrack)
virtual Bool_t finalize()
TString getOutFileName()
const HPidAlgStandCutsPar* getPidParams() const
virtual Bool_t init()
virtual TClass* IsA() const
Bool_t isLepInMdc(Int_t iPid, HPidHitData* pHitData, HPidTrackData* pTrackData)
Bool_t isLepInRich(Int_t iPid, HPidHitData* pHitData)
Bool_t isLepInShower(Int_t iPid, HPidHitData* pHitData, HPidTrackData* pTrackData)
Bool_t isLepInTof(Int_t iPid, HPidHitData* pHitData, HPidTrackData* pTrackData)
Bool_t isLepInTofino(Int_t iPid, HPidHitData* pHitData, HPidTrackData* pTrackData)
Bool_t openOutFile()
HPidAlgStandCuts& operator=(const HPidAlgStandCuts&)
virtual void print() const
virtual Bool_t reinit()
void setAlgorithms(Option_t* opt_algorithms = "RICHRINGCUTS,TOFPVSBETA,TOFINOPVSBETA,SHOWERMAXFVSP,SHOWERSUM0VSP")
void setDefaults()
void setOutFileName(TString outFile)
void setPidParams(HPidAlgStandCutsPar* pP)
virtual void ShowMembers(TMemberInspector& insp, char* parent)
virtual void Streamer(TBuffer& b)
void StreamerNVirtual(TBuffer& b)
Bool_t writeOutFile()
private:
HPidAlgStandCutsPar* pParams
TString sOutFileName name of the output file with ntuple
TFile* pOutFile output file
TNtuple* pNtuple control ntuple
TF1* pMdcDiscriminator Function used to decide upon leptonic id based upon energy loss in mdc
Bool_t bSim check experiment or simulation
Int_t iSimEventFlag simulation event structure flag
Bool_t bRich lepton in Rich
Bool_t bMdc lepton in Mdc
Bool_t bShower lepton in Shower
Bool_t bTofino lepton in Tofino
Bool_t bTof lepton in Tof
Bool_t bUseRichRingCuts use Rich ring condition
Bool_t bUseMdcDeDx use Mdc de/dx condition
Bool_t bUseTofPvsBeta use P vs Beta in Tof region
Bool_t bUseTofinoPvsBeta use P vs Beta in Tofiono region
Bool_t bUseShowerMaxFvsP use Shower maxF vs P algorithm
Bool_t bUseShowerSum0vsP use Shower Sum0 vs P algorithm
Bool_t bUseShowerSumDiffvsP use Shower SumDiff vs P algorithm
HPidAlgStandCuts
Test algorithm for the PID with the standard cuts (windows):
1. RICH-rings (PM, AC, NP, RC)
2. Beta vs Momentum windows in the system 0 and system 1
3. Pre-Shower algorithm (F vs momentum, Sum0 vs momentum)
4. Pre-Shower algorithm (sum2+sum1-sum0 vs momentum)
In case of lepton in the window the pdf vector element is set to 1.
Otherwise, it is set to 0
HPidAlgStandCuts(TString outFileName, Option_t *opt_algorithms)
: HPidAlgorithm("StandCuts", algStandCuts)
constructor
HPidAlgStandCuts(Float_t fWeight,TString outFileName,Option_t *opt_algorithms)
: HPidAlgorithm("StandCuts", algStandCuts, fWeight)
constructor with weight
void setDefaults(void)
default values
void clear(void)
clear
void setAlgorithms(Option_t *opt_algorithms)
Bool_t init(void)
Bool_t openOutFile()
open output ntuple file
Bool_t buildOutNtuple()
bild output ntuple for sim/exp
Bool_t writeOutFile(void)
write output ntuple to file
Bool_t reinit(void)
Check if parameter context corresponds to the
appropriate momentum algorithm index
If not -> return kFALSE (analysis can not be done in such case)
and print error message
Bool_t finalize(void)
write output file
Bool_t isLepInRich(Int_t iPid,HPidHitData *pHitData)
Check RICH ring quality by using standard cuts
Bool_t isLepInShower(Int_t iPid,HPidHitData *pHitData,HPidTrackData *pTrackData)
Check whether lepton fulfills PreShower conditions
F(p)- threshold, Sum0(p)- upper limit
Bool_t isLepInTofino(Int_t iPid,HPidHitData *pHitData, HPidTrackData *pTrackData)
Check whether lepton is in beta vs momentum window (system 0)
Bool_t isLepInMdc(Int_t iPid,HPidHitData *pHitData, HPidTrackData *pTrackData)
Check whether this track has a lepton-like energy loss
Bool_t isLepInTof(Int_t iPid,HPidHitData *pHitData, HPidTrackData *pTrackData)
Check whether lepton is in beta vs momentum window (system 1)
Bool_t fillOutNtuple(Int_t pid,HPidTrackCandSim *pTrack)
fill ntuple in experiment and in simulation
Bool_t calculatePDF(HPidTrackCand *pTrack,
Float_t afReturn[], Short_t &nCatIndex)
Calculate PDF from data (return 1 in case of lepton in the window)
Bool_t calculateCL(HPidTrackCand *pTrack,
Float_t afReturn[], Short_t &nCatIndex)
void print(void) const
Empty function for printing info of the algorithm
Inline Functions
void ~HPidAlgStandCuts()
void setPidParams(HPidAlgStandCutsPar* pP)
const HPidAlgStandCutsPar* getPidParams() const
void setOutFileName(TString outFile)
TString getOutFileName()
TClass* Class()
TClass* IsA() const
void ShowMembers(TMemberInspector& insp, char* parent)
void Streamer(TBuffer& b)
void StreamerNVirtual(TBuffer& b)
HPidAlgStandCuts HPidAlgStandCuts(const HPidAlgStandCuts&)
HPidAlgStandCuts& operator=(const HPidAlgStandCuts&)
Author: Jacek Otwinowski 22/02/2005
Last update: Jacek Otwinowski 15/03/2005
ROOT page - Class index - Class Hierarchy - Top of the page
This page has been automatically generated. If you have any comments or suggestions about the page layout send a mail to ROOT support, or contact the developers with any questions or problems regarding ROOT.