#include <PndDskFLGHitProducerIdeal.h>

Inheritance diagram for PndDskFLGHitProducerIdeal:

Public Member Functions
	PndDskFLGHitProducerIdeal ()

	PndDskFLGHitProducerIdeal (Int_t verbose)

virtual	~PndDskFLGHitProducerIdeal ()

virtual InitStatus	Init ()

virtual void	Exec (Option_t *option)

virtual void	Finish ()

PndDskFLGHit *	AddHit (Int_t trackID, Int_t detectorID, TVector3 position_store, TVector3 momentum_store, Double_t time, Double_t angIn, Double_t thetaC_store, TVector3 Cherenkov_photon, Int_t light_guide, Int_t pixel)

void	SetPersistency (Bool_t val=kTRUE)

Bool_t	GetPersistency ()

Private Member Functions
void	SetDefaultParameters ()

Private Attributes
Int_t	fDetectorID

TVector3	fPosHit

TVector3	fDPosHit

Double_t	fThetaC

Double_t	fErrThetaC

Int_t	fRefIndex

Int_t	f_light_guide

Int_t	f_pixel

TClonesArray *	fDskParticleArray

TClonesArray *	fHitArray

PndGeoDskFLG *	fGeo

Int_t	fVerbose

Int_t	nevents

TRandom *	rand

Detailed Description

PndDskFLGHitProducer.h

Class for producing DSK hits directly from MCPoints

Definition at line 39 of file PndDskFLGHitProducerIdeal.h.

Constructor & Destructor Documentation

PndDskFLGHitProducerIdeal::PndDskFLGHitProducerIdeal ( )

Default constructor

Definition at line 41 of file PndDskFLGHitProducerIdeal.cxx.

References fGeo.

 :PndPersistencyTask("PndDskFLGHitProducerIdeal")
 {
 
   fGeo         = new PndGeoDskFLG();  
 }

PndDskFLGHitProducerIdeal::PndDskFLGHitProducerIdeal ( Int_t verbose )

Constructor with verbosity

Definition at line 51 of file PndDskFLGHitProducerIdeal.cxx.

References fGeo, fVerbose, and verbose.

   :PndPersistencyTask("PndDskFLGHitProducerIdeal")
 {
   fGeo         = new PndGeoDskFLG();
   fVerbose = verbose;  
 }

PndDskFLGHitProducerIdeal::~PndDskFLGHitProducerIdeal ( )

virtual

Destructor

Definition at line 61 of file PndDskFLGHitProducerIdeal.cxx.

References fGeo.

 {
  
   if (fGeo) delete fGeo; 
 }

Member Function Documentation

PndDskFLGHit * PndDskFLGHitProducerIdeal::AddHit	(	Int_t	trackID,
		Int_t	detectorID,
		TVector3	position_store,
		TVector3	momentum_store,
		Double_t	time,
		Double_t	angIn,
		Double_t	thetaC_store,
		TVector3	Cherenkov_photon,
		Int_t	light_guide,
		Int_t	pixel
	)

method AddHit

Adds a DskHit to the HitCollection

Definition at line 203 of file PndDskFLGHitProducerIdeal.cxx.

References fHitArray.

Referenced by Exec().

                                      {
   TClonesArray& clref = *fHitArray;
   Int_t size = clref.GetEntriesFast();
   return new(clref[size]) PndDskFLGHit( trackID,detectorID,position_store,momentum_store,time,
                 angIn,thetaC_store,Cherenkov_photon,light_guide,pixel);
 }

void PndDskFLGHitProducerIdeal::Exec ( Option_t * option )

virtual

Executed task

Definition at line 108 of file PndDskFLGHitProducerIdeal.cxx.

References acos(), AddHit(), Double_t, energy, f_light_guide, f_pixel, fabs(), fDskParticleArray, fGeo, fHitArray, fVerbose, PndDskParticle::GetAngIn(), PndDskParticle::GetMass(), PndDskParticle::GetNPhot(), PndDskParticle::GetThetaC(), nevents, P, Pi, printf(), PndGeoDskFLG::Propagate(), pt(), rand, PndGeoDskFLG::reflect_threshold(), sqrt(), CAMath::Sqrt(), and thetaC.

 {
   if ( ! fHitArray ) Fatal("Exec", "No HitArray");
   fHitArray->Clear();
 
 
   PndDskParticle* pt=NULL;
   nevents++;
    
   if (fVerbose > 1) {
     cout << " ----------------- DSK Hit Producer --------------------" << endl;
     cout <<" Number of input MC points in the plate: "<<fDskParticleArray->GetEntries()<<endl;
   }
 
   // fNHits = 0;
 
   // Loop over PndDskPoints
   for(Int_t j=0; j<fDskParticleArray->GetEntriesFast(); j++) {
     if (fVerbose > 1) printf("\n\n=====> Event No. %d\n", nevents); 
     
     pt = (PndDskParticle*)fDskParticleArray->At(j);
   
     Double_t Px= pt->GetPx();
     Double_t Py= pt->GetPy();
     Double_t Pz= pt->GetPz();
     Double_t P = sqrt(Px*Px + Py*Py +Pz*Pz);
     Double_t mass = pt->GetMass();
     Double_t energy = TMath::Sqrt(P*P + mass*mass); 
 
     Double_t beta;
     if(energy != 0) {
       beta = P/energy;
     }
     else {
       beta = -1.;
       if (fVerbose >1) cout << "Beta not calculated " << endl; 
     }
 
     Double_t thetaC;
     if (fabs(1./(1.47*beta)) > 1. || P == 0. || energy == 0.){
       thetaC = -1.;
     }else{
       thetaC = acos(1/(1.47*beta));
     }
 
     if (thetaC != -1. && beta > 1/1.47){  
 
     //collect that need to be stored
     Int_t trackID = pt->GetTrackID(); 
     Int_t detectorID = 0;
     TVector3 position_store(pt->GetX(), pt->GetY(), pt->GetZ());
     TVector3 momentum_store(Px, Py, Pz);
     Double_t time = 0;
     Double_t angIn = pt->GetAngIn();
     Double_t thetaC_store = pt->GetThetaC();  
     
     //generate Cherenkov photons in a cone
 
     int N_cherenko_photon = pt->GetNPhot();
 
     TVector3 Cone_Pos(pt->GetX()*10, pt->GetY()*10, pt->GetZ()*10);
     TVector3 Cone_Axis(Px, Py, Pz);
     TVector3 Cherenkov_photon(1,0,0); //initialize
     Cherenkov_photon.SetTheta(Cone_Axis.Theta() + thetaC);
     Cherenkov_photon.SetPhi(Cone_Axis.Phi());
     for(int i_photon = 0; i_photon < N_cherenko_photon; i_photon++){
        Double_t rot_angle = 2*(TMath::Pi())*rand->Rndm();
        Cherenkov_photon.Rotate(rot_angle, Cone_Axis); 
 
        //cout<<"Cone_Axis: "<<Cone_Axis.X()<<" "<<Cone_Axis.Y()<<" "<<Cone_Axis.Z()<<endl;
        //cout<<"Cone_Pos: "<<Cone_Pos.X()<<" "<<Cone_Pos.Y()<<" "<<Cone_Pos.Z()<<endl;
        //cout<<"photon_dir: "<<Cherenkov_photon.X()<<" "<<Cherenkov_photon.Y()<<" "<<Cherenkov_photon.Z()<<"  theta: "<<Cherenkov_photon.Theta()<<"  phi: "<<Cherenkov_photon.Phi()<<endl;
        if(Cherenkov_photon.Theta()<0.8 || Cherenkov_photon.Z()<0) continue;
 
        f_light_guide = -99;
        f_pixel = -99; //fixme
        //Propagate(Cone_Pos, Cherenkov_photon);  
 
        if(Cherenkov_photon.Theta()>fGeo->reflect_threshold())  fGeo->Propagate(Cone_Pos, Cherenkov_photon, f_light_guide, f_pixel);
 
        if(f_light_guide != -99)
         AddHit(trackID,detectorID,position_store,momentum_store,time,
                 angIn,thetaC_store,Cherenkov_photon,f_light_guide,f_pixel);
 
     }
 
 
     }
   }
 }

void PndDskFLGHitProducerIdeal::Finish ( )

virtual

Finish task

Definition at line 217 of file PndDskFLGHitProducerIdeal.cxx.

 {
 
   cout << "-I- PndDskFLGHitProducerIdeal: Finish" << endl;
  }

Bool_t PndPersistencyTask::GetPersistency ( )

inlineinherited

Definition at line 32 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency.

Referenced by PndLmdPixelHitProducerFast::GetPersistance(), PndMdtDigitization::Init(), PndMdtHitProducerIdeal::Init(), PndMdtClusterTask::Init(), PndFtsHitProducerRealFast::Init(), PndRichHitProducer::Init(), PndSttHitProducerRealFast::Init(), PndDiscTaskReconstruction::Init(), PndSttHelixHitProducer::Init(), PndDiscTaskPID::Init(), PndIdealTrackFinder::Init(), PndSttMvdGemTracking::Init(), PndMdtTrkProducer::Init(), PndFtsHitProducerRealFull::Init(), PndLmdPixelClusterTask::Init(), PndSttHitProducerRealFull::Init(), PndLmdStripClusterTask::Init(), PndEmcApdHitProducer::Init(), PndMissingPzCleanerTask::Init(), PndEmcMakeRecoHit::Init(), PndEmcMakeClusterOnline::Init(), PndTrackSmearTask::Init(), PndEmcFWEndcapTimebasedWaveforms::Init(), PndSttHitProducerIdeal::Init(), PndEmcFWEndcapDigi::Init(), PndFtsHitProducerIdeal::Init(), PndEmcMakeCluster::Init(), PndMdtPointsToWaveform::Init(), PndDiscTaskDigitization::Init(), PndEmcMakeDigi::Init(), PndSdsTimeWalkCorrTask::Init(), PndLmdPixelHitProducerFast::Init(), PndDrcHitFinder::Init(), PndRichHitFinder::Init(), PndEmcMakeCorr::Init(), PndFtofHitProducerIdeal::Init(), PndEmcHitsToWaveform::Init(), PndSciTDigiTask::Init(), PndDrcHitProducerIdeal::Init(), PndSdsHitProducerIdeal::Init(), PndSciTHitProducerIdeal::Init(), PndRecoMultiKalmanTask2::Init(), PndEmcHitProducer::Init(), PndDrcHitProducerReal::Init(), Init(), PndEmcTmpWaveformToDigi::Init(), PndDrcDigiTask::Init(), PndEmcWaveformToDigi::Init(), PndSttMatchTracks::Init(), PndEmcWaveformToCalibratedDigi::Init(), PndTrkTracking2::Init(), PndSttFindTracks::Init(), PndEmcMultiWaveformToCalibratedDigi::Init(), PndRecoKalmanTask2::Init(), PndDrcTimeDigiTask::Init(), PndEmcExpClusterSplitter::Init(), PndFtsHoughTrackerTask::Init(), PndSdsNoiseProducer::Init(), PndEmcPhiBumpSplitter::Init(), PndSdsIdealRecoTask::Init(), PndSdsHybridHitProducer::Init(), PndRecoMultiKalmanTask::Init(), PndSdsIdealClusterTask::Init(), PndRecoKalmanTask::Init(), PndSdsStripHitProducerDif::Init(), PndGemDigitize::Init(), PndSdsStripHitProducer::Init(), PndGemFindHits::Init(), PndSdsPixelClusterTask::Init(), PndSdsStripClusterTask::Init(), PndMvdGemTrackFinderOnHits::Init(), PndBarrelTrackFinder::Init(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcMakeBump::PndEmcMakeBump(), PndUnassignedHitsTask::RegisterBranches(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndEmcMakeBump::SetStorageOfData(), and PndEmcFullDigiTask::StoreDigi().

32 { return fPersistency; }

PndPersistencyTask::fPersistency

Bool_t fPersistency

Definition: PndPersistencyTask.h:35

InitStatus PndDskFLGHitProducerIdeal::Init ( )

virtual

Initialization of the task

Definition at line 71 of file PndDskFLGHitProducerIdeal.cxx.

References fDskParticleArray, fHitArray, PndPersistencyTask::GetPersistency(), nevents, and rand.

 {
   cout << " ---------- INITIALIZATION ------------" << endl;
   nevents = 0;
   // Get RootManager
   FairRootManager* ioman = FairRootManager::Instance();
   if ( ! ioman ) {
     cout << "-E- PndDskFLGHitProducerIdeal::Init: "
          << "RootManager not instantiated!" << endl;
     return kFATAL;
   }
 
   // Get input array
   fDskParticleArray = (TClonesArray*) ioman->GetObject("DskParticle");
    if ( ! fDskParticleArray ) {
     cout << "-W- PndDskFLGHitProducerIdeal::Init: "
          << "No DskBarPoint array!" << endl;
     return kERROR;
   }
 
    // Create and register output array
    fHitArray = new TClonesArray("PndDskFLGHit");
    ioman->Register("PndDskFLGHit","Dsk",fHitArray, GetPersistency());
 
    //initialize random for producing cherenkov photons
    rand = new TRandom();   
  
    //fVerbose = 2;
    cout << "-I- PndDskFLGHitProducerIdeal: Intialization successfull" << endl;
 
    return kSUCCESS;
    
 }