class AliCaloTrackMCReader: public AliCaloTrackReader
Class for reading data (Kinematics) in order to do prompt gamma or other particle identification and correlations Separates generated particles into charged (CTS) and neutral (PHOS or EMCAL acceptance) Now, it only works with data stored in Kinematics.root and not in filtered Kinematics branch in AODs *-- Author: Gustavo Conesa (LNF-INFN)
Function Members (Methods)
public:
protected:
| virtual void | TObject::DoError(int level, const char* location, const char* fmt, va_list va) const |
| void | TObject::MakeZombie() |
private:
| AliCaloTrackMCReader(const AliCaloTrackMCReader& r) | |
| AliCaloTrackMCReader& | operator=(const AliCaloTrackMCReader& r) |
Data Members
public:
| enum AliCaloTrackReader::inputDataType { | kESD | |
| kAOD | ||
| kMC | ||
| }; | ||
| enum AliCaloTrackReader::detector { | kEMCAL | |
| kPHOS | ||
| kCTS | ||
| kDCAL | ||
| kDCALPHOS | ||
| }; | ||
| enum TObject::EStatusBits { | kCanDelete | |
| kMustCleanup | ||
| kObjInCanvas | ||
| kIsReferenced | ||
| kHasUUID | ||
| kCannotPick | ||
| kNoContextMenu | ||
| kInvalidObject | ||
| }; | ||
| enum TObject::[unnamed] { | kIsOnHeap | |
| kNotDeleted | ||
| kZombie | ||
| kBitMask | ||
| kSingleKey | ||
| kOverwrite | ||
| kWriteDelete | ||
| }; |
protected:
| TList* | AliCaloTrackReader::fAODBranchList | -> List with AOD branches created and needed in analysis |
| TArrayI | AliCaloTrackReader::fAcceptEventsWithBit | Accept events if trigger bit is on |
| Bool_t | AliCaloTrackReader::fAcceptFastCluster | Accept events from fast cluster, exclude these events for LHC11a |
| Bool_t | AliCaloTrackReader::fAcceptOnlyHIJINGLabels | Select clusters or tracks that where generated by HIJING, reject other generators in case of cocktail |
| Bool_t | AliCaloTrackReader::fAccessTrackTOF | Access the track TOF, in case of problems when accessing GetTOFBunchCrossing |
| Bool_t | AliCaloTrackReader::fAnaLED | Analyze LED data only. |
| AliAODJetEventBackground* | AliCaloTrackReader::fBackgroundJets | ! background jets |
| Int_t | AliCaloTrackReader::fBitEGA | Trigger bit on VCaloTrigger for EGA |
| Int_t | AliCaloTrackReader::fBitEJE | Trigger bit on VCaloTrigger for EJE |
| Float_t | AliCaloTrackReader::fCTSPtMax | pT Threshold on charged particles |
| Float_t | AliCaloTrackReader::fCTSPtMin | pT Threshold on charged particles |
| TObjArray* | AliCaloTrackReader::fCTSTracks | -> temporal array with tracks |
| AliCalorimeterUtils* | AliCaloTrackReader::fCaloUtils | Pointer to CalorimeterUtils |
| Int_t | AliCaloTrackReader::fCentralityBin[2] | Minimum and maximum value of the centrality for the analysis |
| TString | AliCaloTrackReader::fCentralityClass | Name of selected centrality class |
| Int_t | AliCaloTrackReader::fCentralityOpt | Option for the returned value of the centrality, possible options 5, 10, 100 |
| Bool_t | AliCaloTrackReader::fCheckFidCut | Do analysis for clusters in defined region |
| Bool_t | AliCaloTrackReader::fComparePtHardAndClusterPt | In MonteCarlo, jet events, reject events with too large cluster energy |
| Bool_t | AliCaloTrackReader::fComparePtHardAndJetPt | In MonteCarlo, jet events, reject fake events with wrong jet energy. |
| Bool_t | AliCaloTrackReader::fCorrectELinearity | Correct cluster linearity, always on |
| TObjArray* | AliCaloTrackReader::fDCALClusters | -> temporal array with DCAL CaloClusters, not needed in the normal case, use just EMCal array with DCal limits |
| Int_t | AliCaloTrackReader::fDataType | Select MC:Kinematics, Data:ESD/AOD, MCData:Both |
| Int_t | AliCaloTrackReader::fDebug | Debugging level |
| TString | AliCaloTrackReader::fDeltaAODFileName | Delta AOD file name |
| Bool_t | AliCaloTrackReader::fDoPileUpEventRejection | Select pile-up events by SPD |
| Bool_t | AliCaloTrackReader::fDoRejectNoTrackEvents | Reject events with no selected tracks in event |
| Bool_t | AliCaloTrackReader::fDoV0ANDEventSelection | Select events depending on V0AND |
| Bool_t | AliCaloTrackReader::fDoVertexBCEventSelection | Select events with vertex on BC=0 or -100 |
| AliVCaloCells* | AliCaloTrackReader::fEMCALCells | ! temporal array with EMCAL CaloCells |
| TObjArray* | AliCaloTrackReader::fEMCALClusters | -> temporal array with EMCAL CaloClusters |
| TString | AliCaloTrackReader::fEMCALClustersListName | Alternative list of clusters produced elsewhere and not from InputEvent |
| Double_t | AliCaloTrackReader::fEMCALParamTimeCutMax[4] | Remove clusters/cells with time larger than parametrized value, in ns |
| Float_t | AliCaloTrackReader::fEMCALParamTimeCutMin[4] | Remove clusters/cells with time smaller than parametrized value, in ns |
| Float_t | AliCaloTrackReader::fEMCALPtMax | pT Threshold on emcal clusters |
| Float_t | AliCaloTrackReader::fEMCALPtMin | pT Threshold on emcal clusters |
| Double_t | AliCaloTrackReader::fEMCALTimeCutMax | Remove clusters/cells with time larger than this value, in ns |
| Double_t | AliCaloTrackReader::fEMCALTimeCutMin | Remove clusters/cells with time smaller than this value, in ns |
| Int_t | AliCaloTrackReader::fEMCalBCEvent[19] | Fill one entry per event if there is a cluster in a given BC |
| Int_t | AliCaloTrackReader::fEMCalBCEventCut[19] | Fill one entry per event if there is a cluster in a given BC, depend on cluster E, acceptance cut |
| Int_t | AliCaloTrackReader::fEventNumber | Event number |
| TString | AliCaloTrackReader::fEventPlaneMethod | Name of event plane method, by default "Q" |
| Bool_t | AliCaloTrackReader::fEventTrigCentral | Event is AliVEvent::kCentral on its name, it should correspond to PbPb |
| Bool_t | AliCaloTrackReader::fEventTrigEMCALL0 | Event is EMCal L0 on its name, it should correspond to AliVEvent::kEMC7, AliVEvent::kEMC1 |
| Bool_t | AliCaloTrackReader::fEventTrigEMCALL1Gamma1 | Event is L1-Gamma, threshold 1 on its name, it should correspond kEMCEGA |
| Bool_t | AliCaloTrackReader::fEventTrigEMCALL1Gamma2 | Event is L1-Gamma, threshold 2 on its name, it should correspond kEMCEGA |
| Bool_t | AliCaloTrackReader::fEventTrigEMCALL1Jet1 | Event is L1-Gamma, threshold 1 on its name, it should correspond kEMCEGA |
| Bool_t | AliCaloTrackReader::fEventTrigEMCALL1Jet2 | Event is L1-Gamma, threshold 2 on its name, it should correspond kEMCEGA |
| Bool_t | AliCaloTrackReader::fEventTrigMinBias | Event is min bias on its name, it should correspond to AliVEvent::kMB, AliVEvent::kAnyInt |
| Bool_t | AliCaloTrackReader::fEventTrigSemiCentral | Event is AliVEvent::kSemiCentral on its name, it should correspond to PbPb |
| Bool_t | AliCaloTrackReader::fEventTriggerAtSE | select triggered event at SE base task or here |
| UInt_t | AliCaloTrackReader::fEventTriggerMask | select this triggerered event |
| AliFiducialCut* | AliCaloTrackReader::fFiducialCut | Acceptance cuts |
| Bool_t | AliCaloTrackReader::fFillCTS | use data from CTS |
| Bool_t | AliCaloTrackReader::fFillDCAL | use data from DCAL, not needed in the normal case, use just EMCal array with DCal limits |
| Bool_t | AliCaloTrackReader::fFillEMCAL | use data from EMCAL |
| Bool_t | AliCaloTrackReader::fFillEMCALCells | use data from EMCAL |
| Bool_t | AliCaloTrackReader::fFillInputBackgroundJetBranch | Flag to use data from background jets |
| Bool_t | AliCaloTrackReader::fFillInputNonStandardJetBranch | Flag to use data from non standard jets |
| Bool_t | AliCaloTrackReader::fFillPHOS | use data from PHOS |
| Bool_t | AliCaloTrackReader::fFillPHOSCells | use data from PHOS |
| TString | AliCaloTrackReader::fFiredTriggerClassName | Name of trigger event type used to do the analysis |
| TString | AliCaloTrackReader::fInputBackgroundJetBranchName | Name of background jet branch |
| AliVEvent* | AliCaloTrackReader::fInputEvent | ! pointer to esd or aod input |
| TString | AliCaloTrackReader::fInputNonStandardJetBranchName | Name of non standard jet branch |
| Bool_t | AliCaloTrackReader::fIsBadCellEvent | Bad cell triggered event flag, any cell in cluster is bad |
| Bool_t | AliCaloTrackReader::fIsBadMaxCellEvent | Bad cell triggered event flag, only max energy cell is bad |
| Bool_t | AliCaloTrackReader::fIsExoticEvent | Exotic trigger event flag |
| Bool_t | AliCaloTrackReader::fIsTriggerMatch | Could match the event to a trigger patch? |
| Bool_t | AliCaloTrackReader::fIsTriggerMatchOpenCut[3] | Could not match the event to a trigger patch?, retry opening cuts |
| Int_t | AliCaloTrackReader::fLastMixedCaloEvent | Temporary container with the last event added to the mixing list for photons |
| Int_t | AliCaloTrackReader::fLastMixedTracksEvent | Temporary container with the last event added to the mixing list for tracks |
| TList** | AliCaloTrackReader::fListMixedCaloEvents | ! Container for photon stored for different events, used in case of own mixing, set in analysis class |
| TList** | AliCaloTrackReader::fListMixedTracksEvents | ! Container for tracks stored for different events, used in case of own mixing, set in analysis class |
| AliMCEvent* | AliCaloTrackReader::fMC | ! Monte Carlo Event Handler |
| UInt_t | AliCaloTrackReader::fMixEventTriggerMask | select this triggerered event for mixing, tipically kMB or kAnyINT |
| AliMixedEvent* | AliCaloTrackReader::fMixedEvent | ! mixed event object. This class is not the owner |
| TLorentzVector | AliCaloTrackReader::fMomentum | ! Temporal TLorentzVector container, avoid declaration of TLorentzVectors per event |
| Int_t | AliCaloTrackReader::fNMCProducedMax | In case of cocktail, select particles in the list with label up to this value |
| Int_t | AliCaloTrackReader::fNMCProducedMin | In case of cocktail, select particles in the list with label from this value |
| Int_t | AliCaloTrackReader::fNMixedEvent | number of events in mixed event buffer |
| Int_t | AliCaloTrackReader::fNNonPileUpClusters | Number of clusters with time below 20 ns |
| Int_t | AliCaloTrackReader::fNPileUpClusters | Number of clusters with time avobe 20 ns |
| Int_t | AliCaloTrackReader::fNPileUpClustersCut | Cut to select event as pile-up |
| TClonesArray* | AliCaloTrackReader::fNonStandardJets | ! temporal array with jets |
| AliAODEvent* | AliCaloTrackReader::fOutputEvent | ! pointer to aod output |
| AliVCaloCells* | AliCaloTrackReader::fPHOSCells | ! temporal array with PHOS CaloCells |
| TObjArray* | AliCaloTrackReader::fPHOSClusters | -> temporal array with PHOS CaloClusters |
| Float_t | AliCaloTrackReader::fPHOSPtMax | pT Threshold on phos clusters |
| Float_t | AliCaloTrackReader::fPHOSPtMin | pT Threshold on phos clusters |
| Double_t | AliCaloTrackReader::fPileUpParamSPD[5] | Parameters to pass to method IsPileupFromSPD: Int_t minContributors, |
| Float_t | AliCaloTrackReader::fPtHardAndClusterPtFactor | Factor between ptHard and cluster pT to reject/accept event. |
| Float_t | AliCaloTrackReader::fPtHardAndJetPtFactor | Factor between ptHard and jet pT to reject/accept event. |
| Bool_t | AliCaloTrackReader::fReadAODMCParticles | Access kine information from filtered AOD MC particles |
| Bool_t | AliCaloTrackReader::fReadStack | Access kine information from stack |
| Bool_t | AliCaloTrackReader::fRecalculateClusters | Correct clusters, recalculate them if recalibration parameters is given |
| Bool_t | AliCaloTrackReader::fRecalculateVertexBC | Recalculate vertex BC from tracks pointing to vertex |
| Bool_t | AliCaloTrackReader::fRejectEMCalTriggerEventsWith2Tresholds | Reject events EG2 also triggered by EG1 or EJ2 also triggered by EJ1 |
| TArrayI | AliCaloTrackReader::fRejectEventsWithBit | Reject events if trigger bit is on |
| Bool_t | AliCaloTrackReader::fRemoveBadTriggerEvents | Remove triggered events because trigger was exotic, bad, or out of BC |
| Bool_t | AliCaloTrackReader::fRemoveLEDEvents | Remove events where LED was wrongly firing - EMCAL LHC11a |
| Bool_t | AliCaloTrackReader::fRemoveUnMatchedTriggers | Analyze events where trigger patch and cluster where found or not |
| Bool_t | AliCaloTrackReader::fSelectEmbeddedClusters | Use only simulated clusters that come from embedding. |
| Bool_t | AliCaloTrackReader::fSelectSPDHitTracks | Ensure that track hits SPD layers |
| TString | AliCaloTrackReader::fTaskName | Name of task that executes the analysis |
| Float_t | AliCaloTrackReader::fTimeStampEventFracMax | Maximum value of time stamp fraction event |
| Float_t | AliCaloTrackReader::fTimeStampEventFracMin | Minimum value of time stamp fraction event |
| Bool_t | AliCaloTrackReader::fTimeStampEventSelect | Select events within a fraction of data taking time |
| Double_t | AliCaloTrackReader::fTimeStampRunMax | Maximum value of time stamp in run |
| Double_t | AliCaloTrackReader::fTimeStampRunMin | Minimum value of time stamp in run |
| Int_t | AliCaloTrackReader::fTrackBCEvent[19] | Fill one entry per event if there is a track in a given BC |
| Int_t | AliCaloTrackReader::fTrackBCEventCut[19] | Fill one entry per event if there is a track in a given BC, depend on track pT, acceptance cut |
| Double_t | AliCaloTrackReader::fTrackDCACut[3] | Remove tracks with DCA larger than cut, parameters of function stored here |
| Int_t | AliCaloTrackReader::fTrackMult | Track multiplicity |
| Float_t | AliCaloTrackReader::fTrackMultEtaCut | Track multiplicity eta cut |
| ULong_t | AliCaloTrackReader::fTrackStatus | Track selection bit, select tracks refitted in TPC, ITS ... |
| Double_t | AliCaloTrackReader::fTrackTimeCutMax | Remove tracks with time larger than this value, in ns |
| Double_t | AliCaloTrackReader::fTrackTimeCutMin | Remove tracks with time smaller than this value, in ns |
| Int_t | AliCaloTrackReader::fTriggerClusterBC | Event triggered by a cluster in BC -5 0 to 5 |
| Int_t | AliCaloTrackReader::fTriggerClusterId | Id of trigger cluster (cluster->GetID()) |
| Int_t | AliCaloTrackReader::fTriggerClusterIndex | Index in clusters array of trigger cluster |
| Bool_t | AliCaloTrackReader::fTriggerClusterTimeRecal | In case cluster already calibrated, do not try to recalibrate even if recalib on in RecoUtils. |
| Float_t | AliCaloTrackReader::fTriggerL0EventThreshold | L0 Threshold to look for triggered events, set outside |
| Float_t | AliCaloTrackReader::fTriggerL1EventThreshold | L1 Threshold to look for triggered events, set in data |
| Bool_t | AliCaloTrackReader::fTriggerL1EventThresholdFix | L1 Threshold is fix and set outside |
| Bool_t | AliCaloTrackReader::fTriggerPatchClusterMatch | Search for the trigger patch and check if associated cluster was the trigger |
| Int_t | AliCaloTrackReader::fTriggerPatchTimeWindow[2] | Trigger patch selection window |
| Bool_t | AliCaloTrackReader::fUseEMCALTimeCut | Do time cut selection |
| Bool_t | AliCaloTrackReader::fUseEventsWithPrimaryVertex | Select events with primary vertex |
| Bool_t | AliCaloTrackReader::fUseParamTimeCut | Use simple or parametrized time cut |
| Bool_t | AliCaloTrackReader::fUseTrackDCACut | Do DCA selection |
| Bool_t | AliCaloTrackReader::fUseTrackTimeCut | Do time cut selection |
| Int_t | AliCaloTrackReader::fV0ADC[2] | Integrated V0 signal |
| Int_t | AliCaloTrackReader::fV0Mul[2] | Integrated V0 Multiplicity |
| Double_t** | AliCaloTrackReader::fVertex | ! vertex array 3 dim for each mixed event buffer |
| Int_t | AliCaloTrackReader::fVertexBC | Vertex BC |
| Bool_t | AliCaloTrackReader::fWriteOutputDeltaAOD | Write the created delta AOD objects into file |
| Float_t | AliCaloTrackReader::fZvtxCut | Cut on vertex position |
private:
| TArrayI* | fChargedParticlesArray | Keep charged particles of this list in calorimeter. |
| Bool_t | fCheckOverlap | Check of overlapped photons from pi0 enter the calorimeter |
| Bool_t | fDecayPi0 | If not decayed, decay pi0 by hand |
| Float_t | fEMCALOverlapAngle | Aperture angle of photons from decay that is not resolved by EMCAL, in radians |
| TLorentzVector | fGamDecayMom1 | ! Gamma decay 1 momentum |
| TLorentzVector | fGamDecayMom2 | ! Gamma decay 2 momentum |
| Int_t | fIndex2ndPhoton | Check overlap of first decay photon already done, internal use. |
| Bool_t | fKeepAllStatus | Do or do not select particles depending on their status code. |
| TLorentzVector | fMomentum | ! momentum |
| TArrayI* | fNeutralParticlesArray | Do not keep neutral particles of this list in calorimeter. |
| Bool_t | fOnlyGeneratorParticles | Use particles only generated by PYTHIA/HERWIG/... and not by the MC tranport G3/G4/FLUKA ... |
| Float_t | fPHOSOverlapAngle | Aperture angle of photons from decay that is not resolved by PHOS, in radians |
| TLorentzVector | fPi0Momentum | ! Pi0 momentum |
| TArrayI* | fStatusArray | Keep particles with status of the list. |
Class Charts
| Inheritance Chart: | ||||||||||||
|
Function documentation
void InitParameters()
void CheckOverlap(Float_t anglethres, Int_t imom, Int_t& iPrimary, Int_t& index, Int_t& pdg)
Check overlap of decay photons
void FillCalorimeters(Int_t& iParticle, TParticle* particle)
Fill AODCaloClusters or TParticles lists of PHOS or EMCAL
Bool_t FillInputEvent(Int_t iEntry, const char* currentFileName)
Fill the event counter and input lists that are needed, called by the analysis maker.
Bool_t KeepParticleWithStatus(Int_t status) const
Check if status is equal to one of the list These particles will be used in analysis.
Bool_t KeepChargedParticles(Int_t pdg) const
Check if pdg is equal to one of the charged particles list These particles will be added to the calorimeters lists.
void MakePi0Decay()
, Double_t &angle) Perform isotropic decay pi0 -> 2 photons fPi0Momentum is pi0 4-momentum (ipnut) fGamDecayMom1 and fGamDecayMom2 are photon 4-momenta (output)
void SetInputOutputMCEvent(AliVEvent* esd, AliAODEvent* aod, AliMCEvent* mc)
Connect the data pointer
Bool_t SkipNeutralParticles(Int_t pdg) const
Check if pdg is equal to one of the neutral particles list These particles will be skipped from analysis.
void SetTrackChargeAndPID(Int_t pdgCode, AliAODTrack* track) const
Give a PID weight for tracks equal to 1 depending on the particle type
void SetCaloClusterPID(Int_t pdgCode, AliVCluster* calo) const
Give a PID weight for CaloClusters equal to 1 depending on the particle type
void AddNeutralParticlesArray(TArrayI& array)
Particle type, status, arrays
{ fNeutralParticlesArray = new TArrayI(array) ; }