FairRoot/PandaRoot
Functions | Variables
anasim.C File Reference

Go to the source code of this file.

Functions

gROOT LoadMacro ("../Tools.C")
 
 LoadPandaStyle ()
 
timer Start ()
 
picture ReplaceAll (".root",".ps")
 
t SetBranchAddress ("MVDPoint",&hit_array)
 
t SetBranchAddress ("MCTrack",&mc_array)
 
hisdedx SetTitle ("MVD MC Points, dE/dx(P);p/GeV / cm;dE/dx / (GeV/cm)")
 
hisLocalXY SetTitle ("Local Hit XY;x_{L} / cm;y_{L} / cm")
 
hisLocalZ SetTitle ("Locel MC Hit Z;z / cm")
 
 for (Int_t j=0;j< nEvents &&j< t->GetEntriesFast();j++)
 
can1 Divide (3, 3)
 
can1 cd (1)
 
 DrawNice2DHisto (hisxy)
 
can1 cd (2)
 
 DrawNice2DHisto (hisrz)
 
can1 cd (3)
 
gPad SetLogy ()
 
hisde DrawCopy ("")
 
can1 cd (4)
 
 DrawNice2DHisto (hisdedx)
 
can1 cd (5)
 
 DrawNice2DHisto (hisLocalXY)
 
can1 cd (6)
 
can1 cd (7)
 
can1 Print (picture.Data())
 
timer Stop ()
 

Variables

TStopwatch timer
 
PndFileNameCreator namecreator ("../data/mvdStrip.root")
 
std::string inFile = namecreator.GetSimFileName(false)
 
TString picture = namecreator.GetSimFileName(false)
 
TFile * f = new TFile(inFile.c_str())
 
TTree * t =(TTree *) f->Get("pndsim")
 
TClonesArray * hit_array =new TClonesArray("PndSdsMCPoint")
 
TClonesArray * mc_array =new TClonesArray("PndMCTrack")
 
TGeoManager * geoMan = (TGeoManager*) gDirectory->Get("FAIRGeom")
 
PndGeoHandlingfGeoH = new PndGeoHandling()
 
TH2D * hisxy = new TH2D("hisxy","MVD MC Points, xy view",200,-15.,15.,200,-15.,15.)
 
TH2D * hisrz = new TH2D("hisrz","MVD MC Points, rz view",200,-20.,20.,200,-15.,25.)
 
TH1D * hisde = new TH1D("hisde","MVD MC Points, Energyloss",200,0.,0.002)
 
TH2D * hisdedx =new TH2D("hisdedx","",200,0.,0.4,200,0.,1.0)
 
TH1D * hismom = new TH1D("hismom","MVD MC Points, momentum",100,0.,1.5)
 
TH2D * hisLocalXY = new TH2D("hisLocalxy","",100,-5.,5.,100,-5.,5.)
 
TH1D * hisLocalZ = new TH1D("localz","",100,-0.005, 0.005)
 
int nEvents = 1000
 
bool verbose = false
 
TVector3 vecs
 
TVector3 veco
 
TVector3 vecFront
 
TVector3 vecBack
 
TVector3 vecP
 
Double_t dx
 
Double_t dE
 
Double_t p
 
Double_t dEdX
 
Double_t tmpMaster [3]
 
Double_t tmpLocal [3]
 
TGeoHMatrix * currentTransMat
 
TString detname
 
TCanvas * can1 = new TCanvas("can1","MCHit view in MVD",0,0,800,800)
 
Double_t rtime = timer.RealTime()
 
Double_t ctime = timer.CpuTime()
 

Function Documentation

can1 cd ( )
can1 cd ( )
can1 cd ( )
can1 cd ( )
can1 cd ( )
can1 cd ( )
can1 cd ( )
can1 Divide ( ,
 
)
hisde DrawCopy ( ""  )
DrawNice2DHisto ( hisxy  )
DrawNice2DHisto ( hisrz  )
DrawNice2DHisto ( hisdedx  )
DrawNice2DHisto ( hisLocalXY  )
for ( )

Definition at line 63 of file anasim.C.

References dx, PndGeoHandling::GetPath(), PndSdsMCPoint::GetXOut(), PndSdsMCPoint::GetYOut(), PndSdsMCPoint::GetZOut(), hit(), i, mcpdg, and vecs.

64  {
65  t->GetEntry(j);
66  if( verbose || 0 == j%1000 ) cout<<"Event No "<<j<<endl;
67  for (Int_t i=0; i<hit_array->GetEntriesFast(); i++)
68  {
69  if(verbose) cout<<"Point No "<<i<<endl;
71  int mcpdg = -1;
72 
73  //PndMCTrack *mctruth = (PndMCTrack*)mc_array->At(hit->GetTrackID());
74  //mcpdg = mctruth->GetPdgCode();
75  //cout<<"mcpdg="<<mcpdg<<endl;
76 
77  vecs.SetXYZ(hit->GetX(), hit->GetY(), hit->GetZ());
78  veco.SetXYZ(hit->GetXOut(),hit->GetYOut(),hit->GetZOut());
79  vecP.SetXYZ(hit->GetPx(),hit->GetPy(),hit->GetPz());
80  dx=(veco-vecs).Mag();
81  dE=hit->GetEnergyLoss();
82  p=vecP.Mag();
83  Int_t layer = Int_t(10.*vecs->Mag());
84  if(verbose) cout<<vecs.x()<<"\t"<<vecs.y()<<"\t"<<vecs.z()<<endl;
85 
86  hisxy->Fill(vecs.x(),vecs.y());
87  hisrz->Fill(vecs.z(),((vecs.y()>0.)?1.:-1.)*vecs.Perp());
88  hisde->Fill(dE);
89  hismom->Fill(p);
90 
91  if(dx!=0)
92  {
93  dEdX=(dE/dx);
94  hisdedx->Fill(p,dEdX);
95  }
96 
97  detname = fGeoH->GetPath( hit->GetDetName());
98  geoMan->cd( detname.Data() );
99  currentTransMat = geoMan->GetCurrentMatrix();
100  tmpMaster[0]=0.5*(veco.x()+vecs.x());
101  tmpMaster[1]=0.5*(veco.y()+vecs.y());
102  tmpMaster[2]=0.5*(veco.z()+vecs.z());
103  currentTransMat->MasterToLocal(tmpMaster,tmpLocal);
104  vecs.SetXYZ(tmpLocal[0],tmpLocal[1],tmpLocal[2]);
105  hisLocalXY->Fill(vecs.x(),vecs.y());
106  hisLocalZ->Fill(vecs.z());
107 
108  }//end for i (points in event)
109  }// end for j (events)
TH2D * hisLocalXY
Definition: anasim.C:48
TVector3 veco
Definition: anasim.C:56
Double_t GetXOut() const
Definition: PndSdsMCPoint.h:81
TH2D * hisdedx
Definition: anasim.C:42
TClonesArray * hit_array
Definition: anasim.C:27
Int_t i
Definition: run_full.C:25
PndGeoHandling * fGeoH
Definition: anasim.C:34
Double_t GetZOut() const
Definition: PndSdsMCPoint.h:83
Double_t tmpMaster[3]
Definition: anasim.C:59
bool verbose
Definition: anasim.C:54
TString GetPath(Int_t shortID)
for a given shortID the path is returned
Double_t p
Definition: anasim.C:58
Double_t dx
Definition: anasim.C:58
Double_t dEdX
Definition: anasim.C:58
TH2D * hisrz
Definition: anasim.C:38
TGeoManager * geoMan
Definition: anasim.C:33
TH1D * hismom
Definition: anasim.C:47
Double_t dE
Definition: anasim.C:58
TVector3 vecs
Definition: anasim.C:56
TString detname
Definition: anasim.C:61
TVector3 vecP
Definition: anasim.C:57
int mcpdg
TH1D * hisde
Definition: anasim.C:39
int hit(Int_t nEvents=0, TString inFile="sim.root", TString parFile="par.root", TString inDigi="digi.root", TString outFile="hit.root", Int_t timeBased=0)
Definition: hit.C:1
Double_t GetYOut() const
Definition: PndSdsMCPoint.h:82
TGeoHMatrix * currentTransMat
Definition: anasim.C:60
Double_t tmpLocal[3]
Definition: anasim.C:59
TH1D * hisLocalZ
Definition: anasim.C:50
TH2D * hisxy
Definition: anasim.C:37
TTree * t
Definition: anasim.C:26
gROOT LoadMacro ( "../Tools.C"  )
LoadPandaStyle ( )
can1 Print ( picture.  Data())
picture ReplaceAll ( ".root"  ,
".ps"   
)
t SetBranchAddress ( "MVDPoint"  ,
hit_array 
)
t SetBranchAddress ( "MCTrack"  ,
mc_array 
)
gPad SetLogy ( )
hisdedx SetTitle ( "MVD MC  Points,
dE/dx(P);p/GeV/cm;dE/dx/(GeV/cm)"   
)
hisLocalXY SetTitle ( "Local Hit XY;x_{L} / cm;y_{L} / cm )
hisLocalZ SetTitle ( "Locel MC Hit Z;z / cm )
timer Start ( )
timer Stop ( )

Variable Documentation

TCanvas* can1 = new TCanvas("can1","MCHit view in MVD",0,0,800,800)

Definition at line 111 of file anasim.C.

Double_t ctime = timer.CpuTime()

Definition at line 138 of file anasim.C.

TGeoHMatrix* currentTransMat

Definition at line 60 of file anasim.C.

Referenced by QAmacro_mvd_ana().

Double_t dEdX
TString detname

Definition at line 58 of file anasim.C.

Referenced by for().

TFile* f = new TFile(inFile.c_str())

Definition at line 25 of file anasim.C.

Definition at line 34 of file anasim.C.

Referenced by anaclust(), ClassImp(), QAmacro_mvd_ana(), ReadHCal(), ReadHits(), ReadLines(), and TrackHits().

TGeoManager* geoMan = (TGeoManager*) gDirectory->Get("FAIRGeom")

Definition at line 33 of file anasim.C.

TH1D* hisde = new TH1D("hisde","MVD MC Points, Energyloss",200,0.,0.002)

Definition at line 39 of file anasim.C.

TH2D* hisdedx =new TH2D("hisdedx","",200,0.,0.4,200,0.,1.0)

Definition at line 42 of file anasim.C.

Referenced by ana_MCOpt(), and ana_MCpid().

TH2D* hisLocalXY = new TH2D("hisLocalxy","",100,-5.,5.,100,-5.,5.)

Definition at line 48 of file anasim.C.

TH1D* hisLocalZ = new TH1D("localz","",100,-0.005, 0.005)

Definition at line 50 of file anasim.C.

TH1D* hismom = new TH1D("hismom","MVD MC Points, momentum",100,0.,1.5)

Definition at line 47 of file anasim.C.

TH2D* hisrz = new TH2D("hisrz","MVD MC Points, rz view",200,-20.,20.,200,-15.,25.)

Definition at line 38 of file anasim.C.

TH2D* hisxy = new TH2D("hisxy","MVD MC Points, xy view",200,-15.,15.,200,-15.,15.)

Definition at line 37 of file anasim.C.

TClonesArray* hit_array =new TClonesArray("PndSdsMCPoint")

Definition at line 27 of file anasim.C.

std::string inFile = namecreator.GetSimFileName(false)

Definition at line 21 of file anasim.C.

TClonesArray* mc_array =new TClonesArray("PndMCTrack")

Definition at line 30 of file anasim.C.

PndFileNameCreator namecreator("../data/mvdStrip.root")
int nEvents = 1000

Definition at line 53 of file anasim.C.

Definition at line 58 of file anasim.C.

Referenced by RhoSimpleVertexSelector::Accept(), RhoGoodPhotonSelector::Accept(), RhoGoodTrackSelector::Accept(), DecayTreeFitter::FitParams::addChiSquare(), PndCAInternal::Allocator< T, alignment >::Alloc(), PndFTSInternal::Allocator< T, alignment >::Alloc(), PndCAInternal::Allocator< T, PndCAFullyCacheLineAligned >::Alloc(), PndFTSInternal::Allocator< T, PndFTSFullyCacheLineAligned >::Alloc(), PndCAInternal::Allocator< T, 0 >::Alloc(), PndFTSInternal::Allocator< T, 0 >::Alloc(), ana_dsdsj_full(), ana_Lambda(), ana_Lambda_fit(), ana_Lambda_fit2(), ana_MCOpt(), ana_MCpid(), ana_pid(), AnalyseSimFile(), Photospp::PhotosParticle::boostAlongZ(), PndRiemannTrack::calcErrorXY1XY2(), PndRiemannTrack::calcIntersection(), PndBarrelTrackFinder::CalcPhi(), calculate_ratio(), PndFTSCAParam::CalculateFieldSlice(), checkMomentumConservationInEvent(), DecayTreeFitter::RecoParticle::chiSquare(), DecayTreeFitter::RecoComposite::chiSquare(), PndTrkTrack::ComputePhi(), PndTrkTools::ComputeSegmentCircleIntersection(), nsL1::SimdAlloc< T >::construct(), KFParticleBase::ConstructGammaBz(), KFParticleBaseSIMD::ConstructGammaBz(), AliKFVertex::ConstructPrimaryVertex(), KFVertex::ConstructPrimaryVertex(), PndTrkLegendreSecTask::CorrectZ(), PndTrkLegendreTask::CorrectZ(), PndTrkLegendreSecTask2::CorrectZ(), Photospp::PhotosBranch::createBranches(), PndGemContFact::createContainer(), PndTutContFact::createContainer(), PndAnaContFact::createContainer(), PndDrcContFact::createContainer(), CbmPlaneContFact::createContainer(), PndPassiveContFact::createContainer(), PndFtsContFact::createContainer(), PndMdtContFact::createContainer(), PndMvdContFact::createContainer(), PndRichContFact::createContainer(), PndSttContFact::createContainer(), PndSensorNameContFact::createContainer(), PndLmdContFact::createContainer(), PndFtofContFact::createContainer(), PndHypContFact::createContainer(), PndHypGeContFact::createContainer(), PndEmcContFact::createContainer(), PndPidContFact::createContainer(), PndSciTContFact::createContainer(), PndGemDigiPar::CreateStations(), createstt(), PndFastSim::cutAndSmear(), nsL1::SimdAlloc< T >::deallocate(), chigen::PythiaChiGen::decay(), PndFsmSimpleTracker::detected(), PndFsmRich::detected(), PndFsmMvd2::detected(), PndFsmMdcFS::detected(), PndFsmMdcTS::detected(), PndFsmMvd::detected(), PndFsmDrcDisc::detected(), PndFsmMdtPid::detected(), PndFsmDrcBarrel::detected(), PndFsmEffTracker::detected(), PndFsmEmcPid::detected(), GFDetPlane::dist(), PndFsmSimpleTracker::dp(), PndFsmDetTemplate::dp(), PndFsmMvd2::dp(), PndFsmMvd::dp(), PndFsmEffTracker::dp(), PndFsmSimpleTracker::dphi(), PndFsmEffTracker::dphi(), PndFTSCADisplay::DrawGBHits(), PndFTSCADisplay::DrawHelix(), PndFTSCADisplay::DrawParticleGlobal(), PndFTSCADisplay::DrawRecoTrack(), PndFsmSimpleTracker::dtheta(), PndFsmEffTracker::dtheta(), PndRichRecoTask::Exec(), PndMvdDigiPixelDraw::Exec(), PndDiscTaskReconstruction::Exec(), PndKFParticleFinderQA::Exec(), PndDiscTaskPID::Exec(), PndFtsDataAccessor::Exec(), PndFtsCATracking::Exec(), PndLLbarAnaTask::Exec(), PndEmcApdHitProducer::Exec(), PndLmdBPRungeKuttaTask::Exec(), PndHypFullIdealAna::Exec(), PndHypFullAna::Exec(), PndHypSimpleAna::Exec(), PndCATracking::Exec(), PndHypMicroIdealWriter::Exec(), PndMvdEventAnaTask::Exec(), PndMvdMSAnaTask::Exec(), PndEmcHitProducer::Exec(), PndHypKalmanTask::Exec(), PndHypDKalmanTask::Exec(), PndSimpleCombinerTask::Exec(), PndEmcPhiBumpSplitter::Exec(), PndCaloDraw::Exec(), DecayTreeFitter::Constraint::filter(), Photospp::PhotosEvent::filterParticles(), KFParticleFinder::Find2DaughterDecay(), PndRiemannTrackFinder::FindTracks(), PndMvdSttGemRiemannTrackFinder::FindTracks(), DecayTreeFitter::Fitter::fitParams(), PndEmcPSAParabolic::FitPeak(), PndEmcPSAParabolicBaseline::FitPeak(), fixForMctester(), G__CbmPlaneDict_202_0_1(), G__CbmPlaneDict_202_0_14(), G__CbmPlaneDict_209_0_1(), G__CbmPlaneDict_209_0_19(), G__CbmPlaneDict_560_0_1(), G__CbmPlaneDict_560_0_2(), G__CbmPlaneDict_560_0_3(), G__CbmPlaneDict_561_0_1(), G__CbmPlaneDict_561_0_2(), G__CbmPlaneDict_561_0_30(), G__CbmPlaneDict_565_0_16(), G__CbmPlaneDict_565_0_2(), G__setup_memvarCbmGeoPlane(), G__setup_memvarCbmGeoPlanePar(), G__setup_memvarCbmPlane(), G__setup_memvarCbmPlaneContFact(), G__setup_memvarCbmPlanePoint(), G__Sizep2memfuncCbmPlaneDict::G__Sizep2memfuncCbmPlaneDict(), Photospp::PH_HEPEVT_Interface::get(), PndEmcErrorMatrix::Get4MomentumErrorMatrix(), get_trap_vertices(), PndEmcPSAParabolicBaseline::GetBaseline(), PndRichCalDb::GetBetaEff(), PndRichCalDb::GetBetaMean(), PndRichCalDb::GetBetaSig(), PndRiemannTrack::getCharge(), Photospp::PhotosHEPEVTParticle::getDaughters(), Photospp::PhotosParticle::getDecayTree(), PndGemDigiPar::GetDetectorIdByName(), PndEmcErrorMatrix::GetErrorP7(), KFParticleBase::GetEta(), KFParticleBaseSIMD::GetEta(), PndMvdEventAnaTask::GetFairHit(), PndForwardTrackFinderTask::getMcId(), KFParticleBase::GetMomentum(), RKTrackRep::getPosMomCov(), PndGemDigiPar::GetSensorByName(), RhoCandList::GetTotalMomentum(), PndMQStraightLineTrackFinder::GetTrackCurvature(), PndHypIdealTrackingTask::GetTrackCurvature(), PndHypIdealTrackFinderTask::GetTrackCurvature(), PndLmdTrackFinderTask::GetTrackCurvature(), PndMvdIdealTrackFinderTask::GetTrackCurvature(), PndStraightLineTrackFinderTask::GetTrackCurvature(), PndLmdTrackFinderCATask::GetTrackCurvature(), PndMQStraightLineTrackFinder::GetTrackDip(), PndHypIdealTrackingTask::GetTrackDip(), PndHypIdealTrackFinderTask::GetTrackDip(), PndLmdTrackFinderTask::GetTrackDip(), PndMvdIdealTrackFinderTask::GetTrackDip(), PndStraightLineTrackFinderTask::GetTrackDip(), PndLmdTrackFinderCATask::GetTrackDip(), PndSttMCPointDraw::GetVector(), PndTrackCandDraw::GetVector(), PndTrkTracking2::GetVolumeCharacteristics(), PndFsmTrack::HelixRep(), PndDrcReco::Init(), main(), makeIni4Vector(), ROOT::new_CbmGeoPlane(), ROOT::new_CbmGeoPlanePar(), ROOT::new_CbmPlane(), ROOT::new_CbmPlaneContFact(), ROOT::new_CbmPlanePoint(), ROOT::newArray_CbmGeoPlane(), ROOT::newArray_CbmGeoPlanePar(), ROOT::newArray_CbmPlane(), ROOT::newArray_CbmPlaneContFact(), ROOT::newArray_CbmPlanePoint(), PndFTSCAMCTrack::operator=(), operator>>(), PndEmcClusterMoments::Phi1(), PndEmcClusterMoments::Phi2(), Photospp::PhotosBranch::PhotosBranch(), PndCAGBTracker::PickUpHits(), plot_radmap(), chigen::PndChiGenExclusive::PndChiGenExclusive(), PndPmtPoormantracks::PoorManTracks(), PoorManTracks(), GFBookkeeping::Print(), PndPidCandidate::PrintOn(), PndEmcPSAMatchedDigiFilter::Process(), PndSttMvdGemTracking::PropagateToGemPlaneAsHelix(), PndRhoTupleQA::qaCand(), chigen::PndChiGenExclusive::ReadEvent(), PndCorrDistGenerator::ReadEvent(), PndHypBupGenerator::ReadEvent(), ReadLines(), KFPTopoReconstructor::ReconstructPrimVertex(), PndTrackCollection::refitAllTracks(), PndTrackCollection::refitHit(), PndFsmCmpDet::respond(), PndFsmCombiDet::respond(), PndFsmRich::respond(), PndFsmStt::respond(), PndFsmSttPid::respond(), PndFsmMvd::respond(), PndFsmTof::respond(), PndFsmEmcPid::respond(), PndFsmMdtPid::respond(), PndFsmDrcBarrel::respond(), PndFsmDrcDisc::respond(), PndFsmEffTracker::respond(), run_ana_eta_c_tpc(), runMC_dpm(), PndStack::SelectTracks(), PndDrcContFact::setAllContainers(), CbmPlaneContFact::setAllContainers(), PndPassiveContFact::setAllContainers(), PndRichContFact::setAllContainers(), PndGemContFact::setAllContainers(), PndFtofContFact::setAllContainers(), PndSciTContFact::setAllContainers(), PndSensorNameContFact::setAllContainers(), chigen::PndChiGen::setChiGen(), GFPlanarHitPolicy::setDetPlane(), PndFastSim::SetFlatCovMatrix(), PndFTSCAMCTrack::SetMCTrack(), Photospp::PhotosHEPEVTEvent::setParticle(), DecayTreeFitter::Projection::setParticle(), PndTrackArrayMerger::SetPersistance(), PndMvdRadDamIonizingTask::SetPersistance(), PndMvdRadDamTask::SetPersistance(), PndMvdIdealRecoTask::SetPersistance(), PndLmdPixelHitProducerFast::SetPersistance(), PndSorterTask::SetPersistance(), PndSdsHitProducerIdeal::SetPersistance(), PndMvdConvertApvTask::SetPersistance(), PndXYResidualTask::SetPersistance(), PndPosCorrectorTask::SetPersistance(), PndTSCorrectorTask::SetPersistance(), PndMapSorterTask::SetPersistance(), PndTimeOffsetFinderTask::SetPersistance(), PndSdsDetector::SetPersistance(), RhoColumn::SetPointer(), chigen::ChiGen::setPythiaChiGen(), PndPidBremCorrector::SetStorageOfData(), PndEmc2DLocMaxFinder::SetStorageOfData(), PndEmcExpClusterSplitter::SetStorageOfData(), PndEmcPhiBumpSplitter::SetStorageOfData(), PndSimpleAnalysis::SetupAnalysis(), RhoBoolColumn::SetValue(), RhoIntColumn::SetValue(), RhoFloatColumn::SetValue(), RhoDoubleColumn::SetValue(), simfast_dpm(), simfast_dpm_cmp(), G__Sizep2memfuncCbmPlaneDict::sizep2memfunc(), SmartAli(), smearMom(), softtrigger_kin5(), softtrigger_toy12(), RhoCalculationTools::StateFromTrajectory(), switch_history_entries_status(), toy_core(), trapmap_barrel(), trapmap_disc(), UserDecayConfig(), VarCorrP(), Vec< T >::Vec(), PndSttHelixTrackFitter::XYFit(), PndSttMvdGemTracking::ZFit(), and Photospp::PhotosHepMCEvent::~PhotosHepMCEvent().

TString picture = namecreator.GetSimFileName(false)

Definition at line 22 of file anasim.C.

Double_t rtime = timer.RealTime()

Definition at line 137 of file anasim.C.

TTree* t =(TTree *) f->Get("pndsim")

Definition at line 26 of file anasim.C.

TStopwatch timer

Definition at line 11 of file anasim.C.

Double_t tmpLocal[3]

Definition at line 59 of file anasim.C.

Referenced by QAmacro_mvd_ana().

Double_t tmpMaster[3]

Definition at line 59 of file anasim.C.

Referenced by QAmacro_mvd_ana().

TVector3 vecBack

Definition at line 57 of file anasim.C.

Referenced by ana_MCpid().

TVector3 vecFront

Definition at line 57 of file anasim.C.

Referenced by ana_MCOpt(), and ana_MCpid().

TVector3 veco

Definition at line 56 of file anasim.C.

TVector3 vecP

Definition at line 57 of file anasim.C.

Referenced by ana_MCOpt(), and ana_MCpid().

TVector3 vecs

Definition at line 56 of file anasim.C.

Referenced by for().

bool verbose = false

Definition at line 54 of file anasim.C.