PndHypMicroIdealWriter Class Reference

#include <PndHypMicroIdealWriter.h>

Inheritance diagram for PndHypMicroIdealWriter:

Public Types
typedef std::map< Int_t, Float_t >	mapper

Public Member Functions
	PndHypMicroIdealWriter ()
	~PndHypMicroIdealWriter ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
virtual void	Finish ()

Protected Member Functions
void	propagate (TLorentzVector &l, TVector3 &p, float charge)

Protected Attributes
TLorentzVector	track1
TLorentzVector	track2

Private Member Functions
virtual void	SetParContainers ()
	ClassDef (PndHypMicroIdealWriter, 1)

Private Attributes
int	evtcnt
TClonesArray *	fTrArray
TClonesArray *	fHitArray
TClonesArray *	fMCTrack
TClonesArray *	fChargedCandidates
TClonesArray *	fNeutralCandidates
TClonesArray *	fMcCandidates
TClonesArray *	fMicroIdealCandidates
TClonesArray *	fEventInfo
bool	fStoreNeutral
bool	fStoreCharged
bool	fStoreMC
Int_t	mvd_hitidx [1000]
Int_t	stt_hitidx [1000]
Int_t	tpc_hitidx [1000]

Detailed Description

Definition at line 28 of file PndHypMicroIdealWriter.h.

Member Typedef Documentation

typedef std::map<Int_t, Float_t> PndHypMicroIdealWriter::mapper

Definition at line 32 of file PndHypMicroIdealWriter.h.

Constructor & Destructor Documentation

PndHypMicroIdealWriter::PndHypMicroIdealWriter

(

)

Default constructor

Definition at line 54 of file PndHypMicroIdealWriter.cxx.

                                               :
  FairTask("FastSim Dump") { 
}

PndHypMicroIdealWriter::~PndHypMicroIdealWriter

(

)

Destructor

Definition at line 60 of file PndHypMicroIdealWriter.cxx.

References fChargedCandidates, fEventInfo, fMcCandidates, fMicroIdealCandidates, and fNeutralCandidates.

                                                {
  if (fChargedCandidates) {fChargedCandidates->Delete(); delete fChargedCandidates;}
  if (fNeutralCandidates) {fNeutralCandidates->Delete(); delete fNeutralCandidates;}
  if (fMcCandidates) {fMcCandidates->Delete(); delete fMcCandidates;}
  if (fMicroIdealCandidates) {fMicroIdealCandidates->Delete(); delete fMicroIdealCandidates;}
  if (fEventInfo) {fEventInfo->Delete(); delete fEventInfo;}

 }

Member Function Documentation

PndHypMicroIdealWriter::ClassDef ( PndHypMicroIdealWriter  , 1    )

private

void PndHypMicroIdealWriter::Exec ( Option_t *  opt )

virtual

Virtual method Exec

****** Extrapolation to primary vertex ******

Definition at line 172 of file PndHypMicroIdealWriter.cxx.

References d, evtcnt, fHitArray, Finish(), fMicroIdealCandidates, fStoreCharged, fTrArray, GFTrack::getCand(), GFTrack::getCardinalRep(), GFTrackCand::getHit(), GFTrackCand::getNHits(), GFAbsTrackRep::getStatusFlag(), GFTrack::getTrackRep(), hit, i, init, jj, mvd_hitidx, p, pos, PndHypPoint::PositionIn(), PndPidCandidate::SetPosition(), and tpc_hitidx.

{
  if ((++evtcnt)%100==0)cout <<"evt: "<<evtcnt<<endl;
  
  // find # in input array
  Int_t nTracks = 0;
  if (fStoreCharged) nTracks=fTrArray->GetEntriesFast();
  
  Int_t nCluster = 0;
  //if (fStoreNeutral) nCluster=fEmcArray->GetEntriesFast();
  
  int nMCTrack=0;
  //if (fStoreMC) nMCTrack=fMCTrack->GetEntriesFast();
  
  //PndStack *fStack=(PndStack*)gMC->GetStack();  
  //int nMCTracks=fStack->GetNtrack();

   // Reset output array
  //if (fChargedCandidates->GetEntriesFast() != 0)  fChargedCandidates->Clear("C");
  //  if (fNeutralCandidates->GetEntriesFast() != 0)  fNeutralCandidates->Clear("C");
  //   if (fMcCandidates->GetEntriesFast() != 0)  fMcCandidates->Clear("C");
  if (fMicroIdealCandidates->GetEntriesFast() != 0) fMicroIdealCandidates->Clear("C");
  // if (fEventInfo->GetEntriesFast() != 0) fEventInfo->Clear("C");
 
  
  //TClonesArray &chrgCandidates  = *fChargedCandidates;
 //  TClonesArray &neutCandidates  = *fNeutralCandidates;
//   TClonesArray &mctracks         = *fMcCandidates;
  TClonesArray &microCandidates = *fMicroIdealCandidates;
 
  //RhoCandList l;
  
  TLorentzVector McSumP4(0,0,0,0);
  TVector3 McAvgVtx(0,0,0);
  
  int nPrimary=0;
  
 
  GFTrack *tr1;
 
  
  // *************************
  // Loop over the charged tracks
  // ************************
   LSLTrackRep* grep=0;
   GFAbsTrackRep* rep=0;
  
  for (Int_t i=0; i<nTracks; i++)
  {
    
    //Int_t chcandsize = chrgCandidates.GetEntriesFast();
    Int_t micsize = microCandidates.GetEntriesFast();
    
    // Get the Tpc Track information
    
    tr1 = (GFTrack *)fTrArray->At(i);  
  
     // check that track has been fitted
    if(tr1->getCardinalRep()->getStatusFlag()!=0){
      std::cout<< "Discarding track. Status flag !=0" <<std::endl;
      continue;
    }
    LSLTrackRep* myrep=dynamic_cast<LSLTrackRep*>(tr1->getTrackRep(0));//tr1->getCardinalRep());
    TVectorD d(6);
    //TVector3 vtx;
    
    d = myrep->getGlobal();
    //TVector3 pos =tr1->getPos();
    //TVector3 pos2 =tr1->getTrackRep(0)->getPos();
    //TVector3 mom =tr1->getTrackRep(0)->getMom();
    
   //  std::cout<< "pos1 "<<pos.x()<<" "<<pos.x()<<" "<<pos.z()<<std::endl;
//     std::cout<< "pos2 "<<pos2.x()<<" "<<pos2.x()<<" "<<pos2.z()<<std::endl;
//     std::cout<< "mom "<<mom.x()<<" "<<mom.x()<<" "<<mom.z()<<std::endl;

    TVector3 vtx(d[0],d[1],d[2]);
   
    TLorentzVector lv;
    lv.SetXYZM(d[3],d[4],d[5],0.13957);
   
    //std::cout<< "pos glob "<<d[0]<<" "<<d[1]<<" "<<d[2]<<std::endl;
//     std::cout<< "mom glob "<<d[3]<<" "<<d[4]<<" "<<d[5]<<std::endl;
    
  
    
    // propagate(lv,vtx,myrep->getCharge());
    
    // create the RhoCandidate (keep this for the time being)
    //RhoCandidate *tcand=new (chrgCandidates[chcandsize]) RhoCandidate(lv,myrep->getCharge());
    
    //tcand->SetPos(vtx);
    
    //set the convariance matrix of tcand
    
    TMatrixD globalCov = myrep->getGlobalCov();
    TMatrixD mat(7,7);
    int ii,jj;
   
    
      for (ii=0;ii<6;ii++) for(jj=0;jj<6;jj++) mat[ii][jj]=globalCov[ii][jj]; 
    
    
    //Extend matrix for energy (with default pion hypothesis)
    double invE = 1./lv.E();
   //  mat[0+3][3+3] = mat[3+3][0+3] = (lv.X()*mat[0+3][0+3]+lv.Y()*mat[0+3][1+3]+lv.Z()*mat[0+3][2+3])*invE;
//     mat[1+3][3+3] = mat[3+3][1+3] = (lv.X()*mat[0+3][1+3]+lv.Y()*mat[1+3][1+3]+lv.Z()*mat[1+3][2+3])*invE;
//     mat[2+3][3+3] = mat[3+3][2+3] = (lv.X()*mat[0+3][2+3]+lv.Y()*mat[1+3][2+3]+lv.Z()*mat[2+3][2+3])*invE;
//     mat[3+3][3+3] = (lv.X()*lv.X()*mat[0+3][0+3]+lv.Y()*lv.Y()*mat[1+3][1+3]+lv.Z()*lv.Z()*mat[2+3][2+3]
//                   +2.0*lv.X()*lv.Y()*mat[0+3][1+3]
//                   +2.0*lv.X()*lv.Z()*mat[0+3][2+3]
//                   +2.0*lv.Y()*lv.Z()*mat[1+3][2+3])*invE*invE;
    
//     mat[3+3][4-4] = mat[4-4][3+3] = (lv.X()*mat[0+3][4-4]+lv.Y()*mat[1+3][4-4]+lv.Z()*mat[2+3][4-4])*invE;
//     mat[3+3][5-4] = mat[5-4][3+3] = (lv.X()*mat[0+3][5-4]+lv.Y()*mat[1+3][5-4]+lv.Z()*mat[2+3][5-4])*invE;
//     mat[3+3][6-4] = mat[6-4][3+3] = (lv.X()*mat[0+3][6-4]+lv.Y()*mat[1+3][6-4]+lv.Z()*mat[2+3][6-4])*invE;
    
    //tcand->SetCov7(mat);
    
    //l.Add(*tcand);
     unsigned int detId, hitId;
    unsigned int numhits=0,mvdhits=0,stthits=0,tpchits=0;
    
    numhits=tr1->getCand().getNHits();
    //cout<<" numhits "<<numhits<<endl;
  // create the PndMicroIdealCandidate
    
    PndPidCandidate *micro=new (microCandidates[micsize]) PndPidCandidate((Int_t)myrep->getCharge(),vtx,lv,mat);//myrep->getCharge()
    
    
    
    for (ii=0;ii<numhits;++ii)
      { 
        tr1->getCand().getHit(ii,detId,hitId);
        //cout<<" numhits loop "<<ii<<" "<<detId<<" "<<hitId<<endl;
        switch (detId)
          {
          case 2:  mvd_hitidx[mvdhits]=hitId; if(mvdhits<1000) mvdhits++; break;
            //case 1:  stt_hitidx[stthits]=hitId; if(stthits<1000) stthits++;break;
          default: tpc_hitidx[tpchits]=hitId; if(tpchits<1000) tpchits++;break;
          }
      }
    micro->SetMvdHitIndexArray(mvdhits, mvd_hitidx);
    //micro->SetSttHitIndexArray(stthits, stt_hitidx);
    //micro->SetTpcHitIndexArray(tpchits, tpc_hitidx);
    

    
    //rep=tr1->getTrackRep(0)->clone();
 
     tr1->getCand().getHit(0,detId,hitId);
     PndHypPoint* hit = (PndHypPoint*) fHitArray->At(hitId);
   
   
    //std::cout<<" hit "<<hit->GetX()<<std::endl;
     TVector3 init;
     if(hit==0)continue;
      hit->PositionIn(init);
    
    
     //DetPlane pinit(init,TVector3(1,0,0),TVector3(0,1,0));
   
    //myrep->extrapolate(pinit);
    // myrep->setReferencePlane(pinit);
   
    TVector3 pos(0,0,-76.5);
    GFDetPlane pfin(pos,TVector3(1,0,0),TVector3(0,1,0));

    TVector3 dist,fin;
    
     TMatrixT<double> state(5,1);
     TMatrixT<double> cov(5,5);
     GFDetPlane p;
     dist=myrep->getPos(pfin) ;
     micro->SetPosition(dist);

     //std::cout<< (micro->GetPosition()).x()<<std::endl;
     //dist = myrep->extrapolateToPoca(pos,state,cov,p);
     //vtx.SetXYZ(dist.x(),dist.y(),dist.z());

    // std::cout<<dist.x()<<" "<<dist.y()<<" "<<dist.z()<<std::endl;
    //std::cout<<fin.x()<<" "<<fin.y()<<" "<<fin.z()<<std::endl;
    // grep=dynamic_cast<LSLTrackRep*>(rep);

//     if(grep!=0){

//       //grep->setPropDir(-1);
//       TMatrixT<double> state(5,1);
//       TMatrixT<double> cov(5,5);
//       DetPlane p;
//       //dist = grep->extrapolateToPoca(pos,state,cov,p);
//       //dist=trk->getPos();
//       //grep->Print();
      
//       vtx.SetXYZ(dist.x(),dist.y(),dist.z());
      
//       //double length=grep->extrapolate(pmed,state);
//       //std::cout<<"length "<<length<<std::endl;
//       //std::cout<<" x poca "<<grep->getChiSqu()<<" "<<dist.x()<<" "<<dist.y()<<" "<<dist.z()<<std::endl;
//       //std::cout<<"error poca "<<sqrt(grep->getCovElem(0, 0))<<" erry "<<sqrt(grep->getCovElem(1, 1))<<std::endl;
      
//     }

  
    
  }
  
  // *************************
  // Loop over the neutral clusters
  // ************************
  
  double calFactor=1.035;
  
  
  
    
    Finish();
    
}

void PndHypMicroIdealWriter::Finish ( )

virtual

Definition at line 393 of file PndHypMicroIdealWriter.cxx.

References fHitArray, and fTrArray.

Referenced by Exec().

 {
   //fChargedCandidates->Delete();
//   //  fNeutralCandidates->Delete(); 
// //    fMcCandidates->Delete(); 
   fTrArray->Delete();
   fHitArray->Delete();
   //fMicroIdealCandidates->Delete();
//    //fEventInfo->Delete(); 
//  //fInvMass->Write();
 }

InitStatus PndHypMicroIdealWriter::Init ( )

virtual

Virtual method Init

Definition at line 73 of file PndHypMicroIdealWriter.cxx.

References evtcnt, fHitArray, fMicroIdealCandidates, fStoreCharged, fStoreMC, fStoreNeutral, and fTrArray.

{
 
  fStoreNeutral=true;
  fStoreCharged=true;
  fStoreMC=true;
 
  //TDatabasePDG *dbpdg=TDatabasePDG::Instance();
  
  //TRho::Instance()->SetPDG(dbpdg);
  
  cout << " Inside the Init function****" << endl;
  
  //FairDetector::Initialize();
  //FairRun* sim = FairRun::Instance();
  //FairRuntimeDb* rtdb=sim->GetRuntimeDb();
  
  // Get RootManager
  FairRootManager* ioman = FairRootManager::Instance();
  if ( ! ioman ) {
    cout << "-E- PndHypMicroIdealWriter::Init: "
         << "RootManager not instantiated!" << endl;
    return kFATAL;
  }
  // Get input array
  fTrArray = (TClonesArray*) ioman->GetObject("Track");
  if ( ! fTrArray) {
    cout << "-W- PndHypMicroIdealWriter::Init: "
         << "No TpcTrack array!" << endl;
    fTrArray=new TClonesArray("Track");
    fStoreCharged=false;
    //return kERROR;
  }
  
  fHitArray = (TClonesArray*) ioman->GetObject("HypPoint");
  if ( ! fHitArray) {
    cout << "-W- PndHypMicroIdealWriter::Init: "
         << "No Hit array!" << endl;
    fHitArray = new TClonesArray("HypPoint");
    fStoreNeutral=false;
    //return kERROR;
  }
  
  /*
  fMCTrack = (TClonesArray*) ioman->GetObject("MCTrack");
  if ( ! fMCTrack) {
    cout << "-W- PndHypMicroIdealWriter::Init: "
         << "No MCTrack array!" << endl;
         fMCTrack = new TClonesArray("MCTrack");
    fStoreMC=false;
    //return kERROR;
  }
*/
  
  // fChargedCandidates = new TClonesArray("RhoCandidate");
//   FairRootManager::Instance()->Register("PndChargedCandidates","FullSim", fChargedCandidates, kTRUE);
  
//   fNeutralCandidates = new TClonesArray("RhoCandidate");
//   FairRootManager::Instance()->Register("PndNeutralCandidates","FullSim", fNeutralCandidates, kTRUE);
  
//   fMcCandidates = new TClonesArray("RhoCandidate");
//   FairRootManager::Instance()->Register("PndMcTracks","FullSim", fMcCandidates, kTRUE);

  fMicroIdealCandidates = new TClonesArray("PndPidCandidate");
  FairRootManager::Instance()->Register("PndMicroIdealCandidates","FullSim", fMicroIdealCandidates, kTRUE);

//   fEventInfo = new TClonesArray("PndEventInfo");
//   FairRootManager::Instance()->Register("PndEventSummary","FullSim", fEventInfo, kTRUE);

  // Create and register output array
  cout << "-I- PndHypMicroIdealWriter: Intialization successfull" << endl;

  //fInvMass = new TH1D("invmass","",100,0.0,4.0);
  
  evtcnt=0;
  
  return kSUCCESS;

}

void PndHypMicroIdealWriter::propagate ( TLorentzVector &  l, TVector3 &  p, float  charge  )

protected

void PndHypMicroIdealWriter::SetParContainers ( )

privatevirtual

Geo file to use Get parameter containers

Definition at line 153 of file PndHypMicroIdealWriter.cxx.

References rtdb.

                                              {

 // Get Base Container
  FairRunAna* ana = FairRunAna::Instance();
  FairRuntimeDb* rtdb=ana->GetRuntimeDb();

  // Get run and runtime database
 //  FairRunAna* run = FairRunAna::Instance();
//   if ( ! run ) Fatal("SetParContainers", "No analysis run");

 //  FairRuntimeDb* db = run->GetRuntimeDb();
//   if ( ! db ) Fatal("SetParContainers", "No runtime database");

 
}

Member Data Documentation

int PndHypMicroIdealWriter::evtcnt

private

book all the histograms

Definition at line 62 of file PndHypMicroIdealWriter.h.

Referenced by Exec(), and Init().

TClonesArray* PndHypMicroIdealWriter::fChargedCandidates

private

Definition at line 74 of file PndHypMicroIdealWriter.h.

Referenced by ~PndHypMicroIdealWriter().

TClonesArray* PndHypMicroIdealWriter::fEventInfo

private

Definition at line 80 of file PndHypMicroIdealWriter.h.

Referenced by ~PndHypMicroIdealWriter().

TClonesArray* PndHypMicroIdealWriter::fHitArray

private

Input array of EmcCandidates

Definition at line 68 of file PndHypMicroIdealWriter.h.

Referenced by Exec(), Finish(), and Init().

TClonesArray* PndHypMicroIdealWriter::fMcCandidates

private

Definition at line 76 of file PndHypMicroIdealWriter.h.

Referenced by ~PndHypMicroIdealWriter().

TClonesArray* PndHypMicroIdealWriter::fMCTrack

private

Input array of MCTracks

Definition at line 71 of file PndHypMicroIdealWriter.h.

TClonesArray* PndHypMicroIdealWriter::fMicroIdealCandidates

private

Definition at line 77 of file PndHypMicroIdealWriter.h.

Referenced by Exec(), Init(), and ~PndHypMicroIdealWriter().

TClonesArray* PndHypMicroIdealWriter::fNeutralCandidates

private

Definition at line 75 of file PndHypMicroIdealWriter.h.

Referenced by ~PndHypMicroIdealWriter().

bool PndHypMicroIdealWriter::fStoreCharged

private

Definition at line 83 of file PndHypMicroIdealWriter.h.

Referenced by Exec(), and Init().

bool PndHypMicroIdealWriter::fStoreMC

private

Definition at line 84 of file PndHypMicroIdealWriter.h.

Referenced by Init().

bool PndHypMicroIdealWriter::fStoreNeutral

private

Definition at line 82 of file PndHypMicroIdealWriter.h.

Referenced by Init().

TClonesArray* PndHypMicroIdealWriter::fTrArray

private

Input array of TpcLheTrack

Definition at line 65 of file PndHypMicroIdealWriter.h.

Referenced by Exec(), Finish(), and Init().

Int_t PndHypMicroIdealWriter::mvd_hitidx[1000]

private

Definition at line 86 of file PndHypMicroIdealWriter.h.

Referenced by Exec().

Int_t PndHypMicroIdealWriter::stt_hitidx[1000]

private

Definition at line 87 of file PndHypMicroIdealWriter.h.

Int_t PndHypMicroIdealWriter::tpc_hitidx[1000]

private

Definition at line 88 of file PndHypMicroIdealWriter.h.

Referenced by Exec().

TLorentzVector PndHypMicroIdealWriter::track1

protected

Definition at line 56 of file PndHypMicroIdealWriter.h.

TLorentzVector PndHypMicroIdealWriter::track2

protected

Definition at line 57 of file PndHypMicroIdealWriter.h.

---

The documentation for this class was generated from the following files:

• PndHypMicroIdealWriter.h
• PndHypMicroIdealWriter.cxx