PndHypStripClusterTask Class Reference

#include <PndHypStripClusterTask.h>

Inheritance diagram for PndHypStripClusterTask:

Public Member Functions
	PndHypStripClusterTask ()
	PndHypStripClusterTask (Double_t chargecut, TString geofile)
virtual	~PndHypStripClusterTask ()
virtual void	SetParContainers ()
virtual InitStatus	Init ()
virtual InitStatus	ReInit ()
virtual void	Exec (Option_t *opt)

Private Member Functions
TVector2	CalcLineCross (TVector2 point1, TVector2 dir1, TVector2 point2, TVector2 dir2)
void	Register ()
void	Reset ()
void	ProduceHits ()
	ClassDef (PndHypStripClusterTask, 3)

Private Attributes
TClonesArray *	fDigiArray
TClonesArray *	fClusterArray
TClonesArray *	fHitArray
Double_t	fRadius
Int_t	fFEcolumns
Int_t	fFErows
TString	fGeoFile
Double_t	fChargeCut
PndHypStripDigiPar *	fDigiPar
PndHypCalcStrip *	fstripcalcTOP
	Digitization Parameters. More...
PndHypCalcStrip *	fstripcalcBOT
PndHypGeoHandling *	fGeoH
PndGeoHypPar *	fGeoPar
	Geometry name handling. More...

Detailed Description

Definition at line 32 of file PndHypStripClusterTask.h.

Constructor & Destructor Documentation

PndHypStripClusterTask::PndHypStripClusterTask

(

)

Default constructor

Definition at line 37 of file PndHypStripClusterTask.cxx.

References fChargeCut, and fGeoFile.

                                               :
  FairTask("HYP Strip Clustertisation Task")
{
  fChargeCut = 1.e6; // this ist really large and shall have no effect
  fGeoFile = "";
}

PndHypStripClusterTask::PndHypStripClusterTask	(	Double_t	chargecut,
		TString	geofile
	)

Definition at line 46 of file PndHypStripClusterTask.cxx.

References chargecut, fChargeCut, and fGeoFile.

                                                                                 :
  FairTask("HYP Strip Clustertisation Task")
{
  fChargeCut = chargecut;
  fGeoFile = geoFile;
}

PndHypStripClusterTask::~PndHypStripClusterTask ( )

virtual

Destructor

Definition at line 56 of file PndHypStripClusterTask.cxx.

References fGeoH.

{
  if(0!=fGeoH)delete fGeoH;
}

Member Function Documentation

TVector2 PndHypStripClusterTask::CalcLineCross ( TVector2  point1, TVector2  dir1, TVector2  point2, TVector2  dir2  )

private

Definition at line 425 of file PndHypStripClusterTask.cxx.

References Double_t, dx, dy, s, x, and y.

Referenced by Exec().

{
  Double_t dx, dy, s,  M, x, y;//t, //[R.K. 01/2017] unused variable
  dx = point2.X() - point1.X();
  dy = point2.Y() - point1.Y();
  
  M  = dir1.X()*dir2.Y() - dir1.Y()*dir2.X();
  
  if(M !=0.)
    {
      s  = dir1.Y()*dx/M - dir1.X()*dy/M;
      x  = point2.X() + dir2.X()*s;
      y  = point2.Y() + dir2.Y()*s;
    }
  else
    {
      std::cout<<"Warning in PndHypStripClusterTask::CalcLineCross(): M=0 setting (x,y) to 0"<<std::endl;
      x=0.;y=0.;
    }
  
  TVector2 result(x,y);
  return result;
}

PndHypStripClusterTask::ClassDef ( PndHypStripClusterTask  , 3    )

private

void PndHypStripClusterTask::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Muss ich den copy construktor selbst schreiben? fehlt was im linkdef?

Definition at line 144 of file PndHypStripClusterTask.cxx.

References PndHypStripClusterBuilder::AddDigi(), b, PndHypCalcStrip::CalcFeChToStrip(), CalcLineCross(), PndHypCalcStrip::CalcStripPointOnLine(), cos(), Double_t, fabs(), fChargeCut, fClusterArray, fDigiArray, fDigiPar, fGeoH, fHitArray, fstripcalcBOT, fstripcalcTOP, fVerbose, PndHypStripClusterBuilder::GetBotClusterIDs(), PndHypStripDigiPar::GetBotPitch(), PndHypDigiStrip::GetChannel(), PndHypDigi::GetCharge(), PndHypDigi::GetDetID(), PndHypDigi::GetDetName(), PndHypDigi::GetFE(), PndHypStripDigiPar::GetOrient(), PndHypStripDigiPar::GetSkew(), PndHypCalcStrip::GetStripDirection(), PndHypStripClusterBuilder::GetTopClusterIDs(), PndHypStripDigiPar::GetTopPitch(), i, PndHypGeoHandling::LocalToMasterErrorsId(), PndHypGeoHandling::LocalToMasterId(), PndHypStripClusterBuilder::SearchClusters(), PndHypHit::SetBotIndex(), sin(), sqrt(), strip, t, and TString.

{
  if (fVerbose > 2)
    std::cout<<"Sarting PndHypStripClusterTask::Exec()"<<std::endl;
  std::vector<PndHypDigiStrip> digiStripArray;
//   std::vector<PndHypClusterCand> clusters;

  // Reset output array
  if ( ! fClusterArray ) Fatal("Exec", "No ClusterArray");
  fClusterArray->Clear();

  if ( ! fHitArray ) Fatal("Exec", "No HitArray");
  fHitArray->Clear();
  
  //Int_t nPoints = fDigiArray->GetEntriesFast(); //[R.K. 01/2017] unused variable

  // load the Clusterfinder
  PndHypStripClusterBuilder clusterbuilder;
 //  PndHypCalcStrip stripcalcTOP(fDigiPar, kTOP);
//   PndHypCalcStrip stripcalcBOT(fDigiPar, kBOTTOM);
  if(fVerbose>1) fDigiPar->Print();
  
  TVector2 topDirection, botDirection ;
  //TGeoVolume* actVolume; //[R.K. 01/2017] unused variable
  //TGeoBBox* actBox; //[R.K. 01/2017] unused variable
  TVector3 sensorDim,localpos, localDpos;
  TString detName;//std::string detName;
  Int_t detID, clindex, botIndex;//;iDigi,  mcindex,topIndex, //[R.K. 01/2017] unused variable
  Double_t mycharge;
  TVector2 meantopPoint, meanbotPoint, onsensorPoint;
  TVector3 hitPos,hitErr;
  Double_t t, b; //local[3], master[3], //[R.K. 01/2017] unused variable

 

  //std::map<std::string,Int_t> firedSensors;
    Int_t strip;
    SensorSide side;
  PndHypDigiStrip* myDigi=0;
  //PndHypCluster* myCandTop=0; //[R.K. 01/2017] unused variable
  //PndHypCluster* myCandBot=0; //[R.K. 01/2017] unused variable

 
  for (Int_t iPoint = 0; iPoint < fDigiArray->GetEntriesFast(); iPoint++)
    { // sort digis by sensor name and stripnumber
    
      myDigi = (PndHypDigiStrip*)(fDigiArray->At(iPoint));
    
      detName = myDigi->GetDetName().Data();
     
      fstripcalcTOP->CalcFeChToStrip(myDigi->GetFE(), myDigi->GetChannel(), strip, side);

      std::cout<<" strip "<< strip<<" side "<<std::endl;
      
      //stripcalcBOT.CalcFeChToStrip(myDigi->GetFE(), myDigi->GetChannel(), stripB, sideB);
      
     
      clusterbuilder.AddDigi(detName.Data(),side,strip,iPoint);
      
    
     
      //if(firedSensors[detName.Data()]) firedSensors[detName.Data()] = myDigi->GetDetID();
    
    }

   


  std::vector< PndHypCluster > clusters;
  std::vector< Int_t > topclusters;// contains index to fClusterArray
  std::vector< Int_t > botclusters;// contains index to fClusterArray
  std::vector< Int_t > oneclustertop;
  std::vector< Int_t > oneclusterbot;
 

    //  -----  search for clusters  -----
    clusters = clusterbuilder.SearchClusters();

    topclusters = clusterbuilder.GetTopClusterIDs();
    botclusters = clusterbuilder.GetBotClusterIDs();
    
    for(std::vector< PndHypCluster >::iterator clit= clusters.begin();
        clit!=clusters.end(); ++clit)
      {
        clindex = fClusterArray->GetEntriesFast();
        //std::cout<<" clindex "<<clindex<<std::endl;
        new((*fClusterArray)[clindex]) PndHypCluster(*clit);
      }

    if(fVerbose > 2)
    {
      std::cout<<"Top Clusters: ";
      for (std::vector< Int_t>::iterator itTop = topclusters.begin();
        itTop!=topclusters.end(); ++itTop)
      {
        std::cout<<*itTop<<" | ";
      }
      std::cout<<std::endl;
      std::cout<<"Bot Clusters: ";
      for (std::vector< Int_t>::iterator itBot = botclusters.begin();
        itBot!=botclusters.end(); ++itBot)
      {
        std::cout<<*itBot<<" | ";
      }
      std::cout<<std::endl;
    }
    //  -----  calculate mean strip numbers  -----

    // -----  merge top/bot clusters to hits  -----
    // loop on clusters from the top side
    for (std::vector< Int_t>::iterator itTop = topclusters.begin();
          itTop!=topclusters.end(); ++itTop)
    {
      // prepare values for the top side
      Double_t topcharge = 0., meantopstrip=0.;

      oneclustertop = clusters[*itTop].GetClusterList();
     
      TString detnametop = ((PndHypDigiStrip*)fDigiArray->At(oneclustertop[0]))->GetDetName();
      
      //std::cout<<"det id top "<<detnametop<<std::endl;
      topDirection = fstripcalcTOP->GetStripDirection();
      botDirection = fstripcalcBOT->GetStripDirection();

     detID = ((PndHypDigiStrip*)fDigiArray->At(oneclustertop[0]))->GetDetID();
     
      //Double_t lastcharge=0; //[R.K. 01/2017] unused variable
      for (std::vector<Int_t>::iterator itTopDigi = oneclustertop.begin();
           itTopDigi != oneclustertop.end(); ++itTopDigi)
        {// I use the temp. variables from the top of this method again
          // calculate the mean charge and stripnumber
          
          myDigi = (PndHypDigiStrip*)fDigiArray->At(*itTopDigi);

          fstripcalcTOP->CalcFeChToStrip(myDigi->GetFE(), myDigi->GetChannel(), strip, side);
          topcharge += myDigi->GetCharge();
          meantopstrip += myDigi->GetCharge() * strip;
          
          if(fVerbose > 2)std::cout<<"FE no: "<<myDigi->GetFE()<<" | channel no: "
                                   <<myDigi->GetChannel()<<" | topstrip no: "
                                   <<strip<<"| detID no: "<<myDigi->GetDetID()<<std::endl;
        }
      
      if(topcharge>0)
        {
          meantopstrip = meantopstrip/topcharge;
          fstripcalcTOP->CalcStripPointOnLine(meantopstrip, meantopPoint);
          // loop on bottom side
          for (std::vector< Int_t>::iterator itBot = botclusters.begin();
              itBot!=botclusters.end(); ++itBot)
            {
              botIndex = *itBot;
              Double_t botcharge = 0., meanbotstrip=0.;
              
              oneclusterbot =( clusters[*itBot]).GetClusterList();
              //std::cout<<" size bot "<<oneclusterbot.size()<<std::endl;
              
              TString detnamebot = ((PndHypDigiStrip*)fDigiArray->At(oneclusterbot[0]))->GetDetName();
              //std::cout<<" bot name "<<detnamebot<<std::endl;
              
              if(detnamebot != detnametop) continue;
              // cout<<"  name detID "<<detnametop<<" bot "<<detnamebot<<" ind "
              //  <<botIndex<<std::endl;
              
              for (std::vector<Int_t>::iterator itBotDigi = oneclusterbot.begin();
                   itBotDigi != oneclusterbot.end(); ++itBotDigi)
                {// I use the temp. variables from the top of this method again
           
                  myDigi = (PndHypDigiStrip*)fDigiArray->At(*itBotDigi);

                  fstripcalcTOP->CalcFeChToStrip(myDigi->GetFE(), myDigi->GetChannel(), 
                                               strip, side);
                  botcharge += myDigi->GetCharge();
                  meanbotstrip += myDigi->GetCharge() * strip;
                  
                  if(fVerbose > 2)std::cout<<"FE no: "<<myDigi->GetFE()<<" | channel no: "
                                           <<myDigi->GetChannel()<<" | botstrip no: "<<strip<<std::endl;
                }
              
              if(botcharge>0)
                {
                  meanbotstrip = meanbotstrip/botcharge;
                  // look if the charges are not too differently
                  if(fVerbose > 2)
                    {
                      std::cout<<"Charges: Ctop = "<<topcharge<<" | Cbot = "<<botcharge
                               <<" | difference bot-top = "<<botcharge-topcharge<<std::endl;
                    }
                  if(fabs(botcharge-topcharge)<fChargeCut)
                    {
                      mycharge = (botcharge + topcharge) / 2.;
                      fstripcalcBOT->CalcStripPointOnLine(meanbotstrip, meanbotPoint);
                      // get the backmapped point
                      onsensorPoint = CalcLineCross(meantopPoint,
                                                    topDirection,meanbotPoint, botDirection);
                      
                      if (fVerbose > 2)
                        {
                          std::cout<<"Strip no. top: "<<meantopstrip
                                   <<" | bot: "<<meanbotstrip
                                   <<" | Mean Top side: ("<<meantopPoint.X()<<" ; "
                                   <<meantopPoint.Y()<<")"<<"\n"<<"Mean Bot side: ("
                                   <<meanbotPoint.X()<<" ; "<<meanbotPoint.Y()<<")"<<"\n"
                                   <<"On Sensor Hit: ("<<onsensorPoint.X()<<" ; "
                                   <<onsensorPoint.Y()<<")"<<std::endl;
                        }
                      std::cout<<" botIndex @final "<<botIndex<<std::endl;
                      
                      localpos.SetXYZ(  onsensorPoint.X(),
                                        onsensorPoint.Y(),
                                        0. );
                      // local[0]=onsensorPoint.X();
                      //local[1]=onsensorPoint.Y();//0.
                      //local[2]=0.;//onsensorPoint.Y();//local coord.X-Z
                      
                      //shift sensor system into middle
                      
                      //local[0]-=sensorDim.X();
                      //local[1]-=sensorDim.Y();
                      //local[2]-=sensorDim.Z();
                      
                      //do the transformation from sensor to lab frame
                      hitPos = fGeoH->LocalToMasterId(localpos,detnametop.Data());

                      //gGeoManager->LocalToMaster(local,master);
                      //hitPos.SetXYZ(master[0],master[1],master[2]);
                      
                      // calculate the errors corresponding to a skewed system!
                      
                      t = fDigiPar->GetTopPitch()*cos(fDigiPar->GetOrient());
                      b = fDigiPar->GetBotPitch()*cos(fDigiPar->GetOrient()+fDigiPar->GetSkew());
                      
                      localDpos.SetX(sqrt((t*t+b*b)/12.));
                      
                      t = fDigiPar->GetTopPitch()*sin(fDigiPar->GetOrient());
                      b = fDigiPar->GetBotPitch()*sin(fDigiPar->GetOrient()+fDigiPar->GetSkew());
                       
                      localDpos.SetY(sqrt((t*t+b*b)/12.));
                       localDpos.SetZ(0.);
                      //x-z c.s local[2]=sqrt((t*t+b*b)/12.);
                      //local[1]=sqrt((t*t+b*b)/12.);//x-z cs.0.
                      //local[2]=0.;
                      //hitErr.SetXYZ(local[0],local[1],0.);
                      hitErr = fGeoH->LocalToMasterErrorsId(localDpos,detnametop.Data());
                      //gGeoManager->LocalToMaster(local,master);
                      //hitErr.SetXYZ(master[0],master[1],master[2]);
                     
                      
                      Int_t i = fHitArray->GetEntriesFast();
                     
                      new((*fHitArray)[i])PndHypHit(detID,detnametop.Data(),
                                                    hitPos,hitErr,*itTop,mycharge,
                                                    oneclusterbot.size()+oneclustertop.size());
                      ((PndHypHit*)((*fHitArray)[i]))->SetBotIndex(*itBot);
                      
                    }else
                    if (fVerbose > 2) std::cout<<"Strip charge contents too differently"<<std::endl;
                } 
            }// loop bot clusters
        }
    }// loop top clusters
    
    //}// end loop itSensors
    
    
    
    
    if (fVerbose > 1)
      std::cout << "ClusterArraySize: " << fClusterArray->GetEntriesFast()
                << " | HitArraySize: " << fHitArray->GetEntriesFast() << " Hits calculated."
                << " out of " <<fDigiArray->GetEntriesFast()<< " Digis"<<std::endl;
  
     
    return;
    
}

InitStatus PndHypStripClusterTask::Init ( )

virtual

Definition at line 91 of file PndHypStripClusterTask.cxx.

References fClusterArray, fDigiArray, fDigiPar, fGeoFile, fGeoH, fHitArray, fstripcalcBOT, fstripcalcTOP, kBOTTOM, kTOP, and SetParContainers().

{
  
  FairRootManager* ioman = FairRootManager::Instance();

  if ( ! ioman ) 
    {
      std::cout << "-E- PndHypStripClusterTask::Init: "
     << "RootManager not instantiated!" << std::endl;
      return kFATAL;
    }

  // Get input array
  fDigiArray = (TClonesArray*) ioman->GetObject("HypStripDigis");

  if ( ! fDigiArray ) 
  {
      std::cout << "-W- PndHypStripClusterTask::Init: "
     << "No HYPDigi array!" << std::endl;
      return kERROR;
  }

  // set output arrays
  fHitArray = new TClonesArray("PndHypHit");
  //ioman->Register("HYPStripClusterHit", "HYP", fHitArray, kTRUE);
  ioman->Register("HYPHit", "HYP", fHitArray, kTRUE);
  
  fClusterArray = new TClonesArray("PndHypCluster");
  ioman->Register("HYPStripClusterCand","HYP",fClusterArray,kTRUE);

  SetParContainers();

  // geo name handling
  if (fGeoFile=="") fGeoFile = ioman->GetInFile()->GetName();
  fGeoH = new PndHypGeoHandling(fGeoFile.Data());

 if ( ! fDigiPar )
   {
     std::cout << "-W- PndMvdStripClusterTask::Init: "
               << "No fDigiParTrap!" << std::endl;
     return kERROR;
   }
 fstripcalcTOP = new PndHypCalcStrip (fDigiPar, kTOP);
 fstripcalcBOT  = new PndHypCalcStrip (fDigiPar, kBOTTOM);

  std::cout << "-I- PndHypStripClusterTask: Initialisation successfull" << std::endl;
  return kSUCCESS;
}

void PndHypStripClusterTask::ProduceHits ( )

private

void PndHypStripClusterTask::Register ( )

private

InitStatus PndHypStripClusterTask::ReInit ( )

virtual

Definition at line 73 of file PndHypStripClusterTask.cxx.

{
  
  InitStatus stat=kERROR;
  return stat;
  
  
  /*FairRunAna* ana = FairRunAna::Instance();
  FairRuntimeDb* rtdb=ana->GetRuntimeDb();
 fGeoPar=(PndGeoHypPar*)(rtdb->getContainer("PndGeoHypPar"));
  return kSUCCESS;*/
 
  /*
  return kSUCCESS;
  */
}

void PndHypStripClusterTask::Reset ( )

private

void PndHypStripClusterTask::SetParContainers ( )

virtual

Virtual method Init

Definition at line 63 of file PndHypStripClusterTask.cxx.

References fDigiPar, and rtdb.

Referenced by Init().

{
  // Get Base Container
  
  FairRunAna* ana = FairRunAna::Instance();
  FairRuntimeDb* rtdb=ana->GetRuntimeDb();
  //fGeoPar = (PndGeoHypPar*)(rtdb->getContainer("PndGeoHypPar"));
  fDigiPar = (PndHypStripDigiPar*)(rtdb->getContainer("PndHypStripDigiPar"));
}

Member Data Documentation

Double_t PndHypStripClusterTask::fChargeCut

private

Definition at line 71 of file PndHypStripClusterTask.h.

Referenced by Exec(), and PndHypStripClusterTask().

TClonesArray* PndHypStripClusterTask::fClusterArray

private

Definition at line 60 of file PndHypStripClusterTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndHypStripClusterTask::fDigiArray

private

Definition at line 59 of file PndHypStripClusterTask.h.

Referenced by Exec(), and Init().

PndHypStripDigiPar* PndHypStripClusterTask::fDigiPar

private

Definition at line 73 of file PndHypStripClusterTask.h.

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

Int_t PndHypStripClusterTask::fFEcolumns

private

Definition at line 68 of file PndHypStripClusterTask.h.

Int_t PndHypStripClusterTask::fFErows

private

Definition at line 69 of file PndHypStripClusterTask.h.

TString PndHypStripClusterTask::fGeoFile

private

Definition at line 70 of file PndHypStripClusterTask.h.

Referenced by Init(), and PndHypStripClusterTask().

PndHypGeoHandling* PndHypStripClusterTask::fGeoH

private

Definition at line 76 of file PndHypStripClusterTask.h.

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

PndGeoHypPar* PndHypStripClusterTask::fGeoPar

private

Geometry name handling.

Definition at line 77 of file PndHypStripClusterTask.h.

TClonesArray* PndHypStripClusterTask::fHitArray

private

Definition at line 61 of file PndHypStripClusterTask.h.

Referenced by Exec(), and Init().

Double_t PndHypStripClusterTask::fRadius

private

Definition at line 67 of file PndHypStripClusterTask.h.

PndHypCalcStrip* PndHypStripClusterTask::fstripcalcBOT

private

Definition at line 75 of file PndHypStripClusterTask.h.

Referenced by Exec(), and Init().

PndHypCalcStrip* PndHypStripClusterTask::fstripcalcTOP

private

Digitization Parameters.

Definition at line 74 of file PndHypStripClusterTask.h.

Referenced by Exec(), and Init().

---

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

• PndHypStripClusterTask.h
• PndHypStripClusterTask.cxx