ERNeuRad.h

/********************************************************************************
 *              Copyright (C) Joint Institute for Nuclear Research              *
 *                                                                              *
 *              This software is distributed under the terms of the             * 
 *         GNU Lesser General Public Licence version 3 (LGPL) version 3,        *  
 *                  copied verbatim in the file "LICENSE"                       *
 ********************************************************************************/

#ifndef ERNeuRad_H
#define ERNeuRad_H

#include "TClonesArray.h"
#include "TLorentzVector.h"
#include "TArrayI.h"

#include "FairVolume.h"

#include "ERDetector.h"
#include "ERMCTrack.h" //for ExpertTrackingStatus

class TClonesArray;
class TLorentzVector;
class TArrayI;

class FairVolume;

class ERNeuRadStep;
class ERNeuRadPoint;

class ERNeuRad : public ERDetector
{

public:

    ERNeuRad();

    ERNeuRad(const char* name, Bool_t active, Int_t verbose);

    virtual ~ERNeuRad();

    /* Modifiers */

    void SetStorePrimarySteps() {fStorePrimarySteps = kTRUE;}

    void SetStoreAllSteps() {fStoreAllSteps = kTRUE;}

    void SetMaxStep(Double_t step) {fStep = step;}

    /*Accessors*/

    virtual TClonesArray* GetCollection(Int_t iColl) const;

public:

    virtual void Initialize();

    virtual void Register();

    virtual Bool_t ProcessHits(FairVolume* vol = 0);

    virtual void EndOfEvent();

    virtual void Print(Option_t *option="") const;

    virtual void Reset();

    virtual void CopyClones(TClonesArray* cl1, TClonesArray* cl2, Int_t offset);

    virtual Bool_t CheckIfSensitive(std::string name);

private:

    void StartNewPoint();

    void FinishNewPoint();

    ERNeuRadPoint* AddPoint();

    ERNeuRadStep* AddFirstStep();

    ERNeuRadStep* AddStep();

    Double_t CurPointLen(TLorentzVector& posIn);

private:

    TClonesArray*  fNeuRadPoints;     
    TClonesArray*  fNeuRadFirstStep;  
    TClonesArray*  fNeuRadSteps;      

    Bool_t fStorePrimarySteps;        
    Bool_t fStoreAllSteps;            
    Int_t   fVerbose;                 

    Float_t fFullEnergy;              
    Float_t fFullLY;                  

    Float_t fStep;                    

    /* current step information*/

    Int_t          fEventID;          
    Int_t          fTrackID;          
    Int_t          fMot0TrackID;      
    Double_t       fMass;             
    TLorentzVector fPosIn;            
    TLorentzVector fPosOut;           
    TLorentzVector fCurPosIn;         
    TVector3       fPosInLocal;       
    TLorentzVector fMomIn;            
    TLorentzVector fMomOut;           
    TLorentzVector fCurMomIn;         
    Double32_t     fTimeIn;           
    Double32_t     fTimeOut;          
    Double32_t     fTrackLength;      
    Double32_t     fELoss;            
    Double_t       fLightYield;       
    Int_t          fFiberNb;          
    Int_t          fPixelNb;          
    Int_t          fModuleNb;         
    Int_t          fStepNb;           
    ExpertTrackingStatus fTrackStatus; 
    TArrayI        fProcessesID;      

    ClassDef(ERNeuRad,1)
};

#endif // ERNeuRad_H