// -------------------------------------------------------------------------
// -----                        R3BTwinMusic header file               -----
// -----                  Created 26/03/12  by J.L. Rodriguez          -----
// -----                                                               -----
// -------------------------------------------------------------------------

/**  R3BTwinMusic.h
 **/

#ifndef R3BTwinMusic_H
#define R3BTwinMusic_H

#include "R3BDetector.h"

#include "TLorentzVector.h"

class TClonesArray;
class R3BTwinMusicPoint;
class FairVolume;
class TGeoRotation;

class R3BTwinMusic : public R3BDetector
{

 public:

  /** Default constructor **/
  R3BTwinMusic();


  /** Standard constructor.
   *@param name    detetcor name
   *@param active  sensitivity flag
   **/
  R3BTwinMusic(const char* name, Bool_t active);


  /** Destructor **/
  virtual ~R3BTwinMusic();


  /** Virtual method ProcessHits
   **
   ** Defines the action to be taken when a step is inside the
   ** active volume. Creates a R3BTwinMusicPoint and adds it to the
   ** collection.
   *@param vol  Pointer to the active volume
   **/
  virtual Bool_t ProcessHits(FairVolume* vol = 0);


  /** Virtual method BeginEvent
   **
   ** If verbosity level is set, print hit collection at the
   ** end of the event and resets it afterwards.
   **/

  virtual void BeginEvent();

  /** Virtual method EndOfEvent
   **
   ** If verbosity level is set, print hit collection at the
   ** end of the event and resets it afterwards.
   **/

  virtual void EndOfEvent();



  /** Virtual method Register
   **
   ** Registers the hit collection in the ROOT manager.
   **/
  virtual void Register();


  /** Accessor to the hit collection **/
  virtual TClonesArray* GetCollection(Int_t iColl) const;


  /** Virtual method Print
   **
   ** Screen output of hit collection.
   **/
  virtual void Print() const;


  /** Virtual method Reset
   **
   ** Clears the hit collection
   **/
  virtual void Reset();


  /** Virtual method CopyClones
   **
   ** Copies the hit collection with a given track index offset
   *@param cl1     Origin
   *@param cl2     Target
   *@param offset  Index offset
   **/
  virtual void CopyClones(TClonesArray* cl1, TClonesArray* cl2,
        Int_t offset);


  /** Virtaul method Construct geometry
   **
   ** Constructs the TwinMusic geometry
   **/
  virtual void ConstructGeometry();

  virtual void Initialize();
  virtual void SetSpecialPhysicsCuts();
  void SetEnergyCutOff( Double_t cutE ){fCutE = cutE;}
  Double_t  GetEnergyCutOff ( ) {return fCutE;}


  private:

    /** Track information to be stored until the track leaves the
  active volume. **/
    Int_t          fTrackID;           //!  track index
    Int_t          fVolumeID;          //!  volume id
    Int_t          fPdgId;             //!  Pdg 
    TLorentzVector fPosIn, fPosOut;    //!  position
    TLorentzVector fMomIn, fMomOut;    //!  momentum
    Double32_t     fTime;              //!  time
    Double32_t     fLength;            //!  length
    Double32_t     fELoss;             //!  energy loss
    Double32_t     fCutE;              //!  energy loss

    Int_t          fPosIndex;          //!
    TClonesArray*  ftwinCollection;    //!  The hit collection
    Bool_t         kGeoSaved;          //!
    TList *flGeoPar;                   //!
    Int_t fGeoVersion;                 //!

    Int_t fDetectorType[19];           //!
    
    /** Private method AddHit
     **
     ** Adds a TwinMusicPoint to the HitCollection
     **/
    R3BTwinMusicPoint* AddHit(Int_t trackID, Int_t detID, Int_t type, Int_t pdg, TVector3 posIn,
      TVector3 pos_out, TVector3 momIn, 
      TVector3 momOut, Double_t time, 
      Double_t length, Double_t eLoss);


    /** Private method ResetParameters
     **
     ** Resets the private members for the track parameters
     **/
    void ResetParameters();

    Int_t  GetDetectorType(Int_t volID);



    ClassDef(R3BTwinMusic,1);

};

inline Int_t R3BTwinMusic::GetDetectorType(Int_t volID) {
Int_t type=-1;

for (Int_t i=0;i<20;i++ ){
    if (volID==fDetectorType[i]) {
	type=i+1; 
	return (type);
    }
}
return type;
}


inline void R3BTwinMusic::ResetParameters() {
  fTrackID = fVolumeID = 0;
  fPosIn.SetXYZM(0.0, 0.0, 0.0, 0.0);
  fPosOut.SetXYZM(0.0, 0.0, 0.0, 0.0);
  fMomIn.SetXYZM(0.0, 0.0, 0.0, 0.0);
  fMomOut.SetXYZM(0.0, 0.0, 0.0, 0.0);
  fTime = fLength = fELoss = 0;
  fPosIndex = 0;
};


#endif